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.展开更多
To study wave-current actions on 3-D bodies a time-domain numerical model was established using a higher-order boundary element method(HOBEM).By assuming small flow velocities,the velocity potential could be expressed...To study wave-current actions on 3-D bodies a time-domain numerical model was established using a higher-order boundary element method(HOBEM).By assuming small flow velocities,the velocity potential could be expressed for linear and higher order components by perturbation expansion.A 4th-order Runge-Kutta method was applied for time marching.An artificial damping layer was adopted at the outer zone of the free surface mesh to dissipate scattering waves.Validation of the numerical method was carried out on run-up,wave exciting forces,and mean drift forces for wave-currents acting on a bottom-mounted vertical cylinder.The results were in close agreement with the results of a frequency-domain method and a published time-domain method.The model was then applied to compute wave-current forces and run-up on a Seastar mini tension-leg platform.展开更多
This study examined the impact of the leading-edge sweep angle on the vibration characteristics of a marine cycloidal propeller(MCP)blade during different ship maneuvering motions using a coupled three-dimensional bou...This study examined the impact of the leading-edge sweep angle on the vibration characteristics of a marine cycloidal propeller(MCP)blade during different ship maneuvering motions using a coupled three-dimensional boundary element method(BEM)and finite element method(FEM)approach.Through this approach,the study captured the interaction between hydrodynamics and structural dynamics,providing a comprehensive understanding of the response of the swept MCP blade.The following ship maneuvers were analyzed:bollard pull,crabbing,crash stop,cruising,and turning circle.During MCP operation,each blade undergoes one oscillation about its own longitudinal axis for each rotation of the horizontal propeller disc.The face and back of the propeller blade interchange during each oscillation.Consequently,the propeller blades are subjected to higher fluctuations in loading because of changes in the angle of attack and inflow velocity at each time instant.This results in complex and unstable fluid dynamics at the blade location.Variations in the sweep angle can profoundly influence the performance of the blade by altering the hydrodynamic loads and structural responses.The impact of the sweep angle is depicted through changes in the displacement,velocity,twisting angle,twisting moment,and von Mises stress of the blade.Furthermore,because of the load fluctuations on the blade,fatigue and load variations in each disc revolution must be considered during the design of cycloidal propellers.Therefore,a preliminary fatigue assessment for each maneuver was conducted.The research provides valuable information into the behavior of swept MCP blades under various loading conditions.展开更多
Metasurface design often requires solving field distributions across varying structural parameters and frequencies,where neural operators offer a promising avenue for fast prediction.However,conventional neural operat...Metasurface design often requires solving field distributions across varying structural parameters and frequencies,where neural operators offer a promising avenue for fast prediction.However,conventional neural operators have problems with degradation of the accuracy in multi-scale structural analysis.In this work,we propose a Generative Residual Enhanced Neural Operator(GRE-NO)framework that introduces a generative residual network to model the systematic bias of the main predictor.The core model retains the DeepONet architecture with both branch and trunk networks implemented using Fourier Neural Operators,combining strong generalization and efficient global representation.To handle the complexity of unbounded acoustic scattering problems,we integrate the Boundary Element Method(BEM)into data modeling and field computation,which reduces the problem dimensionality and enables training with samples at the 104 scale.Numerical experiments on some 2D and 3D acoustic metasurface problems demonstrate that the developed GRE-NO achieves excellent accuracy in results with relative errors under 1%in this study,outperforming conventional neural networks in accuracy of prediction.展开更多
This paper presents a comparative study of the solid weighting functions within the text of the modified partially saturated method(MPSM),which is an effective fluid-solid boundary condition in the lattice Boltzmann-d...This paper presents a comparative study of the solid weighting functions within the text of the modified partially saturated method(MPSM),which is an effective fluid-solid boundary condition in the lattice Boltzmann-discrete element coupling method(LBM-DEM).In its original form,the solid weighting function is τ-dependent.Previous studies have shown that the computational drag is viscosity-dependent when using the τ-dependent solid weighting function to solve fluid-particle interactions.To address this issue,two modified solid weighting functions,namely,the higher-order function and the solid-coverage function,are proposed.Nevertheless,the literature lacks a comparison of these functions,especially for viscosity dependence.In this study,the solid weighting functions are implemented and tested through two benchmark multiphase configurations,i.e.,a sphere settling between two parallel plates and the‘drafting,kissing and tumbling’of two settling spheres.The computational accuracy,viscosity dependence and convergence of the two modified functions are validated and compared against the original τ-dependent function.The LBM-DEM-MPSM formulation is then applied to study the settling behavior of a particle pack with varying solid fractions in a narrow fracture,which highlights the potential of employing the LBM-DEM-MPSM approach to a broader range of fluid-particle systems.展开更多
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.展开更多
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.展开更多
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.展开更多
In this paper, we developed the theory and algorithm of an elastic one-way boundary element method(BEM) and a corresponding hybrid elastic thin-slab propagator for earth media with sharp boundaries between strong co...In this paper, we developed the theory and algorithm of an elastic one-way boundary element method(BEM) and a corresponding hybrid elastic thin-slab propagator for earth media with sharp boundaries between strong contrast media. This approach can takes the advantage of accurate boundary condition of BEM and completely overcomes the weak contrast limitation of the perturbationtheory based one-way operator approach. The one-way BEM is a smooth boundary approximation, which avoids huge matrix operations in exact full BEM. In addition, the one-way BEM can model the primary-only transmitted and reflected waves and therefore is a valuable tool in elastic imaging and inversion. Through numerical tests for some simple models,we proved the validity and efficiency of the proposed method.展开更多
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.展开更多
In the adjoint-state method, the forward-propagated source wavefield and the backward-propagated receiver wavefield must be available simultaneously either for seismic imaging in migration or for gradient calculation ...In the adjoint-state method, the forward-propagated source wavefield and the backward-propagated receiver wavefield must be available simultaneously either for seismic imaging in migration or for gradient calculation in inversion. A feasible way to avoid the excessive storage demand is to reconstruct the source wavefield backward in time by storing the entire history of the wavefield in perfectly matched layers. In this paper, we make full use of the elementwise global property of the Laplace operator of the spectral element method (SEM) and propose an efficient source wavefield reconstruction method at the cost of storing the wavefield history only at single boundary layer nodes. Numerical experiments indicate that the accuracy of the proposed method is identical to that of the conventional method and is independent of the order of the Lagrange polynomials, the element type, and the temporal discretization method. In contrast, the memory-saving ratios of the conventional method versus our method is at least N when using either quadrilateral or hexahedron elements, respectively, where N is the order of the Lagrange polynomials used in the SEM. A higher memorysaving ratio is achieved with triangular elements versus quadrilaterals. The new method is applied to reverse time migration by considering the Marmousi model as a benchmark. Numerical results demonstrate that the method is able to provide the same result as the conventional method but with about 1/25 times lower storage demand. With the proposed wavefield reconstruction method, the storage demand is dramatically reduced;therefore, in-core memory storage is feasible even for large-scale three-dimensional adjoint inversion problems.展开更多
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.展开更多
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.展开更多
Due to the low permeability of tight reservoirs,throats play a significant role in controlling fluid flow.Although many studies have been conducted to investigate fluid flow in throats in the microscale domain,compara...Due to the low permeability of tight reservoirs,throats play a significant role in controlling fluid flow.Although many studies have been conducted to investigate fluid flow in throats in the microscale domain,comparatively fewer works have been devoted to study the effect of adsorption boundary layer(ABL)in throats based on the digital rock method.By considering an ABL,we investigate its effects on fluid flow.We build digital rock model based on computed tomography technology.Then,microscopic pore structures are extracted with watershed segmentation and pore geometries are meshed through Delaunay triangulation approach.Finally,using the meshed digital simulation model and finite element method,we investigate the effects of viscosity and thickness of ABL on microscale flow.Our results demonstrate that viscosity and thickness of ABL are major factors that significantly hinder fluid flow in throats.展开更多
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.展开更多
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.展开更多
The stochastic boundary element method is developed to analyze elasticity problems with random material and/or geometrical parameters and ran- domly perturbed boundaries. Based on the first-order Taylor series expansi...The stochastic boundary element method is developed to analyze elasticity problems with random material and/or geometrical parameters and ran- domly perturbed boundaries. Based on the first-order Taylor series expansion, the boundary integration equations concerning the mean and deviation of the displace- ments are derived, respectively. It is found that the randomness of material param- eters is equivalent to a random body force, so the mean and covariance matrices of unknown boundary displacements and tractions can be obtained. Furthermore, the mean and covariance of displacements and stresses at inner points can also be obtained. Numerical examples show that the proposed stochastic boundary element method gives satisfactory solutions, as compared with those obtained by theoretical analysis or other numerical methods.展开更多
This paper applies a machine learning technique to find a general and efficient numerical integration scheme for boundary element methods.A model based on the neural network multi-classification algorithmis constructe...This paper applies a machine learning technique to find a general and efficient numerical integration scheme for boundary element methods.A model based on the neural network multi-classification algorithmis constructed to find the minimum number of Gaussian quadrature points satisfying the given accuracy.The constructed model is trained by using a large amount of data calculated in the traditional boundary element method and the optimal network architecture is selected.The two-dimensional potential problem of a circular structure is tested and analyzed based on the determined model,and the accuracy of the model is about 90%.Finally,by incorporating the predicted Gaussian quadrature points into the boundary element analysis,we find that the numerical solution and the analytical solution are in good agreement,which verifies the robustness of the proposed method.展开更多
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.展开更多
基金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 Natural Science Foundation of China under (Grant No.107 72040,50709005 and 50921001)the Major National Science and Technology Projects of China under (Grant No.2008ZX05026-02)the Open Fund of State Key Laboratory of Ocean Engineering
摘要To study wave-current actions on 3-D bodies a time-domain numerical model was established using a higher-order boundary element method(HOBEM).By assuming small flow velocities,the velocity potential could be expressed for linear and higher order components by perturbation expansion.A 4th-order Runge-Kutta method was applied for time marching.An artificial damping layer was adopted at the outer zone of the free surface mesh to dissipate scattering waves.Validation of the numerical method was carried out on run-up,wave exciting forces,and mean drift forces for wave-currents acting on a bottom-mounted vertical cylinder.The results were in close agreement with the results of a frequency-domain method and a published time-domain method.The model was then applied to compute wave-current forces and run-up on a Seastar mini tension-leg platform.
摘要This study examined the impact of the leading-edge sweep angle on the vibration characteristics of a marine cycloidal propeller(MCP)blade during different ship maneuvering motions using a coupled three-dimensional boundary element method(BEM)and finite element method(FEM)approach.Through this approach,the study captured the interaction between hydrodynamics and structural dynamics,providing a comprehensive understanding of the response of the swept MCP blade.The following ship maneuvers were analyzed:bollard pull,crabbing,crash stop,cruising,and turning circle.During MCP operation,each blade undergoes one oscillation about its own longitudinal axis for each rotation of the horizontal propeller disc.The face and back of the propeller blade interchange during each oscillation.Consequently,the propeller blades are subjected to higher fluctuations in loading because of changes in the angle of attack and inflow velocity at each time instant.This results in complex and unstable fluid dynamics at the blade location.Variations in the sweep angle can profoundly influence the performance of the blade by altering the hydrodynamic loads and structural responses.The impact of the sweep angle is depicted through changes in the displacement,velocity,twisting angle,twisting moment,and von Mises stress of the blade.Furthermore,because of the load fluctuations on the blade,fatigue and load variations in each disc revolution must be considered during the design of cycloidal propellers.Therefore,a preliminary fatigue assessment for each maneuver was conducted.The research provides valuable information into the behavior of swept MCP blades under various loading conditions.
基金funded by National Natural Science Foundation of China(Grant number 12372198).
摘要Metasurface design often requires solving field distributions across varying structural parameters and frequencies,where neural operators offer a promising avenue for fast prediction.However,conventional neural operators have problems with degradation of the accuracy in multi-scale structural analysis.In this work,we propose a Generative Residual Enhanced Neural Operator(GRE-NO)framework that introduces a generative residual network to model the systematic bias of the main predictor.The core model retains the DeepONet architecture with both branch and trunk networks implemented using Fourier Neural Operators,combining strong generalization and efficient global representation.To handle the complexity of unbounded acoustic scattering problems,we integrate the Boundary Element Method(BEM)into data modeling and field computation,which reduces the problem dimensionality and enables training with samples at the 104 scale.Numerical experiments on some 2D and 3D acoustic metasurface problems demonstrate that the developed GRE-NO achieves excellent accuracy in results with relative errors under 1%in this study,outperforming conventional neural networks in accuracy of prediction.
基金supported by the National Natural Science Foundation of China (Grant No. 52304025)the Heilongjiang Provincial Natural Science Foundation of China (Grant No. YQ2024E010)the National Natural Science Foundation of China (Grant Nos. U23A20596 and U24B6004)
摘要This paper presents a comparative study of the solid weighting functions within the text of the modified partially saturated method(MPSM),which is an effective fluid-solid boundary condition in the lattice Boltzmann-discrete element coupling method(LBM-DEM).In its original form,the solid weighting function is τ-dependent.Previous studies have shown that the computational drag is viscosity-dependent when using the τ-dependent solid weighting function to solve fluid-particle interactions.To address this issue,two modified solid weighting functions,namely,the higher-order function and the solid-coverage function,are proposed.Nevertheless,the literature lacks a comparison of these functions,especially for viscosity dependence.In this study,the solid weighting functions are implemented and tested through two benchmark multiphase configurations,i.e.,a sphere settling between two parallel plates and the‘drafting,kissing and tumbling’of two settling spheres.The computational accuracy,viscosity dependence and convergence of the two modified functions are validated and compared against the original τ-dependent function.The LBM-DEM-MPSM formulation is then applied to study the settling behavior of a particle pack with varying solid fractions in a narrow fracture,which highlights the potential of employing the LBM-DEM-MPSM approach to a broader range of fluid-particle systems.
摘要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 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 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 National Scientific Foundation of China with Grant No. 41774067
摘要In this paper, we developed the theory and algorithm of an elastic one-way boundary element method(BEM) and a corresponding hybrid elastic thin-slab propagator for earth media with sharp boundaries between strong contrast media. This approach can takes the advantage of accurate boundary condition of BEM and completely overcomes the weak contrast limitation of the perturbationtheory based one-way operator approach. The one-way BEM is a smooth boundary approximation, which avoids huge matrix operations in exact full BEM. In addition, the one-way BEM can model the primary-only transmitted and reflected waves and therefore is a valuable tool in elastic imaging and inversion. Through numerical tests for some simple models,we proved the validity and efficiency of the proposed method.
基金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.
基金financial support for this work contributed by the National Key Research and Development Program of China (grant numbers 2016YFC0600101 and 2016YFC 0600201)the National Natural Science Foundation of China (grant numbers 41874065, 41604076, 41674102, 41674095, 41522401, 41574082, and 41774097)
摘要In the adjoint-state method, the forward-propagated source wavefield and the backward-propagated receiver wavefield must be available simultaneously either for seismic imaging in migration or for gradient calculation in inversion. A feasible way to avoid the excessive storage demand is to reconstruct the source wavefield backward in time by storing the entire history of the wavefield in perfectly matched layers. In this paper, we make full use of the elementwise global property of the Laplace operator of the spectral element method (SEM) and propose an efficient source wavefield reconstruction method at the cost of storing the wavefield history only at single boundary layer nodes. Numerical experiments indicate that the accuracy of the proposed method is identical to that of the conventional method and is independent of the order of the Lagrange polynomials, the element type, and the temporal discretization method. In contrast, the memory-saving ratios of the conventional method versus our method is at least N when using either quadrilateral or hexahedron elements, respectively, where N is the order of the Lagrange polynomials used in the SEM. A higher memorysaving ratio is achieved with triangular elements versus quadrilaterals. The new method is applied to reverse time migration by considering the Marmousi model as a benchmark. Numerical results demonstrate that the method is able to provide the same result as the conventional method but with about 1/25 times lower storage demand. With the proposed wavefield reconstruction method, the storage demand is dramatically reduced;therefore, in-core memory storage is feasible even for large-scale three-dimensional adjoint inversion problems.
基金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.
基金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.
基金National Natural Science Foundation of China(No.51674280,51774308,51704033,51722406,51950410591)Shandong Provincial Natural Science Foundation(ZR2019JQ21,JQ201808)+3 种基金the Fundamental Research Funds for the Central Universities(No.20CX02113A)National Science and Technology Major Project(2016ZX05014-000407)Program for Changjiang Scholars and Innovative Research Team in University(IRT_16R69)PetroChina Innovation Foundation(No.2018D-5007-0210)。
摘要Due to the low permeability of tight reservoirs,throats play a significant role in controlling fluid flow.Although many studies have been conducted to investigate fluid flow in throats in the microscale domain,comparatively fewer works have been devoted to study the effect of adsorption boundary layer(ABL)in throats based on the digital rock method.By considering an ABL,we investigate its effects on fluid flow.We build digital rock model based on computed tomography technology.Then,microscopic pore structures are extracted with watershed segmentation and pore geometries are meshed through Delaunay triangulation approach.Finally,using the meshed digital simulation model and finite element method,we investigate the effects of viscosity and thickness of ABL on microscale flow.Our results demonstrate that viscosity and thickness of ABL are major factors that significantly hinder fluid flow in throats.
基金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.
基金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.
摘要The stochastic boundary element method is developed to analyze elasticity problems with random material and/or geometrical parameters and ran- domly perturbed boundaries. Based on the first-order Taylor series expansion, the boundary integration equations concerning the mean and deviation of the displace- ments are derived, respectively. It is found that the randomness of material param- eters is equivalent to a random body force, so the mean and covariance matrices of unknown boundary displacements and tractions can be obtained. Furthermore, the mean and covariance of displacements and stresses at inner points can also be obtained. Numerical examples show that the proposed stochastic boundary element method gives satisfactory solutions, as compared with those obtained by theoretical analysis or other numerical methods.
基金The authors thank the financial support of National Natural Science Foundation of China(NSFC)under Grant(No.11702238).
摘要This paper applies a machine learning technique to find a general and efficient numerical integration scheme for boundary element methods.A model based on the neural network multi-classification algorithmis constructed to find the minimum number of Gaussian quadrature points satisfying the given accuracy.The constructed model is trained by using a large amount of data calculated in the traditional boundary element method and the optimal network architecture is selected.The two-dimensional potential problem of a circular structure is tested and analyzed based on the determined model,and the accuracy of the model is about 90%.Finally,by incorporating the predicted Gaussian quadrature points into the boundary element analysis,we find that the numerical solution and the analytical solution are in good agreement,which verifies the robustness of the proposed method.
基金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.