Fast High Order Algorithm for Three-Dimensional Helmholtz Equation Involving Impedance Boundary Condition with Large Wave Numbers

Acoustic fields with impedance boundary conditions have high engineering applications, such as noise control and evaluation of sound insulation materials, and can be approximated by three-dimensional Helmholtz boundary value problems. Finite difference method is widely applied to solving these problems due to its ease of use. However, when the wave number is large, the pollution effects are still a major difficulty in obtaining accurate numerical solutions. We develop a fast algorithm for solving three-dimensional Helmholtz boundary problems with large wave numbers. The boundary of computational domain is discrete based on high-order compact difference scheme. Using the properties of the tensor product and the discrete Fourier sine transform method, the original problem is solved by splitting it into independent small tridiagonal subsystems. Numerical examples with impedance boundary conditions are used to verify the feasibility and accuracy of the proposed algorithm. Results demonstrate that the algorithm has a fourth- order convergence in and -norms, and costs less CPU calculation time and random access memory.

Three-dimensional parabolic equation model for seismo-acoustic propagation: Theoretical development and preliminary numerical implementation

A three-dimensional（3D） parabolic equation（PE） model for sound propagation in a seismo-acoustic waveguide is developed in Cartesian coordinates, with x, y, and z representing the marching direction, the longitudinal direction, and the depth direction, respectively. Two sets of 3D PEs for horizontally homogenous media are derived by rewriting the 3D elastic motion equations and simultaneously choosing proper dependent variables. The numerical scheme is for now restricted to the y-independent bathymetry. Accuracy of the numerical scheme is validated, and its azimuthal limitation is analyzed. In addition, effects of horizontal refraction in a wedge-shaped waveguide and another waveguide with a polyline bottom are illustrated. Great efforts should be made in future to provide this model with the ability to handle arbitrarily irregular fluid-elastic interfaces.

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.

FOURIER-CHEBYSHEV PSEUDOSPECTRAL METHOD FOR THREE-DIMENSIONAL VORTICITY EQUATION WITH UNILATERALLY PERIODIC BOUNDARY CONDITION

A Fourier-Chebyshev pseudospectral scheme is proposed for three-dimensionalvorticily equation with unilaterally periodic boundary condition. The generalized stability and convergence are analysed. The numerical results are presented.

INTERFACIAL CRACK ANALYSIS IN THREE-DIMENSIONAL TRANSVERSELY ISOTROPIC BI-MATERIALS BY BOUNDARY INTEGRAL EQUATION METHOD

The integral-differential equations for three-dimensional planar interfacial cracks of arbitrary shape in transversely isotropic bimaterials were derived by virtue of the Somigliana identity and the fundamental solutions, in which the displacement discontinuities across the crack faces are the unknowns to be determined. The interface is parallel to both the planes of isotropy. The singular behaviors of displacement and stress near the crack border were analyzed and the stress singularity indexes were obtained by integral equation method. The stress intensity factors were expressed in terms of the displacement discontinuities. In the non-oscillatory case, the hyper-singular boundary integral-differential equations were reduced to hyper-singular boundary integral equations similar to those of homogeneously isotropic materials.

OpenMP-Based PCG Solver for Three-Dimensional Heat Equation

As one of the most important mathematics-physics equations, heat equation has been widely used in engineering area and computing science research. Large-scale heat problems are difficult to solve due to computational intractability. The parallelization of heat equation is available to improve the simulation model efficiency. In order to solve the three-dimensional heat problems more rapidly, the OpenMP was adopted to parallelize the preconditioned conjugate gradient （PCG） algorithm in this paper. A numerical experiment on the three-dimensional heat equation model was carried out on a computer with four cores. Based on the test results, it is found that the execution time of the original serial PCG program is about 1.71 to 2.81 times of the parallel PCG program executed with different number of threads. The experiment results also demonstrate the available performance of the parallel PCG algorithm based on OpenMP in terms of solution quality and computational performance.

Numerical Solutions of Three-Dimensional Coupled Burgers’ Equations by Using Some Numerical Methods

In this paper, we found the numerical solution of three-dimensional coupled Burgers’ Equations by using more efficient methods: Laplace Adomian decomposition method, Laplace transform homotopy perturbation method, variational iteration method, variational iteration decomposition method and variational iteration homotopy perturbation method. Example is examined to validate the efficiency and accuracy of these methods and they reduce the size of computation without the restrictive assumption to handle nonlinear terms and it gives the solutions rapidly.

A New Fourth Order Difference Approximation for the Solution of Three-dimensional Non-linear Biharmonic Equations Using Coupled Approach

This paper deals with a new higher order compact difference scheme, which is, O(h4) using coupled approach on the 19-point 3D stencil for the solution of three dimensional nonlinear biharmonic equations. At each internal grid point, the solution u(x,y,z) and its Laplacian Δ4u are obtained. The resulting stencil algo-rithm is presented and hence this new algorithm can be easily incorporated to solve many problems. The present discretization allows us to use the Dirichlet boundary conditions only and there is no need to discretize the derivative boundary conditions near the boundary. We also show that special treatment is required to handle the boundary conditions. Convergence analysis for a model problem is briefly discussed. The method is tested on three problems and compares very favourably with the corresponding second order approximation which we also discuss using coupled approach.

A NEW HIGH-ORDER ACCURACY EXPLICIT DIFFERENCE SCHEME FOR SOLVING THREE-DIMENSIONAL PARABOLIC EQUATIONS

