Time-Domain Higher-Order Boundary Element Method for Simulating High Forward-Speed Ship Motions in Waves

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 Boundary Element Method for Ordinary State-Based Peridynamics

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.

Generalized nth-Order Perturbation Method Based on Loop Subdivision Surface Boundary Element Method for Three-Dimensional Broadband Structural Acoustic Uncertainty Analysis

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.

Enriched Constant Elements in the Boundary Element Method for Solving 2D Acoustic Problems at Higher Frequencies

The boundary element method(BEM)is a popular method for solving acoustic wave propagation problems,especially those in exterior domains,owing to its ease in handling radiation conditions at infinity.However,BEM models must meet the requirement of 6–10 elements per wavelength,using the conventional constant,linear,or quadratic elements.Therefore,a large storage size of memory and long solution time are often needed in solving higher-frequency problems.In this work,we propose two new types of enriched elements based on conventional constant boundary elements to improve the computational efficiency of the 2D acoustic BEM.The first one uses a plane wave expansion,which can be used to model scattering problems.The second one uses a special plane wave expansion,which can be used tomodel radiation problems.Five examples are investigated to showthe advantages of the enriched elements.Compared with the conventional constant elements,the new enriched elements can deliver results with the same accuracy and in less computational time.This improvement in the computational efficiency is more evident at higher frequencies(with the nondimensional wave numbers exceeding 100).The paper concludes with the potential of our proposed enriched elements and plans for their further improvement.

FINITE ELEMENT-ARTIFICIAL TRANSMITTING BOUNDARY METHOD FOR WAVE SCATTERING FROM IRREGULAR CYLINDER

The finite element artificial transmitting boundary method is employed here to treat the near field scattering of a cylindrical wave from an irregular cylinder. A comparison is made between this method and the analytical one. And then examples are given to demonstrate the solution of several problems of the irregular object scattering. The method can not only produce clear physical pictures, but can efficiently handle many complicated scattering problems.

Identification and Balancing of Flexible Rotors by Boundary Element Method

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.

STOCHASTIC BOUNDARY ELEMENT METHODS FOR 3D PROBLEMS WITH BODY FORCES AND ITS APPLICATION IN RELIABILITY OF TURBINE DISKS

The stochastic boundary element method(SBEM)is developed in this paper for 3D problems with body forces and reliability analysis of engineering structures.The integral equations of SBEM are established by the approach of partial derivation with respect to stochastic variables,considering the yield limit,rotation speeds and material density to be the fundamental stochastic variables.Through analyzing a numerical example and a turbo-disk of an aeroengine,the results show that the method developed is successful.

A Multi-Domain Boundary Element Method to Analyze the Reflection and Transmission of Oblique Waves From Double Porous Thin Walls

In the present paper, we examine the performance of an efficient type of wave-absorbing porous marine structure under the attack of regular oblique waves by using a Multi-Domain Boundary Element Method(MDBEM). The structure consists of two perforated vertical thin barriers creating what can be called a wave absorbing chamber system. The barriers are surface piercing, thereby eliminating wave overtopping. The problem of the interaction of obliquely incident linear waves upon a pair of perforated barriers is first formulated in the context of linear diffraction theory. The resulting boundary integral equation, which is matched with far-field solutions presented in terms of analytical series with unknown coefficients, as well as the appropriate boundary conditions at the free surface, seabed, and barriers, is then solved numerically using MDBEM. Dissipation of the wave energy due to the presence of the perforated barriers is represented by a simple yet effective relation in terms of the porosity parameter appropriate for thin perforated walls. The results are presented in terms of reflection and transmission coefficients. The effects of the incident wave angles, relative water depths, porosities, depths of the walls, and other major parameters of interest are explored.

A stable implicit nodal integration-based particle finite element method(N-PFEM)for modelling saturated soil dynamics

