THE RIEMANN PROBLEM WITH DELTA INITIAL DATA FOR THE ONE-DIMENSIONAL CHAPLYGIN GAS EQUATIONS

In this article, we study the Riemann problem with delta initial data for the one-dimensional Chaplygin gas equations. Under the generalized Rankine-Hugoniot conditions and the entropy condition, we constructively obtain the global existence of generalized solutions that explicitly exhibit four kinds of different structures. Moreover, we obtain the stability of generalized solutions by making use of the perturbation of the initial data.

An RKDG finite element method for the one-dimensional inviscid compressible gas dynamics equations in a Lagrangian coordinate

In this paper,Runge-Kutta Discontinuous Galerkin（RKDG） finite element method is presented to solve the onedimensional inviscid compressible gas dynamic equations in a Lagrangian coordinate.The equations are discretized by the DG method in space and the temporal discretization is accomplished by the total variation diminishing Runge-Kutta method.A limiter based on the characteristic field decomposition is applied to maintain stability and non-oscillatory property of the RKDG method.For multi-medium fluid simulation,the two cells adjacent to the interface are treated differently from other cells.At first,a linear Riemann solver is applied to calculate the numerical ？ux at the interface.Numerical examples show that there is some oscillation in the vicinity of the interface.Then a nonlinear Riemann solver based on the characteristic formulation of the equation and the discontinuity relations is adopted to calculate the numerical ？ux at the interface,which suppresses the oscillation successfully.Several single-medium and multi-medium fluid examples are given to demonstrate the reliability and efficiency of the algorithm.

ENTROPY CONSISTENCY OF LARGE TIME STEP SCHEMES FOR ISENTROPIC EQUATIONS OF GAS DYNAMICS

In this paper, We show for isentropic equations of gas dynamics with adiabatic exponent gamma=3 that approximations of weak solutions generated by large time step Godunov's scheme or Glimm's scheme give entropy solution in the limit if Courant number is less than or equal to 1.

Numerical Computation for Time-fractional Gas Dynamics Equations by Fractional Reduced Differential Transforms Method

The present article is concerned with the implementation of a recent semi-analytical method referred to as fractional reduced differential transform method （FRDTM） for computation of approximate solution of time-fractional gas dynamics equation （TFGDE） arising in shock fronts. In this approach, the fractional derivative is described in the Caputo sense. Four numeric experiments have been carried out to confirm the validity and the efficiency of the method. It is found that the exact or a closed approximate analytical solution of a fractional nonlinear differential equations arising in allied science and engineering can be obtained easily. Moreover, due to its small size of calculation contrary to the other analytical approaches while dealing with a complex and tedious physical problems arising in various branches of natural sciences and engineering, it is very easy to implement.

Gas seepage equation of deep mined coal seams and its application

In order to obtain a gas seepage law of deep mined coal seams, according to the properties of coalbed methane seepage in in-situ stress and geothermal temperature fields, the gas seepage equation of deep mined coal seams with the Klinkenberg effect was obtained by confirming the coatbed methane permeability in in-situ stress and geothermal temperature fields. Aimed at the condition in which the coal seams have or do not have an outcrop and outlet on the ground, the application of the gas seepage equation of deep mined coal seams in in-situ stress and geothermal temperature fields on the gas pressure calculation of deep mined coal seams was investigated. The comparison between calculated and measured results indicates that the calculation method of gas pressure, based on the gas seepage equation of deep mined coal seams in in-situ stress and geothermal temperature fields can accu- rately be identical with the measured values and theoretically perfect the calculation method of gas pressure of deep mined coal seams.

The pressure compensation technology of deep-sea sampling based on the real gas state equation

Compressed gas is usually used for the pressure compensation of the deep-sea pressure-maintaining sampler.The pressure and volume of the recovered fluid sample are highly related to the precharged gas. To better understand the behavior of the gas under high pressure, we present a new real gas state equation based on the compression factor Z which was derived from experimental data. Then theoretical calculation method of the pressure and volume of the sample was introduced based on this empirical gas state equation. Finally, the proposed calculation method was well verified by the high-pressure vessel experiment of the sampler under 115 MPa.

Study on the mechanism of interaction for coal and methane gas

Although two moulds for methane gas in coal with the free state and adsorption state have been popularly considered, the derivation between the real methane gas state equation in coal and the perfect gas state equation has been fuzzily considered and the mechanism of interaction for coal aromatics and methane gas molecules has not been understood. Then these problems have been discussed in this paper applied the principle of statistical thermo mechanics and quantum chemistry as well as based on the numerical calculating of experiential data in quantum chemistry. Therefore, it is revealed by research results that the experience state equation for real methane gas in coal, which is put forward in this paper, is closer to actual situation and the interaction process for methane gas adsorption on the surface of coal aromatics can be formulated by Morse potential function. Furthermore it is most stable through this research that the structural mould for methane gas molecule adsorption on the surface of coal nuclear with one gas molecule on top of another aromatics in regular triangle cone has been understood, and it is a physical adsorption for methane gas adsorption with single layer molecule on the surface of coal nuclear.

