Gas kinetic flux solver based finite volume weighted essentially non-oscillatory scheme for inviscid compressible flows

A high-order gas kinetic flux solver(GKFS)is presented for simulating inviscid compressible flows.The weighted essentially non-oscillatory(WENO)scheme on a uniform mesh in the finite volume formulation is combined with the circular function-based GKFS(C-GKFS)to capture more details of the flow fields with fewer grids.Different from most of the current GKFSs,which are constructed based on the Maxwellian distribution function or its equivalent form,the C-GKFS simplifies the Maxwellian distribution function into the circular function,which ensures that the Euler or Navier-Stokes equations can be recovered correctly.This improves the efficiency of the GKFS and reduces its complexity to facilitate the practical application of engineering.Several benchmark cases are simulated,and good agreement can be obtained in comparison with the references,which demonstrates that the high-order C-GKFS can achieve the desired accuracy.

Three-Dimensional Lattice Boltzmann Model for High-Speed Compressible Flows

A highly efficient three-dimensional （31）） Lattice Boltzmann （LB） model for high-speed compressible flows is proposed. This model is developed from the original one by Kataoka and Tsutahara [Phys. Rev. E 69 （2004） 056702]. The convection term is discretized by the Non-oscillatory, containing No free parameters and Dissipative （NND） scheme, which effectively damps oscillations at discontinuities. To be more consistent with the kinetic theory of viscosity and to further improve the numerical stability, an additional dissipation term is introduced. Model parameters are chosen in such a way that the von Neumann stability criterion is satisfied. The new model is validated by well-known benchmarks, （i） Riemann problems, including the problem with Lax shock tube and a newly designed shock tube problem with high Mach number; （ii） reaction of shock wave on droplet or bubble. Good agreements are obtained between LB results and exact ones or previously reported solutions. The model is capable of simulating flows from subsonic to supersonic and capturing jumps resulted from shock waves.

Flux Limiter Lattice Boltzmann for Compressible Flows

In this paper, a new flux limiter scheme with the splitting technique is successfully incorporated into a multiple-relaxation-time lattice Boltzmann （LB） model for shacked compressible flows. The proposed flux limiter scheme is efficient in decreasing the artificial oscillations and numerical diffusion around the interface. Due to the kinetic nature, some interface problems being difficult to handle at the macroscopic level can be modeled more naturally through the LB method. Numerical simulations for the Richtmyer-Meshkov instability show that with the new model the computed interfaces are smoother and more consistent with physical analysis. The growth rates of bubble and spike present a satisfying agreement with the theoretical predictions and other numerical simulations.

Flux Limiter Lattice Boltzmann Scheme Approach to Compressible Flows with Flexible Specific-Heat Ratio and Prandtl Number

We further develop the lattice Boltzmann （LB） model [Physica A 382 （2007） 502] for compressible flows from two aspects. Firstly, we modify the Bhatnagar Gross Krook （BGK） collision term in the LB equation, which makes the model suitable for simulating flows with different Prandtl numbers. Secondly, the flux limiter finite difference （FLFD） scheme is employed to calculate the convection term of the LB equation, which makes the unphysical oscillations at discontinuities be effectively suppressed and the numerical dissipations be significantly diminished. The proposed model is validated by recovering results of some well-known benchmarks, including （i） The thermal Couette flow; （ii） One- and two-dlmenslonal FLiemann problems. Good agreements are obtained between LB results and the exact ones or previously reported solutions. The flexibility, together with the high accuracy of the new model, endows the proposed model considerable potential for tracking some long-standing problems and for investigating nonlinear nonequilibrium complex systems.

An h-Adaptivity DG Method on Locally Curved Tetrahedral Mesh for Solving Compressible Flows

For the numerical simulation of compressible flows,normally different mesh sizes are expected in different regions.For example,smaller mesh sizes are required to improve the local numerical resolution in the regions where the physical variables vary violently(for example,near the shock waves or in the boundary layers)and larger elements are expected for the regions where the solution is smooth.h-adaptive mesh has been widely used for complex flows.However,there are two difficulties when employing h-adaptivity for high-order discontinuous Galerkin(DG)methods.First,locally curved elements are required to precisely match the solid boundary,which significantly increases the difficulty to conduct the"refining"and"coarsening"operations since the curved information has to be maintained.Second,h-adaptivity could break the partition balancing,which would significantly affect the efficiency of parallel computing.In this paper,a robust and automatic h-adaptive method is developed for high-order DG methods on locally curved tetrahedral mesh,for which the curved geometries are maintained during the h-adaptivity.Furthermore,the reallocating and rebalancing of the computational loads on parallel clusters are conducted to maintain the parallel efficiency.Numerical results indicate that the introduced h-adaptive method is able to generate more reasonable mesh according to the structure of flow-fields.

