Numerical simulations were conducted on a 10-blade Sevik rotor ingesting wake downstream of two turbulence-generating grids.These simulations were based on implicit large-eddy simulation(ILES)and the boundary data imm...Numerical simulations were conducted on a 10-blade Sevik rotor ingesting wake downstream of two turbulence-generating grids.These simulations were based on implicit large-eddy simulation(ILES)and the boundary data immersion method(BDIM)for compressible flows,which were solved using a fully self-programmed Fortran code.Results show that the predicted thrust spectrum aligns closely with the experimental measurements.In addition,it captures the thrust dipole directivity of the noise around the rotating propeller due to random pressure pulsations on the blades,as well as the flow structures simultaneously.Furthermore,the differences in the statistical characteristics,flow structures,and low-frequency broadband thrust spectra due to different turbulence levels were investigated.This analysis indicates that the interaction between the upstream,which is characterized by a lower turbulence level and a higher turbulent length of scale,and the rotating propeller results in a lower amplitude in force spectra and a slight increase in the scale of tip vortices.展开更多
The diffuse-interface immersed boundary method(IBM)possesses excellent capabilities for simulating flows around complex geometries and moving boundaries.In this method,the flow field is solved on a fixed Cartesian mes...The diffuse-interface immersed boundary method(IBM)possesses excellent capabilities for simulating flows around complex geometries and moving boundaries.In this method,the flow field is solved on a fixed Cartesian mesh,while the solid boundary is discretized into a series of Lagrangian points immersed in the flow field.The boundary condition is implemented by introducing a force term into the momentum equation,and the interaction between the immersed boundary and the fluid domain is achieved via an interpolation process.Over the past decades,the diffuse-interface IBM has gained popularity and spawned many variants,effectively handling a wide range of flow problems from isothermal to thermal flows,from laminar to turbulent flows,and from complex geometries to fluidstructure interaction scenarios.This paper first outlines the basic principles of the diffuse-interface IBM,then highlights recent advancements achieved by the authors’research group,and finally shows the method’s excellent numerical performance and wide applicability through several case studies involving complex moving boundary problems.展开更多
This study explores a sensitivity analysis method based on the boundary element method(BEM)to address the computational complexity in acoustic analysis with ground reflection problems.The advantages of BEM in acoustic...This study explores a sensitivity analysis method based on the boundary element method(BEM)to address the computational complexity in acoustic analysis with ground reflection problems.The advantages of BEM in acoustic simulations and its high computational cost in broadband problems are examined.To improve efficiency,a Taylor series expansion is applied to decouple frequency-dependent terms in BEM.Additionally,the SecondOrder Arnoldi(SOAR)model order reduction method is integrated to reduce computational costs and enhance numerical stability.Furthermore,an isogeometric sensitivity boundary integral equation is formulated using the direct differentiation method,incorporating Cauchy principal value integrals and Hadamard finite part integrals to handle singularities.The proposed method improves the computational efficiency,and the acoustic sensitivity analysis provides theoretical support for further acoustic structure optimization.展开更多
Hydraulic fracturing plays a critical role in enhancing shale gas production in deep shale reservoirs.Conventional hydraulic fracturing simulation methods rely on prefabricated grids,which can be hindered by the chall...Hydraulic fracturing plays a critical role in enhancing shale gas production in deep shale reservoirs.Conventional hydraulic fracturing simulation methods rely on prefabricated grids,which can be hindered by the challenge of being computationally overpowered.This study proposes an efficient fracturing simulator to analyze fracture morphology during hydraulic fracturing processes in deep shale gas reservoirs.The simulator integrates the boundary element displacement discontinuity method and the finite volume method to model the fluid-solid coupling process by employing a pseudo-3D fracture model to calculate the fracture height.In particular,the Broyden iteration method was introduced to improve the computational efficiency and model robustness;it achieved a 46.6%reduction in computation time compared to the Newton-Raphson method.The influences of horizontal stress differences,natural fracture density,and natural fracture angle on the modified zone of the reservoir were simulated,and the following results were observed.(1)High stress difference reservoirs have smaller stimulated reservoir area than low stress difference reservoirs.(2)A higher natural fracture angle resulted in larger modification zones at low stress differences,while the effect of a natural fracture angle at high stress differences was not significant.(3)High-density and long natural fracture zones played a significant role in enhancing the stimulated reservoir area.These findings are critical for comprehending the impact of geological parameters on deep shale reservoirs.展开更多
The boundary knot method(BKM)is a simple boundary-type meshless method.Due to the use of non-singular general solutions rather than singular fundamental solutions,BKM does not need to consider the artificial boundary....The boundary knot method(BKM)is a simple boundary-type meshless method.Due to the use of non-singular general solutions rather than singular fundamental solutions,BKM does not need to consider the artificial boundary.Therefore,this method has the merits of purely meshless,easy to program,high solution accuracy and so on.In this paper,we investigate the effectiveness of the BKM for solving Helmholtz-type problems under various conditions through a series of novel numerical experiments.The results demonstrate that the BKM is efficient and achieves high computational accuracy for problems with smooth or continuous boundary conditions.However,when applied to discontinuous boundary problems,the method exhibits significant numerical instability,potentially leading to substantial deviations in the computed results.Finally,three potential improvement strategies are proposed to mitigate this limitation.展开更多
The traditional topology optimization method of continuum structure generally uses quadrilateral elements as the basic mesh.This approach often leads to jagged boundary issues,which are traditionally addressed through...The traditional topology optimization method of continuum structure generally uses quadrilateral elements as the basic mesh.This approach often leads to jagged boundary issues,which are traditionally addressed through post-processing,potentially altering the mechanical properties of the optimized structure.A topology optimization method of Movable Morphable Smooth Boundary(MMSB)is proposed based on the idea of mesh adaptation to solve the problem of jagged boundaries and the influence of post-processing.Based on the ICM method,the rational fraction function is introduced as the filtering function,and a topology optimization model with the minimum weight as the objective and the displacement as the constraint is established.A triangular mesh is utilized as the base mesh in this method.The mesh is re-divided in the optimization process based on the contour line,and a smooth boundary parallel to the contour line is obtained.Numerical examples demonstrate that the MMSB method effectively resolves the jagged boundary issues,leading to enhanced structural performance.展开更多
The identification of the traction acting on a portion of the surface of an anisotropic solid is very important in structural health monitoring and optimal design of structures. The traction can be determined using in...The identification of the traction acting on a portion of the surface of an anisotropic solid is very important in structural health monitoring and optimal design of structures. The traction can be determined using inverse methods in which displacement or strain measurements are taken at several points on the body. This paper presents an inverse method based on the method of fundamental solutions for the traction identification problem in two-dimensional anisotropic elasticity. The method of fundamental solutions is an efficient boundary-type meshless method widely used for analyzing various problems. Since the problem is linear, the sensitivity analysis is simply performed by solving the corresponding direct problem several times with different loads. The effects of important parameters such as the number of measurement data, the position of the measurement points, the amount of measurement error, and the type of measurement, i.e., displacement or strain, on the results are also investigated. The results obtained show that the presented inverse method is suitable for the problem of traction identification. It can be concluded from the results that the use of strain measurements in the inverse analysis leads to more accurate results than the use of displacement measurements. It is also found that measurement points closer to the boundary with unknown traction provide more reliable solutions. Additionally, it is found that increasing the number of measurement points increases the accuracy of the inverse solution. However, in cases with a large number of measurement points, further increasing the number of measurement data has little effect on the results.展开更多
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.展开更多
This study analyzed 10 years of sounding data from Yongxing Island to characterize temperature,humidity,and wind profiles in the atmospheric boundary layer(ABL).Our key findings are as follows.(1)Relative humidity inv...This study analyzed 10 years of sounding data from Yongxing Island to characterize temperature,humidity,and wind profiles in the atmospheric boundary layer(ABL).Our key findings are as follows.(1)Relative humidity inversions(RIs)were the most frequent events(99.8%),followed by surface-based temperature inversions(SBIs,63%),elevated temperature inversions(EIs,50%),and low-level jets(LLJs,41%).(2)Higher near-surface temperature leads to a decrease in EI height,thickness,and intensity.(3)When EIs occurred,RIs were consistently observed beneath the base of the EIs.(4)The intensity of the RIs correlated negatively with surface humidity,whereas the intensity of the LLJs correlated positively with surface winds except in summer.(5)The boundary layer heights(BLHs)are best estimated using the potential temperature gradient method when EIs occur,otherwise,the Richardson number method is the best method.In contrast,the relative humidity and temperature gradient methods exhibit high dispersion or uncertainty.(6)The BLHs showed seasonal variability,peaking in winter(350–1450 m).These findings advance our understanding of ABL dynamics and BLH estimation over tropical islands.展开更多
The weakly compressible smooth particle hydrodynamics(WCSPH)model is studied to address the boundary pressure instability of the SPH method,resulting in the development of the SPH method with improved dynamic boundary...The weakly compressible smooth particle hydrodynamics(WCSPH)model is studied to address the boundary pressure instability of the SPH method,resulting in the development of the SPH method with improved dynamic boundary conditions.This method employs the‘fan’search method for free surface detection,effectively identifying cavity interface particles with diameters smaller than the support domain’s radius,thereby indirectly enhancing the algorithm’s accuracy.On this basis,an improved dynamic boundary condition is proposed by updating the boundary particle pressure calculation scheme to achieve a more stable and continuous pressure field,thereby effectively preventing particles from penetrating the boundary.The SPH method with improved dynamic boundary conditions is used to simulate typical high-speed impact problems such as wedge entry and dam break.The simulation results are in good agreement with the experimental data and other numerical results.展开更多
In order to solve the problem of the variable coefficient ordinary differen-tial equation on the bounded domain,the Lagrange interpolation method is used to approximate the exact solution of the equation,and the error...In order to solve the problem of the variable coefficient ordinary differen-tial equation on the bounded domain,the Lagrange interpolation method is used to approximate the exact solution of the equation,and the error between the numerical solution and the exact solution is obtained,and then compared with the error formed by the difference method,it is concluded that the Lagrange interpolation method is more effective in solving the variable coefficient ordinary differential equation.展开更多
A hybrid approach based on the immersed boundary method(IBM)is developed for computation of flow-induced sound around moving bodies.In this method,a high-fidelity direct numerical simulation(DNS)solver is used to simu...A hybrid approach based on the immersed boundary method(IBM)is developed for computation of flow-induced sound around moving bodies.In this method,a high-fidelity direct numerical simulation(DNS)solver is used to simulate the incompressible flow field.The sound field is predicted by discretizing acoustic perturbation equations(APEs)with dispersionrelation-preserving space scheme and low-dispersion and low-dissipation Runge-Kutta time integration.A sharp-interface IBM based on ghost-cell is implemented for present two-step DNS-APE approach to deal with complex moving bodies with Cartesian grids.The present method is validated through simulations of sound generation caused by flow past a rotating cylinder,an oscillating cylinder,and tandem oscillating and stationary cylinders.The sound generated by typical kinds of complicated bio-inspired locomotions,i.e.,flapping flight by wings of varied shapes and collective undulatory swimming in tandem,are investigated using present method.The results demonstrate potential of the hybrid approach in addressing flowinduced sound generation and propagation with complex moving boundaries in a fluid medium,especially for the sound characteristics of bio-mimetic flows,which might shed lights on investigations on bio-acoustics,ethology of complex animal system,and related bio-mimetic design for quietness.展开更多
Applying numerical simulation technology to investigate fluid-solid interaction involving complex curved bound-aries is vital in aircraft design,ocean,and construction engineering.However,current methods such as Latti...Applying numerical simulation technology to investigate fluid-solid interaction involving complex curved bound-aries is vital in aircraft design,ocean,and construction engineering.However,current methods such as Lattice Boltzmann(LBM)and the immersion boundary method based on solid ratio(IMB)have limitations in identifying custom curved boundaries.Meanwhile,IBM based on velocity correction(IBM-VC)suffers from inaccuracies and numerical instability.Therefore,this study introduces a high-accuracy curve boundary recognition method(IMB-CB),which identifies boundary nodes by moving the search box,and corrects the weighting function in LBM by calculating the solid ratio of the boundary nodes,achieving accurate recognition of custom curve boundaries.In addition,curve boundary image and dot methods are utilized to verify IMB-CB.The findings revealed that IMB-CB can accurately identify the boundary,showing an error of less than 1.8%with 500 lattices.Also,the flow in the custom curve boundary and aerodynamic characteristics of the NACA0012 airfoil are calculated and compared to IBM-VC.Results showed that IMB-CB yields lower lift and drag coefficient errors than IBM-VC,with a 1.45%drag coefficient error.In addition,the characteristic curve of IMB-CB is very stable,whereas that of IBM-VC is not.For the moving boundary problem,LBM-IMB-CB with discrete element method(DEM)is capable of accurately simulating the physical phenomena of multi-moving particle flow in complex curved pipelines.This research proposes a new curve boundary recognition method,which can significantly promote the stability and accuracy of fluid-solid interaction simulations and thus has huge applications in engineering.展开更多
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.展开更多
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.展开更多
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.展开更多
How to utilize existing flow control mechanisms to make profiled end wall design more flexible,efficient,and physical is a meaningful challenge.This study presents a three-dimensional inverse method for profiled end w...How to utilize existing flow control mechanisms to make profiled end wall design more flexible,efficient,and physical is a meaningful challenge.This study presents a three-dimensional inverse method for profiled end wall design to achieve the application of flow control mechanisms.The predetermined pressure distribution on the end wall is reached by modifying the end wall geometry during flow field calculation.A motion velocity model is derived from the normal momentum equation of the moving no-slip boundary to modify the end wall geometry.A Reynolds-Averaged Navier-Stokes(RANS)solver based on the Semi-Implicit Method for Pressure Linked Equations(SIMPLE)algorithm is adopted to simulate the flow field.Based on the mechanism understanding obtained through numerical optimization results,this study adopts the inverse method to redesign an optimized end wall in a compressor cascade.The results indicate that the redesigned end wall exhibits better loss reduction,reducing the overall total pressure loss by 5.5%,whereas the optimized end wall reduces it by 3%.The inverse method allows the imposition of desired influences on the end wall flow without constructing a database,making it highly flexible,efficient,and physical.展开更多
In this study,a wavelet multi-resolution interpolation Galerkin method(WMIGM)is proposed to solve linear singularly perturbed boundary value problems.Unlike conventional wavelet schemes,the proposed algorithm can be r...In this study,a wavelet multi-resolution interpolation Galerkin method(WMIGM)is proposed to solve linear singularly perturbed boundary value problems.Unlike conventional wavelet schemes,the proposed algorithm can be readily extended to special node generation techniques,such as the Shishkin node.Such a wavelet method allows a high degree of local refinement of the nodal distribution to efficiently capture localized steep gradients.All the shape functions possess the Kronecker delta property,making the imposition of boundary conditions as easy as that in the finite element method.Four numerical examples are studied to demonstrate the validity and accuracy of the proposedwavelet method.The results showthat the use ofmodified Shishkin nodes can significantly reduce numerical oscillation near the boundary layer.Compared with many other methods,the proposed method possesses satisfactory accuracy and efficiency.The theoretical and numerical results demonstrate that the order of theε-uniform convergence of this wavelet method can reach 5.展开更多
基金Supported by the National Key R&D Program of China(2022YFB3303500).
文摘Numerical simulations were conducted on a 10-blade Sevik rotor ingesting wake downstream of two turbulence-generating grids.These simulations were based on implicit large-eddy simulation(ILES)and the boundary data immersion method(BDIM)for compressible flows,which were solved using a fully self-programmed Fortran code.Results show that the predicted thrust spectrum aligns closely with the experimental measurements.In addition,it captures the thrust dipole directivity of the noise around the rotating propeller due to random pressure pulsations on the blades,as well as the flow structures simultaneously.Furthermore,the differences in the statistical characteristics,flow structures,and low-frequency broadband thrust spectra due to different turbulence levels were investigated.This analysis indicates that the interaction between the upstream,which is characterized by a lower turbulence level and a higher turbulent length of scale,and the rotating propeller results in a lower amplitude in force spectra and a slight increase in the scale of tip vortices.
基金partially supported by the National Natural Science Foundation of China(Nos.92271103,12202191)。
文摘The diffuse-interface immersed boundary method(IBM)possesses excellent capabilities for simulating flows around complex geometries and moving boundaries.In this method,the flow field is solved on a fixed Cartesian mesh,while the solid boundary is discretized into a series of Lagrangian points immersed in the flow field.The boundary condition is implemented by introducing a force term into the momentum equation,and the interaction between the immersed boundary and the fluid domain is achieved via an interpolation process.Over the past decades,the diffuse-interface IBM has gained popularity and spawned many variants,effectively handling a wide range of flow problems from isothermal to thermal flows,from laminar to turbulent flows,and from complex geometries to fluidstructure interaction scenarios.This paper first outlines the basic principles of the diffuse-interface IBM,then highlights recent advancements achieved by the authors’research group,and finally shows the method’s excellent numerical performance and wide applicability through several case studies involving complex moving boundary problems.
基金supported by the Shanxi Scholarship Council of China(Grant No.2023-036)the Natural Science Foundation of Shanxi Province(Grant No.202303021222020).
文摘This study explores a sensitivity analysis method based on the boundary element method(BEM)to address the computational complexity in acoustic analysis with ground reflection problems.The advantages of BEM in acoustic simulations and its high computational cost in broadband problems are examined.To improve efficiency,a Taylor series expansion is applied to decouple frequency-dependent terms in BEM.Additionally,the SecondOrder Arnoldi(SOAR)model order reduction method is integrated to reduce computational costs and enhance numerical stability.Furthermore,an isogeometric sensitivity boundary integral equation is formulated using the direct differentiation method,incorporating Cauchy principal value integrals and Hadamard finite part integrals to handle singularities.The proposed method improves the computational efficiency,and the acoustic sensitivity analysis provides theoretical support for further acoustic structure optimization.
文摘Hydraulic fracturing plays a critical role in enhancing shale gas production in deep shale reservoirs.Conventional hydraulic fracturing simulation methods rely on prefabricated grids,which can be hindered by the challenge of being computationally overpowered.This study proposes an efficient fracturing simulator to analyze fracture morphology during hydraulic fracturing processes in deep shale gas reservoirs.The simulator integrates the boundary element displacement discontinuity method and the finite volume method to model the fluid-solid coupling process by employing a pseudo-3D fracture model to calculate the fracture height.In particular,the Broyden iteration method was introduced to improve the computational efficiency and model robustness;it achieved a 46.6%reduction in computation time compared to the Newton-Raphson method.The influences of horizontal stress differences,natural fracture density,and natural fracture angle on the modified zone of the reservoir were simulated,and the following results were observed.(1)High stress difference reservoirs have smaller stimulated reservoir area than low stress difference reservoirs.(2)A higher natural fracture angle resulted in larger modification zones at low stress differences,while the effect of a natural fracture angle at high stress differences was not significant.(3)High-density and long natural fracture zones played a significant role in enhancing the stimulated reservoir area.These findings are critical for comprehending the impact of geological parameters on deep shale reservoirs.
基金Supported by the Key Scientific Research Plan of Colleges and Universities in Henan Province(23B140006)。
文摘The boundary knot method(BKM)is a simple boundary-type meshless method.Due to the use of non-singular general solutions rather than singular fundamental solutions,BKM does not need to consider the artificial boundary.Therefore,this method has the merits of purely meshless,easy to program,high solution accuracy and so on.In this paper,we investigate the effectiveness of the BKM for solving Helmholtz-type problems under various conditions through a series of novel numerical experiments.The results demonstrate that the BKM is efficient and achieves high computational accuracy for problems with smooth or continuous boundary conditions.However,when applied to discontinuous boundary problems,the method exhibits significant numerical instability,potentially leading to substantial deviations in the computed results.Finally,three potential improvement strategies are proposed to mitigate this limitation.
基金supported by the National Natural Science Foundation of China(Grant 12472113).
文摘The traditional topology optimization method of continuum structure generally uses quadrilateral elements as the basic mesh.This approach often leads to jagged boundary issues,which are traditionally addressed through post-processing,potentially altering the mechanical properties of the optimized structure.A topology optimization method of Movable Morphable Smooth Boundary(MMSB)is proposed based on the idea of mesh adaptation to solve the problem of jagged boundaries and the influence of post-processing.Based on the ICM method,the rational fraction function is introduced as the filtering function,and a topology optimization model with the minimum weight as the objective and the displacement as the constraint is established.A triangular mesh is utilized as the base mesh in this method.The mesh is re-divided in the optimization process based on the contour line,and a smooth boundary parallel to the contour line is obtained.Numerical examples demonstrate that the MMSB method effectively resolves the jagged boundary issues,leading to enhanced structural performance.
基金funded by Vice Chancellor of Research at Shiraz University(grant 3GFU2M1820).
文摘The identification of the traction acting on a portion of the surface of an anisotropic solid is very important in structural health monitoring and optimal design of structures. The traction can be determined using inverse methods in which displacement or strain measurements are taken at several points on the body. This paper presents an inverse method based on the method of fundamental solutions for the traction identification problem in two-dimensional anisotropic elasticity. The method of fundamental solutions is an efficient boundary-type meshless method widely used for analyzing various problems. Since the problem is linear, the sensitivity analysis is simply performed by solving the corresponding direct problem several times with different loads. The effects of important parameters such as the number of measurement data, the position of the measurement points, the amount of measurement error, and the type of measurement, i.e., displacement or strain, on the results are also investigated. The results obtained show that the presented inverse method is suitable for the problem of traction identification. It can be concluded from the results that the use of strain measurements in the inverse analysis leads to more accurate results than the use of displacement measurements. It is also found that measurement points closer to the boundary with unknown traction provide more reliable solutions. Additionally, it is found that increasing the number of measurement points increases the accuracy of the inverse solution. However, in cases with a large number of measurement points, further increasing the number of measurement data has little effect on the results.
基金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.
基金National Key Research and Development Program of China(2023YFC3008002)National Natural Science Foundation of China(U21A6001,42075059)Key Laboratory of Guangdong Province(2020B1212060025)。
文摘This study analyzed 10 years of sounding data from Yongxing Island to characterize temperature,humidity,and wind profiles in the atmospheric boundary layer(ABL).Our key findings are as follows.(1)Relative humidity inversions(RIs)were the most frequent events(99.8%),followed by surface-based temperature inversions(SBIs,63%),elevated temperature inversions(EIs,50%),and low-level jets(LLJs,41%).(2)Higher near-surface temperature leads to a decrease in EI height,thickness,and intensity.(3)When EIs occurred,RIs were consistently observed beneath the base of the EIs.(4)The intensity of the RIs correlated negatively with surface humidity,whereas the intensity of the LLJs correlated positively with surface winds except in summer.(5)The boundary layer heights(BLHs)are best estimated using the potential temperature gradient method when EIs occur,otherwise,the Richardson number method is the best method.In contrast,the relative humidity and temperature gradient methods exhibit high dispersion or uncertainty.(6)The BLHs showed seasonal variability,peaking in winter(350–1450 m).These findings advance our understanding of ABL dynamics and BLH estimation over tropical islands.
基金supported by the National Natural Science Foundation of China(Grant No.52071094).
文摘The weakly compressible smooth particle hydrodynamics(WCSPH)model is studied to address the boundary pressure instability of the SPH method,resulting in the development of the SPH method with improved dynamic boundary conditions.This method employs the‘fan’search method for free surface detection,effectively identifying cavity interface particles with diameters smaller than the support domain’s radius,thereby indirectly enhancing the algorithm’s accuracy.On this basis,an improved dynamic boundary condition is proposed by updating the boundary particle pressure calculation scheme to achieve a more stable and continuous pressure field,thereby effectively preventing particles from penetrating the boundary.The SPH method with improved dynamic boundary conditions is used to simulate typical high-speed impact problems such as wedge entry and dam break.The simulation results are in good agreement with the experimental data and other numerical results.
文摘In order to solve the problem of the variable coefficient ordinary differen-tial equation on the bounded domain,the Lagrange interpolation method is used to approximate the exact solution of the equation,and the error between the numerical solution and the exact solution is obtained,and then compared with the error formed by the difference method,it is concluded that the Lagrange interpolation method is more effective in solving the variable coefficient ordinary differential equation.
基金supported by the National Natural Science Foundation of China(Grant Nos.52301380 and 11602277).
文摘A hybrid approach based on the immersed boundary method(IBM)is developed for computation of flow-induced sound around moving bodies.In this method,a high-fidelity direct numerical simulation(DNS)solver is used to simulate the incompressible flow field.The sound field is predicted by discretizing acoustic perturbation equations(APEs)with dispersionrelation-preserving space scheme and low-dispersion and low-dissipation Runge-Kutta time integration.A sharp-interface IBM based on ghost-cell is implemented for present two-step DNS-APE approach to deal with complex moving bodies with Cartesian grids.The present method is validated through simulations of sound generation caused by flow past a rotating cylinder,an oscillating cylinder,and tandem oscillating and stationary cylinders.The sound generated by typical kinds of complicated bio-inspired locomotions,i.e.,flapping flight by wings of varied shapes and collective undulatory swimming in tandem,are investigated using present method.The results demonstrate potential of the hybrid approach in addressing flowinduced sound generation and propagation with complex moving boundaries in a fluid medium,especially for the sound characteristics of bio-mimetic flows,which might shed lights on investigations on bio-acoustics,ethology of complex animal system,and related bio-mimetic design for quietness.
基金WJD,JYZ,CLC,ZX,and ZGY were supported by the National Natural Science Foundation of China(Grant Number 51705143)the Education Department of Hunan Province(Grant Number 22B0464)the Postgraduate Scientific Research Innovation Project of Hunan Province(Grant Number QL20230249).
文摘Applying numerical simulation technology to investigate fluid-solid interaction involving complex curved bound-aries is vital in aircraft design,ocean,and construction engineering.However,current methods such as Lattice Boltzmann(LBM)and the immersion boundary method based on solid ratio(IMB)have limitations in identifying custom curved boundaries.Meanwhile,IBM based on velocity correction(IBM-VC)suffers from inaccuracies and numerical instability.Therefore,this study introduces a high-accuracy curve boundary recognition method(IMB-CB),which identifies boundary nodes by moving the search box,and corrects the weighting function in LBM by calculating the solid ratio of the boundary nodes,achieving accurate recognition of custom curve boundaries.In addition,curve boundary image and dot methods are utilized to verify IMB-CB.The findings revealed that IMB-CB can accurately identify the boundary,showing an error of less than 1.8%with 500 lattices.Also,the flow in the custom curve boundary and aerodynamic characteristics of the NACA0012 airfoil are calculated and compared to IBM-VC.Results showed that IMB-CB yields lower lift and drag coefficient errors than IBM-VC,with a 1.45%drag coefficient error.In addition,the characteristic curve of IMB-CB is very stable,whereas that of IBM-VC is not.For the moving boundary problem,LBM-IMB-CB with discrete element method(DEM)is capable of accurately simulating the physical phenomena of multi-moving particle flow in complex curved pipelines.This research proposes a new curve boundary recognition method,which can significantly promote the stability and accuracy of fluid-solid interaction simulations and thus has huge applications in engineering.
基金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 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 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.
基金Supported by the National Natural Science Foundation of China(No.52376021).
文摘How to utilize existing flow control mechanisms to make profiled end wall design more flexible,efficient,and physical is a meaningful challenge.This study presents a three-dimensional inverse method for profiled end wall design to achieve the application of flow control mechanisms.The predetermined pressure distribution on the end wall is reached by modifying the end wall geometry during flow field calculation.A motion velocity model is derived from the normal momentum equation of the moving no-slip boundary to modify the end wall geometry.A Reynolds-Averaged Navier-Stokes(RANS)solver based on the Semi-Implicit Method for Pressure Linked Equations(SIMPLE)algorithm is adopted to simulate the flow field.Based on the mechanism understanding obtained through numerical optimization results,this study adopts the inverse method to redesign an optimized end wall in a compressor cascade.The results indicate that the redesigned end wall exhibits better loss reduction,reducing the overall total pressure loss by 5.5%,whereas the optimized end wall reduces it by 3%.The inverse method allows the imposition of desired influences on the end wall flow without constructing a database,making it highly flexible,efficient,and physical.
基金supported by the National Natural Science Foundation of China (No.12172154)the 111 Project (No.B14044)+1 种基金the Natural Science Foundation of Gansu Province (No.23JRRA1035)the Natural Science Foundation of Anhui University of Finance and Economics (No.ACKYC20043).
文摘In this study,a wavelet multi-resolution interpolation Galerkin method(WMIGM)is proposed to solve linear singularly perturbed boundary value problems.Unlike conventional wavelet schemes,the proposed algorithm can be readily extended to special node generation techniques,such as the Shishkin node.Such a wavelet method allows a high degree of local refinement of the nodal distribution to efficiently capture localized steep gradients.All the shape functions possess the Kronecker delta property,making the imposition of boundary conditions as easy as that in the finite element method.Four numerical examples are studied to demonstrate the validity and accuracy of the proposedwavelet method.The results showthat the use ofmodified Shishkin nodes can significantly reduce numerical oscillation near the boundary layer.Compared with many other methods,the proposed method possesses satisfactory accuracy and efficiency.The theoretical and numerical results demonstrate that the order of theε-uniform convergence of this wavelet method can reach 5.
