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



  • 上次的问题解决后,面临这个新问题,最近也找了三角函数的定义,但都是关于边界条件的定义,没有明确的类似案例。
    总的公式仍然是: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];
    

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



  • @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变量了?

登录后回复
 

与 CFD中文网 的连接断开,我们正在尝试重连,请耐心等待