VANISHING VISCOSITY FOR NON-ISENTROPIC GAS DYNAMICS WITH INTERACTING SHOCKS

In this paper, we study the vanishing viscosity limit for non-isentropic gas dy- namics with interacting shocks. Given any entropy solution of non-isentropic gas dynamics which consists of two different families of shocks interacting at some positive time, we show that such solution is the vanishing viscosity limit of a family of smooth global solutions for a viscous system of conservation law. We remark that, after the interacting time, not only shocks but also contact discontinuity are generated.

EVOLUTION AND INTERACTION OF δ-WAVES IN THE ZERO-PRESSURE GAS DYNAMICS SYSTEM

Evolution and interaction of plane waves of the multidimensional zero-pressure gas dynamics system leads to the study of the corresponding one dimensional system.In this paper,we study the initial value problem for one dimensional zero-pressure gas dynamics system.Here the first equation is the Burgers equation and the second one is the continuity equation.We consider the solution with initial data in the space of bounded Borel measures.First we prove a general existence result in the algebra of generalized functions of Colombeau.Then we study in detail special solutions withδ-measures as initial data.We study interaction of waves originating from initial data concentrated on two point sources and interaction with classical shock/rarefaction waves.This gives an understanding of plane-wave interactions in the multidimensional case.We use the vanishing viscosity method in our analysis as this gives the physical solution.

Model of gas exchange dynamics for modified-atmosphere packages containing fresh produce

A model for modified-atmosphere packaging （MAP） systems containing fruits and vegetables was developed.The computer simulation was performed to predict the gas mass concentrations inside the packages and was successfully verified by experiments with yellow peaches at 5,15 and 25 ℃ using two types of packaging films.A Michaelis-Menten type respiration model with noncompetitive inhibition mechanism due to CO2 was adopted while the respiration rates were measured with an improved permeable system method suitable for either steady or unsteady state.The applicability of the model in the design of MAP systems was demonstrated with a calculation to evaluate film specification and equilibrium concentrations of O2 and CO2 in the package containing yellow peaches.

Analysis of CH_(4) and H_(2) Adsorption on Heterogeneous Shale Surfaces Using aMolecular Dynamics Approach

Determining the adsorption of shale gas on complex surfaces remains a challenge in molecular simulation studies.Difficulties essentially stem from the need to create a realistic shale structure model in terms of mineral heterogeneityand multiplicity.Moreover,precise characterization of the competitive adsorption of hydrogen andmethane in shale generally requires the experimental determination of the related adsorptive capacity.In thisstudy,the adsorption of adsorbates,methane(CH_(4)),and hydrogen(H_(2))on heterogeneous shale surface modelsof Kaolinite,Orthoclase,Muscovite,Mica,C_(60),and Butane has been simulated in the frame of a moleculardynamic’s numerical technique.The results show that these behaviors are influenced by pressure and potentialenergy.On increasing the pressure from 500 to 2000 psi,the sorption effect for CH_(4)significantly increasesbut shows a decline at a certain stage(if compared to H_(2)).The research findings also indicate that raw shalehas a higher capacity to adsorb CH_(4)compared to hydrogen.However,in shale,this difference is negligible.

High density gas state at water/graphite interface studied by molecular dynamics simulation

In this paper molecular dynamics simulations are performed to study the accumulation behaviour of N2 and H2 at water/graphite interface under ambient temperature and pressure. It finds that both N2 and H2 molecules can accumulate at the interface and form one of two states according to the ratio of gas molecules number to square of graphite surface from our simulation results： gas films （pancake-like） for a larger ratio and nanobubbles for a smaller ratio. In addition, we discuss the stabilities of nanobubbles at different environment temperatures. Surprisingly, it is found that the density of both kinds of gas states can be greatly increased, even comparable with that of the liquid N2 and liquid H2. The present results are expected to be helpful for the understanding of the stable existence of gas film （pancake-like） and nanobubbles.

Modeling and scenario prediction of a natural gas demand system based on a system dynamics method

