With the development of industrial activities,global warming has accelerated due to excessive emission of CO_(2).Enhanced Geothermal System(EGS)utilizes deep geothermal heat for power generation.Although porous medium...With the development of industrial activities,global warming has accelerated due to excessive emission of CO_(2).Enhanced Geothermal System(EGS)utilizes deep geothermal heat for power generation.Although porous medium theory is commonly employed to model geothermal reservoirs in EGS,Hot Dry Rock(HDR)presents a challenge as it consists of impermeable granite with zero porosity,potentially distorting the physical interpretation.To address this,the Lattice Boltzmann Method(LBM)is employed to simulate CO_(2)flow within geothermal reservoirs and the Finite Volume Method(FVM)to solve the energy conservation equation for temperature distribution.This combined method of LBM and FVM is imple-mented using MATLAB.The results showed that the Reynolds numbers(Re)of 3,000 and 8,000 lead to higher heat extraction rates from geothermal reservoirs.However,higher Re values may accelerate thermal breakthrough,posing challenges to EGS operation.Meanwhile,non-equilibrium of density in fractures becomes more pronounced during the system's life cycle,with non-Darcy's law becoming significant at Re values of 3,000 and 8,000.Density stratification due to buoyancy effects significantly impacts temperature distribution within geothermal reservoirs,with buoyancy effects at Re=100 under gravitational influence being noteworthy.Larger Re values(3,000 and 8,000)induce stronger forced convection,leading to more uniform density distribution.The addition of proppant negatively affects heat transfer performance in geothermal reservoirs,especially in single fractures.Practical engineering considerations should determine the quantity of proppant through detailed numerical simulations.展开更多
In this paper,we consider the maximal positive definite solution of the nonlinear matrix equation.By using the idea of Algorithm 2.1 in ZHANG(2013),a new inversion-free method with a stepsize parameter is proposed to ...In this paper,we consider the maximal positive definite solution of the nonlinear matrix equation.By using the idea of Algorithm 2.1 in ZHANG(2013),a new inversion-free method with a stepsize parameter is proposed to obtain the maximal positive definite solution of nonlinear matrix equation X+A^(*)X|^(-α)A=Q with the case 0<α≤1.Based on this method,a new iterative algorithm is developed,and its convergence proof is given.Finally,two numerical examples are provided to show the effectiveness of the proposed method.展开更多
Marine thin plates are susceptible to welding deformation owing to their low structural stiffness.Therefore,the efficient and accurate prediction of welding deformation is essential for improving welding quality.The t...Marine thin plates are susceptible to welding deformation owing to their low structural stiffness.Therefore,the efficient and accurate prediction of welding deformation is essential for improving welding quality.The traditional thermal elastic-plastic finite element method(TEP-FEM)can accurately predict welding deformation.However,its efficiency is low because of the complex nonlinear transient computation,making it difficult to meet the needs of rapid engineering evaluation.To address this challenge,this study proposes an efficient prediction method for welding deformation in marine thin plate butt welds.This method is based on the coupled temperature gradient-thermal strain method(TG-TSM)that integrates inherent strain theory with a shell element finite element model.The proposed method first extracts the distribution pattern and characteristic value of welding-induced inherent strain through TEP-FEM analysis.This strain is then converted into the equivalent thermal load applied to the shell element model for rapid computation.The proposed method-particularly,the gradual temperature gradient-thermal strain method(GTG-TSM)-achieved improved computational efficiency and consistent precision.Furthermore,the proposed method required much less computation time than the traditional TEP-FEM.Thus,this study lays the foundation for future prediction of welding deformation in more complex marine thin plates.展开更多
Applying the standard Galerkin finite element method for solving flow problems in porous media encounters some difficulties such as numerical oscillation at the shock front and discontinuity of the velocity field on e...Applying the standard Galerkin finite element method for solving flow problems in porous media encounters some difficulties such as numerical oscillation at the shock front and discontinuity of the velocity field on element faces.Discontinuity of velocity field leads this method not to conserve mass locally.Moreover,the accuracy and stability of a solution is highly affected by a non-conservative method.In this paper,a three dimensional control volume finite element method is developed for twophase fluid flow simulation which overcomes the deficiency of the standard finite element method,and attains high-orders of accuracy at a reasonable computational cost.Moreover,this method is capable of handling heterogeneity in a very rational way.A fully implicit scheme is applied to temporal discretization of the governing equations to achieve an unconditionally stable solution.The accuracy and efficiency of the method are verified by simulating some waterflooding experiments.Some representative examples are presented to illustrate the capability of the method to simulate two-phase fluid flow in heterogeneous porous media.展开更多
A 3D compressible nonhydrostatic dynamic core based on a three-point multi-moment constrained finite-volume (MCV) method is developed by extending the previous 2D nonhydrostatic atmospheric dynamics to 3D on a terrain...A 3D compressible nonhydrostatic dynamic core based on a three-point multi-moment constrained finite-volume (MCV) method is developed by extending the previous 2D nonhydrostatic atmospheric dynamics to 3D on a terrainfollowing grid. The MCV algorithm defines two types of moments: the point-wise value (PV) and the volume-integrated average (VIA). The unknowns (PV values) are defined at the solution points within each cell and are updated through the time evolution formulations derived from the governing equations. Rigorous numerical conservation is ensured by a constraint on the VIA moment through the flux form formulation. The 3D atmospheric dynamic core reported in this paper is based on a three-point MCV method and has some advantages in comparison with other existing methods, such as uniform third-order accuracy, a compact stencil, and algorithmic simplicity. To check the performance of the 3D nonhydrostatic dynamic core, various benchmark test cases are performed. All the numerical results show that the present dynamic core is very competitive when compared to other existing advanced models, and thus lays the foundation for further developing global atmospheric models in the near future.展开更多
To evaluate the water storage and project the future evolution of glaciers, the ice-thickness of glaciers is an essential input. However, direct measurements of ice thickness are labo- rious, not feasible everywhere, ...To evaluate the water storage and project the future evolution of glaciers, the ice-thickness of glaciers is an essential input. However, direct measurements of ice thickness are labo- rious, not feasible everywhere, and necessarily restricted to a small number of glaciers. In this article, we develop a simple method to estimate the ice-thickness along flow-line of mountain glaciers. Different from the traditional method based on shallow ice approximation （SIA）, which gives a relationship be- tween ice thickness, surface slope, and yield stress of glaciers, the improved method considers and pre- sents a simple way to calibrate the influence of valley wall on ice discharge. The required inputs are the glacier surface topography and outlines. This shows the potential of the method for estimating the ice-thickness distribution and volume of glaciers without using of direct thickness measurements.展开更多
In order to obtain the accurate sonic conductance of solenoid valves, the polytropic exponent is used in the data processing of the discharge method for measuring the flow rate characteristics of pneumatic components....In order to obtain the accurate sonic conductance of solenoid valves, the polytropic exponent is used in the data processing of the discharge method for measuring the flow rate characteristics of pneumatic components. Three data processing principles are first introduced, and then the discharge pressure data obtained from the measured solenoid valve are processed to obtain the sonic conductance with three methods： the partial polytropic exponent method, the complete polytropic exponent method and the adiabatic method. By comparison of the obtained results it indicates that the complete polytropic exponent method is the most accurate. However, the partial polytropic exponent method is of a high applicable value, because it is easy and simple to measure and the obtained results are relatively accurate.展开更多
The finite volume method has been successfully applied in several engineering fields and has shown outstanding performance in fluid dynamics simulation. In this paper, the general framework for the simulation ofnear-w...The finite volume method has been successfully applied in several engineering fields and has shown outstanding performance in fluid dynamics simulation. In this paper, the general framework for the simulation ofnear-wellbore systems using the finite volume method is described. The mathematical model and the numerical model developed by the authors are presented and discussed. A radial geometry in the vertical plane was implemented so as to thoroughly describe near-wellbore phenomena. The model was then used to simulate injection tests in an oil reservoir through a horizontal well and proved very powerful to correctly reproduce the transient pressure behavior. The reason for this is the robustness of the method, which is independent of the gridding options because the discretization is performed in the physical space. The model is able to describe the phenomena taking place in the reservoir even in complex situations, i.e. in the presence of heterogeneities and permeability barriers, demonstrating the flexibility of the finite volume method when simulating non-conventional tests. The results are presented in comparison with those obtained with the finite difference numerical approach and with analytical methods, if possible.展开更多
The finite volume method （FVM） has many advantages in 2-D shallow water numerical simulation. In this study, the finite volume method is used with unstructured triangular grids to simulate the tidal currents. The Ro...The finite volume method （FVM） has many advantages in 2-D shallow water numerical simulation. In this study, the finite volume method is used with unstructured triangular grids to simulate the tidal currents. The Roe scheme is applied in the calculation of the intercell numerical flux, and the MUSCL method is introduced to improve its accuracy. The time integral is a two-step scheme of forecast and revision. For the verification of the present method, the Stoker＇s problem is calculated and the result is compared with the mathematically analytic solutions. The comparison indicates that the method is feasible. A sea area of a port is used as an example to test the method established here. The result shows that the present computational method is satisfactory, and it could be applied to the engineering fields.展开更多
A method to simulate processes of forging and subsequent heat treatment of an axial symmetric rod is formulated in eulerian description and the feasibility is investigated. This method uses finite volume mushes for t...A method to simulate processes of forging and subsequent heat treatment of an axial symmetric rod is formulated in eulerian description and the feasibility is investigated. This method uses finite volume mushes for troching material deformation and an automatically refined facet surface to accurately trace the free surface of the deforming material.In the method,the deforming work piece flows through fixed finite volume meshes using eulerian formulation to describe the conservation laws,Fixed finite volume meshing is particularly suitable for large three-dimensional deformation such as forging because remeshing techniques are not required, which are commonly considered to be the main bottelencek in the ssimulations of large defromation by using the finite element method,By means of this finite volume method, an approach has been developed in the framework of 'metallo-thermo-mechanics' to simulate metallic structure, temperature and stress/strain coupled in the heat treatment process.In a first step of simulation, the heat treatment solver is limited in small deformation hypothesis,and un- coupled with forging. The material is considered as elastic-plastic and takes into account of strain, strain rate and temperature effects on the yield stress.Heat generation due to deformation,heat con- duction and thermal stress are considered.Temperature - dependent phase transformation,stress-in- duced phase transformation,latent heat,transformation stress and strain are included.These ap- proaches are implemented into the commerical commercial computer program MSC/SuperForge and a verification example with experimental date is given as comparison.展开更多
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).展开更多
This paper presents alternating direction finite volume element methods for three-dimensional parabolic partial differential equations and gives four computational schemes, one is analogous to Douglas finite differenc...This paper presents alternating direction finite volume element methods for three-dimensional parabolic partial differential equations and gives four computational schemes, one is analogous to Douglas finite difference scheme with second-order splitting error, the other two schemes have third-order splitting error, and the last one is an extended LOD scheme. The L2 norm and H1 semi-norm error estimates are obtained for the first scheme and second one, respectively. Finally, two numerical examples are provided to illustrate the efficiency and accuracy of the methods.展开更多
The mixed covolume method for the regularized long wave equation is devel- oped and studied. By introducing a transfer operator γh, which maps the trial function space into the test function space, and combining the ...The mixed covolume method for the regularized long wave equation is devel- oped and studied. By introducing a transfer operator γh, which maps the trial function space into the test function space, and combining the mixed finite element with the finite volume method, the nonlinear and linear Euler fully discrete mixed covolume schemes are constructed, and the existence and uniqueness of the solutions are proved. The optimal error estimates for these schemes are obtained. Finally, a numerical example is provided to examine the efficiency of the proposed schemes.展开更多
Fractional-in-space Allen-Cahn equation containing a very strong nonlinear source term and small perturbation shows metastability and a quartic double well potential.Using a finite volume unstructured triangular mesh ...Fractional-in-space Allen-Cahn equation containing a very strong nonlinear source term and small perturbation shows metastability and a quartic double well potential.Using a finite volume unstructured triangular mesh method, the present paper solves the twodimensional fractional-in-space Allen-Cahn equation with homogeneous Neumann boundary condition on different irregular domains. The efficiency of the method is presented through numerical computation of the two-dimensional fractional-in-space Allen-Cahn equation on different domains.展开更多
This paper proposes a hybrid vertex-centered fi- nite volume/finite element method for solution of the two di- mensional （2D） incompressible Navier-Stokes equations on unstructured grids. An incremental pressure fra...This paper proposes a hybrid vertex-centered fi- nite volume/finite element method for solution of the two di- mensional （2D） incompressible Navier-Stokes equations on unstructured grids. An incremental pressure fractional step method is adopted to handle the velocity-pressure coupling. The velocity and the pressure are collocated at the node of the vertex-centered control volume which is formed by join- ing the centroid of cells sharing the common vertex. For the temporal integration of the momentum equations, an im- plicit second-order scheme is utilized to enhance the com- putational stability and eliminate the time step limit due to the diffusion term. The momentum equations are discretized by the vertex-centered finite volume method （FVM） and the pressure Poisson equation is solved by the Galerkin finite el- ement method （FEM）. The momentum interpolation is used to damp out the spurious pressure wiggles. The test case with analytical solutions demonstrates second-order accuracy of the current hybrid scheme in time and space for both veloc- ity and pressure. The classic test cases, the lid-driven cavity flow, the skew cavity flow and the backward-facing step flow, show that numerical results are in good agreement with the published benchmark solutions.展开更多
Level set methods are widely used for predicting evolutions of complex free surface topologies,such as the crystal and crack growth,bubbles and droplets deformation,spilling and breaking waves,and two-phase flow pheno...Level set methods are widely used for predicting evolutions of complex free surface topologies,such as the crystal and crack growth,bubbles and droplets deformation,spilling and breaking waves,and two-phase flow phenomena.This paper presents a characteristic level set equation which is derived from the two-dimensional level set equation by using the characteristic-based scheme.An explicit finite volume element method is developed to discretize the equation on triangular grids.Several examples are presented to demonstrate the performance of the proposed method for calculating interface evolutions in time.The proposed level set method is also coupled with the Navier-Stokes equations for two-phase immiscible incompressible flow analysis with surface tension.The Rayleigh-Taylor instability problem is used to test and evaluate the effectiveness of the proposed scheme.展开更多
In this paper, a high accuracy finite volume element method is presented for two-point boundary value problem of second order ordinary differential equation, which differs from the high order generalized difference me...In this paper, a high accuracy finite volume element method is presented for two-point boundary value problem of second order ordinary differential equation, which differs from the high order generalized difference methods. It is proved that the method has optimal order error estimate O(h3) in H1 norm. Finally, two examples show that the method is effective.展开更多
The fluid-coupled discrete element method(DEM)and the constant volume method as two types of discrete modeling methods for fundamental study of undrained responses of granular materials,have been discussed by many res...The fluid-coupled discrete element method(DEM)and the constant volume method as two types of discrete modeling methods for fundamental study of undrained responses of granular materials,have been discussed by many researchers.The fluid-coupled DEM,which couples the motions of discrete particles with pore fluid movements,is theoretically robust although it requires a large amount of computation time.As a substitution for the complex fluid-coupled DEM,the constant volume method simulates an undrained condition for a saturated granular material by simply preserving the total volume of a granular assembly without considering interactions between fluids and particles;hence,the validity of its results is questionable.In this paper,the undrained behaviors of granular assemblies simulated using the aforementioned two methods are compared.Based on a comparison of both macroscopic and microscopic responses given by the two methods,it is demonstrated that the constant volume method may reasonably simulate the responses of a loose saturated granular material with very coarse grains,which has a high permeability,and thus a rapid pore pressure equalization.However,it is ineffective in simulating the responses of a loose material with fine components due to its failure to capture the process of a slow dissipation of the excess pore pressure among the individual pores.With regard to the dense material adopted,similar behaviors at the early and intermediate shearing stages given by the two methods are displayed.展开更多
1 Introduction Hydraulic fracturing is a technique for increasing permeability in oil and gas resource development,grouting reinforcement in mine management,and geostress measurement.For the purpose of enhancing hydra...1 Introduction Hydraulic fracturing is a technique for increasing permeability in oil and gas resource development,grouting reinforcement in mine management,and geostress measurement.For the purpose of enhancing hydraulic fracturing in horizontal wells,oriented perforating methods have been developed(Kurdi,2018;Michael and Gupta,2020a;Yan et al.,2020).Fluid is injected into the rock through perforations,which increases fluid pressure within rock and decreases rock temperature.Then,the rock around the perforation is fractured.Therefore,fracture initiation pressure is intimately connected to the reservoir’s physical and mechanical properties,geo-stress,and temperature(Morgan and Aral,2015).展开更多
Numerical simulations are performed to examine the packing behavior of human red blood cells(RBCs). A combined ?nite-discrete element method(FDEM) is utilized, in which the RBCs are modeled as no-friction and no-adhes...Numerical simulations are performed to examine the packing behavior of human red blood cells(RBCs). A combined ?nite-discrete element method(FDEM) is utilized, in which the RBCs are modeled as no-friction and no-adhesion solid bodies. The packed volume and the void ratio of a large number of randomly packed RBCs are clari?ed,and the effects of the RBC shape, the mesh size, the cell number, and the container size are investigated. The results show that the packed human RBCs with normal shape have a void ratio of 28.45%, which is slightly higher than that of the ?at or thick cells used in this study. Such information is bene?cial to the further understanding on the geometric features of human RBCs and the research on RBC simulations.展开更多
基金supported by the Hebei Province Graduate Innovation Funding Project(CXZZBS2022029).
文摘With the development of industrial activities,global warming has accelerated due to excessive emission of CO_(2).Enhanced Geothermal System(EGS)utilizes deep geothermal heat for power generation.Although porous medium theory is commonly employed to model geothermal reservoirs in EGS,Hot Dry Rock(HDR)presents a challenge as it consists of impermeable granite with zero porosity,potentially distorting the physical interpretation.To address this,the Lattice Boltzmann Method(LBM)is employed to simulate CO_(2)flow within geothermal reservoirs and the Finite Volume Method(FVM)to solve the energy conservation equation for temperature distribution.This combined method of LBM and FVM is imple-mented using MATLAB.The results showed that the Reynolds numbers(Re)of 3,000 and 8,000 lead to higher heat extraction rates from geothermal reservoirs.However,higher Re values may accelerate thermal breakthrough,posing challenges to EGS operation.Meanwhile,non-equilibrium of density in fractures becomes more pronounced during the system's life cycle,with non-Darcy's law becoming significant at Re values of 3,000 and 8,000.Density stratification due to buoyancy effects significantly impacts temperature distribution within geothermal reservoirs,with buoyancy effects at Re=100 under gravitational influence being noteworthy.Larger Re values(3,000 and 8,000)induce stronger forced convection,leading to more uniform density distribution.The addition of proppant negatively affects heat transfer performance in geothermal reservoirs,especially in single fractures.Practical engineering considerations should determine the quantity of proppant through detailed numerical simulations.
基金Supported in part by Natural Science Foundation of Guangxi(2023GXNSFAA026246)in part by the Central Government's Guide to Local Science and Technology Development Fund(GuikeZY23055044)in part by the National Natural Science Foundation of China(62363003)。
文摘In this paper,we consider the maximal positive definite solution of the nonlinear matrix equation.By using the idea of Algorithm 2.1 in ZHANG(2013),a new inversion-free method with a stepsize parameter is proposed to obtain the maximal positive definite solution of nonlinear matrix equation X+A^(*)X|^(-α)A=Q with the case 0<α≤1.Based on this method,a new iterative algorithm is developed,and its convergence proof is given.Finally,two numerical examples are provided to show the effectiveness of the proposed method.
基金Supported by the National Natural Science Foundation of China under Grant No.51975138the High-Tech Ship Scientific Research Project from the Ministry of Industry and Information Technology under Grant No.CJ05N20the National Defense Basic Research Project under Grant No.JCKY2023604C006.
文摘Marine thin plates are susceptible to welding deformation owing to their low structural stiffness.Therefore,the efficient and accurate prediction of welding deformation is essential for improving welding quality.The traditional thermal elastic-plastic finite element method(TEP-FEM)can accurately predict welding deformation.However,its efficiency is low because of the complex nonlinear transient computation,making it difficult to meet the needs of rapid engineering evaluation.To address this challenge,this study proposes an efficient prediction method for welding deformation in marine thin plate butt welds.This method is based on the coupled temperature gradient-thermal strain method(TG-TSM)that integrates inherent strain theory with a shell element finite element model.The proposed method first extracts the distribution pattern and characteristic value of welding-induced inherent strain through TEP-FEM analysis.This strain is then converted into the equivalent thermal load applied to the shell element model for rapid computation.The proposed method-particularly,the gradual temperature gradient-thermal strain method(GTG-TSM)-achieved improved computational efficiency and consistent precision.Furthermore,the proposed method required much less computation time than the traditional TEP-FEM.Thus,this study lays the foundation for future prediction of welding deformation in more complex marine thin plates.
基金Iranian Offshore Oil Company (IOOC) for financial support of this work
文摘Applying the standard Galerkin finite element method for solving flow problems in porous media encounters some difficulties such as numerical oscillation at the shock front and discontinuity of the velocity field on element faces.Discontinuity of velocity field leads this method not to conserve mass locally.Moreover,the accuracy and stability of a solution is highly affected by a non-conservative method.In this paper,a three dimensional control volume finite element method is developed for twophase fluid flow simulation which overcomes the deficiency of the standard finite element method,and attains high-orders of accuracy at a reasonable computational cost.Moreover,this method is capable of handling heterogeneity in a very rational way.A fully implicit scheme is applied to temporal discretization of the governing equations to achieve an unconditionally stable solution.The accuracy and efficiency of the method are verified by simulating some waterflooding experiments.Some representative examples are presented to illustrate the capability of the method to simulate two-phase fluid flow in heterogeneous porous media.
基金supported by the National Key Research and Development Program of China (Grant Nos. 2017YFC1501901 and 2017YFA0603901)the Beijing Natural Science Foundation (Grant No. JQ18001)
文摘A 3D compressible nonhydrostatic dynamic core based on a three-point multi-moment constrained finite-volume (MCV) method is developed by extending the previous 2D nonhydrostatic atmospheric dynamics to 3D on a terrainfollowing grid. The MCV algorithm defines two types of moments: the point-wise value (PV) and the volume-integrated average (VIA). The unknowns (PV values) are defined at the solution points within each cell and are updated through the time evolution formulations derived from the governing equations. Rigorous numerical conservation is ensured by a constraint on the VIA moment through the flux form formulation. The 3D atmospheric dynamic core reported in this paper is based on a three-point MCV method and has some advantages in comparison with other existing methods, such as uniform third-order accuracy, a compact stencil, and algorithmic simplicity. To check the performance of the 3D nonhydrostatic dynamic core, various benchmark test cases are performed. All the numerical results show that the present dynamic core is very competitive when compared to other existing advanced models, and thus lays the foundation for further developing global atmospheric models in the near future.
基金supported by the National Basic Research Program of China (No. 2007CB411501)the Knowledge Innovation Project of the Chinese Academy of Sciences (No. KZCX2-EW-311)+1 种基金the National Natural Science Foundation of China (Nos. 91025012, J0930003/J0109)the Project for Outstanding Young Scientists of the National Natural Science Foundation of China (No. 40121101)
文摘To evaluate the water storage and project the future evolution of glaciers, the ice-thickness of glaciers is an essential input. However, direct measurements of ice thickness are labo- rious, not feasible everywhere, and necessarily restricted to a small number of glaciers. In this article, we develop a simple method to estimate the ice-thickness along flow-line of mountain glaciers. Different from the traditional method based on shallow ice approximation （SIA）, which gives a relationship be- tween ice thickness, surface slope, and yield stress of glaciers, the improved method considers and pre- sents a simple way to calibrate the influence of valley wall on ice discharge. The required inputs are the glacier surface topography and outlines. This shows the potential of the method for estimating the ice-thickness distribution and volume of glaciers without using of direct thickness measurements.
文摘In order to obtain the accurate sonic conductance of solenoid valves, the polytropic exponent is used in the data processing of the discharge method for measuring the flow rate characteristics of pneumatic components. Three data processing principles are first introduced, and then the discharge pressure data obtained from the measured solenoid valve are processed to obtain the sonic conductance with three methods： the partial polytropic exponent method, the complete polytropic exponent method and the adiabatic method. By comparison of the obtained results it indicates that the complete polytropic exponent method is the most accurate. However, the partial polytropic exponent method is of a high applicable value, because it is easy and simple to measure and the obtained results are relatively accurate.
文摘The finite volume method has been successfully applied in several engineering fields and has shown outstanding performance in fluid dynamics simulation. In this paper, the general framework for the simulation ofnear-wellbore systems using the finite volume method is described. The mathematical model and the numerical model developed by the authors are presented and discussed. A radial geometry in the vertical plane was implemented so as to thoroughly describe near-wellbore phenomena. The model was then used to simulate injection tests in an oil reservoir through a horizontal well and proved very powerful to correctly reproduce the transient pressure behavior. The reason for this is the robustness of the method, which is independent of the gridding options because the discretization is performed in the physical space. The model is able to describe the phenomena taking place in the reservoir even in complex situations, i.e. in the presence of heterogeneities and permeability barriers, demonstrating the flexibility of the finite volume method when simulating non-conventional tests. The results are presented in comparison with those obtained with the finite difference numerical approach and with analytical methods, if possible.
基金This paper was supported bythe Natural Science Foundation of Shandong Province (Grant No.y2004f13)
文摘The finite volume method （FVM） has many advantages in 2-D shallow water numerical simulation. In this study, the finite volume method is used with unstructured triangular grids to simulate the tidal currents. The Roe scheme is applied in the calculation of the intercell numerical flux, and the MUSCL method is introduced to improve its accuracy. The time integral is a two-step scheme of forecast and revision. For the verification of the present method, the Stoker＇s problem is calculated and the result is compared with the mathematically analytic solutions. The comparison indicates that the method is feasible. A sea area of a port is used as an example to test the method established here. The result shows that the present computational method is satisfactory, and it could be applied to the engineering fields.
文摘A method to simulate processes of forging and subsequent heat treatment of an axial symmetric rod is formulated in eulerian description and the feasibility is investigated. This method uses finite volume mushes for troching material deformation and an automatically refined facet surface to accurately trace the free surface of the deforming material.In the method,the deforming work piece flows through fixed finite volume meshes using eulerian formulation to describe the conservation laws,Fixed finite volume meshing is particularly suitable for large three-dimensional deformation such as forging because remeshing techniques are not required, which are commonly considered to be the main bottelencek in the ssimulations of large defromation by using the finite element method,By means of this finite volume method, an approach has been developed in the framework of 'metallo-thermo-mechanics' to simulate metallic structure, temperature and stress/strain coupled in the heat treatment process.In a first step of simulation, the heat treatment solver is limited in small deformation hypothesis,and un- coupled with forging. The material is considered as elastic-plastic and takes into account of strain, strain rate and temperature effects on the yield stress.Heat generation due to deformation,heat con- duction and thermal stress are considered.Temperature - dependent phase transformation,stress-in- duced phase transformation,latent heat,transformation stress and strain are included.These ap- proaches are implemented into the commerical commercial computer program MSC/SuperForge and a verification example with experimental date is given as comparison.
基金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).
文摘This paper presents alternating direction finite volume element methods for three-dimensional parabolic partial differential equations and gives four computational schemes, one is analogous to Douglas finite difference scheme with second-order splitting error, the other two schemes have third-order splitting error, and the last one is an extended LOD scheme. The L2 norm and H1 semi-norm error estimates are obtained for the first scheme and second one, respectively. Finally, two numerical examples are provided to illustrate the efficiency and accuracy of the methods.
基金supported by the National Natural Science Fundation of China (No. 11061021)the Science Research of Inner Mongolia Advanced Education (Nos. NJ10006, NJ10016, and NJZZ12011)the National Science Foundation of Inner Mongolia (Nos. 2011BS0102 and 2012MS0106)
文摘The mixed covolume method for the regularized long wave equation is devel- oped and studied. By introducing a transfer operator γh, which maps the trial function space into the test function space, and combining the mixed finite element with the finite volume method, the nonlinear and linear Euler fully discrete mixed covolume schemes are constructed, and the existence and uniqueness of the solutions are proved. The optimal error estimates for these schemes are obtained. Finally, a numerical example is provided to examine the efficiency of the proposed schemes.
基金Supported by the National Natural Science Foundation of China(11105040,61773153)Supported by the Foundation of Henan Educational Committee(18B110003,15A110015)+1 种基金Supported by the Excellent Young Scientific Talents Cultivation Foundation of Henan University(yqpy20140037)Supported by the Science and Technology Program of Henan Province(162300410061)
文摘Fractional-in-space Allen-Cahn equation containing a very strong nonlinear source term and small perturbation shows metastability and a quartic double well potential.Using a finite volume unstructured triangular mesh method, the present paper solves the twodimensional fractional-in-space Allen-Cahn equation with homogeneous Neumann boundary condition on different irregular domains. The efficiency of the method is presented through numerical computation of the two-dimensional fractional-in-space Allen-Cahn equation on different domains.
基金supported by the Natural Science Foundation of China (11061021)the Program of Higher-level talents of Inner Mongolia University (SPH-IMU,Z200901004)the Scientific Research Projection of Higher Schools of Inner Mongolia(NJ10016,NJ10006)
文摘This paper proposes a hybrid vertex-centered fi- nite volume/finite element method for solution of the two di- mensional （2D） incompressible Navier-Stokes equations on unstructured grids. An incremental pressure fractional step method is adopted to handle the velocity-pressure coupling. The velocity and the pressure are collocated at the node of the vertex-centered control volume which is formed by join- ing the centroid of cells sharing the common vertex. For the temporal integration of the momentum equations, an im- plicit second-order scheme is utilized to enhance the com- putational stability and eliminate the time step limit due to the diffusion term. The momentum equations are discretized by the vertex-centered finite volume method （FVM） and the pressure Poisson equation is solved by the Galerkin finite el- ement method （FEM）. The momentum interpolation is used to damp out the spurious pressure wiggles. The test case with analytical solutions demonstrates second-order accuracy of the current hybrid scheme in time and space for both veloc- ity and pressure. The classic test cases, the lid-driven cavity flow, the skew cavity flow and the backward-facing step flow, show that numerical results are in good agreement with the published benchmark solutions.
基金King Mongkut’s University of Technology North Bangkok (KMUTNB)the Office of the Higher Education Commission (OHEC)the National Metal and Materials Technology Center (MTEC) for supporting this research work
文摘Level set methods are widely used for predicting evolutions of complex free surface topologies,such as the crystal and crack growth,bubbles and droplets deformation,spilling and breaking waves,and two-phase flow phenomena.This paper presents a characteristic level set equation which is derived from the two-dimensional level set equation by using the characteristic-based scheme.An explicit finite volume element method is developed to discretize the equation on triangular grids.Several examples are presented to demonstrate the performance of the proposed method for calculating interface evolutions in time.The proposed level set method is also coupled with the Navier-Stokes equations for two-phase immiscible incompressible flow analysis with surface tension.The Rayleigh-Taylor instability problem is used to test and evaluate the effectiveness of the proposed scheme.
基金heprojectissupportedbyNNSFofChina (No .1 9972 0 39) .
文摘In this paper, a high accuracy finite volume element method is presented for two-point boundary value problem of second order ordinary differential equation, which differs from the high order generalized difference methods. It is proved that the method has optimal order error estimate O(h3) in H1 norm. Finally, two examples show that the method is effective.
基金The authors would like to express their acknowledgments to Dr.Timothy Topper for his insightful com。
文摘The fluid-coupled discrete element method(DEM)and the constant volume method as two types of discrete modeling methods for fundamental study of undrained responses of granular materials,have been discussed by many researchers.The fluid-coupled DEM,which couples the motions of discrete particles with pore fluid movements,is theoretically robust although it requires a large amount of computation time.As a substitution for the complex fluid-coupled DEM,the constant volume method simulates an undrained condition for a saturated granular material by simply preserving the total volume of a granular assembly without considering interactions between fluids and particles;hence,the validity of its results is questionable.In this paper,the undrained behaviors of granular assemblies simulated using the aforementioned two methods are compared.Based on a comparison of both macroscopic and microscopic responses given by the two methods,it is demonstrated that the constant volume method may reasonably simulate the responses of a loose saturated granular material with very coarse grains,which has a high permeability,and thus a rapid pore pressure equalization.However,it is ineffective in simulating the responses of a loose material with fine components due to its failure to capture the process of a slow dissipation of the excess pore pressure among the individual pores.With regard to the dense material adopted,similar behaviors at the early and intermediate shearing stages given by the two methods are displayed.
基金This work is supported by the National Natural Science Foundation of China(Nos.51890914 and 52179119)the Natural Science Foundation of Shandong Province(No.ZR2019MEE001)the Open Research Fund of Hunan Provincial Key Laboratory of Hydropower Development Key Technology(No.PKLHD202001),China。
文摘1 Introduction Hydraulic fracturing is a technique for increasing permeability in oil and gas resource development,grouting reinforcement in mine management,and geostress measurement.For the purpose of enhancing hydraulic fracturing in horizontal wells,oriented perforating methods have been developed(Kurdi,2018;Michael and Gupta,2020a;Yan et al.,2020).Fluid is injected into the rock through perforations,which increases fluid pressure within rock and decreases rock temperature.Then,the rock around the perforation is fractured.Therefore,fracture initiation pressure is intimately connected to the reservoir’s physical and mechanical properties,geo-stress,and temperature(Morgan and Aral,2015).
基金Project supported by the Engineering and Physical Sciences Research Council(EPSRC)Turbulence Consortium Grant(No.EP/G069581/1)the Marie Curie International Incoming Fellowship(No.PIIF-GA-253453)
文摘Numerical simulations are performed to examine the packing behavior of human red blood cells(RBCs). A combined ?nite-discrete element method(FDEM) is utilized, in which the RBCs are modeled as no-friction and no-adhesion solid bodies. The packed volume and the void ratio of a large number of randomly packed RBCs are clari?ed,and the effects of the RBC shape, the mesh size, the cell number, and the container size are investigated. The results show that the packed human RBCs with normal shape have a void ratio of 28.45%, which is slightly higher than that of the ?at or thick cells used in this study. Such information is bene?cial to the further understanding on the geometric features of human RBCs and the research on RBC simulations.
