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 compli-cated 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


Introduction
The solution of wave propagation phenomena in spatially discrete method for the numerical simulation always requires the elimination of the spurious events. The unnecessary spurious events have been produced by the boundaries of the numerical grids. In numerical simulation, a finite region of space is covered by the numerical mesh. As a result, these excessive boundary effects events arise in the system. The boundary effects have been appeared usually as reflection. These unwanted events are always inessential to the real physical events. As a result, the elimination of these events is desired in the study.
In recent years, study says that the explicit finite difference schemes are being used widely in finding the approximate solution for wave propagation problems.  [1]. The technique described here appears particularly advantageous in case of only long time integrations and the primary limitation of the technique is the propagation medium that must be uniform at the boundary. Appelö, Hagstrom, and Kreiss (2006) presented a model with perfectly matched layer boundary condition with the limitation of some special cases [2]. Benito, Urena, Gavete, Salete, and Muelas (2013) pointed out perfectly matched layer absorbing boundary condition in unbounded domain and shown the avoided reflection over the grid point from edges [3]. Clayton and Enquist (1977) described the elimination of boundary effects [4]. The findings remain some major difficulties in cooperating with corner points of the rectangular domain. Costen and Bérenger (2012) has shown the implementation of absorbing boundary condition in corner regions of computational domain [5]. Gächter and Grote (2003) derived an exact non-reflecting boundary condition which is investigated for three space dimensions [6]. Givoli and Neta (2003) incorporated a finite difference scheme and proposed a non-reflecting boundary scheme for time-dependent wave problems in unbounded domains [7]. The boundary effects along the boundary discussed here caused by normal and shear stress modules. Grote and Keller (2000) demonstrated finite difference method and showed high improvement in accuracy about the implementation of absorbing boundary conditions over standard methods [8]. Hagstrom, Mar-Or, and Givoli (2008) introduced artificial absorbing boundary condition to two dimensional wave guides [9]. Hall and Wang (2009) discussed both reflecting and absorbing boundary conditions in two dimensional models of heterogeneous, anisotropic and fractured media [10]. The artificial boundary is able to reduce the edges reflection but no significance approach is analyzed. Lysmer and Kuhlemeyer (1969) proposed a possible solution concentrating on viscous damping situation [11]. The proposed method can describe the energy propagation in case of outward direction only. However, this method almost reduces reflected compressional waves in boundary. It also works over a long range of incident angles. But, the method is not able to eliminate reflected shear waves completely. Mur (1981) presented absorbing boundary condition for maximum elimination of boundary effect [12]. We overviewed the Mur's absorbing boundary condition as more efficient and accurate in implementing boundary effect. Reynolds (1978) developed a finite difference models that shows unwanted reflections are produced from the edges of boundary [13]. This happens due to use of boundary conditions modeled by Dirichlet or Neumann.
The model was used for generating synthetic seismograms. Reynold in his study shows boundary conditions and the reduction of the edge reflection in a significant manner as well. Smith (1974) proposed a method which was claimed to work perfectly for all incident angles [14]. In this method, the simulation is ex-  [16].
In this study, a systematic method is presented through applying the set of absorbing boundary conditions. These conditions are noticed to be able to eliminate the unnecessary boundary events. The approached method is applied to the two dimensional wave equation in solving initial boundary value problems. The applied absorbing boundary conditions would be able to reduce or absorb the reflecting effect by the boundaries thorough the outward-moving energy.

Model and Method
In this study, two dimensional wave equation is considered to solve and analysis the boundary effects in the space-time domain is a rectangular spatial domain. The complete initial-boundary value problem is defined by where the four sides of the rectangle The initial condition (2) represents the initial shape of the incidence waves while another initial condition (3) specifies initial velocity of waves. Furthermore, in solving PDEs requires boundary conditions. The Equation (4) and (5) represents fixed boundary in the rectangular domain. Here, the fixed boundaries reflect the incidence waves at the boundary. Later, the investigation is carried out to eliminate the reflection from the boundary by introducing the absorbing boundary conditions. Waves in one-dimension move along the line. The reflection is occurred here by some particular change or disturbance in the line and consequently waves move backwards. Contrary, when we go through the higher dimensions to know the reflection effects and to know the wave disturbance occurred initially in some localized region and spreads it out is far from obvious. In this study, we initiate absorbing boundary conditions and shows how waves is being absorbed rather than reflected in the boundaries of the rectangular region. The phenomena can be accounted to be investigated like falling a pebble or other solid substance into still water create an outward forwarding circle of ripples. The boundary effects of the produced ripples are minimized in the study by introducing explicit finite difference method.

