CFD中文网

@wang123sd 不是啊😁

@李东岳 在 OpenFOAM中移动接触线问题（接触角）的边界条件设置 中说：

我去，竟然来了！！

别人提问的话，我这边能收到邮件呢。工作太忙，顾不上逛帖子了。

另外的话呢，对于像OpenFOAM和fluent这类使用CSF模型来计算气液两相流的话，低速问题时虚拟流动（spurious current）影响较大，做学术研究的话精度恐怕不够。建议用Basilisk去做，精度很高。但是工程问题的话，商软或者OpenFOAM是可以的。其他看到的吧友可以借鉴一下。
我目前已经不做单纯的界面流动问题了，研究方向稍微有了变化。

@闹市闲民 再给推荐一本粘性流体力学教材，个人亲测食用最佳:chitang: 7d7a70c8-f0cf-4676-b460-7866a9178ba8-image.png

我这个是一个圆柱，分成1/26，成为周期对称模型。中间是空气，不能省略。

@dabiu 在 关于icoFOAM中通量phi的更新 中说：

icoFOAM中更新通量phi采用的如下语句（注：我理解这个语句）

if (piso.finalNonOrthogonalIter()) { phi = phiHbyA - pEqn.flux(); }

为什么不在下面这条语句的后面，直接用phi = fvc::flux(U)来更新通量，这样更新跟上面那条语句有什么区别吗？用phi=fvc::flux(U)会有什么问题吗？

U = HbyA - rAU*fvc::grad(p); U.correctBoundaryConditions();

初学icoFOAM，求各位高手解答，OpenFOAM用的版本为OpenFOAM 7.

我觉得这里的phi是经过Rhie-Chow插值修正过的面通量。

截断误差首项阶次一样

我推了一下，CD的更小。不一样啊。

捕获.PNG

@李东岳 明白了，谢谢李老师耐心细致的解答。

谢谢您的解答，在openFoam v6版本里icoUncoupledKinematicParcelFoam这个求解器里，有pairCollision(软球模型),我现在尝试用软球模型先试试。
sprayFoam里面是有O’Rourke碰撞的，但我做的东西和化学反应没啥关联，在后面做的过程中如果pairCollision没啥进展，可以借鉴您的建议把O’Rourke移植过去。

#include "udf.h"
DEFINE_EC_RATE(user_ec_rate, f, fthread, r, V, current, didV, Eeq)
{
double alpha_a = r->alpha_a, alpha_c = r->alpha_c;
double io = r->io
double T = F_T(f,fthread);
double arg1, arg2;
cxboolean tafelmethod = r->tafelmethod;
int i;
double eta;
if (tafelmethod)
{
alpha_a = 2.303 * UNIVERSAL_GAS_CONSTANT * 298.15 /(alpha_a * FARADAY_CONSTANT);
alpha_c = 2.303 * UNIVERSAL_GAS_CONSTANT * 298.15 /(alpha_c * FARADAY_CONSTANT);
}
Eeq = r->Eeq;
eta = V - Eeq;
for(i = 0; i<r->n_reactants; i++)
if( ABS( r->exp_reactant[i] ) > SMALL_S )
{
int ni = r->reactant[i];
io = pow((F_YI(f,fthread,ni)/MAX(r->yi_ref[ni],SMALL) + 1.0e-20), r->exp_reactant[i]);
}
for(i = 0; i<r->n_products; i++)
if( ABS( r->exp_product[i] ) > SMALL_S )
{
int ni = r->product[i];
io = pow((F_YI(f,fthread,ni)/MAX(r->yi_ref[ni],SMALL) + 1.0e-20), r->exp_product[i]);
}
arg1 = FARADAY_CONSTANT / (UNIVERSAL_GAS_CONSTANTT);
arg2 = arg1eta;
current = io( exp( arg2alpha_a ) - exp( -arg2alpha_c ) );
didV = io( arg1alpha_aexp( arg2alpha_a ) + arg1alpha_cexp( -arg2alpha_c ) );
/* If multiple electrochemical reactions are used, you can define rate for each reaction
using the following if-statement /
/
if (STREQ(r->name, "reaction-1"))
{
...
}
else if (STREQ(r->name, "reaction-2"))
{
...
}
*/
}

@cfd-china 这个问题，有意思。去知乎上搜一下

哈哈，多谢老铁帮吹流弊！你们都是真流弊，我是吹流弊 :140:

@李东岳 收到，谢谢李老师。

我导收了一个湘潭大学的搅拌器的英文sci。实验+cfd，用openfoam搞的。我针对算例设置这一段提出了1个不痛不痒的问题，不能一个不提啊。目前已返回我导。估计明后天就返回编辑吧？哈哈。我导没有抵制，不过提了2页问题...anyway，问题不大。:sunglasses:

@alex 看到你的网格了，感谢分享，感觉划分的非常不错，好厉害！学习了！

网格质量如何？

编译方式不一样？

@cfd-china 门槛啊，真高！

RSM 怎么样？ 对雷诺数没有的限制

提一下，目前可以通过蒙特卡洛求解、矩方法求解、最小二乘法求解、加权残数法求解。不知通过傅里叶变换、拉普拉斯变换求解是否有这个可能？

你查一下这方面的实验论文吧？
壁面速度低，容易造成管壁聚集