A compressible multiphase flow model for violent aerated wave impact problems
A compressible multiphase flow model for violent aerated wave impact problems
复制标题
DOI：
10.1098/rspa.2014.0542
复制
发表时间：
2014-12
期刊：
影响因子：
--
通讯作者：
Zhihua Ma;D. Causon;Ling Qian;C. Mingham;H. Gu;P. M. Ferrer
中科院分区：
文献类型：
--
作者：
Zhihua Ma;D. Causon;Ling Qian;C. Mingham;H. Gu;P. M. Ferrer
本文着重于数值模拟的波浪冲击事件下的空气截留和曝气的影响。底层流动模型将离散水波视为具有均匀材料特性的空气和水的可压缩混合物。相应的数学方程是基于多相流模型，该模型建立在质量、动量和能量守恒定律以及气相体积分数平流方程的基础上。采用基于单调上游中心守恒律重构格式的高阶有限体积格式离散控制方程的积分形式。通过HLLC近似黎曼求解器的网格单元面的数值通量估计。采用三阶总变差递减龙格-库塔格式获得时间精确解。该模型为处理空气和水气混合物的可压缩性问题提供了一种有效的方法。几个测试用例已计算使用本方法，包括重力引起的液体活塞，自由下降的水柱在一个封闭的坦克，水-空气激波管，砰的一声平板到仍然纯净和充气的水和暴跌波冲击在垂直的墙壁。所得结果与实验、精确解及其它数值计算结果吻合良好。这表明了目前的方法来解决更一般的波-空气-结构相互作用问题的潜力。
This paper focuses on the numerical modelling of wave impact events under air entrapment and aeration effects. The underlying flow model treats the dispersed water wave as a compressible mixture of air and water with homogeneous material properties. The corresponding mathematical equations are based on a multiphase flow model which builds on the conservation laws of mass, momentum and energy as well as the gas-phase volume fraction advection equation. A high-order finite volume scheme based on monotone upstream-centred schemes for conservation law reconstruction is used to discretize the integral form of the governing equations. The numerical flux across a mesh cell face is estimated by means of the HLLC approximate Riemann solver. A third-order total variation diminishing Runge–Kutta scheme is adopted to obtain a time-accurate solution. The present model provides an effective way to deal with the compressibility of air and water–air mixtures. Several test cases have been calculated using the present approach, including a gravity-induced liquid piston, free drop of a water column in a closed tank, water–air shock tubes, slamming of a flat plate into still pure and aerated water and a plunging wave impact at a vertical wall. The obtained results agree well with experiments, exact solutions and other numerical computations. This demonstrates the potential of the current method to tackle more general wave–air–structure interaction problems.