CFD中文网

    CFD中文网

    • 登录
    • 搜索
    • 最新

    codedFixedValue这样写对吗?

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

      请教各位老师,在文献中看到,射流入口用了双曲正切函数分布,我用codedFixedValue这样写对吗?
      我编译可以通过,但是 plot over line看了下分布,和直接给速度定值喷射好像没有差别。

      文献定义

      文献中的速度分布

          inletJet
          {
              type            codedFixedValue;
              value           uniform (0 0 0); 
      
              name            tanhInletJet;    
      
              code
              #{
                  const fvPatch& boundaryPatch = patch();
                  const vectorField& Cf = boundaryPatch.Cf();
                  vectorField& field = *this;
      
                  scalar Uj_y = 1250, r_0 = 0.001, theta_0 = 0.00009;   
      
                  forAll(Cf, faceI)
                  {
                      const scalar x = Cf[faceI][0]; 
                      const scalar z = Cf[faceI][2]; 
                      const scalar r = sqrt(x*x+z*z); 
      
                      field[faceI] = vector(0,Uj_y*0.5*(1-tanh(1/(4*theta_0)*((r/r_0)-(r_0/r)))),0); 
                  }
              #};
      
              // codeInclude
              //#{
              //    #include "fvCFD.H"
              //    #include <cmath>
              //    #include <iostream>
              //#};
      
              // codeOptions
              //#{
              //    -I$(LIB_SRC)/finiteVolume/lnInclude
              //    -I$(LIB_SRC)/meshTools/lnInclude 
              //#};
          }
      
      W 1 条回复 最后回复 回复 引用
      • W
        wangfei9088 @尚善若水 最后由 wangfei9088 编辑

        @尚善若水 代码没错,是公式(3.2)的原因。分析一下:
        当 r=r0 时,tanh[..]=0, w/Wj=0.5;当r=0时,tanh[...]=-1,w/Wj=1。这个曲线是可以用excel或matlab画出的。
        我觉得这个公式对应的应该是图中Elliptic (major axis)这条曲线,实质上也就是你说的和直接给速度定值喷射一样。要得到Circular曲线也可以,把公式(3.2)中的r除以2即可。

        尚 1 条回复 最后回复 回复 引用
        • 尚
          尚善若水 @wangfei9088 最后由 编辑

          @wangfei9088 公式应该没问题,它这个处理好像就是在径向有一个微弱的速度衰减。

          1 条回复 最后回复 回复 引用
          • 梦中飞翔的阿涛
            梦中飞翔的阿涛 最后由 编辑

            大佬,我刚刚发了一个帖子,能不能帮我看一下,关于煤粉燃烧的。之前你帮我回答过相关问题,感谢了。:xinxin:

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