MPS与GPU结合数值模拟LNG液舱晃荡

陈翔; 万德成

doi:10.6052/0459-1879-18-410

MPS与GPU结合数值模拟LNG液舱晃荡

doi: 10.6052/0459-1879-18-410

陈翔,
万德成^2), ,

上海交通大学船舶海洋与建筑工程学院, 海洋工程国家重点实验室, 高新船舶与深海开发装备协同创新中心,上海 200240

基金项目: 1)国家自然科学基金 (51879159, 51490675, 11432009, 51579145),长江学者奖励计划(T2014099),上海高校特聘教授(东方学者)岗位跟踪计划(2013022),上海市优秀学术带头人计划(17XD1402300),工信部数值水池创新专项课题(2016-23/09)资助项目.

详细信息

通讯作者:
万德成

中图分类号: O35,U633;
计量
- 文章访问数: 979
- HTML全文浏览量: 116
- PDF下载量: 195
- 被引次数: 0
出版历程
- 收稿日期: 2018-12-05
- 刊出日期: 2019-05-18

NUMERICAL SIMULATION OF LIQUID SLOSHING IN LNG TANK USING GPU-ACCELERATED MPS METHOD

Xiang Chen,
Decheng Wan^{2)
, ,}

Collaborative Innovation Center for Advanced Ship and Deep-Sea Exploration, State Key Laboratory of Ocean Engineering, School of Naval Architecture, Ocean and Civil Engineering, Shanghai Jiao Tong University, Shanghai 200240, China

摘要

摘要: 液舱晃荡是一种在外部激励作用下部分装载的液舱内液体的波动现象,它会对液舱结构强度和运输船舶稳性产生危害.移动粒子半隐式法(moving particle semi-implicit,MPS)是一种典型的无网格粒子类方法,可以有效地模拟剧烈的液舱晃荡问题.但MPS方法存在计算效率低的缺点,难以模拟大规模三维问题,而GPU并行加速技术已广泛应用于科学计算领域.因此,本文将MPS方法与GPU并行加速技术相结合,采用CUDA程序语言编写,自主开发了MPSGPU-SJTU求解器,对三维液化天然气(liquefiednatural gas, LNG)型液舱晃荡进行了数值模拟.通过三种不同粒子间距的数值模拟,验证了求解器的收敛性,其中最大计算粒子数达到了200多万.与其他研究结果相比,MPSGPU-SJTU求解器能够准确地预测壁面砰击压力,并且捕捉晃荡过程中自由面的大幅度变形和强非线性破碎现象.相比CPU求解器的计算时间,GPU并行加速技术可以大幅度地减小计算时长,提高MPS方法的计算效率.本文将LNG型液舱与方型液舱的晃荡进行对比,结果表明在高充液率下LNG型液舱可以有效地减小晃荡幅值和壁面砰击压力.但在中低充液率下,LNG型液舱则会加剧晃荡,自由面呈现明显的三维特征.本文还进一步研究了水和LNG两种不同介质的液舱晃荡现象,数值模拟结果表明二者的流场基本相似,砰击压力则正比于液体密度.
- 无网格粒子法 /
- 移动粒子半隐式法 /
- 并行加速技术 /
- 液舱晃荡
Abstract: Liquid sloshing is a common phenomenon induced in partially filled tanks under external excitations, which may destroy the tank structure and vessel stability. Moving particle semi-implicit (MPS) method is a typical meshfree method which can effectively simulate violent liquid sloshing problem. However, the low computational efficiency of MPS is the bottleneck of its application in large-scale three-dimensional problems. In the past years, GPU parallel acceleration technique has been widely used in the field of scientific computing. In this work, GPU parallel acceleration technique is introduced into MPS method and an in-house solver MPSGPU-SJTU is developed by using CUDA language. Then this solver is used to simulate 3-D liquid sloshing in liquefied natural gas (LNG) tank. The convergent validation of particle spacing is carried out to verify the accuracy of present solver. The maximum particle number of simulation model is over two million particles. MPSGPU-SJTU solver can accurately predict the impact pressures by comparing with other results. In addition, the violent flow phenomena such as large deformation and nonlinear fragmentation of free surface can be observed in these simulations. The comparison of computation time between GPU and CPU solvers demonstrates GPU parallel acceleration technique can significantly reduce the computation time and improve the computational efficiency of MPS. The phenomena of liquid sloshing in LNG tank and rectangular tank are compared. The results show that LNG tank can reduce the sloshing amplitude and impact pressure in high filling level. However, the sloshing is more violent and the free surface presents three-dimensional feature in LNG tank with middle and low filling level. Finally, the investigation of the effect of different fluids such as water and LNG on sloshing phenomena is also conducted in this paper. It shows that the flow fields of both liquids are almost similar and the impact pressure is proportional to the liquid density.
- meshfree particle method /
- MPS method /
- parallel acceleration technique /
- liquid sloshing

HTML全文

参考文献(1)

