Numerical modelling for Fluid-structure interaction
3 likes1,517 views
(1) Numerical modeling of fluid-structure interaction (FSI) requires solving complex multi-physics problems at the fluid-solid interface. (2) Two common approaches to FSI modeling are body-fitted meshes, where the mesh conforms to the solid boundary, and unfitted/immersed methods, which use a fixed background mesh. (3) Coupling strategies for solving the fluid and solid sub-problems can be monolithic, solving them simultaneously, or staggered, where they are solved separately. Staggered methods are computationally more efficient but can suffer from added-mass instabilities.
Numerical modelling for Fluid-structure interaction
- 1. Numerical modelling for Fluid-structure interaction EGEM07 – Fluid-structure interaction Dr Wulf G. Dettmer Dr Chennakesava Kadapa Swansea University, UK.
- 2. Table of contents (1)Introduction (2)Aspects of numerical modelling for FSI (3)Body-fitted Vs Unfitted/immersed methods (4)Monolithic Vs Staggered schemes (5)A stabilised immersed framework for FSI
- 3. Introduction to FSI Interactions of fluid and solid A multi-physics phenomenon Abundant in nature – Almost every life form Occurs in many areas of engineering – Aerospace: Aircraft, parachutes, rockets – Civil: Bridges, dams, cable/roof structures – Mechanical: Automobiles, turbines, pumps – Naval: Ships, off-shore structures, submarines
- 4. Governing equations Fluid: (Eulerian) Solid: (Lagrangian) Interface: Kinematic condition: Equilibrium condition:
- 5. Aspects of numerical modelling Solid solver Fluid solver Coupling Can we solve all the FSI problems if we use the best available solvers for fluid and solid sub-problems?
- 6. Aspects of numerical modelling Solid solver Fluid solver Coupling Can we solve all the FSI problems if we use the best available solvers for fluid and solid sub-problems? No. But, why?
- 7. Aspects of numerical modelling Solid solver Fluid solver Coupling Can we solve all the FSI problems if we use the best available solvers for fluid and solid sub-problems? No. But, why? The devil is at the interface.
- 8. Caution! If someone tells you that his/her scheme/tool can solve a FSI problem without actually looking at the problem, then it is highly likely that he/she is lying.
- 9. Important properties of numerical schemes for FSI (1) Existence Does the tool have FSI capability? (2) Robustness For a reasonable time step, does the scheme work without crashing? (3) Accuracy How accurate is the solution? (4) Efficiency What is the amount of time required?
- 10. Body-fitted meshes Meshes aligned with the solid boundary Finite Element or Finite Volume schemes for the fluid problem How to deal with moving solids?
- 11. Body-fitted meshes When the solid moves Surrounding fluid mesh also moves Arbitrary Lagrangian-Eulerian (ALE) formulation for the fluid For small displacements mesh deformation schemes For large displacements re-meshing techniques
- 12. Body-fitted meshes Advantages Efficient and accurate for simple problems Well established Available in commercial software Disadvantages Mesh generation is cumbersome Require sophisticated re-meshing algorithms Complicated and inefficient in 3D Difficulty in capturing topological changes
- 13. Unfitted/immersed methods ● Solids immersed/embedded on fixed grids Advantages No need for body-fitted meshes No need for re-meshing Ideal for multi-phase flows, fracture Complex FSI problems can be solved Disadvantages Needs to develop a fluid solver Majority of the schemes are only 1st order accurate in time Very limited availability in commercial software
- 14. Integration in time Only implicit schemes are considered Fluid: 1st order - Backward Euler 2nd order – Crank-Nicolson/Trapezoidal, Generalised-alpha, BDF2 Solid: 1st order - Backward Euler 2nd order - Crank-Nicolson/Trapezoidal, Generalised-alpha
- 15. Coupling strategies Monolithic Vs Staggered ✔ Spatial discretisation ✔ Temporal discretisation
- 16. Governing equations Fluid: (Eulerian) Solid: (Lagrangian) Interface: Kinematic condition: Equilibrium condition:
- 17. Coupling Data transfer between fluid and solid Types of techniques Dirichlet-Neumann (body-fitted, unfitted) Robin-Robin (body-fitted, unfitted) Body-force (standard Immersed methods) We consider Dirichlet-Neumann The most intuitive and physical Velocity boundary condition on the Fluid Force boundary condition on the Solid
- 18. Monolithic schemes • Fixed-point or Newton-Raphson • Advantages ➔ No added-mass instabilities ➔ 2nd order accuracy in time is possible • Disadvantages ➔ Need to develop customised solvers ➔ Computationally expensive ➔ Difficult to linearise ➔ Convergence issues
- 19. Staggered schemes • Solve solid and fluid separately • Advantages ➔ Computationally appealing ➔ Existing solvers can be used ● Disadvantages ➔ Added-mass instabilities ➔ Difficult to get 2nd order accurate schemes for FSI with flexible structures in the presence of significant added-mass ➔ Efficiency and accuracy decrease with the increase in added mas
- 20. Summary of FSI schemes Body-fitted Unfitted Monolithic Staggered Commerical software No added-mass issue Expensive Efficient Easiest of all Added-mass issues No added-mass issue Complicated Expensive Efficient Relatively easy Many applications Added-mass issue
- 21. What is added mass issue? Instability arising when 1) The density of the solid is close to or less than that of the fluid Blood flow through arteries 2) When the structure is very thin Shell structures 3) When the structure is highly flexible Roof membranes, parachutes
- 22. A model problem for FSI d d Dettmer, W. G. and Peric, D. A new staggered scheme for fluid-structure interaction, IJNME, 93, 1-22, 2013.
- 23. A stabilised immersed framework for FSI Combines the state-of-the-art Hierarchical b-splines SUPG/PSPG stabilisation for the fluid Ghost-penalty stabilisation for cut-cells Solid-Solid contact Staggered solution schemes Wide variety of applications
- 24. B-Splines and hierarchical refinement - spatial discretisations for unfitted meshes
- 25. B-Splines
- 29. B-Splines • Nice mathematical properties • Tensor product nature • Partition of unity • Higher-order continuities across element boundaries • Always positive • No hanging nodes • Ease of localised refinements • Efficient programming techniques and data structures
- 30. Sample meshes
- 31. Formulation Incompressible Navier-Stokes Variational formulation Time integration: Backward Euler (O(dt)) and Generalised-alpha (O(dt^2))
- 34. Sedimentation of multiple particles
- 35. Model turbines
- 36. Ball check valve
- 37. Relief valve in 3D
- 38. References (1) W. G. Dettmer and D. Perić. A new staggered scheme for fluid-structure interaction, IJNME, 93, 1-22, 2013. (2) W. G. Dettmer, C. Kadapa, D. Perić, A stabilised immersed boundary method on hierarchical b-spline grids, CMAME, Vol. 311, pp. 415-437, 2016. (3) C. Kadapa, W. G. Dettmer, D. Perić, A stabilised immersed boundary method on hierarchical b-spline grids for fluid-rigid body interaction with solid-solid contact, CMAME, Vol. 318, pp. 242-269, 2017. (4) Y. Bazilevs, K. Takizawa, T. E. Tezduyar, Computational Fluid-Structure Interaction: Methods and Applications, Wiley, 2013.