drift-flux model中相对速度计算问题
-
fluent中提供了不同的滑移速度的计算方法，我想尝试将其移植到OpenFOAM中去，在处理滑移速度的时候遇到了问题：相对速度计算时需要用到相对雷诺数，相对雷诺数计算时需要用相对速度，有点套娃的感觉。这里在计算时需要使用迭代计算的方法么，有没有资料可以参考呢~ 先行谢过各位大佬~
-
谢谢李老师，我先试试~
-
在将这个方程植入到OpenFOAM中之后，出现了明显的相对速度的震荡，惯性加速度最大能到10^5m/s2。植入时直接使用- fvc::ddt(U) - (U & fvc::grad(U))进行计算，使用的时1D的气液激波管问题做的验证，此时无滑移时的速度为20m/s左右，
给定滑移速度时计算的速度场没有震荡，使用的是OpenFOAM中的general方法。
想了解下出现这种问题的可能原因是什么呢~ -
把源代码放上来，我用的是v2312版本的OpenFOAM，tutorials里面有参考算例
compressibleDriftFoam.zip -
这个模型是不是存在一个问题，就是滑移速度在计算时除了和加速度相关，与雷诺数相关，似乎与别的参数没有什么关系。这就会导致随着加速度的增加，滑移速度急剧增加，这里是否还缺少能量上的约束以避免出现不合理的结果呢
-
我在计算的课题是欠膨胀气液两相流内的激波，想看看滑移速度的影响，有没有可能计算的场景的瞬态特征导致了计算的不稳定啊~
-
请教个问题，就是再滑移速度模型植入时，出现了类似于棋盘一样的相对速度Udm，想问下这种问题该怎么处理呢~
求解器更换成了twoPhaseMixingCentralFoam（Kraposhin Matvey的求解器https://github.com/unicfdlab/hybridCentralSolvers），在其中的YEqn中添加了滑移速度Udm的通量，但是出现了Udm的棋盘一样的变化，不清楚问题出在了哪里，请大佬指条明路~
Udm植入时使用了固定大小的Udc，Udc的方向由加速度方向决定（认为气液两相惯性的差异导致速度的差异）。
这个算例在孔内是50%质量含气率的气液两相，入口压力10MPa，温度383K，（这是个测试，并非真实场景）
10/11