Gas-kinetic numerical method for solving mesoscopic velocity distribution function equation

A gas-kinetic numerical method for directly solving the mesoscopic velocity distribution function equation is presented and applied to the study of three-dimensional complex flows and micro-channel flows covering various flow regimes. The unified velocity distribution function equation describing gas transport phenomena from rarefied transition to continuum flow regimes can be presented on the basis of the kinetic Boltzmann-Shakhov model equation. The gas-kinetic finite-difference schemes for the velocity distribution function are constructed by developing a discrete velocity ordinate method of gas kinetic theory and an unsteady time-splitting technique from computational fluid dynamics. Gas-kinetic boundary conditions and numerical modeling can be established by directly manipulating on the mesoscopic velocity distribution function. A new Gauss-type discrete velocity numerical integra- tion method can be developed and adopted to attack complex flows with different Mach numbers. HPF paral- lel strategy suitable for the gas-kinetic numerical method is investigated and adopted to solve three-dimensional complex problems. High Mach number flows around three-dimensional bodies are computed preliminarilywith massive scale parallel. It is noteworthy and of practical importance that the HPF parallel algorithm for solving three-dimensional complex problems can be effectively developed to cover various flow regimes. On the other hand, the gas-kinetic numerical method is extended and used to study micro-channel gas flows including the classical Couette flow, the Poiseuillechannel flow and pressure-driven gas flows in twodimensional short micro-channels. The numerical experience shows that the gas-kinetic algorithm may be a powerful tool in the numerical simulation of microscale gas flows occuring in the Micro-Electro-Mechanical System （MEMS）.

An improved theoretical procedure for the pore-size analysis of activated carbon by gas adsorption

Amorphous carbon materials play a vital role in adsorbed natural gas(ANG) storage. One of the key issues in the more prevalent use of ANG is the limited adsorption capacity, which is primarily determined by the porosity and surface characteristics of porous materials. To identify suitable adsorbents, we need a reliable computational tool for pore characterization and, subsequently, quantitative prediction of the adsorption behavior. Within the framework of adsorption integral equation(AIE), the pore-size distribution(PSD) is sensitive to the adopted theoretical models and numerical algorithms through isotherm fitting. In recent years, the classical density functional theory(DFT) has emerged as a common choice to describe adsorption isotherms for AIE kernel construction. However,rarely considered is the accuracy of the mean-field approximation(MFA) commonly used in commercial software. In this work, we calibrate four versions of DFT methods with grand canonical Monte Carlo(GCMC) molecular simulation for the adsorption of CH_4 and CO_2 gas in slit pores at 298 K with the pore width varying from 0.65 to 5.00 nm and pressure from 0.2 to 2.0 MPa. It is found that a weighted-density approximation proposed by Yu(WDA-Yu) is more accurate than MFA and other non-local DFT methods. In combination with the trapezoid discretization of AIE, the WDA-Yu method provides a faithful representation of experimental data, with the accuracy and stability improved by 90.0% and 91.2%, respectively, in comparison with the corresponding results from MFA for fitting CO_2 isotherms. In particular, those distributions in the feature pore width range(FPWR)are proved more representative for the pore-size analysis. The new theoretical procedure for pore characterization has also been tested with the methane adsorption capacity in seven activated carbon samples.

Adaptive Moving Mesh Central-Upwind Schemes for Hyperbolic System of PDEs:Applications to Compressible Euler Equations and Granular Hydrodynamics

We introduce adaptive moving mesh central-upwind schemes for one-and two-dimensional hyperbolic systems of conservation and balance laws.The proposed methods consist of three steps.First,the solution is evolved by solving the studied system by the second-order semi-discrete central-upwind scheme on either the one-dimensional nonuniform grid or the two-dimensional structured quadrilateral mesh.When the evolution step is complete,the grid points are redistributed according to the moving mesh differential equation.Finally,the evolved solution is projected onto the new mesh in a conservative manner.The resulting adaptive moving mesh methods are applied to the one-and two-dimensional Euler equations of gas dynamics and granular hydrodynamics systems.Our numerical results demonstrate that in both cases,the adaptive moving mesh central-upwind schemes outperform their uniform mesh counterparts.

Influence of Different Equations of State on Simulation Results of Supercritical CO_(2) Centrifugal Compressor

