38,553 research outputs found
An Explicit Finite-Difference Scheme for Simulation of Moving Particles
We present an explicit finite-difference scheme for direct simulation of the motion of solid particles in a fluid. The method is based on a second order MacCormack finite-difference solver for the flow, and Newton’s equations for the particles. The fluid is modeled with fully compressible mass and momentum balances; the technique is intended to be used at moderate particle Reynolds number. Several examples are shown, including a single stationary circular particle in a uniform flow between two moving walls, a particle dropped in a stationary fluid at particle Reynolds number of 20, the drafting, kissing, and tumbling of two particles, and 100 particles falling in a closed box
Analytic-Numerical Modeling and Investigation of Nanostructures' Dynamics on Material Surfaces After Laser Irradiation
The work is devoted to modeling the formation and behavior of solid nano-sized particles on the surface of materials. In the simulation, it is assumed that the main processing technology of surface nanostructures is laser irradiation, which causes the Brownian motion of nanoparticles, due mainly to thermal fluctuations: if the temperature around the nanoparticles is uniformly distributed, the time average of the Brownian fluctuations is zero; however, if there is a temperature gradient around the nanoparticles, the thermal fluctuations affect the nanoparticle in different ways from different sides, and there is a force like the thermophoretic force, biasing the average position of the nanoparticle. When building a 1D model of the formation and flow behavior of nanoparticles, three important assumptions are introduced: the impact of nanoparticles on the process of irradiation is negligible; the impact of nanoparticles on each other as compared to the effect of laser irradiation on them is also negligible; and nanoparticles after laser irradiation can move both forward and backward and at every fixed period of time, moving the nanoparticles does not impose any steric constraints. Under the above assumptions, a 1D continuous model is built, implicit and explicit finite difference schemes to solve it are developed; their convergence and order of convergence are studied; an output condition ensuring the stability of the explicit difference scheme is obtained, the unconditional stability of implicit difference scheme is proven, and software for computer implementation of some of the obtained analytical and numerical results developed.
A Simulation Method to Resolve Hydrodynamic Interactions in Colloidal Dispersions
A new computational method is presented to resolve hydrodynamic interactions
acting on solid particles immersed in incompressible host fluids. In this
method, boundaries between solid particles and host fluids are replaced with a
continuous interface by assuming a smoothed profile. This enabled us to
calculate hydrodynamic interactions both efficiently and accurately, without
neglecting many-body interactions. The validity of the method was tested by
calculating the drag force acting on a single cylindrical rod moving in an
incompressible Newtonian fluid. This method was then applied in order to
simulate sedimentation process of colloidal dispersions.Comment: 7pages, 7 figure
Diagnosing numerical Cherenkov instabilities in relativistic plasma simulations based on general meshes
Numerical Cherenkov radiation (NCR) or instability is a detrimental effect
frequently found in electromagnetic particle-in-cell (EM-PIC) simulations
involving relativistic plasma beams. NCR is caused by spurious coupling between
electromagnetic-field modes and multiple beam resonances. This coupling may
result from the slow down of poorly-resolved waves due to numerical (grid)
dispersion and from aliasing mechanisms. NCR has been studied in the past for
finite-difference-based EM-PIC algorithms on regular (structured) meshes with
rectangular elements. In this work, we extend the analysis of NCR to
finite-element-based EM-PIC algorithms implemented on unstructured meshes. The
influence of different mesh element shapes and mesh layouts on NCR is studied.
Analytic predictions are compared against results from finite-element-based
EM-PIC simulations of relativistic plasma beams on various mesh types.Comment: 31 pages, 20 figure
Variational Formulation of Macro-Particle Models for Electromagnetic Plasma Simulations
A variational method is used to derive a self-consistent macro-particle model
for relativistic electromagnetic kinetic plasma simulations. Extending earlier
work [E. G. Evstatiev and B. A. Shadwick, J. Comput. Phys., vol. 245, pp.
376-398, 2013], the discretization of the electromagnetic Low Lagrangian is
performed via a reduction of the phase-space distribution function onto a
collection of finite-sized macro-particles of arbitrary shape and
discretization of field quantities onto a spatial grid. This approach may be
used with both lab frame coordinates or moving window coordinates; the latter
can greatly improve computational efficiency for studying some types of
laser-plasma interactions. The primary advantage of the variational approach is
the preservation of Lagrangian symmetries, which in our case leads to energy
conservation and thus avoids difficulties with grid heating. Additionally, this
approach decouples particle size from grid spacing and relaxes restrictions on
particle shape, leading to low numerical noise. The variational approach also
guarantees consistent approximations in the equations of motion and is amenable
to higher order methods in both space and time. We restrict our attention to
the 1-1/2 dimensional case (one coordinate and two momenta). Simulations are
performed with the new models and demonstrate energy conservation and low
noise.Comment: IEEE Transaction on Plasma Science (TPS) Special Issue: Plenary and
Invited Papers of the Pulsed Power and Plasma Science Conference (PPPS 2013
- …