能问下关于湍流方程是否可以添加源项吗?
-
最近想添加一个湍动能的源项,不知道是否可以直接在k方程添加源项呢?
-
可以用 fvm::susp 添加。
-
@wwzhao 准备怼of湍流模型的调用机制,导致扫了一下K-equation 是可以添加源项的,但是具体的还不是很清楚,比如说如何引入源项(因为源项是通过某个数学模型计算得到的)。
-
@mohui 打开源码,找到kEqn,添加源项,完成。
-
@wwzhao 谢谢,我看到这个方程了,现在我在想如何引进源项。
-
@mohui 你可以把源项 IOobject 注册到内存中,再在湍流模型里通过 lookupObject 查找得到源项。
-
@wwzhao 好的,谢谢。现在稍微了解到3.0以后版本的湍流模型的编译较为复杂,需要仔细看看。
-
用pimpleFoam的话可以直接在fvoptions里添加momentum source项。
-
@zym604 动量方程和湍动能方程是两个不同的概念,fvOptions只能在动量方程中使用。因为只有动量方程的源项中加了
fvOptions(U)
,而湍动能方程没有。 -
@wwzhao 啊,sorry,没看清题目。
-
@mohui 也不复杂吧,比如对于kepsilon,用incompressible turbulence模型,
1 在ras目录下直接在kepsilon做更改,
2 回到incompressible目录下,编译
3 回到求解器,在编译一下
ok 还算挺方便的
当然也可以添加新的文件 mykepsilon,需要在incompressible目录下修改某子文件夹下的 xxxmodels 添加两句代码,就可以了 -
@wwzhao 其实很多模型里面已经有ksource项了,当然默认下这项为零,可以直接将这个source给激活就行了。虽然我没用过
-
of5.0里面已经有fvOption了,
tmp<fvScalarMatrix> epsEqn ( fvm::ddt(alpha, rho, epsilon_) + fvm::div(alphaRhoPhi, epsilon_) - fvm::laplacian(alpha*rho*DepsilonEff(), epsilon_) == C1_*alpha()*rho()*G*epsilon_()/k_() - fvm::SuSp(((2.0/3.0)*C1_ - C3_)*alpha()*rho()*divU, epsilon_) - fvm::Sp(C2_*alpha()*rho()*epsilon_()/k_(), epsilon_) + epsilonSource() + fvOptions(alpha, rho, epsilon_) );
或者复制模仿realizableKE的格式,赋值于epsilonSource() 或者 kSource().
论坛登录问题反馈可联系 li.dy@dyfluid.com