OpenFOAM与python的耦合

zym604

我想用python生成一个湍流模型的替代模型，放到OpenFOAM中，让OpenFOAM用我的python模型算动量方程中的雷诺应力项。请大家给我一些建议：
我打算修改pisoFoam.C

        fvVectorMatrix UEqn
        (
            fvm::ddt(U)
          + fvm::div(phi, U)
          + turbulence->divDevReff(U)
        );

改成

        fvVectorMatrix UEqn
        (
            fvm::ddt(U)
          + fvm::div(phi, U)
          - fvm::laplacian(nu+mynut, U)
        );

mynut是我算出来的一个tensorField，跟LRR里面一样处理。

但是我不确定是不是只改这里就行了，后面有句：

    turbulence->correct();

李东岳

turbulence->correct();这句可以删掉不管，别的不需要改动了。

zym604

@李东岳多谢！我可以问问turbulence->correct()这句是干什么的吗？

李东岳

turbulence->correct()要执行correct()函数，这个函数里面执行其他函数。。。

为何要用python，你用python算的湍流？

zym604

@李东岳是的，我用python写了个回归函数，根据DNS给出的数据做回归，然后想往OpenFOAM里带。

D.Benjamin

@东岳请问李老师，turbulence->correct()中的correct()函数具体是怎么实施的？源码在哪里可以找到？谢谢