In this study,we present a novel nodal integration-based particle finite element method(N-PFEM)designed for the dynamic analysis of saturated soils.Our approach incorporates the nodal integration technique into a generalised Hellinger-Reissner(HR)variational principle,creating an implicit PFEM formulation.To mitigate the volumetric locking issue in low-order elements,we employ a node-based strain smoothing technique.By discretising field variables at the centre of smoothing cells,we achieve nodal integration over cells,eliminating the need for sophisticated mapping operations after re-meshing in the PFEM.We express the discretised governing equations as a min-max optimisation problem,which is further reformulated as a standard second-order cone programming(SOCP)problem.Stresses,pore water pressure,and displacements are simultaneously determined using the advanced primal-dual interior point method.Consequently,our numerical model offers improved accuracy for stresses and pore water pressure compared to the displacement-based PFEM formulation.Numerical experiments demonstrate that the N-PFEM efficiently captures both transient and long-term hydro-mechanical behaviour of saturated soils with high accuracy,obviating the need for stabilisation or regularisation techniques commonly employed in other nodal integration-based PFEM approaches.This work holds significant implications for the development of robust and accurate numerical tools for studying saturated soil dynamics.

Numerical Simulation for Water Entry of a Wedge at Varying Speed by a High Order Boundary Element Method

A high order boundary element method was developed for the complex velocity potential problem. The method ensures not only the continuity of the potential at the nodes of each element but also the velocity. It can be applied to a variety of velocity potential problems. The present paper, however, focused on its application to the problem of water entry of a wedge with varying speed. The continuity of the velocity achieved herein is particularly important for this kind of nonlinear free surface flow problem, because when the time stepping method is used, the free surface is updated through the velocity obtained at each node and the accuracy of the velocity is therefore crucial. Calculation was made for a case when the distance S that the wedge has travelled and time t follow the relationship s=Dtα, where D and α are constants, which is found to lead to a self similar flow field when the effect due to gravity is ignored.

Three-Dimensional Boundary Element Method Applied to Nonlinear Wave Transformation

For higher accuracy in simulating the transformation of three dimensional waves, in consideration of the advantages of constant panels and linear elements, a combined boundary elements is applied in this research. The method can be used to remove the transverse vibration due to the accumulation of computational errors. A combined boundary condition of sponge layer and Sommerfeld radiation condition is used to remove the reflected waves from the computing domain. By following the water particle on the water surface, the third order Stokes wave transform is simulated by the numerical wave flume technique. The computed results are in good agreement with theoretical ones.

Introducing an effective coherence function to generate non-uniform ground motion on topographic site using time-domain boundary element method

In this study,a comprehensive parametric analysis was performed on non-uniform excitation of V-shaped topography using the boundary element method in time domain.For this purpose,wave scattering analysis was carried out on a topography subjected to the SV-wave for different predominant frequencies and shape ratios.Based on the numerical results,new coherence and time delay functions are proposed to generate non-uniform ground motion for topographic irregularities.The efficiency and accuracy of the proposed functions for real engineering problems are indicated by comparison with observations reported in previous literature.

BOUNDARY ELEMENT METHOD FOR MOVING AND ROLLING CONTACT OF 2D ELASTIC BODIES WITH DEFECTS

A scheme of boundary element method for moving contact of two-dimensional elastic bodies using conforming discretization is presented. Both the displacement and the traction boundary conditions are satisfied on the contacting region in the sense of discretization. An algorithm to deal with the moving of the contact boundary on a larger possible contact region is presented. The algorithm is generalized to rolling contact problem as well. Some numerical examples of moving and rolling contact of 2D elastic bodies with or without friction, including the bodies with a hole-type defect, are given to show the effectiveness and the accuracy of the presented schemes.

Isogeometric Boundary Element Analysis for 2D Transient Heat Conduction Problem with Radial Integration Method

This paper presents an isogeometric boundary element method(IGABEM)for transient heat conduction analysis.The Non-Uniform Rational B-spline(NURBS)basis functions,which are used to construct the geometry of the structures,are employed to discretize the physical unknowns in the boundary integral formulations of the governing equations.Bezier extraction technique is employed to accelerate the evaluation of NURBS basis functions.We adopt a radial integration method to address the additional domain integrals.The numerical examples demonstrate the advantage of IGABEM in dimension reduction and the seamless connection between CAD and numerical analysis.

