3,687 research outputs found
A Lagrange multiplier method for a Stokes-Biot fluid-poroelastic structure interaction model
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
Discontinuous Galerkin Methods for Mass Transfer through Semi-Permeable Membranes
A discontinuous Galerkin (dG) method for the numerical solution of
initial/boundary value multi-compartment partial differential equation (PDE)
models, interconnected with interface conditions, is presented and analysed.
The study of interface problems is motivated by models of mass transfer of
solutes through semi-permeable membranes. More specifically, a model problem
consisting of a system of semilinear parabolic advection-diffusion-reaction
partial differential equations in each compartment, equipped with respective
initial and boundary conditions, is considered. Nonlinear interface conditions
modelling selective permeability, congestion and partial reflection are applied
to the compartment interfaces. An interior penalty dG method is presented for
this problem and it is analysed in the space-discrete setting. The a priori
analysis shows that the method yields optimal a priori bounds, provided the
exact solution is sufficiently smooth. Numerical experiments indicate agreement
with the theoretical bounds and highlight the stability of the numerical method
in the advection-dominated regime
- ā¦