Artificial boundary condition for Klein-Gordon equation by constructing mechanics structure

An innovative local artificial boundary condition is proposed to numerically solve the Cauchy problem of the Klein-Gordon equation in an unbounded domain.Initially,the equation is considered as the axial wave prop-agation in a bar supported on a spring foundation.The numerical model is then truncated by replacing the half-infinitely long bar with an equivalent mechanical structure.The effective frequency-dependent stiffness of the half-infinitely long bar is expressed as the sum of rational terms using Pade approximation.For each term,a corresponding substructure composed of dampers and masses is constructed.Finally,the equivalent mechan-ical structure is obtained by parallelly connecting these substructures.The proposed approach can be easily implemented within a standard finite element framework by incorporating additional mass points and damper elements.Numerical examples show that with just a few extra degrees of freedom,the proposed approach effec-tively suppresses artificial reflections at the truncation boundary and exhibits first-order convergence.

New artificial boundary condition for saturated soil foundations

Anew artificial boundary model based on multi-directional transmitting and viscous-spring artificial boundary theories is proposed to absorb stress waves in a saturated soil foundation in dynamic analysis. Since shear waves （S-waves） are the same in a saturated soil foundation and a single-phase medium foundation, a tangential visco-elastic boundary condition for a single-phase medium foundation can also be used for saturated soil foundations. Thus, the purpose of the artificial boundary proposed in this paper is primarily to absorb two types of P-waves in a saturated soil foundation. The main idea is that the stress of the P-waves in the saturated soil foundation is decomposed into two types. The first type of stress, δra＇ is absorbed by the first artificial boundary. The second type of stress, δrb, is balanced by the stress generated by the second artificial boundary. Ultimately, both types of P-waves （fast-P-waves and slow-P-waves） are absorbed by the artificial boundary model proposed in this paper. In particular, note that the fast-P-waves and slow-P-waves are absorbed at the position of the first boundary. Thus, the artificial boundary model proposed herein can simultaneously absorb P-fast waves, P-slow waves and shear waves. Finally, a numerical example is given to examine the proposed artificial boundary model, and the results show that it is very accurate.

3D viscous-spring artificial boundary in time domain

After a brief review of studies on artificial boundaries in dynamic soil-structure interaction, a three-dimensional viscous-spring artificial boundary （VSAB） in the time domain is developed in this paper. First, the 3D VSAB equations in the normal and tangential directions are derived based on the elastic wave motion theory. Secondly, a numerical simulation technique of wave motion equations along with the VSAB condition in the time domain is studied. Finally, numerical examples of some classical elastic wave motion problems are presented and the results are compared with the associated theoretical solutions, demonstrating that high precision and adequate stability can be achieved by using the proposed 3D VSAB. The proposed 3D VSAB can be conveniently incorporated in the general finite element program, which is commonly used to study dynamic soil-structure interaction problems.

BOUNDARY ELEMENT ANALYSIS OF CONTACT PROBLEMS USING ARTIFICIAL BOUNDARY NODE APPROACH

An improved version of the regular boundary element method, the artificial boundary node approach, is derived. A simple contact algorithm is designed and implemented into the direct boundary element, regular boundary element and artificial boundary node approaches. The exisiting and derived approaches are tested using some case studies. The results of the artificial boundary node approach are compared with those of the existing boundary element program, the regular element approach, ANSYS and analytical solution whenever possible. The results show the effectiveness of the artificial boundary node approach for a wider range of boundary offsets.

Artificial boundary conditions for Euler–Bernoulli beam equation

In a semi-discretized Euler-Bernoulli beam equa- tion, the non-nearest neighboring interaction and large span of temporal scales for wave propagations pose challenges to the effectiveness and stability for artificial boundary treat- ments. With the discrete equation regarded as an atomic lattice with a three-atom potential, two accurate artificial boundary conditions are first derived here. Reflection co- efficient and numerical tests illustrate the capability of the proposed methods. In particular, the time history treatment gives an exact boundary condition, yet with sensitivity to nu- merical implementations. The ALEX （almost EXact） bound- ary condition is numerically more effective.

Comparison of perfectly matched layer and multi-transmitting formula artificial boundary condition based on hybrid finite element formulation

The theory of perfectly matched layer （PML） artificial boundary condition （ABC）, which is characterized by absorption any wave motions with arbitrary frequency and arbitrarily incident angle, is introduced. The construction process of PML boundary based on elastodynamic partial differential equation （PDE） system is developed. Combining with velocity-stress hybrid finite element formulation, the applicability of PML boundary is investigated and the numerical reflection of PML boundary is estimated. The reflectivity of PML and multi-transmitting formula （MTF） boundary is then compared based on body wave and surface wave simulations. The results show that although PML boundary yields some reflection, its absorption performance is superior to MTF boundary in the numerical simulations of near-fault wave propagation, especially in comer and large angle grazing incidence situations. The PML boundary does not arise any unstable phenomenon and the stability of PML boundary is better than MTF boundary in hybrid finite element method. For a specified problem and analysis tolerance, the computational efficiency of PML boundary is only a little lower than MTF boundary.

