Investigating the Performance of Twin Marine Propellers in Different Ship Wake Fields Using an Unsteady Viscous and Inviscid Solver

Abstract

In this study, the performance of a twin-screw propeller under the influence of the wake field of a fully appended ship was investigated using a coupled Reynolds-averaged Navier—Stokes (RANS)/boundary element method (BEM) code. The unsteady BEM is an efficient approach to predicting propeller performance. By applying the time-stepping method in the BEM solver, the trailing vortex sheet pattern of the propeller can be accurately captured at each time step. This is the main innovation of the coupled strategy. Furthermore, to ascertain the effect of the wake field of the ship with acceptable accuracy, a RANS solver was developed. A finite volume method was used to discretize the Navier—Stokes equations on fully unstructured grids. To simulate ship motions, the volume of the fluid method was applied to the RANS solver. The validation of each solver (BEM/RANS) was separately performed, and the results were compared with experimental data. Ultimately, the BEM and RANS solvers were coupled to estimate the performance of a twin-screw propeller, which was affected by the wake field of the fully appended hull. The proposed model was applied to a twin-screw oceanography research vessel. The results demonstrated that the presented model can estimate the thrust coefficient of a propeller with good accuracy as compared to an experimental self-propulsion test. The wake sheet pattern of the propeller in open water (uniform flow) was also compared with the propeller in a real wake field.

Keywords

Twin propeller / Reynolds-averaged Navier—Stokes (RANS) / Boundary element method (BEM) / Time-stepping method (TSM) / Wake sheet pattern / Effective wake field

References

[1] Benek J, Steger J, Dougherty FC (1983) A flexible grid embedding technique with application to the Euler equations. 6th Computational Fluid Dynamics Conference, Danvers

[2] Carlton J (2018) Marine propellers and propulsion. Butterworth-Heinemann

[3] Durante D, Dubbioso G, Testa C. Simplified hydrodynamic models for the analysis of marine propellers in a wake-field. Journal of Hydrodynamics, 2013, 25(6): 954-965

[4] Gaggero S, Dubbioso G, Villa D, Muscari R, Viviani M. Propeller modeling approaches for off—design operative conditions. Ocean Engineering, 2019, 178: 283-305

[5] Gaggero S, Villa D, Viviani M (2014) An investigation on the discrepancies between RANSE and BEM approaches for the prediction of marine propeller unsteady performances in strongly non-homogeneous wakes. International Conference on Offshore Mechanics and Arctic Engineering, American Society of Mechanical Engineers

[6] Gaggero S, Villa D, Viviani M. An extensive analysis of numerical ship self-propulsion prediction via a coupled BEM/RANS approach. Applied Ocean Research, 2017, 66: 55-78

[7] Gaskell P, Lau A. Curvature-compensated convective transport: SMART, a new boundedness-preserving transport algorithm. International Journal for Numerical Methods in Fluids, 1988, 8(6): 617-641

[8] Greve M, Wöckner-Kluwe K, Abdel-Maksoud M, Rung T. Viscous-inviscid coupling methods for advanced marine propeller applications. International Journal of Rotating Machinery, 2012, 2012: 743060

[9] Henshaw WD, Schwendeman DW. Moving overlapping grids with adaptive mesh refinement for high-speed reactive and non-reactive flow. Journal of Computational Physics, 2006, 216(2): 744-779

[10] Hess J, Smith A. Calculation of non-lifting potential flow about arbitrary three dimensional bodies. J ship Res, 1962, 8(4): 22-44

[11] Hess JL (1972) Calculation of potential flow about arbitrary three-dimensional lifting bodies: Douglas Aircraft Co Long Beach CA

[12] Hess JL, Valarezo WO. Calculation of steady flow about propellers using a surface panel method. Journal of Propulsion and Power, 1985, 1(6): 470-476

[13] ITTC-Recommended Procedures and Guidelines (2011) Performance Prediction Method 7.5-02-03-01.4: 9

[14] Jasak H. Error analysis and estimation for the finite volume method with applications to fluid flows, 1996, London: University of London

[15] Katz J, Plotkin A (2001) Low-speed aerodynamics. Cambridge University Press

