Fluid–structure interaction

From Wikipedia, the free encyclopedia
  (Redirected from Fluid-structure interaction)
Jump to: navigation, search

Fluid–structure interaction (FSI) is the interaction of some movable or deformable structure with an internal or surrounding fluid flow.[1] Fluid–structure interactions can be stable or oscillatory. In oscillatory interactions, the strain induced in the solid structure causes it to move such that the source of strain is reduced, and the structure returns to its former state only for the process to repeat.

Propagation of a pressure wave through an incompressible fluid in a flexible tube

Examples[edit]

Fluid–structure interactions are a crucial consideration in the design of many engineering systems, e.g. aircraft and bridges. Failing to consider the effects of oscillatory interactions can be catastrophic, especially in structures comprising materials susceptible to fatigue. Tacoma Narrows Bridge (1940), the first Tacoma Narrows Bridge, is probably one of the most infamous examples of large-scale failure. Aircraft wings and turbine blades can break due to FSI oscillations. Fluid–structure interaction has to be taken into account for the analysis of aneurysms in large arteries and artificial heart valves. A reed actually produces sound because the system of equations governing its dynamics has oscillatory solutions. The dynamic of reed valves used in two strokes engines and compressors is governed by FSI. The act of "blowing a raspberry" is another such example.

Analysis[edit]

Fluid–structure interaction problems and multiphysics problems in general are often too complex to solve analytically and so they have to be analyzed by means of experiments or numerical simulation. Research in the fields of computational fluid dynamics and computational structural dynamics is still ongoing but the maturity of these fields enables numerical simulation of fluid-structure interaction. Two main approaches exist for the simulation of fluid–structure interaction problems:

  • Monolithic approach: the equations governing the flow and the displacement of the structure are solved simultaneously, with a single solver
  • Partitioned approach: the equations governing the flow and the displacement of the structure are solved separately, with two distinct solvers

The monolithic approach requires a code developed for this particular combination of physical problems whereas the partitioned approach preserves software modularity because an existing flow solver and structural solver are coupled. Moreover, the partitioned approach facilitates solution of the flow equations and the structural equations with different, possibly more efficient techniques which have been developed specifically for either flow equations or structural equations. On the other hand, development of stable and accurate coupling algorithm is required in partitioned simulations.

Numerical simulation[edit]

The Newton–Raphson method or a different fixed-point iteration can be used to solve FSI problems. Methods based on Newton–Raphson iteration are used in both the monolithic [2][3][4] and the partitioned [5][6] approach. These methods solve the nonlinear flow equations and the structural equations in the entire fluid and solid domain with the Newton–Raphson method. The system of linear equations within the Newton–Raphson iteration can be solved without knowledge of the Jacobian with a matrix-free iterative method, using a finite difference approximation of the Jacobian-vector product.

Whereas Newton–Raphson methods solve the flow and structural problem for the state in the entire fluid and solid domain, it is also possible to reformulate an FSI problem as a system with only the degrees of freedom in the interface’s position as unknowns. This domain decomposition condenses the error of the FSI problem into a subspace related to the interface.[7] The FSI problem can hence be written as either a root finding problem or a fixed point problem, with the interface’s position as unknowns.

Interface Newton–Raphson methods solve this root-finding problem with Newton–Raphson iterations, e.g. with an approximation of the Jacobian from a linear reduced-physics model.[8][9] The interface quasi-Newton method with approximation for the inverse of the Jacobian from a least-squares model couples a black-box flow solver and structural solver [10] by means of the information that has been gathered during the coupling iterations. This technique is based on the interface block quasi-Newton technique with an approximation for the Jacobians from least-squares models which reformulates the FSI problem as a system of equations with both the interface’s position and the stress distribution on the interface as unknowns. This system is solved with block quasi-Newton iterations of the Gauss–Seidel type and the Jacobians of the flow solver and structural solver are approximated by means of least-squares models.[11]

The fixed-point problem can be solved with fixed-point iterations, also called (block) Gauss–Seidel iterations,[6] which means that the flow problem and structural problem are solved successively until the change is smaller than the convergence criterion. However, the iterations converge slowly if at all, especially when the interaction between the fluid and the structure is strong due to a high fluid/structure density ratio or the incompressibility of the fluid.[12] The convergence of the fixed point iterations can be stabilized and accelerated by Aitken relaxation and steepest descent relaxation, which adapt the relaxation factor in each iteration based on the previous iterations.[13]

If the interaction between the fluid and the structure is weak, only one fixed-point iteration is required within each time step. These so-called staggered or loosely coupled methods do not enforce the equilibrium on the fluid–structure interface within a time step but they are suitable for the simulation of aeroelasticity with a heavy and rather stiff structure. Several studies have analyzed the stability of partitioned algorithms for the simulation of fluid-structure interaction[12][14] [15] .[16][17][18]

See also[edit]

Open source codes[edit]

Academic Codes[edit]

Commercial Codes[edit]

