A high-order in space spectral-element methodology for the solution of a strongly coupled fluid-structure interaction (FSI) problem is developed. A methodology is based on a partitioned solution of incompressible fluid equations on body-fitted grids, and nonlinearly-elastic solid deformation equations coupled via a fixed-point iteration approach with Aitken relaxation. A comprehensive verification strategy of the developed methodology is presented, including h-, p-and temporal refinement studies. An expected order of convergence is demonstrated first separately for the corresponding fluid and solid solvers, followed by a self-convergence study on a coupled FSI problem (self-convergence refers to a convergence to a reference solution obtained with the same solver at higher resolution). To this end, a new three-dimensional fluid-structure interaction benchmark is proposed for a verification of the FSI codes, which consists of a fluid flow in a channel with one rigid and one flexible wall. It is shown that, due to a consistent problem formulation, including initial and boundary conditions, a high-order spatial convergence on a fully coupled FSI problem can be demonstrated. Finally, a developed framework is applied successfully to a Direct Numerical Simulation of a turbulent flow in a channel interacting with a compliant wall, where the fluid-structure interface is fully resolved.
more »
« less
Fluid-structure interaction simulations with a LES filtering approach in solids4Foam
Abstract The goal of this paper is to test solids4Foam, the fluid-structure interaction (FSI) toolbox developed for foam-extend (a branch of OpenFOAM), and assess its flexibility in handling more complex flows. For this purpose, we consider the interaction of an incompressible fluid described by a Leray model with a hyperelastic structure modeled as a Saint Venant-Kirchho material. We focus on a strongly coupled, partitioned fluid-structure interaction (FSI) solver in a finite volume environment, combined with an arbitrary Lagrangian-Eulerian approach to deal with the motion of the fluid domain. For the implementation of the Leray model, which features a nonlinear differential low-pass filter, we adopt a three-step algorithm called Evolve-Filter-Relax. We validate our approach against numerical data available in the literature for the 3D cross flow past a cantilever beam at Reynolds number 100 and 400.
more »
« less
- Award ID(s):
- 1953535
- PAR ID:
- 10333506
- Date Published:
- Journal Name:
- Communications in Applied and Industrial Mathematics
- Volume:
- 12
- Issue:
- 1
- ISSN:
- 2038-0909
- Page Range / eLocation ID:
- 13 to 28
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
null (Ed.)Abstract Hydroplaning is a phenomenon that occurs when a layer of water between the tire and pavement pushes the tire upward. The tire detaches from the pavement, preventing it from providing sufficient forces and moments for the vehicle to respond to driver control inputs such as breaking, accelerating, and steering. This work is mainly focused on the tire and its interaction with the pavement to address hydroplaning. Using a tire model that is validated based on results found in the literature, fluid–structure interaction (FSI) between the tire-water-road surfaces is investigated through two approaches. In the first approach, the coupled Eulerian–Lagrangian (CEL) formulation was used. The drawback associated with the CEL method is the laminar assumption and that the behavior of the fluid at length scales smaller than the smallest element size is not captured. To improve the simulation results, in the second approach, an FSI model incorporating finite element methods (FEMs) and the Navier–Stokes equations for a two-phase flow of water and air, and the shear stress transport k–ω turbulence model, was developed and validated, improving the prediction of real hydroplaning scenarios. With large computational and processing requirements, a grid dependence study was conducted for the tire simulations to minimize the mesh size yet retain numerical accuracy. The improved FSI model was applied to hydroplaning speed and cornering force scenarios.more » « less
-
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
-
As a step towards addressing a scarcity of references on this topic, we compared the Eulerian and Lagrangian Computational Fluid Dynamics (CFD) approaches for the solution of free-surface and Fluid–Solid Interaction (FSI) problems. The Eulerian approach uses the Finite Element Method (FEM) to spatially discretize the Navier–Stokes equations. The free surface is handled via the volume-of-fluid (VOF) and the level-set (LS) equations; an Immersed Boundary Method (IBM) in conjunction with the Nitsche’s technique were applied to resolve the fluid–solid coupling. For the Lagrangian approach, the smoothed particle hydrodynamics (SPH) method is the meshless discretization technique of choice; no additional equations are needed to handle free-surface or FSI coupling. We compared the two approaches for a flow around cylinder. The dam break test was used to gauge the performance for free-surface flows. Lastly, the two approaches were compared on two FSI problems—one with a floating rigid body dropped into the fluid and one with an elastic gate interacting with the flow. We conclude with a discussion of the robustness, ease of model setup, and versatility of the two approaches. The Eulerian and Lagrangian solvers used in this study are open-source and available in the public domain.more » « less
-
Human phonation involves the flow-induced vibrations of the vocal folds (VFs) that result from the interaction with airflow through the larynx. Most voice dysfunctions correspond with the fluid–structure interaction (FSI) features as well as the local changes in perfusion within the VF tissue. This study aims to develop a multiphysics computational framework to simulate the interstitial fluid flow dynamics in vibrating VFs using a biphasic description of the tissue and FSI methodology. The integration of FSI and a permeable VF model presents a novel approach to capture phonation physics' complexity and investigate VF tissue's porous nature. The glottal airflow is modeled by the unsteady, incompressible Navier–Stokes equations, and the Brinkman equation is employed to simulate the flow through the saturated porous medium of the VFs. The computational model provides a prediction of tissue deformation metrics and pulsatile glottal flow, in addition to the interstitial fluid velocity and flow circulation within the porous structure. Furthermore, the model is used to characterize the effects of variation in subglottal lung pressure and VF permeability coefficient by conducting parametric studies. Subsequent investigations to quantify the relationships between these input variables, flow perfusion, pore pressure, and vibration amplitude are presented. A linear relationship is found between the vibration amplitude, pore pressure, and filtration flow with subglottal pressure, whereas a nonlinear dependence between the filtration velocity and VF permeability coefficient is detected. The outcomes highlight the importance of poroelasticity in phonation models.more » « less
An official website of the United States government