Finite Difference Discretization for 2D Wave Equation
The finite differences approximation is an extensively used method in solving partial differential equations. Considering, d space dimension and constant wave velocity , the wave equation can be written as ( ] Consider 2D wave equation with some free region which can be discretized as For constant mesh spacing let x y ∆ =∆ . The numerical scheme as follows ( )

Numerical Solution
In general, it is very difficult in most case to obtain the exact solution of partial differential equation. In some cases, though analytical solution can be obtained but that is in very complex form. Therefore, one needs to obtain the approximate solution of PDE. Finite difference method is one of the standard methods that have been used to find the approximate solution in solving partial differential equations. In this section, the graphical representation of wave propagation in two dimensions using the developed numerical scheme is presented for different time-step selection and investigated the boundary effects. The time-dependent numerical simulation for two dimension wave equation can be described by the above numerical results. Here, the four numerical results demonstrated by using the developed finite difference scheme. We consider the constant propagation speed 1 c = for each case. The temporal domain was considered as [0, 10]. The propagation of waves created from a source towards the boundary is shown here for 1st sec in Figure 1(a), 2nd sec in Figure 1(b), 5th sec in Figure 1(c) and finally at last 10th sec in Figure 1(d) respectively.

Stability Analysis of Numerical Scheme
The stability of numerical scheme depends on CFL condition. The CFL stands for the developed numerical scheme as c t x ∆ ∆ . Optimization of the plane sampling intervals ( ) , x y ∆ ∆ depends on some criteria while we select time sampling interval that also depends on some criteria as well. The numerical scheme shows Journal of Applied Mathematics and Physics Because the CFL value cannot be taken as negative value in the implementation, 1 2 c t x ∆ ≤ ∆ must be carried out to obtain the stable solution.
Here, Figure 2(a) shows the rate of propagation of waves towards the boundaries with different angle of incidence. The CFL condition does not exceed the limit criteria as we choose the CFL as 0.5 which is less than 1 2 . Contrary, Figure 2(b) shows no visible wave is being propagated here as CFL condition exceeds here. In this case, we select CFL value as 0.8 which is greater than 1 2 . So, the stability criteria 1 2 c t x ∆ ≤ ∆ is failed which results the instable graph. The investigation satisfied the stability criteria properly.

Reflection of Waves in Boundaries
The wave can be undergoing by reflection or be absorbed by the boundary. In the wave propagation model described in Equation (1)   Here, the Dirichlet boundary condition is applied and the reflections of waves in boundaries have been noticed. We investigate the effect of spurious reflection for different time step selection. We carried out the investigation for 10 sec. We noticed that, at 7 sec the waves hit the boundary and reflected backward which is shown in Figure 3(a). The reflection is growing significantly with time. Figure   3(b) shows the reflection at 8sec while Figure 3(c) gives the results of reflection at 9 sec. Finally, we experienced the maximum reflection at the boundary at 10sec in Figure 3(d).

Absorbing Boundary Condition for 2D Wave Equation
In this section we consider a finite two-dimensional plane. We reform the un-   For absorption across the boundary, the absorbing boundary conditions are The discretized form of the described absorbing condition is implemented in finite computational domain. The result is illustrated using MATLAB and investigated how the reflection is reduced or eliminated. The boundary conditions presented here are more stable and it is said to be computationally efficient as well. The numerical experiments and boundary effect of absorbing boundary condition is demonstrated in this study graphically. Here, we have presented the effect of non-reflecting and absorbing boundary condition for different time steps selection. After creating the waves from a source, the results of last four seconds have been shown here in this article. We see that, the waves are not reflected in each case as the propagation goes on. In Figure 4(a), we see that at 7 sec of the computation, the waves are being absorbed by the boundaries. Then, gradually in Figures 4(b)-(d), we have chosen the time selection for 8 sec, 9 sec and 10 sec respectively and found that all waves are being absorbed by the boundaries. The spurious reflections are being eliminated significantly here is noticed.

Conclusion
In numerical simulation for seismic wave fields, absorbing boundary condition plays a significant role in reducing boundary reflection. It is often essential to introduce artificial boundaries to limit the area of computation. Elimination of reflection effects by the boundaries may lead the controlling of different devastating disaster. The absorbing boundary condition and its effect are investigated in this paper. The applied boundary condition not only guarantees stable difference approximation, but also minimizes the spurious reflections that occur at the boundaries.