Navigation

    CFD中文网

    CFD中文网

    • Login
    • Search
    • 最新

    RANS动量方程中雷诺应力表达的问题

    OpenFOAM
    2
    4
    704
    Loading More Posts
    • Oldest to Newest
    • Newest to Oldest
    • Most Votes
    Reply
    • Reply as topic
    Log in to reply
    This topic has been deleted. Only users with topic management privileges can see it.
    • F
      fangyuanaza last edited by

      RANS中张量形式动量方程为:
      798afe64-8bcc-426b-80a0-a9086108325d-image.png
      在rhoSimpleFoam求解器中对应为UEqn.H
      368d0d2c-93e1-428a-9448-a00f28ef6d11-image.png
      由于压力梯度项一般不直接放入UEqn中,而是对动量方程进行修正,因此turbulence->divDevRhoReff(U) 对应07bd23db-8f4c-4d77-b1e5-afb046e5e712-image.png 前一项可用拉普拉斯符号表示,后一项实际上就是雷诺应力。
      问题一:我发现turbulence->divDevRhoReff在/home/dyfluid/OpenFOAM/OpenFOAM-7/src/TurbulenceModels/turbulenceModels/路径下的linearViscousStress/linearViscousStress.C及ReynoldsStress/ReynoldsStress.C中均有体现,如果我要修改动量方程中的雷诺应力项,那么我两个都需要修改?
      问题二:turbulence->divDevRhoReff对应的代码是
      335942ce-0263-449f-a98d-54fbf388ee5b-image.png 这里的第一项也就是雷诺应力项表达式为
      04392fb8-9aaf-4137-9af7-8349a57fee31-image.png
      这为什么与Boussinesq假定(如下)不一致呢?
      a574197b-04cf-4ac4-bbc2-2f7402831d04-image.png
      Boussinesq假定的第一项b2d0bd9f-1d34-4c49-8230-66db725d6e7e-image.png为什么在代码中没有体现呢?

      F 1 Reply Last reply Reply Quote
      • F
        fangyuanaza @fangyuanaza last edited by

        @fangyuanaza https://caefn.com/openfoam/diffusion-term1 我找到了这方面的解析,说是这一项放到了压力梯度中,但是我还是没有找到

        1 Reply Last reply Reply Quote
        • 李东岳
          李东岳 管理员 last edited by

          没有显性的操作,只不过是存在压力p里面,但是没有显性的说p=p+1/3k之类。所以,其实求解出来的p,是p+1/3k

          CFD课程 改成线上了 http://dyfluid.com/class.html
          CFD高性能服务器 http://dyfluid.com/servers.html

          F 1 Reply Last reply Reply Quote
          • F
            fangyuanaza @李东岳 last edited by

            @东岳 懂了~ 谢谢智慧的李老师~~

            1 Reply Last reply Reply Quote
            • First post
              Last post

            CFD中文网 | 东岳流体 | 京ICP备15017992号-2