THE APPROXIMATIONS OF THE EXACT BOUNDARY CONDITION AT AN ARTIFICIAL BOUNDARY FOR LINEARIZED INCOMPRESSIBLE VISCOUS FLOWS

We consider the linearized incompressible Navier-Stokes (Oseen) equations in a flat channel. A sequence of approximations to the exact boundary condition at an artificial boundary is derived. Then the original problem is reduced to a boundary value problem in a bounded domain, which is well-posed. A finite element approximation on the bounded domain is given, furthermore the error estimate of the finite element approximation is obtained. Numerical example shows that our artificial boundary conditions are very effective.

Mixed Finite Element Method and Higher-Order Local Artificial Boundary Conditions for Exterior 3-D Poisson Equation

The mixed finite element method is used to solve the exterior Poisson equations with higher-order local artificial boundary conditions in 3-D space. New unknowns are introduced to reduce the order of the derivatives of the unknown to two. The result is an equivalent mixed variational problem which was solved using bilinear finite elements. The primary advantage is that special finite elements are not needed on the adjacent layer of the artificial boundary for the higher-order derivatives. Error estimates are obtained for some local artificial boundary conditions with prescibed orders. A numerical example demonstrates the effectiveness of this method.

Artificial Boundary Conditions for Nonlocal Heat Equations on Unbounded Domain

This paper is concerned with numerical approximations of a nonlocal heat equation define on an infinite domain.Two classes of artificial boundary conditions(ABCs)are designed,namely,nonlocal analog Dirichlet-to-Neumann-type ABCs(global in time)and high-order Pad´e approximate ABCs(local in time).These ABCs reformulate the original problem into an initial-boundary-value(IBV)problem on a bounded domain.For the global ABCs,we adopt a fast evolution to enhance computational efficiency and reduce memory storage.High order fully discrete schemes,both second-order in time and space,are given to discretize two reduced problems.Extensive numerical experiments are carried out to show the accuracy and efficiency of the proposed methods.

APPROXIMATION, STABILITY AND FAST EVALUATION OF EXACT ARTIFICIAL BOUNDARY CONDITION FOR THE ONE-DIMENSIONAL HEAT EQUATION

In this paper we consider the numerical solution of the one-dimensional heat equation on unbounded domains. First an exact semi-discrete artificial boundary condition is derived by discretizing the time variable with the Crank-Nicolson method. The semi-discretized heat equation equipped with this boundary condition is then proved to be unconditionally stable, and its solution is shown to have second-order accuracy. In order to reduce the computational cost, we develop a new fast evaluation method for the convolution operation involved in the exact semi-discrete artificial boundary condition. A great advantage of this method is that the unconditional stability held by the semi-discretized heat equation is preserved. An error estimate is also given to show the dependence of numerical errors on the time step and the approximation accuracy of the convolution kernel. Finally, a simple numerical example is presented to validate the theoretical results.

ARTIFICIAL BOUNDARY METHOD FOR BURGERS＇ EQUATION USING NONLINEAR BOUNDARY CONDITIONS

This paper discusses the numerical solution of Burgers＇ equation on unbounded domains. Two artificial boundaries are introduced and boundary conditions are obtained on the artificial boundaries, which are in nonlinear forms. Then the original problem is reduced to an equivalent problem on a bounded domain. Finite difference method is applied to the reduced problem, and some numerical examples are given to show the effectiveness of the new approach.

THE GLOBAL ARTIFICIAL BOUNDARY CONDITIONS FOR NUMERICAL SIMULATIONS OF THE 3D FLOW AROUND A SUBMERGED BODY

We consider the numerical approximations of the three-dimensional steady potential flow around a body moving in a liquid of finite constant depth at constant speed and distance below a free surface in a channel. One vertical side is introduced as the upstream artificial boundary and two vertical sides are introduced as the downstream artificial boundaries. On the artificial boundaries, a sequence of high-order global artificial boundary conditions are given. Then the original problem is reduced to a problem defined on a finite computational domain, which is equivalent to a variational problem. After solving the variational problem by the finite element method, we obtain the numerical approximation of the original problem. The numerical examples show that the artificial boundary conditions given in this paper are very effective.

ARTIFICIAL BOUNDARY METHOD FOR THE THREE-DIMENSIONAL EXTERIOR PROBLEM OF ELASTICITY

The exact boundary condition on a spherical artificial boundary is derived for the three-dimensional exterior problem of linear elasticity in this paper. After this boundary condition is imposed on the artificial boundary, a reduced problem only defined in a bounded domain is obtained. A series of approximate problems with increasing accuracy can be derived if one truncates the series term in the variational formulation, which is equivalent to the reduced problem. An error estimate is presented to show how the error depends on the finite element discretization and the accuracy of the approximate problem. In the end, a numerical example is given to demonstrate the performance of the proposed method.

