The moving morphable component(MMC)topology optimization method,as a typical explicit topology optimization method,has been widely concerned.In the MMC topology optimization framework,the surrogate material model is m...The moving morphable component(MMC)topology optimization method,as a typical explicit topology optimization method,has been widely concerned.In the MMC topology optimization framework,the surrogate material model is mainly used for finite element analysis at present,and the effectiveness of the surrogate material model has been fully confirmed.However,there are some accuracy problems when dealing with boundary elements using the surrogate material model,which will affect the topology optimization results.In this study,a boundary element reconstruction(BER)model is proposed based on the surrogate material model under the MMC topology optimization framework to improve the accuracy of topology optimization.The proposed BER model can reconstruct the boundary elements by refining the local meshes and obtaining new nodes in boundary elements.Then the density of boundary elements is recalculated using the new node information,which is more accurate than the original model.Based on the new density of boundary elements,the material properties and volume information of the boundary elements are updated.Compared with other finite element analysis methods,the BER model is simple and feasible and can improve computational accuracy.Finally,the effectiveness and superiority of the proposed method are verified by comparing it with the optimization results of the original surrogate material model through several numerical examples.展开更多
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).展开更多
In this work,the Hierarchical Quadrature Element Method(HQEM)formulation of geometrically exact shells is proposed and applied for geometrically nonlinear analyses of both isotropic and laminated shells.The stress res...In this work,the Hierarchical Quadrature Element Method(HQEM)formulation of geometrically exact shells is proposed and applied for geometrically nonlinear analyses of both isotropic and laminated shells.The stress resultant formulation is developed within the HQEM framework,consequently significantly simplifying the computations of residual force and stiffness matrix.The present formulation inherently avoids shear and membrane locking,benefiting from its high-order approximation property.Furthermore,HQEM’s independent nodal distribution capability conveniently supports local p-refinement and flexibly facilitates mesh generation in various structural configurations through the combination of quadrilateral and triangular elements.Remarkably,in lateral buckling analysis,the HQEM outperforms the weak-form quadrilateral element(QEM)in accuracy with identical nodal degrees of freedom(three displacements and two rotations).Under high-load nonlinear response,the QEM exhibits a maximum relative deviation of approximately 9.5%from the reference,while the HQEM remains closely aligned with the benchmark results.In addition,for the cantilever beam under tip moment,HQEM produces virtually no out-of-plane deviation,compared to a slight deviation of 0.00001 with QEM,confirming its superior numerical reliability.In summary,the method demonstrates high accuracy,superior convergence,and robustness in handling large rotations and complex post-buckling behaviors across a series of benchmark problems.展开更多
The analysis of the dynamics of surface girders is of great importance in the design of engineering structures such as steel welded bridge plane girders or concrete plate-column structures.This work is an extension of...The analysis of the dynamics of surface girders is of great importance in the design of engineering structures such as steel welded bridge plane girders or concrete plate-column structures.This work is an extension of the classical deterministic problem of free vibrations of thin(Kirchhoff)plates.Themain aim of this work is the study of stochastic eigenvibrations of thin(Kirchhoff)elastic plates resting on internal continuous and column supports by the Boundary Element Method(BEM).This work is a continuation of previous research related to the random approach in plate analysis using the BEM.The static fundamental solution(Green’s function)is applied,coupled with a nonsingular formulation of the boundary and domain integral equations.These are derived using a modified and simplified formulation of the boundary conditions,inwhich there is no need to introduce theKirchhoff forces on a plate boundary.The role of the Kirchhoff corner forces is played by the boundary elements placed close to a single corner.Internal column or linear continuous supports are introduced using the Bezine technique,where the additional collocation points are introduced inside a plate domain.This allows for significant simplification of the BEM computational algorithm.An application of the polynomial approximations in the Least Squares Method(LSM)recovery of the structural response is done.The probabilistic analysis will employ three independent computational approaches:semi-analytical method(SAM),stochastic perturbation technique(SPT),and Monte-Carlo simulations.Numerical investigations include the fundamental eigenfrequencies of an elastic,thin,homogeneous,and isotropic plate.展开更多
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.展开更多
Hydraulic fracturing plays a critical role in enhancing shale gas production in deep shale reservoirs.Conventional hydraulic fracturing simulation methods rely on prefabricated grids,which can be hindered by the chall...Hydraulic fracturing plays a critical role in enhancing shale gas production in deep shale reservoirs.Conventional hydraulic fracturing simulation methods rely on prefabricated grids,which can be hindered by the challenge of being computationally overpowered.This study proposes an efficient fracturing simulator to analyze fracture morphology during hydraulic fracturing processes in deep shale gas reservoirs.The simulator integrates the boundary element displacement discontinuity method and the finite volume method to model the fluid-solid coupling process by employing a pseudo-3D fracture model to calculate the fracture height.In particular,the Broyden iteration method was introduced to improve the computational efficiency and model robustness;it achieved a 46.6%reduction in computation time compared to the Newton-Raphson method.The influences of horizontal stress differences,natural fracture density,and natural fracture angle on the modified zone of the reservoir were simulated,and the following results were observed.(1)High stress difference reservoirs have smaller stimulated reservoir area than low stress difference reservoirs.(2)A higher natural fracture angle resulted in larger modification zones at low stress differences,while the effect of a natural fracture angle at high stress differences was not significant.(3)High-density and long natural fracture zones played a significant role in enhancing the stimulated reservoir area.These findings are critical for comprehending the impact of geological parameters on deep shale reservoirs.展开更多
To address the issue that traditional finite element methods cannot fully consider the semi-infinite earth strata and have lower solution accuracy,a new equivalent force model for induced deformation during oil and ga...To address the issue that traditional finite element methods cannot fully consider the semi-infinite earth strata and have lower solution accuracy,a new equivalent force model for induced deformation during oil and gas reservoir development is derived from the perspective of semi-infinite strata.A brand-new volume boundary element numerical method solution has been developed and verified and tested.The influences of internal flow and flow boundary of the reservoir on strata deformation are equivalent to the impacts on strata deformation when external forces act at the interior and boundary of the reservoir,respectively.Calculation methods for the flow equivalent force and boundary equivalent force are provided.The deformation solution at any point in the strata can be obtained through the convolution of flow equivalent forces,boundary equivalent forces and Green’s functions.After discretization,the deformation solution at any point in the strata can be obtained by multiplying the grid boundary equivalent forces,grid flow equivalent forces with their corresponding grid boundary sources and grid volume sources respectively,and then summing them up.This numerical method is termed the Volumetric Boundary Element Method(VBEM).Compared with traditional commercial simulators,VBEM fully considers the effects of reservoir flow boundaries,pore pressure gradient fields within the reservoir,and fluid mass changes within pores on formation deformation.It eliminates the need for meshing outside the reservoir,achieves significantly improved solution accuracy,and provides a new technical framework for simulating deformation induced by reservoir development.展开更多
The finite element artificial transmitting boundary method is employed here to treat the near field scattering of a cylindrical wave from an irregular cylinder. A comparison is made between this method and the analy...The finite element artificial transmitting boundary method is employed here to treat the near field scattering of a cylindrical wave from an irregular cylinder. A comparison is made between this method and the analytical one. And then examples are given to demonstrate the solution of several problems of the irregular object scattering. The method can not only produce clear physical pictures, but can efficiently handle many complicated scattering problems.展开更多
In this paper, an identification method to estimate the unbalances is introduced, which is based on the boundary element method (BEM). By using the vibration response measured at some points on the flexible rotor the ...In this paper, an identification method to estimate the unbalances is introduced, which is based on the boundary element method (BEM). By using the vibration response measured at some points on the flexible rotor the unbalances can be identified conveniently. Therefore, the rotor can be balanced without test runs.展开更多
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.展开更多
In the present paper, we examine the performance of an efficient type of wave-absorbing porous marine structure under the attack of regular oblique waves by using a Multi-Domain Boundary Element Method(MDBEM). The str...In the present paper, we examine the performance of an efficient type of wave-absorbing porous marine structure under the attack of regular oblique waves by using a Multi-Domain Boundary Element Method(MDBEM). The structure consists of two perforated vertical thin barriers creating what can be called a wave absorbing chamber system. The barriers are surface piercing, thereby eliminating wave overtopping. The problem of the interaction of obliquely incident linear waves upon a pair of perforated barriers is first formulated in the context of linear diffraction theory. The resulting boundary integral equation, which is matched with far-field solutions presented in terms of analytical series with unknown coefficients, as well as the appropriate boundary conditions at the free surface, seabed, and barriers, is then solved numerically using MDBEM. Dissipation of the wave energy due to the presence of the perforated barriers is represented by a simple yet effective relation in terms of the porosity parameter appropriate for thin perforated walls. The results are presented in terms of reflection and transmission coefficients. The effects of the incident wave angles, relative water depths, porosities, depths of the walls, and other major parameters of interest are explored.展开更多
A high order boundary element method was developed for the complex velocity potential problem. The method ensures not only the continuity of the potential at the nodes of each element but also the velocity. It can be ...A high order boundary element method was developed for the complex velocity potential problem. The method ensures not only the continuity of the potential at the nodes of each element but also the velocity. It can be applied to a variety of velocity potential problems. The present paper, however, focused on its application to the problem of water entry of a wedge with varying speed. The continuity of the velocity achieved herein is particularly important for this kind of nonlinear free surface flow problem, because when the time stepping method is used, the free surface is updated through the velocity obtained at each node and the accuracy of the velocity is therefore crucial. Calculation was made for a case when the distance S that the wedge has travelled and time t follow the relationship s=Dtα, where D and α are constants, which is found to lead to a self similar flow field when the effect due to gravity is ignored.展开更多
For higher accuracy in simulating the transformation of three dimensional waves, in consideration of the advantages of constant panels and linear elements, a combined boundary elements is applied in this research. The...For higher accuracy in simulating the transformation of three dimensional waves, in consideration of the advantages of constant panels and linear elements, a combined boundary elements is applied in this research. The method can be used to remove the transverse vibration due to the accumulation of computational errors. A combined boundary condition of sponge layer and Sommerfeld radiation condition is used to remove the reflected waves from the computing domain. By following the water particle on the water surface, the third order Stokes wave transform is simulated by the numerical wave flume technique. The computed results are in good agreement with theoretical ones.展开更多
In this study,a comprehensive parametric analysis was performed on non-uniform excitation of V-shaped topography using the boundary element method in time domain.For this purpose,wave scattering analysis was carried o...In this study,a comprehensive parametric analysis was performed on non-uniform excitation of V-shaped topography using the boundary element method in time domain.For this purpose,wave scattering analysis was carried out on a topography subjected to the SV-wave for different predominant frequencies and shape ratios.Based on the numerical results,new coherence and time delay functions are proposed to generate non-uniform ground motion for topographic irregularities.The efficiency and accuracy of the proposed functions for real engineering problems are indicated by comparison with observations reported in previous literature.展开更多
A scheme of boundary element method for moving contact of two-dimensional elastic bodies using conforming discretization is presented. Both the displacement and the traction boundary conditions are satisfied on the co...A scheme of boundary element method for moving contact of two-dimensional elastic bodies using conforming discretization is presented. Both the displacement and the traction boundary conditions are satisfied on the contacting region in the sense of discretization. An algorithm to deal with the moving of the contact boundary on a larger possible contact region is presented. The algorithm is generalized to rolling contact problem as well. Some numerical examples of moving and rolling contact of 2D elastic bodies with or without friction, including the bodies with a hole-type defect, are given to show the effectiveness and the accuracy of the presented schemes.展开更多
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.展开更多
This paper presents an isogeometric boundary element method(IGABEM)for transient heat conduction analysis.The Non-Uniform Rational B-spline(NURBS)basis functions,which are used to construct the geometry of the structu...This paper presents an isogeometric boundary element method(IGABEM)for transient heat conduction analysis.The Non-Uniform Rational B-spline(NURBS)basis functions,which are used to construct the geometry of the structures,are employed to discretize the physical unknowns in the boundary integral formulations of the governing equations.Bezier extraction technique is employed to accelerate the evaluation of NURBS basis functions.We adopt a radial integration method to address the additional domain integrals.The numerical examples demonstrate the advantage of IGABEM in dimension reduction and the seamless connection between CAD and numerical analysis.展开更多
A cardiac vector model is presented and verified, and then the forward problem for cardiac magnetic fields and electric potential are discussed based on this model and the realistic human torso volume conductor model,...A cardiac vector model is presented and verified, and then the forward problem for cardiac magnetic fields and electric potential are discussed based on this model and the realistic human torso volume conductor model, including lungs. A torso-cardiac vector model is used for a 12-lead electrocardiographic （ECG） and magneto-cardiogram （MCG） simulation study by using the boundary element method （BEM）. Also, we obtain the MCG wave picture using a compound four-channel HTc.SQUID system in a magnetically shielded room. By comparing the simulated results and experimental results, we verify the cardiac vector model and then do a preliminary study of the forward problem of MCG and ECG. Therefore, the results show that the vector model is reasonable in cardiac electrophysiology.展开更多
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.展开更多
基金supported by the Science and Technology Research Project of Henan Province(242102241055)the Industry-University-Research Collaborative Innovation Base on Automobile Lightweight of“Science and Technology Innovation in Central Plains”(2024KCZY315)the Opening Fund of State Key Laboratory of Structural Analysis,Optimization and CAE Software for Industrial Equipment(GZ2024A03-ZZU).
文摘The moving morphable component(MMC)topology optimization method,as a typical explicit topology optimization method,has been widely concerned.In the MMC topology optimization framework,the surrogate material model is mainly used for finite element analysis at present,and the effectiveness of the surrogate material model has been fully confirmed.However,there are some accuracy problems when dealing with boundary elements using the surrogate material model,which will affect the topology optimization results.In this study,a boundary element reconstruction(BER)model is proposed based on the surrogate material model under the MMC topology optimization framework to improve the accuracy of topology optimization.The proposed BER model can reconstruct the boundary elements by refining the local meshes and obtaining new nodes in boundary elements.Then the density of boundary elements is recalculated using the new node information,which is more accurate than the original model.Based on the new density of boundary elements,the material properties and volume information of the boundary elements are updated.Compared with other finite element analysis methods,the BER model is simple and feasible and can improve computational accuracy.Finally,the effectiveness and superiority of the proposed method are verified by comparing it with the optimization results of the original surrogate material model through several numerical examples.
基金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 National Natural Science Foundation of China(Grant Nos.12472194,12002018,11972004,11772031,11402015).
文摘In this work,the Hierarchical Quadrature Element Method(HQEM)formulation of geometrically exact shells is proposed and applied for geometrically nonlinear analyses of both isotropic and laminated shells.The stress resultant formulation is developed within the HQEM framework,consequently significantly simplifying the computations of residual force and stiffness matrix.The present formulation inherently avoids shear and membrane locking,benefiting from its high-order approximation property.Furthermore,HQEM’s independent nodal distribution capability conveniently supports local p-refinement and flexibly facilitates mesh generation in various structural configurations through the combination of quadrilateral and triangular elements.Remarkably,in lateral buckling analysis,the HQEM outperforms the weak-form quadrilateral element(QEM)in accuracy with identical nodal degrees of freedom(three displacements and two rotations).Under high-load nonlinear response,the QEM exhibits a maximum relative deviation of approximately 9.5%from the reference,while the HQEM remains closely aligned with the benchmark results.In addition,for the cantilever beam under tip moment,HQEM produces virtually no out-of-plane deviation,compared to a slight deviation of 0.00001 with QEM,confirming its superior numerical reliability.In summary,the method demonstrates high accuracy,superior convergence,and robustness in handling large rotations and complex post-buckling behaviors across a series of benchmark problems.
基金funded by research grant OPUS no.2021/41/B/ST8/02432 entitled Probabilistic entropy in engineering computations sponsored by The National Science Center in Polandthe Institute of Structural Analysis of Poznan University of Technology in the framework of the internal research grant 0411/SBAD/0010.
文摘The analysis of the dynamics of surface girders is of great importance in the design of engineering structures such as steel welded bridge plane girders or concrete plate-column structures.This work is an extension of the classical deterministic problem of free vibrations of thin(Kirchhoff)plates.Themain aim of this work is the study of stochastic eigenvibrations of thin(Kirchhoff)elastic plates resting on internal continuous and column supports by the Boundary Element Method(BEM).This work is a continuation of previous research related to the random approach in plate analysis using the BEM.The static fundamental solution(Green’s function)is applied,coupled with a nonsingular formulation of the boundary and domain integral equations.These are derived using a modified and simplified formulation of the boundary conditions,inwhich there is no need to introduce theKirchhoff forces on a plate boundary.The role of the Kirchhoff corner forces is played by the boundary elements placed close to a single corner.Internal column or linear continuous supports are introduced using the Bezine technique,where the additional collocation points are introduced inside a plate domain.This allows for significant simplification of the BEM computational algorithm.An application of the polynomial approximations in the Least Squares Method(LSM)recovery of the structural response is done.The probabilistic analysis will employ three independent computational approaches:semi-analytical method(SAM),stochastic perturbation technique(SPT),and Monte-Carlo simulations.Numerical investigations include the fundamental eigenfrequencies of an elastic,thin,homogeneous,and isotropic plate.
基金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.
文摘Hydraulic fracturing plays a critical role in enhancing shale gas production in deep shale reservoirs.Conventional hydraulic fracturing simulation methods rely on prefabricated grids,which can be hindered by the challenge of being computationally overpowered.This study proposes an efficient fracturing simulator to analyze fracture morphology during hydraulic fracturing processes in deep shale gas reservoirs.The simulator integrates the boundary element displacement discontinuity method and the finite volume method to model the fluid-solid coupling process by employing a pseudo-3D fracture model to calculate the fracture height.In particular,the Broyden iteration method was introduced to improve the computational efficiency and model robustness;it achieved a 46.6%reduction in computation time compared to the Newton-Raphson method.The influences of horizontal stress differences,natural fracture density,and natural fracture angle on the modified zone of the reservoir were simulated,and the following results were observed.(1)High stress difference reservoirs have smaller stimulated reservoir area than low stress difference reservoirs.(2)A higher natural fracture angle resulted in larger modification zones at low stress differences,while the effect of a natural fracture angle at high stress differences was not significant.(3)High-density and long natural fracture zones played a significant role in enhancing the stimulated reservoir area.These findings are critical for comprehending the impact of geological parameters on deep shale reservoirs.
基金Supported by the China National Natural Science Foundation Project(52274048)Beijing Natural Science Foundation Project(3222037)。
文摘To address the issue that traditional finite element methods cannot fully consider the semi-infinite earth strata and have lower solution accuracy,a new equivalent force model for induced deformation during oil and gas reservoir development is derived from the perspective of semi-infinite strata.A brand-new volume boundary element numerical method solution has been developed and verified and tested.The influences of internal flow and flow boundary of the reservoir on strata deformation are equivalent to the impacts on strata deformation when external forces act at the interior and boundary of the reservoir,respectively.Calculation methods for the flow equivalent force and boundary equivalent force are provided.The deformation solution at any point in the strata can be obtained through the convolution of flow equivalent forces,boundary equivalent forces and Green’s functions.After discretization,the deformation solution at any point in the strata can be obtained by multiplying the grid boundary equivalent forces,grid flow equivalent forces with their corresponding grid boundary sources and grid volume sources respectively,and then summing them up.This numerical method is termed the Volumetric Boundary Element Method(VBEM).Compared with traditional commercial simulators,VBEM fully considers the effects of reservoir flow boundaries,pore pressure gradient fields within the reservoir,and fluid mass changes within pores on formation deformation.It eliminates the need for meshing outside the reservoir,achieves significantly improved solution accuracy,and provides a new technical framework for simulating deformation induced by reservoir development.
文摘The finite element artificial transmitting boundary method is employed here to treat the near field scattering of a cylindrical wave from an irregular cylinder. A comparison is made between this method and the analytical one. And then examples are given to demonstrate the solution of several problems of the irregular object scattering. The method can not only produce clear physical pictures, but can efficiently handle many complicated scattering problems.
文摘In this paper, an identification method to estimate the unbalances is introduced, which is based on the boundary element method (BEM). By using the vibration response measured at some points on the flexible rotor the unbalances can be identified conveniently. Therefore, the rotor can be balanced without test runs.
文摘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.
文摘In the present paper, we examine the performance of an efficient type of wave-absorbing porous marine structure under the attack of regular oblique waves by using a Multi-Domain Boundary Element Method(MDBEM). The structure consists of two perforated vertical thin barriers creating what can be called a wave absorbing chamber system. The barriers are surface piercing, thereby eliminating wave overtopping. The problem of the interaction of obliquely incident linear waves upon a pair of perforated barriers is first formulated in the context of linear diffraction theory. The resulting boundary integral equation, which is matched with far-field solutions presented in terms of analytical series with unknown coefficients, as well as the appropriate boundary conditions at the free surface, seabed, and barriers, is then solved numerically using MDBEM. Dissipation of the wave energy due to the presence of the perforated barriers is represented by a simple yet effective relation in terms of the porosity parameter appropriate for thin perforated walls. The results are presented in terms of reflection and transmission coefficients. The effects of the incident wave angles, relative water depths, porosities, depths of the walls, and other major parameters of interest are explored.
文摘A high order boundary element method was developed for the complex velocity potential problem. The method ensures not only the continuity of the potential at the nodes of each element but also the velocity. It can be applied to a variety of velocity potential problems. The present paper, however, focused on its application to the problem of water entry of a wedge with varying speed. The continuity of the velocity achieved herein is particularly important for this kind of nonlinear free surface flow problem, because when the time stepping method is used, the free surface is updated through the velocity obtained at each node and the accuracy of the velocity is therefore crucial. Calculation was made for a case when the distance S that the wedge has travelled and time t follow the relationship s=Dtα, where D and α are constants, which is found to lead to a self similar flow field when the effect due to gravity is ignored.
基金National Natural Science Foundation of China(No.49876026)
文摘For higher accuracy in simulating the transformation of three dimensional waves, in consideration of the advantages of constant panels and linear elements, a combined boundary elements is applied in this research. The method can be used to remove the transverse vibration due to the accumulation of computational errors. A combined boundary condition of sponge layer and Sommerfeld radiation condition is used to remove the reflected waves from the computing domain. By following the water particle on the water surface, the third order Stokes wave transform is simulated by the numerical wave flume technique. The computed results are in good agreement with theoretical ones.
文摘In this study,a comprehensive parametric analysis was performed on non-uniform excitation of V-shaped topography using the boundary element method in time domain.For this purpose,wave scattering analysis was carried out on a topography subjected to the SV-wave for different predominant frequencies and shape ratios.Based on the numerical results,new coherence and time delay functions are proposed to generate non-uniform ground motion for topographic irregularities.The efficiency and accuracy of the proposed functions for real engineering problems are indicated by comparison with observations reported in previous literature.
基金The project supported by the National Natural Science Foundation of China (19772025)
文摘A scheme of boundary element method for moving contact of two-dimensional elastic bodies using conforming discretization is presented. Both the displacement and the traction boundary conditions are satisfied on the contacting region in the sense of discretization. An algorithm to deal with the moving of the contact boundary on a larger possible contact region is presented. The algorithm is generalized to rolling contact problem as well. Some numerical examples of moving and rolling contact of 2D elastic bodies with or without friction, including the bodies with a hole-type defect, are given to show the effectiveness and the accuracy of the presented schemes.
基金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.
基金funded by National Natural Science Foundation of China(NSFC)under Grant Nos.11702238,51904202,and 11902212Nanhu Scholars Program for Young Scholars of XYNU.
文摘This paper presents an isogeometric boundary element method(IGABEM)for transient heat conduction analysis.The Non-Uniform Rational B-spline(NURBS)basis functions,which are used to construct the geometry of the structures,are employed to discretize the physical unknowns in the boundary integral formulations of the governing equations.Bezier extraction technique is employed to accelerate the evaluation of NURBS basis functions.We adopt a radial integration method to address the additional domain integrals.The numerical examples demonstrate the advantage of IGABEM in dimension reduction and the seamless connection between CAD and numerical analysis.
基金supported by the State Key Development Program for Basic Research of China (Grant No. 2011CBA00106)the National Natural Science Foundation of China (Grant Nos. 10674006, 81171421, and 61101046)the National High Technology Research and Development Program of China (Grant No. 2007AA03Z238)
文摘A cardiac vector model is presented and verified, and then the forward problem for cardiac magnetic fields and electric potential are discussed based on this model and the realistic human torso volume conductor model, including lungs. A torso-cardiac vector model is used for a 12-lead electrocardiographic （ECG） and magneto-cardiogram （MCG） simulation study by using the boundary element method （BEM）. Also, we obtain the MCG wave picture using a compound four-channel HTc.SQUID system in a magnetically shielded room. By comparing the simulated results and experimental results, we verify the cardiac vector model and then do a preliminary study of the forward problem of MCG and ECG. Therefore, the results show that the vector model is reasonable in cardiac electrophysiology.
基金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.
