A new quadrilateral edge element method is proposed and analyzed for Maxwell equations.This proposed method is based on Duan-Liang quadrilateral element(Math.Comp.73(2004),pp.1–18).When applied to the eigenvalue prob...A new quadrilateral edge element method is proposed and analyzed for Maxwell equations.This proposed method is based on Duan-Liang quadrilateral element(Math.Comp.73(2004),pp.1–18).When applied to the eigenvalue problem,the method is spectral-correct and spurious-free.Stability and error estimates are obtained,including the interpolation error estimates and the error estimates between the finite element solution and the exact solution.The method is suitable for singular solution as well as smooth solution,and consequently,the method is valid for nonconvex domains which may have a number of reentrant corners.Of course,the method is suitable for arbitrary quadrilaterals(under the usual shape-regular condition).展开更多
This study explores a sensitivity analysis method based on the boundary element method(BEM)to address the computational complexity in acoustic analysis with ground reflection problems.The advantages of BEM in acoustic...This study explores a sensitivity analysis method based on the boundary element method(BEM)to address the computational complexity in acoustic analysis with ground reflection problems.The advantages of BEM in acoustic simulations and its high computational cost in broadband problems are examined.To improve efficiency,a Taylor series expansion is applied to decouple frequency-dependent terms in BEM.Additionally,the SecondOrder Arnoldi(SOAR)model order reduction method is integrated to reduce computational costs and enhance numerical stability.Furthermore,an isogeometric sensitivity boundary integral equation is formulated using the direct differentiation method,incorporating Cauchy principal value integrals and Hadamard finite part integrals to handle singularities.The proposed method improves the computational efficiency,and the acoustic sensitivity analysis provides theoretical support for further acoustic structure optimization.展开更多
Sudden and unforeseen seismic failures of coal mine overburden(OB)dump slopes interrupt mining operations,cause loss of lives and delay the production of coal.Consideration of the spatial heterogeneity of OB dump mate...Sudden and unforeseen seismic failures of coal mine overburden(OB)dump slopes interrupt mining operations,cause loss of lives and delay the production of coal.Consideration of the spatial heterogeneity of OB dump materials is imperative for an adequate evaluation of the seismic stability of OB dump slopes.In this study,pseudo-static seismic stability analyses are carried out for an OB dump slope by considering the material parameters obtained from an insitu field investigation.Spatial heterogeneity is simulated through use of the random finite element method(RFEM)and the random limit equilibrium method(RLEM)and a comparative study is presented.Combinations of horizontal and vertical spatial correlation lengths were considered for simulating isotropic and anisotropic random fields within the OB dump slope.Seismic performances of the slope have been reported through the probability of failure and reliability index.It was observed that the RLEM approach overestimates failure probability(P_(f))by considering seismic stability with spatial heterogeneity.The P_(f)was observed to increase with an increase in the coefficient of variation of friction angle of the dump materials.Further,it was inferred that the RLEM approach may not be adequately applicable for assessing the seismic stability of an OB dump slope for a horizontal seismic coefficient that is more than or equal to 0.1.展开更多
The stochastic boundary element method(SBEM)is developed in this paper for 3D problems with body forces and reliability analysis of engineering structures.The integral equations of SBEM are established by the approach...The stochastic boundary element method(SBEM)is developed in this paper for 3D problems with body forces and reliability analysis of engineering structures.The integral equations of SBEM are established by the approach of partial derivation with respect to stochastic variables,considering the yield limit,rotation speeds and material density to be the fundamental stochastic variables.Through analyzing a numerical example and a turbo-disk of an aeroengine,the results show that the method developed is successful.展开更多
To analyze the differences in the transport and distribution of different types of proppants and to address issues such as the short effective support of proppant and poor placement in hydraulically intersecting fract...To analyze the differences in the transport and distribution of different types of proppants and to address issues such as the short effective support of proppant and poor placement in hydraulically intersecting fractures,this study considered the combined impact of geological-engineering factors on conductivity.Using reservoir production parameters and the discrete elementmethod,multispherical proppants were constructed.Additionally,a 3D fracture model,based on the specified conditions of the L block,employed coupled(Computational Fluid Dynamics)CFD-DEM(Discrete ElementMethod)for joint simulations to quantitatively analyze the transport and placement patterns of multispherical proppants in intersecting fractures.Results indicate that turbulent kinetic energy is an intrinsic factor affecting proppant transport.Moreover,the efficiency of placement and migration distance of low-sphericity quartz sand constructed by the DEM in the main fracture are significantly reduced compared to spherical ceramic proppants,with a 27.7%decrease in the volume fraction of the fracture surface,subsequently affecting the placement concentration and damaging fracture conductivity.Compared to small-angle fractures,controlling artificial and natural fractures to expand at angles of 45°to 60°increases the effective support length by approximately 20.6%.During hydraulic fracturing of gas wells,ensuring the fracture support area and post-closure conductivity can be achieved by controlling the sphericity of proppants and adjusting the perforation direction to control the direction of artificial fractures.展开更多
Magnetization configurations were calculated under various magnetic fields for nanocrystalline Pr-Fe-B permanent magnets by micromagnetic finite element method.According to the configurations during demagnetization pr...Magnetization configurations were calculated under various magnetic fields for nanocrystalline Pr-Fe-B permanent magnets by micromagnetic finite element method.According to the configurations during demagnetization process, the mechanism of magnetization reversal was analyzed.For the Pr2Fe14B with 10 nm grains or its composite with 10vol.% α-Fe, the coercivity was determined by nucleation of reversed domain that took place at grain boundaries.However, for Pr2Fe14B with 30 nm grains, coercivity was controlled by pinning of the nucle-ated domain.For Pr2Fe14B/α-Fe with 30vol.% α-Fe, the demagnetization behavior was characterized by continuous reversal of α-Fe moment.展开更多
In this paper,two-grid immersed finite element (IFE) algorithms are proposed and analyzed for semi-linear interface problems with discontinuous diffusion coefficients in two dimension.Because of the advantages of fini...In this paper,two-grid immersed finite element (IFE) algorithms are proposed and analyzed for semi-linear interface problems with discontinuous diffusion coefficients in two dimension.Because of the advantages of finite element (FE) formulation and the simple structure of Cartesian grids,the IFE discretization is used in this paper.Two-grid schemes are formulated to linearize the FE equations.It is theoretically and numerically illustrated that the coarse space can be selected as coarse as H =O(h^1/4)(or H =O(h^1/8)),and the asymptotically optimal approximation can be achieved as the nonlinear schemes.As a result,we can settle a great majority of nonlinear equations as easy as linearized problems.In order to estimate the present two-grid algorithms,we derive the optimal error estimates of the IFE solution in the L^p norm.Numerical experiments are given to verify the theorems and indicate that the present two-grid algorithms can greatly improve the computing efficiency.展开更多
The past decade has witnessed the substantial growth in research interests and progress on the subject of coupled hydro-mechanical processes in rocks and soils,driven mainly by the surge of research in unconventional ...The past decade has witnessed the substantial growth in research interests and progress on the subject of coupled hydro-mechanical processes in rocks and soils,driven mainly by the surge of research in unconventional hydrocarbon reservoirs and associated hazards.Many coupling techniques have been developed to include the effects of fluid flow in the discrete element method(DEM),and the techniques have been applied to a variety of geomechanical problems.Although these coupling methods have been successfully applied in various engineering fields,no single fluid/DEM coupling method is universal due to the complexity of engineering problems and the limitations of the numerical methods.For researchers and engineers,the key to solve a specific problem is to select the most appropriate fluid/DEM coupling method among these modeling technologies.The purpose of this paper is to give a comprehensive review of fluid flow/DEM coupling methods and relevant research.Given their importance,the availability or unavailability of best practice guidelines is outlined.The theoretical background and current status of DEM are introduced first,and the principles,applications,and advantages and disadvantages of different fluid flow/DEM coupling methods are discussed.Finally,a summary with speculation on future development trends is given.展开更多
This paper discusses the forward and inverse problem for cardiac magnetic fields and electric potentials. A torso-heart model established by boundary element method （BEM） is used for studying the distributions of ca...This paper discusses the forward and inverse problem for cardiac magnetic fields and electric potentials. A torso-heart model established by boundary element method （BEM） is used for studying the distributions of cardiac magnetic fields and electric potentials. Because node-to-node and triangle-to-triangle BEM can lead to discrepant field distributions, their properties and influences are compared. Then based on constructed torso-heart model and supposed current source functional model-current dipole array, the magnetic and electric imaging by optimal constrained linear inverse method are applied at the same time. Through figure and reconstructing parameter comparison, though the magnetic current dipole array imaging possesses better reconstructing effect, however node-to-node BEM and triangleto-triangle BEM make little difference to magnetic and electric imaging.展开更多
In this article, on the basis of two-level discretizations and multiscale finite element method, two kinds of finite element algorithms for steady Navier-Stokes problem are presented and discussed. The main technique ...In this article, on the basis of two-level discretizations and multiscale finite element method, two kinds of finite element algorithms for steady Navier-Stokes problem are presented and discussed. The main technique is first to use a standard finite element discretization on a coarse mesh to approximate low frequencies, then to apply the simple and Newton scheme to linearize discretizations on a fine grid. At this process, multiscale finite element method as a stabilized method deals with the lowest equal-order finite element pairs not satisfying the inf-sup condition. Under the uniqueness condition, error analyses for both algorithms are given. Numerical results are reported to demonstrate the effectiveness of the simple and Newton scheme.展开更多
On triangle or quadrilateral meshes, two finite element methods are proposed for solving the Reissner-Mindlin plate problem either by augmenting the Galerkin formulation or modifying the plate-thickness. In these meth...On triangle or quadrilateral meshes, two finite element methods are proposed for solving the Reissner-Mindlin plate problem either by augmenting the Galerkin formulation or modifying the plate-thickness. In these methods, the transverse displacement is approximated by conforming （bi）linear macroelements or （bi）quadratic elements, and the rotation by conforming （bi）linear elements. The shear stress can be locally computed from transverse displacement and rotation. Uniform in plate thickness, optimal error bounds are obtained for the transverse displacement, rotation, and shear stress in their natural norms. Numerical results are presented to illustrate the theoretical results.展开更多
This article presents a micro-structure tensor enhanced elasto-plastic finite element(FE)method to address strength anisotropy in three-dimensional(3D)soil slope stability analysis.The gravity increase method(GIM)is e...This article presents a micro-structure tensor enhanced elasto-plastic finite element(FE)method to address strength anisotropy in three-dimensional(3D)soil slope stability analysis.The gravity increase method(GIM)is employed to analyze the stability of 3D anisotropic soil slopes.The accuracy of the proposed method is first verified against the data in the literature.We then simulate the 3D soil slope with a straight slope surface and the convex and concave slope surfaces with a 90turning corner to study the 3D effect on slope stability and the failure mechanism under anisotropy conditions.Based on our numerical results,the end effect significantly impacts the failure mechanism and safety factor.Anisotropy degree notably affects the safety factor,with higher degrees leading to deeper landslides.For concave slopes,they can be approximated by straight slopes with suitable boundary conditions to assess their stability.Furthermore,a case study of the Saint-Alban test embankment A in Quebec,Canada,is provided to demonstrate the applicability of the proposed FE model.展开更多
Wellbore breakout is one of the critical issues in drilling due to the fact that the related problems result in additional costs and impact the drilling scheme severely.However,the majority of such wellbore breakout a...Wellbore breakout is one of the critical issues in drilling due to the fact that the related problems result in additional costs and impact the drilling scheme severely.However,the majority of such wellbore breakout analyses were based on continuum mechanics.In addition to failure in intact rocks,wellbore breakouts can also be initiated along natural discontinuities,e.g.weak planes and fractures.Furthermore,the conventional models in wellbore breakouts with uniform distribution fractures could not reflect the real drilling situation.This paper presents a fully coupled hydro-mechanical model of the SB-X well in the Tarim Basin,China for evaluating wellbore breakouts in heavily fractured rocks under anisotropic stress states using the distinct element method(DEM)and the discrete fracture network(DFN).The developed model was validated against caliper log measurement,and its stability study was carried out by stress and displacement analyses.A parametric study was performed to investigate the effects of the characteristics of fracture distribution(orientation and length)on borehole stability by sensitivity studies.Simulation results demonstrate that the increase of the standard deviation of orientation when the fracture direction aligns parallel or perpendicular to the principal stress direction aggravates borehole instability.Moreover,an elevation in the average fracture length causes the borehole failure to change from the direction of the minimum in-situ horizontal principal stress(i.e.the direction of wellbore breakouts)towards alternative directions,ultimately leading to the whole wellbore failure.These findings provide theoretical insights for predicting wellbore breakouts in heavily fractured rocks.展开更多
The demagnetization curves were calculated using micromagnetic finite-element method for nanocomposite Pr 2Fe 14B/α-Fe permanent magnets with precipitate-typed microstructure. Due to intergrain exchange coupling,...The demagnetization curves were calculated using micromagnetic finite-element method for nanocomposite Pr 2Fe 14B/α-Fe permanent magnets with precipitate-typed microstructure. Due to intergrain exchange coupling, both remanence enhancement and a single magnetic phase behavior in demagnetization curve were found. For the samples with the hard phase as the precipitate and the soft one as the matrix, a coercivity μ 0H c of 0.78 T, a remanence J r of 1.18 T and a large energy product (BH) max of 200 kJ·m -3 are obtained for the sample with hard grain size being 23 nm. While, for the sample with the soft phase as the precipitate, μ 0H c of 0.95 T, J r of 1.24 T and (BH) max of 240 kJ·m -3 are obtained for the sample with soft grain size being 10 nm. The calculating results were compared with the experimental Pr 8Fe 87B 5 ribbons. The dependence of remanence and coercivity on the microstructure was discussed extensively.展开更多
A unified analysis is presented for the stabilized methods including the pres- sure projection method and the pressure gradient local projection method of conforming and nonconforming low-order mixed finite elements f...A unified analysis is presented for the stabilized methods including the pres- sure projection method and the pressure gradient local projection method of conforming and nonconforming low-order mixed finite elements for the stationary Navier-Stokes equa- tions. The existence and uniqueness of the solution and the optimal error estimates are proved.展开更多
In this paper,a new strategy for a sub-element-based shock capturing for discontinuous Galerkin(DG)approximations is presented.The idea is to interpret a DG element as a col-lection of data and construct a hierarchy o...In this paper,a new strategy for a sub-element-based shock capturing for discontinuous Galerkin(DG)approximations is presented.The idea is to interpret a DG element as a col-lection of data and construct a hierarchy of low-to-high-order discretizations on this set of data,including a first-order finite volume scheme up to the full-order DG scheme.The dif-ferent DG discretizations are then blended according to sub-element troubled cell indicators,resulting in a final discretization that adaptively blends from low to high order within a single DG element.The goal is to retain as much high-order accuracy as possible,even in simula-tions with very strong shocks,as,e.g.,presented in the Sedov test.The framework retains the locality of the standard DG scheme and is hence well suited for a combination with adaptive mesh refinement and parallel computing.The numerical tests demonstrate the sub-element adaptive behavior of the new shock capturing approach and its high accuracy.展开更多
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.展开更多
Chaotic vibrations of flexible non-linear Euler-Bernoulli beams subjected to harmonic load and with various boundary conditions（symmetric and non-symmetric）are studied in this work.Reliability of the obtained result...Chaotic vibrations of flexible non-linear Euler-Bernoulli beams subjected to harmonic load and with various boundary conditions（symmetric and non-symmetric）are studied in this work.Reliability of the obtained results is verified by the finite difference method（FDM）and the finite element method（FEM）with the Bubnov-Galerkin approximation for various boundary conditions and various dynamic regimes（regular and non-regular）.The influence of boundary conditions on the Euler-Bernoulli beams dynamics is studied mainly,dynamic behavior vs.control parameters { ωp,q0 } is reported,and scenarios of the system transition into chaos are illustrated.展开更多
This paper presents a novel framework for stochastic analysis of linear elastic fracture problems.Monte Carlo simulation(MCs)is adopted to address the multi-dimensional uncertainties,whose computation cost is reduced ...This paper presents a novel framework for stochastic analysis of linear elastic fracture problems.Monte Carlo simulation(MCs)is adopted to address the multi-dimensional uncertainties,whose computation cost is reduced by combination of Proper Orthogonal Decomposition(POD)and the Radial Basis Function(RBF).In order to avoid re-meshing and retain the geometric exactness,isogeometric boundary element method(IGABEM)is employed for simulation,in which the Non-Uniform Rational B-splines(NURBS)are employed for representing the crack surfaces and discretizing dual boundary integral equations.The stress intensity factors(SIFs)are extracted by M integral method.The numerical examples simulate several cracked structures with various uncertain parameters such as load effects,materials,geometric dimensions,and the results are verified by comparison with the analytical solutions.展开更多
Since 1930, the analysis of slope stability is done according to the limit equilibrium approach. Several methods were developed of which certain remain applicable because of their simplicity. However, major disadvanta...Since 1930, the analysis of slope stability is done according to the limit equilibrium approach. Several methods were developed of which certain remain applicable because of their simplicity. However, major disadvantages of these methods are (1) they do not take into account the soil behavior and (2) the complex cases cannot be studied with precision. The use of the finite elements in calculations of stability has to overcome the weakness of the traditional methods. An analysis of stability was applied to a slope, of complex geometry, composed of alternating sandstone and marls using finite elements and limit equilibrium methods. The calculation of the safety factors did not note any significant difference between the two approaches. Various calculations carried out illustrate perfectly benefits that can be gained from modeling the behavior by the finite elements method. In the finite elements analysis, the shape of deformations localization in the slope is nearly circular and confirms the shape of the failure line which constitutes the basic assumption of the analytical methods. The integration of the constitutive laws of soils and the use of field’s results tests in finite elements models predict the failure mode, to better approach the real behavior of slope soil formations and to optimize its reinforcement.展开更多
基金supported by the National Natural Science Foundation of China(12401482)the second author was supported by the National Natural Science Foundation of China(12371371,12261160361,11971366)supported by the Open Research Fund of Hubei Key Laboratory of Computational Science,Wuhan University.
文摘A new quadrilateral edge element method is proposed and analyzed for Maxwell equations.This proposed method is based on Duan-Liang quadrilateral element(Math.Comp.73(2004),pp.1–18).When applied to the eigenvalue problem,the method is spectral-correct and spurious-free.Stability and error estimates are obtained,including the interpolation error estimates and the error estimates between the finite element solution and the exact solution.The method is suitable for singular solution as well as smooth solution,and consequently,the method is valid for nonconvex domains which may have a number of reentrant corners.Of course,the method is suitable for arbitrary quadrilaterals(under the usual shape-regular condition).
基金supported by the Shanxi Scholarship Council of China(Grant No.2023-036)the Natural Science Foundation of Shanxi Province(Grant No.202303021222020).
文摘This study explores a sensitivity analysis method based on the boundary element method(BEM)to address the computational complexity in acoustic analysis with ground reflection problems.The advantages of BEM in acoustic simulations and its high computational cost in broadband problems are examined.To improve efficiency,a Taylor series expansion is applied to decouple frequency-dependent terms in BEM.Additionally,the SecondOrder Arnoldi(SOAR)model order reduction method is integrated to reduce computational costs and enhance numerical stability.Furthermore,an isogeometric sensitivity boundary integral equation is formulated using the direct differentiation method,incorporating Cauchy principal value integrals and Hadamard finite part integrals to handle singularities.The proposed method improves the computational efficiency,and the acoustic sensitivity analysis provides theoretical support for further acoustic structure optimization.
基金the financial support provided by MHRD,Govt.of IndiaCoal India Limited for providing financial assistance for the research(Project No.CIL/R&D/01/73/2021)the partial financial support provided by the Ministry of Education,Government of India,under SPARC project(Project No.P1207)。
文摘Sudden and unforeseen seismic failures of coal mine overburden(OB)dump slopes interrupt mining operations,cause loss of lives and delay the production of coal.Consideration of the spatial heterogeneity of OB dump materials is imperative for an adequate evaluation of the seismic stability of OB dump slopes.In this study,pseudo-static seismic stability analyses are carried out for an OB dump slope by considering the material parameters obtained from an insitu field investigation.Spatial heterogeneity is simulated through use of the random finite element method(RFEM)and the random limit equilibrium method(RLEM)and a comparative study is presented.Combinations of horizontal and vertical spatial correlation lengths were considered for simulating isotropic and anisotropic random fields within the OB dump slope.Seismic performances of the slope have been reported through the probability of failure and reliability index.It was observed that the RLEM approach overestimates failure probability(P_(f))by considering seismic stability with spatial heterogeneity.The P_(f)was observed to increase with an increase in the coefficient of variation of friction angle of the dump materials.Further,it was inferred that the RLEM approach may not be adequately applicable for assessing the seismic stability of an OB dump slope for a horizontal seismic coefficient that is more than or equal to 0.1.
文摘The stochastic boundary element method(SBEM)is developed in this paper for 3D problems with body forces and reliability analysis of engineering structures.The integral equations of SBEM are established by the approach of partial derivation with respect to stochastic variables,considering the yield limit,rotation speeds and material density to be the fundamental stochastic variables.Through analyzing a numerical example and a turbo-disk of an aeroengine,the results show that the method developed is successful.
基金funded by the project of the Major Scientific and Technological Projects of CNOOC in the 14th Five-Year Plan(No.KJGG2022-0701)the CNOOC Research Institute(No.2020PFS-03).
文摘To analyze the differences in the transport and distribution of different types of proppants and to address issues such as the short effective support of proppant and poor placement in hydraulically intersecting fractures,this study considered the combined impact of geological-engineering factors on conductivity.Using reservoir production parameters and the discrete elementmethod,multispherical proppants were constructed.Additionally,a 3D fracture model,based on the specified conditions of the L block,employed coupled(Computational Fluid Dynamics)CFD-DEM(Discrete ElementMethod)for joint simulations to quantitatively analyze the transport and placement patterns of multispherical proppants in intersecting fractures.Results indicate that turbulent kinetic energy is an intrinsic factor affecting proppant transport.Moreover,the efficiency of placement and migration distance of low-sphericity quartz sand constructed by the DEM in the main fracture are significantly reduced compared to spherical ceramic proppants,with a 27.7%decrease in the volume fraction of the fracture surface,subsequently affecting the placement concentration and damaging fracture conductivity.Compared to small-angle fractures,controlling artificial and natural fractures to expand at angles of 45°to 60°increases the effective support length by approximately 20.6%.During hydraulic fracturing of gas wells,ensuring the fracture support area and post-closure conductivity can be achieved by controlling the sphericity of proppants and adjusting the perforation direction to control the direction of artificial fractures.
基金supported by the National Natural Science Foundation of China (10574156)
文摘Magnetization configurations were calculated under various magnetic fields for nanocrystalline Pr-Fe-B permanent magnets by micromagnetic finite element method.According to the configurations during demagnetization process, the mechanism of magnetization reversal was analyzed.For the Pr2Fe14B with 10 nm grains or its composite with 10vol.% α-Fe, the coercivity was determined by nucleation of reversed domain that took place at grain boundaries.However, for Pr2Fe14B with 30 nm grains, coercivity was controlled by pinning of the nucle-ated domain.For Pr2Fe14B/α-Fe with 30vol.% α-Fe, the demagnetization behavior was characterized by continuous reversal of α-Fe moment.
基金Project supported by the National Natural Science Foundation of China(Nos.11671157 and11826212)
文摘In this paper,two-grid immersed finite element (IFE) algorithms are proposed and analyzed for semi-linear interface problems with discontinuous diffusion coefficients in two dimension.Because of the advantages of finite element (FE) formulation and the simple structure of Cartesian grids,the IFE discretization is used in this paper.Two-grid schemes are formulated to linearize the FE equations.It is theoretically and numerically illustrated that the coarse space can be selected as coarse as H =O(h^1/4)(or H =O(h^1/8)),and the asymptotically optimal approximation can be achieved as the nonlinear schemes.As a result,we can settle a great majority of nonlinear equations as easy as linearized problems.In order to estimate the present two-grid algorithms,we derive the optimal error estimates of the IFE solution in the L^p norm.Numerical experiments are given to verify the theorems and indicate that the present two-grid algorithms can greatly improve the computing efficiency.
基金supported by the National Natural Science Foundation of China (Grant Nos. 41772286 and 42077247)the Fundamental Research Funds for the Central Universities, China
文摘The past decade has witnessed the substantial growth in research interests and progress on the subject of coupled hydro-mechanical processes in rocks and soils,driven mainly by the surge of research in unconventional hydrocarbon reservoirs and associated hazards.Many coupling techniques have been developed to include the effects of fluid flow in the discrete element method(DEM),and the techniques have been applied to a variety of geomechanical problems.Although these coupling methods have been successfully applied in various engineering fields,no single fluid/DEM coupling method is universal due to the complexity of engineering problems and the limitations of the numerical methods.For researchers and engineers,the key to solve a specific problem is to select the most appropriate fluid/DEM coupling method among these modeling technologies.The purpose of this paper is to give a comprehensive review of fluid flow/DEM coupling methods and relevant research.Given their importance,the availability or unavailability of best practice guidelines is outlined.The theoretical background and current status of DEM are introduced first,and the principles,applications,and advantages and disadvantages of different fluid flow/DEM coupling methods are discussed.Finally,a summary with speculation on future development trends is given.
基金Project supported by the State Key Development Program for Basic Research of China (Grant No. 2006CB601007)the National Natural Science Foundation of China (Grant No. 10674006)the National High Technology Research and Development Program of China (Grant No. 2007AA03Z238)
文摘This paper discusses the forward and inverse problem for cardiac magnetic fields and electric potentials. A torso-heart model established by boundary element method （BEM） is used for studying the distributions of cardiac magnetic fields and electric potentials. Because node-to-node and triangle-to-triangle BEM can lead to discrepant field distributions, their properties and influences are compared. Then based on constructed torso-heart model and supposed current source functional model-current dipole array, the magnetic and electric imaging by optimal constrained linear inverse method are applied at the same time. Through figure and reconstructing parameter comparison, though the magnetic current dipole array imaging possesses better reconstructing effect, however node-to-node BEM and triangleto-triangle BEM make little difference to magnetic and electric imaging.
文摘In this article, on the basis of two-level discretizations and multiscale finite element method, two kinds of finite element algorithms for steady Navier-Stokes problem are presented and discussed. The main technique is first to use a standard finite element discretization on a coarse mesh to approximate low frequencies, then to apply the simple and Newton scheme to linearize discretizations on a fine grid. At this process, multiscale finite element method as a stabilized method deals with the lowest equal-order finite element pairs not satisfying the inf-sup condition. Under the uniqueness condition, error analyses for both algorithms are given. Numerical results are reported to demonstrate the effectiveness of the simple and Newton scheme.
基金supported by NSFC(11571266,91430106,11171168,11071132)NSFC-RGC(China-Hong Kong)(11661161017)
文摘On triangle or quadrilateral meshes, two finite element methods are proposed for solving the Reissner-Mindlin plate problem either by augmenting the Galerkin formulation or modifying the plate-thickness. In these methods, the transverse displacement is approximated by conforming （bi）linear macroelements or （bi）quadratic elements, and the rotation by conforming （bi）linear elements. The shear stress can be locally computed from transverse displacement and rotation. Uniform in plate thickness, optimal error bounds are obtained for the transverse displacement, rotation, and shear stress in their natural norms. Numerical results are presented to illustrate the theoretical results.
基金supported by the National Natural Science Foundation of China(Grant Nos.51890912,51979025 and 52011530189).
文摘This article presents a micro-structure tensor enhanced elasto-plastic finite element(FE)method to address strength anisotropy in three-dimensional(3D)soil slope stability analysis.The gravity increase method(GIM)is employed to analyze the stability of 3D anisotropic soil slopes.The accuracy of the proposed method is first verified against the data in the literature.We then simulate the 3D soil slope with a straight slope surface and the convex and concave slope surfaces with a 90turning corner to study the 3D effect on slope stability and the failure mechanism under anisotropy conditions.Based on our numerical results,the end effect significantly impacts the failure mechanism and safety factor.Anisotropy degree notably affects the safety factor,with higher degrees leading to deeper landslides.For concave slopes,they can be approximated by straight slopes with suitable boundary conditions to assess their stability.Furthermore,a case study of the Saint-Alban test embankment A in Quebec,Canada,is provided to demonstrate the applicability of the proposed FE model.
基金supported by National Natural Science Foundation of China(Grant Nos.52074312 and 52211530097)CNPC Science and Technology Innovation Foundation(Grant No.2021DQ02-0505).
文摘Wellbore breakout is one of the critical issues in drilling due to the fact that the related problems result in additional costs and impact the drilling scheme severely.However,the majority of such wellbore breakout analyses were based on continuum mechanics.In addition to failure in intact rocks,wellbore breakouts can also be initiated along natural discontinuities,e.g.weak planes and fractures.Furthermore,the conventional models in wellbore breakouts with uniform distribution fractures could not reflect the real drilling situation.This paper presents a fully coupled hydro-mechanical model of the SB-X well in the Tarim Basin,China for evaluating wellbore breakouts in heavily fractured rocks under anisotropic stress states using the distinct element method(DEM)and the discrete fracture network(DFN).The developed model was validated against caliper log measurement,and its stability study was carried out by stress and displacement analyses.A parametric study was performed to investigate the effects of the characteristics of fracture distribution(orientation and length)on borehole stability by sensitivity studies.Simulation results demonstrate that the increase of the standard deviation of orientation when the fracture direction aligns parallel or perpendicular to the principal stress direction aggravates borehole instability.Moreover,an elevation in the average fracture length causes the borehole failure to change from the direction of the minimum in-situ horizontal principal stress(i.e.the direction of wellbore breakouts)towards alternative directions,ultimately leading to the whole wellbore failure.These findings provide theoretical insights for predicting wellbore breakouts in heavily fractured rocks.
文摘The demagnetization curves were calculated using micromagnetic finite-element method for nanocomposite Pr 2Fe 14B/α-Fe permanent magnets with precipitate-typed microstructure. Due to intergrain exchange coupling, both remanence enhancement and a single magnetic phase behavior in demagnetization curve were found. For the samples with the hard phase as the precipitate and the soft one as the matrix, a coercivity μ 0H c of 0.78 T, a remanence J r of 1.18 T and a large energy product (BH) max of 200 kJ·m -3 are obtained for the sample with hard grain size being 23 nm. While, for the sample with the soft phase as the precipitate, μ 0H c of 0.95 T, J r of 1.24 T and (BH) max of 240 kJ·m -3 are obtained for the sample with soft grain size being 10 nm. The calculating results were compared with the experimental Pr 8Fe 87B 5 ribbons. The dependence of remanence and coercivity on the microstructure was discussed extensively.
基金supported by the National Natural Science Foundation of China(Nos.11271273 and 11271298)
文摘A unified analysis is presented for the stabilized methods including the pres- sure projection method and the pressure gradient local projection method of conforming and nonconforming low-order mixed finite elements for the stationary Navier-Stokes equa- tions. The existence and uniqueness of the solution and the optimal error estimates are proved.
文摘In this paper,a new strategy for a sub-element-based shock capturing for discontinuous Galerkin(DG)approximations is presented.The idea is to interpret a DG element as a col-lection of data and construct a hierarchy of low-to-high-order discretizations on this set of data,including a first-order finite volume scheme up to the full-order DG scheme.The dif-ferent DG discretizations are then blended according to sub-element troubled cell indicators,resulting in a final discretization that adaptively blends from low to high order within a single DG element.The goal is to retain as much high-order accuracy as possible,even in simula-tions with very strong shocks,as,e.g.,presented in the Sedov test.The framework retains the locality of the standard DG scheme and is hence well suited for a combination with adaptive mesh refinement and parallel computing.The numerical tests demonstrate the sub-element adaptive behavior of the new shock capturing approach and its high accuracy.
文摘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.
文摘Chaotic vibrations of flexible non-linear Euler-Bernoulli beams subjected to harmonic load and with various boundary conditions（symmetric and non-symmetric）are studied in this work.Reliability of the obtained results is verified by the finite difference method（FDM）and the finite element method（FEM）with the Bubnov-Galerkin approximation for various boundary conditions and various dynamic regimes（regular and non-regular）.The influence of boundary conditions on the Euler-Bernoulli beams dynamics is studied mainly,dynamic behavior vs.control parameters { ωp,q0 } is reported,and scenarios of the system transition into chaos are illustrated.
基金The authors thank the financial support of National Natural Science Foundation of China(NSFC)under Grant(Nos.51904202,11902212,11901578).
文摘This paper presents a novel framework for stochastic analysis of linear elastic fracture problems.Monte Carlo simulation(MCs)is adopted to address the multi-dimensional uncertainties,whose computation cost is reduced by combination of Proper Orthogonal Decomposition(POD)and the Radial Basis Function(RBF).In order to avoid re-meshing and retain the geometric exactness,isogeometric boundary element method(IGABEM)is employed for simulation,in which the Non-Uniform Rational B-splines(NURBS)are employed for representing the crack surfaces and discretizing dual boundary integral equations.The stress intensity factors(SIFs)are extracted by M integral method.The numerical examples simulate several cracked structures with various uncertain parameters such as load effects,materials,geometric dimensions,and the results are verified by comparison with the analytical solutions.
文摘Since 1930, the analysis of slope stability is done according to the limit equilibrium approach. Several methods were developed of which certain remain applicable because of their simplicity. However, major disadvantages of these methods are (1) they do not take into account the soil behavior and (2) the complex cases cannot be studied with precision. The use of the finite elements in calculations of stability has to overcome the weakness of the traditional methods. An analysis of stability was applied to a slope, of complex geometry, composed of alternating sandstone and marls using finite elements and limit equilibrium methods. The calculation of the safety factors did not note any significant difference between the two approaches. Various calculations carried out illustrate perfectly benefits that can be gained from modeling the behavior by the finite elements method. In the finite elements analysis, the shape of deformations localization in the slope is nearly circular and confirms the shape of the failure line which constitutes the basic assumption of the analytical methods. The integration of the constitutive laws of soils and the use of field’s results tests in finite elements models predict the failure mode, to better approach the real behavior of slope soil formations and to optimize its reinforcement.