A new magneto-cardiogram study using a vector model with a virtual heart and the boundary element method

A cardiac vector model is presented and verified, and then the forward problem for cardiac magnetic fields and electric potential are discussed based on this model and the realistic human torso volume conductor model, including lungs. A torso-cardiac vector model is used for a 12-lead electrocardiographic （ECG） and magneto-cardiogram （MCG） simulation study by using the boundary element method （BEM）. Also, we obtain the MCG wave picture using a compound four-channel HTc.SQUID system in a magnetically shielded room. By comparing the simulated results and experimental results, we verify the cardiac vector model and then do a preliminary study of the forward problem of MCG and ECG. Therefore, the results show that the vector model is reasonable in cardiac electrophysiology.

Forward and inverse problem for cardiac magnetic field and electric potential using two boundary element methods

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.

Monte Carlo Simulation of Fractures Using Isogeometric Boundary Element Methods Based on POD-RBF

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.

A Combined Shape and Topology Optimization Based on Isogeometric Boundary Element Method for 3D Acoustics

A combined shape and topology optimization algorithm based on isogeometric boundary element method for 3D acoustics is developed in this study.The key treatment involves using adjoint variable method in shape sensitivity analysis with respect to non-uniform rational basis splines control points,and in topology sensitivity analysis with respect to the artificial densities of sound absorption material.OpenMP tool in Fortran code is adopted to improve the efficiency of analysis.To consider the features and efficiencies of the two types of optimization methods,this study adopts a combined iteration scheme for the optimization process to investigate the simultaneous change of geometry shape and distribution of material to achieve better noise control.Numerical examples,such as sound barrier,simple tank,and BeTSSi submarine,are performed to validate the advantage of combined optimization in noise reduction,and to demonstrate the potential of the proposed method for engineering problems.

An improved boundary element-free method (IBEFM) for two-dimensional potential problems

The interpolating moving least-squares （IMLS） method is discussed first in this paper. And the formulae of the IMLS method obtained by Lancaster are revised. Then on the basis of the boundary element-free method （BEFM）, combining the boundary integral equation （BIE） method with the IMLS method, the improved boundary element-free method （IBEFM） for two-dimensional potential problems is presented, and the corresponding formulae of the IBEFM are obtained. In the BEFM, boundary conditions are applied directly, but the shape function in the MLS does not satisfy the property of the Kronecker ~ function. This is a problem of the BEFM, and must be solved theoretically. In the IMLS method, when the shape function satisfies the property of the Kronecker 5 function, then the boundary conditions, in the meshless method based on the IMLS method, can be applied directly. Then the IBEFM, based on the IMLS method, is a direct meshless boundary integral equation method in which the basic unknown quantity is the real solution of the nodal variables, and the boundary conditions can be applied directly and easily, thus it gives a greater computational precision. Some numerical examples are presented to demonstrate the method.

Application of dual reciprocity boundary element method to predict acoustic attenuation characteristics of marine engine exhaust silencers

In marine engine exhaust silencing systems, the presence of exhaust gas flow influences the sound propagation inside the systems and the acoustic attenuation performance of silencers. In order to investigate the effects of three-dimensional gas flow and acoustic damping on the acoustic attenuation characteristics of marine engine exhaust silencers, a dual reciprocity boundary element method (DRBEM) was developed. The acoustic governing equation in three-dimensional potential flow was derived first, and then the DRBEM numerical procedure is given. Compared to the conventional boundary element method (CBEM), the DRBEM considers the second order terms of flow Mach number in the acoustic governing equation, so it is suitable for the cases with higher Mach number subsonic flow. For complex exhaust silencers, it is difficult to apply the single-domain boundary element method, so a substructure approach based on the dual reciprocity boundary element method is presented. The experiments for measuring transmission loss of silencers are conducted, and the experimental setup and measurements are explained. The transmission loss of a single expansion chamber silencer with extended inlet and outlet were predicted by DRBEM and compared with the measurements. The good agreements between predictions and measurements are observed, which demonstrated that the derived acoustic governing equation and the DRBEM numerical procedure in the present study are correct.