In this paper, a new three-level explicit difference scheme with high-order accuracy is proposed for solving three-dimensional parabolic equations. The stability condition is r = Delta t/Delta x(2) = Delta t/Delta gamma(2) = Delta t/Delta z(2) less than or equal to 1/4, and the truncation error is O(Delta t(2) + Delta x(4)).

THE BOUNDARY INTEGRAL METHOD FOR THE HELMHOLTZ EQUATION WITH CRACKS INSIDE A BOUNDED DOMAIN

We consider a kind of scattering problem by a crack F that is buried in a bounded domain D, and we put a point source inside the domain D. This leads to a mixed boundary value problem to the Helmholtz equation in the domain D with a crack Г. Both sides of the crack F are given Dirichlet-impedance boundary conditions, and different boundary condition （Dirichlet, Neumann or Impedance boundary condition） is set on the boundary of D. Applying potential theory, the problem can be reformulated as a system of boundary integral equations. We establish the existence and uniqueness of the solution to the system by using the Fredholm theory.

Finite-Difference Solution of the Helmholtz Equation Based on Two Domain Decomposition Algorithms

In this paper, wave simulation with the finite difference method for the Helmholtz equation based on the domain decomposition method is investigated. The method solves the problem by iteratively solving subproblems defined on smaller subdomains. Two domain decomposition algorithms both for nonoverlapping and overlapping methods are described. More numerical computations including the benchmark Marmousi model show the effectiveness of the proposed algorithms. This method can be expected to be used in the full-waveform inversion in the future.

ON THE UNIQUENESS OF BOUNDARY INTEGRAL EQUATION FOR THE EXTERIOR HELMHOLTZ PROBLEM

From the point of view of energy analysis, the cause that the uniqueness of the boundary integral equation induced from the exterior Helmholtz problem does not hold is investigated in this paper. It is proved that the Sommerfeld's condition at the infinity is changed so that it is suitable not only for the radiative wave but also for the absorptive wave when we use the boundary integral equation to describe the exterior Helmholtz problem. There fore, the total energy of the system is conservative. The mathematical dealings to guarantee the uniqueness are discussed based upon this explanation.

A Fourth-Order Modified Method for the Cauchy Problem of the Modified Helmholtz Equation

This paper is concerned with the Cauchy problem for the modified Helmholtz equation in an infinite strip domain 0<x≤1,y∈R.The Cauchy data at x = 0 is given and the solution is then sought for the interval 0<x≤1.This problem is highly ill-posed and the solution(if it exists) does not depend continuously on the given data. In this paper,we propose a fourth-order modified method to solve the Cauchy problem. Convergence estimates are presented under the suitable choices of regularization parameters and the a priori assumption on the bounds of the exact solution.Numerical implementation is considered and the numerical examples show that our proposed method is effective and stable.

Numerical Solution for the Helmholtz Equation with Mixed Boundary Condition

We consider the numerical solution for the Helmholtz equation in R^2 with mixed boundary conditions.The solvability of this mixed boundary value problem is estab- lished by the boundary integral equation method.Based on the Green formula,we express the solution in terms of the boundary data.The key to the numerical real- ization of this method is the computation of weakly singular integrals.Numerical performances show the validity and feasibility of our method.The numerical schemes proposed in this paper have been applied in the realization of probe method for inverse scattering problems.

Semi-Analytical Solution of the 1D Helmholtz Equation, Obtained from Inversion of Symmetric Tridiagonal Matrix

An interesting semi-analytic solution is given for the Helmholtz equation. This solution is obtained from a rigorous discussion of the regularity and the inversion of the tridiagonal symmetric matrix. Then, applications are given, showing very good accuracy. This work provides also the analytical inverse of the skew-symmetric tridiagonal matrix.

The Best Finite-Difference Scheme for the Helmholtz Equation

The best finite-difference scheme for the Helmholtz equation is suggested. A method of solving obtained finite-difference scheme is developed. The efficiency and accuracy of method were tested on several examples.

On Vector Helmholtz Equation with a Coupling Boundary Condition

The Helmholtz equation is sometimes supplemented by conditions that include the specification of the boundary value of the divergence of the unknown.In this paper, we study the vector Helmholtz problem in domains of both C^(1,1)and Lipschitz.We es- tablish a rigorous variational analysis such as equivalence,existence and uniqueness. And we propose finite element approximations based on the uncoupled solutions.Fi- nally we present a convergence analysis and error estimates.

An Optimal Sixth-order Finite Difference Scheme for the Helmholtz Equation in One-dimension

In this paper, we present an optimal 3-point finite difference scheme for solving the 1D Helmholtz equation. We provide a convergence analysis to show that the scheme is sixth-order in accuracy. Based on minimizing the numerical dispersion, we propose a refined optimization rule for choosing the scheme’s weight parameters. Numerical results are presented to demonstrate the efficiency and accuracy of the optimal finite difference scheme.

Fundamental Solution of Dirichlet Boundary Value Problem of Axisymmetric Helmholtz Equation

Fundamental solution of Dirichlet boundary value problem of axisymmetric Helmholtz equation is constructed via modi?ed Bessel function of the second kind, which uni?ed the formulas of fundamental solution of Helmholtz equation, elliptic type Euler-Poisson-Darboux equation and Laplace equation in any dimensional space.