An official website of the United States government
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
Analysis of a 3D nonlinear moving boundary problem describing fluid-mesh-sell interaction
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
- Award ID(s):
- 1853340
- PAR ID:
- 10148592
- Date Published:
- Journal Name:
- Transactions of the American Mathematical Society
- ISSN:
- 2330-0000
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract We consider the finite element approximation of a coupled fluid‐structure interaction (FSI) system, which comprises a three‐dimensional (3D) Stokes flow and a two‐dimensional (2D) fourth‐order Euler–Bernoulli or Kirchhoff plate. The interaction of these parabolic and hyperbolic partial differential equations (PDE) occurs at the boundary interface which is assumed to be fixed. The vertical displacement of the plate dynamics evolves on the flat portion of the boundary where the coupling conditions are implemented via the matching velocities of the plate and fluid flow, as well as the Dirichlet boundary trace of the pressure. This pressure term also acts as a coupling agent, since it appears as a forcing term on the flat, elastic plate domain. Our main focus in this work is to generate some numerical results concerning the approximate solutions to the FSI model. For this, we propose a numerical algorithm that sequentially solves the fluid and plate subsystems through an effective decoupling approach. Numerical results of test problems are presented to illustrate the performance of the proposed method.more » « less
-
Underwater explosion poses a significant threat to the structural integrity of ocean vehicles and platforms. Accurate prediction of the dynamic loads from an explosion and the resulting structural response is crucial to ensuring safety without overconservative design. When the distance between the explosive charge and the structure is relatively small (i.e., near-field explosion), the dynamics of the gaseous explosion product, i.e., the “bubble”, comes into play, rendering a multiphysics problem that features the interaction of the bubble, the surrounding liquid water, and the solid structure. The problem is highly nonlinear, as it involves shock waves, large deformation, yielding, contact, and possibly fracture. This paper investigates the two-way interaction between the cyclic expansion and collapse of an explosion bubble and the deformation of a thin-walled elastoplastic cylindrical shell in its vicinity. Intuitively, when a shock wave impinges on a thin cylindrical shell, the shell would collapse in the direction of shock propagation. However, some recent laboratory experiments have shown that under certain conditions the shell collapsed in a counter-intuitive mode in which the direction of collapse is perpendicular to that of shock propagation. In other words, the nearest point on the structural surface moved towards the explosion charge, despite being impacted by a compressive shock. This paper focuses on replicating this phenomenon through numerical simulation and elucidating the underlying mechanisms. A recently developed computational framework (“FIVER”) coupling a nonlinear finite element structural dynamics solver and a finite volume compressible fluid dynamics solver is used to complete this study. The solver utilizes an embedded boundary method to track the wetted surface of the structure (i.e. the fluid-structure interface), which is capable of handling large structural deformation and topological changes (e.g., fracture). The solver also adopts the level set method for tracking the bubble surface (i.e. the liquid-gas interface). The fluid-structure and liquid-gas interface conditions are enforced by constructing and solving one-dimensional multi-material Riemann problems, which naturally accommodates the propagation of shock waves across the interfaces. In this paper, mesh refinement study is made to examine the sensitivity of the results to various meshing parameters. The results show that the intermediate level of refinement is appropriate in terms of both the accuracy and the computation costs. Next, the deformation history of both the bubble and the structure are presented and analyzed to provide a detailed view of the counter-intuitive collapse mode mentioned above. We show that timewise, the structural collapse spans multiple cycles of bubble oscillation. Additional details about the time-histories of fluid pressure, structure displacement, and bubble size are presented to elucidate this dynamic bubble-structure interaction and the resulting structural failure.more » « less
-
This paper introduces the 3D Peskin problem: a two-dimensional elastic membrane immersed in a three-dimensional steady Stokes flow. We obtain the equations that model this free boundary problem and show that they admit a boundary integral reduction, providing an evolution equation for the elastic interface. We consider general nonlinear elastic laws, i.e. the fully nonlinear Peskin problem, and prove that the problem is well-posed in low-regularity Hölder spaces. Moreover, we prove that the elastic membrane becomes smooth instantly in time.more » « less
-
Aeromechanics of highly flexible flapping wings is a complex nonlinear fluid–structure interaction problem and, therefore, cannot be analyzed using conventional linear aeroelasticity methods. This paper presents a standalone coupled aeroelastic framework for highly flexible flapping wings in hover for micro air vehicle (MAV) applications. The MAV-scale flapping wing structure is modeled using fully nonlinear beam and shell finite elements. A potential-flow-based unsteady aerodynamic model is then coupled with the structural model to generate the coupled aeroelastic framework. Both the structural and aerodynamic models are validated independently before coupling. Instantaneous lift force and wing deflection predictions from the coupled aeroelastic simulations are compared with the force and deflection measurements (using digital image correlation) obtained from in-house flapping wing experiments at both moderate (13 Hz) and high (20 Hz) flapping frequencies. Coupled trim analysis is then performed by simultaneously solving wing response equations and vehicle trim equations until trim controls, wing elastic response, inflow and circulation converge all together. The dependence of control inputs on weight and center of gravity (cg) location of the vehicle is studied for the hovering flight case.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- The DOI is not currently available.