A GHOST FLUID BASED FRONT TRACKING METHOD FOR MULTIMEDIUM COMPRESSIBLE FLOWS

Recent years the modify ghost fluid method （MGFM） and the real ghost fluid method （RGFM） based on Riemann problem have been developed for multimedium compressible flows. According to authors, these methods have only been used with the level set technique to track the interface. In this paper, we combine the MCFM and the RGFM respectively with front tracking method, for which the fluid interfaces are explicitly tracked by connected points. The method is tested with some one-dimensional problems, and its applicability is also studied. Furthermore, in order to capture the interface more accurately, especially for strong shock impacting on interface, a shock monitor is proposed to determine the initial states of the Riemann problem. The present method is applied to various one- dimensional problems involving strong shock-interface interaction. An extension of the present method to two dimension is also introduced and preliminary results are given.

Gas flow characteristics of argon inductively coupled plasma and advections of plasma species under incompressible and compressible flows

In this work, incompressible and compressible flows of background gas are characterized in argon inductively coupled plasma by using a fluid model, and the respective influence of the two flows on the plasma properties is specified. In the incompressible flow, only the velocity variable is calculated, while in the compressible flow, both the velocity and density variables are calculated. The compressible flow is more realistic; nevertheless, a comparison of the two types of flow is convenient for people to investigate the respective role of velocity and density variables. The peripheral symmetric profile of metastable density near the chamber sidewall is broken in the incompressible flow. At the compressible flow, the electron density increases and the electron temperature decreases. Meanwhile, the metastable density peak shifts to the dielectric window from the discharge center, besides for the peripheral density profile distortion, similar to the incompressible flow.The velocity profile at incompressible flow is not altered when changing the inlet velocity, whereas clear peak shift of velocity profile from the inlet to the outlet at compressible flow is observed as increasing the gas flow rate. The shift of velocity peak is more obvious at low pressures for it is easy to compress the rarefied gas. The velocity profile variations at compressible flow show people the concrete residing processes of background molecule and plasma species in the chamber at different flow rates. Of more significance is it implied that in the usual linear method that people use to calculate the residence time, one important parameter in the gas flow dynamics, needs to be rectified. The spatial profile of pressure simulated exhibits obvious spatial gradient. This is helpful for experimentalists to understand their gas pressure measurements that are always taken at the chamber outlet. At the end, the work specification and limitations are listed.

A multiple-relaxation-time lattice Boltzmann method for high-speed compressible flows

This paper presents a coupling compressible model of the lattice Boltzmann method. In this model, the multiplerelaxation-time lattice Boltzmann scheme is used for the evolution of density distribution functions, whereas the modified single-relaxation-time （SRT） lattice Boltzmann scheme is applied for the evolution of potential energy distribution functions. The governing equations are discretized with the third-order Monotone Upwind Schemes for scalar conservation laws finite volume scheme. The choice of relaxation coefficients is discussed simply. Through the numerical simulations, it is found that compressible flows with strong shocks can be well simulated by present model. The numerical results agree well with the reference results and are better than that of the SRT version.

Simulation of thermoacoustic waves by a pressure-based algorithm for compressible flows

A modified SIMPLEC method which can solve compressible flows at low Mach number is introduced and used to study thermoacoustic waves induced by a rapid change of temperature at a solid wall and alternating- direction flows generated by thermoacoustic effects in a ta- pered resonator. The results indicate that the algorithm adopted in this paper can be used for calculating com- pressible flows and thermoacoustic waves. It is found that the pressure and velocity in the resonator behave as stand- ing waves, and the tapered resonator can suppress high- frequency harmonic waves as observed in a cylindrical res- onator.

Shocklets in compressible flows

The mechanism of shocklets is studied theoretically and numerically for the stationary fluid, uniform compressible flow, and boundary layer flow. The conditions that trigger shock waves for sound wave, weak discontinuity, and Tollmien-Schlichting （T-S） wave in compressible flows are investigated. The relations between the three types of waves and shocklets are further analyzed and discussed. Different stages of the shocklet formation process are simulated. The results show that the three waves in compressible flows will transfer to shocklets only when the initial disturbance amplitudes are greater than the certain threshold values. In compressible boundary layers, the shocklets evolved from T-S wave exist only in a finite region near the surface instead of the whole wavefront.

An Eulerian SPH method with WENO reconstruction for compressible and incompressible flows

