CFD中文网

    CFD中文网

    • 登录
    • 搜索
    • 最新

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

    OpenFOAM
    2
    4
    841
    正在加载更多帖子
    • 从旧到新
    • 从新到旧
    • 最多赞同
    回复
    • 在新帖中回复
    登录后回复
    此主题已被删除。只有拥有主题管理权限的用户可以查看。
    • F
      fangyuanaza 最后由 编辑

      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 条回复 最后回复 回复 引用
      • F
        fangyuanaza @fangyuanaza 最后由 编辑

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

        1 条回复 最后回复 回复 引用
        • 李东岳
          李东岳 管理员 最后由 编辑

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

          线上CFD课程开始报名:http://www.dyfluid.com/class.html

          CFD高性能服务器 http://dyfluid.com/servers.html

          F 1 条回复 最后回复 回复 引用
          • F
            fangyuanaza @李东岳 最后由 编辑

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

            1 条回复 最后回复 回复 引用
            • First post
              Last post