9,269 research outputs found
A Nonlinear Multigrid Steady-State Solver for Microflow
We develop a nonlinear multigrid method to solve the steady state of
microflow, which is modeled by the high order moment system derived recently
for the steady-state Boltzmann equation with ES-BGK collision term. The solver
adopts a symmetric Gauss-Seidel iterative scheme nested by a local Newton
iteration on grid cell level as its smoother. Numerical examples show that the
solver is insensitive to the parameters in the implementation thus is quite
robust. It is demonstrated that expected efficiency improvement is achieved by
the proposed method in comparison with the direct time-stepping scheme
Kinetic Solvers with Adaptive Mesh in Phase Space
An Adaptive Mesh in Phase Space (AMPS) methodology has been developed for
solving multi-dimensional kinetic equations by the discrete velocity method. A
Cartesian mesh for both configuration (r) and velocity (v) spaces is produced
using a tree of trees data structure. The mesh in r-space is automatically
generated around embedded boundaries and dynamically adapted to local solution
properties. The mesh in v-space is created on-the-fly for each cell in r-space.
Mappings between neighboring v-space trees implemented for the advection
operator in configuration space. We have developed new algorithms for solving
the full Boltzmann and linear Boltzmann equations with AMPS. Several recent
innovations were used to calculate the discrete Boltzmann collision integral
with dynamically adaptive mesh in velocity space: importance sampling,
multi-point projection method, and the variance reduction method. We have
developed an efficient algorithm for calculating the linear Boltzmann collision
integral for elastic and inelastic collisions in a Lorentz gas. New AMPS
technique has been demonstrated for simulations of hypersonic rarefied gas
flows, ion and electron kinetics in weakly ionized plasma, radiation and light
particle transport through thin films, and electron streaming in
semiconductors. We have shown that AMPS allows minimizing the number of cells
in phase space to reduce computational cost and memory usage for solving
challenging kinetic problems
Solving the Boltzmann equation deterministically by the fast spectral method : application to gas microflows
Based on the fast spectral approximation to the Boltzmann collision operator, we present an accurate and efficient deterministic numerical method for solving the Boltzmann equation. First, the linearised Boltzmann equation is solved for Poiseuille and thermal creep flows, where the influence of different molecular models on the mass and heat flow rates is assessed, and the Onsager-Casimir relation at the microscopic level for large Knudsen numbers is demonstrated. Recent experimental measurements of mass flow rates along a rectangular tube with large aspect ratio are compared with numerical results for the linearised Boltzmann equation. Then, a number of two-dimensional micro flows in the transition and free molecular flow regimes are simulated using the nonlinear Boltzmann equation. The influence of the molecular model is discussed, as well as the applicability of the linearised Boltzmann equation. For thermally driven flows in the free molecular regime, it is found that the magnitudes of the flow velocity are inversely proportional to the Knudsen number. The streamline patterns of thermal creep flow inside a closed rectangular channel are analysed in detail: when the Knudsen number is smaller than a critical value, the flow pattern can be predicted based on a linear superposition of the velocity profiles of linearised Poiseuille and thermal creep flows between parallel plates. For large Knudsen numbers, the flow pattern can be determined using the linearised Poiseuille and thermal creep velocity profiles at the critical Knudsen number. The critical Knudsen number is found to be related to the aspect ratio of the rectangular channel
Three-dimensional Boltzmann-Hydro code for core-collapse in massive stars I. special relativistic treatments
We propose a novel numerical method for solving multi-dimensional, special
relativistic Boltzmann equations for neutrinos coupled to hydrodynamics
equations. It is meant to be applied to simulations of core-collapse
supernovae. We handle special relativity in a non-conventional way, taking
account of all orders of v/c. Consistent treatment of advection and collision
terms in the Boltzmann equations is the source of difficulties, which we
overcome by employing two different energy grids: Lagrangian remapped and
laboratory fixed grids. We conduct a series of basic tests and perform a
one-dimensional simulation of core-collapse, bounce and shock-stall for a
15M_{sun} progenitor model with a minimum but essential set of microphysics. We
demonstrate in the latter simulation that our new code is capable of handling
all phases in core-collapse supernova. For comparison, a non-relativistic
simulation is also conducted with the same code, and we show that they produce
qualitatively wrong results in neutrino transfer. Finally, we discuss a
possible incorporation of general relativistic effects in our method.Comment: 25 pages, 22 figures, submitted to Ap
Numerical Regularized Moment Method of Arbitrary Order for Boltzmann-BGK Equation
We introduce a numerical method for solving Grad's moment equations or
regularized moment equations for arbitrary order of moments. In our algorithm,
we do not need explicitly the moment equations. As an instead, we directly
start from the Boltzmann equation and perform Grad's moment method \cite{Grad}
and the regularization technique \cite{Struchtrup2003} numerically. We define a
conservative projection operator and propose a fast implementation which makes
it convenient to add up two distributions and provides more efficient flux
calculations compared with the classic method using explicit expressions of
flux functions. For the collision term, the BGK model is adopted so that the
production step can be done trivially based on the Hermite expansion. Extensive
numerical examples for one- and two-dimensional problems are presented.
Convergence in moments can be validated by the numerical results for different
number of moments.Comment: 33 pages, 13 figure
A Unified Gas-kinetic Scheme for Continuum and Rarefied Flows IV: full Boltzmann and Model Equations
Fluid dynamic equations are valid in their respective modeling scales. With a
variation of the modeling scales, theoretically there should have a continuous
spectrum of fluid dynamic equations. In order to study multiscale flow
evolution efficiently, the dynamics in the computational fluid has to be
changed with the scales. A direct modeling of flow physics with a changeable
scale may become an appropriate approach. The unified gas-kinetic scheme (UGKS)
is a direct modeling method in the mesh size scale, and its underlying flow
physics depends on the resolution of the cell size relative to the particle
mean free path. The cell size of UGKS is not limited by the particle mean free
path. With the variation of the ratio between the numerical cell size and local
particle mean free path, the UGKS recovers the flow dynamics from the particle
transport and collision in the kinetic scale to the wave propagation in the
hydrodynamic scale.
The previous UGKS is mostly constructed from the evolution solution of
kinetic model equations. This work is about the further development of the UGKS
with the implementation of the full Boltzmann collision term in the region
where it is needed. The central ingredient of the UGKS is the coupled treatment
of particle transport and collision in the flux evaluation across a cell
interface, where a continuous flow dynamics from kinetic to hydrodynamic scales
is modeled. The newly developed UGKS has the asymptotic preserving (AP)
property of recovering the NS solutions in the continuum flow regime, and the
full Boltzmann solution in the rarefied regime. In the mostly unexplored
transition regime, the UGKS itself provides a valuable tool for the flow study
in this regime. The mathematical properties of the scheme, such as stability,
accuracy, and the asymptotic preserving, will be analyzed in this paper as
well
- …