Experimental study on the size effect on the equation of state of concretes under shock loading

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.

Basic quantities of the equation of state in isospin asymmetric nuclear matter

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))

A compositional model for CO_(2) ooding including CO_(2) equilibria between water and oil using the Peng–Robinson equation of state with the Wong–Sandler mixing rule

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.

Validation for equation of state in wide regime:Copper as prototype

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.

Derivation of Martin-Hou Equation of State from Hard-particle Perturbation Theory

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.

Thermal equation of state for lattice Boltzmann gases

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.

EQUATION OF STATE CALCULATIONS FOR HOT, DENSE MATTER AT ARBITRARY DENSITIES AND TEMPERATURES

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.

Numerical method of the Riemann problem for two-dimensional multi-fluid flows with general equation of state

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.

Thermal equation of state of a natural kyanite up to 8.55 GPa and 1273 K

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.

Coupled thermo-hydro-mechanical simulation of CO2 enhanced gas recovery with an extended equation of state module for TOUGH2MP-FLAC3D

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.

Joule-Thomson inversion curve prediction by using equation of state

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.

Unreacted equation of states of typical energetic materials under static compression:A review

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 Brueckner-Hartree-Fock Equation of State for Nuclear Matter and Neutron Skin

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.

CUBIC CHAIN-OF-ROTATORS EQUATION OF STATE——LOCAL COMPOSITION VERSION

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.

An insight into the role of the association equations of states in gas hydrate modeling:a review

Encouraged by the wide spectrum of novel applications of gas hydrates,e.g.,energy recovery,gas separation,gas storage,gas transportation,water desalination,and hydrogen hydrate as a green energy resource,as well as CO2 capturing,many scientists have focused their attention on investigating this important phenomenon.Of course,from an engineering viewpoint,the mathematical modeling of gas hydrates is of paramount importance,as anticipation of gas hydrate stability conditions is effective in the design and control of industrial processes.Overall,the thermodynamic modeling of gas hydrate can be tackled as an equilibration of three phases,i.e.,liquid,gas,and solid hydrate.The inseparable component in all hydrate systems,water,is highly polar and non-ideal,necessitating the use of more advanced equation of states(EoSs) that take into account more intermolecular forces for thermodynamic modeling of these systems.Motivated by the ever-increasing number of publications on this topic,this study aims to review the application of associating EoSs for the thermodynamic modeling of gas hydrates.Three most important hydrate-based models available in the literature including the van der Waals-Platteeuw(vdW-P) model,Chen-Guo model,and Klauda-Sandler model coupled with and SAFT EoSs were investigated and compared with cubic EoSs.It was concluded that the CPA and SAFT EoSs gave very accurate results for hydrate systems as they take into account the association interactions,which are very crucial in gas hydrate systems in which water,methanol,glycols,and other types of associating compounds are available.Moreover,it was concluded that the CPA EoS is easier to use than the SAFT-type EoSs and our suggestion for the gas hydrate systems is the CPA EoS.

Calculating detonation performance of explosives by VLWR thermodynamics code introduced with universal VINET equation of state

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.

Perturbation and variational approach for the equation of state for hard-sphere and Lennard-Jones fluids

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.

New insight into prediction of phase behavior of natural gas hydrate by different cubic equations of state coupled with various mixing rules

Progress in hydrate thermodynamic study necessitates robust and fast models to be incorporated in reservoir simulation softwares. However, numerous models presented in the literature makes selection of the best,proper predictive model a cumbersome task. It is of industrial interest to make use of cubic equations of state(EOS) for modeling hydrate equilibria. In this regard, this study focuses on evaluation of three common EOSs including Peng–Robinson, Soave–Redlich–Kwong and Valderrama–Patel–Teja coupled with van der Waals and Platteeuw theory to predict hydrate P–T equilibrium of a real natural gas sample. Each EOS was accompanied with three mixing rules, including van der Waals(vd W),Avlonitis non-density dependent(ANDD) and general nonquadratic(GNQ). The prediction of cubic EOSs was in sufficient agreement with experimental data and with overall AARD% of less than unity. In addition, PR plus ANDD proved to be the most accurate model in this study for prediction of hydrate equilibria with AARD% of 0.166.It was observed that the accuracy of cubic EOSs studied in this paper depends on mixing rule coupled with them,especially at high-pressure conditions. Lastly, the present study does not include any adjustable parameter to be correlated with hydrate phase equilibrium data.

Structure phase transformation and equation of state of cerium metal under pressures up to 51 GPa

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.

Equation of state for warm dense lithium： A first principles investigation

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.