Publications 2015
Research papers, software, lecture notes

No. 2015.08
Higher order moment equations for rarefied gas mixtures
V. K. Gupta and M. Torrilhon
Subject: Kinetic theory, binary mixture, Grad’s moment method, 26-moment equations

Abstract

The fully nonlinear Grad’s N×26-moment (N×G26) equations for a mixture of N monatomic-inert-ideal gases made up of Maxwell molecules are derived. The boundary conditions for these equations are derived by using Maxwell’s accommodation model for each component in the mixture. The linear stability analysis is performed to show that the 2×G26 equations for a binary gas mixture of Maxwell molecules are linearly stable. The derived equations are used to study the heat flux problem for binary gas mixtures confined between parallel plates having different temperatures.

Reference

Proc. A. 471 (2015), no. 2173, 20140754, 20 pp.

No. 2015.07
Convergence study of moment approximations for boundary value problems of the Boltzmann-BGK equation
M. Torrilhon
Subject: Kinetic gas theory, Hermite expansion, channel flow

Abstract

The accuracy of moment equations as approximations of kinetic gas theory is studied for four different boundary value problems. The kinetic setting is given by the BGK equation linearized around a globally constant Maxwellian using one space dimension and a three-dimensional velocity space. The boundary value problems include Couette and Poiseuille flow as well as heat conduction between walls and heat conduction based on a locally varying heating source. The polynomial expansion of the distribution function allows for different moment theories of which two popular families are investigated in detail. Furthermore, optimal approximations for a given number of variables are studied empirically. The paper focuses on approximations with relatively low number of variables which allows to draw conclusions in particular about specific moment theories like the regularized 13-moment equations.

Reference

Commun. Comput. Phys. 18 (2015), no. 3, 529–557

No. 2015.06
Approximation of the linearized Boltzmann collision operator for hard-sphere and inverse-power-law model
Z. Cai and M. Torrilhon
Subject: Linearized collision model, hard sphere, inverse power law, series expansion

Abstract

A sequence of approximate linear collision models for hard-sphere and inverse-power-law gases is introduced. These models are obtained by expanding the linearized Boltzmann collision operator into series, and a practical algorithm is proposed for evaluating the coefficients in the series. The sequence is proven to be convergent to the linearized Boltzmann operator, and it established a connection between the Shakhov model and the linearized collision model. The convergence is demonstrated by solving the spatially homogeneous Boltzmann equation. By observing the magnitudes of the coefficients, simpler models are developed through removing small entries in the coefficient matrices.

Reference

J. Comput. Phys. 295 (2015), 617–643

No. 2015.05
Comparison of relaxation phenomena in binary gas-mixtures of Maxwell molecules and hard spheres
V. K. Gupta and M. Torrilhon
Subject: Kinetic theory, binary mixture, moment equations, relaxation, maxwell molecule, hard spheres

Abstract

The strategy for computing the Boltzmann collision integrals for gaseous mixtures is presented and bestowed to compute the fully non-linear Boltzmann collision integrals for hard sphere gas-mixtures. The Boltzmann collision integrals associated with the first 26 moments of each constituent in a gas-mixture are presented. Moreover, the Boltzmann collision integrals are exploited to study the relaxation phenomena of diffusion velocities, stresses and heat fluxes in binary gas-mixtures of Maxwell molecules and hard spheres.

Reference

Comput. Math. Appl. 70 (2015), no. 1, 73–88

No. 2015.04
On Singular Closures for the 5-Moment System in Kinetic Gas Theory
R. P. Schaerer and M. Torrilhon
Subject: Non-equilibrium gas dynamics, hyperbolic moment equations, kinetic theory, maximum entropy closures

Abstract

Moment equations provide a flexible framework for the approximation of the Boltzmann equation in kinetic gas theory. While moments up to second order are sufficient for the description of equilibrium processes, the inclusion of higher order moments, such as the heat flux vector, extends the validity of the Euler equations to non-equilibrium gas flows in a natural way.