[1]	Faltinsen OM.A numerical nonlinear method of sloshing in tanks with two dimensional flow. Journal of Ship Research, 1978, 22(3): 193-202
[2]	Nakayama T, Washizu K.Nonlinear analysis of liquid motion in a container subjected to forced pitching oscillation. International Journal for Numerical Methods in Engineering, 1980, 15(8): 1207-1220
[3]	Nakayama T, Washizu K.The boundary element method applied to the analysis of two-dimensional nonlinear sloshing problems. International Journal for Numerical Methods in Engineering, 1981, 17(11): 1631-1646
[4]	Kim Y, Shin YS, Lee KH.Numerical study on slosh-induced impact pressures on three-dimensional prismatic tanks. Applied Ocean Research, 2004, 26(5): 213-226
[5]	Zhu RQ, Fang ZY, Zhang ZG, et al.Level-set method for predicting impact pressure induced by violent sloshing in a tank. Journal of Ship Mechanics, 2008, 12(3): 344-351
[6]	Xue MA, Lin PZ.Numerical study of ring baffle effects on reducing violent liquid sloshing. Computers & Fluids, 2011, 52: 116-129
[7]	朱跃, 姜胜耀, 杨星团等. 粒子法中压力振荡的机理研究. 力学学报, 2018, 50(3): 688-698
[7]	(Zhu Yue, Jiang Shengyao, Yang Xingtuan, et al.Mechanism analysis of pressure oscillation in particle method. Chinese Journal of Theoretical and Applied Mechanics, 2018, 50(3): 688-698 (in Chinese))
[8]	马文涛. 二维弹性力学问题的光滑无网格伽辽金法.力学学报, 2018, 50(5): 1115-1124
[8]	(Ma Wentao.A smoothed meshfree Galerkin method for 2D elasticity problem. Chinese Journal of Theoretical and Applied Mechanics, 2018, 50(5): 1115-1124 (in Chinese))
[9]	李艾伦, 傅卓佳, 李柏纬等. 含肿瘤皮肤组织传热分析的广义有限差分法. 力学学报, 2018, 50(5): 1198-1205
[9]	(Li Ailun, Fu Zhuojia, Li Powei, et al.Generalized finite difference method for bioheat transfer analysis on skin tissue with tumors. Chinese Journal of Theoretical and Applied Mechanics, 2018, 50(5): 1198-1205 (in Chinese))
[10]	Shao JR, Li HQ, Liu GR, et al. An improved SPH method for modeling liquid sloshing dynamics. Computers & Structures, 2012, 100-101: 18-26
[11]	Shao JR, Li SM, Li ZR, et al.A comparative study of different baffles on mitigating liquid sloshing in a rectangular tank due to a horizontal excitation. Engineering Computations, 2015, 32(4): 1172-1190
[12]	Zhang ZL, Liu MB.A decoupled finite particle method for modeling incompressible flows with free surfaces. Applied Mathematical Modelling, 2018, 60: 606-633
[13]	Rueda GE, Tsukamoto MM, Medeiros HF, et al. Validation study of MPS (moving particle semi-implicit method) for sloshing and damage stability analysis//Proceedings of the ASME 2008 27th International Conference on Offshore Mechanics and Arctic Engineering, Estoril, Portugal, 2008-6-15-20. 2008 ASME: 483-489
[14]	潘徐杰, 张怀新. 移动粒子半隐式法晃荡模拟中的压力震荡现象研究. 水动力学研究与进展, 2008, 23(4): 453-463
[14]	(Pan Xujie, Zhang Huaixin.A study on the oscillations appearing in pressure calculation for sloshing simulation by using moving-particle semi-implicit method. Chinese Journal of Hydrodynamics, 2008, 23(4): 453-463 (in Chinese))
[15]	Lee BH, Park JC, Kim MH. Two-dimensional vessel-motion/liquid-sloshing interactions and impact loads by using a particle method//Proceedings of the ASME 2010 29th International Conference on Ocean, Offshore and Arctic Engineering, Shanghai, China, 2010-6-6-11. 2010 ASME: 679-686
[16]	Tsukamoto MM, Cheng LY, Nishimoto K.Analytical and numerical study of the effects of an elastically-linked body on sloshing. Computers & Fluids, 2011, 49(1): 1-21
[17]	Zhang YX, Wan DC.Comparative study of MPS method and level-set method for sloshing flows. Journal of Hydrodynamics, 2014, 26(4): 577-585
[18]	杨亚强. 基于MPS方法的液舱晃荡数值模拟与分析. [硕士论文]. 上海: 上海交通大学, 2016
[18]	(Yang Yaqiang.Numerical investigation of liquid sloshing by MPS method. [Master Thesis]. Shanghai: Shanghai Jiao Tong University, 2016 (in Chinese))
[19]	Zhang YL, Wan DC.MPS-FEM coupled method for sloshing flows in an elastic tank. Ocean Engineering, 2018, 152: 416-427
[20]	Zhang YL, Chen X, Wan DC.MPS-FEM coupled method for the comparison study of liquid sloshing flows interacting with rigid and elastic baffles. Applied Mathematics and Mechanics, 2016, 37(12): 1359-1377
[21]	Wen X, Wan DC. Numerical simulation of three-layer-liquid sloshing by multiphase mps method//Proceedings of the ASME 2018 37th International Conference on Ocean, Offshore and Arctic Engineering, Madrid, Spain, 2018-6-17-22. 2018 ASME, OMAE2018-78387
[22]	Wen X, Wan DC, Chen G. Multiphase MPS method for two-layer-liquid sloshing flows in oil-water separators//Proceedings of the Twenty-eighth (2018) International Ocean and Polar Engineering Conference, Sapporo, Japan, 2018-6-10-15. 2018 by the International Society of Offshore and Polar Engineers (ISOPE): 859-866
[23]	Hori C, Gotoh H, Ikari H, et al.GPU-acceleration for moving particle semi-implicit method. Computers & Fluids, 2011, 51: 174-183
[24]	Zhu XS, Cheng L, Lu L, et al.Implementation of the moving particle semi-implicit method on GPU. Science China Physics, Mechanics and Astronomy, 2011, 54: 523-532
[25]	Kakuda K, Nagashima T, Hayashi Y, et al.Particle-based fluid flow simulations on GPGPU using CUDA. Computer Modeling in Engineering & Sciences, 2012, 88: 17-28
[26]	Kakuda K, Nagashima T, Hayashi Y, et al.Three dimensional fluid flow simulations using GPU-based particle method. Computer Modeling in Engineering & Sciences, 2013, 93: 363-376
[27]	李海洲, 唐振远, 万德成. 三维自由面流动模拟中GPU并行计算技术. 海洋工程, 2016, 34(5): 20-29
[27]	(Li Haizhou. Tang Zhenyuan, Wan Decheng.Application of GPU acceleration techniques in 3D violent flow. The Ocean Engineering, 2016, 34(5): 20-29 (in Chinese))
[28]	Chen X, Wan DC. Numerical simulation of three-dimensional violent free surface flows by GPU-based MPS method. International Journal of Computational Methods, 2018, 15(8): 1843012-1-20
[29]	Koshizuka S, Oka Y.Moving-particle semi-implicit method for fragmentation of incompressible fluid. Nuclear Science and Engineering, 1996, 123(3): 421-434
[30]	Zhang YX, Wan DC.Numerical simulation of liquid sloshing in low-filling tank by MPS. Journal of Hydrodynamics, 2012, 27: 101-107
[31]	Tanaka M, Masunaga T.Stabilization and smoothing of pressure in MPS method by quasi-compressibility. Journal of Computational Physics, 2010, 229(11): 4279-4290
[32]	Lee BH, Park JC, Kim MH, et al.Step-by-step improvement of MPS method in simulating violent free-surface motions and impact loads. Computer Methods in Applied Mechanics and Engineering, 2011, 200(9-12): 1113-1125
[33]	NIVIDIA. CUDA Toolkit Documentation v9.2.88. , 2018
[34]	张雨新, 万德成. MPS方法在三维溃坝问题中的应用. 中国科学: 物理学力学天文学, 2011, 41(2): 140-154
[34]	(Zhang Yuxin, Wan Decheng.Application of MPS in 3D dam breaking flows. Scientia Sinica Phys, Mech & Astron, 2011, 41(2): 140-154 (in Chinese))
[35]	Tang ZY, Wan DC, Chen G, et al.Numerical simulation of 3D violent free-surface flows by multi-resolution MPS method. Journal of Ocean Engineering and Marine Energy, 2016, 2: 355-364
[36]	Rao CP, Wan DC.Numerical study of the wave-induced slamming force on the elastic plate based on MPS-FEM coupled method. Journal of Hydrodynamics, 2018, 30(1): 70-78
[37]	Wen X, Wan DC. Numerical simulation of rayleigh--taylor instability by multiphase MPS method. International Journal of Computational Methods, 2018, 15(3): 1846005-1-1846005-12
[38]	蔡忠华. 液货船液舱晃荡问题研究. [博士论文]. 上海: 上海交通大学, 2012
[38]	(Cai Zhonghua.Study on the sloshing problems of liquid cargo tanks. [PhD Thesis]. Shanghai: Shanghai Jiao Tong University, 2012 (in Chinese))
[39]	陆志妹, 范佘明, 朱仁传 .基于CFD的不同介质液舱晃荡比较分析//第二十五届全国水动力学研讨会暨第十二届全国水动力学学术会议, 浙江, 舟山, 2013-9-13: 665-672
[39]	(Lu Zhimei, Fan Sheming, Zhu Renchuan. Comparative analysis of sloshing with different fluid materials based CFD//Proceedings of the 25th National Conference on Hydrodynamics & 12th National Congress on Hydrodynamics, Zhou Shan, Zhejiang, 2013-9-13: 665-672 (in Chinese))