Based on the study of the relationship between structure and feedback of China’s natural gas demand system, this paper establishes a system dynamics model. In order to simulate the total demand and consumption structure of natural gas in China, we set up seven scenarios by changing some of the parameters of the model. The results showed that the total demand of natural gas would increase steadily year by year and reach in the range from 3600 to 4500 billion cubic meters in 2035. Furthermore, in terms of consumption structure, urban gas consumption would still be the largest term, followed by the gas consumption as industrial fuel, gas power generation and natural gas chemical industry. In addition, compared with the population growth, economic development still plays a dominant role in the natural gas demand growth, the impact of urbanization on urban gas consumption is significant, and the promotion of natural gas utilization technology can effectively reduce the total consumption of natural gas.

GLOBAL L~∞ SOLUTIONS TO SYSTEM OF ISENTROPIC GAS DYNAMICS IN A DIVERGENT NOZZLE WITH FRICTION

In this article, we study the global L^∞ entropy solutions for the Cauchy problem of system of isentropic gas dynamics in a divergent nozzle with a friction. Especially when the adiabatic exponent γ=3, we apply for the maximum principle to obtain the L^∞ estimates w(ρ^δ,ε, u^δ,ε)≤ B(t) and z(ρ^δ,ε, u^δ,ε)≤ B(t) for the viscosity solutions (ρ^δ,ε, u^δ,ε), where B(t) is a nonnegative bounded function for any finite time t. This work, in the special case γ≥ 3, extends the previous works, which provided the global entropy solutions for the same Cauchy problem with the restriction w(ρ^δ,ε, u^δ,ε)≤ 0 or z(ρ^δ,ε, u^δ,ε)≤ 0.

CONVERGENCE OF THE LAX-FRIEDRICHS SCHEME FOR ISENTROPIC GAS DYNAMICS (Ⅰ)

A compactness frame of the Lax-Friedrichs scheme for the equations of gas dynamics is obtained by using some embedding theorems and an analysis of the difference scheme and the weak entropy.

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.

Desired Dynamic Equation for Primary Frequency Modulation Control of Gas Turbines

Gas turbines play core roles in clean energy supply and the construction of comprehensive energy systems.The control performance of primary frequency modulation of gas turbines has a great impact on the frequency control of the power grid.However,there are some control difficulties in the primary frequency modulation control of gas turbines,such as the coupling effect of the fuel control loop and speed control loop,slow tracking speed,and so on.To relieve the abovementioned difficulties,a control strategy based on the desired dynamic equation proportional integral(DDE-PI)is proposed in this paper.Based on the parameter stability region,a parameter tuning procedure is summarized.Simulation is carried out to address the ease of use and simplicity of the proposed tuning method.Finally,DDE-PI is applied to the primary frequency modulation system of an MS6001B heavy-duty gas turbine.The simulation results indicate that the gas turbine with the proposed strategy can obtain the best control performance with a strong ability to deal with system uncertainties.The proposed method shows good engineering application potential.

Natural gas density under extremely high pressure and high temperature: Comparison of molecular dynamics simulation with corresponding state model

This work applied molecular dynamics(MD)simulation to calculate densities of natural gas mixtures at extremely high pressure(>138 MPa)and high temperature(>200℃)conditions(x HPHT)to bridge the knowledge and technical gaps between experiments and classical theories.The experimental data are scarce at these conditions which are also out of assumptions for classical predictive correlations,such as the Dranchuk&Abou-Kassem(DAK)equation of state(EOS).Force fields of natural gas components were carefully chosen from literatures and the simulation results are validated with experimental data.The largest relative error is 2.67%for pure hydrocarbons,2.99%for C1/C3 mixture,7.85%for C1/C4 mixture,and 8.47%for pure H2S.These satisfactory predictions demonstrate that the MD simulation approach is reliable to predict natural-and acid-gases thermodynamic properties.The validated model is further used to generate data for the study of the EOS with pressure up to 276 MPa and temperature up to 573 K.Our results also reveal that the Dranchuk&Abou-Kassem(DAK)EOS is capable of predicting natural gas compressibility to a satisfactory accuracy at x HPHT conditions,which extends the confidence range of the DAK EOS.

A SUPPLEMENT TO THE PAPERS CONVERGENCE OF THE LAX-FRIEDRICHS SCHEME FOR ISENTROPIC GAS DYNAMICS (Ⅱ)-(Ⅲ)

In the papers (Ⅱ)—(Ⅲ), we studied the global existence of the solutions for the isentropic equation of gas dynamics

STRONG SHOCKS IN NONISENTROPIC GAS DYNAMICS

