相关论文: Robust preconditioners for perturbed saddle-point …
We derive parameter-robust quasi-optimal error estimates for mixed finite element methods for the nonlinear Darcy--Forchheimer equations with mixed boundary conditions. Using the framework of operator preconditioning, we also design…
In many applications, one wants to model physical systems consisting of two different physical processes in two different domains that are coupled across a common interface. A crucial challenge is then that the solutions of the two…
We consider the dynamic Biot model (see [Biot, M. A. J. Appl. Phys. 33, 1482--1498 (1962)]) describing the interaction between fluid flow and solid deformation including wave propagation phenomena in both the liquid and solid phases of a…
We consider the problem of iteratively solving large and sparse double saddle-point systems arising from the stationary Stokes-Darcy equations in two dimensions, discretized by the Marker-and-Cell (MAC) finite difference method. We analyze…
The goal of this paper is to propose preconditioners for the system of linear equations that arises from a discretization of fourth order elliptic problems using spectral element methods. These preconditioners are constructed using…
In this paper, we propose a distributed computing approach to solving large-scale robust stability problems on the simplex. Our approach is to formulate the robust stability problem as an optimization problem with polynomial variables and…
We consider a coupled model of free-flow and porous medium flow, governed by stationary Stokes and Darcy flow, respectively. The coupling between the two systems is enforced by introducing a single variable representing the normal flux…
In this paper, we propose and analyze a mixed formulation for the Kelvin-Voigt-Brinkman-Forchheimer equations for unsteady viscoelastic flows in porous media. Besides the velocity and pressure, our approach introduces the vorticity as a…
A second-order-accurate finite volume method, hybridized by blending an extended double-flux algorithm and a traditionally conservative scheme, is developed. In this scheme, hybrid convective fluxes as well as hybrid interpolation…
Simulating infiltration in porous media using Richards' equation remains computationally challenging due to its parabolic structure and nonlinear coefficients. While a wide range of numerical methods for differential equations have been…
Computational formulations for large strain, polyconvex, nearly incompressible elasticity have been extensively studied, but research on enhancing solution schemes that offer better tradeoffs between accuracy, robustness, and computational…
This study examines the dynamics of the third body in an elliptic restricted three-body problem (ERTBP) framework, taking into account perturbations from radiation pressure, oblateness, and elongation of the primary bodies, as well as…
We consider a saddle point formulation for a sixth order partial differential equation and its finite element approximation, for two sets of boundary conditions. We follow the Ciarlet-Raviart formulation for the biharmonic problem to…
This work develops optimal preconditioners for the discrete H(curl) and H(div) problems on two-dimensional surfaces by nodal auxiliary space preconditioning [R. Hiptmair, J. Xu: SIAM J. Numer. Anal. \textbf{45}, 2483-2509 (2007)]. In…
In this paper we propose a variant of the substructuring preconditioner for solving three-dimensional elliptic-type equations with strongly discontinuous coefficients. In the proposed preconditioner, we use the simplest coarse solver…
The relaxed physical factorization (RPF) preconditioner is a recent algorithm allowing for the efficient and robust solution to the block linear systems arising from the three-field displacement-velocity-pressure formulation of coupled…
Based on the needs of convergence proofs of preconditioned proximal point methods, we introduce notions of partial strong submonotonicity and partial (metric) subregularity of set-valued maps. We study relationships between these two…
In this paper, we propose a descent method for composite optimization problems with linear operators. Specifically, we first design a structure-exploiting preconditioner tailored to the linear operator so that the resulting preconditioned…
The discretization of Cahn-Hilliard equation with obstacle potential leads to a block 2 by 2 non-linear system, where the p1, 1q block has a non-linear and non-smooth term. Recently a globally convergent Newton Schur method was proposed for…
We study mixed finite element methods for the linearized rotating shallow water equations with linear drag and forcing terms. By means of a strong energy estimate for an equivalent second-order formulation for the linearized momentum, we…