如何给时间空间二阶精度两步映射法的求解器添加紊流模型和VOF？

[OItoCFD]

各位老师和同学，
我用Admas-bashforth格式写了一个新的求解器，用的两步映射法，没用PISO或pimple，简要代码如下：

        fvVectorMatrix UEqn
        (
          fvm::ddt(U)
          + 1.5*fvc::div(phi, U) - 0.5*fvc::div(phi.oldTime(), U.oldTime())
          == 1.5*fvm::laplacian(nu, U) - 0.5*fvc::laplacian(nu, U)
          + 0.5*fvc::grad(p)
        );

        fvScalarMatrix pEqn
        (
            fvm::laplacian(p) == 2.0/3.0*fvc::div(phi)/dt
        );

请问，

1. 如何添加紊流模型？
   直接按照原始求解器，建立紊流模型的指针后，获取nuEff(), 把nu替换成nuEff()，这样是正确的吗？
2. 如何添加VOF？
   我还没尝试，不知道直接仿照interFoam，把该包括的子文件包括进去，Ueqn加上重建的那一部分，是不是可行，有没有什么误区。

谢谢各位！

[李东岳]

自己添加VOF比较费劲，这个要做到跟原始求解器一样的稳健，需要弄好多东西。最好是从interFoam开始，把piso改成你的压力速度耦合求解。不过要注意phi的更新。

[OItoCFD]

@李东岳 好的谢谢老师，那我添加紊流模型，直接nu用nuEff()代替可以吗，我就调用下LES，不用RANS。

[李东岳]

@oitocfd 可以

[OItoCFD]

@李东岳 谢谢！