No abstract

No abstract

No abstract

To investigate dynamic behaviors of monocharged particle systems, a direct truncation (DT) method and a hybrid particle-cell (HPC) method are implemented into the discrete element method coupled with computational fluid dynamics (DEM-CFD) with defined cutoff distances. The DT method only considers electrostatic interactions between particles within the cutoff distance while the HPC method computes electrostatic interactions in the entire computational domain. The deposition process of monocharged particles in a container in air was simulated using the developed DEM-CFD. It was found that using the DT method, the macrostructure, evolution of granular temperature, and radial distribution function of the particle system were sensitive to the specified cutoff distance. In contrast, using the HPC method, these results were independent of the specified cutoff distance, as expected. This implies that, although electrostatic interactions between particles with large separation distances are weak, they should be considered in DEM-CFD for accurate modeling of charged particle systems.

Previous studies of the stable atmospheric boundary layer using techniques of nonlinear dynamical systems have shown that the equations support multiple solutions in certain parameter spaces. When geostrophic speed is used as a bifurcation parameter, two stable equilibria are found-a warm solution corresponding to the high-wind regime where the surface layer of the atmosphere stays coupled to the outer layer, and a cold solution corresponding to the low-wind, decoupled case. Between the stable equilibria is an unstable region where multiple solutions exist. The bifurcation diagram is a classic S shape with the foldback region showing the multiple solutions. These studies were carried out using a simple two-layer model of the atmosphere with a fairly complete surface energy budget. This allowed the dynamical analysis to be carried out on a coupled set of four ordinary differential equations. The present paper extends this work by examining additional bifurcation parameters and, more importantly, analyzing a set of partial differential equations with full vertical dependence. Simple mathematical representations of classical problems in dynamical analysis often exhibit interesting behavior, such as multiple solutions, that is not retained in the behavior of more complete representations. In the present case the S-shaped bifurcation diagram remains with only slight variations from the two-layer model. For the parameter space in the foldback region, the evolution of the boundary layer may be dramatically affected by the initial conditions at sunset. An eigenvalue analysis carried out to determine whether the system might support pure limit-cycle behavior showed that purely complex eigenvalues are not found. Thus, any cyclic behavior is likely to be transient.