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.展开更多
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.展开更多
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.展开更多
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 computational framework for parachute inflation is developed based on the immersed boundary/finite element approach within the open-source IBAMR library.The fluid motion is solved by Peskin's diffuse-interface i...A computational framework for parachute inflation is developed based on the immersed boundary/finite element approach within the open-source IBAMR library.The fluid motion is solved by Peskin's diffuse-interface immersed boundary(IB)method,which is attractive for simulating moving-boundary flows with large deformations.The adaptive mesh refinement technique is employed to reduce the computational cost while retain the desired resolution.The dynamic response of the parachute is solved with the finite element approach.The canopy and cables of the parachute system are modeled with the hyperelastic material.A tether force is introduced to impose rigidity constraints for the parachute system.The accuracy and reliability of the present framework is validated by simulating inflation of a constrained square plate.Application of the present framework on several canonical cases further demonstrates its versatility for simulation of parachute inflation.展开更多
The boundary element method(BEM)is a popular method for solving acoustic wave propagation problems,especially those in exterior domains,owing to its ease in handling radiation conditions at infinity.However,BEM models...The boundary element method(BEM)is a popular method for solving acoustic wave propagation problems,especially those in exterior domains,owing to its ease in handling radiation conditions at infinity.However,BEM models must meet the requirement of 6–10 elements per wavelength,using the conventional constant,linear,or quadratic elements.Therefore,a large storage size of memory and long solution time are often needed in solving higher-frequency problems.In this work,we propose two new types of enriched elements based on conventional constant boundary elements to improve the computational efficiency of the 2D acoustic BEM.The first one uses a plane wave expansion,which can be used to model scattering problems.The second one uses a special plane wave expansion,which can be used tomodel radiation problems.Five examples are investigated to showthe advantages of the enriched elements.Compared with the conventional constant elements,the new enriched elements can deliver results with the same accuracy and in less computational time.This improvement in the computational efficiency is more evident at higher frequencies(with the nondimensional wave numbers exceeding 100).The paper concludes with the potential of our proposed enriched elements and plans for their further improvement.展开更多
The peridynamics(PD),as a promising nonlocal continuum mechanics theory,shines in solving discontinuous problems.Up to now,various numerical methods,such as the peridynamic mesh-free particlemethod(PD-MPM),peridynamic...The peridynamics(PD),as a promising nonlocal continuum mechanics theory,shines in solving discontinuous problems.Up to now,various numerical methods,such as the peridynamic mesh-free particlemethod(PD-MPM),peridynamic finite element method(PD-FEM),and peridynamic boundary element method(PD-BEM),have been proposed.PD-BEM,in particular,outperforms other methods by eliminating spurious boundary softening,efficiently handling infinite problems,and ensuring high computational accuracy.However,the existing PD-BEM is constructed exclusively for bond-based peridynamics(BBPD)with fixed Poisson’s ratio,limiting its applicability to crack propagation problems and scenarios involving infinite or semi-infinite problems.In this paper,we address these limitations by introducing the boundary element method(BEM)for ordinary state-based peridynamics(OSPD-BEM).Additionally,we present a crack propagationmodel embeddedwithin the framework ofOSPD-BEM to simulate crack propagations.To validate the effectiveness of OSPD-BEM,we conduct four numerical examples:deformation under uniaxial loading,crack initiation in a double-notched specimen,wedge-splitting test,and threepoint bending test.The results demonstrate the accuracy and efficiency of OSPD-BEM,highlighting its capability to successfully eliminate spurious boundary softening phenomena under varying Poisson’s ratios.Moreover,OSPDBEMsignificantly reduces computational time and exhibits greater consistencywith experimental results compared to PD-MPM.展开更多
In this paper,a generalized nth-order perturbation method based on the isogeometric boundary element method is proposed for the uncertainty analysis of broadband structural acoustic scattering problems.The Burton-Mill...In this paper,a generalized nth-order perturbation method based on the isogeometric boundary element method is proposed for the uncertainty analysis of broadband structural acoustic scattering problems.The Burton-Miller method is employed to solve the problem of non-unique solutions that may be encountered in the external acoustic field,and the nth-order discretization formulation of the boundary integral equation is derived.In addition,the computation of loop subdivision surfaces and the subdivision rules are introduced.In order to confirm the effectiveness of the algorithm,the computed results are contrasted and analyzed with the results under Monte Carlo simulations(MCs)through several numerical examples.展开更多
The solution to heat transfer problems in two-dimensional heterogeneous media is attended based on the scaled boundary finite element method(SBFEM)coupled with equilibrated basis functions(EqBFs).The SBFEM reduces the...The solution to heat transfer problems in two-dimensional heterogeneous media is attended based on the scaled boundary finite element method(SBFEM)coupled with equilibrated basis functions(EqBFs).The SBFEM reduces the model order by scaling the boundary solution onto the inner element.To this end,tri-lateral elements are emanated from a scaling center,followed by the development of a semi-analytical solution along the radial direction and a finite element solution along the circumferential/boundary direction.The discretization is thus limited to the boundaries of the model,and the semi-analytical radial solution is found through the solution of an eigenvalue problem,which restricts the methods’applicability to heterogeneous media.In this research,we first extracted the SBFEM formulation considering the heterogeneity of the media.Then,we replaced the semi-analytical radial solution with the EqBFs and removed the eigenvalue solution step from the SBFEM.The varying coefficients of the partial differential equation(PDE)resulting from the heterogeneity of the media are replaced by a finite series in the radial and circumferential directions of the element.A weighted residual approach is applied to the radial equation.The equilibrated radial solution series is used in the new formulation of the SBFEM.展开更多
The hydrodynamic performance of a high forward-speed ship in obliquely propagating waves is numerically examined to assess both free motions and wave field in comparison with a low forward-speed ship.This numerical mo...The hydrodynamic performance of a high forward-speed ship in obliquely propagating waves is numerically examined to assess both free motions and wave field in comparison with a low forward-speed ship.This numerical model is based on the time-domain potential flow theory and higher-order boundary element method,where an analytical expression is completely expanded to determine the base-unsteady coupling flow imposed on the moving condition of the ship.The ship in the numerical model may possess different advancing speeds,i.e.stationary,low speed,and high speed.The role of the water depth,wave height,wave period,and incident wave angle is analyzed by means of the accurate numerical model.It is found that the resonant motions of the high forward-speed ship are triggered by comparison with the stationary one.More specifically,a higher forward speed generates a V-shaped wave region with a larger elevation,which induces stronger resonant motions corresponding to larger wave periods.The shoaling effect is adverse to the motion of the low-speed ship,but is beneficial to the resonant motion of the high-speed ship.When waves obliquely propagate toward the ship,the V-shaped wave region would be broken due to the coupling effect between roll and pitch motions.It is also demonstrated that the maximum heave motion occurs in beam seas for stationary cases but occurs in head waves for high speeds.However,the variation of the pitch motion with period is hardly affected by wave incident angles.展开更多
Due to the wide application of closely spaced multi-well horizontal pads for developing unconventional gas reservoirs,interference between wells becomes a significant concern.Communication between wells mainly occurs ...Due to the wide application of closely spaced multi-well horizontal pads for developing unconventional gas reservoirs,interference between wells becomes a significant concern.Communication between wells mainly occurs through natural fractures.However,previous studies have found that interwell communication through natural fractures is varied,and non-communication also appears in the mid and late stages of production due to natural fracture closure.This study proposes a boundary element method for coupling multi-connected regions for the first time.Using this method,we coupled multiple flow fields to establish dual-well models with various connectivity conditions of the stimulated reservoir volume(SRV)region.These models also take into consideration of adsorption and desorption mechanism of natural gas as well as the impact of fracturing fluid retention.The study found that when considering the non-communication of SRV regions between multi-well horizontal pads,the transient behavior of the targeted well exhibits a transitional flow stage occurring before the well interference flow stage.In addition,sensitivity analysis shows that the well spacing and production regime,as well as the connectivity conditions of the SRV region,affect the timing of interwell interference.Meanwhile,the productivity of the two wells,reservoir properties,and fracturing operations affect the intensity of interwell interference.展开更多
In low-permeability geothermal reservoirs,hydro-shearing of pre-existing natural fractures plays a crucial role in improving connectivity between injection and production wells,thereby enhancing heat extraction effici...In low-permeability geothermal reservoirs,hydro-shearing of pre-existing natural fractures plays a crucial role in improving connectivity between injection and production wells,thereby enhancing heat extraction efficiency.This process increases fracture conductivity through dilation caused by injectioninduced slip;however,it also carries the risk of inducing seismic events,posing significant challenges for geothermal operations.This study employs a coupled hydro-mechanical numerical model based on the boundary element method to simulate hydro-shearing under two distinct fluid injection scenarios:(1)monotonic injection and(2)cyclic injection regulated by a traffic light system(TLS).The model assesses the effectiveness of these injection regimes in enhancing fracture conductivity while mitigating seismic hazards.Results indicate that monotonic injection frequently triggers a cascade of seismic events,disrupting pressure and stress distributions on nearby faults and resulting in complex seismic and aseismic interactions.In contrast,TLS-regulated cyclic injection,when carefully managed,promotes stable slip behavior and improves fracture conductivity.This approach proves particularly effective over extended durations during the simultaneous stimulation of two parallel faults.However,in multi-stage stimulation scenariosdwhere natural fractures are stimulated sequentiallydTLS-based cyclic injection,while more efficient at enhancing conductivity,may increase seismicity risk with prolonged application,thereby limiting its safe operational window.展开更多
Under external disturbances,the shear mechanical responses and debonding failure mechanisms at anisotropic interfaces of anchoring system composed of multiphase media are inherently difficult to characterize due to th...Under external disturbances,the shear mechanical responses and debonding failure mechanisms at anisotropic interfaces of anchoring system composed of multiphase media are inherently difficult to characterize due to the concealment nature of interfacial interactions.This study establishes an equivalent shear model for a bolt-resin-rock anchoring system and conducts direct shear tests under dynamic normal load(DNL)boundary from both laboratory experiments and discrete element method(DEM)simulations.The research investigates the influence of normal dynamic load amplitude(An)and rock type on shear strength parameters,elucidating the evolutionary characteristics and underlying mechanisms of shear load and normal displacement fluctuations induced by cyclic normal loading,with maximum shear load decreasing by 36.81%to 46.94%as An increases from 10%to 70%when rock type varies from coal to limestone.Through analysis of strain field evolution,the critical impact of rock type on localization of shear failure surface is revealed,with systematic summarization of differentiated wear characteristics,failure modes,and key controlling factors associated with shear failure surface.Mesoscopic investigations enabled by DEM simulations uncover the nonuniform distribution of contact force chains within the material matrix and across the anisotropic interfaces under various DNL boundaries,clarify rock type dependent crack propagation pathways,and quantitatively assess the damage extent of shear failure surface,with the anisotropic interface damage factor increasing from 34.9%to 56.6%as An rises from 10%to 70%,and decreasing from 49.6%to 23.4%as rock type varies from coal to limestone.展开更多
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.展开更多
The matrix expression for the 3 D transient dynamic boundary integral equation in Laplace transform space is obtained and the degenerative element method has been implemented to treat the kernel function over the sin...The matrix expression for the 3 D transient dynamic boundary integral equation in Laplace transform space is obtained and the degenerative element method has been implemented to treat the kernel function over the singular element. In the computer program BEMTDY the Koizumi′s numerical inversion method is used and three examples of the 3 D vibrated foundation under harmonic forces and the influence with both adjacent foundations are studied.展开更多
In this work, the Direct Numerical Simulation(DNS) and Oil-Film Interferometry(OFI)technique are used to investigate the hypersonic boundary layer transition induced by single and double roughness elements at Mach num...In this work, the Direct Numerical Simulation(DNS) and Oil-Film Interferometry(OFI)technique are used to investigate the hypersonic boundary layer transition induced by single and double roughness elements at Mach number 5. For single roughness, the DNS results showed that both horseshoe vortices and hairpin vortices caused by shear layer instability can affect the boundary layer instability. The generation of the near-wall unstable structure is the key point of boundary layer transition behind the roughness element. At the downstream of the roughness element, the interaction between horseshoe vortices and hairpin vortices will spread in the spanwise direction.For double roughness elements, the effect of the spacing between roughness elements on the transition is studied. It is found that the case of higher spacing between roughness elements is more effective for inducing transition than the lower one. The interaction between two adjacent roughness elements can suppress the evolution of horseshoe vortices in the downstream and trigger the instability of shear layer. Thus, the transition will be suppressed accordingly.展开更多
The interpolating moving least-squares （IMLS） method is discussed first in this paper. And the formulae of the IMLS method obtained by Lancaster are revised. Then on the basis of the boundary element-free method （B...The interpolating moving least-squares （IMLS） method is discussed first in this paper. And the formulae of the IMLS method obtained by Lancaster are revised. Then on the basis of the boundary element-free method （BEFM）, combining the boundary integral equation （BIE） method with the IMLS method, the improved boundary element-free method （IBEFM） for two-dimensional potential problems is presented, and the corresponding formulae of the IBEFM are obtained. In the BEFM, boundary conditions are applied directly, but the shape function in the MLS does not satisfy the property of the Kronecker ~ function. This is a problem of the BEFM, and must be solved theoretically. In the IMLS method, when the shape function satisfies the property of the Kronecker 5 function, then the boundary conditions, in the meshless method based on the IMLS method, can be applied directly. Then the IBEFM, based on the IMLS method, is a direct meshless boundary integral equation method in which the basic unknown quantity is the real solution of the nodal variables, and the boundary conditions can be applied directly and easily, thus it gives a greater computational precision. Some numerical examples are presented to demonstrate the method.展开更多
This study examines the design provisions of the Chinese GB 50011-2010 code for seismic design of buildings for the special boundary elements of T-shaped reinforced concrete walls and proposes an improved design metho...This study examines the design provisions of the Chinese GB 50011-2010 code for seismic design of buildings for the special boundary elements of T-shaped reinforced concrete walls and proposes an improved design method. Comparison of the design provisions of the GB 50011-2010 code and those of the American code ACI 318-14 indicates a possible deficiency in the T-shaped wall design provisions in GB 50011-2010. A case study of a typical T-shaped wall designed in accordance with GB 50011-2010 also indicates the insufficient extent of the boundary element at the non-flange end and overly conservative design of the flange end boundary element. Improved designs for special boundary elements ofT-shaped walls are developed using a displacement-based method. The proposed design formulas produce a longer boundary element at the non-flange end and a shorter boundary element at the flange end, relative to those of the GB 50011-2010 provisions. Extensive numerical analysis indicates that T-shaped walls designed using the proposed formulas develop inelastic drift of 0.01 for both cases of the flange in compression and in tension.展开更多
基金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.
文摘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 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.
基金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 Open Project of Key Laboratory of Aerospace EDLA,CASC(No.EDL19092208)。
文摘A computational framework for parachute inflation is developed based on the immersed boundary/finite element approach within the open-source IBAMR library.The fluid motion is solved by Peskin's diffuse-interface immersed boundary(IB)method,which is attractive for simulating moving-boundary flows with large deformations.The adaptive mesh refinement technique is employed to reduce the computational cost while retain the desired resolution.The dynamic response of the parachute is solved with the finite element approach.The canopy and cables of the parachute system are modeled with the hyperelastic material.A tether force is introduced to impose rigidity constraints for the parachute system.The accuracy and reliability of the present framework is validated by simulating inflation of a constrained square plate.Application of the present framework on several canonical cases further demonstrates its versatility for simulation of parachute inflation.
基金the National Natural Science Foundation of China(https://www.nsfc.gov.cn/,Project No.11972179)the Natural Science Foundation of Guangdong Province(http://gdstc.gd.gov.cn/,No.2020A1515010685)the Department of Education of Guangdong Province(http://edu.gd.gov.cn/,No.2020ZDZX2008).
文摘The boundary element method(BEM)is a popular method for solving acoustic wave propagation problems,especially those in exterior domains,owing to its ease in handling radiation conditions at infinity.However,BEM models must meet the requirement of 6–10 elements per wavelength,using the conventional constant,linear,or quadratic elements.Therefore,a large storage size of memory and long solution time are often needed in solving higher-frequency problems.In this work,we propose two new types of enriched elements based on conventional constant boundary elements to improve the computational efficiency of the 2D acoustic BEM.The first one uses a plane wave expansion,which can be used to model scattering problems.The second one uses a special plane wave expansion,which can be used tomodel radiation problems.Five examples are investigated to showthe advantages of the enriched elements.Compared with the conventional constant elements,the new enriched elements can deliver results with the same accuracy and in less computational time.This improvement in the computational efficiency is more evident at higher frequencies(with the nondimensional wave numbers exceeding 100).The paper concludes with the potential of our proposed enriched elements and plans for their further improvement.
基金supported by the National Key R&D Program of China(2020YFA0710500).
文摘The peridynamics(PD),as a promising nonlocal continuum mechanics theory,shines in solving discontinuous problems.Up to now,various numerical methods,such as the peridynamic mesh-free particlemethod(PD-MPM),peridynamic finite element method(PD-FEM),and peridynamic boundary element method(PD-BEM),have been proposed.PD-BEM,in particular,outperforms other methods by eliminating spurious boundary softening,efficiently handling infinite problems,and ensuring high computational accuracy.However,the existing PD-BEM is constructed exclusively for bond-based peridynamics(BBPD)with fixed Poisson’s ratio,limiting its applicability to crack propagation problems and scenarios involving infinite or semi-infinite problems.In this paper,we address these limitations by introducing the boundary element method(BEM)for ordinary state-based peridynamics(OSPD-BEM).Additionally,we present a crack propagationmodel embeddedwithin the framework ofOSPD-BEM to simulate crack propagations.To validate the effectiveness of OSPD-BEM,we conduct four numerical examples:deformation under uniaxial loading,crack initiation in a double-notched specimen,wedge-splitting test,and threepoint bending test.The results demonstrate the accuracy and efficiency of OSPD-BEM,highlighting its capability to successfully eliminate spurious boundary softening phenomena under varying Poisson’s ratios.Moreover,OSPDBEMsignificantly reduces computational time and exhibits greater consistencywith experimental results compared to PD-MPM.
基金sponsored by the Graduate Student Research and Innovation Fund of Xinyang Normal University under No.2024KYJJ012.
文摘In this paper,a generalized nth-order perturbation method based on the isogeometric boundary element method is proposed for the uncertainty analysis of broadband structural acoustic scattering problems.The Burton-Miller method is employed to solve the problem of non-unique solutions that may be encountered in the external acoustic field,and the nth-order discretization formulation of the boundary integral equation is derived.In addition,the computation of loop subdivision surfaces and the subdivision rules are introduced.In order to confirm the effectiveness of the algorithm,the computed results are contrasted and analyzed with the results under Monte Carlo simulations(MCs)through several numerical examples.
文摘The solution to heat transfer problems in two-dimensional heterogeneous media is attended based on the scaled boundary finite element method(SBFEM)coupled with equilibrated basis functions(EqBFs).The SBFEM reduces the model order by scaling the boundary solution onto the inner element.To this end,tri-lateral elements are emanated from a scaling center,followed by the development of a semi-analytical solution along the radial direction and a finite element solution along the circumferential/boundary direction.The discretization is thus limited to the boundaries of the model,and the semi-analytical radial solution is found through the solution of an eigenvalue problem,which restricts the methods’applicability to heterogeneous media.In this research,we first extracted the SBFEM formulation considering the heterogeneity of the media.Then,we replaced the semi-analytical radial solution with the EqBFs and removed the eigenvalue solution step from the SBFEM.The varying coefficients of the partial differential equation(PDE)resulting from the heterogeneity of the media are replaced by a finite series in the radial and circumferential directions of the element.A weighted residual approach is applied to the radial equation.The equilibrated radial solution series is used in the new formulation of the SBFEM.
基金supported by the National Natural Science Foundation of China(Grant Nos.52271278 and 52111530137)the Natural Science Foundation of Jiangsu Province(Grant No.SBK2022020579)the Newton Advanced Fellowships by the Royal Society(Grant No.NAF\R1\180304).
文摘The hydrodynamic performance of a high forward-speed ship in obliquely propagating waves is numerically examined to assess both free motions and wave field in comparison with a low forward-speed ship.This numerical model is based on the time-domain potential flow theory and higher-order boundary element method,where an analytical expression is completely expanded to determine the base-unsteady coupling flow imposed on the moving condition of the ship.The ship in the numerical model may possess different advancing speeds,i.e.stationary,low speed,and high speed.The role of the water depth,wave height,wave period,and incident wave angle is analyzed by means of the accurate numerical model.It is found that the resonant motions of the high forward-speed ship are triggered by comparison with the stationary one.More specifically,a higher forward speed generates a V-shaped wave region with a larger elevation,which induces stronger resonant motions corresponding to larger wave periods.The shoaling effect is adverse to the motion of the low-speed ship,but is beneficial to the resonant motion of the high-speed ship.When waves obliquely propagate toward the ship,the V-shaped wave region would be broken due to the coupling effect between roll and pitch motions.It is also demonstrated that the maximum heave motion occurs in beam seas for stationary cases but occurs in head waves for high speeds.However,the variation of the pitch motion with period is hardly affected by wave incident angles.
基金supported by the National Science Fund for Excellent Young Scholars(No.52222402)State Key Program of National Natural Science Foundation of China(No.U23A2022)+7 种基金State Key Program of National Natural Science Foundation of China(No.52234003)Sichuan Science and Technology Program(No.2022JDJQ0009)National Natural Science Foundation of China(No.52074235)Science and Technology Cooperation Project of the CNPC-SWPU Innovation Alliance(Nos.2020CX020202 and 2020CX030202)Shale Gas industry Development Institute of Sichuan Province111 Project(No.D18016)China Postdoctoral Science Foundation(No.2022M722637)the Science Foundation of Sichuan Province(No.2022NSFSC0186)。
文摘Due to the wide application of closely spaced multi-well horizontal pads for developing unconventional gas reservoirs,interference between wells becomes a significant concern.Communication between wells mainly occurs through natural fractures.However,previous studies have found that interwell communication through natural fractures is varied,and non-communication also appears in the mid and late stages of production due to natural fracture closure.This study proposes a boundary element method for coupling multi-connected regions for the first time.Using this method,we coupled multiple flow fields to establish dual-well models with various connectivity conditions of the stimulated reservoir volume(SRV)region.These models also take into consideration of adsorption and desorption mechanism of natural gas as well as the impact of fracturing fluid retention.The study found that when considering the non-communication of SRV regions between multi-well horizontal pads,the transient behavior of the targeted well exhibits a transitional flow stage occurring before the well interference flow stage.In addition,sensitivity analysis shows that the well spacing and production regime,as well as the connectivity conditions of the SRV region,affect the timing of interwell interference.Meanwhile,the productivity of the two wells,reservoir properties,and fracturing operations affect the intensity of interwell interference.
基金the financial support of the Helmholtz Association's Initiative and Networking Fund for the Helmholtz Young Investigator Group ARES(Contract number VHNG-1516).
文摘In low-permeability geothermal reservoirs,hydro-shearing of pre-existing natural fractures plays a crucial role in improving connectivity between injection and production wells,thereby enhancing heat extraction efficiency.This process increases fracture conductivity through dilation caused by injectioninduced slip;however,it also carries the risk of inducing seismic events,posing significant challenges for geothermal operations.This study employs a coupled hydro-mechanical numerical model based on the boundary element method to simulate hydro-shearing under two distinct fluid injection scenarios:(1)monotonic injection and(2)cyclic injection regulated by a traffic light system(TLS).The model assesses the effectiveness of these injection regimes in enhancing fracture conductivity while mitigating seismic hazards.Results indicate that monotonic injection frequently triggers a cascade of seismic events,disrupting pressure and stress distributions on nearby faults and resulting in complex seismic and aseismic interactions.In contrast,TLS-regulated cyclic injection,when carefully managed,promotes stable slip behavior and improves fracture conductivity.This approach proves particularly effective over extended durations during the simultaneous stimulation of two parallel faults.However,in multi-stage stimulation scenariosdwhere natural fractures are stimulated sequentiallydTLS-based cyclic injection,while more efficient at enhancing conductivity,may increase seismicity risk with prolonged application,thereby limiting its safe operational window.
基金support from the National Natural Science Foundation of China(Nos.51504247,52174092,51904290,and 52074259)the Natural Science Foundation of Jiangsu Province,China(No.BK20220157)+1 种基金the Fundamental Research Funds for the Central Universities,China(No.2022YCPY0202)the China University of Mining and Technology(CUMT)Open Sharing Fund for Large-scale Instruments and Equipment(No.DYGX-2025-47)is gratefully acknowledged.
文摘Under external disturbances,the shear mechanical responses and debonding failure mechanisms at anisotropic interfaces of anchoring system composed of multiphase media are inherently difficult to characterize due to the concealment nature of interfacial interactions.This study establishes an equivalent shear model for a bolt-resin-rock anchoring system and conducts direct shear tests under dynamic normal load(DNL)boundary from both laboratory experiments and discrete element method(DEM)simulations.The research investigates the influence of normal dynamic load amplitude(An)and rock type on shear strength parameters,elucidating the evolutionary characteristics and underlying mechanisms of shear load and normal displacement fluctuations induced by cyclic normal loading,with maximum shear load decreasing by 36.81%to 46.94%as An increases from 10%to 70%when rock type varies from coal to limestone.Through analysis of strain field evolution,the critical impact of rock type on localization of shear failure surface is revealed,with systematic summarization of differentiated wear characteristics,failure modes,and key controlling factors associated with shear failure surface.Mesoscopic investigations enabled by DEM simulations uncover the nonuniform distribution of contact force chains within the material matrix and across the anisotropic interfaces under various DNL boundaries,clarify rock type dependent crack propagation pathways,and quantitatively assess the damage extent of shear failure surface,with the anisotropic interface damage factor increasing from 34.9%to 56.6%as An rises from 10%to 70%,and decreasing from 49.6%to 23.4%as rock type varies from coal to limestone.
文摘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.
文摘The matrix expression for the 3 D transient dynamic boundary integral equation in Laplace transform space is obtained and the degenerative element method has been implemented to treat the kernel function over the singular element. In the computer program BEMTDY the Koizumi′s numerical inversion method is used and three examples of the 3 D vibrated foundation under harmonic forces and the influence with both adjacent foundations are studied.
基金the China Scholarship Council(CSC)the Aeronautics Science Foundation of China(No.20163252037)+2 种基金the China Postdoctoral Science Foundation(No.2017M610325)the Natural Science Foundation of Jiangsu Province of China(No.BK20170771)the Fundamental Research Funds for the Central Universities of China(No.NP2017202)for their support
文摘In this work, the Direct Numerical Simulation(DNS) and Oil-Film Interferometry(OFI)technique are used to investigate the hypersonic boundary layer transition induced by single and double roughness elements at Mach number 5. For single roughness, the DNS results showed that both horseshoe vortices and hairpin vortices caused by shear layer instability can affect the boundary layer instability. The generation of the near-wall unstable structure is the key point of boundary layer transition behind the roughness element. At the downstream of the roughness element, the interaction between horseshoe vortices and hairpin vortices will spread in the spanwise direction.For double roughness elements, the effect of the spacing between roughness elements on the transition is studied. It is found that the case of higher spacing between roughness elements is more effective for inducing transition than the lower one. The interaction between two adjacent roughness elements can suppress the evolution of horseshoe vortices in the downstream and trigger the instability of shear layer. Thus, the transition will be suppressed accordingly.
基金Project supported by the National Natural Science Foundation of China (Grant No 10871124)Innovation Program of Shanghai Municipal Education Commission (Grant No 09ZZ99)Shanghai Leading Academic Discipline Project (Grant No J50103)
文摘The interpolating moving least-squares （IMLS） method is discussed first in this paper. And the formulae of the IMLS method obtained by Lancaster are revised. Then on the basis of the boundary element-free method （BEFM）, combining the boundary integral equation （BIE） method with the IMLS method, the improved boundary element-free method （IBEFM） for two-dimensional potential problems is presented, and the corresponding formulae of the IBEFM are obtained. In the BEFM, boundary conditions are applied directly, but the shape function in the MLS does not satisfy the property of the Kronecker ~ function. This is a problem of the BEFM, and must be solved theoretically. In the IMLS method, when the shape function satisfies the property of the Kronecker 5 function, then the boundary conditions, in the meshless method based on the IMLS method, can be applied directly. Then the IBEFM, based on the IMLS method, is a direct meshless boundary integral equation method in which the basic unknown quantity is the real solution of the nodal variables, and the boundary conditions can be applied directly and easily, thus it gives a greater computational precision. Some numerical examples are presented to demonstrate the method.
基金International Science&Technology Cooperation Program of China under Grant No.2014DFA70950Tsinghua University Initiative Scientific Research Program under Grant No.2012THZ02-1National Natural Science Foundation of China under Grant No.91315301
文摘This study examines the design provisions of the Chinese GB 50011-2010 code for seismic design of buildings for the special boundary elements of T-shaped reinforced concrete walls and proposes an improved design method. Comparison of the design provisions of the GB 50011-2010 code and those of the American code ACI 318-14 indicates a possible deficiency in the T-shaped wall design provisions in GB 50011-2010. A case study of a typical T-shaped wall designed in accordance with GB 50011-2010 also indicates the insufficient extent of the boundary element at the non-flange end and overly conservative design of the flange end boundary element. Improved designs for special boundary elements ofT-shaped walls are developed using a displacement-based method. The proposed design formulas produce a longer boundary element at the non-flange end and a shorter boundary element at the flange end, relative to those of the GB 50011-2010 provisions. Extensive numerical analysis indicates that T-shaped walls designed using the proposed formulas develop inelastic drift of 0.01 for both cases of the flange in compression and in tension.