Unfortunately, the classical closure theory proposed by Grad leads to moment equations, which suffer not only from a restricted hyperbolicity region but are also affected by non-physical sub-shocks in the continuous shock-structure problem if the shock velocity exceeds a critical value. Amore recently suggested closure theory based on the maximum entropy principle yields symmetric hyperbolic moment equations. However, if moments higher than second order are included, the computational demand of this closure can be overwhelming. Additionally, it was shown for the 5-moment system that the closing flux becomes singular on a subset of moments including the equilibrium state.

Motivated by recent promising results of closed-form, singular closures based on the maximum entropy approach, we study regularized singular closures that become singular on a subset of moments when the regularizing terms are removed. In order to study some implications of singular closures, we use a recently proposed explicit closure for the 5-moment equations. We show that this closure theory results in a hyperbolic system that can mitigate the problem of sub-shocks independent of the shock wave velocity and handle strongly non-equilibrium gas flows.

Reference

Commun. Comput. Phys. 17 (2015), no. 2, 371–400

No. 2015.03
Deterministic model for the transport of energetic particles: Application in the electron radiotherapy
J. Caron, J.-L. Feugeas, B. Dubroca, G. Kantor, C. Dejean, G. Birindelli, T. Pichard, P. Nicolaï, E. d’Humières, M. Frank, and V. Tikhonchuk
Subject: Radiotherapy treatment planning system, deterministic algorithm, electron dosimetry, Monte Carlo algorithm

Abstract

A new deterministic method for calculating the dose distribution in the electron radiotherapy field is presented. The aim of this work was to validate our model by comparing it with the Monte Carlo simulation toolkit, GEANT4. A comparison of the longitudinal and transverse dose deposition profiles and electron distributions in homogeneous water phantoms showed a good accuracy of our model for electron transport, while reducing the calculation time by a factor of 50. Although the Bremsstrahlung effect is not yet implemented in our model, we propose here a method that solves the Boltzmann kinetic equation and provides a viable and efficient alternative to the expensive Monte Carlo modeling.

Reference

Physica Medica 31 (2015) 912–921

No. 2015.02
The nonclassical Boltzmann equation and diffusion-based approximations to the Boltzmann equation
M. Frank, K. Krycki, E. W. Larsen, and R. Vasques
Subject: Boltzmann equation, diffusion approximation, nonclassical transport, Monte Carlo

Abstract

We show that several diffusion-based approximations (classical diffusion or $SP_1$, $SP_2$, $SP_3$) to the linear Boltzmann equation can (for an infinite, homogeneous medium) be represented exactly by a nonclassical transport equation. As a consequence, we indicate a method to solve these diffusion-based approximations to the Boltzmann equation via Monte Carlo methods, with only statistical errors—no truncation errors.

Reference

SIAM J. Appl. Math. 75 (2015), no. 3, 1329–1345

Download

arXiv:1412.1000

No. 2015.01
Flux splitting for stiff equations: a notion on stability
J. Schütz, S. Noelle
Subject: IMEX finite volume, asymptotic preserving, flux splitting, modified equation, stability analysis

Abstract

For low Mach number flows, there is a strong recent interest in the development and analysis of IMEX (implicit/explicit) schemes, which rely on a splitting of the convective flux into stiff and nonstiff parts. A key ingredient of the analysis is the so-called Asymptotic Preserving property, which guarantees uniform consistency and stability as the Mach number goes to zero. While many authors have focused on asymptotic consistency, we study asymptotic stability in this paper: does an IMEX scheme allow for a CFL number which is independent of the Mach number? We derive a stability criterion for a general linear hyperbolic system. In the decisive eigenvalue analysis, the advective term, the upwind diffusion and a quadratic term stemming from the truncation in time all interact in a subtle way. As an application, we show that a new class of splittings based on characteristic decomposition, for which the commutator vanishes, avoids the deterioration of the time step which has sometimes been observed in the literature.

Reference

J. Sci. Comput. 64 (2015), no. 2, 522–540

Download

arXiv:1412.1595