4 research outputs found
Parameter-robust discretization and preconditioning of Biot's consolidation model
Biot's consolidation model in poroelasticity has a number of applications in
science, medicine, and engineering. The model depends on various parameters,
and in practical applications these parameters ranges over several orders of
magnitude. A current challenge is to design discretization techniques and
solution algorithms that are well behaved with respect to these variations. The
purpose of this paper is to study finite element discretizations of this model
and construct block diagonal preconditioners for the discrete Biot systems. The
approach taken here is to consider the stability of the problem in non-standard
or weighted Hilbert spaces and employ the operator preconditioning approach. We
derive preconditioners that are robust with respect to both the variations of
the parameters and the mesh refinement. The parameters of interest are small
time-step sizes, large bulk and shear moduli, and small hydraulic conductivity.Comment: 24 page
Three-field block preconditioners for models of coupled magma/mantle dynamics
For a prescribed porosity, the coupled magma/mantle flow equations can be
formulated as a two-field system of equations with velocity and pressure as
unknowns. Previous work has shown that while optimal preconditioners for the
two-field formulation can be obtained, the construction of preconditioners that
are uniform with respect to model parameters is difficult. This limits the
applicability of two-field preconditioners in certain regimes of practical
interest. We address this issue by reformulating the governing equations as a
three-field problem, which removes a term that was problematic in the two-field
formulation in favour of an additional equation for a pressure-like field. For
the three-field problem, we develop and analyse new preconditioners and we show
numerically that they are optimal in terms of problem size and less sensitive
to model parameters, compared to the two-field preconditioner. This extends the
applicability of optimal preconditioners for coupled mantle/magma dynamics into
parameter regimes of physical interest