Convergence proof of the DSMC method and the Gas-Kinetic Unified Algorithm for the Boltzmann equation

This paper investigates the convergence proof of the Direct Simulation Monte Carlo(DSMC) method and the Gas-Kinetic Unified Algorithm in simulating the Boltzmann equation.It can be shown that the particle velocity distribution function obtained by the DSMC method converges to a modified form of the Boltzmann equation,which is the equation of the gas-kinetic unified algorithm to directly solve the molecular velocity distribution function.Their convergence is derived through mathematical treatment.The collision frequency is presented using various molecular models and the local equilibrium distribution function is obtained by Enskog expansion using the converged equation of the DSMC method.These two expressions agree with those used in the unified algorithm.Numerical validation of the converging consistency between these two approaches is illustrated by simulating the pressure driven Poiseuille flow in the slip transition flow regime and the two-dimensional and three-dimensional flows around a circular cylinder and spherical-cone reentry body covering the whole flow regimes from low speed micro-channel flow to high speed non-equilibrium aerothermodynamics.

Gas-kinetic unified algorithm for computable modeling of Boltzmann equation and application to aerothermodynamics for falling disintegration of uncontrolled Tiangong-No.1 spacecraft

How to solve the hypersonic aerothermodynamics around large-scale uncontrolled spacecraft during falling disintegrated process from outer space to earth,is the key to resolve the problems of the uncontrolled Tiangong-No.1 spacecraft reentry crash.To study aerodynamics of spacecraft reentry covering various flow regimes,a Gas-Kinetic Unified Algorithm(GKUA)has been presented by computable modeling of the collision integral of the Boltzmann equation over tens of years.On this basis,the rotational and vibrational energy modes are considered as the independent variables of the gas molecular velocity distribution function,a kind of Boltzmann model equation involving in internal energy excitation is presented by decomposing the collision term of the Boltzmann equation into elastic and inelastic collision terms.Then,the gas-kinetic numerical scheme is constructed to capture the time evolution of the discretized velocity distribution functions by developing the discrete velocity ordinate method and numerical quadrature technique.The unified algorithm of the Boltzmann model equation involving thermodynamics non-equilibrium effect is presented for the whole range of flow regimes.The gas-kinetic massive parallel computing strategy is developed to solve the hypersonic aerothermodynamics with the processor cores 500~45,000 at least 80%parallel efficiency.To validate the accuracy of the GKUA,the hypersonic flows are simulated including the reentry Tiangong-1 spacecraft shape with the wide range of Knudsen numbers of 220~0.00005 by the comparison of the related results from the DSMC and N-S coupled methods,and the low-density tunnel experiment etc.For uncontrolling spacecraft falling problem,the finite-element algorithm for dynamic thermalforce coupling response is presented,and the unified simulation of the thermal structural response and the hypersonic flow field is tested on the Tiangong-1 shape under reentry aerodynamic environment.Then,the forecasting analysis platform of end-of-life largescale spacecraft flying track is established on the basis of ballistic computation combined with reentry aerothermodynamics and deformation failure/disintegration.

Study on the Splitting Methods for Separable Convex Optimization in a Unified Algorithmic Framework

It is well recognized the convenience of converting the linearly constrained convex optimization problems to a monotone variational inequality.Recently,we have proposed a unified algorithmic framework which can guide us to construct the solution methods for solving these monotone variational inequalities.In this work,we revisit two full Jacobian decomposition of the augmented Lagrangian methods for separable convex programming which we have studied a few years ago.In particular,exploiting this framework,we are able to give a very clear and elementary proof of the convergence of these solution methods.

A Gas-Kinetic Unified Algorithm for Non-Equilibrium Polyatomic Gas Flows Covering Various Flow Regimes

In this paper,a gas-kinetic unified algorithm(GKUA)is developed to investigate the non-equilibrium polyatomic gas flows covering various regimes.Based on the ellipsoidal statistical model with rotational energy excitation,the computable modelling equation is presented by unifying expressions on the molecular collision relaxing parameter and the local equilibrium distribution function.By constructing the corresponding conservative discrete velocity ordinate method for this model,the conservative properties during the collision procedure are preserved at the discrete level by the numerical method,decreasing the computational storage and time.Explicit and implicit lower-upper symmetric Gauss-Seidel schemes are constructed to solve the discrete hyperbolic conservation equations directly.Applying the new GKUA,some numerical examples are simulated,including the Sod Riemann problem,homogeneous flow rotational relaxation,normal shock structure,Fourier and Couette flows,supersonic flows past a circular cylinder,and hypersonic flow around a plate placed normally.The results obtained by the analytic,experimental,direct simulation Monte Carlo method,and other measurements in references are compared with the GKUA results,which are in good agreement,demonstrating the high accuracy of the present algorithm.Especially,some polyatomic gas non-equilibrium phenomena are observed and analysed by solving the Boltzmann-type velocity distribution function equation covering various flow regimes.