Adopting the classical theory of hydrocodes,the constitutive relations of concretes are separated into an equation of state(EoS)which describes the volumetric behavior of concrete material and a strength model which d...Adopting the classical theory of hydrocodes,the constitutive relations of concretes are separated into an equation of state(EoS)which describes the volumetric behavior of concrete material and a strength model which depicts the shear properties of concrete.The experiments on the EoS of concrete is always challenging due to the technical difficulties and equipment limitations,especially for the specimen size effect on the EoS.Although some researchers investigate the shock properties of concretes by fly-plate impact tests,the specimens used in their tests are usually in one size.In this paper,the fly-plate impact tests on concrete specimens with different sizes are performed to investigate the size effect on the shock properties of concrete materials.The mechanical background of the size effect on the shock properties are revealed,which is related to the lateral rarefaction effect and the deviatoric stress produced in the specimen.According to the tests results,the modified EoS considering the size effect on the shock properties of concrete are proposed,which the bulk modulus of concrete is unpredicted by up to 20% if size effects are not accounted for.展开更多
Based on the Hugenholtz-Van Hove theorem six basic quantities of the EoS in isospin asymmetric nuclear matter are expressed in terms of the nucleon kinetic energy t(k),the isospin symmetric and asymmetric parts of the...Based on the Hugenholtz-Van Hove theorem six basic quantities of the EoS in isospin asymmetric nuclear matter are expressed in terms of the nucleon kinetic energy t(k),the isospin symmetric and asymmetric parts of the single-nucleon potentials U_(0)(ρ,k)and U_(sym,i)(ρ,k).The six basic quantities include the quadratic symmetry energy E_(sym,2)(ρ),the quartic symmetry energy E_(sym,4)(ρ),their corresponding density slopes L_(2)(ρ)and L_(4)(ρ),and the incompressibility coefficients K_(2)(ρ)and K_(4)(ρ).By using four types of well-known effective nucleon-nucleon interaction models,namely the BGBD,MDI,Skyrme,and Gogny forces,the density-and isospin-dependent properties of these basic quantities are systematically calculated and their values at the saturation density q_(0)are explicitly given.The contributions to these quantities from t(k)U_(0)(ρ,k),and U_(sym,i)(ρ,k)are also analyzed at the norma nuclear density q_(0).It is clearly shown that the first-order asymmetric term U_(sym,1)(ρ,k)(also known as the symmetry potential in the Lane potential)plays a vital role in determining the density dependence of the quadratic symmetry energy E_(sym,2)(ρ).It is also shown that the contributions from the high-order asymmetric parts of the single-nucleon potentials(U_(sym,i)(ρ,k)with i>1)cannot be neglected in the calculations of the other five basic quantities Moreover,by analyzing the properties of asymmetric nuclear matter at the exact saturation densityρ_(sat)(δ),the corresponding quadratic incompressibility coefficient is found to have a simple empirical relation K_(sat,2)=K_(2)(ρ_(0))-4.14L_(2)(ρ_(0))展开更多
This paper presents a three-dimensional, three-phase compositional model considering CO2 phase equilibrium between water and oil. In this model, CO2 is mutually soluble in aqueous and hydrocarbon phases, while other c...This paper presents a three-dimensional, three-phase compositional model considering CO2 phase equilibrium between water and oil. In this model, CO2 is mutually soluble in aqueous and hydrocarbon phases, while other components, except water,exist in hydrocarbon phase. The Peng–Robinson(PR) equation of state and the Wong–Sandler mixing rule with non-random two-liquid parameters are used to calculate CO2 fugacity in the aqueous phase. One-dimensional and three-dimensional CO2 flooding examples show that a significant amount of injected CO2 is dissolved in water. Our simulation shows 7% of injected CO2 can be dissolved in the aqueous phase, which delays oil recovery by 4%. The gas rate predicted by the model is smaller than the conventional model as long as water is undersaturated by CO2, which can be considered as 'lost' in the aqueous phase. The model also predicts that the delayed oil can be recovered after the gas breakthrough, indicating that delayed oil is hard to recover in field applications. A three-dimensional example reveals that a highly stratified reservoir causes uneven displacement and serious CO2 breakthrough. If mobility control measures like water alternating gas are undertaken, the solubility e ects will be more pronounced than this example.展开更多
In this paper we introduce the wide regime equation of state(WEOS)developed in Institute of Applied Physics and Computational Mathematics(IAPCM).A semi-empirical model of the WEOS is given by a thermodynamically compl...In this paper we introduce the wide regime equation of state(WEOS)developed in Institute of Applied Physics and Computational Mathematics(IAPCM).A semi-empirical model of the WEOS is given by a thermodynamically complete potential of the Helmholtz free energy which combines several theoretical models and has some adjustable parameters calibrated via some experimental and theoretical data.The validation methods of the equation of state in wide regime are presented using copper as a prototype.The results of the WEOS are well consistent with the available theoretical and experimental data,including ab initio cold curve under compression,isotherm,Hugoniot,off-Hugoniot and sound velocity data.It enhances our confidence in the accuracy of the WEOS,which is very important for the validation and verification of equation of state in high temperature and pressure technology.展开更多
In this paper, the Martin-Hou equation of state is derived by using a power series representation of radial distribution function and an analytic representation of multi-section potential based on the Barker-Henderso...In this paper, the Martin-Hou equation of state is derived by using a power series representation of radial distribution function and an analytic representation of multi-section potential based on the Barker-Henderson hard-particle perturbation theory including high-order terms. In the derivation, a theoretical form of Martin-Hou equation was obtained. It had a similar form and the same capability to predict P-V-T properties as the Martin-Hou equation and no additional data were required for evaluating the constants. The characteristic constants of the theoretical expression have certain relationships with the molecular parameters.展开更多
The Galilean invariance and the induced thermo-hydrodynamics of the lattice Boltzmann Bhatnagar-Gross-Krook model are proposed together with their rigorous theoretical background. From the viewpoint of group invarianc...The Galilean invariance and the induced thermo-hydrodynamics of the lattice Boltzmann Bhatnagar-Gross-Krook model are proposed together with their rigorous theoretical background. From the viewpoint of group invariance, recovering the Galilean invariance for the isothermal lattice Boltzmann Bhatnagar-Gross-Krook equation (LBGKE) induces a new natural thermal-dynamical system, which is compatible with the elementary statistical thermodynamics.展开更多
Within the approximations of spherical lattice cell, central-field, and relativistic Fermi statis- tics, an algorithm with average atom model is presented to calculate the electronic energy levels and equation of stat...Within the approximations of spherical lattice cell, central-field, and relativistic Fermi statis- tics, an algorithm with average atom model is presented to calculate the electronic energy levels and equation of state for hot and dense matter at arbitrary densities and temperatures. Choosing Zink's analytical potential as initial potential, we have solved the Dirac-Slater equation which satisfies the Weigner-Seitz boundary condition. The electronic energy bands are not taken into account. Tak- ing energy level degeneracy as a continuous function of density, we have considered the pressure ionization effects for highly dense matter. Results for ^(13)Al atom are shown.展开更多
Based on the classical Roe method, we develop an interface capture method according to the general equation of state, and extend the single-fluid Roe method to the two-dimensional (2D) multi-fluid flows, as well as ...Based on the classical Roe method, we develop an interface capture method according to the general equation of state, and extend the single-fluid Roe method to the two-dimensional (2D) multi-fluid flows, as well as construct the continuous Roe matrix for the whole flow field. The interface capture equations and fluid dynamic conservative equations are coupled together and solved by using any high-resolution schemes that usually suit for the single-fluid flows. Some numerical examples are given to illustrate the solution of 1D and 2D multi-fluid Riemann problems.展开更多
The thermal equation of state of a natural kyanite has been investigated with a DIA-type,cubic-anvil apparatus(SAM85)combined with an energy-dispersive synchrotron X-ray radiation technique up to 8.55 GPa and 1273 K.N...The thermal equation of state of a natural kyanite has been investigated with a DIA-type,cubic-anvil apparatus(SAM85)combined with an energy-dispersive synchrotron X-ray radiation technique up to 8.55 GPa and 1273 K.No phase transition was observed in the studied pressure-temperature(P-T)range.The Le Bail full profile refinement technique was used to derive the unit-cell parameters.By fixing the bulk modulus K 0 as 196 GPa and its pressure derivative K0 as 4,our P-V(volume)-T data were fitted to the high temperature BircheMurnaghan equation of state.The obtained parameters for the kyanite are:V_(0)=294.05(9)Å^(3),a=2.53(11)×10^(-5)K^(-1) and(K/T)P=-0.021(8)GPa∙K^(-1).These parameters have been combined with other experimentally-measured thermodynamic data for the relevant phases to calculate the P-T locus of the reaction kyanite¼stishoviteþcorundum.With this thermodynamically constrained phase boundary,previous high-pressure phase equi-librium experimental studies with the multi-anvil press have been evaluated.展开更多
As one of the most important ways to reduce the greenhouse gas emission,carbon dioxide(CO2)enhanced gas recovery(CO2-EGR) is attractive since the gas recovery can be enhanced simultaneously with CO2sequestration.B...As one of the most important ways to reduce the greenhouse gas emission,carbon dioxide(CO2)enhanced gas recovery(CO2-EGR) is attractive since the gas recovery can be enhanced simultaneously with CO2sequestration.Based on the existing equation of state(EOS) module of TOUGH2 MP,extEOS7C is developed to calculate the phase partition of H2O-CO2-CH4-NaCl mixtures accurately with consideration of dissolved NaCI and brine properties at high pressure and temperature conditions.Verifications show that it can be applied up to the pressure of 100 MPa and temperature of 150℃.The module was implemented in the linked simulator TOUGH2MP-FLAC3 D for the coupled hydro-mechanical simulations.A simplified three-dimensional(3D)1/4 model(2.2 km×1 km×1 km) which consists of the whole reservoir,caprock and baserock was generated based on the geological conditions of a gas field in the North German Basin.The simulation results show that,under an injection rate of 200,000 t/yr and production rate of 200,000 sm3/d,CO2breakthrough occurred in the case with the initial reservoir pressure of 5 MPa but did not occur in the case of 42 MPa.Under low pressure conditions,the pressure driven horizontal transport is the dominant process;while under high pressure conditions,the density driven vertical flow is dominant.Under the considered conditions,the CO2-EGR caused only small pressure changes.The largest pore pressure increase(2 MPa) and uplift(7 mm) occurred at the caprock bottom induced by only CO2injection.The caprock had still the primary stress state and its integrity was not affected.The formation water salinity and temperature variations of ±20℃ had small influences on the CO2-EGR process.In order to slow down the breakthrough,it is suggested that CO2-EGR should be carried out before the reservoir pressure drops below the critical pressure of CO2.展开更多
In this paper five equations of state are tested for checking their ability to predict the Joule-Thomson inversion curve.These five equations of state are:Mohsennia-Modarres-Mansoori(MMM),Ji-Lemp(JL),modified Soave-Re...In this paper five equations of state are tested for checking their ability to predict the Joule-Thomson inversion curve.These five equations of state are:Mohsennia-Modarres-Mansoori(MMM),Ji-Lemp(JL),modified Soave-Redlich-Kwang(SRK)equation of state by Graboski(MSRK1),modified SRK equation of state by Peneloux and Rauzy(MSRK2),and modified Peng-Robinson (PR)equation of state by Rauzy(PRmr).The investigated equations of state give good prediction of the low-temperature branch of the inversion curve,except for MMM equation of state.The high-temperature branch and the peak of the inversion curve have been observed,in general,to be sensitive to the applied equation of state.The values of the maximum inversion temperature and maximum inversion pressure are calculated for each component used in this work.展开更多
The unreacted equation of state(EOS) of energetic materials is an important thermodynamic relationship to characterize their high pressure behaviors and has practical importance. The previous experimental and theore...The unreacted equation of state(EOS) of energetic materials is an important thermodynamic relationship to characterize their high pressure behaviors and has practical importance. The previous experimental and theoretical works on the equation of state of several energetic materials including nitromethane, 1,3,5-trinitrohexahydro-1,3,5-triazine(RDX),1,3,5,7-tetranitro-1,3,5,7-tetrazacyclooctane(HMX), hexanitrostilbene(HNS), hexanitrohexaazaisowurtzitane(HNIW or CL-20), pentaerythritol tetranitrate(PETN), 2,6-diamino-3,5-dinitropyrazine-1-oxide(LLM-105), triamino-trinitrobenzene(TATB), 1,1-diamino-2,2-dinitroethene(DADNE or FOX-7), and trinitrotoluene(TNT) are reviewed in this paper. The EOS determined from hydrostatic and non-hydrostatic compressions are discussed and compared. The theoretical results based on ab initio calculations are summarized and compared with the experimental data.展开更多
The equation of state for nuclear matter is presented within the Brueckner Hartree-Fock (BHF) scheme, by using the realistic Argonne VI8 or Bonn B two-nucleon potentials plus their corresponding microscopic three-nu...The equation of state for nuclear matter is presented within the Brueckner Hartree-Fock (BHF) scheme, by using the realistic Argonne VI8 or Bonn B two-nucleon potentials plus their corresponding microscopic three-nucleon forces. It is then applied to calculate the properties of finite nuclei within a simple liquid-drop model, and we compare the calculated volume, surface, and Coulomb parameters with the empirical ones from the liquid drop model. Nuclear density distributions and charge radii in good agreement with the experimental data are obtained~ and we predict the neutron skin thickness of various nuclei.展开更多
For further improving the representation of mixture VLE data,the local composition version of CCORequation of state has been developed and tested on 42 sets low-pressure and high-pressure as well as polarand nonpolar ...For further improving the representation of mixture VLE data,the local composition version of CCORequation of state has been developed and tested on 42 sets low-pressure and high-pressure as well as polarand nonpolar VLE data.The data reduction results were compared with conventional quadratic mixing ruleand activity coefficient method.The comparison with quadratic mixing rule showed that the local composition version significantly im-proved the data fitting of polar systems,especially for those highly nonideal mixtures where quadratic mixingrule failed to fit satisfactorily.The comparison with the well-known activity coefficient method——Hayden-O’Connell-Wilson model,indicated that this new version gave,in general,better fit to those low-pressure strongly polar systems,which traditionally has to be treated by activity coefficient approach.展开更多
Thermodynamic calculation is the theoretical basis for the study of initiation and detonation,as well as the prerequisite for forecasting the detonation performance of unknown explosives.Based on the VLWR(Virial-Wu)th...Thermodynamic calculation is the theoretical basis for the study of initiation and detonation,as well as the prerequisite for forecasting the detonation performance of unknown explosives.Based on the VLWR(Virial-Wu)thermodynamic code,this paper introduced the universal solid equation of state(EOS)VINET.In order to truly reflect the compressibility of nanocarbon under the extremely high-temperature and high-pressure environment in detonation,an SVM(support vector machine)was utilized to optimize the input parameters of carbon.The detonation performance of several explosives with different densities was calculated by the optimized universal EOS,and the results show that the thermodynamic code coupled with the universal solid EOS VINET can predict the detonation performance parameters of explosives well.To investigate the application of the thermodynamic code with the improved VINET EOS in the working capacity of explosives,the interrelationship between pressure P-particle velocity u and pressure P-volume V were computed for the detonation products of TNT and HMX-based PBX(HMX:binder:insensitive agent=95:4.3:0.7)in the CJ isentropic state.A universal curve proposed by Cooper was used to compared the computed isentropic state,where the ratio of pressure to CJ state were plotted against the ratio of velocity to CJ state.The parameters of the JWL(Jones-Wilkins-Lee)EOS for detonation products were obtained by fitting the P-V curve.The cylinder tests of TNT and HMX-based PBX were numerically simulated using the LS-DYNA,it is verified that,within a certain range,the improved algorithm has superiority in describing the working capacity of explosives.展开更多
The present work uses the concept of a scaled particle along with the perturbation and variation approach, to develop an equation of state (EOS) for a mixture of hard sphere (HS), Lennar-Jones (L J) fluids. A su...The present work uses the concept of a scaled particle along with the perturbation and variation approach, to develop an equation of state (EOS) for a mixture of hard sphere (HS), Lennar-Jones (L J) fluids. A suitable flexible functional form for the radial distribution function G(R) is assumed for the mixture, with R as a variable. The function G(R) has an arbitrary parameter m and a different equation of state can be obtained with a suitable choice of m. For m = 0.75 and m = 0.83 results are close to molecular dynamics (MD) result for pure HS and LJ fluid respectively.展开更多
This study presents high pressure phase transitions and equation of states of cerium under pressures up to 51 GPa at room temperature. The angle-dispersive x-ray diffraction experiments are carried out using a high en...This study presents high pressure phase transitions and equation of states of cerium under pressures up to 51 GPa at room temperature. The angle-dispersive x-ray diffraction experiments are carried out using a high energy synchrotron x-ray source. The bulk moduli of high pressure phases of cerium are calculated using the Birch-Mumaghan equation. We discuss and correct several previous controversial conclusions, which are caused by the measurement accuracy or personal explanation. The c/a axial ratio of e-Ce has a maximum value at about 29 GPa, i.e., c/a ≈ 1.690.展开更多
The quantum molecular dynamics based on the density functional theory has been adopted to simulate the equation of state for the shock compressed lithium. In contrary to some earlier experimental measurement and theor...The quantum molecular dynamics based on the density functional theory has been adopted to simulate the equation of state for the shock compressed lithium. In contrary to some earlier experimental measurement and theoretical simulation,there is not any evidence of the ‘kink' in the Hugoniot curve in our accurate simulation. Throughout the shock compression process, only a simple solid-to-liquid melting behavior is demonstrated, instead of complicated solid–solid phase transitions. Moreover, the x-ray absorption near-edge spectroscopy has been predicted as a feasible way to diagnose the structural evolution of warm dense lithium in this density region.展开更多
In this paper, we employ the concept of probability for creating a cavity with diameter d in fluid along with the perturbation and variation approach, and develop an equation of state (EOS) for a hard sphere (HS) ...In this paper, we employ the concept of probability for creating a cavity with diameter d in fluid along with the perturbation and variation approach, and develop an equation of state (EOS) for a hard sphere (HS) and Lennard Jones (L J) fluids. A suitable axiomatic form for surface tension S(r) is assumed for the pure fluid, with r as a variable. The function S(r) has an arbitrary parameter rn. S(r) = A + B(d/r)/[1 + m(d/r)]. We use the condition in terms of radial distribution function G(λd,η) containing the self-consistent parameter λ and the condition of continuity at r = d/2 to determine A and B. A different EOS can be obtained with a suitable choice of rn and the EOS has a lower root-mean-square deviation than that of Barker-Henderson BH2 for LJ fluids.展开更多
The thermodynamic properties of the ε phase of solid oxygen are studied by using the analytic mean field approach (AMFP). Analytic expressions for the Helmholtz free energy, internal energy and equation of state of...The thermodynamic properties of the ε phase of solid oxygen are studied by using the analytic mean field approach (AMFP). Analytic expressions for the Helmholtz free energy, internal energy and equation of state of solid oxygen have been derived based on the multi-exponential potential. The formulism for the case of double-exponential (DE) model is applied to the ε phase of solid oxygen. Its four potential parameters are determined through fitting the experimental compression data of the ε phase of solid oxygen. Numerical results of the pressure dependence of the volume calculated by using the AMFP are in good agreement with the original experimental data. This suggests that the AMFP is a useful approach to study the thermodynamic properties of the ε phase of solid oxygen. Furthermore, we predict the variation of the volume, lattice parameters and intermolecular distances with pressure, and some thermodynamic quantities versus volume, at several higher temperatures.展开更多
基金supported by the National Natural Science Foundation of China[Grant Nos.51938011 and 51908405]Australian Research Council。
文摘Adopting the classical theory of hydrocodes,the constitutive relations of concretes are separated into an equation of state(EoS)which describes the volumetric behavior of concrete material and a strength model which depicts the shear properties of concrete.The experiments on the EoS of concrete is always challenging due to the technical difficulties and equipment limitations,especially for the specimen size effect on the EoS.Although some researchers investigate the shock properties of concretes by fly-plate impact tests,the specimens used in their tests are usually in one size.In this paper,the fly-plate impact tests on concrete specimens with different sizes are performed to investigate the size effect on the shock properties of concrete materials.The mechanical background of the size effect on the shock properties are revealed,which is related to the lateral rarefaction effect and the deviatoric stress produced in the specimen.According to the tests results,the modified EoS considering the size effect on the shock properties of concrete are proposed,which the bulk modulus of concrete is unpredicted by up to 20% if size effects are not accounted for.
基金supported by the National Natural Science Foundation of China(No.11822503)。
文摘Based on the Hugenholtz-Van Hove theorem six basic quantities of the EoS in isospin asymmetric nuclear matter are expressed in terms of the nucleon kinetic energy t(k),the isospin symmetric and asymmetric parts of the single-nucleon potentials U_(0)(ρ,k)and U_(sym,i)(ρ,k).The six basic quantities include the quadratic symmetry energy E_(sym,2)(ρ),the quartic symmetry energy E_(sym,4)(ρ),their corresponding density slopes L_(2)(ρ)and L_(4)(ρ),and the incompressibility coefficients K_(2)(ρ)and K_(4)(ρ).By using four types of well-known effective nucleon-nucleon interaction models,namely the BGBD,MDI,Skyrme,and Gogny forces,the density-and isospin-dependent properties of these basic quantities are systematically calculated and their values at the saturation density q_(0)are explicitly given.The contributions to these quantities from t(k)U_(0)(ρ,k),and U_(sym,i)(ρ,k)are also analyzed at the norma nuclear density q_(0).It is clearly shown that the first-order asymmetric term U_(sym,1)(ρ,k)(also known as the symmetry potential in the Lane potential)plays a vital role in determining the density dependence of the quadratic symmetry energy E_(sym,2)(ρ).It is also shown that the contributions from the high-order asymmetric parts of the single-nucleon potentials(U_(sym,i)(ρ,k)with i>1)cannot be neglected in the calculations of the other five basic quantities Moreover,by analyzing the properties of asymmetric nuclear matter at the exact saturation densityρ_(sat)(δ),the corresponding quadratic incompressibility coefficient is found to have a simple empirical relation K_(sat,2)=K_(2)(ρ_(0))-4.14L_(2)(ρ_(0))
基金financially supported by National Natural Science Foundation of China(U1762101)National Science and Technology Major Projects(2017ZX05069)
文摘This paper presents a three-dimensional, three-phase compositional model considering CO2 phase equilibrium between water and oil. In this model, CO2 is mutually soluble in aqueous and hydrocarbon phases, while other components, except water,exist in hydrocarbon phase. The Peng–Robinson(PR) equation of state and the Wong–Sandler mixing rule with non-random two-liquid parameters are used to calculate CO2 fugacity in the aqueous phase. One-dimensional and three-dimensional CO2 flooding examples show that a significant amount of injected CO2 is dissolved in water. Our simulation shows 7% of injected CO2 can be dissolved in the aqueous phase, which delays oil recovery by 4%. The gas rate predicted by the model is smaller than the conventional model as long as water is undersaturated by CO2, which can be considered as 'lost' in the aqueous phase. The model also predicts that the delayed oil can be recovered after the gas breakthrough, indicating that delayed oil is hard to recover in field applications. A three-dimensional example reveals that a highly stratified reservoir causes uneven displacement and serious CO2 breakthrough. If mobility control measures like water alternating gas are undertaken, the solubility e ects will be more pronounced than this example.
基金supported by the National Natural Science Foundation of China(Nos.10804011,11176002).
文摘In this paper we introduce the wide regime equation of state(WEOS)developed in Institute of Applied Physics and Computational Mathematics(IAPCM).A semi-empirical model of the WEOS is given by a thermodynamically complete potential of the Helmholtz free energy which combines several theoretical models and has some adjustable parameters calibrated via some experimental and theoretical data.The validation methods of the equation of state in wide regime are presented using copper as a prototype.The results of the WEOS are well consistent with the available theoretical and experimental data,including ab initio cold curve under compression,isotherm,Hugoniot,off-Hugoniot and sound velocity data.It enhances our confidence in the accuracy of the WEOS,which is very important for the validation and verification of equation of state in high temperature and pressure technology.
基金Zhejiang Provincial Natural Science Foundation of China!(No. 298013)
文摘In this paper, the Martin-Hou equation of state is derived by using a power series representation of radial distribution function and an analytic representation of multi-section potential based on the Barker-Henderson hard-particle perturbation theory including high-order terms. In the derivation, a theoretical form of Martin-Hou equation was obtained. It had a similar form and the same capability to predict P-V-T properties as the Martin-Hou equation and no additional data were required for evaluating the constants. The characteristic constants of the theoretical expression have certain relationships with the molecular parameters.
基金Project supported by the National Natural Science Foundation of China (Grant Nos 90816013 and 10572083)Shanghai Leading Academic Discipline Project,China (Grant No Y0103)
文摘The Galilean invariance and the induced thermo-hydrodynamics of the lattice Boltzmann Bhatnagar-Gross-Krook model are proposed together with their rigorous theoretical background. From the viewpoint of group invariance, recovering the Galilean invariance for the isothermal lattice Boltzmann Bhatnagar-Gross-Krook equation (LBGKE) induces a new natural thermal-dynamical system, which is compatible with the elementary statistical thermodynamics.
文摘Within the approximations of spherical lattice cell, central-field, and relativistic Fermi statis- tics, an algorithm with average atom model is presented to calculate the electronic energy levels and equation of state for hot and dense matter at arbitrary densities and temperatures. Choosing Zink's analytical potential as initial potential, we have solved the Dirac-Slater equation which satisfies the Weigner-Seitz boundary condition. The electronic energy bands are not taken into account. Tak- ing energy level degeneracy as a continuous function of density, we have considered the pressure ionization effects for highly dense matter. Results for ^(13)Al atom are shown.
文摘Based on the classical Roe method, we develop an interface capture method according to the general equation of state, and extend the single-fluid Roe method to the two-dimensional (2D) multi-fluid flows, as well as construct the continuous Roe matrix for the whole flow field. The interface capture equations and fluid dynamic conservative equations are coupled together and solved by using any high-resolution schemes that usually suit for the single-fluid flows. Some numerical examples are given to illustrate the solution of 1D and 2D multi-fluid Riemann problems.
基金This work is financially supported by the Natural Science Foundation of China(Grant No.41273072)by the Strategic Priority Research Program(B)of the Chinese Academy of Sciences(XDB18030602).
文摘The thermal equation of state of a natural kyanite has been investigated with a DIA-type,cubic-anvil apparatus(SAM85)combined with an energy-dispersive synchrotron X-ray radiation technique up to 8.55 GPa and 1273 K.No phase transition was observed in the studied pressure-temperature(P-T)range.The Le Bail full profile refinement technique was used to derive the unit-cell parameters.By fixing the bulk modulus K 0 as 196 GPa and its pressure derivative K0 as 4,our P-V(volume)-T data were fitted to the high temperature BircheMurnaghan equation of state.The obtained parameters for the kyanite are:V_(0)=294.05(9)Å^(3),a=2.53(11)×10^(-5)K^(-1) and(K/T)P=-0.021(8)GPa∙K^(-1).These parameters have been combined with other experimentally-measured thermodynamic data for the relevant phases to calculate the P-T locus of the reaction kyanite¼stishoviteþcorundum.With this thermodynamically constrained phase boundary,previous high-pressure phase equi-librium experimental studies with the multi-anvil press have been evaluated.
基金funded by the National Natural Science Foundation of China(Grant No.NSFC51374147)the German Society for Petroleum and Coal Science and Technology(Grant No.DGMK680-4)
文摘As one of the most important ways to reduce the greenhouse gas emission,carbon dioxide(CO2)enhanced gas recovery(CO2-EGR) is attractive since the gas recovery can be enhanced simultaneously with CO2sequestration.Based on the existing equation of state(EOS) module of TOUGH2 MP,extEOS7C is developed to calculate the phase partition of H2O-CO2-CH4-NaCl mixtures accurately with consideration of dissolved NaCI and brine properties at high pressure and temperature conditions.Verifications show that it can be applied up to the pressure of 100 MPa and temperature of 150℃.The module was implemented in the linked simulator TOUGH2MP-FLAC3 D for the coupled hydro-mechanical simulations.A simplified three-dimensional(3D)1/4 model(2.2 km×1 km×1 km) which consists of the whole reservoir,caprock and baserock was generated based on the geological conditions of a gas field in the North German Basin.The simulation results show that,under an injection rate of 200,000 t/yr and production rate of 200,000 sm3/d,CO2breakthrough occurred in the case with the initial reservoir pressure of 5 MPa but did not occur in the case of 42 MPa.Under low pressure conditions,the pressure driven horizontal transport is the dominant process;while under high pressure conditions,the density driven vertical flow is dominant.Under the considered conditions,the CO2-EGR caused only small pressure changes.The largest pore pressure increase(2 MPa) and uplift(7 mm) occurred at the caprock bottom induced by only CO2injection.The caprock had still the primary stress state and its integrity was not affected.The formation water salinity and temperature variations of ±20℃ had small influences on the CO2-EGR process.In order to slow down the breakthrough,it is suggested that CO2-EGR should be carried out before the reservoir pressure drops below the critical pressure of CO2.
文摘In this paper five equations of state are tested for checking their ability to predict the Joule-Thomson inversion curve.These five equations of state are:Mohsennia-Modarres-Mansoori(MMM),Ji-Lemp(JL),modified Soave-Redlich-Kwang(SRK)equation of state by Graboski(MSRK1),modified SRK equation of state by Peneloux and Rauzy(MSRK2),and modified Peng-Robinson (PR)equation of state by Rauzy(PRmr).The investigated equations of state give good prediction of the low-temperature branch of the inversion curve,except for MMM equation of state.The high-temperature branch and the peak of the inversion curve have been observed,in general,to be sensitive to the applied equation of state.The values of the maximum inversion temperature and maximum inversion pressure are calculated for each component used in this work.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.11174045 and 11404050)
文摘The unreacted equation of state(EOS) of energetic materials is an important thermodynamic relationship to characterize their high pressure behaviors and has practical importance. The previous experimental and theoretical works on the equation of state of several energetic materials including nitromethane, 1,3,5-trinitrohexahydro-1,3,5-triazine(RDX),1,3,5,7-tetranitro-1,3,5,7-tetrazacyclooctane(HMX), hexanitrostilbene(HNS), hexanitrohexaazaisowurtzitane(HNIW or CL-20), pentaerythritol tetranitrate(PETN), 2,6-diamino-3,5-dinitropyrazine-1-oxide(LLM-105), triamino-trinitrobenzene(TATB), 1,1-diamino-2,2-dinitroethene(DADNE or FOX-7), and trinitrotoluene(TNT) are reviewed in this paper. The EOS determined from hydrostatic and non-hydrostatic compressions are discussed and compared. The theoretical results based on ab initio calculations are summarized and compared with the experimental data.
基金Supported by the National Natural Science Foundation of China under Grant Nos 11075037 and 11475045the Scientific Research Foundation for the Returned Overseas Chinese Scholars of the Ministry of Education of China+2 种基金the Fundamental Research Funds for the Central Universities of Chinathe Shanghai Leading Academic Discipline Project under Grant No B107the 'NewCompStar',COST Action MP1304
文摘The equation of state for nuclear matter is presented within the Brueckner Hartree-Fock (BHF) scheme, by using the realistic Argonne VI8 or Bonn B two-nucleon potentials plus their corresponding microscopic three-nucleon forces. It is then applied to calculate the properties of finite nuclei within a simple liquid-drop model, and we compare the calculated volume, surface, and Coulomb parameters with the empirical ones from the liquid drop model. Nuclear density distributions and charge radii in good agreement with the experimental data are obtained~ and we predict the neutron skin thickness of various nuclei.
文摘For further improving the representation of mixture VLE data,the local composition version of CCORequation of state has been developed and tested on 42 sets low-pressure and high-pressure as well as polarand nonpolar VLE data.The data reduction results were compared with conventional quadratic mixing ruleand activity coefficient method.The comparison with quadratic mixing rule showed that the local composition version significantly im-proved the data fitting of polar systems,especially for those highly nonideal mixtures where quadratic mixingrule failed to fit satisfactorily.The comparison with the well-known activity coefficient method——Hayden-O’Connell-Wilson model,indicated that this new version gave,in general,better fit to those low-pressure strongly polar systems,which traditionally has to be treated by activity coefficient approach.
基金the financial support from the National Natural Science Foundation of China(No.11902298)the National Key Research and Development Program of China(No.2017YFC0804701)。
文摘Thermodynamic calculation is the theoretical basis for the study of initiation and detonation,as well as the prerequisite for forecasting the detonation performance of unknown explosives.Based on the VLWR(Virial-Wu)thermodynamic code,this paper introduced the universal solid equation of state(EOS)VINET.In order to truly reflect the compressibility of nanocarbon under the extremely high-temperature and high-pressure environment in detonation,an SVM(support vector machine)was utilized to optimize the input parameters of carbon.The detonation performance of several explosives with different densities was calculated by the optimized universal EOS,and the results show that the thermodynamic code coupled with the universal solid EOS VINET can predict the detonation performance parameters of explosives well.To investigate the application of the thermodynamic code with the improved VINET EOS in the working capacity of explosives,the interrelationship between pressure P-particle velocity u and pressure P-volume V were computed for the detonation products of TNT and HMX-based PBX(HMX:binder:insensitive agent=95:4.3:0.7)in the CJ isentropic state.A universal curve proposed by Cooper was used to compared the computed isentropic state,where the ratio of pressure to CJ state were plotted against the ratio of velocity to CJ state.The parameters of the JWL(Jones-Wilkins-Lee)EOS for detonation products were obtained by fitting the P-V curve.The cylinder tests of TNT and HMX-based PBX were numerically simulated using the LS-DYNA,it is verified that,within a certain range,the improved algorithm has superiority in describing the working capacity of explosives.
文摘The present work uses the concept of a scaled particle along with the perturbation and variation approach, to develop an equation of state (EOS) for a mixture of hard sphere (HS), Lennar-Jones (L J) fluids. A suitable flexible functional form for the radial distribution function G(R) is assumed for the mixture, with R as a variable. The function G(R) has an arbitrary parameter m and a different equation of state can be obtained with a suitable choice of m. For m = 0.75 and m = 0.83 results are close to molecular dynamics (MD) result for pure HS and LJ fluid respectively.
基金supported by the National Natural Science Foundation of China(Grant No.NSAF.U1330115)the National Major Scientific Instrument and Equipment Development Project of China(Grant No.2012YQ130234)
文摘This study presents high pressure phase transitions and equation of states of cerium under pressures up to 51 GPa at room temperature. The angle-dispersive x-ray diffraction experiments are carried out using a high energy synchrotron x-ray source. The bulk moduli of high pressure phases of cerium are calculated using the Birch-Mumaghan equation. We discuss and correct several previous controversial conclusions, which are caused by the measurement accuracy or personal explanation. The c/a axial ratio of e-Ce has a maximum value at about 29 GPa, i.e., c/a ≈ 1.690.
基金supported by the National Natural Science Foundation of China(Grant Nos.11474034 and 11675024)the Foundation for Development of Science and Technology of China Academy of Engineering Physics(Grant Nos.2015B0102020 and 2015B0102022)the Science Challenge Project(Grant No.TZ2016005)
文摘The quantum molecular dynamics based on the density functional theory has been adopted to simulate the equation of state for the shock compressed lithium. In contrary to some earlier experimental measurement and theoretical simulation,there is not any evidence of the ‘kink' in the Hugoniot curve in our accurate simulation. Throughout the shock compression process, only a simple solid-to-liquid melting behavior is demonstrated, instead of complicated solid–solid phase transitions. Moreover, the x-ray absorption near-edge spectroscopy has been predicted as a feasible way to diagnose the structural evolution of warm dense lithium in this density region.
文摘In this paper, we employ the concept of probability for creating a cavity with diameter d in fluid along with the perturbation and variation approach, and develop an equation of state (EOS) for a hard sphere (HS) and Lennard Jones (L J) fluids. A suitable axiomatic form for surface tension S(r) is assumed for the pure fluid, with r as a variable. The function S(r) has an arbitrary parameter rn. S(r) = A + B(d/r)/[1 + m(d/r)]. We use the condition in terms of radial distribution function G(λd,η) containing the self-consistent parameter λ and the condition of continuity at r = d/2 to determine A and B. A different EOS can be obtained with a suitable choice of rn and the EOS has a lower root-mean-square deviation than that of Barker-Henderson BH2 for LJ fluids.
基金supported by the Joint Foundation of National Natural Science Foundation of China and China Academy of Engineering Physics (Grant No 10476007)the Program for New Century Excellent Talents in University (Grant No NCET-05-0799)the Program for Excellent Talents of University of Electronic Science and Technology (Grant No 23601008)
文摘The thermodynamic properties of the ε phase of solid oxygen are studied by using the analytic mean field approach (AMFP). Analytic expressions for the Helmholtz free energy, internal energy and equation of state of solid oxygen have been derived based on the multi-exponential potential. The formulism for the case of double-exponential (DE) model is applied to the ε phase of solid oxygen. Its four potential parameters are determined through fitting the experimental compression data of the ε phase of solid oxygen. Numerical results of the pressure dependence of the volume calculated by using the AMFP are in good agreement with the original experimental data. This suggests that the AMFP is a useful approach to study the thermodynamic properties of the ε phase of solid oxygen. Furthermore, we predict the variation of the volume, lattice parameters and intermolecular distances with pressure, and some thermodynamic quantities versus volume, at several higher temperatures.