CFD中文网

Cloud1990

@cloud1990 在 VOF求解自由水面问题，求解空气区域吗？中说：

Left: Correct jet shape predicted by TruVOF technique used in FLOW-3D. Right: Incorrect jet shape predicted by pseudo-VOF technique used by other CFD codes.

Left: FLOW-3D‘s TruVOF technique predicts jet impingement on wall and some outflow. Right: Pseudo-VOF methods don’t predict realistic jetting of fluid on side walls.

从结果上看，flow3D 所谓的TruVOF是有问题的，air velocity 和　air resistance = 0?　这个可能比较适合粗网格计算，作为一个tricky的存在．但理论上，这个TruVOF是无法解决wind wave的问题．

Cloud1990

Cloud1990

但是建议做个sourceMask，能量方程可以写成
\begin{equation}
LHS=-\frac{\partial p}{\partial t}\times (sourceMask-1)-S_c\times sourceMask
\end{equation}

Cloud1990

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()．

Cloud1990

@李东岳在 turbulence里k设置的问题中说：

你这个可以通过codedfixedvalue来实现，在论坛里有代码你找找这个关键词。我现在手机写不了代码，找不到再回帖。

还有一个方法上是用timeVaryingMappedFixedValue这个boundary condition。

参见tutorials/incompressible/simpleFoam/pitzDailyExptInlet。

需要在constant里建立一个boundaryData的文件，一般只要给定0 时刻分布，就相当于给定constant分布。

boundaryData中点的坐标和具体数值可用python或matlab自行生成。

必须注意的是，boundaryData里面文件的具体格式，不同版本OpenFOAM可能不同。

CFD中文网

Cloud1990

帖子