A new closure for Grad's 13 moment equations is presented that adds terms of Super-Burnett order to the balances of pressure deviator and heat flux vector. The additional terms are derived from equations for higher moments by means of the distribution function for 13 moments. The resulting system of equations contains the Burnett and Super-Burnett equations when expanded in a series in the Knudsen number. However, other than the Burnett and Super-Burnett equations, the new set of equations is linearly stable for all wavelengths and frequencies. Dispersion relation and damping for the new equations agree better with experimental data than those for the Navier-Stokes-Fourier equations, or the original 13 moments system. The new equations also allow the description of Knudsen boundary layers.

Recently a new system of field equations for the accurate description of flows in rarefied gases, called regularized 13-moment equations, was obtained by means of a hybrid gas kinetic approach. The first part of this paper discusses the relationship of the new system to classical high-order theories like the Burnett and super-Burnett equations as well as to modified models like the augmented and regularized Burnett equations. In the second part, shock structure calculations with the new theory are presented and compared to direct-simulation Monte Carlo (DSMC) solutions and to solutions of the Burnett models. Owing to additional higher-order dissipation in the system, the profiles are smooth for any Mach number, in contrast to the results of Grad's 13-moment case. The results show reliable quantitative agreement with DSMC simulations for Mach numbers up to M 0 ≈ 3.0. The agreement is better for Maxwell molecules than for hard spheres. The results of the augmented Burnett equations are comparable, but these equations are shown to be spatially unstable. Additionally, a validiation procedure for the new equations is presented by investigating the positivity of Grad's distribution function.

This article discusses the development of continuum models to describe processes in gases in which the particle collisions cannot maintain thermal equilibrium. Such a situation typically is present in rarefied or diluted gases, for flows in microscopic settings, or in general whenever the Knudsen number-the ratio between the mean free path of the particles and a macroscopic length scale-becomes significant. The continuum models are based on the stochastic description of the gas by Boltzmann's equation in kinetic gas theory. With moment approximations, extended fluid dynamic equations can be derived, such as the regularized 13-moment equations. Moment equations are introduced in detail, and typical results are reviewed for channel flow, cavity flow, and flow past a sphere in the low-Mach number setting for which both evolution equations and boundary conditions are well established. Conversely, nonlinear, high-speed processes require special closures that are still under development. Current approaches are examined, along with the challenge of computing shock wave profiles based on continuum equations. 429

No abstract

Boundary conditions are the major obstacle in simulations based on advanced continuum models of rarefied and micro-flows of gases. In this paper we present a theory how to combine the regularized 13-moment-equations derived from Boltzmann's equation with boundary conditions obtained from Maxwell's kinetic accommodation model. While for the linear case these kinetic boundary conditions suffice, we need additional conditions in the non-linear case. They are provided by the bulk solutions obtained after properly transforming the equations while keeping their asymptotic accuracy with respect to Boltzmann's equation.After finding a suitable set of boundary conditions and equations, a numerical method for generic shear flow problems is formulated. Several test simulations demonstrate the stable and oscillation-free performance of the new approach.

No abstract

Shallow flow models are used for a large number of applications including weather forecasting, open channel hydraulics and simulation-based natural hazard assessment. In these applications the shallowness of the process motivates depth-averaging. While the shallow flow formulation is advantageous in terms of computational efficiency, it also comes at the price of losing vertical information such as the flow's velocity profile. This gives rise to a model error, which limits the shallow flow model's predictive power and is often not explicitly quantifiable. We propose the use of vertical moments to overcome this problem. The shallow moment approximation preserves information on the vertical flow structure while still making use of the simplifying framework of depth-averaging. In this article, we derive a generic shallow flow moment system of arbitrary order starting from a set of balance laws, which has been reduced by scaling arguments. The derivation is based on a fully vertically resolved reference model with the vertical coordinate mapped onto the unit interval. We specify the shallow flow moment hierarchy for kinematic and Newtonian flow conditions and present 1D numerical results for shallow moment systems up to third order. Finally, we assess their performance with respect to both the standard shallow flow equations as well as with respect to the vertically resolved reference model. Our results show that depending on the parameter regime, e.g. friction and slip, shallow moment approximations significantly reduce the model error in shallow flow regimes and have a lot of potential to increase the predictive power of shallow flow models, while keeping them computationally cost efficient.

By a further study of the mechanism of the hyperbolic regularization of the moment system for Boltzmann equation proposed in [Z. Cai, Y. Fan, R. Li, Comm. Math. Sci. 11(2): 547-571, 2013], we point out that the key point is treating the time and space derivative in the same way. Based on this understanding, a uniform framework to derive globally hyperbolic moment systems from kinetic equations using an operator projection method is proposed. The framework is so concise and clear that it can be treated as an algorithm with four inputs to derive hyperbolic moment system by routine calculations. Almost all existing globally hyperbolic moment system can be included in the framework, as well as some new moment system including globally hyperbolic regularized versions of Grad ordered moment system and a multidimensional extension of the quadrature-based moment system.Comment: 32 pages, 2 figure