Monolithic coupling of the pressure and rigid body motion equations in computational marine hydrodynamics

Abstract

In Fluid Structure Interaction (FSI) problems encountered in marine hydrodynamics, the pressure field and the velocity of the rigid body are tightly coupled. This coupling is traditionally resolved in a partitioned manner by solving the rigid body motion equations once per nonlinear correction loop, updating the position of the body and solving the fluid flow equations in the new configuration. The partitioned approach requires a large number of nonlinear iteration loops per time–step. In order to enhance the coupling, a monolithic approach is proposed in Finite Volume (FV) framework, where the pressure equation and the rigid body motion equations are solved in a single linear system. The coupling is resolved by solving the rigid body motion equations once per linear solver iteration of the pressure equation, where updated pressure field is used to calculate new forces acting on the body, and by introducing the updated rigid body boundary velocity in to the pressure equation. In this paper the monolithic coupling is validated on a simple 2D heave decay case. Additionally, the method is compared to the traditional partitioned approach (i.e. “strongly coupled” approach) in terms of computational efficiency and accuracy. The comparison is performed on a seakeeping case in regular head waves, and it shows that the monolithic approach achieves similar accuracy with fewer nonlinear correctors per time–step. Hence, significant savings in computational time can be achieved while retaining the same level of accuracy.

Keywords

monolithic coupling / pressure equation / rigid body motion / computational fluid dynamics / marine hydrodynamics / seakeeping

References

[1] Bna S, Manservisi S, Aulisa E. A multilevel domain decomposition solver for monolithic fluid-structure interaction problems. 11th International Conference of Numerical Analysis and Applied Mathematics 2013, 2013, 1558: 871-874

[2] Castiglione T, Stern F, Bova S, Kandasamy M. Numerical investigation of the seakeeping behavior of a catamaran advancing in regular head waves. Ocean Engineering, 2011, 38: 1806-1822

[3] Demirdžić I, Perić M. Space conservation law in finite volume calculations of fluid flow. International Journal for Numerical Methods in Fluids, 1988, 8(9): 1037-1050

[4] Eken A, Sahin M. A parallel monolithic algorithm for the numerical simulation of large-scale fluid structure interaction problems. International Journal for Numerical Methods in Fluids, 2016, 80(12): 687-714

[5] Farah P, Vuong AT, Wall WA, Popp A. Volumetric coupling approaches for multiphysics simulations on nonmatching meshes. International Journal for Numerical Methods in Engineering, 2016, 108(12): 1550-1576

[6] Hachem E, Feghali S, Codina R, Coupez T. Immersed stress method for fluidstructure interaction using anisotropic mesh adaptation. International Journal for Numerical Methods in Engineering, 2013, 94(9): 805-825

[7] Heil M, Hazel AL, Boyle J. Solvers for large-displacement fluid-structure interaction problems: segregated versus monolithic approaches. Computational Mechanics, 2008, 43(1): 91-101

[8] Hu Z, Tang WY, Xue HX, Zhang XY. A simple-based monolithic implicit method for strong-coupled fluid-structure interaction problems with free surfaces. Computer Methods in Applied Mechanics and Engineering, 2016, 299: 90-115

[9] Irons BM, Tuck RC. A version of the Aitken accelerator for computer iteration. International Journal for Numerical Methods in Engineering, 1969, 1: 275-277

[10] Jasak H. Error analysis and estimation for the finite volume method with applications to fluid flows, 1996, London: Imperial College of Science, Technology & Medicine

[11] Jasak H, Vukčević V, Gatin I. Numerical simulation of wave loads on static offshore structures. CFD for Wind and Tidal Offshore Turbines, 2015, 95-105

[12] Jog CS, Pal RK. A monolithic strategy for fluid-structure interaction problems. International Journal for Numerical Methods in Engineering, 2011, 85(4): 429-460

[13] Langer U, Yang HD. Robust and efficient monolithic fluid-structure-interaction solvers. International Journal for Numerical Methods in Engineering, 2016, 108(4): 303-325

[14] Larsson L, Stern F, Visonneau M, Hirata N, Hino T, Kim J. Tokyo 2015: A Workshop on CFD in Ship Hydrodynamics, 2015, Tokyo, Japan: NMRI (National Maritime Research Institute)

[15] Legay A, Zilian A, Janssen C. A rheological interface model and its space-time finite element formulation for fluid-structure interaction. International Journal for Numerical Methods in Engineering, 2011, 86(6): 667-687

[16] Miyata H, Orihara H, Sato Y. Nonlinear ship waves and computational fluid dynamics. Proceedings of the Japan Academy Series B—Physical and Biological Sciences, 2014, 90: 278-300

[17] Orihara H, Miyata H. Evaluation of added resistance in regular incident waves by computational fluid dynamics motion simulation using an overlapping grid system. Journal of Marine Science and Technology, 2003, 8: 47-60

[18] Press WH, Teukolsky SA, Vetterling WT, Flannery BP. Numerical Recipes in C++: The Art of Scientific Computing, 2002, Cambridge: Cambridge University Press

[19] Robinson-Mosher A, Schroeder C, Fedkiw R. A symmetric positive definite formulation for monolithic fluid structure interaction. Journal of Computational Physics, 2011, 230(4): 1547-1566

[20] Saad Y. Iterative methods for sparse linear systems, 2003

[21] Simonsen CD, Otzen JF, Joncquez S, Stern F. EFD and CFD for KCS heaving and pitching in regular head waves. Journal of Marine Science and Technology, 2013, 18: 435-459

[22] Tezdogan T, Demirel YK, Kellett P, Khorasanchi M, Incecik A, Turan O. Full-scale unsteady RANS CFD simulations of ship behaviour and performance in head seas due to slow steaming. Ocean Engineering, 2015, 97: 186-206

[23] Vukčević V. Numerical modelling of coupled potential and viscous flow for marine applications, 2016, University of Zagreb: Faculty of Mechanical Engineering and Naval Architecture

[24] Vukčević V, Jasak H. Seakeeping validation and verification using decomposition model based on embedded free surface method. Tokyo 2015: A Workshop on CFD in Ship Hydrodynamics, 2015

[25] Vukčević V, Jasak H. Validation and verification of decomposition model based on embedded free surface method for oblique wave seakeeping simulations. Tokyo 2015: A Workshop on CFD in Ship Hydrodynamics, 2015

[26] Vukčević V, Jasak H, Gatin I, Malenica S. Seakeeping sensitivity studies using the decomposition CFD model based on the ghost fluid method. Proceedings of the 31st Symposium on Naval Hydrodynamics, 2016

[27] Vukčević V, Jasak H, Malenica S. Decomposition model for naval hydrodynamic applications, Part I: Computational method. Ocean Eng., 2016, 121: 37-46

[28] Vukčević V, Jasak H, Malenica S. Decomposition model for naval hydrodynamic applications, Part II: Verification and validation. Ocean Eng., 2016, 121: 76-88

[29] Wu CS, Zhou DC, Gao L, Miao QM. CFD computation of ship motions and added resistance for a high speed trimaran in regular head waves. International Journal of Naval Architecture and Ocean Engineering, 2011, 3: 105-110

[30] Yang P, Xiang J, Fang F, Pavlidis D, Latham JP, Pain CC. Modelling of fluid-structure interaction with multiphase viscous flows using an immersed-body method. Journal of Computational Physics, 2016, 321: 571-592