A frequency-dependent absorbing boundary condition for numerically solving u-U elastic wave equations in layered and fluid-saturated porous media

Guoliang Zhang,Mi Zhao,Xiuli Du,Xiaoling Zhang
DOI: https://doi.org/10.1016/j.soildyn.2020.106189
IF: 4
2020-08-01
Soil Dynamics and Earthquake Engineering
Abstract:An accurate absorbing boundary condition (ABC) is proposed to incorporate in the frequency-domain finite element simulation for numerically solving u-U elastic wave equations in horizontally multilayer fluid-saturated porous media. The u-U elastic wave equations are first discretized only along the depth direction in which the material properties are heterogeneous, analytical method using in the remaining coordinate direction. A general eigenvalue problem on the horizontal wavenumber is then solved with the frequency as known parameter. Dynamic stiffness on the artificial boundary of finite domain is finally obtained as ABC based on the solution to the above eigenvalue problem. The proposed ABC is mathematically derived without introducing any approximate assumptions, and the only approximation comes from the finite-element discretization of the fluid-saturated porous medium in the depth direction. The proposed ABC can also be coupled seamlessly with the finite element method (FEM). Numerical examples of wave radiation problems are given to demonstrate the effectiveness of the proposed ABC and its coupling with FEM. Finally, a practical application of proposed ABC to assess the topography effects of the fluid-saturated porous basin-type site is presented to show that the proposed ABC can also be used to solve wave scattering problems in geotechnical earthquake engineering.
geosciences, multidisciplinary,engineering, geological
What problem does this paper attempt to address?
The problem that this paper attempts to solve is how to impose an accurate absorbing boundary condition (ABC) on the artificial boundary of a finite domain when performing numerical simulation of wave propagation in fluid - saturated porous media, in order to simulate the radiation damping effect of an infinite domain. Specifically, the paper proposes a new frequency - dependent absorbing boundary condition for solving the numerical solution problem of the elastic wave equation in horizontally multi - layered fluid - saturated porous media. This method aims to improve the accuracy of numerical simulation, reduce the influence of spurious reflected waves, and thus more realistically simulate the wave propagation process, especially in the application of seismic engineering. The paper mentions that traditional absorbing boundary conditions have made a lot of progress in dealing with wave propagation problems in single - phase media, but there are relatively few studies on fluid - saturated porous media. Most of the existing methods are extended from the absorbing boundary conditions of single - phase media. Although these methods are easy to implement, they are usually approximate and require the establishment of a large enough finite - element model to dissipate spurious reflected waves. Therefore, the paper proposes a new accurate global absorbing boundary condition. By solving an eigenvalue problem, the dynamic stiffness on the artificial boundary is obtained, and then the absorbing boundary condition is constructed. This method not only improves the accuracy of numerical simulation but also can be seamlessly combined with the finite - element method and is suitable for frequency - domain analysis. The paper also verifies the effectiveness of the proposed absorbing boundary condition through numerical examples and shows its application in evaluating the topographic effect of fluid - saturated porous basins, proving the potential of this method in solving wave - scattering problems in geotechnical engineering.