CFD中文网

请问各位前辈，有没有界面蒸发源项的模型啊，找了好多都是壁面成核的，因为要做的工况是在气液界面处产生相变，而不是壁面加热导致相变的，看了一些文章好像找到没有one-group下界面蒸发的模型

各位前辈，大佬们，最近在学习中遇到了几个问题希望各位前辈指教一下：
目前在大多数文献中我看到气液两相流中的曳力模型系数都是以下形式：
$k=\frac{3}{4}\alpha_d\rho_c\frac{C_D}{d_d}\left|U_c-U_d\right|$
但是最近学习文献的时候在核口的一些CFD代码文章比如SIMMER-III，AFDM，还有一些其他程序里出现了另外一种形式：
$k=A+B*\left|U_c-U_d\right|$
B项一般是跟$C_d$相关，在SIMMER-III中其形式为：
$B=\frac{1}{8}a_i\rho_cC_D,a_i是交界面面积$
而A项给出的解释一般是comes from Stokes law for viscous flow，形式一般是：
$A=\frac{3\alpha_c\mu_c}{d_d}$
比较好奇这种形式的曳力系数是怎么得到的，尤其是A项，在文献里也没看到相关的引用。
是否是由于这些代码一般是用来计算液态金属-气体两相流，因此比较特殊呢？
希望各位大佬不吝赐教，非常感谢！

李老师在另一篇帖子有说过
https://cfd-china.com/post/11859

我记得setField是在求解器运行前用的吧，直接在setFieldDict指定好box以及相关变量，然后在算例文件下blockMesh setField xxxFoam这样用吧，如果单纯需要初始化场应该不用修改求解器来着吧，认识有限，说错了希望大佬们指正

这个是我想求解的方程形式，这段时间看了点相关东西，请问各位前辈，foam-extend里的blockMatrix能实现我想要的求解方式吗，然后就是曳力项应该OF里应该怎么写呢？是fvm : : Sp ( K, Ug ) - fvm : : Sp ( K, Ul )吗？希望各位大佬不吝赐教啊

请教各位大佬一个问题，两相流，在动量方程这块想处理相间作用力，OF我目前了解到的求解器是采用显式处理，

，但是我想做成隐式处理，弄成两个类似于a1Ul+b1Ug=c1,a2Ul+b2Ug=c2的矩阵,消元后再解线性方程组，这种有实现的可能吗，还是说这种应该做成迭代求解的模式呢？非常感谢!

已经没问题了

各位大佬，我想请教下关于EOS的事，我老师想让我弄一个算液态水密度的函数，我之前根据IAPWS97写了一个，但是发现在临界区算出来的密度跟实际好像差了很多，然后又看了下P-R eos,请问各位有这种根据p T计算水密度的相关资料推荐下吗

@李东岳 您的意思是在EEqn这里求解出的温度，他就默认是最终的正确温度是吗？然后在pEqn这里只是对质量速度压力密度这些的迭代求解，这里就相当于buoyantPimpleFoam里的比较EOS得到的rho和连续性方程的rho这个过程。

春节期间没啥事，或许我可以抽时间把openfoam中的算法拔出来。

期待！感谢李老师！

@李东岳 我的理解里，传质的过程应该伴随着能量的变化吧，求出正确的p后才能确定对应的温度这些？
能量方程右侧也有传质的源项，如果传质没求解的话怎么得到正确的温度呢
我的认识非常有限，如果错了请李老师批评指正！

@李东岳
reactingTwoPhaseEulerFoam求解器的相变求解过程
李老师，我想请问下在org这边的multiphaseEuler，pEqn里既然考虑了传质，但是为啥找不到对能量进行的修正呢，从代码里似乎只能看到在EEqn求解后进行了thermo.correct()，是我理解有问题吗
还有就是能否请您大概提示下这个求解器的求解过程或者pEqn的迭代过程，万分感激！

@李东岳 谢谢李老师，我先去找点相关的文章看一看吧

各位老师，前辈，我想学习下multiphaseEuler,或者是reactingTwoPhaseEulerFoam的算法，请问各位前辈有什么相关的文章可以推荐下吗？非常感谢
阅读了xiaopingqiu老师的博客，也看了东岳老师的reactingTwoPhaseEulerFoam的解析，有点疑惑,东岳老师写的连续性方程

ddt(alpha,rho) + div(alphaRhoPhi) = 0

如果考虑传质的话，右侧为什么没有源项呢？

@zy-gong 好像确实可以，太惭愧了，总是问这种低级问题，学习之路任重而道远
谢谢李老师了

还是说这种也可以直接用lookupobject呢？

@李东岳
抱歉李老师，我重新组织下语言
以我有限的认识，因为求解器用的是heRhoThermo,无论是new()还是correct()似乎最后都会调用heRhoThermo::calculate()来计算热物性
而我想植入的一个气相的eos，需要的参数是液相的温度和气相的温度，大概就是rho(Tgas,Tliquid)，植入eos的方法我大概已经掌握了，目前的想法是建立一个新的thermoType，newThermo的calculate()函数是用两个温度场作为参数来计算rho的，但是我目前不知道怎么在thermo里获取另一相的温度
假设液相的温度，我可以在多项欧拉求解器里通过fluid.phases()[0].T()来访问，那在newThermo里应该怎么做才能读到液相的温度场呢？
感谢李老师!

目前在往OF里植入一些模型，因为在使用多项欧拉求解器，想植入一个EOS是根据气液两相的温度来计算气相的rho
现在想的方法是基于原有的heRhoThermo,弄一个user_lib，改一下原有的calculate函数，把rho(p，T)，换成rho(T,T)，但目前碰到的问题是我不知道该怎么在我的thermo里声明fluid.phase()[0].T()这个场
请问各位前辈，这种应该怎么写呢？或者各位前辈有没有更好的方法
感谢各位前辈的指教，万分感激

@李东岳 感谢李老师

@李东岳 李老师，我现在用.boundaryField()[patchi][facei]对边界进行赋值，总是提示我read-only，我需要用什么特定边界条件吗，或者是我的语句不对？

@zy-gong 哦哦我反应过来了，用mesh.oner()就可以了，打扰李老师了