Search

In this paper, we propose a modification of an acoustic-transport operator splitting Lagrange-projection method for simulating compressible flows with gravity. The original method involves two steps that respectively account for acoustic and transport effects. Our work proposes a simple modification of the transport step, and the resulting modified scheme turns out to be a flux-splitting method. This new numerical method is less computationally expensive in the low-Mach regime, more memory efficient, and easier to implement than the original one. We prove stability properties for this new scheme by showing that under classical CFL conditions, the method is positivity preserving for mass, energy and entropy satisfying. The flexible flux-splitting structure of the method enables straightforward extensions of the method to multi-dimensional problems (with respect to space) and high-order discretizations that are presented in this work. We also propose an interpretation of the flux-splitting solver as a relaxation approximation. Both the stability and the accuracy of the new method are tested against one-dimensional and two-dimensional numerical experiments that involve highly compressible flows and low-Mach regimes. • The all-regime acoustic/transport splitting solver proposed in [16] can be modified into a flux-splitting, one-step version. • The resulting method is less computationally expensive, more memory efficient, and easier to implement than the original one. • Stability properties are derived by interpreting the scheme as a convex combination of the advection/pressure steps. • This recasting allows an easy high-order extension as the new method is a standard flux-based solver. • Following the lines of [32] , we equip our scheme with the well-balanced property for gravity source terms. [ABSTRACT FROM AUTHOR]

In this paper we deal with the extension of the Fast Kinetic Scheme (FKS) (Dimarco and Loubère, 2013 [26] ) originally constructed for solving the BGK equation, to the more challenging case of the Boltzmann equation. The scheme combines a robust and fast method for treating the transport part based on an innovative Lagrangian technique supplemented with conservative fast spectral schemes to treat the collisional operator by means of an operator splitting approach. This approach along with several implementation features related to the parallelization of the algorithm permits to construct an efficient simulation tool which is numerically tested against exact and reference solutions on classical problems arising in rarefied gas dynamic. We present results up to the 3 D × 3 D case for unsteady flows for the Variable Hard Sphere model which may serve as benchmark for future comparisons between different numerical methods for solving the multidimensional Boltzmann equation. For this reason, we also provide for each problem studied details on the computational cost and memory consumption as well as comparisons with the BGK model or the limit model of compressible Euler equations. [ABSTRACT FROM AUTHOR]

Abstract: Several plasma devices, including Hall and ion thrusters, operate by ionizing a low density neutral gas for which the mean free path between collisions of gas molecules is greater than typical device dimensions. In general, the discrete-particle algorithms used to calculate the neutral gas ignore velocity changes due to collisions between gas molecules. However, particle algorithms are a source of unphysical statistical noise that may detract from the study of the plasma physics, the prime purpose of most simulations. In this paper we present a new neutral gas algorithm for use in plasma simulation codes that exploits the fact that very few collisions change the velocity of neutral gas molecules. The algorithm assumes that the particle velocity distribution function for neutrals emitted from a given surface remains unchanged except for a scale factor that reflects the loss of neutrals to ionization. The sources of neutrals may be gas inlets, and isotropic, thermally accommodated, gas molecules coming off chamber surfaces including recombined ions. The algorithm is implemented in two dimensions (R–Z) with emitting surfaces represented as surfaces of revolution. The advantage of this algorithm over the conventional particle approach is the absence of statistical noise. [ABSTRACT FROM AUTHOR]

Abstract: We present a new time-accurate algorithm for the explicit numerical integration of the compressible Euler equations of gas dynamics. This technique is based on the discrete-event simulation (DES) methodology for nonlinear flux-conservative PDEs [Y.A. Omelchenko, H. Karimabadi, Self-adaptive time integration of flux-conservative equations with sources, J. Comput. Phys. 216 (1) (2006) 179–194]. DES enables adaptive distribution of CPU resources in accordance with local time scales of the underlying numerical solution. It distinctly stands apart from multiple (local) time-stepping algorithms in that it requires neither selecting a global synchronization time step nor pre-determining a sequence of time-integration operations for individual parts of a heterogeneous numerical system. In this paper we extend the DES methodology in three important directions: (i) we apply DES to a system of coupled gas dynamics equations discretized via a central-upwind scheme [A. Kurganov, E. Tadmor, New high-resolution central schemes for nonlinear conservation laws and convection–diffusion equations, J. Comput. Phys. 160 (2000) 241–282; A. Kurganov, S. Noelle, G. Petrova, Semidiscrete central-upwind schemes for hyperbolic conservation laws and Hamilton–Jacobi equations, SIAM J. Sci. Comput. 23 (3) (2001) 707–740]; (ii) we introduce a new Preemptive Event Processing (PEP) technique, which automatically enforces synchronous execution of events with sufficiently close update times; (iii) we significantly improve the accuracy of the previous algorithm [Y.A. Omelchenko, H. Karimabadi, Self-adaptive time integration of flux-conservative equations with sources, J. Comput. Phys. 216 (1) (2006) 179–194] by applying locally second-order-in-time flux-conserving corrections to the solution obtained with the forward Euler scheme. The performance of the new technique is demonstrated in a series of one-dimensional gas dynamics test problems by comparing numerical solutions obtained in event-driven and equivalent time-stepping simulations. [Copyright &y& Elsevier]

The present paper highlights the importance of management of the discontinuous reconstruction in the kinetic schemes for gasdynamic equation systems. Firstly, it is revealed by the analysis of the gas kinetic-BGK scheme [JCP 171 (2001) 289] that a continuous reconstruction created from a discontinuous one is a key to the successful kinetic schemes. When it is applied to a well-resolved region, the numerical flux that takes account of the collision effect correctly becomes Lax–Wendroff-like. When it is applied to an unresolved region, such as a shock layer, an appreciable numerical dissipation, which contributes to the suppression of spurious oscillations, is produced. Secondly, new kinetic schemes for the compressible Navier–Stokes (Euler) equations are developed by using the key. The numerical flux of one of the schemes is computed by using the splitting algorithm, where the effect of the molecular collision is directly taken into account and the undesirable error of the splitting algorithm in the case where the time step is much larger than the mean free time is avoided by a simple modification of the initial data. Although a discontinuous reconstruction is employed in the approximation of the initial data, the continuity is automatically taken into account in the dominant part of the numerical flux. The other schemes are the extensions of the Lax–Wendroff-type scheme to the case of the key reconstruction. Thirdly, the performance of these schemes is tested. It is demonstrated that they work as shock capturing schemes and yield fine boundary-layer profiles with a reasonable number of cells, such as 10 cells in the layer. [Copyright &y& Elsevier]