Three-dimensional acoustic propagation model for shallow waters based on an indirect boundary element method

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.

Alternating Direction Finite Volume Element Methods for Three-Dimensional Parabolic Equations

This paper presents alternating direction finite volume element methods for three-dimensional parabolic partial differential equations and gives four computational schemes, one is analogous to Douglas finite difference scheme with second-order splitting error, the other two schemes have third-order splitting error, and the last one is an extended LOD scheme. The L2 norm and H1 semi-norm error estimates are obtained for the first scheme and second one, respectively. Finally, two numerical examples are provided to illustrate the efficiency and accuracy of the methods.

CHEBYSHEV PSEUDOSPECTRAL-HYBRID FINITE ELEMENT METHOD FOR THREE-DIMENSIONAL VORTICITY EQUATION

In this paper,Chebyshev pseudospectral-finite element schemes are proposed for solving three dimensional vorticity equation.Some approximation results in nonisotropic Sobolev spaces are given.The generalized stability and the convergence are proved strictly.The numerical results show the advantages of this method.The technique in this paper is also applicable to other three-dimensional nonlinear problems in fluid dynamics.

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.

HIGH ACCURACY FINITE VOLUME ELEMENT METHOD FOR TWO-POINT BOUNDARY VALUE PROBLEM OF SECOND ORDER ORDINARY DIFFERENTIAL EQUATIONS

In this paper, a high accuracy finite volume element method is presented for two-point boundary value problem of second order ordinary differential equation, which differs from the high order generalized difference methods. It is proved that the method has optimal order error estimate O(h3) in H1 norm. Finally, two examples show that the method is effective.

Analysis of multiple interfacial cracks in three-dimensional bimaterials using hypersingular integro-differential equation method

By using the concept of finite-part integral, a set of hypersingular integro-differential equations for multiple interracial cracks in a three-dimensional infinite bimaterial subjected to arbitrary loads is derived. In the numerical analysis, unknown displacement discontinuities are approximated with the products of the fundamental density functions and power series. The fundamental functions are chosen to express a two-dimensional interface crack rigorously. As illustrative examples, the stress intensity factors for two rectangular interface cracks are calculated for various spacing, crack shape and elastic constants. It is shown that the stress intensity factors decrease with the crack spacing.

Diagonal form fast multipole boundary element method for 2D acoustic problems based on Burton-Miller boundary integral equation formulation and its applications

This paper describes formulation and implementation of the fast multipole boundary element method （FMBEM） for 2D acoustic problems. The kernel function expansion theory is summarized, and four building blocks of the FMBEM are described in details. They are moment calculation, moment to moment translation, moment to local translation, and local to local translation. A data structure for the quad-tree construction is proposed which can facilitate implementation. An analytical moment expression is derived, which is more accurate, stable, and efficient than direct numerical computation. Numerical examples are presented to demonstrate the accuracy and efficiency of the FMBEM, and radiation of a 2D vibration rail mode is simulated using the FMBEM.

THE COUPLING OF BOUNDARY ELEMENT AND FINITE ELEMENT METHODS FOR THE EXTERIOR NONSTATIONARY NAVIER-STOKES EQUATIONS

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.

FINITE ELEMENT METHOD FOR SOLVING TWO-DIMENSIONAL DIFFUSION-REACTION EQUATIONS OF BOUNDARY LAYER TYPE IN POROUS CATALYST PELLET

In this paper,finite element method(FEM)is used to solve two-dimensional diffu-sion-reaction equations of boundary layer type.This kind of equations are usually too complicatedand diffcult to be solved by applying the traditional methods used in chemical engineering becauseof the steep gradients of concentration and temperature.But,these difficulties are easy to be over-comed when the FEM is used.The integraded steps of solving this kind of problems by the FEMare presented in this paper.By applying the FEM to the two actual examples,the conclusion can bereached that the FEM has the advantages of simplicity and good accuracy.

SINGULAR INTEGRAL EQUATIONS AND BOUNDARY ELEMENT METHOD OF CRACKS IN THERMALLY STRESSED PLANAR SOLIDS

Using the method of the boundary integral equation, a set of singular integral equations of the hear transfer problems and the thermo-elastic problems of a crack embedded in a two-dimensional finite body is derived, and then,its numerical method is proposed by the numerical method of the singular integral equations combined with boundary element method. Moreover, the singular nature of temperature gradient field near the crack front is proved by the main-part analysis method of the singular integral equation, and the singular temperature gradients are exactly obtained. Finally, several typical examples calculated.

Simulating regular wave propagation over a submerged bar by boundary element method model and Boussinesq equation model