Supercritical CO_(2)(SCO_(2))Brayton cycle has received more and more attention in the field of power generation due to its high cycle efficiency and compact structure.SCO_(2) compressor is the core component of the cycle,and the improvement of its performance is the key to improving the efficiency of the entire cycle.However,the operation of the SCO_(2) compressor near the critical point has brought many design and operation problems.Based on the Reynolds Averaged Navier-Stokes(RANS)model,the performance and flow field of SCO_(2) centrifugal compressors based on different CO_(2) working fluid models are numerically investigated in this paper.The stability and convergence of the compressor steady-state simulation are also discussed.The results show that the fluid based on the Span-Wanger(SW)equation can obtain a more ideal compressor performance curve and capture a more accurate flow field structure,while the CO_(2) ideal gas is not suitable for the calculation of SCO_(2) centrifugal compressors.But its flow field can be used as the initial flow field for numerical calculation of centrifugal compressor based on CO_(2) real gas.

Direct modeling for computational fluid dynamics

All fluid dynamic equations are valid under their modeling scales, such as the particle mean free path and mean collision time scale of the Boltzmann equation and the hydrodynamic scale of the Navier-Stokes （NS） equations. The current computational fluid dynamics （CFD） focuses on the numerical solution of partial differential equations （PDEs）, and its aim is to get the accurate solution of these governing equations. Under such a CFD practice, it is hard to develop a unified scheme that covers flow physics from kinetic to hydrodynamic scales continuously because there is no such governing equation which could make a smooth transition from the Boltzmann to the NS modeling. The study of fluid dynamics needs to go beyond the traditional numer- ical partial differential equations. The emerging engineering applications, such as air-vehicle design for near-space flight and flow and heat transfer in micro-devices, do require fur- ther expansion of the concept of gas dynamics to a larger domain of physical reality, rather than the traditional dis- tinguishable governing equations. At the current stage, the non-equilibrium flow physics has not yet been well explored or clearly understood due to the lack of appropriate tools. Unfortunately, under the current numerical PDE approach, it is hard to develop such a meaningful tool due to the absence of valid PDEs. In order to construct multiscale and multiphysics simulation methods similar to the modeling process of con- structing the Boltzmann or the NS governing equations, the development of a numerical algorithm should be based on the first principle of physical modeling. In this paper, instead of following the traditional numerical PDE path, we introduce direct modeling as a principle for CFD algorithm develop- ment. Since all computations are conducted in a discretized space with limited cell resolution, the flow physics to be mod- eled has to be done in the mesh size and time step scales. Here, the CFD is more or less a direct construction of dis- crete numerical evolution equations, where the mesh size and time step will play dynamic roles in the modeling process. With the variation of the ratio between mesh size and local particle mean free path, the scheme will capture flow physics from the kinetic particle transport and collision to the hydro- dynamic wave propagation. Based on the direct modeling, a continuous dynamics of flow motion will be captured in the unified gas-kinetic scheme. This scheme can be faithfully used to study the unexplored non-equilibrium flow physics in the transition regime.

A Well-Balanced Kinetic Scheme for Gas Dynamic Equations under Gravitational Field

In this paper,a well-balanced kinetic scheme for the gas dynamic equations under gravitational field is developed.In order to construct such a scheme,the physical process of particles transport through a potential barrier at a cell interface is considered,where the amount of particle penetration and reflection is evaluated according to the incident particle velocity.This work extends the approach of Perthame and Simeoni for the shallow water equations[Calcolo,38(2001),pp.201-231]to the Euler equations under gravitational field.For an isolated system,this scheme is probably the only well-balanced method which can precisely preserve an isothermal steady state solution under time-independent gravitational potential.A few numerical examples are used to validate the above approach.

AnAll-Regime Lagrange-Projection Like Scheme for the Gas Dynamics Equations on Unstructured Meshes

We propose an all regime Lagrange-Projection like numerical scheme for the gas dynamics equations.By all regime,we mean that the numerical scheme is able to compute accurate approximate solutions with an under-resolved discretization with respect to the Mach number M,i.e.such that the ratio between the Mach number M and the mesh size or the time step is small with respect to 1.The key idea is to decouple acoustic and transport phenomenon and then alter the numerical flux in the acoustic approximation to obtain a uniform truncation error in term of M.This modified scheme is conservative and endowed with good stability properties with respect to the positivity of the density and the internal energy.A discrete entropy inequality under a condition on the modification is obtained thanks to a reinterpretation of the modified scheme in the Harten Lax and van Leer formalism.A natural extension to multi-dimensional problems discretized over unstructured mesh is proposed.Then a simple and efficient semi implicit scheme is also proposed.The resulting scheme is stable under a CFL condition driven by the(slow)material waves and not by the(fast)acoustic waves and so verifies the all regime property.Numerical evidences are proposed and show the ability of the scheme to deal with tests where the flow regime may vary from low to high Mach values.

