A numerical study of three-dimensional liquid sloshing in tanks
A numerical study of three-dimensional liquid sloshing in tanks
复制标题
罐内三维液体晃动的数值研究
DOI：
10.1016/j.jcp.2007.12.006
复制
发表时间：
2008-04-01
影响因子：
4.1
通讯作者：
Lin, Pengzhi
中科院分区：
文献类型：
--
作者：
Liu, Dongming;Lin, Pengzhi
本文建立了一个数值模型NEWTANK (numerical Wave TANK)，用于研究具有破碎自由表面的三维非线性液体晃动。该数值模型求解了两相流的空间平均Navier-Stokes方程，该方程建立在具有任意六自由度运动的非惯性参照系上。采用大涡模拟（LES）方法，利用Smagorinsky亚网格尺度（SGS）闭合模型模拟湍流效应。数值求解采用两步投影法，辅以Bi-CGSTAB技术求解过滤后压力场的压力泊松方程。采用二阶精确流体体积法（VOF）对变形和破碎的自由表面进行跟踪。对矩形槽内二维和三维非线性液体晃动进行了室内实验研究。本文还建立了浪涌和摇摆耦合激励下三维液体晃动的线性解析解。本文首先利用已有的无粘流体和粘性流体二维晃动的解析解和实验数据对数值模型进行了验证。将验证进一步扩展到三维液体晃动。当激励幅值较小时，数值结果与解析解吻合。当激励幅值较大时，当晃动高度非线性时，数值计算结果与解析解存在较大差异，但解析解与实验结果吻合较好。最后，对六自由度激励下具有破碎自由表面的剧烈液体晃动进行了模拟和讨论。(C) 2007爱思唯尔公司版权所有。
A numerical model NEWTANK (Numerical Wave TANK) has been developed to study three-dimensional (3-D) non-linear liquid sloshing with broken free surfaces. The numerical model solves the spatially averaged Navier-Stokes equations, which are constructed on a non-inertial reference frame having arbitrary six degree-of-freedom (DOF) of motions, for two-phase flows. The large-eddy-simulation (LES) approach is adopted to model the turbulence effect by using the Smagorinsky sub-grid scale (SGS) closure model. The two-step projection method is employed in the numerical solutions, aided by the Bi-CGSTAB technique to solve the pressure Poisson equation for the filtered pressure field. The second-order accurate volume-of-fluid (VOF) method is used to track the distorted and broken free surface. Laboratory experiments are conducted for both 2-D and 3-D non-linear liquid sloshing in a rectangular tank. A linear analytical solution of 3-D liquid sloshing under the coupled surge and sway excitation is also developed in this study. The numerical model is first validated against the available analytical solution and experimental data for 2-D liquid sloshing of both inviscid and viscous fluids. The validation is further extended to 3-D liquid sloshing. The numerical results match with the analytical solution when the excitation amplitude is small. When the excitation amplitude is large where sloshing becomes highly non-linear, large discrepancies are developed between the numerical results and the analytical solutions, the former of which, however, agree well with the experimental data. Finally, as a demonstration, a violent liquid sloshing with broken free surfaces under six DOF excitations is simulated and discussed. (C) 2007 Elsevier Inc. All rights reserved.