While Eulerian smoothed particle hydrodynamics(SPH)method has received increasing attention in scientific and industrial communities owing to its high spatial accuracy,it exhibits excessive numerical dissipation due to the fact that the flux is derived in particle pair pattern.In this paper,we adopt a one-dimensional weighted essentially non-oscillatory(WENO)reconstruction to reduce the numerical dissipation and improve the overall accuracy particularly in capturing the contact discontinuity.The underlying principle is to construct a 4-point stencil along the interacting line of each particle pair and then the WENO scheme is applied to reconstruct the initial states of the Riemann problem which determines the flow flux.A set of benchmark tests for both compressible and incompressible flows are studied to investigate the accuracy,robustness and versatility of the proposed Eulerian SPH method with the WENO reconstruction(ESPH-WENO).

ENERGY CONSERVATION FOR THE WEAK SOLUTIONS TO THE 3D COMPRESSIBLE NEMATIC LIQUID CRYSTAL FLOW

In this paper,we establish some regularity conditions on the density and velocity fields to guarantee the energy conservation of the weak solutions for the three-dimensional compressible nematic liquid crystal flow in the periodic domain.

Improvement and application of wall function boundary condition for high-speed compressible flows

In order to develop a wall function boundary condition for high-speed flows so as to reduce the grid-dependence of the simula- tion for the skin friction and heat flux, a research was performed to improve the compressible wall function boundary condition proposed by Nichols. Values of parameters in the velocity law-of-the-wall were revised according to numerical experiments and the expression of temperature law-of-the-wall was modified based on theoretical analysis and numerical simulation. Be- sides, the formula of the heat conduction term in near-wall region was derived so that the coupling between the wall function boundary condition and CFD code was realized more accurately. Whereafter, the application study of the modified wall func- tion was carried out. The numerical case of supersonic turbulent boundary layer on a flat plate illustrated that the modified wall function produces reasonable results of skin friction and heat flux, and profiles of velocity, temperature and turbulent eddy viscosity for coarse grids with the initial wall spacing of y＋〈400, and that the modifications to the original wall function can obviously improve the simulation precision. As for the application of separation flows, it was found from the numerical cases of supersonic cavity flow and hypersonic axisymmetric compression comer that the compressible velocity law-of-the-wall originally established based on the fully-developed attached turbulent boundary layer approximately holds in the near-wall re- gion inside the separation flows, which ensures that reliable skin friction and heat flux can be given by the wall function inside the separation flows, while for the region near separation and reattachment points, the wall function gives results with a rela- tively large error, because the velocity law-of-the-wall used in the wall function takes on obvious deviation from the real ve- locity profiles near the separation and reattachment points.

Lattice Boltzmann modeling and simulation of compressible flows

In this mini-review we summarize the progress of Lattice Boltzmann （LB） modeling and simulating compressible flows in our group in recent years. Main contents include （i） Single-Relaxation-Time （SRT） LB model supplemented by additional viscosity, （ii） Multiple-Relaxation-Time （MRT） LB model, and （iii） LB study on hydrodynamic instabilities. The former two belong to improvements of physical modeling and the third belongs to simulation or application. The SRT-LB model sup- plemented by additional viscosity keeps the original framework of Lattice Bhatnagar-Gross Krook （LBGK）. So, it is easier and more convenient for previous SRT-LB users. The MRT-LB is a com- pletely new framework for physical modeling. It significantly extends the range of LB applications. The cost is longer computational time. The developed SRT-LB and MRT-LB are complementary from the sides of convenience and applicability.

A unified theory for gas dynamics and aeroacoustics in viscous compressible flows.PartⅡ.Sources on solid boundary

This work attempts to extend the fundamental theory for classic gas dynamics to viscous compressible flow,of which aeroacoustics will naturally be a special branch.As a continuation of Part I.Unbounded fluid(Mao et al.,2022),this paper studies the source of longitudinal field at solid boundary,caused by the on-wall kinematic and viscous dynamic coupling of longitudinal and transverse processes.We find that at this situation the easiest choice for the two independent thermodynamic variables is the dimensionless pressure P and temperature T.The two-level structure of boundary dynamics of longitudinal field is obtained by applying the continuity equation and its normal derivative to the surface.We show that the boundary dilatation flux represents faithfully the boundary production of vortex sound and entropy sound,and the mutual generation mechanism of the longitudinal and transverse fields on the boundary does not occur symmetrically"at the samc level,but appears along a zigzag route.At the first level,it is the pressure gradient that generates vorticity unidirectionally;while at the second level,it is the vorticity that generates dilatation unidirectionally.

Development and Comparative Studies of Three Non-free Parameter Lattice Boltzmann Models for Simulation of Compressible Flows

