Ilona Ambartsumyan

Ilona Ambartsumyan, Eldar Khattatov, Ivan Yotov et al.

We study a finite element computational model for solving the coupled problem arising in the interaction between a free fluid and a fluid in a poroelastic medium. The free fluid is governed by the Stokes equations, while the flow in the poroelastic medium is modeled using the Biot poroelasticity system. Equilibrium and kinematic conditions are imposed on the interface. A mixed Darcy formulation is employed, resulting in continuity of flux condition of essential type. A Lagrange multiplier method is employed to impose weakly this condition. A stability and error analysis is performed for the semi-discrete continuous-in-time and the fully discrete formulations. A series of numerical experiments is presented to confirm the theoretical convergence rates and to study the applicability of the method to modeling physical phenomena and the sensitivity of the model with respect to its parameters.

Ilona Ambartsumyan, Vincent J. Ervin, Truong Nguyen et al.

We develop and analyze a model for the interaction of a quasi-Newtonian free fluid with a poroelastic medium. The flow in the fluid region is described by the nonlinear Stokes equations and in the poroelastic medium by the nonlinear quasi-static Biot model. Equilibrium and kinematic conditions are imposed on the interface. We establish existence and uniqueness of a solution to the weak formulation and its semidiscrete continuous-in-time finite element approximation. We present error analysis, complemented by numerical experiments.

Ilona Ambartsumyan, Eldar Khattatov, ChangQing Wang et al.

Three algorithms are developed for uncertainty quantification in modeling coupled Stokes and Darcy flows. The porous media may consist of multiple regions with different properties. The permeability is modeled as a non-stationary stochastic variable, with its log represented as a sum of local Karhunen-Loève (KL) expansions. The problem is approximated by stochastic collocation on either tensor-product or sparse grids, coupled with a multiscale mortar mixed finite element method for the spatial discretization. A non-overlapping domain decomposition algorithm reduces the global problem to a coarse scale mortar interface problem, which is solved by an iterative solver, for each stochastic realization. In the traditional implementation, each subdomain solves a local Dirichlet or Neumann problem in every interface iteration. To reduce this cost, two additional algorithms based on deterministic or stochastic multiscale flux basis are introduced. The basis consists of the local flux (or velocity trace) responses from each mortar degree of freedom. It is computed by each subdomain independently before the interface iteration begins. The use of the multiscale flux basis avoids the need for subdomain solves on each iteration. The deterministic basis is computed at each stochastic collocation and used only at this realization. The stochastic basis is formed by further looping over all local realizations of a subdomain's KL region before the stochastic collocation begins. It is reused over multiple realizations. Numerical tests are presented to illustrate the performance of the three algorithms, with the stochastic multiscale flux basis showing significant savings in computational cost.

Ilona Ambartsumyan, Eldar Khattatov, Jeonghun Lee et al.

We develop higher order multipoint flux mixed finite element (MFMFE) methods for solving elliptic problems on quadrilateral and hexahedral grids that reduce to cell-based pressure systems. The methods are based on a new family of mixed finite elements, which are enhanced Raviart-Thomas spaces with bubbles that are curls of specially chosen polynomials. The velocity degrees of freedom of the new spaces can be associated with the points of tensor-product Gauss-Lobatto quadrature rules, which allows for local velocity elimination and leads to a symmetric and positive definite cell-based system for the pressures. We prove optimal $k$-th order convergence for the velocity and pressure in their natural norms, as well as $(k+1)$-st order superconvergence for the pressure at the Gauss points. Moreover, local postprocessing gives a pressure that is superconvergent of order $(k+1)$ in the full $L^2$-norm. Numerical results illustrating the validity of our theoretical results are included.