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.

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.

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.

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.

APPLICATION OF THE PML ABSORBING BOUNDARY CONDITION TO FD-TD SIMULATION

The key problem of finite-difference time-domain (FD-TD) method is the skillful application of special conditions on the boundaries of the computational domain. A new technique named Perfectly Matched Layer(PML) yields a robust Absorbing Boundary Condition(ABC) independent of the angle of incidence and the frequency of outgoing waves. In this paper, the principle of the PML technique is briefly presented. Then some problems in the application and their settlements are discussed emphatically. Finally three numerical tests and a measured result are devoted to examine the accuracy and effectiveness of this approach.

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.

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 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.

High-Accurate Transparent Boundary Conditions for Time-Dependent Quantum Wave Packet Method

Complex absorbing potential is usually required in a time-dependent wave packet method to accomplish the calculation in a truncated region.Usually it works effectively but becomes inefficient when the wave function involves translational energy of broad range,particularly involving ultra-low energy.In this work,a new transparent boundary condition(TBC)is proposed for the time-dependent wave packet method.It in principle is of spectral accuracy when typical discrete variable representations are applied.The prominent merit of the new TBC is that its accuracy is insensitive to the translational energy distribution of the wave function,in contrast with the complex absorbing potential.Application of the new TBC is given to one-dimensional particle wave packet scatterings from a barrier with a potential well,which supports resonances states.

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.

发展方程的完美匹配边界层模型

在有界区域上使用数值方法模拟无界区域上的偏微分方程时,要求在有界区域边界上添加人工边界条件。构造的人工边界条件必须能消除边界上不必要的反射,完美匹配边界层是非常有效的人工吸收边界条件之一。本研究采用复值坐标拉伸的方法,在不分离变量的前提下对发展方程推导出了完美匹配边界层模型。通过平面波分析可知,波在初始求解区域边界上无反射,且在有界的完美匹配边界层上振幅是指数衰减的。

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.