Numerical Methods for Fluid-Structure Interaction …

1 Commun. Comput. : , No. 2, pp. 337-377 August 2012 REVIEWARTICLEN umerical Methods for fluid -StructureInteraction A ReviewGene Hou1, Jin Wang2, and Anita Layton31 Department of Mechanical and Aerospace Engineering, Old Dominion University,Norfolk, VA 23529, of Mathematics and Statistics, Old Dominion University, Norfolk, VA23529, of Mathematics, Duke University, Durham, NC 27708, 29 December 2010; Accepted (in revised version) 29 April 2011 Available online 20 February interactions between incompressible fluid flows and immersed struc-tures are nonlinear multi-physics phenomena that have applications to a wide range ofscientific and engineering disciplines.

2 In this article, wereview representative numeri-cal Methods based on conforming and non-conforming meshes that are currently avail-able for computing Fluid-Structure Interaction problems, with an emphasis on some ofthe recent developments in the field. A goal is to categorize the selected Methods andassess their accuracy and efficiency. We discuss challengesfaced by researchers in thisfield, and we emphasize the importance of interdisciplinaryeffort for advancing thestudy in Fluid-Structure subject classifications: 65-02, 65Z05, 74F10 Key words: Fluid-Structure Interaction , conforming and non-conforming meshes, immersed Introduction3382 FSI problem formulation3413 Conforming-mesh methods3434 FSI computation using immersed methods3545 Discussion366 Corresponding Layton) 2012 Global-Science Press338G.

3 Hou, J. Wang and A. Layton / Commun. Comput. Phys.,12(2012), pp. 337-3771 IntroductionIn Fluid-Structure Interaction (FSI) problems, one or more solid structures interact with aninternal or surrounding fluid flow. FSI problems play prominent roles in many scientificand engineering fields, yet a comprehensive study of such problems remains a challengedue to their strong nonlinearity and multidisciplinary nature (Chakrabarti 2005, Dowelland Hall 2001, Morand and Ohayon 1995). For most FSI problems,analytical solutions tothe model equations are impossible to obtain, whereas laboratory experiments are limitedin scope; thus to investigate the fundamental physics involved in the complex interactionbetween fluids and solids, Numerical simulations may be recent advances of computer technology, simulations of scientific and engineer-ing systems have become increasingly sophisticated and complicated.

4 For example, thespeed requirement of a planing boat hull has advanced to sucha degree and with sucha speed that has outpaced the availability of testing data and existing design equations(Weymouth et al. 2006, 2008). To fill the technological gap, anefficient Numerical algo-rithm can be used to investigate in detail the Interaction between water waves and themotion of the boat. Such an investigation is typically multidisciplinary. In this example,the performance of the boat is a result of the Interaction between water hydrodynamicsand structural dynamics.

5 Other FSI applications include, but are not limited to, sedimen-tation (Mucha et al. 2004, Tornberg and Shelley 2004, Wang andLayton 2009), particleassembly (Liu et al. 2006), aerodynamics (Haase 2001, Zhang, Jiang and Ye 2007), turbu-lence (Kaligzin and Iaccarino 2003, Yang and Balaras 2006), complex flows in irregular do-mains (Fadlun et al. 2000, Udaykumar et al. 1996, 2001), electro-hydrodynamics (Hoburgand Melcher 1976), magneto-hydrodynamic flows (Grigoriadiset al. 2009), biofluid andbio-mechanics (such as cell aggregation and deformation, blood-heart Interaction , innerear fluid dynamics, jellyfish swimming, sperm motility, cilliary beating, etc.)

6 The Numerical procedures to solve these FSI problems may be broadly classified intotwo approaches: themonolithic approachand thepartitioned approach. It is understoodthat the distinction between the monolithic and partitioned approaches may be vieweddifferently by researchers from different fields. In this paper, we intend to define thesetwo approaches from the engineering application point of view. Fig. 1 illustrates thesolution procedures of the monolithic and partitioned monolithic approach (Hubner et al. 2004, Michler et al. 2004, Ryzhakov et ) treats the fluid and structure dynamics in the same mathematical framework toform a single system equation for the entire problem, which is solved simultaneously bya unified algorithm.

7 The interfacial conditions are implicitin the solution procedure. Thisapproach can potentially achieve better accuracy for a multidisciplinary problem, but itmay require substantially more resources and expertise to develop and maintain such aspecialized code. In contrast, the partitioned approach treats the fluid and the structureas two computational fields which can be solved separately with their respective meshdiscretization and Numerical algorithm. The interfacial conditions are used explicitly tocommunicate information between the fluid and structure solutions.

8 A motivation ofG. Hou, J. Wang and A. Layton / Commun. Comput. Phys.,12(2012), pp. 337-377339 Figure 1: Schematic of the monolithic approach (a) and the partitioned approach (b) for Fluid-Structure inter-actions, whereSfandSsdenote the fluid and structure solutions, later approach is to integrate available disciplinary ( , fluidic and structural) algo-rithms and reduce the code development time by taking advantage of the legacy codesor Numerical algorithms that have been validated and used for solving many complicatedfluid or structural problems.

9 As a result, a successful partitioned method can solve a FSIproblem with sophisticated fluid and structural physics. The challenge of this approachis, however, to coordinate the disciplinary algorithms to achieve accurate and efficientfluid- structure Interaction solution with minimal code modification. Particularly, the in-terface location that divides the fluid and the structure domains is not known a priori andusually changed in time; thus, the partitioned approach requires the tracking of the newinterface location and its related quantities, which can becumbersome and general classification of the FSI solution procedures is based upon the treat-ment of meshes: theconforming mesh methodsandnon-conforming mesh Methods .

10 The con-forming mesh Methods consider the interface conditions as physical boundary condi-tions, which treat the interface location as part of the solution, and requires meshes thatconform to the interface. Owing to the movement and/or deformation of the solid struc-ture, re-meshing (or mesh-updating) is needed as the solution is advanced. On the otherhand, the non-conforming mesh Methods treat the boundary location and the related in-terface conditions as constraints imposed on the model equations so that non-conformingmeshes can be employed.