2,465 research outputs found
Inertial Coupling Method for particles in an incompressible fluctuating fluid
We develop an inertial coupling method for modeling the dynamics of
point-like 'blob' particles immersed in an incompressible fluid, generalizing
previous work for compressible fluids. The coupling consistently includes
excess (positive or negative) inertia of the particles relative to the
displaced fluid, and accounts for thermal fluctuations in the fluid momentum
equation. The coupling between the fluid and the blob is based on a no-slip
constraint equating the particle velocity with the local average of the fluid
velocity, and conserves momentum and energy. We demonstrate that the
formulation obeys a fluctuation-dissipation balance, owing to the
non-dissipative nature of the no-slip coupling. We develop a spatio-temporal
discretization that preserves, as best as possible, these properties of the
continuum formulation. In the spatial discretization, the local averaging and
spreading operations are accomplished using compact kernels commonly used in
immersed boundary methods. We find that the special properties of these kernels
make the discrete blob a particle with surprisingly physically-consistent
volume, mass, and hydrodynamic properties. We develop a second-order
semi-implicit temporal integrator that maintains discrete
fluctuation-dissipation balance, and is not limited in stability by viscosity.
Furthermore, the temporal scheme requires only constant-coefficient Poisson and
Helmholtz linear solvers, enabling a very efficient and simple FFT-based
implementation on GPUs. We numerically investigate the performance of the
method on several standard test problems...Comment: Contains a number of corrections and an additional Figure 7 (and
associated discussion) relative to published versio
A comparison of Eulerian and Lagrangian transport and non-linear reaction algorithms
When laboratory-measured chemical reaction rates are used in simulations at the field-scale, the models typically overpredict the apparent reaction rates. The discrepancy is primarily due to poorer mixing of chemically distinct waters at the larger scale. As a result, realistic field-scale predictions require accurate simulation of the degree of mixing between fluids. The Lagrangian particle-tracking (PT) method is a now-standard way to simulate the transport of conservative or sorbing solutes. The method’s main advantage is the absence of numerical dispersion (and its artificial mixing) when simulating advection. New algorithms allow particles of different species to interact in nonlinear (e.g., bimolecular) reactions. Therefore, the PT methods hold a promise of more accurate field-scale simulation of reactive transport because they eliminate the masking effects of spurious mixing due to advection errors inherent in grid-based methods. A hypothetical field-scale reaction scenario is constructed and run in PT and Eulerian (finite-volume/finite-difference) simulators. Grid-based advection schemes considered here include 1st- to 3rd-order spatially accurate total-variation-diminishing flux-limiting schemes, both of which are widely used in current transport/reaction codes. A homogeneous velocity field in which the Courant number is everywhere unity, so that the chosen Eulerian methods incur no error when simulating advection, shows that both the Eulerian and PT methods can achieve convergence in the L1 (integrated concentration) norm, but neither shows stricter pointwise convergence. In this specific case with a constant dispersion coefficient and bimolecular reaction A+B¿P, the correct total amount of product is 0.221MA0, where MA0 is the original mass of reactant A. When the Courant number drops, the grid-based simulations can show remarkable errors due to spurious over- and under-mixing. In a heterogeneous velocity field (keeping the same constant and isotropic dispersion), the PT simulations show an increased reaction total from 0.221MA0 to 0.372MA0 due to fluid deformation, while the 1st-order Eulerian simulations using ˜ 106 cells (with a classical grid Peclet number ¿x/aL of 10) have total product of 0.53MA0, or approximately twice as much additional reaction due to advection error. The 3rd-order TVD algorithm fares better, with total product of 0.394MA0, or about 1.14 times the increased reaction total. A very strict requirement on grid Peclet numbers for Eulerian simulations will be required for realistic reactions because of their nonlinear nature. We analytically estimate the magnitude of the effect for the end-member cases of very fast and very slow reactions and show that in either case, the mass produced is proportional to View the MathML source where Pe is the Peclet number. Therefore, extra mass is produced according to View the MathML source where the dispersion includes any numerical dispersion error. We test two PT methods, one that kills particles upon reaction and another that decrements a particle’s mass. For the bimolecular reaction studied here, the computational demands of the particle-killing methods are much smaller than, and the particle-number-preserving algorithm are on par with, the fastest Eulerian methods.Peer ReviewedPostprint (author's final draft
A Hybrid Godunov Method for Radiation Hydrodynamics
From a mathematical perspective, radiation hydrodynamics can be thought of as
a system of hyperbolic balance laws with dual multiscale behavior (multiscale
behavior associated with the hyperbolic wave speeds as well as multiscale
behavior associated with source term relaxation). With this outlook in mind,
this paper presents a hybrid Godunov method for one-dimensional radiation
hydrodynamics that is uniformly well behaved from the photon free streaming
(hyperbolic) limit through the weak equilibrium diffusion (parabolic) limit and
to the strong equilibrium diffusion (hyperbolic) limit. Moreover, one finds
that the technique preserves certain asymptotic limits. The method incorporates
a backward Euler upwinding scheme for the radiation energy density and flux as
well as a modified Godunov scheme for the material density, momentum density,
and energy density. The backward Euler upwinding scheme is first-order accurate
and uses an implicit HLLE flux function to temporally advance the radiation
components according to the material flow scale. The modified Godunov scheme is
second-order accurate and directly couples stiff source term effects to the
hyperbolic structure of the system of balance laws. This Godunov technique is
composed of a predictor step that is based on Duhamel's principle and a
corrector step that is based on Picard iteration. The Godunov scheme is
explicit on the material flow scale but is unsplit and fully couples matter and
radiation without invoking a diffusion-type approximation for radiation
hydrodynamics. This technique derives from earlier work by Miniati & Colella
2007. Numerical tests demonstrate that the method is stable, robust, and
accurate across various parameter regimes.Comment: accepted for publication in Journal of Computational Physics; 61
pages, 15 figures, 11 table
A mass-conserving sparse grid combination technique with biorthogonal hierarchical basis functions for kinetic simulations
The exact numerical simulation of plasma turbulence is one of the assets and
challenges in fusion research. For grid-based solvers, sufficiently fine
resolutions are often unattainable due to the curse of dimensionality. The
sparse grid combination technique provides the means to alleviate the curse of
dimensionality for kinetic simulations. However, the hierarchical
representation for the combination step with the state-of-the-art hat functions
suffers from poor conservation properties and numerical instability.
The present work introduces two new variants of hierarchical multiscale basis
functions for use with the combination technique: the biorthogonal and full
weighting bases. The new basis functions conserve the total mass and are shown
to significantly increase accuracy for a finite-volume solution of constant
advection. Further numerical experiments based on the combination technique
applied to a semi-Lagrangian Vlasov--Poisson solver show a stabilizing effect
of the new bases on the simulations
- …