A high order cell-centered semi-Lagrangian scheme for multi-dimensional kinetic simulations of neutral gas flows

DOI10.1016/J.JCP.2012.01.008MaRDI QIDQ417861zbMATH OpenOpenAlexFDO

Authors Y. Güçlü, W. N. G. Hitchon

Publication date 14 May 2012

Published in Journal of Computational Physics (Search for Journal in Brave)

Full work available at URL https://arxiv.org/abs/1305.5477

zbMATH Keywords

Boltzmann equationadvectionconvected schemeneutral gas kinetics

Mathematics Subject Classification ID

Particle methods and lattice-gas methods (76M28) Rarefied gas flows, Boltzmann equation in fluid mechanics (76P05)

Abstract: The term `Convected Scheme' (CS) refers to a family of algorithms, most usually applied to the solution of Boltzmann's equation, which uses a method of characteristics in an integral form to project an initial cell forward to a group of final cells. As such the CS is a `forward-trajectory' semi-Lagrangian scheme. For multi-dimensional simulations of neutral gas flows, the cell-centered version of this semi-Lagrangian (CCSL) scheme has advantages over other options due to its implementation simplicity, low memory requirements, and easier treatment of boundary conditions. The main drawback of the CCSL-CS to date has been its high numerical diffusion in physical space, because of the 2${}^{extnd}$ order remapping that takes place at the end of each time step. By means of a Modified Equation Analysis, it is shown that a high order estimate of the remapping error can be obtained a priori, and a small correction to the final position of the cells can be applied upon remapping, in order to achieve full compensation of this error. The resulting scheme is 4${}^{extth}$ order accurate in space while retaining the desirable properties of the CS: it is conservative and positivity-preserving, and the overall algorithm complexity is not appreciably increased. Two monotone (i.e. non-oscillating) versions of the fourth order CCSL-CS are also presented: one uses a common flux-limiter approach; the other uses a non-polynomial reconstruction to evaluate the derivatives of the density function. The method is illustrated in simple one- and two-dimensional examples, and a fully 3D solution of the Boltzmann equation describing expansion of a gas into vacuum through a cylindrical tube.

Recommendations

• Arbitrarily high order convected scheme solution of the Vlasov-Poisson system
• A truly forward semi-Lagrangian WENO scheme for the Vlasov-Poisson system
• scientific article; zbMATH DE number 6453528
• High order conservative semi-Lagrangian scheme for the BGK model of the Boltzmann equation
• A class of low dissipative schemes for solving kinetic equations

Cites work

• scientific article; zbMATH DE number 4014334 (Why is no real title available?)
• scientific article; zbMATH DE number 914383 (Why is no real title available?)
• A Model for Collision Processes in Gases. I. Small Amplitude Processes in Charged and Neutral One-Component Systems
• Accelerated solution of the Boltzmann equation
• Antidiffusive Velocities for Multipass Donor Cell Advection
• Convective scheme solution of the Boltzmann transport equation for nanoscale semiconductor devices
• Convergence Analysis for a Class of High-Order Semi-Lagrangian Advection Schemes
• Direct simulation of gas flows at the molecular level
• Flux-corrected transport. I: SHASTA, a fluid transport algorithm that works
• High Resolution Schemes Using Flux Limiters for Hyperbolic Conservation Laws
• High resolution schemes for hyperbolic conservation laws
• Kinetic description of flow past a micro-plate
• Limiter-Free Third Order Logarithmic Reconstruction
• Local Piecewise Hyperbolic Reconstruction of Numerical Fluxes for Nonlinear Scalar Conservation Laws
• On sub-linear convergence for linearly degenerate waves in capturing schemes
• Study of rarefied shear flow by the discrete velocity method
• The modified equation approach to the stability and accuracy analysis of finite-difference methods
• The multidimensional positive definite advection transport algorithm: Nonoscillatory option
• Towards the ultimate conservative difference scheme. II: Monotonicity and conservation combined in a second-order scheme

Cited in

(12)

• A class of low dissipative schemes for solving kinetic equations
• Fast kinetic scheme: efficient MPI parallelization strategy for 3D Boltzmann equation
• Implicit discontinuous Galerkin method for the Boltzmann equation
• Towards an ultra efficient kinetic scheme. II: The high order case
• A performance comparison of semi-Lagrangian discontinuous Galerkin and spline based Vlasov solvers in four dimensions
• Quantum algorithm for the collisionless Boltzmann equation
• A spectral-Lagrangian Boltzmann solver for a multi-energy level gas
• Arbitrarily high order convected scheme solution of the Vlasov-Poisson system
• Noiseless Vlasov-Poisson simulations with linearly transformed particles
• Towards an ultra efficient kinetic scheme. part III: high-performance-computing
• A truly forward semi-Lagrangian WENO scheme for the Vlasov-Poisson system
• An efficient numerical method for solving the Boltzmann equation in multidimensions

Describes a project that uses

Uses Software

• SHASTA

This page was built for publication: A high order cell-centered semi-Lagrangian scheme for multi-dimensional kinetic simulations of neutral gas flows

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q417861)