1,659 research outputs found
A momentum-conserving, consistent, Volume-of-Fluid method for incompressible flow on staggered grids
The computation of flows with large density contrasts is notoriously
difficult. To alleviate the difficulty we consider a consistent mass and
momentum-conserving discretization of the Navier-Stokes equation.
Incompressible flow with capillary forces is modelled and the discretization is
performed on a staggered grid of Marker and Cell type. The Volume-of-Fluid
method is used to track the interface and a Height-Function method is used to
compute surface tension. The advection of the volume fraction is performed
using either the Lagrangian-Explicit / CIAM (Calcul d'Interface Affine par
Morceaux) method or the Weymouth and Yue (WY) Eulerian-Implicit method. The WY
method conserves fluid mass to machine accuracy provided incompressiblity is
satisfied which leads to a method that is both momentum and mass-conserving. To
improve the stability of these methods momentum fluxes are advected in a manner
"consistent" with the volume-fraction fluxes, that is a discontinuity of the
momentum is advected at the same speed as a discontinuity of the density. To
find the density on the staggered cells on which the velocity is centered, an
auxiliary reconstruction of the density is performed. The method is tested for
a droplet without surface tension in uniform flow, for a droplet suddenly
accelerated in a carrying gas at rest at very large density ratio without
viscosity or surface tension, for the Kelvin-Helmholtz instability, for a
falling raindrop and for an atomizing flow in air-water conditions
Eulerian and modified Lagrangian approaches to multi-dimensional condensation and collection
Turbulence is argued to play a crucial role in cloud droplet growth. The
combined problem of turbulence and cloud droplet growth is numerically
challenging. Here, an Eulerian scheme based on the Smoluchowski equation is
compared with two Lagrangian superparticle (or su- perdroplet) schemes in the
presence of condensation and collection. The growth processes are studied
either separately or in combination using either two-dimensional turbulence, a
steady flow, or just gravitational acceleration without gas flow. Good
agreement between the differ- ent schemes for the time evolution of the size
spectra is observed in the presence of gravity or turbulence. Higher moments of
the size spectra are found to be a useful tool to characterize the growth of
the largest drops through collection. Remarkably, the tails of the size spectra
are reasonably well described by a gamma distribution in cases with gravity or
turbulence. The Lagrangian schemes are generally found to be superior over the
Eulerian one in terms of computational performance. However, it is shown that
the use of interpolation schemes such as the cloud-in-cell algorithm is
detrimental in connection with superparticle or superdroplet approaches.
Furthermore, the use of symmetric over asymmetric collection schemes is shown
to reduce the amount of scatter in the results.Comment: 36 pages, 17 figure
A minimal model for acoustic forces on Brownian particles
We present a generalization of the inertial coupling (IC) [Usabiaga et al. J.
Comp. Phys. 2013] which permits the resolution of radiation forces on small
particles with arbitrary acoustic contrast factor. The IC method is based on a
Eulerian-Lagrangian approach: particles move in continuum space while the fluid
equations are solved in a regular mesh (here we use the finite volume method).
Thermal fluctuations in the fluid stress, important below the micron scale, are
also taken into account following the Landau-Lifshitz fluid description. Each
particle is described by a minimal cost resolution which consists on a single
small kernel (bell-shaped function) concomitant to the particle. The main role
of the particle kernel is to interpolate fluid properties and spread particle
forces. Here, we extend the kernel functionality to allow for an arbitrary
particle compressibility. The particle-fluid force is obtained from an imposed
no-slip constraint which enforces similar particle and kernel fluid velocities.
This coupling is instantaneous and permits to capture the fast, non-linear
effects underlying the radiation forces on particles. Acoustic forces arise
either because an excess in particle compressibility (monopolar term) or in
mass (dipolar contribution) over the fluid values. Comparison with theoretical
expressions show that the present generalization of the IC method correctly
reproduces both contributions. Due to its low computational cost, the present
method allows for simulations with many particles using a standard Graphical
Processor Unit (GPU)
Variational approach to low-frequency kinetic-MHD in the current coupling scheme
Hybrid kinetic-MHD models describe the interaction of an MHD bulk fluid with
an ensemble of hot particles, which is described by a kinetic equation. When
the Vlasov description is adopted for the energetic particles, different
Vlasov-MHD models have been shown to lack an exact energy balance, which was
recently recovered by the introduction of non-inertial force terms in the
kinetic equation. These force terms arise from fundamental approaches based on
Hamiltonian and variational methods. In this work we apply Hamilton's
variational principle to formulate new current-coupling kinetic-MHD models in
the low-frequency approximation (i.e. large Larmor frequency limit). More
particularly, we formulate current-coupling hybrid schemes, in which energetic
particle dynamics are expressed in either guiding-center or gyrocenter
coordinates.Comment: v3.0. 30 page
- …