This paper at first shows the details of finite volume-based lattice Boltzmann method(FV-LBM)for simulation of compressible flows with shock waves.In the FV-LBM,the normal convective flux at the interface of a cell is evaluated by using one-dimensional compressible lattice Boltzmann model,while the tangential flux is calculated using the same way as used in the conventional Euler solvers.The paper then presents a platform to construct one-dimensional compressible lattice Boltzmann model for its use in FV-LBM.The platform is formed from the conservation forms of moments.Under the platform,both the equilibrium distribution functions and lattice velocities can be determined,and therefore,non-free parameter model can be developed.The paper particularly presents three typical non-free parameter models,D1Q3,D1Q4 and D1Q5.The performances of these three models for simulation of compressible flows are investigated by a brief analysis and their application to solve some one-dimensional and two-dimensional test problems.Numerical results showed that D1Q3 model costs the least computation time and D1Q4 and D1Q5 models have the wider application range of Mach number.From the results,it seems that D1Q4 model could be the best choice for the FVLBM simulation of hypersonic flows.

HIGH-ORDER I-STABLE CENTERED DIFFERENCE SCHEMES FOR VISCOUS COMPRESSIBLE FLOWS

In this paper we present high-order I-stable centered difference schemes for the numerical simulation of viscous compressible flows. Here I-stability refers to time discretizations whose linear stability regions contain part of the imaginary axis. This class of schemes has a numerical stability independent of the cell-Reynolds number Re, thus allows one to simulate high Reynolds number flows with relatively larger Re, or coarser grids for a fixed Re. On the other hand, Re cannot be arbitrarily large if one tries to obtain adequate numerical resolution of the viscous behavior. We investigate the behavior of high-order I-stable schemes for Burgers' equation and the compressible Navier-Stokes equations. We demonstrate that, for the second order scheme, Re ≤ 3 is an appropriate constraint for numerical resolution of the viscous profile, while for the fourth-order schemes the constraint can be relaxed to Re ≤ 6.0ur study indicates that the fourth order scheme is preferable: better accuracy, higher resolution, and larger cell-Reynolds numbers.

Multi-scale Equations for Compressible Turbulent Flows

The short-range property of interactions between scales in the compressible turbulent flow was examined. An estimation of the short-range scale scope and some formulae for the short-range eddy stress and heat transfer etc. were given. A concept of resonant-range interactions between extremely contiguous scales was introduced and some formulae for the resonant-range eddy stress and heat transfer etc. were also given. Multi-scale equations for the compressible turbulent flows were presented. The multi-scale equations are approximately closed and do not contain any empirical constants. The compressibility effects on turbulence are determined by the Farve averaged variables and the nonlinear relationships between the Farve- and physical-averaged variables.

Fourier Collocation and Reduced Basis Methods for Fast Modeling of Compressible Flows

A projection-based reduced order model(ROM)based on the Fourier collocation method is proposed for compressible flows.The incorporation of localized artificial viscosity model and filtering is pursued to enhance the robustness and accuracy of the ROM for shock-dominated flows.Furthermore,for Euler systems,ROMs built on the conservative and the skew-symmetric forms of the governing equation are compared.To ensure efficiency,the discrete empirical interpolation method(DEIM)is employed.An alternative reduction approach,exploring the sparsity of viscosity is also investigated for the viscous terms.A number of one-and two-dimensional benchmark cases are considered to test the performance of the proposed models.Results show that stable computations for shock-dominated cases can be achieved with ROMs built on both the conservative and the skew-symmetric forms without additional stabilization components other than the viscosity model and filtering.Under the same parameters,the skew-symmetric form shows better robustness and accuracy than its conservative counterpart,while the conservative form is superior in terms of efficiency.

Artificial neural network-based subgrid-scale models for LES of compressible turbulent channel flow

Fully connected neural networks(FCNNs)have been developed for the closure of subgrid-scale(SGS)stress and SGS heat flux in large-eddy simulations of compressible turbulent channel flow.The FCNNbased SGS model trained using data with Mach number Ma=3.0 and Reynolds number Re=3000 was applied to situations with different Mach numbers and Reynolds numbers.The input variables of the neural network model were the filtered velocity gradients and temperature gradients at a single spatial grid point.The a priori test showed that the FCNN model had a correlation coefficient larger than 0.91 and a relative error smaller than 0.43,with much better reconstructions of SGS unclosed terms than the dynamic Smagorinsky model(DSM).In a posteriori test,the behavior of the FCNN model was marginally better than that of the DSM in predicting the mean velocity profiles,mean temperature profiles,turbulent intensities,total Reynolds stress,total Reynolds heat flux,and mean SGS flux of kinetic energy,and outperformed the Smagorinsky model.