A sharp interface method for compressible liquid-vapor flow with phase transition and surface tension
A sharp interface method for compressible liquid-vapor flow with phase transition and surface tension
复制标题
DOI：
10.1016/j.jcp.2017.02.001
复制
发表时间：
2015-10
期刊：
影响因子：
--
通讯作者：
S. Fechter;C. Munz;C. Rohde;C. Zeiler
中科院分区：
文献类型：
--
作者：
S. Fechter;C. Munz;C. Rohde;C. Zeiler
在可压缩区域内的非等温液-汽流动的数值近似是一项困难的任务，因为相界面处的复杂物理效应可以控制全局流动行为。我们提出了一个尖锐的接口的方法，将接口作为一个冲击波一样的不连续性。通过仅在体区域中使用流动求解器和接近界面的虚流体方法，避免了流体相的任何混合。体区域数值解的耦合状态由界面上局部两相黎曼问题的解确定。黎曼解通过在相边界处施加适当的跳跃条件来解释相关的物理现象。各种各样的界面效果都可以用一种逻辑上一致的方式来处理。这包括通过相变的表面张力或质量/能量转移。此外，界面的局部法向速度，这是需要计算的界面的时间演化，由黎曼解。界面跟踪本身是基于水平集方法的，本文的重点是描述两阶段黎曼求解器及其在锐界面方法中的使用。一维问题的选择，以验证该方法。最后，一个摆动液滴和冲击液滴在两个维度的相互作用的三维模拟。在这两个问题中，相变和表面张力决定了整体行为。
The numerical approximation of non-isothermal liquid–vapor flow within the compressible regime is a difficult task because complex physical effects at the phase interfaces can govern the global flow behavior. We present a sharp interface approach which treats the interface as a shock-wave like discontinuity. Any mixing of fluid phases is avoided by using the flow solver in the bulk regions only, and a ghost-fluid approach close to the interface. The coupling states for the numerical solution in the bulk regions are determined by the solution of local two-phase Riemann problems across the interface. The Riemann solution accounts for the relevant physics by enforcing appropriate jump conditions at the phase boundary. A wide variety of interface effects can be handled in a thermodynamically consistent way. This includes surface tension or mass/energy transfer by phase transition. Moreover, the local normal speed of the interface, which is needed to calculate the time evolution of the interface, is given by the Riemann solution. The interface tracking itself is based on a level-set method.The focus in this paper is the description of the two-phase Riemann solver and its usage within the sharp interface approach. One-dimensional problems are selected to validate the approach. Finally, the three-dimensional simulation of a wobbling droplet and a shock droplet interaction in two dimensions are shown. In both problems phase transition and surface tension determine the global bulk behavior.