Your browser does not seem to support JavaScript. As a result, your viewing experience will be diminished, and you have been placed in read-only mode.
Please download a browser that supports JavaScript, or enable it if it's disabled (i.e. NoScript).
reactingMultiphaseEulerFOam可以代替multiphaseEulerFoam
@金石为开 假期出去了，才看到。我没研究过升力模型，不过我感觉你可以找现成的升力模型加到动量方程里，或者加到压力泊松方程里。
@东岳 reactingMultiphaseEulerFoam求解器之前尝试用过，自带的算例全是层流的，没用湍流模型，我添加湍流模型后湍流方程矩阵求解不收敛，尝试多次没成功后，就没用了～如果有哪位用reactingMultiphaseEulerFOam求解器选择LES模型运行成功，请分享下，谢谢～
@zwl 加湍流的时候要防止相场接近于0，可以把初始场设置一个小的相场值。LES那面可以试试smagorinsky，不需要求解湍流变量
@zwl 你好，我粗略的拜读了一下大作，我想问一下drag force里面的当地离散项的直径是怎么计算的，论文说取为50mm，是没有计算吗？
@yfclark 是50微米，根据文献实验数据取值，这个目前我还没见到有做动态直径的，而且这个直径还很少有相应理论计算，如果你见到相关理论，请分享一下～
@zwl 我是才接触雾化仿真，我个人有个不成熟的想法，假设离散的液滴和气泡是球形，用曲率转化得到等效直径，这个想法存在一些问题：离散项并非都是球形，特别是在稠密区的液核附近，出现了液核上脱落液丝，液丝破碎为液滴，这些几何结构并不是球形的，即便是在稀疏区，液滴和气泡也有可能震荡变形，这样得到的直径只是等效直径，算得准不准不好说，另外要识别离散项的位置，如果用液核表面的曲率得到一个等效直径是错位的，怎样识别是离散项还是液核表面也是一个问题
曲率转化似乎也不太对
@东岳 为什么会加湍流时要防止相场接近于0？还请东岳老师明示，这是哪部分数值处理出现的问题？如果物理情况确实是有相场等于零的情况，那这个求解器就不能选择LES/kEqn模型来求解湍流了，对吗？
@zwl 防止奇异问题。多相湍流方程和动量方程都是一样的，都需要处理，在我之前那个CES里面讨论过。
如果物理情况确实是有相场等于零的情况，那这个求解器就不能选择LES/kEqn模型来求解湍流了，对吗？
可以，需要特殊处理下，OpenFOAM已经给你处理好了啊。只不过Smagorinsky不存在这个问题不需要特殊处理
@东岳 那为什么reactingEulerFoam选LES/kEqn还是不能用于有相场为0的条件？您说的of已经给我处理好了是指？我正想去看个究竟，谢谢。
@东岳 @zwl 我有个问题，E-E-VOF混合模型里面表面张力不应该是“微观”的力吗，曳力和虚拟质量力应该是积分后的“宏观”的力，那么曳力和虚拟质量力应该包含了表面张力效应吧，把它们写在同一个方程里面，不就相当于表面张力算了两次吗？我没有想明白，请老师解答一下。
@东岳 我最近用reactingMultiphaseEulerFoam做算例，也是设置的时候某一相为0，湍流模型一直出现divide by zero，找了半天没发现错误，更改之后搞定了，谢谢东岳老师
