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 paper applied the isogeometric boundary element method(IGABEM)to thermoelastic problems.The Non-Uniform Rational B-splines(NURBS)used to construct geometric models are employed to discretize the boundary integral ...The paper applied the isogeometric boundary element method(IGABEM)to thermoelastic problems.The Non-Uniform Rational B-splines(NURBS)used to construct geometric models are employed to discretize the boundary integral formulation of the governing equation.Due to the existence of thermal stress,the domain integral term appears in the boundary integral equation.We resolve this problem by incorporating radial integration method into IGABEM which converts the domain integral to the boundary integral.In this way,IGABEM can maintain its advantages in dimensionality reduction and more importantly,seamless integration of CAD and numerical analysis based on boundary representation.The algorithm is verified by numerical examples.展开更多
This paper presents an elasto-viscoplastic consistent tangent operator (CTO) based boundary element formulation, and application for calculation of path-domain independentJ integrals (extension of the classicalJ integ...This paper presents an elasto-viscoplastic consistent tangent operator (CTO) based boundary element formulation, and application for calculation of path-domain independentJ integrals (extension of the classicalJ integrals) in nonlinear crack analysis. When viscoplastic deformation happens, the effective stresses around the crack tip in the nonlinear region is allowed to exceed the loading surface, and the pure plastic theory is not suitable for this situation. The concept of consistency employed in the solution of increment viscoplastic problem, plays a crucial role in preserving the quadratic rate asymptotic convergence of iteractive schemes based on Newton's method. Therefore, this paper investigates the viscoplastic crack problem, and presents an implicit viscoplastic algorithm using the CTO concept in a boundary element framework for path-domain independentJ integrals. Applications are presented with two numerical examples for viscoplastic crack problems andJ integrals.展开更多
This work presents some numerical aspects of isogeometric boundary element methods(IGABEM).The behavior of hyper-singular and nearly-singular integration is first explored on the distorted NURBS surface.Several numeri...This work presents some numerical aspects of isogeometric boundary element methods(IGABEM).The behavior of hyper-singular and nearly-singular integration is first explored on the distorted NURBS surface.Several numerical treatments are proposed to enhance the quadrature in the framework of isogeometric analysis.Then a numerical implementation of IGABEM on the trimmed NURBS is detailed.Based on this idea,the surface crack problem is modeled incorporation with the phantom element method.The proposed method allows the crack to intersect with the boundary of the body while preserving the original parametrization of the NURBS-based CAD geometry.展开更多
In this paper, the domain integral of the form of Poisson equation is translatedinto complete boundary integral by the fundamental solution of higher-order Laplaceoperator, the dimensions of the problem can be contrac...In this paper, the domain integral of the form of Poisson equation is translatedinto complete boundary integral by the fundamental solution of higher-order Laplaceoperator, the dimensions of the problem can be contracted into one. The numericalexamples for Stokes equations show that this method is efficient.展开更多
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.展开更多
In this work, the vortex methods for Euler equations with initial boundary value problem is considered, Poisson equations are solved using boundary element methods which can be seen to require less operations to compu...In this work, the vortex methods for Euler equations with initial boundary value problem is considered, Poisson equations are solved using boundary element methods which can be seen to require less operations to compute the velocity field from the vorticity by Chorin([6]). We prove that the rate of convergence of the boundary element schemes can be independent of the vortex blob parameters.展开更多
In this work,an acoustic topology optimizationmethod for structural surface design covered by porous materials is proposed.The analysis of acoustic problems is performed using the isogeometric boundary elementmethod.T...In this work,an acoustic topology optimizationmethod for structural surface design covered by porous materials is proposed.The analysis of acoustic problems is performed using the isogeometric boundary elementmethod.Taking the element density of porousmaterials as the design variable,the volume of porousmaterials as the constraint,and the minimum sound pressure or maximum scattered sound power as the design goal,the topology optimization is carried out by solid isotropic material with penalization(SIMP)method.To get a limpid 0–1 distribution,a smoothing Heaviside-like function is proposed.To obtain the gradient value of the objective function,a sensitivity analysis method based on the adjoint variable method(AVM)is proposed.To find the optimal solution,the optimization problems are solved by the method of moving asymptotes(MMA)based on gradient information.Numerical examples verify the effectiveness of the proposed topology optimization method in the optimization process of two-dimensional acoustic problems.Furthermore,the optimal distribution of sound-absorbingmaterials is highly frequency-dependent and usually needs to be performed within a frequency band.展开更多
A higher order boundary element method(HOBEM)is presented for inviscid flow passing cylinders in bounded or unbounded domain.The traditional boundary integral equation is established with respect to the velocity poten...A higher order boundary element method(HOBEM)is presented for inviscid flow passing cylinders in bounded or unbounded domain.The traditional boundary integral equation is established with respect to the velocity potential and its normal derivative.In present work,a new integral equation is derived for the tangential velocity.The boundary is discretized into higher order elements to ensure the continuity of slope at the element nodes.The velocity potential is also expanded with higher order shape functions,in which the unknown coefficients involve the tangential velocity.The expansion then ensures the continuities of the velocity and the slope of the boundary at element nodes.Through extensive comparison of the results for the analytical solution of cylinders,it is shown that the present HOBEM is much more accurate than the conventional BEM.展开更多
This work has a two-fold purpose.On the one hand,the theoretical formulation of a three-dimensional(3D)acoustic propagation model for shallow waters with a constant sound speed is presented,based on the boundary eleme...This work has a two-fold purpose.On the one hand,the theoretical formulation of a three-dimensional(3D)acoustic propagation model for shallow waters with a constant sound speed is presented,based on the boundary element method(BEM),which uses a half-space Green function instead of the more conventional free-space Green function.On the other hand,a numerical implementation is illustrated to explore the formulation in simple idealized cases,controlled by a few parameters,which provides necessary tests for the accuracy and performance of the model.The half-space Green's function,which has been previously used in scattering and diffraction,adds terms to the usual expressions of the integral operators without altering their continuity properties.Verifications against the wavenumber integration solution of the Pekeris waveguide suggest that the model allows an adequate prediction for the acoustic field.Likewise,numerical experiments in relation to the necessary mesh size for the description of the water-marine sediment interface lead to the conclusion that a transmission loss prediction with acceptable accuracy can be obtained with the use of a limited mesh around the desired evaluation region.展开更多
Consider a time-harmonic electromagnetic plane wave incident on a biperiodic structure in R^3. The periodic structure separates two homogeneous regions. The medium inside the structure is chiral and nonhomogeneous. In...Consider a time-harmonic electromagnetic plane wave incident on a biperiodic structure in R^3. The periodic structure separates two homogeneous regions. The medium inside the structure is chiral and nonhomogeneous. In this paper, variational formulations coupling finite element methods in the chiral medium with a method of integral equations on the periodic interfaces are studied. The well-posedness of the continuous and discretized problems is established. Uniform convergence for the coupling variational approximations of the model problem is obtained.展开更多
In this paper, we represent a new numerical method for solving the nonstationary Stokes equations in an unbounded domain. The technique consists in coupling the boundary integral and finite element methods. The variat...In this paper, we represent a new numerical method for solving the nonstationary Stokes equations in an unbounded domain. The technique consists in coupling the boundary integral and finite element methods. The variational formulation and well posedness of the coupling method are obtained. The convergence and optimal estimates for the approximation solution are provided.展开更多
We investigate time domain boundary element methods for the wave equation in R3, with a view towards sound emission problems in computational acoustics. The Neumann problem is reduced to a time dependent integral equa...We investigate time domain boundary element methods for the wave equation in R3, with a view towards sound emission problems in computational acoustics. The Neumann problem is reduced to a time dependent integral equation for the hypersingular operator, and we present a priori and a posteriori error estimates for conforming Galerkin approxima- tions in the more general case of a screen. Numerical experiments validate the convergence of our boundary element scheme and compare it with the numerical approximations ob- tained from an integral equation of the second kind. Computations in a half-space illustrate the influence of the reflection properties of a flat street.展开更多
The present article is concerned with the numerical solution of boundary integral e- quations by an adaptive wavelet boundary element method. This method approximates the solution with a computational complexity that ...The present article is concerned with the numerical solution of boundary integral e- quations by an adaptive wavelet boundary element method. This method approximates the solution with a computational complexity that is proportional to the solution＇s best N-term approximation. The focus of this article is on algorithmic issues which includes the crucial building blocks and details about the efficient implementation. By numerical examples for the Laplace equation and the Helmholtz equation, solved for different geome- tries and right-hand sides, we validate the feasibility and efficiency of the adaptive wavelet boundary element method.展开更多
A stochastic boundary element method (SBEM) is developed for 3-Dproblems with body forces. The integral equations of SBEM are established by the approach of partial derivation with respect to stochastic variables, con...A stochastic boundary element method (SBEM) is developed for 3-Dproblems with body forces. The integral equations of SBEM are established by the approach of partial derivation with respect to stochastic variables, considering the strengthlimit, rotation speeds and material density to be the fundamental stochastic variables.The method developed is applied to analyzing the strength reliability of the turbo diskof an aero-engine.展开更多
In this paper, an identification method to estimate the unbalances is introduced, which is based on the boundary element method (BEM). By using the vibration response measured at some points on the flexible rotor the ...In this paper, an identification method to estimate the unbalances is introduced, which is based on the boundary element method (BEM). By using the vibration response measured at some points on the flexible rotor the unbalances can be identified conveniently. Therefore, the rotor can be balanced without test runs.展开更多
In the present paper, we examine the performance of an efficient type of wave-absorbing porous marine structure under the attack of regular oblique waves by using a Multi-Domain Boundary Element Method(MDBEM). The str...In the present paper, we examine the performance of an efficient type of wave-absorbing porous marine structure under the attack of regular oblique waves by using a Multi-Domain Boundary Element Method(MDBEM). The structure consists of two perforated vertical thin barriers creating what can be called a wave absorbing chamber system. The barriers are surface piercing, thereby eliminating wave overtopping. The problem of the interaction of obliquely incident linear waves upon a pair of perforated barriers is first formulated in the context of linear diffraction theory. The resulting boundary integral equation, which is matched with far-field solutions presented in terms of analytical series with unknown coefficients, as well as the appropriate boundary conditions at the free surface, seabed, and barriers, is then solved numerically using MDBEM. Dissipation of the wave energy due to the presence of the perforated barriers is represented by a simple yet effective relation in terms of the porosity parameter appropriate for thin perforated walls. The results are presented in terms of reflection and transmission coefficients. The effects of the incident wave angles, relative water depths, porosities, depths of the walls, and other major parameters of interest are explored.展开更多
基金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.
基金This study was funded by the National Natural Science Foundation of China(NSFC)(Grant Nos.11702238,51904202 and 11902212)and Nanhu Scholars Program for Young Scholars of XYNU.
文摘The paper applied the isogeometric boundary element method(IGABEM)to thermoelastic problems.The Non-Uniform Rational B-splines(NURBS)used to construct geometric models are employed to discretize the boundary integral formulation of the governing equation.Due to the existence of thermal stress,the domain integral term appears in the boundary integral equation.We resolve this problem by incorporating radial integration method into IGABEM which converts the domain integral to the boundary integral.In this way,IGABEM can maintain its advantages in dimensionality reduction and more importantly,seamless integration of CAD and numerical analysis based on boundary representation.The algorithm is verified by numerical examples.
基金The project supported by National Natural Science Foundation of China(9713008)Zhejiang Natural Science Foundation Special Funds No. RC.9601
文摘This paper presents an elasto-viscoplastic consistent tangent operator (CTO) based boundary element formulation, and application for calculation of path-domain independentJ integrals (extension of the classicalJ integrals) in nonlinear crack analysis. When viscoplastic deformation happens, the effective stresses around the crack tip in the nonlinear region is allowed to exceed the loading surface, and the pure plastic theory is not suitable for this situation. The concept of consistency employed in the solution of increment viscoplastic problem, plays a crucial role in preserving the quadratic rate asymptotic convergence of iteractive schemes based on Newton's method. Therefore, this paper investigates the viscoplastic crack problem, and presents an implicit viscoplastic algorithm using the CTO concept in a boundary element framework for path-domain independentJ integrals. Applications are presented with two numerical examples for viscoplastic crack problems andJ integrals.
基金National Natural Science Foundation of China(NSFC)under Grant(No.51904202).
文摘This work presents some numerical aspects of isogeometric boundary element methods(IGABEM).The behavior of hyper-singular and nearly-singular integration is first explored on the distorted NURBS surface.Several numerical treatments are proposed to enhance the quadrature in the framework of isogeometric analysis.Then a numerical implementation of IGABEM on the trimmed NURBS is detailed.Based on this idea,the surface crack problem is modeled incorporation with the phantom element method.The proposed method allows the crack to intersect with the boundary of the body while preserving the original parametrization of the NURBS-based CAD geometry.
文摘In this paper, the domain integral of the form of Poisson equation is translatedinto complete boundary integral by the fundamental solution of higher-order Laplaceoperator, the dimensions of the problem can be contracted into one. The numericalexamples for Stokes equations show that this method is efficient.
基金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.
文摘In this work, the vortex methods for Euler equations with initial boundary value problem is considered, Poisson equations are solved using boundary element methods which can be seen to require less operations to compute the velocity field from the vorticity by Chorin([6]). We prove that the rate of convergence of the boundary element schemes can be independent of the vortex blob parameters.
基金sponsored by Natural Science Foundation of Henan under Grant No.222300420498.
文摘In this work,an acoustic topology optimizationmethod for structural surface design covered by porous materials is proposed.The analysis of acoustic problems is performed using the isogeometric boundary elementmethod.Taking the element density of porousmaterials as the design variable,the volume of porousmaterials as the constraint,and the minimum sound pressure or maximum scattered sound power as the design goal,the topology optimization is carried out by solid isotropic material with penalization(SIMP)method.To get a limpid 0–1 distribution,a smoothing Heaviside-like function is proposed.To obtain the gradient value of the objective function,a sensitivity analysis method based on the adjoint variable method(AVM)is proposed.To find the optimal solution,the optimization problems are solved by the method of moving asymptotes(MMA)based on gradient information.Numerical examples verify the effectiveness of the proposed topology optimization method in the optimization process of two-dimensional acoustic problems.Furthermore,the optimal distribution of sound-absorbingmaterials is highly frequency-dependent and usually needs to be performed within a frequency band.
基金financially supported by the National Natural Science Foundation of China (Grant Nos.52271276,52271319,and 52201364)the Natural Science Foundation of Jiangsu Province (Grant No.BK20201006)。
文摘A higher order boundary element method(HOBEM)is presented for inviscid flow passing cylinders in bounded or unbounded domain.The traditional boundary integral equation is established with respect to the velocity potential and its normal derivative.In present work,a new integral equation is derived for the tangential velocity.The boundary is discretized into higher order elements to ensure the continuity of slope at the element nodes.The velocity potential is also expanded with higher order shape functions,in which the unknown coefficients involve the tangential velocity.The expansion then ensures the continuities of the velocity and the slope of the boundary at element nodes.Through extensive comparison of the results for the analytical solution of cylinders,it is shown that the present HOBEM is much more accurate than the conventional BEM.
文摘This work has a two-fold purpose.On the one hand,the theoretical formulation of a three-dimensional(3D)acoustic propagation model for shallow waters with a constant sound speed is presented,based on the boundary element method(BEM),which uses a half-space Green function instead of the more conventional free-space Green function.On the other hand,a numerical implementation is illustrated to explore the formulation in simple idealized cases,controlled by a few parameters,which provides necessary tests for the accuracy and performance of the model.The half-space Green's function,which has been previously used in scattering and diffraction,adds terms to the usual expressions of the integral operators without altering their continuity properties.Verifications against the wavenumber integration solution of the Pekeris waveguide suggest that the model allows an adequate prediction for the acoustic field.Likewise,numerical experiments in relation to the necessary mesh size for the description of the water-marine sediment interface lead to the conclusion that a transmission loss prediction with acceptable accuracy can be obtained with the use of a limited mesh around the desired evaluation region.
基金the NSF grants DMS-0604790the NSF grants CCF-0514078+2 种基金the NSF grants EAR-0724527the ONR grant N000140210365the National Science Foundation of China grant 10428105
文摘Consider a time-harmonic electromagnetic plane wave incident on a biperiodic structure in R^3. The periodic structure separates two homogeneous regions. The medium inside the structure is chiral and nonhomogeneous. In this paper, variational formulations coupling finite element methods in the chiral medium with a method of integral equations on the periodic interfaces are studied. The well-posedness of the continuous and discretized problems is established. Uniform convergence for the coupling variational approximations of the model problem is obtained.
文摘In this paper, we represent a new numerical method for solving the nonstationary Stokes equations in an unbounded domain. The technique consists in coupling the boundary integral and finite element methods. The variational formulation and well posedness of the coupling method are obtained. The convergence and optimal estimates for the approximation solution are provided.
文摘We investigate time domain boundary element methods for the wave equation in R3, with a view towards sound emission problems in computational acoustics. The Neumann problem is reduced to a time dependent integral equation for the hypersingular operator, and we present a priori and a posteriori error estimates for conforming Galerkin approxima- tions in the more general case of a screen. Numerical experiments validate the convergence of our boundary element scheme and compare it with the numerical approximations ob- tained from an integral equation of the second kind. Computations in a half-space illustrate the influence of the reflection properties of a flat street.
文摘The present article is concerned with the numerical solution of boundary integral e- quations by an adaptive wavelet boundary element method. This method approximates the solution with a computational complexity that is proportional to the solution＇s best N-term approximation. The focus of this article is on algorithmic issues which includes the crucial building blocks and details about the efficient implementation. By numerical examples for the Laplace equation and the Helmholtz equation, solved for different geome- tries and right-hand sides, we validate the feasibility and efficiency of the adaptive wavelet boundary element method.
文摘A stochastic boundary element method (SBEM) is developed for 3-Dproblems with body forces. The integral equations of SBEM are established by the approach of partial derivation with respect to stochastic variables, considering the strengthlimit, rotation speeds and material density to be the fundamental stochastic variables.The method developed is applied to analyzing the strength reliability of the turbo diskof an aero-engine.
文摘In this paper, an identification method to estimate the unbalances is introduced, which is based on the boundary element method (BEM). By using the vibration response measured at some points on the flexible rotor the unbalances can be identified conveniently. Therefore, the rotor can be balanced without test runs.
文摘In the present paper, we examine the performance of an efficient type of wave-absorbing porous marine structure under the attack of regular oblique waves by using a Multi-Domain Boundary Element Method(MDBEM). The structure consists of two perforated vertical thin barriers creating what can be called a wave absorbing chamber system. The barriers are surface piercing, thereby eliminating wave overtopping. The problem of the interaction of obliquely incident linear waves upon a pair of perforated barriers is first formulated in the context of linear diffraction theory. The resulting boundary integral equation, which is matched with far-field solutions presented in terms of analytical series with unknown coefficients, as well as the appropriate boundary conditions at the free surface, seabed, and barriers, is then solved numerically using MDBEM. Dissipation of the wave energy due to the presence of the perforated barriers is represented by a simple yet effective relation in terms of the porosity parameter appropriate for thin perforated walls. The results are presented in terms of reflection and transmission coefficients. The effects of the incident wave angles, relative water depths, porosities, depths of the walls, and other major parameters of interest are explored.
