Navigation

    CFD中文网

    CFD中文网

    • Login
    • Search
    • 最新

    在某指定网格加入源项,源项为随时间变化的三角函数

    OpenFOAM
    3
    3
    2100
    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.
    • V
      VamosZhi last edited by 李东岳

      上次的问题解决后,面临这个新问题,最近也找了三角函数的定义,但都是关于边界条件的定义,没有明确的类似案例。
      总的公式仍然是:0_1518392340331_Screenshot-2018-1-31 DTU-1 pdf.png

      只是公式8中的源项换了,换成了:
      0_1518392417894_Screenshot-2018-2-11 DTU-1 pdf.png
      我已经写好程序去找到这个网格,

      label cellI;
      label cellIob;
      DynamicList<label> p1I;
      scalar value(1000);
      vector p1loc(1,0.5,0.5);
      forAll(mesh.cells(),cellI)
      {
            if(mag(U.mesh().C()[cellI]-p1loc)<=value)
            {
              cellIob=cellI; 
              value=mag(U.mesh().C()[cellI]-p1loc);
            }
      }
      p1I.append(cellIob);
      Info<<"location"<<U.mesh().C()[cellIob];
      

      但是如何在定义三角函数源项以及在网格中添加该源项仍然难。如有了解,希望提示一下。谢谢!

      zym604 1 Reply Last reply Reply Quote
      • zym604
        zym604 @VamosZhi last edited by

        @vamoszhi 0_1518409145264_770044e1-808d-47f0-883e-7fa553c8a2a0-image.png
        0_1518409156721_e340af6d-62dc-44f8-99de-0810a1f9a20f-image.png
        我一开始没看明白,是因为你这个公式里面的符号跟一般书上的不一样。

        1. 看上去你这里的u'反而是average量,U_i才是湍动量。是么?
        2. 公式7的雷诺应力tau是不是少了一项,2/3Kdelta_ij?
          0_1518409365573_3956194b-bbe9-4faf-91c6-d4ee33b901e2-image.png
        3. 你的意思是只在某个cell上套用公式9,其他cell继续用公式8么?
        4. 我对能量方程不是很熟,公式8里面有P和p',这就是两个变量了,你是不是还用了一些别的closure law来表达P变量了?
        1 Reply Last reply Reply Quote
        • Cloud1990
          Cloud1990 last edited by 李东岳

          但是建议做个sourceMask,能量方程可以写成
          \begin{equation}
          LHS=-\frac{\partial p}{\partial t}\times (sourceMask-1)-S_c\times sourceMask
          \end{equation}

          1 Reply Last reply Reply Quote
          • First post
            Last post

          CFD中文网 | 东岳流体 | 京ICP备15017992号-2
          论坛登录问题反馈可联系 li.dy@dyfluid.com