References[edit]

  1. ^ Bungartz, Hans-Joachim; Schäfer, Michael, eds. (2006). Fluid-structure Interaction: Modelling, Simulation, Optimization. Springer-Verlag. ISBN 3-540-34595-7. 
  2. ^ M. Heil (2004). "An efficient solver for the fully coupled solution of large-displacement fluid-structure interaction problems". Computer Methods in Applied Mechanics and Engineering 193: 1–23. Bibcode:2004CMAME.193....1H. doi:10.1016/j.cma.2003.09.006. 
  3. ^ K.-J. Bathe; H. Zhang (2004). "Finite element developments for general fluid flows with structural interactions". International Journal for Numerical Methods in Engineering 60: 213–232. Bibcode:2004IJNME..60..213B. doi:10.1002/nme.959. 
  4. ^ J. Hron, S. Turek (2006). H.-J. Bungartz; M. Schäfer, eds. A monolithic FEM/multigrid solver for ALE formulation of fluid-structure interaction with application in biomechanics. Lecture Notes in Computational Science and Engineering. Fluid–Structure Interaction – Modelling, Simulation, Optimisation. Springer-Verlag. pp. 146–170. ISBN 3-540-34595-7. 
  5. ^ H. Matthies, J. Steindorf (2003). "Partitioned strong coupling algorithms for fluid-structure interaction". Computers and Structures 81: 805–812. doi:10.1016/S0045-7949(02)00409-1. 
  6. ^ a b H. Matthies, R. Niekamp, J. Steindorf (2006). "Algorithms for strong coupling procedures". Computer Methods in Applied Mechanics and Engineering 195: 2028–2049. Bibcode:2006CMAME.195.2028M. doi:10.1016/j.cma.2004.11.032. 
  7. ^ C. Michler, E. van Brummelen, R. de Borst (2006). "Error-amplification analysis of subiteration-preconditioned GMRES for fluid-structure interaction". Computer Methods in Applied Mechanics and Engineering 195: 2124–2148. Bibcode:2006CMAME.195.2124M. doi:10.1016/j.cma.2005.01.018. 
  8. ^ J.-F. Gerbeau, M. Vidrascu (2003). "A quasi-Newton algorithm based on a reduced model for fluid-structure interaction problems in blood flows". ESAIM: Mathematical Modelling and Numerical Analysis 37 (4): 631–648. doi:10.1051/m2an:2003049. 
  9. ^ J.-F. Gerbeau, M. Vidrascu, P. Frey (2005). "Fluid-structure interaction in blood flows on geometries based on medical imaging". Computers and Structures 83 (2-3): 155–165. doi:10.1016/j.compstruc.2004.03.083. 
  10. ^ J. Degroote, K.-J. Bathe, J. Vierendeels (2009). "Performance of a new partitioned procedure versus a monolithic procedure in fluid–structure interaction". Computers and Structures 87: 793. doi:10.1016/j.compstruc.2008.11.013. 
  11. ^ J. Vierendeels, L. Lanoye, J. Degroote, P. Verdonck (2007). "Implicit coupling of partitioned fluid-structure interaction problems with reduced order models". Computers and Structures 85 (11–14): 970–976. doi:10.1016/j.compstruc.2006.11.006. 
  12. ^ a b P. Causin, J.-F. Gerbeau, F. Nobile (2005). "Added-mass effect in the design of partitioned algorithms for fluid-structure problems". Computer Methods in Applied Mechanics and Engineering 194 (42-44): 4506–4527. Bibcode:2005CMAME.194.4506C. doi:10.1016/j.cma.2004.12.005. 
  13. ^ U. Küttler, W. Wall (2008). "Fixed-point fluid-structure interaction solvers with dynamic relaxation". Computational Mechanics 43 (1): 61–72. Bibcode:2008CompM..43...61K. doi:10.1007/s00466-008-0255-5. 
  14. ^ J. Degroote, P. Bruggeman, R. Haelterman, J. Vierendeels (2008). "Stability of a coupling technique for partitioned solvers in FSI applications". Computers and Structures 86 (23–24): 2224–2234. doi:10.1016/j.compstruc.2008.05.005. 
  15. ^ R. Jaiman, X. Jiao, P. Geubelle, E. Loth (2006). "Conservative load transfer along curved fluid-solid interface with non-matching meshes". Journal of Computational Physics 218 (1): 372–397. Bibcode:2006JCoPh.218..372J. doi:10.1016/j.jcp.2006.02.016. 
  16. ^ J. Vierendeels, K. Dumont, E. Dick, P. Verdonck (2005). "Analysis and stabilization of fluid-structure interaction algorithm for rigid-body motion". AIAA Journal 43 (12): 2549–2557. Bibcode:2005AIAAJ..43.2549V. doi:10.2514/1.3660. 
  17. ^ Christiane Förster, Wolfgang A. Wall, Ekkehard Ramm (2006). P. Wesseling, E. Oñate, J. Périaux, ed. "The artificial added mass effect in sequential staggered fluid-structure interaction algorithms". European Conference on Computational Fluid Dynamics ECCOMAS CFD 2006. The Netherlands. 
  18. ^ Christiane Förster, Wolfgang A. Wall, Ekkehard Ramm (2007). "Artificial added mass instabilities in sequential staggered coupling of nonlinear structures and incompressible viscous flows". Computer Methods in Applied Mechanics and Engineering 196: 1278–1293. Bibcode:2007CMAME.196.1278F. doi:10.1016/j.cma.2006.09.002.