Power/Speed Prediction of Ice-Going Ships Based on DEM-SPH Fluid-Structure Interaction Model

doi:10.19423/j.cnki.31-1561/u.2024.01.007

Abstract

Abstract: The navigation of ships in sea ice is affected by the resistance of the ice-water environment, which is an important factor in the research of the power/speed of ice ships. The ship sailing in ice areas has been simulated by a fluid-structure interaction model based on the discrete element method (DEM) and smoothed particle hydrodynamics (SPH) method, to reasonably analyze the power/speed performance of ice ships. It obtains the resistance of the ship at different speeds, and the thrust and torque of the propeller, as well as the speed of revolution required for navigation at constant speed. A DEM-SPH model for the coupling of the hull structure, sea ice and seawater is established by directly calculating the interaction force between the structure and the fluid through the fitting of the relative motion between the SPH particles and the fixed particle boundaries, in order to study the fluid structure interaction between the hull structure, sea ice and seawater. Then the DEM-SPH coupling model is used to simulate the resistance and propulsion of the ships in ice regions. The resistance and thrust of the ships are matched through the fitting method with consideration of the influence of the resistance increment caused by the wake field, to predict the speed of revolution required for the self-propelled ship at a specific speed. The resistances of the DTMB 5415 ship model in ice floes and level ice and the propeller thrusts at different speeds of revolution are also simulated to predict the speed of revolution required for the ship sailing at specific speeds under different operation conditions. The results show that the fluid-structure interaction between the ship-ice and propeller-ice can be well simulated by the DEM-SPH coupling model, fully describing the drag effect of the ship hull and the wake field on the sea ice. A numerical simulation-based prediction of the speed of revolution for the self-propelled ship is performed through the simulations of resistance and thrust, and the fitting analysis of the forced force. It can provide foundation for further experimental verification and engineering applications.

Key words: ice-going ships, DEM-SPH fluid-structure interaction model, discrete element method, ship resistance, propeller thrust

摘要： 船舶在冰区海域中航行会受到冰水环境阻力的影响,是冰区船舶快速性研究中的重要影响因素。为合理分析冰区船舶的快速性能,该文采用基于离散元（discrete element method, DEM）和光滑粒子流体动力学（smoothed particle hydrodynamics, SPH）方法的流固耦合模型模拟船舶冰区航行过程,获得不同航速下的船舶阻力和推进力,进一步计算出螺旋桨的推力、扭矩以及定速航行所需的螺旋桨转速等参数。为研究船体结构、海冰与海水之间的流固耦合作用,文中通过SPH粒子与固定粒子边界相对运动的拟合项直接计算固体与流体之间的相互作用力,建立船体结构、海冰与海水耦合的DEM-SPH模型,并基于该模型分别对船舶在冰区的航行阻力和推进力进行模拟,通过拟合的方式匹配航行阻力和推进力,并考虑尾部流场导致的船体阻力增额,从而预报船舶在特定航速下实现自航所需的螺旋桨转速。此外,文中还模拟了DTMB 5415船模在浮冰区和层冰区中航行的阻力和不同螺旋桨转速下的推力,对船模在不同工况下实现特定航速航行所需的螺旋桨转速进行了预报。计算结果表明：DEM-SPH耦合模型对船-冰、桨-冰作用中的流固耦合过程模拟效果出色,可完整描述船体及尾部伴流场对海冰的拖曳作用;通过文中所述阻力-推力模拟算例及强制力的拟合分析,所形成的基于数值模拟方法的船舶自航下螺旋桨转速预报,可为进一步的试验验证和工程应用推广奠定基础。

关键词: 冰区船舶, DEM-SPH流固耦合, 离散元方法, 船体阻力, 螺旋桨推力

CLC Number:

WU Jie, LIU Lu, TIAN Yukui, JI Shunying. Power/Speed Prediction of Ice-Going Ships Based on DEM-SPH Fluid-Structure Interaction Model[J]. Ship & Boat, 2024, 35(01): 70-83.

吴捷, 刘璐, 田于逵, 季顺迎. 基于DEM-SPH流固耦合的冰区船舶快速性预报[J]. 船舶, 2024, 35(01): 70-83.

References

