Absorbing Boundary Conditions for Simulating Water Waves Near Solid Bodies

The objective of this paper is to present a new method for designing absorbing or non-reflective boundary conditions (ABC) or (NRBC), illustrated by the case study of the modelling of a solid body in water, specifically the capillary gravity waves generated by its motion at the surface. The study analyses the flow of an inviscid, barotropic, and compressible fluid around the stationary solid body. The dynamic behaviour of the fluid is analysed using a two-dimensional coupled Neumann-Kelvin model extended with capillarity and inertia terms. For computational purposes, it is necessary to truncate the unbounded spatial domain with artificial boundaries and then introduce appropriate absorbing boundary conditions. The propagation of short wavelength waves in a convective fluid medium with significant differences in properties between the interior and the surface of the fluid presents a number of difficulties in the design of these conditions. The results are illustrated numerically and commented upon.

Application of the double absorbing boundary condition in seismic modeling

We apply the newly proposed double absorbing boundary condition（DABC）（Hagstrom et al., 2014） to solve the boundary reflection problem in seismic finite-difference（FD） modeling. In the DABC scheme, the local high-order absorbing boundary condition is used on two parallel artificial boundaries, and thus double absorption is achieved. Using the general 2D acoustic wave propagation equations as an example, we use the DABC in seismic FD modeling, and discuss the derivation and implementation steps in detail. Compared with the perfectly matched layer（PML）, the complexity decreases, and the stability and fl exibility improve. A homogeneous model and the SEG salt model are selected for numerical experiments. The results show that absorption using the DABC is considerably improved relative to the Clayton–Engquist boundary condition and nearly the same as that in the PML.

Perfectly matched layer-absorbing boundary condition for finite-element time-domain modeling of elastic wave equations

The perfectly matched layer （PML） is a highly efficient absorbing boundary condition used for the numerical modeling of seismic wave equation. The article focuses on the application of this technique to finite-element time-domain numerical modeling of elastic wave equation. However, the finite-element time-domain scheme is based on the second- order wave equation in displacement formulation. Thus, the first-order PML in velocity-stress formulation cannot be directly applied to this scheme. In this article, we derive the finite- element matrix equations of second-order PML in displacement formulation, and accomplish the implementation of PML in finite-element time-domain modeling of elastic wave equation. The PML has an approximate zero reflection coefficients for bulk and surface waves in the finite-element modeling of P-SV and SH wave propagation in the 2D homogeneous elastic media. The numerical experiments using a two-layer model with irregular topography validate the efficiency of PML in the modeling of seismic wave propagation in geological models with complex structures and heterogeneous media.

Improved absorbing boundary condition based on linear interpolation for ADI-FDTD method

With the linear interpolation method, an improved absorbing boundary condition（ABC）is introduced and derived, which is suitable for the alternating-direction-implicit finite- difference time-domain （ADI-FDTD） method. The reflection of the ABC caused by both the truncated error and the phase velocity error is analyzed. Based on the phase velocity estimation and the nonuniform cell, two methods are studied and then adopted to improve the performance of the ABC. A calculation case of a rectangular waveguide which is a typical dispersive transmission line is carried out using the ADI-FDTD method with the improved ABC for evaluation. According to the calculated case, the comparison is given between the reflection coefficients of the ABC with and without the velocity estimation and also the comparison between the reflection coefficients of the ABC with and without the nonuniform processing. The reflection variation of the ABC under different time steps is also analyzed and the acceptable worsening will not obscure the improvement on the absorption. Numerical results obviously show that efficient improvement on the absorbing performance of the ABC is achieved based on these methods for the ADI-FDTD.

The implementation of an improved NPML absorbing boundary condition in elastic wave modeling

In elastic wave forward modeling, absorbing boundary conditions （ABC） are used to mitigate undesired reflections from the model truncation boundaries. The perfectly matched layer （PML） has proved to be the best available ABC. However, the traditional splitting PML （SPML） ABC has some serious disadvantages： for example, global SPML ABCs require much more computing memory, although the implementation is easy. The implementation of local SPML ABCs also has some difficulties, since edges and corners must be considered. The traditional non-splitting perfectly matched layer （NPML） ABC has complex computation because of the convolution. In this paper, based on non-splitting perfectly matched layer （NPML） ABCs combined with the complex frequency-shifted stretching function （CFS）, we introduce a novel numerical implementation method for PML absorbing boundary conditions with simple calculation equations, small memory requirement, and easy programming.

Hybrid absorbing boundary condition for threedimensional elastic wave modeling

Edge reflections are inevitable in numerical modeling of seismic wavefields, and they are usually attenuated by absorbing boundary conditions. However, the commonly used perfectly matched layer （PML） boundary condition requires special treatment for the absorbing zone, and in three-dimensional （3D） modeling, it has to split each variable into three corresponding variables, which increases the computing time and memory storage. In contrast, the hybrid absorbing boundary condition （HABC） has the advantages such as ease of implementation, less computation time, and near-perfect absorption; it is thus able to enhance the computational efficiency of 3D elastic wave modeling. In this study, a HABC is developed from two-dimensional （2D） modeling into 3D modeling based on the I st Higdon one way wave equations, and a HABC is proposed that is suitable for a 3D elastic wave numerical simulation. Numerical simulation results for a homogenous model and a complex model indicate that the proposed HABC method is more effective and has better absorption than the traditional PML method.