In nonisentropic gas dynamics, with general equations of state, strong cornpressive shocks satisfying the Liu E-condition may violate the Second Law of thermodynamics.

THE VACUUM IN NONISENTROPIC GAS DYNAMICS

We investigate the vacuum in nonisentropic gas dynamics in one space vari- able, with the most general equation of states allowed by thermodynamics. We recall physical constraints on the equations of state and give explicit and easily checkable conditions under which vacuums occur in the solution of the Riemann problem. We then present a class of models for which the Riemann problem admits unique global solutions without vacuums.

Computational fluid dynamics simulation of gas dispersion in complex facilities using Kit Fox field experiments:Validation and statistical evaluation

Gas release and its dispersion is a major concern in chemical industries.In order to manage and mitigate the risk of gas dispersion and its consequences,it is necessary to predict gas dispersion behavior and its concentration at various locations upon emission.Therefore,models and commercial packages such as Phast and ALOHA have been developed.Computational fluid dynamics(CFD)can be a useful tool to simulate gas dispersion in complex areas and conditions.The validation of the models requires the employment of the experimental data from filed and wind tunnel experiments.It appears that the use of the experimental data to validate the CFD method that only includes certain monitor points and not the entire domain can lead to unreliable results for the intended areas of concern.In this work,some of the trials of the Kit Fox field experiment,which provided a wide-range database for gas dispersion,were simulated by CFD.Various scenarios were considered with different mesh sizes,physical conditions,and types of release.The results of the simulations were surveyed in the whole domain.The data matching each scenario was varied by the influence of the dominant displacement force(wind or diffusivity).Furthermore,the statistical parameters suggested for the heavy gas dispersion showed a dependency on the lower band of gas concentration.Therefore,they should be used with precaution.Finally,the results and computation cost of the simulation could be affected by the chosen scenario,the location of the intended points,and the release type.

A NOTE ON THE EXISTENCE AND NONEXISTENCE OF GLOBALLY BOUNDED CLASSICAL SOLUTIONS FOR NONISENTROPIC GAS DYNAMICS

Existence of globally bounded classical solution for nonisentropic gas dynamics system has long been studied, especially in the case of polytropic gas. In [4], Liu claimed that sufficient condition has been established. However, the authors find that the argument he used is not true in general. In this article, the authors give a counter example of his argument. Hence, his claim is not valid. The authors believe that it is difficult to impose general conditions on the initial data to obtain globally bounded classical solution.

MOLECULAR DYNAMICS SIMULATIONS OF FILLED AND EMPTY CAGE-LIKE WATER CLUSTERS IN LIQUID WATER AND THEIR SIGNIFICANCE TO GAS HYDRATE FORMATION MECHANISMS

Molecular dynamics simulations are performed to observe the evolutions of 512 and 51262 cage-like water clusters filled with or without a methane molecule immersed in bulk liquid water at 250 K and 230 K. The lifetimes of these clusters are calculated according to their Lindemann index δ (t) using the criteria of δ≥0.07. For both the filled and empty clusters, we find the dynamics of bulk water determines the lifetimes of cage-like water clusters, and that the lifetime of 512 62 cage-like cluster is the same as that of 512 cage-like cluster. Although the methane molecule indeed makes the filled cage-like cluster more stable than the empty one, the empty cage-like cluster still has chance to be long-lived compared with the filled clusters. These observations support the labile cluster hypothesis on the formation mechanisms of gas hydrates.

SOME EXACT SOLUTIONS OF 3-DIMENSIONAL ZERO-PRESSURE GAS DYNAMICS SYSTEM

The 3-dimensional zero-pressure gas dynamics system appears in the modeling for the large scale structure formation in the universe. The aim of this paper is to construct spherically symmetric solutions to the system. The radial component of the velocity and density satisfy a simpler one dimensional problem. First we construct explicit solutions of this one dimensional case with initial and boundary conditions. Then we get special radial solutions with different behaviours at the origin.

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.

CONVERGENCE OF THE APPROXIMATE SOLUTIONS TO ISENTROPIC GAS DYNAMICS

This paper gives four pairs of entropies (η_i, q_i) (i=1, 2, 3, 4) to the isentropic gas dynamics equations ρ_t+(ρu)_x=0 (ρu)_t+(ρu^2+p(ρ))_x=0 p(ρ)=k^2ρ~γ,1<γ<3。 when all the function equations are satisfied