A method of solving the stiffness problem in Biot's poroelastic equations using a staggered high-order finite-difference
Zhao Hai-Bo(赵海波)a)†, Wang Xiu-Ming(王秀明)a)b), and Chen Hao(陈浩)a)
a Institute of Acoustics, Chinese Academy of Sciences, Beijing 100080, China; b Australian Resources Research Centre, Australian Commonwealth Scientific and Industrial Research Organization, Bentley WA 6102, Australia
Abstract In modelling elastic wave propagation in a porous medium, when the ratio between the fluid viscosity and the medium permeability is comparatively large, the stiffness problem of Biot's poroelastic equations will be encountered. In the paper, a partition method is developed to solve the stiffness problem with a staggered high-order finite-difference. The method splits the Biot equations into two systems. One is stiff, and solved analytically, the other is nonstiff, and solved numerically by using a high-order staggered-grid finite-difference scheme. The time step is determined by the staggered finite-difference algorithm in solving the nonstiff equations, thus a coarse time step may be employed. Therefore, the computation efficiency and computational stability are improved greatly. Also a perfect by matched layer technology is used in the split method as absorbing boundary conditions. The numerical results are compared with the analytical results and those obtained from the conventional staggered-grid finite-difference method in a homogeneous model, respectively. They are in good agreement with each other. Finally, a slightly more complex model is investigated and compared with related equivalent model to illustrate the good performance of the staggered-grid finite-difference scheme in the partition method.
Received: 13 April 2006
Revised: 29 April 2006
Accepted manuscript online:
Fund: Project supported by the ``100 Talents Project" of the Chinese Academy of Sciences and the Major Program of the National Natural Science Foundation of China (Grant No 10534040).
Cite this article:
Zhao Hai-Bo(赵海波), Wang Xiu-Ming(王秀明), and Chen Hao(陈浩) A method of solving the stiffness problem in Biot's poroelastic equations using a staggered high-order finite-difference 2006 Chinese Physics 15 2819
Altmetric calculates a score based on the online attention an article receives. Each coloured thread in the circle represents a different type of online attention. The number in the centre is the Altmetric score. Social media and mainstream news media are the main sources that calculate the score. Reference managers such as Mendeley are also tracked but do not contribute to the score. Older articles often score higher because they have had more time to get noticed. To account for this, Altmetric has included the context data for other articles of a similar age.