A note for the mechanism of high-frequency oscillation instability resulted from absorbing boundary conditions

In this paper the explanation of the mechanism of high-frequency oscillation instability resulted from absorbing boundary conditions is further improved. And we analytically prove the proposition that for one dimensional discrete model of elastic wave motion, the module of reflection factor will be greater than 1 in high frequency band when artificial wave velocity is greater than 1.5 times the ratio of discrete space step to discrete time step. Based on the proof, the frequency band in which instability occurs is discussed in detail, showing such high-frequency waves are meaningless for the numerical simulation of wave motion.

Absorbing Boundary Conditions for Hyperbolic Systems

This paper deals with absorbing boundary conditions for hyperbolic systems in one and two space dimensions.We prove the strict well-posedness of the resulting initial boundary value problem in 1D.Afterwards we establish the GKS-stability of the corresponding Lax-Wendroff-type finite difference scheme.Hereby,we have to extend the classical proofs,since the(discretized) absorbing boundary conditions do not fit the standard form of boundary conditions for hyperbolic systems.

A Novel Absorbing Boundary Condition for the Frequency-Dependent Finite-Difference Time-Domain Method

A new absorbing boundary condition (ABC) for frequency dependent finite difference time domain algorithm for the arbitrary dispersive media is presented. The concepts of the digital systems are introduced to the (FD) 2TD method. On the basis of digital filter designing and vector algebra, the absorbing boundary condition under arbitrary angle of incidence are derived. The transient electromagnetic problems in two dimensions and three dimensions are calculated and the validity of the ABC is verified.

A truncated implicit high-order finite-difference scheme combined with boundary conditions

In this paper, first we calculate finite-difference coefficients of implicit finite- difference methods （IFDM） for the first and second-order derivatives on normal grids and first- order derivatives on staggered grids and find that small coefficients of high-order IFDMs exist. Dispersion analysis demonstrates that omitting these small coefficients can retain approximately the same order accuracy but greatly reduce computational costs. Then, we introduce a mirrorimage symmetric boundary condition to improve IFDMs accuracy and stability and adopt the hybrid absorbing boundary condition （ABC） to reduce unwanted reflections from the model boundary. Last, we give elastic wave modeling examples for homogeneous and heterogeneous models to demonstrate the advantages of the proposed scheme.

Time Dependent Wave Propagation Modeling Using Finite Difference Scheme of 2D Wave Equation Based on Absorbing and Reflecting Boundaries

Boundary procedure is an important phenomenon in numerical simulation. To reduce or eliminate the spurious reflections significantly which is occurred in boundary is a challenging and vital approach. The appropriate artificial numerical boundaries can be applied to eliminate the effect of unnecessary spurious reflections in case of the numerical simulations of wave propagation phenomena problems. Typically, to reduce the artificial reflections, the absorbing boundary conditions are necessary. In this paper, we overview and investigate the appropriate typical absorbing boundary conditions and analyzed the boundary effect of two dimensional wave equation numerically. Reflections over the wide-ranging incident angles are complicated to eliminate, but the absorbing boundary conditions that we have applied are computationally cost efficient, easy to apply and able to reduce reflections significantly. For numerical solution, finite difference method is applied to develop numerical scheme using 2D wave equation. Using the developed numerical scheme, we obtain the numerical solution of the governing equation as an initial boundary value problem and realize the qualitative behavior of the solution in infinite space. The finite difference numerical scheme has been investigated by developing MATLAB programming language code. Numerical results have been discussed and analyzed with presenting different qualitative behavior of the numerical scheme. The accuracy and efficiency of the numerical scheme has been illustrated. The stability analysis was discussed and verified stability condition. Using the numerical scheme and absorbing boundary conditions, the boundary effects and absorption of spurious reflection of boundary have been demonstrated.

STABILITY FOR IMPOSING ABSORBING BOUNDARY CONDITIONS IN THE FINITE ELEMENT SIMULATION OF ACOUSTIC WAVE PROPAGATION＊

It is well-known that artificial boundary conditions are crucial for the efficient and accurate computations of wavefields on unbounded domains. In this paper, we investigate stability analysis for the wave equation coupled with the first and the second order absorbing boundary conditions. The computational scheme is also developed. The approach allows the absorbing boundary conditions to be naturally imposed, which makes it easier for us to construct high order schemes for the absorbing boundary conditions. A thirdorder Lagrange finite element method with mass lumping is applied to obtain the spatial discretization of the wave equation. The resulting scheme is stable and is very efficient since no matrix inversion is needed at each time step. Moreover, we have shown both abstract and explicit conditional stability results for the fully-discrete schemes. The results are helpful for designing computational parameters in computations. Numerical computations are illustrated to show the efficiency and accuracy of our method. In particular, essentially no boundary reflection is seen at the artificial boundaries.