[16] Kim D, Choi H. A second-order time-accurate finite volume method for unsteady incompressible flow on hybrid unstructured grids. Journal of Computational Physics, 2000, 162(2): 411-428

[17] Leonard B. The ULTIMATE conservative difference scheme applied to unsteady one-dimensional advection. Computer Methods in Applied Mechanics and Engineering, 1991, 88(1): 17-74

[18] Meakin R, Suhs N (1989) Unsteady aerodynamic simulation of multiple bodies in relative motion. 9th Computational Fluid Dynamics Conference

[19] Najafi S, Abbaspoor M. Numerical investigation of flow pattern and hydrodynamic forces of submerged marine propellers using unsteady boundary element method. Proceedings of the Institution of Mechanical Engineers, Part M: Journal of Engineering for the Maritime Environment, 2019, 233(1): 67-79

[20] Najafi S, Abbaspour M. Numerical study of propulsion performance in swimming fish using boundary element method. Journal of the Brazilian Society of Mechanical Sciences and Engineering, 2017, 39(2): 443-455

[21] Panahi R, Jahanbakhsh E, Seif MS. Development of a VoF-fractional step solver for floating body motion simulation. Applied Ocean Research, 2006, 28(3): 171-181

[22] Panahi R, Jahanbakhsh E, Seif MS. Towards simulation of 3D nonlinear high-speed vessels motion. Ocean Engineering, 2009, 36(3–4): 256-265

[23] Politis G. Unsteady wake rollup modeling using a mollifier based filtering technique. Dev. Appl. Ocean. Eng, 2016, 5(1): 1-28

[24] Politis GK. Simulation of unsteady motion of a propeller in a fluid including free wake modeling. Engineering Analysis with Boundary Elements, 2004, 28(6): 633-653

[25] Politis GK. Unsteady rollup modeling for wake-adapted propellers using a time-stepping method. Journal of Ship Research, 2005, 49(3): 216-231

[26] Politis GK. Application of a BEM time stepping algorithm in understanding complex unsteady propulsion hydrodynamic phenomena. Ocean Engineering, 2011, 38(4): 699-711

[27] Pourmostafa M, Ghadimi P. Applying boundary element method to simulate a high-skew Controllable Pitch Propeller with different hub diameters for preliminary design purposes. Cogent Engineering, 2020, 7(1): 1805857

[28] Pourmostafa M, Ghadimi P. Unsteady 2D and 3D Navier-Stokes Solver with application of multigrid scheme to pressure poisson fractional step on arbitrary unstructured grids in various applications with emphasis on ship motion. Mathematical Problems in Engineering, 2020, 2020: 1-28

[29] Queutey P, Deng GB, Guilmineau E, Salvatore F (2013) A comparison between full RANSE and coupled RANSE-BEM approaches in ship propulsion performance prediction. International Conference on Offshore Mechanics and Arctic Engineering, OMAE 2013–10366

[30] Rao ZQ, Yang CJ. Numerical prediction of effective wake field for a submarine based on a hybrid approach and an RBF interpolation. Journal of Hydrodynamics, 2017, 29(4): 691-701

[31] Rijpkema D, Starke B, Bosschers J (2013) Numerical simulation of propeller-hull interaction and determination of the effective wake field using a hybrid RANS-BEM approach. 3rd International Symposium on Marine Propulsors, Launceston, Australia

[32] RV50 Test Report (2011) Hydrodynamic model tests. Test Report No. 2512/01, Model No.2512, Vienna Model Basin (www.sva.at): 43

[33] Ubbink O, Issa R. A method for capturing sharp fluid interfaces on arbitrary meshes. Journal of Computational Physics, 1999, 153(1): 26-50

[34] Wang Y, Abdel-Maksoud M, Song B. Convergence of different wake alignment methods in a panel code for steady-state flows. Journal of Marine Science and Technology, 2016, 21(4): 567-578

[35] Zhu Q, Wolfgang M, Yue D, Triantafyllou M. Three-dimensional flow structures and vorticity control in fish-like swimming. Journal of Fluid Mechanics, 2002, 468: 1-28