CFD中文网

各位前辈，大佬们，最近在学习中遇到了几个问题希望各位前辈指教一下：
目前在大多数文献中我看到气液两相流中的曳力模型系数都是以下形式：
$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项，在文献里也没看到相关的引用。
是否是由于这些代码一般是用来计算液态金属-气体两相流，因此比较特殊呢？
希望各位大佬不吝赐教，非常感谢！:xiexie:

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

:135: :135:

是的。

openfoam这面PU的求解框架大体一样。

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

@李东岳 :xiexie: 感谢李老师

@李东岳 :135: 不好意思，李老师，再打扰您一下，我还有一个问题，我目前使用的是v10版本的multiphaseEulerFoam，我在dragModel里可以通过interface_.dispersed().d()这种来调用直径参数，现在我想在求解器的.c程序里调用它，应该用什么方式呢
我注意到您再小代码里提到过以下代码

meanDiameter { type coded; libs ("libutilityFunctionObjects.so"); name error; codeExecute #{ const volScalarField& d = mesh().lookupObject<volScalarField>("d.alpha.oil"); scalar d32 = d.weightedAverage(mesh().V()).value(); if (Pstream::master()) { std::ofstream file; file.open ("d32", std::ofstream::out | std::ofstream::app); file << mesh().time().timeName() << " " << d32 << "\n"; file.close(); } #}; }

这段我是否能用呢

请问各位大佬，是否有phaseProperties设置教程呢？

of9中多相欧拉的算例，在constant/phaseProperties 中都要设置一个type；
已知的可选设置有:
basicMultiphaseSystem;
interfaceCompositionPhaseChangeMultiphaseSystem;
interfaceCompositionPhaseChangePopulationBalanceMultiphaseSystem;
populationBalanceMultiphaseSystem;
thermalPhaseChangeMultiphaseSystem;
thermalPhaseChangePopulationBalanceMultiphaseSystem.

目前比较关心的是什么时候使用interfaceCompositionPhaseChangeMultiphaseSystem，什么时候使用thermalPhaseChangeMultiphaseSystem呢？这两者具体有什么不同呢？

还有就是blending相里，算例有时会设置传热和曳力的离散/连续分数设置，有时又不会，这是与传热和曳力的模型有关吗
万分感谢:xiexie:

@李东岳
已经看到了，非常感谢您:xiexie:

@李东岳 @bestucan
感谢两位老师，我大概明白了，马上就去详细看一看:xiexie: