A joint absorbing boundary for the multiple-relaxation-time lattice Boltzmann method in seismic acoustic wavefield modeling

Conventional seismic wave forward simulation generally uses mathematical means to solve the macroscopic wave equation,and then obtains the corresponding seismic wavefield.Usually,when the subsurface structure is finely constructed and the continuity of media is poor,this strategy is difficult to meet the requirements of accurate wavefield calculation.This paper uses the multiple-relaxation-time lattice Boltzmann method(MRT-LBM)to conduct the seismic acoustic wavefield simulation and verify its computational accuracy.To cope with the problem of severe reflections at the truncated boundaries,we analogize the viscous absorbing boundary and perfectly matched layer(PML)absorbing boundary based on the single-relaxation-time lattice Boltzmann(SRT-LB)equation to the MRT-LB equation,and further,propose a joint absorbing boundary through comparative analysis.We give the specific forms of the modified MRT-LB equation loaded with the joint absorbing boundary in the two-dimensional(2D)and three-dimensional(3D)cases,respectively.Then,we verify the effects of this absorbing boundary scheme on a 2D homogeneous model,2D modified British Petroleum(BP)gas-cloud model,and 3D homogeneous model,respectively.The results reveal that by comparing with the viscous absorbing boundary and PML absorbing boundary,the joint absorbing boundary has the best absorption performance,although it is a little bit complicated.Therefore,this joint absorbing boundary better solves the problem of truncated boundary reflections of MRT-LBM in simulating seismic acoustic wavefields,which is pivotal to its wide application in the field of exploration seismology.

Hybrid absorbing boundary condition based on transmitting boundary and its application in 3D fractional viscoacoustic modeling

An accurate numerical simulation for wave equations is essential for understanding of wave propagation in the earth's interior as well as full waveform inversion and reverse time migration. However, due to computational cost and hardware capability limitations, numerical simulations are often performed within a finite domain. Thus, an adequate absorbing boundary condition (ABC) is indispensable for obtaining accurate numerical simulation results. In this study, we develop a hybrid ABC based on a transmitting boundary, which is referred to as THABC, to eliminate artificial boundary reflections in 3D second-order fractional viscoacoustic numerical simulations. Furthermore, we propose an adaptive weighted coefficient to reconcile the transmitting and viscoacoustic wavefields in THABC. Through several numerical examples, we determine that the proposed THABC approach is characterized by the following benefits. First, with the same number of absorbing layers, THABC exhibits a better ability in eliminating boundary reflection than traditional ABC schemes. Second, THABC is more effective in computation, since it only requires the wavefields at the current and last time steps to solve the transmitting formula within the absorbing layers. Benefiting from a simple but effective combination between the transmitting equation and the second-order wave equation, our scheme performs well in the 3D fractional Laplacian viscoacoustic numerical simulation.

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.

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.

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.

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.

Well-Posedness of Initial Boundary Value ProblemsAssociated with Absorbing Boundary Conditions

In this paper we get one-way wave equations by using pseudo-differential operator theory,and present a set of absorbing boundary conditions based on the higher order aPProximations of oneway wave equations. An integral identity is the key point of the approximation. Also, we have provedthe well-posedness of the initial boundary value Problems related to our absorbing boundary conditionsconstructed in this artical.

Absorbing Boundary Conditions for Solving N-Dimensional Stationary Schr¨odinger Equations with Unbounded Potentials and Nonlinearities

We propose a hierarchy of novel absorbing boundary conditions for the onedimensional stationary Schr¨odinger equation with general(linear and nonlinear)potential.The accuracy of the new absorbing boundary conditions is investigated numerically for the computation of energies and ground-states for linear and nonlinear Schr¨odinger equations.It turns out that these absorbing boundary conditions and their variants lead to a higher accuracy than the usual Dirichlet boundary condition.Finally,we give the extension of these ABCs to N-dimensional stationary Schr¨odinger equations.

Analysis of High-Order Absorbing Boundary Conditions for the Schrodinger Equation

The paper is concerned with the numerical solution of Schr¨odinger equations on an unbounded spatial domain.High-order absorbing boundary conditions for one-dimensional domain are derived,and the stability of the reduced initial boundary value problem in the computational interval is proved by energy estimate.Then a second order finite difference scheme is proposed,and the convergence of the scheme is established as well.Finally,numerical examples are reported to confirm our error estimates of the numerical methods.

An Exact Absorbing Boundary Condition for the Schr¨odinger Equation With Sinusoidal Potentials at Infinity

In this paper we study numerical issues related to the Schr ¨odinger equationwith sinusoidal potentials at infinity. An exact absorbing boundary condition in a formof Dirichlet-to-Neumann mapping is derived. This boundary condition is based on ananalytical expression of the logarithmic derivative of the Floquet solution toMathieu’sequation, which is completely new to the author’s knowledge. The implementationof this exact boundary condition is discussed, and a fast evaluation method is used toreduce the computation burden arising from the involved half-order derivative operator.Some numerical tests are given to showthe performance of the proposed absorbingboundary conditions.

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.

An efficient absorbing boundary for finite-difference time-domain field modelling in acoustics

A highly efficient absorbing boundary condition suitable for use in the finitedtherence time-domain (FDTD) modelling of acoustic fields is presented in this paper. The new method seeks a least square esthoate of a transfer matrix for field components near truncating boundaries by matrir pseud-inversion. The proposed absorbing boundary is considerably more effective than most ekisting ones. The method is also computationally econondcal and robust.The performance of the new method is shown by numerical experiments on a point-source radiation problem, a wedge dimaction problem, and a scattering problem in which a plane wave is scattered by a circular cylinder.

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.

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

Analysis of high-frequency local coupling instability induced by multi- transmitting formula - P-SV wave simulation in a 2D waveguide

Local coupling instability will occur when the numerical scheme of absorbing boundary condition and that of the field wave equation allow energies to spontaneously enter into the computational domain. That is, the two schemes support common wave solutions with group velocity pointed into the computation domain. The key to eliminate local coupling instability is to avoid such wave solutions. For lumped-mass finite element simulation of P-SV wave motion in a 2D waveguide, an approach for stable implementation of high order multi-transmitting formula is provided. With a uniform rectangular mesh, it is proven and validated that high-freqaency local coupling instability can be eliminated by setting the ratio of the element size equal to or greater than x/2 times the ratio of the P wave velocity to the S wave velocity. These results can be valuable for dealing instability problems induced by other absorbing boundary conditions.

An improved convolution perfectly matched layer for elastic second-order wave equation

A convolution perfectly matched layer(CPML)can efficiently absorb boundary reflection in numerical simulation.However,the CPML is suitable for the first-order elastic wave equation and is difficult to apply directly to the second-order elastic wave equation.In view of this,based on the first-order CPML absorbing boundary condition,we propose a new CPML(NCPML)boundary which can be directly applied to the second-order wave equation.We first systematically extend the first-order CPML technique into second-order wave equations,neglecting the space-varying characteristics of the partial damping coefficient in the complex-frequency domain,avoiding the generation of convolution in the time domain.We then transform the technique back to the time domain through the inverse Fourier transform.Numerical simulation indicates that the space-varying characteristics of the attenuation factor have little influence on the absorption effect and increase the memory at the same time.A number of numerical examples show that the NCPML proposed in this study is effective in simulating elastic wave propagation,and this algorithm is more efficient and requires less memory allocation than the conventional PML absorbing boundary.