6,027 research outputs found
The remapped particle-mesh advection scheme
We describe the remapped particle-mesh method, a new mass-conserving method
for solving the density equation which is suitable for combining with
semi-Lagrangian methods for compressible flow applied to numerical weather
prediction. In addition to the conservation property, the remapped
particle-mesh method is computationally efficient and at least as accurate as
current semi-Lagrangian methods based on cubic interpolation. We provide
results of tests of the method in the plane, results from incorporating the
advection method into a semi-Lagrangian method for the rotating shallow-water
equations in planar geometry, and results from extending the method to the
surface of a sphere
Strong and auxiliary forms of the semi-Lagrangian method for incompressible flows
We present a review of the semi-Lagrangian method for advection-diusion and incompressible Navier-Stokes equations discretized with high-order methods. In particular, we compare the strong form where the departure points are computed directly via backwards integration with the auxiliary form where an auxiliary advection equation is solved instead; the latter is also referred to as Operator Integration Factor Splitting (OIFS) scheme. For intermediate size of time steps the auxiliary form is preferrable but for large time steps only the strong form is stable
A Moving Frame Algorithm for High Mach Number Hydrodynamics
We present a new approach to Eulerian computational fluid dynamics that is
designed to work at high Mach numbers encountered in astrophysical hydrodynamic
simulations. The Eulerian fluid conservation equations are solved in an
adaptive frame moving with the fluid where Mach numbers are minimized. The
moving frame approach uses a velocity decomposition technique to define local
kinetic variables while storing the bulk kinetic components in a smoothed
background velocity field that is associated with the grid velocity.
Gravitationally induced accelerations are added to the grid, thereby minimizing
the spurious heating problem encountered in cold gas flows. Separately tracking
local and bulk flow components allows thermodynamic variables to be accurately
calculated in both subsonic and supersonic regions. A main feature of the
algorithm, that is not possible in previous Eulerian implementations, is the
ability to resolve shocks and prevent spurious heating where both the preshock
and postshock Mach numbers are high. The hybrid algorithm combines the high
resolution shock capturing ability of the second-order accurate Eulerian TVD
scheme with a low-diffusion Lagrangian advection scheme. We have implemented a
cosmological code where the hydrodynamic evolution of the baryons is captured
using the moving frame algorithm while the gravitational evolution of the
collisionless dark matter is tracked using a particle-mesh N-body algorithm.
The MACH code is highly suited for simulating the evolution of the IGM where
accurate thermodynamic evolution is needed for studies of the Lyman alpha
forest, the Sunyaev-Zeldovich effect, and the X-ray background. Hydrodynamic
and cosmological tests are described and results presented. The current code is
fast, memory-friendly, and parallelized for shared-memory machines.Comment: 19 pages, 5 figure
Direct Numerical Simulation of Complex Multi-Fluid Flows Using a Combined Volume of Fluid and Immersed Boundary Method
In this paper a simulation model is presented for the Direct Numerical Simulation (DNS) of complex multi-fluid flows in which simultaneously (moving) deformable (drops or bubbles) and non-deformable (moving) elements (particles) are present, possibly with the additional presence of free surfaces. Our model combines the VOF model developed by van Sint Annaland et al. (2005) and the Immersed Boundary (IB) model developed by van der Hoef et al. (2006). The Volume of Fluid (VOF) part features i) an interface reconstruction technique based on piecewise linear interface representation ii) a three-dimensional version of the CSF model of Brackbill et al. (1992). The Immersed Boundary (IB) part incorporates both particle-fluid and particle-particle interaction via a Direct Forcing Method (DFM) and a hard sphere Discrete Particle (DP) approach. In our model a fixed (Eulerian) grid is utilized to solve the Navier-Stokes equations for the entire computational domain. The no-slip condition at the surface of the moving particles is enforced via a momentum source term which only acts in the vicinity of the particle surface. For the enforcement of the no-slip condition Lagrangian force points are used which are distributed evenly over the surface of the particle. Dissipative particle-particle and/or particle-wall collisions are accounted via a hard sphere DP approach (Hoomans et al., 1996) using a three-parameter particle-particle interaction model accounting for normal and tangential restitution and tangential friction. The capabilities of the hybrid VOF-IB model are demonstrated with a number of examples in which complex topological changes in the interface are encountered
Constrained hyperbolic divergence cleaning in smoothed particle magnetohydrodynamics with variable cleaning speeds
We present an updated constrained hyperbolic/parabolic divergence cleaning
algorithm for smoothed particle magnetohydrodynamics (SPMHD) that remains
conservative with wave cleaning speeds which vary in space and time. This is
accomplished by evolving the quantity instead of . Doing so
allows each particle to carry an individual wave cleaning speed, , that
can evolve in time without needing an explicit prescription for how it should
evolve, preventing circumstances which we demonstrate could lead to runaway
energy growth related to variable wave cleaning speeds. This modification
requires only a minor adjustment to the cleaning equations and is trivial to
adopt in existing codes. Finally, we demonstrate that our constrained
hyperbolic/parabolic divergence cleaning algorithm, run for a large number of
iterations, can reduce the divergence of the field to an arbitrarily small
value, achieving to machine precision.Comment: 23 pages, 16 figures, accepted for publication in Journal of
Computational Physic
- …