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

1. 我不知道这句有什么用，按照我的需求，是不是要把这句删掉？
2. 还有什么其他的地方需要改吗？

[李东岳]

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

[zym604]

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

[李东岳]

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

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

[zym604]

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

[D.Benjamin]

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