Numerical models based on the boundary element method and Boussinesq equation are used to simulate the wave transform over a submerged bar for regular waves.In the boundary-element-method model the linear element is used,and the integrals are computed by analytical formulas.The Boussinesq-equation model is the well-known FUNWAVE from the University of Delaware.We compare the numerical free surface displacements with the laboratory data on both gentle slope and steep slope,and find that both the models simulate the wave transform well.We further compute the agreement indexes between the numerical result and laboratory data,and the results support that the boundary-element-method model has a stable good performance,which is due to the fact that its governing equation has no restriction on nonlinearity and dispersion as compared with Boussinesq equation.

A Conceptual Numerical Model of the Wave Equation Using the Complex Variable Boundary Element Method

In this work, a conceptual numerical solution of the two-dimensional wave partial differential equation (PDE) is developed by coupling the Complex Variable Boundary Element Method (CVBEM) and a generalized Fourier series. The technique described in this work is suitable for modeling initial-boundary value problems governed by the wave equation on a rectangular domain with Dirichlet boundary conditions and an initial condition that is equal on the boundary to the boundary conditions. The new numerical scheme is based on the standard approach of decomposing the global initial-boundary value problem into a steady-state component and a time-dependent component. The steady-state component is governed by the Laplace PDE and is modeled with the CVBEM. The time-dependent component is governed by the wave PDE and is modeled using a generalized Fourier series. The approximate global solution is the sum of the CVBEM and generalized Fourier series approximations. The boundary conditions of the steady-state component are specified as the boundary conditions from the global BVP. The boundary conditions of the time-dependent component are specified to be identically zero. The initial condition of the time-dependent component is calculated as the difference between the global initial condition and the CVBEM approximation of the steady-state solution. Additionally, the generalized Fourier series approximation of the time-dependent component is fitted so as to approximately satisfy the derivative of the initial condition. It is shown that the strong formulation of the wave PDE is satisfied by the superposed approximate solutions of the time-dependent and steady-state components.

SOLUTION TO THE FORM OF POlSSON EQUATiONBY THE BOUNDARY ELEMENT METHODS

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.

APPLICATION OF THE DUAL RECIPROCITY BOUNDARY ELEMENT METHOD TO SOLUTION OF NONLINEAR DIFFERENTIAL EQUATION

In this paper the dual reciprocity boundary element method is employed to solve nonlinear differential equation del(2) u + u + epsilon u(3) = b. Results obtained in an example have a good agreement with those by FEM and show the applicability and simplicity of dual reciprocity method (DRM) in solving nonlinear differential equations.

Boundary Element Method with Non-overlapping Domain Decomposition for Diffusion Equation

A boundary element method based on non-overlapping domain decomposition method to solve the time-dependent diffusion equations is presented. The time-dependent fundamental solution is used in the formulation of boundary integrals and the time integration process always restarts from the initial time condition. The process of replacing the interface values, which needs a summation of boundary integrals related to the boundary values at previous time steps can be treated in parallel iterative procedure. Numerical experiments demonstrate that the implementation of the present algorithm is efficient.

BOUNDARY INTEGRAL EQUATIONS FOR THE BENDING PROBLEM OF PLATES ON TWO-PARAMETER FOUNDATION

By means of Fourier integral transformation of generalized function, the fundamental solution for the bending problem of plates on two-parameter foundation is derived in this paper, and the fundamental solution is expanded into a uniformly convergent series. On the basis of the above work, two boundary integral equations which are suitable to arbitrary shapes and arbitrary boundary conditions are established by means of the Rayleigh-Green identity. The content of the paper provides the powerful theories for the application of BEM in this problem.

A Direct Implementation of a Modified Boundary Integral Formulation for the Extended Fisher-Kolmogorov Equation

This study is concerned with the numerical approximation of the extended Fisher-Kolmogorov equation with a modified boundary integral method. A key aspect of this formulation is that it relaxes the domain-driven approach of a typical boundary element (BEM) technique. While its discretization keeps faith with the second order accurate BEM formulation, its implementation is element-based. This leads to a local solution of all integral equation and their final assembly into a slender and banded coefficient matrix which is far easier to manipulate numerically. This outcome is much better than working with BEM’s fully populated coefficient matrices resulting from a numerical encounter with the problem domain especially for nonlinear, transient, and heterogeneous problems. Faithful results of high accuracy are achieved when the results obtained herein are compared with those available in literature.

Error Analysis of A New Higher Order Boundary Element Method for A Uniform Flow Passing Cylinders

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.

A Boundary Integral Equation Approach for Boundary Problem of Laplace Equation

Using the second Green formula, the boundary problem of Laplace equation satisfied by potential function of static electric field is transformed to the problem of the boundary integral equation, and then a boundary integral equation approach is established by partitioning boundary using linear boundary element.