PML Absorbing Boundary Condition for Seismic Numerical Modeling by Convolutional Differentiator in Fluid-Saturated Porous Media

The perfectly matched layer（PML） was first introduced by Berenger as an absorbing boundary condition for electromagnetic wave propagation.In this article,a method is developed to ex-tend the PML to simulating seismic wave propagation in fluid-saturated porous medium.This non-physical boundary is used at the computational edge of a Forsyte polynomial convolutional differenti-ator（FPCD） algorithm as an absorbing boundary condition to truncate unbounded media.The incor-poration of PML in Biot＇s equations is given.Numerical results show that the PML absorbing bound-ary condition attenuates the outgoing waves effectively and eliminates the reflections adequately.

HIGH-ORDER LOCAL ABSORBING BOUNDARY CONDITIONS FOR HEAT EQUATION IN UNBOUNDED DOMAINS

With the development of numerical methods the numerical computations require higher and higher accuracy. This paper is devoted to the high-order local absorbing boundary conditions （ABCs） for heat equation. We proved that the coupled system yields a stable problem between the obtained high-order local ABCs and the partial differential equation in the computational domain. This method has been used widely in wave propagation models only recently. We extend the spirit of the methodology to parabolic ones, which will become a basis to design the local ABCs for a class of nonlinear PDEs. Some numerical tests show that the new treatment is very efficient and tractable.

Implementation of absorbing boundary conditions in dynamic simulation of the material point method

Outgoing waves arising from high-velocity impacts between soil and structure can be reflected by the conventional truncated boundaries.Absorbing boundary conditions(ABCs),to attenuate the energy of the outward waves,are necessary to ensure the proper representation of the kinematic field and the accurate quantification of impact forces.In this paper,damping layer and dashpot ABCs are implemented in the material point method(MPM)with slight adjustments.Benchmark scenarios of different dynamic problems are modelled with the ABCs configured.Feasibility of the ABCs is assessed through the velocity fluctuations at specific observation points and the impact force fluctuations on the structures.The impact forces predicted by the MPM with ABCs are verified by comparison with those estimated using a computational fluid dynamics approach.

Software for Acoustic Design

We present the СATEС software, which implements the solution to the problems of computational acoustics. The software is based on the use of the superelement method and finite element modeling algorithms, in-cluding hydrodynamic noise. The paper presents the main possibilities of software for solving acoustic design problems. .

Practical Absorbing Boundary Conditions forWave Propagation on Arbitrary Domain

This paper presents an absorbing boundary conditions(ABCs)for wave propagations on arbitrary computational domains.The purpose of ABCs is to eliminate the unwanted spurious reflection at the artificial boundaries and minimize the finite size effect.Traditional methods are usually complicate in theoretical derivation and implementation and work only for very limited types of boundary geometry.In contrast to other existing methods,our emphasis is placed on the ease of implementation.In particular,we propose a method for which the implementation can be done by fitting or learning from the simulation data in a larger domain,and it is insensitive to the geometry and space dimension of the computational domain.Furthermore,a stability criterion is imposed to ensure the stability of the proposed ABC.Numerical results are presented to demonstrate the effectiveness of our method.

STABILITY ANALYSIS OF FINITE ELEMENT METHODS FORTHE ACOUSTIC WAVE EQUATION WITH ABSORBINGBOUNDARY CONDITIONS (PART I)

In Part I and Part II of this paper initial-boundary value problems of the acoustic wave equation with absorbing boundary conditions are considered. Their finite element-finite difference computational schemes are proposed. The stability of the schemes is discussed and the corresponding stability conditions are given. Part I and Part II concern the first- and the second-order absorbing boundary conditions, respectively. Finally, numerical results are presented in Part II to show the correctness of theoretical analysis. (Author abstract) 7 Refs.

THE LONG WAVE OSCILLATIONS IN HARBORS WITH ABSORBING BOUNDARY CONDITIONS

This paper employs finite element method to solve shallow water equations with absorbing boundary conditions(the third kind,mixed boundary conditions).It is of practical importance in the cases that the land boundaries of the coastal area are made of porous medium allowing sea water flow in or out.The absorbing boundary conditions are treated as natural boundary conditions in wave equation finite element model.The numerical results for rectangu- lar and quarterly annular harbors indicate that the numerical solutions agree very well with ana- lytic solutions,which are also given in this paper.It is found that the land boundary absorbabili- ty may be significant to long wave oscillations,such as tidal waves in coastal harbors.

Solution of the time-dependent Schrdinger equation with absorbing boundary conditions

The performances of absorbing boundary conditions （ABCs） in four widely used finite difference time domain （FDTD） methods, i.e. explicit, implicit, explicit staggered-time, and Chebyshev methods, for solving the time-dependent Schrodinger equation are assessed and compared. The computation efficiency for each approach is also evaluated. A typical evolution problem of a single Gaussian wave packet is chosen to demonstrate the performances of the four methods combined with ABCs. It is found that ABCs perfectly eliminate reflection in implicit and explicit staggered-time methods. However, small reflection still exists in explicit and Chebyshev methods even though ABCs are applied.