Fluid/solid modeling and coupling
Mathematical models are chosen to describe the fluid dynamics on one hand and the solid dynamics on the other hand. For the fluid, this choice depends on two nondimensional numbers: the Mach number which determines whether the flow is compressible or incompressible, and the Reynolds number, which indicates whether the flow is laminar or turbulent. In the current state of the project, low Mach number flows are considered, so that the fluid motion is described as an incompressible viscous flow, and low and high Reynolds number flow cases are adressed. To describe the structural behaviour, two classes of models are considered: rigid solids mounted on springs, as classicaly done in aeroelasticity or civil engineering studies, and deformable structures made of hyperelastic materials.

Fluid modeling
We consider viscous Newtonian fluids and describe their motion in an Eulerian framework with the NavierStokes equations. Models of increasing complexity are used, depending on the physical parameters.
 Incompressible laminar flows
When the Reynolds number Re and Mach number Ma are small, the flow is incompressible and described by the laminar NavierStokes equations. The conservative form of the momentum equation is
where u_{E}(x_{t},t) is the Eulerian velocity field, i.e. the velocity of the fluid at the spatial position x_{t} of the fluid domain. In fluid/structure interaction problems, the fluid domain moves as a result of the solid motion and thus the spatial coordinate x_{t} is time dependent, as indicated by the subscript. This time dependence of the spatial domain is also indicated in the spatial derivative operator with the subscript t. For a Newtonian incompressible viscous fluid, the Cauchy stress tensor is defined as
Mathematically speaking, the pressure field p_{E} is a Lagrangian multiplier associated to the divergencefree constraint
 Incompressible turbulent flows
When the Mach number Ma is small but the Reynolds number Re is large, the flow is incompressible and turbulent. The flow velocity is then decomposed as the sum of a mean component and a fluctuating component, i.e.
The mean component is then governed by the incompressible Reynolds Averaged Navier Stokes (RANS) equations
which differs from the NavierSDtokes equations by an additional stress term. The socalled turbulence Reynolds stresses depend quadratically on the fluctuating velocity field, which is unknown. To close the problem, the turbulence Reynolds stresses are modelled as a function of the mean velocty and pressure component, solely. Using he Boussinesq approximation, which introduces an eddy viscosity field, the unsteady Reynolds Averaged NavierStokes (RANS) equations write
where is the turbulence eddy viscosity, determined with a turbulence model. The SpalartAllmaras turbulence model is a oneequation model that solves a transport equation for a surrogate viscosity.
 Compressible turbulent flows
When the Reynolds number Re and Mach number Ma are large, the flow is compressible and turbulent. It is then described by the compressible unsteady Reynolds Averaged NavierStokes (RANS) equations
\mut represents the turbulent viscosity to be computed through some RANS turbulence model (e.g. SpalartAllmaras).

Solid modeling
We introduce here the notation for modelling rigid structures mounted on spring and the structural model used for the nonlinear elastic material.
 Springmounted rigid cylinder
Add a scketch of the srpingmounted clyinder (similar to what Johann did)
The evolution equation of the vertical displacement of the mass center is a damped harmonic oscillator forced by the ycomponent of the aerodynamical force
where is the vertical displacement of the cylinder's mass center, is the structural frequency, defined as the square root of the spring's stiffness over the cylinder's mass, is the (structural) damping coefficient and is the solidtofluid denisty ratio.
 Translational and torsional springmounted rigid plates
A very simple modelisation of a moving solid consists in considering a rotating and translating rigd body. For example, when a rigid solid is allowed to translate in the vertical direction and rotate around some elastic axis, simple equations of motion based on harmonic oscillators are obtained :
is the lift force (positive upwards)
is the moment of fluid forces evaluated at the elastic center (chosen positive nosedown here)
is the static unbalance. Note that if the center of gravity is not aligned with the elastic axis, the heaving and pitching degrees of freedom are inertially coupled through the static unbalance term.
The following set of nondimensional numbers can be used to describe this system :
They represent respectively the reduced flow velocity (based on pitching characteristc timescale), the ratio of heaving to pitching frequencies, the solidtofluid density ratio and the nondimensional position of the center of gravity.
Ref : E. Dowell, H. Curtiss, R. Scanlan, A modern Course in Aeroelasticity, Fifth Revised and Enlarged Edition, Springer, 2015
 Linear and nonlinear elastic solids
The evolution of an elastic structure is defined by the motion of its material points. The position _{}of any material point in the actual configuration is defined with respect to its position in a reference configuration. In solid mechanics, a natural choice for this reference configuration is the spatial domain occupied by the solid when no external stresses are applied. We will denote the spatial position configuration in this reference configuration. Other choices for the reference configuration could be made, especially when investigating the stability analysis.
The material displacement is a Lagrangian field defined as the difference between the actual position and the reference position, i.e.
and the solid dynamics is described by the evolution equation of the displacement field
where is the solid density, is the deformation gradient of the solid defined as
and is the second PiolaKirchhoff stress tensor. The second PiolaKirchhoff stress is classically used in solid mechanics since it is a measure of the force acting on an element of area in the reference configuration, as opposed to the Cauchy stress tensor, classically used in fluid mechanics, since it is a measure of the force acting on an element of area in the actual configuration. The Cauchy and second PiolaKirchhoff stress tensors are related through the formula
where , the determinant of the deformation gradient
measures the volumic change of the solid and is equal to unity in case of an incompressible solid.
In the AEROFLEX project, we mainly consider hyperelastic materials modelled with the Saint VenantKirchhoff law, which relates the second PiolaKirchhoff stress tensor to the displacement field
via the GreenLagrangian strain tensor defined as
This is a model valid for small deformations and large displacements of an elastic solid. In case of small displacements, the quadratic (nonlinear) term in the GreenLagrangian stress tensor is negligible and the equation governing the solid dynamics reduces to the classical linear elasticity equation

Fluid/Elastic solid coupling
The fluid and solid dynamics are naturally described by variables of different nature. The fluid velocity is an Eulerian variable that describes the velocity of the fluid at a fixed spatial position in space, while the solid displacement is a Lagrangian variable that describes the displacement of a material point in a fixed reference configuration. The complexity in the modelling of fluid/structure interaction problem arises from the coupling between the Eulerian fluid and Lagrangian solid description.
Physicallt, the coupling of the fluid and solid dynamics coupling occurs at the fluid/structure interface through boundary conditions, that express
 the continuity of the fluid and solid velocities
 the continuity of the fluid and solid stresses
To cope with the coupling of the Eulerian and Lagrangian description, various strategies have emerged in the literature dealing with the numerical resolution of the unsteady nonlinear fluid/solid equation. They can be classified into two class of methods: the conforming methods and the nonconforming methods. In the present project, the coupling of the nonlinear fluid/solid equations is described either within an Arbitrary Lagrangian Eulerian (ALE) framework (conforming method) or with a Fictitious Domain approach (nonconforming method).
References:
G. Hou, J. Wang and A. Layton, Numerical Methods for Fluid Structure interaction – A review, Communication in Computational Physics, Vol 12, No 2, pp. 337377