skip to main content


The NSF Public Access Repository (NSF-PAR) system and access will be unavailable from 5:00 PM ET until 11:00 PM ET on Friday, June 21 due to maintenance. We apologize for the inconvenience.

This content will become publicly available on December 1, 2025

Title: Fluid-poroviscoelastic structure interaction problem with nonlinear coupling
We prove the existence of a weak solution to a fluid-structure interaction (FSI) problem between the flow of an incompressible, viscous fluid modeled by the Navier-Stokes equations, and a poroviscoelastic medium modeled by the Biot equations. The two are nonlinearly coupled over an interface with mass and elastic energy, modeled by a reticular plate equation, which is transparent to fluid flow. The existence proof is constructive, consisting of two steps. First, the existence of a weak solution to a regularized problem is shown. Next, a weak-classical consistency result is obtained, showing that the weak solution to the regularized problem converges, as the regularization parameter approaches zero, to a classical solution to the original problem, when such a classical solution exists. While the assumptions in the first step only require the Biot medium to be poroelastic, the second step requires additional regularity, namely, that the Biot medium is poroviscoelastic. This is the first weak solution existence result for an FSI problem with nonlinear coupling involving a Biot model for poro(visco)elastic media.  more » « less
Award ID(s):
Author(s) / Creator(s):
Publisher / Repository:
Date Published:
Journal Name:
Jorunal de Mathematiques Pures and Appliques
Not yet availabe
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. Abstract. We consider a nonlinear, moving boundary, fluid-structure interaction problem between a time dependent incompressible, viscous fluid flow, and an elastic structure composed of a cylindrical shell supported by a mesh of elastic rods. The fluid flow is modeled by the time-dependent Navier- Stokes equations in a three-dimensional cylindrical domain, while the lateral wall of the cylinder is modeled by the two-dimensional linearly elastic Koiter shell equations coupled to a one-dimensional system of conservation laws defined on a graph domain, describing a mesh of curved rods. The mesh supported shell allows displacements in all three spatial directions. Two-way coupling based on kinematic and dynamic coupling conditions is assumed between the fluid and composite structure, and between the mesh of curved rods and Koiter shell. Problems of this type arise in many ap- plications, including blood flow through arteries treated with vascular prostheses called stents. We prove the existence of a weak solution to this nonlinear, moving boundary problem by using the time discretization via Lie operator splitting method combined with an Arbitrary Lagrangian-Eulerian approach, and a non-trivial extension of the Aubin-Lions-Simon compactness result to problems on moving domains. 
    more » « less
  3. Abstract

    This work focuses on modeling the interaction between an incompressible, viscous fluid and a poroviscoelastic material. The fluid flow is described using the time‐dependent Stokes equations, and the poroelastic material using the Biot model. The viscoelasticity is incorporated in the equations using a linear Kelvin–Voigt model. We introduce two novel, noniterative, partitioned numerical schemes for the coupled problem. The first method uses the second‐order backward differentiation formula (BDF2) for implicit integration, while treating the interface terms explicitly using a second‐order extrapolation formula. The second method is the Crank–Nicolson and Leap‐Frog (CNLF) method, where the Crank–Nicolson method is used to implicitly advance the solution in time, while the coupling terms are explicitly approximated by the Leap‐Frog integration. We show that the BDF2 method is unconditionally stable and uniformly stable in time, while the CNLF method is stable under a CFL condition. Both schemes are validated using numerical simulations. Second‐order convergence in time is observed for both methods. Simulations over a longer period of time show that the errors in the solution remain bounded. Cases when the structure is poroviscoelastic and poroelastic are included in numerical examples.

    more » « less
  4. In this paper we introduce a constructive approach to study well-posedness of solutions to stochastic uid-structure interaction with stochastic noise. We focus on a benchmark problem in stochastic uidstructure interaction, and prove the existence of a unique weak solution in the probabilistically strong sense. The benchmark problem consists of the 2D time-dependent Stokes equations describing the ow of an incompressible, viscous uid interacting with a linearly elastic membrane modeled by the 1D linear wave equation. The membrane is stochastically forced by the time-dependent white noise. The uid and the structure are linearly coupled. The constructive existence proof is based on a time-discretization via an operator splitting approach. This introduces a sequence of approximate solutions, which are random variables. We show the existence of a subsequence of approximate solutions which converges, almost surely, to a weak solution in the probabilistically strong sense. The proof is based on uniform energy estimates in terms of the expectation of the energy norms, which are the backbone for a weak compactness argument giving rise to a weakly convergent subsequence of probability measures associated with the approximate solutions. Probabilistic techniques based on the Skorohod representation theorem and the Gyongy-Krylov lemma are then employed to obtain almost sure convergence of a subsequence of the random approximate solutions to a weak solution in the probabilistically strong sense. The result shows that the deterministic benchmark FSI model is robust to stochastic noise, even in the presence of rough white noise in time. To the best of our knowledge, this is the  rst well-posedness result for stochastic uid-structure interaction. 
    more » « less
  5. Abstract We introduce and analyze a partially augmented fully mixed formulation and a mixed finite element method for the coupled problem arising in the interaction between a free fluid and a poroelastic medium. The flows in the free fluid and poroelastic regions are governed by the Navier–Stokes and Biot equations, respectively, and the transmission conditions are given by mass conservation, balance of fluid force, conservation of momentum and the Beavers–Joseph–Saffman condition. We apply dual-mixed formulations in both domains, where the symmetry of the Navier–Stokes and poroelastic stress tensors is imposed in an ultra-weak and weak sense. In turn, since the transmission conditions are essential in the fully mixed formulation, they are imposed weakly by introducing the traces of the structure velocity and the poroelastic medium pressure on the interface as the associated Lagrange multipliers. Furthermore, since the fluid convective term requires the velocity to live in a smaller space than usual, we augment the variational formulation with suitable Galerkin-type terms. Existence and uniqueness of a solution are established for the continuous weak formulation, as well as a semidiscrete continuous-in-time formulation with nonmatching grids, together with the corresponding stability bounds and error analysis with rates of convergence. Several numerical experiments are presented to verify the theoretical results and illustrate the performance of the method for applications to arterial flow and flow through a filter. 
    more » « less