Exact Artificial Boundary Condition for the Poisson Equation in the Simulation of the 2D Schrodinger-Poisson System

We study the computation of ground states and time dependent solutions of the Schr¨odinger-Poisson system(SPS)on a bounded domain in 2D(i.e.in two space dimensions).On a disc-shaped domain,we derive exact artificial boundary conditions for the Poisson potential based on truncated Fourier series expansion inθ,and propose a second order finite difference scheme to solve the r-variable ODEs of the Fourier coefficients.The Poisson potential can be solved within O(M NlogN)arithmetic operations where M,N are the number of grid points in r-direction and the Fourier bases.Combined with the Poisson solver,a backward Euler and a semi-implicit/leap-frog method are proposed to compute the ground state and dynamics respectively.Numerical results are shown to confirm the accuracy and efficiency.Also we make it clear that backward Euler sine pseudospectral(BESP)method in[33]can not be applied to 2D SPS simulation.

Artificial Boundary Method for Calculating Ship Wave Resistance

The calculation of wave resistance for a ship moving at constant speed near a free surface is considered. This wave resistance is calculated with a linearized steady potential model. To deal with the unboundedness of the physical domain in the potential flow problem, we introduce one vertical side as an artificial upstream boundary and two vertical sides as the artificial downstream boundaries. On the artificial boundaries, a sequence of high-order global artificial boundary conditions are given. Then the potential flow problem is reduced to a problem defined on a finite computational domain, which is equivalent to a variational problem. The solution of the variational problem by the finite element method gives the numerical approximation of the potential flow around the ship, which was used to calculate the wave resistance. The numerical examples show the accuracy and efficiency of the proposed numerical scheme.

THE ELLIPSOID ARTIFICIAL BOUNDARY METHOD FOR THREE-DIMENSIONAL UNBOUNDED DOMAINS

The artificial boundary method is applied to solve three-dimensional exterior problems. Two kind of rotating ellipsoids are chosen as the artificial boundaries and the exact artificial boundary conditions are derived explicitly in terms of an infinite series. Then the well-posedness of the coupled variational problem is obtained. It is found that error estimates derived depend on the mesh size, truncation term and the location of the artificial boundary. Three numerical examples are presented to demonstrate the effectiveness and accuracy of the proposed method.

Artificial Boundary Conditions for Time-Fractional Telegraph Equation

In this paper,we study the numerical solution of the time-fractional telegraph equation on the unbounded domain.We first introduce the artificial boundariesГ±to get a finite computational domain.On the artificial boundariesГ±,we use the Laplace transform to construct the exact artificial boundary conditions(ABCs)to reduce the original problem to an initial-boundary value problem on a bounded domain.In addition,we propose a finite difference scheme based on the L_(1−2)formule for the Caputo fractional derivative in time direction and the central difference scheme for the spatial directional derivative to solve the reduced problem.In order to reduce the effect of unsmoothness of the solution at the initial moment,we use a fine mesh and low-order interpolation to discretize the solution near t=0.Finally,some numerical results show the efficiency and reliability of the ABCs and validate our theoretical results.

ARTIFICIAL BOUNDARY CONDITIONS FOR "VORTEX IN CELL" METHOD

This paper mainly designs artificial boundary conditions for 'vortex in cell' method in solving two-dimensional incompressible inviscid fluid under two conditions: one is with periodical initial value in one direction and the other with compact supported initial value. To mimic the vortex motion, Euler equation is transformed into vorticity-stream function and the technique of vortex in cell is applied incorporating with the artificial boundary conditions.

A POSTERIORI ERROR ESTIMATES OF A NON-CONFORMING FINITE ELEMENT METHOD FOR PROBLEMS WITH ARTIFICIAL BOUNDARY CONDITIONS

An a posteriori error estimator is obtained for a nonconforming finite element approximation of a linear elliptic problem, which is derived from a corresponding unbounded domain problem by applying a nonlocal approximate artificial boundary condition. Our method can be easily extended to obtain a class of a posteriori error estimators for various conforming and nonconforming finite element approximations of problems with different artificial boundary conditions. The reliability and efficiency of our a posteriori error estimator are rigorously proved and are verified by numerical examples.

Split Local Artificial Boundary Conditions for the Two-Dimensional Sine-Gordon Equation on R^(2)

In this paper the numerical solution of the two-dimensional sine-Gordon equation is studied.Split local artificial boundary conditions are obtained by the operator splitting method.Then the original problem is reduced to an initial boundary value problem on a bounded computational domain,which can be solved by the finite differencemethod.Several numerical examples are provided to demonstrate the effectiveness and accuracy of the proposed method,and some interesting propagation and collision behaviors of the solitary wave solutions are observed.