An Energy-Based Unit for the Thermodynamic Temperature

The article suggests to replace the conventional unit of thermodynamic temperature with one based on the energy unit joule by including the gas constant into the temperature definition. The suggestion may be seen as a contribution to the ongoing efforts to redefine base units. The gas constant includes the Boltzmann constant which both will be abrogated and molar heat capacity and entropy will become dimensionless pure numbers. The suggestion has no impact on thermodynamic theory, but will make thermodynamic relations more translucent and easier to grab for students.

Fifth-Order A-WENO Schemes Based on the Adaptive Diffusion Central-Upwind Rankine-Hugoniot Fluxes

We construct new fifth-order alternative WENO(A-WENO)schemes for the Euler equations of gas dynamics.The new scheme is based on a new adaptive diffusion centralupwind Rankine-Hugoniot(CURH)numerical flux.The CURH numerical fluxes have been recently proposed in[Garg et al.J Comput Phys 428,2021]in the context of secondorder semi-discrete finite-volume methods.The proposed adaptive diffusion CURH flux contains a smaller amount of numerical dissipation compared with the adaptive diffusion central numerical flux,which was also developed with the help of the discrete RankineHugoniot conditions and used in the fifth-order A-WENO scheme recently introduced in[Wang et al.SIAM J Sci Comput 42,2020].As in that work,we here use the fifth-order characteristic-wise WENO-Z interpolations to evaluate the fifth-order point values required by the numerical fluxes.The resulting one-and two-dimensional schemes are tested on a number of numerical examples,which clearly demonstrate that the new schemes outperform the existing fifth-order A-WENO schemes without compromising the robustness.

THE ASYMPTOTIC PRESERVING UNIFIED GAS KINETIC SCHEME FOR GRAY RADIATIVE TRANSFER EQUATIONS ON DISTORTED QUADRILATERAL MESHES

In this paper,we consider the multi-dimensional asymptotic preserving unified gas kinetic scheme for gray radiative transfer equations on distorted quadrilateral meshes.Different from the former scheme [J.Comput.Phys.285（2015）,265-279] on uniform meshes,in this paper,in order to obtain the boundary fluxes based on the framework of unified gas kinetic scheme（UGKS）,we use the real multi-dimensional reconstruction for the initial data and the macro-terms in the equation of the gray transfer equations.We can prove that the scheme is asymptotic preserving,and especially for the distorted quadrilateral meshes,a nine-point scheme [SIAM J.SCI.COMPUT.30（2008）,1341-1361] for the diffusion limit equations is obtained,which is naturally reduced to standard five-point scheme for the orthogonal meshes.The numerical examples on distorted meshes are included to validate the current approach.

Recovering Navier–Stokes Equations from Asymptotic Limits of the Boltzmann Gas Mixture Equation

This paper is devoted to the derivation of macroscopic fluid dynamics from the Boltzmann mesoscopic dynamics of a binary mixture of hard-sphere gas particles.Specifically the hydrodynamics limit is performed by employing different time and space scalings.The paper shows that,depending on the magnitude of the parameters which define the scaling,the macroscopic quantities(number density,mean velocity and local temperature)are solutions of the acoustic equation,the linear incompressible Euler equation and the incompressible Navier–Stokes equation.The derivation is formally tackled by the recent moment method proposed by[C.Bardos,et al.,J.Stat.Phys.63(1991)323]and the results generalize the analysis performed in[C.Bianca,et al.,Commun.Nonlinear Sci.Numer.Simulat.29(2015)240].

Non-Local Theory of Bose-Einstein Condensate and Stopping of Light

The problem of an adequate description of the transport processes in Bose-Einstein condensates (CBE), including space-temporal evolution of CBE in a gravitational field is considered. The full nonlocal system for the CBE evolution is delivered including particular case and analytical solutions. We show (analytically) that a black hole can evolve in the Bose-Einstein condensate (CBE) regime. At the same time, there are modes in which black hole flickering occurs. Quantization of the black holes flickering is discovered. The corresponding nonlocal hydrodynamic equations indicated for fermions gas.

The Riemann Problem for Chaplygin Gas Flow in a Duct with Discontinuous Cross-Section

The fluid flows in a variable cross-section duct are nonconservative because of the source term.Recently,the Riemann problem and the interactions of the elementary waves for the compressible isentropic gas in a variable cross-section duct were studied.In this paper,the Riemann problem for Chaplygin gas flow in a duct with discontinuous cross-section is studied.The elementary waves include rarefaction waves,shock waves,delta waves and stationary waves.