[1] 黄焱, 孙剑桥, 田育丰. 极地船舶冰区航行性能的试验预报技术发展现状[J]. 船舶, 2023, 34(1): 87-97.
[2] CARRICA P M, CASTRO A M, STERN F.Self-propulsion computations using a speed controller and a discretized propeller with dynamic overset grids[J]. Journal of Marine Science and Technology, 2010, 15(4): 316-330.
[3] XIE C, ZHOU L, DING S, et al.Experimental and numerical investigation on self-propulsion performance of polar merchant ship in brash ice channel[J]. Ocean Engineering, 2023, 269: 113424.
[4] YANG B, SUN Z, ZHANG G, et al.Numerical estimation of ship resistance in broken ice and investigation on the effect of floe geometry[J]. Marine Structures, 2021, 75: 102867.
[5] YANG D, LIU L, JI S.Numerical analysis of interaction between sea ice and propeller based on coupled DEM-FEM model[J]. Ocean Engineering, 2023, 268: 113469.
[6] LU W, LUBBAD R, LØSET S. Simulating ice-sloping structure interactions with the cohesive element method[J]. Journal of Offshore Mechanics & Arctic Engineering, 2014, 136(3): 16.
[7] PRADANA M R, QIAN X.Bridging local parameters with global mechanical properties in bonded discrete elements for ice load prediction on conical structures[J]. Cold Regions Science and Technology, 2020, 173: 102960.
[8] ZHANG Y, ZHANG G, TAO L, et al.Study and discussion on computational efficiency of ice-structure interaction by peridynamic[J]. Journal of Marine Science and Engineering, 2023, 11(6): 1154.
[9] ZHAO X, SHEN H H, CHENG S.Modeling ocean wave propagation under sea ice covers[J]. Acta Mechanica Sinica, 2015, 31(1): 1-15.
[10] 刘璐, 尹振宇, 季顺迎. 船舶与海洋平台结构冰荷载的高性能扩展多面体离散元方法[J]. 力学学报, 2019, 51(6): 1720-1739.
[11] 季顺迎, 田于逵. 基于多介质、多尺度离散元方法的冰载荷数值冰水池[J]. 力学学报, 2021, 53(9): 2427-2453.
[12] JI S, YANG D.Ice loads and ice-induced vibrations of offshore wind turbine based on coupled DEM-FEM simulations[J]. Ocean Engineering, 2022, 243: 110197.
[13] HUANG L F, REN K, LI M H, et al.Fluid-structure interaction of a large ice sheet in waves[J]. Ocean Engineering, 2019, 182: 102-111.
[14] LUO W Z, JIANG D P, WU T C, et al.Numerical simulation of an ice-strengthened bulk carrier in brash ice channel[J]. Ocean Engineering, 2020, 196: 106830.
[15] 周光正, 葛蔚, 李静海. 传统光滑粒子动力学方法的适用性分析[J]. 科学通报, 2013, 58(15): 1414-1421.
[16] 王志超, 李大鸣. 基于SPH-DEM流-固耦合算法的滑坡涌浪模拟[J]. 岩土力学, 2017, 38(4): 1226-1232.
[17] TANG Y, JIANG Q, ZHOU C.A Lagrangian-based SPH-DEM model for fluid-solid interaction with free surface flow in two dimensions[J]. Applied Mathematical Modelling, 2018, 62: 436-460.
[18] ANTUONO M, COLAGROSSI A, MARRONEA S, et al.Propagation of gravity waves through an SPH scheme with numerical diffusive terms[J]. Computer Physics Communications, 2011, 182: 866-877.
[19] SUN X, SAKAI M, YAMADA Y.Three-dimensional simulation of a solid-liquid flow by the DEM-SPH method[J]. Journal of Computational Physics, 2013, 248(5): 147-176.
[20] EITZLMAYR A, KOSCHER G, KHINAST J.A novel method for modeling of complex wall geometries in smoothed particle hydrodynamics[J]. Computer Physics Communications, 2014, 185(10): 2436-2448.
[21] 刘璐, 曹晶, 张志刚, 等. 冰区航行中船体结构冰压力分布特性的离散元分析[J]. 船舶力学, 2021, 25(4): 453-461.
[22] POTYONDY D O, CUNDALL P A.A bonded-particle model for rock[J]. International Journal of Rock Mechanics & Mining Sciences, 2004, 41(8): 1329-1364.
[23] LONG X, JI S Y, WANG Y F.Validation of micro-parameters in discrete element modeling of sea ice failure process[J]. Particulate Science and Technology, 2019, 37(5): 546-555.
[24] LONG X, LIU S W, JI S Y.Discrete element modelling of relationship between ice breaking length and ice load on conical structure[J]. Ocean Engineering, 2020, 201: 107152.
[25] NOMERITAE D E, GRIMALDI S, et al.Explicit incompressible SPH algorithm for free-surface flow modelling: a comparison with weakly compressible schemes[J]. Advances in Water Resources, 2016, 97: 156-167.
[26] MORIKAWA D, SENADHEERA H, ASAI M.Explicit incompressible smoothed particle hydrodynamics in a multi-GPU environment for large-scale simulations[J]. Computational Particle Mechanics, 2021, 8(3): 493-510.
[27] LIU G R, LIU M B.Smoothed particle hydrodynamics: a meshfree particle method[M]. Singapore: World Scientific Publishing, 2003.
[28] CLEARY P W, MONAGHAN J J.Conduction modelling using smoothed particle hydrodynamics[J]. Journal of Computational Physics, 1999, 148(1): 227-264.
[29] SHAO S, LO E Y.Incompressible SPH method for simulating newtonian and non-newtonian flows with a free surface[J]. Advances in Water Resources, 2003, 26(7): 787-800.
[30] RAFIEE A, THIAGARAJAN K P.An SPH projection method for simulating fluid-hypoelastic structure interaction[J]. Computer Methods in Applied Mechanics & Engineering, 2009, 198: 2785-2795.
[31] NI X, SHENG J, FENG W.Simulation of free-surface flow using the smoothed particle hydrodynamics (SPH) method with radiation open boundary conditions[J]. Journal of Atmospheric and Oceanic Technology, 2016, 33(11): 2435-2460.
[32] VERBRUGGHE T, DOMÍNGUEZ J M, ALTOMARE C, et al. Non-linear wave generation and absorption using open boundaries within DualSPHysics[J]. Computer Physics Communications, 2019, 240: 46-59.
[33] HOLMES D W, PIVONKA P.Novel pressure inlet and outlet boundary conditions for smoothed particle hydrodynamics, applied to real problems in porous media flow[J]. Journal of Computational Physics, 2021, 429: 110029.
[34] EITZLMAYR A, KHINAST J.Co-rotating twin-screw extruders: detailed analysis of conveying elements based on smoothed particle hydrodynamics. part 1: hydrodynamics[J]. Chemical Engineering Science, 2015, 134: 861-879.
[35] LIU L, WU J, JI S Y.DEM-SPH coupling method for the interaction between irregularly shaped granular materials and fluids[J]. Powder Technology, 2022, 400: 117249.