This paper is concerned with numerical algorithms for Biot model. By introducing an intermediate variable, the classical 2-field Biot model is written into a 3-field formulation. Based on such a 3-field formulation, we propose a coupled algorithm, some time-extrapolation based decoupled algorithms, and an iterative decoupled algorithm. Our focus is the analysis of the iterative decoupled algorithm. It is shown that the convergence of the iterative decoupled algorithm requires no extra assumptions on physical parameters or stabilization parameters. Numerical experiments are provided to demonstrate the accuracy and efficiency of the proposed method.
more »
« less
This content will become publicly available on May 1, 2024
An Iterative Decoupled Algorithm with Unconditional Stability for Biot Model
This paper is concerned with numerical algorithms for Biot model. By introducing an intermediate variable, the classical 2-field Biot model is written into a 3-field formulation. Based on such a 3-field formulation, we propose a coupled algorithm, some time-extrapolation based decoupled algorithms, and an iterative decoupled algorithm. Our focus is the analysis of the iterative decoupled algorithm. It is shown that the convergence of the iterative decoupled algorithm requires no extra assumptions on physical parameters or stabilization parameters. Numerical experiments are provided to demonstrate the accuracy and efficiency of the proposed method.
more »
« less
- Award ID(s):
- 1700328
- NSF-PAR ID:
- 10406600
- Date Published:
- Journal Name:
- Mathematics of computation
- Volume:
- 92
- Issue:
- 341
- ISSN:
- 1088-6842
- Page Range / eLocation ID:
- 1087-1108
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
In this paper, we develop parameter-robust numerical algorithms for Biot model and apply the algorithms in brain edema simulations. By introducing an intermediate variable, we derive a multiphysics reformulation of the Biot model. Based on the reformulation, the Biot model is viewed as a generalized Stokes subproblem combining with a reaction–diffusion subproblem. Solving the two subproblems together or separately leads to a coupled or a decoupled algorithm. We conduct extensive numerical experiments to show that the two algorithms are robust with respect to the key physical parameters. The algorithms are applied to study the brain swelling caused by abnormal accumulation of cerebrospinal fluid in injured areas. The effects of the key physical parameters on brain swelling are carefully investigated. It is observed that the permeability has the biggest influence on intracranial pressure (ICP) and tissue deformation; the Young’s modulus and the Poisson ratio do not affect the maximum value of ICP too much but have big influence on the tissue deformation and the developing speed of brain swelling.more » « less
-
In this paper, we aim at solving the Biot model under stabilized finite element discretizations. To solve the resulting generalized saddle point linear systems, some iterative methods are proposed and compared. In the first method, we apply the GMRES algorithm as the outer iteration. In the second method, the Uzawa method with variable relaxation parameters is employed as the outer iteration method. In the third approach, Uzawa method is treated as a fixed-point iteration, the outer solver is the so-called Anderson acceleration. In all these methods, the inner solvers are preconditioners for the generalized saddle point problem. In the preconditioners, the Schur complement approximation is derived by using Fourier analysis approach. These preconditioners are implemented exactly or inexactly. Extensive experiments are given to justify the performance of the proposed preconditioners and to compare all the algorithms.more » « less
-
We develop a mixed finite element method for the coupled problem arising in the interaction between a free fluid governed by the Stokes equations and flow in deformable porous medium modeled by the Biot system of poroelasticity. Mass conservation, balance of stress, and the Beavers–Joseph–Saffman condition are imposed on the interface. We consider a fully mixed Biot formulation based on a weakly symmetric stress-displacement-rotation elasticity system and Darcy velocity-pressure flow formulation. A velocity-pressure formulation is used for the Stokes equations. The interface conditions are incorporated through the introduction of the traces of the structure velocity and the Darcy pressure as Lagrange multipliers. Existence and uniqueness of a solution are established for the continuous weak formulation. Stability and error estimates are derived for the semi-discrete continuous-in-time mixed finite element approximation. Numerical experiments are presented to verify the theoretical results and illustrate the robustness of the method with respect to the physical parameters.more » « less
-
more » « less
SUMMARY The advent of fast sensing technologies allow for real-time model updates in many applications where the model parameters are uncertain. Once the observations are collected, Bayesian algorithms offer a pathway for real-time inversion (a.k.a. model parameters/inputs update) because of the flexibility of the Bayesian framework against non-uniqueness and uncertainties. However, Bayesian algorithms rely on the repeated evaluation of the computational models and deep learning (DL) based proxies can be useful to address this computational bottleneck. In this paper, we study the effects of the approximate nature of the deep learned models and associated model errors during the inversion of borehole electromagnetic (EM) measurements, which are usually obtained from logging while drilling. We rely on the iterative ensemble smoothers as an effective algorithm for real-time inversion due to its parallel nature and relatively low computational cost. The real-time inversion of EM measurements is used to determine the subsurface geology and properties, which are critical for real-time adjustments of the well trajectory (geosteering). The use of deep neural network (DNN) as a forward model allows us to perform thousands of model evaluations within seconds, which is very useful to quantify uncertainties and non-uniqueness in real-time. While significant efforts are usually made to ensure the accuracy of the DL models, it is widely known that the DNNs can contain some type of model-error in the regions not covered by the training data, which are unknown and training specific. When the DL models are utilized during inversion of EM measurements, the effects of the model-errors could manifest themselves as a bias in the estimated input parameters and as a consequence might result in a low-quality geosteering decision. We present numerical results highlighting the challenges associated with the inversion of EM measurements while neglecting model-error. We further demonstrate the utility of a recently proposed flexible iterative ensemble smoother in reducing the effect of model-bias by capturing the unknown model-errors, thus improving the quality of the estimated subsurface properties for geosteering operation. Moreover, we describe a procedure for identifying inversion multimodality and propose possible solutions to alleviate it in real-time.
