OpenFOAM与python的耦合



  • 我想用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();这句可以删掉不管,别的不需要改动了。



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



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

    为何要用python,你用python算的湍流?



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



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


Log in to reply
 

CFD中文网 2016 - 2020 | 京ICP备15017992号-2