Rubber-like materials that are commonly used in structural applications are modelled using hyperelastic material models.Most of the hyperelastic materials are nearly incompressible,which poses challenges,i.e.,volumetr...Rubber-like materials that are commonly used in structural applications are modelled using hyperelastic material models.Most of the hyperelastic materials are nearly incompressible,which poses challenges,i.e.,volumetric locking during numerical modelling.There exist many formulations in the context of the finite element method,among which the mixed displacementpressure formulation is robust.However,such a displacement-pressure formulation is less explored in meshfree methods,which mitigates the problem associated with mesh distortion during large deformation.This work addresses this issue of alleviating volumetric locking in the element-free Galerkin method(EFGM),which is one of the popular meshfree methods.A two-field mixed variational formulation using the perturbed Lagrangian approach within the EFGM framework is proposed for modelling nearly incompressible hyperelastic material models,such as Neo-Hookean and Mooney-Rivlin.Taking advantage of the meshless nature of the EFGM,this work introduces a unique approach by randomly distributing pressure nodes across the geometry,following specific guidelines.A wide spectrum of problems involving bending,tension,compression,and contact is solved using two approaches of the proposed displacement-pressure node formulation involving regular and irregular pressure node distribution.It is observed that both approaches give accurate results compared to the reference results,though the latter offers flexibility in the pressure nodal distribution.展开更多
Natural gas hydrate(NGH),considered as a type of premium energy alternative to conventional hydrocarbons,has been broadly studied.The estimate of the total NGH resources in the world has decreased by more than 90%sinc...Natural gas hydrate(NGH),considered as a type of premium energy alternative to conventional hydrocarbons,has been broadly studied.The estimate of the total NGH resources in the world has decreased by more than 90%since the first evaluation in 1973.Geographic and geophysical conditions of the South China Sea(SCS)are favorable for the formation of NGH,which has been proved by drilling results up to date.The recoverability of the NGH in the SCS has been confirmed by the production tests using both vertical and horizontal wells.Since 2001,35 estimates of NGH resources in the SCS have been made,with relatively stable results varying between 600 and 900×109 ton oil equivalent.In these estimations,the volumetric method was commonly adopted,but the geological conditions,the migration-accumulation mechanisms of NGH,and the practical recoverability were not considered.These estimates cannot be regarded as evaluated resources according to the international resource evaluation standards,but are at most about prospective gas content of NGH,thus inefficient for guiding explorations and developments.To solve these problems,this study divides the past NGH surveys in the SCS into seven stages,acquires key geological parameters of every stage based on previous studies and analogy with other areas,evaluates the NGH resources of these seven stages by using the volumetric method,then adopts a new trend-analysis method to simulate the downward trend of these estimates,and finally predicts the NGH resources in the SCS at 2025 and 2030.The downward trend is because of the continuous improvement of NGH understanding over time,which is consistent with the trend of global NGH estimates.At the present stage(from 2019 to 2021),the average technically recoverable resource(ATRR)is 7.0×10^(12)m^(3),and the estimates of 2025 and 2030 ATRR are 6.46×10^(12)m^(3) and 4.01×10^(12)m^(3)respectively,with a difference of less than 40%.Therefore,it can be inferred that the ATRR of NGH in the SCS is between 4.0 and 6.5×10^(12)m^(3),with an average of 5.25×10^(12)m^(3).展开更多
Porous materials present significant advantages for absorbing radioactive isotopes in nuclear waste streams.To improve absorption efficiency in nuclear waste treatment,a thorough understanding of the diffusion-advecti...Porous materials present significant advantages for absorbing radioactive isotopes in nuclear waste streams.To improve absorption efficiency in nuclear waste treatment,a thorough understanding of the diffusion-advection process within porous structures is essential for material design.In this study,we present advancements in the volumetric lattice Boltzmann method(VLBM)for modeling and simulating pore-scale diffusion-advection of radioactive isotopes within geopolymer porous structures.These structures are created using the phase field method(PFM)to precisely control pore architectures.In our VLBM approach,we introduce a concentration field of an isotope seamlessly coupled with the velocity field and solve it by the time evolution of its particle population function.To address the computational intensity inherent in the coupled lattice Boltzmann equations for velocity and concentration fields,we implement graphics processing unit(GPU)parallelization.Validation of the developed model involves examining the flow and diffusion fields in porous structures.Remarkably,good agreement is observed for both the velocity field from VLBM and multiphysics object-oriented simulation environment(MOOSE),and the concentration field from VLBM and the finite difference method(FDM).Furthermore,we investigate the effects of background flow,species diffusivity,and porosity on the diffusion-advection behavior by varying the background flow velocity,diffusion coefficient,and pore volume fraction,respectively.Notably,all three parameters exert an influence on the diffusion-advection process.Increased background flow and diffusivity markedly accelerate the process due to increased advection intensity and enhanced diffusion capability,respectively.Conversely,increasing the porosity has a less significant effect,causing a slight slowdown of the diffusion-advection process due to the expanded pore volume.This comprehensive parametric study provides valuable insights into the kinetics of isotope uptake in porous structures,facilitating the development of porous materials for nuclear waste treatment applications.展开更多
This standard specifies the method summary, reagents, apparatus, sampling, procedure, test results calculation and permissible tolerance of the determination of magnesium oxide by CyDTA volumetric method.
Formulation and numerical evaluation of a novel twice-interpolation finite element method （TFEM） is presented for solid mechanics problems. In this method, the trial function for Galerkin weak form is constructed th...Formulation and numerical evaluation of a novel twice-interpolation finite element method （TFEM） is presented for solid mechanics problems. In this method, the trial function for Galerkin weak form is constructed through two stages of consecutive interpolation. The primary interpolation follows exactly the same procedure of standard FEM and is further reproduced according to both nodal values and averaged nodal gradients obtained from primary interpolation. The trial functions thus constructed have continuous nodal gradients and contain higher order polynomial without increasing total freedoms. Several benchmark examples and a real dam problem are used to examine the TFEM in terms of accuracy and convergence. Compared with standard FEM, TFEM can achieve significantly better accuracy and higher convergence rate, and the continuous nodal stress can be obtained without any smoothing operation. It is also found that TFEM is insensitive to the quality of the elemental mesh. In addition, the present TFEM can treat the incompressible material without any modification.展开更多
This paper presents a quasi-static implicit generalized interpolation material point method(i GIMP)with B-bar approach for large deformation geotechnical problems.The i GIMP algorithm is an extension of the implicit m...This paper presents a quasi-static implicit generalized interpolation material point method(i GIMP)with B-bar approach for large deformation geotechnical problems.The i GIMP algorithm is an extension of the implicit material point method(iMPM).The global stiffness matrix is formed explicitly and the Newton-Raphson iterative method is used to solve the equilibrium equations.Where possible,the implementation procedure closely follows standard finite element method(FEM)approaches to allow easy conversion of other FEM codes.The generalized interpolation function is assigned to eliminate the inherent cell crossing noise within conventional MPM.For the first time,the B-bar approach is used to overcome volumetric locking in standard GIMP method for near-incompressible non-linear geomechanics.The proposed i GIMP was tested and compared with i MPM and analytical solutions via a 1 D column compression problem.Results highlighted the superiority of the i GIMP approach in reducing stress oscillations,thereby improving computational accuracy.Then,elasto-plastic slope stabilities and rigid footing problems were considered,further illustrating the ability of the proposed method to overcome volumetric locking due to incompressibility.Results showed that the proposed i GIMP with B-bar approach can be used to simulate geotechnical problems with large deformations.展开更多
The development of a general discrete element method for irregularly shaped particles is the core issue of the simulation of the dynamic behavior of granular materials.The general energy-conserving contact theory is u...The development of a general discrete element method for irregularly shaped particles is the core issue of the simulation of the dynamic behavior of granular materials.The general energy-conserving contact theory is used to establish a universal discrete element method suitable for particle contact of arbitrary shape.In this study,three dimentional(3D)modeling and scanning techniques are used to obtain a triangular mesh representation of the true particles containing typical concave particles.The contact volumebased energy-conserving model is used to realize the contact detection between irregularly shaped particles,and the contact force model is refined and modified to describe the contact under real conditions.The inelastic collision processes between the particles and boundaries are simulated to verify the robustness of the modified contact force model and its applicability to the multi-point contact mode.In addition,the packing process and the flow process of a large number of irregular particles are simulated with the modified discrete element method(DEM)to illustrate the applicability of the method of complex problems.展开更多
Vibration fatigue is one of the main failure modes of blade.The vibration fatigue life of blade is scattered caused by manufacture error,material property dispersion and external excitation randomness.A new vibration ...Vibration fatigue is one of the main failure modes of blade.The vibration fatigue life of blade is scattered caused by manufacture error,material property dispersion and external excitation randomness.A new vibration fatigue probabilistic life prediction model(VFPLPM)and a prediction method are proposed in this paper.Firstly,as one-dimensional volumetric method(ODVM)only considers the principle calculation direction,a three-dimensional space vector volumetric method(TSVVM)is proposed to improve fatigue life prediction accuracy for actual threedimensional engineering structure.Secondly,based on the two volumetric methods(ODVM and TSVVM),the material C-P-S-N fatigue curve model(CFCM)and the maximum entropy quantile function model(MEQFM),VFPLPM is established to predict the vibration fatigue probabilistic life of blade.The VFPLPM is combined with maximum stress method(MSM),ODVM and TSVVM to estimate vibration fatigue probabilistic life of blade simulator by finite element simulation,and is verified by vibration fatigue test.The results show that all of the three methods can predict the vibration fatigue probabilistic life of blade simulator well.VFPLPM &TSVVM method has the highest computational accuracy for considering stress gradient effect not only in the principle calculation direction but also in other space vector directions.展开更多
Measurement of the volume of gas adsorbed per unit mass of coal with increasing pressure at a constant temperature produces an isotherm that describes the gas storage capacity of this type of coal. The accurate testin...Measurement of the volume of gas adsorbed per unit mass of coal with increasing pressure at a constant temperature produces an isotherm that describes the gas storage capacity of this type of coal. The accurate testing and interpretation of coal sorption isotherm plays an important role in the areas of coal mine methane drainage, coalbed methane （CBM） reservoir resource assessment, enhanced coalbed methane （ECBM） recovery, as well as the carbon dioxide （CO2） sequestration in deep coal seams or similar geological formations. Different coal sorption isotherm testing apparatus and associated calculation methods are critically reviewed and presented in this paper. These include both volumetric and gravimetric based methods, as well as experimental sorption tests with confining stress and direction sorption methods. The volumetric techniques utilise experimental apparatus with sample cell and injection pump and that with both sample cell and reference cell. Whilst the gravimetric approachesinclude methods with sample cell and suspension magnetic balance and that with both sample cell and reference cell. Different testing methods are compared and discussed in this study. A unique in-house-built coal sorption isotherm testing apparatus at the University of Wollongong was presented together with the calculation method, procedures and experimental results. The isotherm results can be calculated by both Soave-Redlich-Kwong （SRK） equation and calibration cure methods which can be used directly to convert the volume of adsorbed gas in different test conditions to standard condition （NTP）.展开更多
基金supported by the DST-SERB and VSSC,ISRO of the project titled“Functionality Enhancement through Design and Development of Advanced Finite Element Algorithms for STR tools”under IMPRINT.IIC(IMP/2019/000276)scheme.
文摘Rubber-like materials that are commonly used in structural applications are modelled using hyperelastic material models.Most of the hyperelastic materials are nearly incompressible,which poses challenges,i.e.,volumetric locking during numerical modelling.There exist many formulations in the context of the finite element method,among which the mixed displacementpressure formulation is robust.However,such a displacement-pressure formulation is less explored in meshfree methods,which mitigates the problem associated with mesh distortion during large deformation.This work addresses this issue of alleviating volumetric locking in the element-free Galerkin method(EFGM),which is one of the popular meshfree methods.A two-field mixed variational formulation using the perturbed Lagrangian approach within the EFGM framework is proposed for modelling nearly incompressible hyperelastic material models,such as Neo-Hookean and Mooney-Rivlin.Taking advantage of the meshless nature of the EFGM,this work introduces a unique approach by randomly distributing pressure nodes across the geometry,following specific guidelines.A wide spectrum of problems involving bending,tension,compression,and contact is solved using two approaches of the proposed displacement-pressure node formulation involving regular and irregular pressure node distribution.It is observed that both approaches give accurate results compared to the reference results,though the latter offers flexibility in the pressure nodal distribution.
基金financially supported by the CAS consultation project“South China Sea Oil and Gas Comprehensive Development Strategy”(2019-ZW11-Z-035)the National Basic Research Program of China(973 Program)(2006CB202300,2011CB201100)the National High-Tech R&D Program of China(863 Program)(2013AA092600)。
文摘Natural gas hydrate(NGH),considered as a type of premium energy alternative to conventional hydrocarbons,has been broadly studied.The estimate of the total NGH resources in the world has decreased by more than 90%since the first evaluation in 1973.Geographic and geophysical conditions of the South China Sea(SCS)are favorable for the formation of NGH,which has been proved by drilling results up to date.The recoverability of the NGH in the SCS has been confirmed by the production tests using both vertical and horizontal wells.Since 2001,35 estimates of NGH resources in the SCS have been made,with relatively stable results varying between 600 and 900×109 ton oil equivalent.In these estimations,the volumetric method was commonly adopted,but the geological conditions,the migration-accumulation mechanisms of NGH,and the practical recoverability were not considered.These estimates cannot be regarded as evaluated resources according to the international resource evaluation standards,but are at most about prospective gas content of NGH,thus inefficient for guiding explorations and developments.To solve these problems,this study divides the past NGH surveys in the SCS into seven stages,acquires key geological parameters of every stage based on previous studies and analogy with other areas,evaluates the NGH resources of these seven stages by using the volumetric method,then adopts a new trend-analysis method to simulate the downward trend of these estimates,and finally predicts the NGH resources in the SCS at 2025 and 2030.The downward trend is because of the continuous improvement of NGH understanding over time,which is consistent with the trend of global NGH estimates.At the present stage(from 2019 to 2021),the average technically recoverable resource(ATRR)is 7.0×10^(12)m^(3),and the estimates of 2025 and 2030 ATRR are 6.46×10^(12)m^(3) and 4.01×10^(12)m^(3)respectively,with a difference of less than 40%.Therefore,it can be inferred that the ATRR of NGH in the SCS is between 4.0 and 6.5×10^(12)m^(3),with an average of 5.25×10^(12)m^(3).
基金supported as part of the Center for Hierarchical Waste Form Materials,an Energy Frontier Research Center funded by the U.S.Department of Energy,Office of Science,Basic Energy Sciences under Award No.DE-SC0016574.
文摘Porous materials present significant advantages for absorbing radioactive isotopes in nuclear waste streams.To improve absorption efficiency in nuclear waste treatment,a thorough understanding of the diffusion-advection process within porous structures is essential for material design.In this study,we present advancements in the volumetric lattice Boltzmann method(VLBM)for modeling and simulating pore-scale diffusion-advection of radioactive isotopes within geopolymer porous structures.These structures are created using the phase field method(PFM)to precisely control pore architectures.In our VLBM approach,we introduce a concentration field of an isotope seamlessly coupled with the velocity field and solve it by the time evolution of its particle population function.To address the computational intensity inherent in the coupled lattice Boltzmann equations for velocity and concentration fields,we implement graphics processing unit(GPU)parallelization.Validation of the developed model involves examining the flow and diffusion fields in porous structures.Remarkably,good agreement is observed for both the velocity field from VLBM and multiphysics object-oriented simulation environment(MOOSE),and the concentration field from VLBM and the finite difference method(FDM).Furthermore,we investigate the effects of background flow,species diffusivity,and porosity on the diffusion-advection behavior by varying the background flow velocity,diffusion coefficient,and pore volume fraction,respectively.Notably,all three parameters exert an influence on the diffusion-advection process.Increased background flow and diffusivity markedly accelerate the process due to increased advection intensity and enhanced diffusion capability,respectively.Conversely,increasing the porosity has a less significant effect,causing a slight slowdown of the diffusion-advection process due to the expanded pore volume.This comprehensive parametric study provides valuable insights into the kinetics of isotope uptake in porous structures,facilitating the development of porous materials for nuclear waste treatment applications.
文摘This standard specifies the method summary, reagents, apparatus, sampling, procedure, test results calculation and permissible tolerance of the determination of magnesium oxide by CyDTA volumetric method.
基金supported by the National Natural Science Foundation of China(50474053,50475134 and 50675081)the 863 project (2007AA042142)
文摘Formulation and numerical evaluation of a novel twice-interpolation finite element method （TFEM） is presented for solid mechanics problems. In this method, the trial function for Galerkin weak form is constructed through two stages of consecutive interpolation. The primary interpolation follows exactly the same procedure of standard FEM and is further reproduced according to both nodal values and averaged nodal gradients obtained from primary interpolation. The trial functions thus constructed have continuous nodal gradients and contain higher order polynomial without increasing total freedoms. Several benchmark examples and a real dam problem are used to examine the TFEM in terms of accuracy and convergence. Compared with standard FEM, TFEM can achieve significantly better accuracy and higher convergence rate, and the continuous nodal stress can be obtained without any smoothing operation. It is also found that TFEM is insensitive to the quality of the elemental mesh. In addition, the present TFEM can treat the incompressible material without any modification.
基金the National Natural Science Foundation of China(Nos.41807223 and 51908175)the Fundamental Research Funds for the Central Universities(No.B210202096)+1 种基金the Natural Science Foundation of Guangdong Province(No.2018A030310346)the Water Conservancy Science and Technology Innovation Project of Guangdong Province(No.2020-11),China。
文摘This paper presents a quasi-static implicit generalized interpolation material point method(i GIMP)with B-bar approach for large deformation geotechnical problems.The i GIMP algorithm is an extension of the implicit material point method(iMPM).The global stiffness matrix is formed explicitly and the Newton-Raphson iterative method is used to solve the equilibrium equations.Where possible,the implementation procedure closely follows standard finite element method(FEM)approaches to allow easy conversion of other FEM codes.The generalized interpolation function is assigned to eliminate the inherent cell crossing noise within conventional MPM.For the first time,the B-bar approach is used to overcome volumetric locking in standard GIMP method for near-incompressible non-linear geomechanics.The proposed i GIMP was tested and compared with i MPM and analytical solutions via a 1 D column compression problem.Results highlighted the superiority of the i GIMP approach in reducing stress oscillations,thereby improving computational accuracy.Then,elasto-plastic slope stabilities and rigid footing problems were considered,further illustrating the ability of the proposed method to overcome volumetric locking due to incompressibility.Results showed that the proposed i GIMP with B-bar approach can be used to simulate geotechnical problems with large deformations.
基金supported by the National Key Research and Development Program of China(2018YFA0605902)the National Natural Science Foundation of China(42176241 and 11872136).
文摘The development of a general discrete element method for irregularly shaped particles is the core issue of the simulation of the dynamic behavior of granular materials.The general energy-conserving contact theory is used to establish a universal discrete element method suitable for particle contact of arbitrary shape.In this study,three dimentional(3D)modeling and scanning techniques are used to obtain a triangular mesh representation of the true particles containing typical concave particles.The contact volumebased energy-conserving model is used to realize the contact detection between irregularly shaped particles,and the contact force model is refined and modified to describe the contact under real conditions.The inelastic collision processes between the particles and boundaries are simulated to verify the robustness of the modified contact force model and its applicability to the multi-point contact mode.In addition,the packing process and the flow process of a large number of irregular particles are simulated with the modified discrete element method(DEM)to illustrate the applicability of the method of complex problems.
基金supported by the Aviation Science Foundation of China(No.20150252003)
文摘Vibration fatigue is one of the main failure modes of blade.The vibration fatigue life of blade is scattered caused by manufacture error,material property dispersion and external excitation randomness.A new vibration fatigue probabilistic life prediction model(VFPLPM)and a prediction method are proposed in this paper.Firstly,as one-dimensional volumetric method(ODVM)only considers the principle calculation direction,a three-dimensional space vector volumetric method(TSVVM)is proposed to improve fatigue life prediction accuracy for actual threedimensional engineering structure.Secondly,based on the two volumetric methods(ODVM and TSVVM),the material C-P-S-N fatigue curve model(CFCM)and the maximum entropy quantile function model(MEQFM),VFPLPM is established to predict the vibration fatigue probabilistic life of blade.The VFPLPM is combined with maximum stress method(MSM),ODVM and TSVVM to estimate vibration fatigue probabilistic life of blade simulator by finite element simulation,and is verified by vibration fatigue test.The results show that all of the three methods can predict the vibration fatigue probabilistic life of blade simulator well.VFPLPM &TSVVM method has the highest computational accuracy for considering stress gradient effect not only in the principle calculation direction but also in other space vector directions.
文摘Measurement of the volume of gas adsorbed per unit mass of coal with increasing pressure at a constant temperature produces an isotherm that describes the gas storage capacity of this type of coal. The accurate testing and interpretation of coal sorption isotherm plays an important role in the areas of coal mine methane drainage, coalbed methane （CBM） reservoir resource assessment, enhanced coalbed methane （ECBM） recovery, as well as the carbon dioxide （CO2） sequestration in deep coal seams or similar geological formations. Different coal sorption isotherm testing apparatus and associated calculation methods are critically reviewed and presented in this paper. These include both volumetric and gravimetric based methods, as well as experimental sorption tests with confining stress and direction sorption methods. The volumetric techniques utilise experimental apparatus with sample cell and injection pump and that with both sample cell and reference cell. Whilst the gravimetric approachesinclude methods with sample cell and suspension magnetic balance and that with both sample cell and reference cell. Different testing methods are compared and discussed in this study. A unique in-house-built coal sorption isotherm testing apparatus at the University of Wollongong was presented together with the calculation method, procedures and experimental results. The isotherm results can be calculated by both Soave-Redlich-Kwong （SRK） equation and calibration cure methods which can be used directly to convert the volume of adsorbed gas in different test conditions to standard condition （NTP）.
