CFD中文网

    CFD中文网

    • 登录
    • 搜索
    • 最新

    SRF求解器的使用

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

      最近在做一个旋转的后台阶流动,如下图所示,旋转的位置在台阶下角落处,我想采用SRF求解器来计算,但是对于其设定有一些疑惑。
      该算例中,流体是沿台阶前的直管流入的,旋转角速度是1000多rpm,我想问的是在设定Urel中进口速度时,该采用什么样的边界条件,是采用SRFvelocity还是fixedvalue?现在我采用fixedValue,计算结果完全没有旋转的迹象
      还有那些地方需要注意,谢谢!
      23685536-cfef-4ee3-b816-9e9de8006aa1-image.png

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

        SRFvelocity,给定的值采用U来给定,然后SRFvelocity会计算Urel。你这个算例几何比较有意思,你试试看,调节不好我给你调

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

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

        benqing 1 条回复 最后回复 回复 引用
        • benqing
          benqing 最后由 编辑

          前辈,我不太明白是什么意思,在这个求解器里面,0文件里面的速度必须是Urel,不能给定U,而且我看资料的话,SRFvelocity给定的相当于是绝对速度,我这算例中,应该是给定相对速度吧,这个case我现在是在of230中算的,版本有影响吗?

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

            @东岳 如果采用SRFVelocity,那里面的relative是设置为yes还是no啊

                inlet
                {
                    type            SRFVelocity;
                    inletValue      uniform ( 5.6 0 0 );
                    relative        no;
                    value           uniform ( 0 0 0 );
                }
            
            1 条回复 最后回复 回复 引用
            • 李东岳
              李东岳 管理员 最后由 编辑

              看这个描述。如果是yes,那你的Urel就是你给的值,如果是yes,那你就给出绝对速度,然后Urel会进行一定的计算后得出真实值。

              你的应该是No

              Description
                  Velocity condition to be used in conjunction with the single
                  rotating frame (SRF) model (see: SRFModel class)
                  Given the free stream velocity in the absolute frame, the condition
                  applies the appropriate rotation transformation in time and space to
                  determine the local velocity.
                  The optional \c relative flag switches the behaviour of the patch
                  such that:
                      - relative = yes: inlet velocity applied 'as is':
                      \f[
                          U_p = U_{in}
                      \f]
                      - relative = no : SRF velocity is subtracted from the inlet velocity:
                      \f[
                          U_p = U_{in} - U_{p,srf}
                      \f]
                  where
                  \vartable
                      U_p     = patch velocity [m/s]
                      U_{in}  = user-specified inlet velocity
                      U_{p,srf} = SRF velocity
                  \endvartable
                  \heading Patch usage
                  \table
                      Property     | Description             | Required    | Default value
                      inletValue   | inlet velocity          | yes         |
                      relative     | inletValue relative motion to the SRF? | yes     |
                  \endtable
                  Example of the boundary condition specification:
                  \verbatim
                  myPatch
                  {
                      type            SRFVelocity;
                      inletValue      uniform (0 0 0);
                      relative        yes;
                      value           uniform (0 0 0);    // initial value
                  }
                  \endverbatim
              

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

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

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

                @东岳

                Given the free stream velocity in the absolute frame,
                

                这个意思不是给定的是绝对坐标系中速度吗?但是我的算例中,入口是随着系统在旋转的,所以采用SRFVelocity合理吗

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

                  如果进口的$U_{in}$是向右1m/s,凑巧你那个位置的SRF $U_{p,srf}$也是1m/s,那你的Urel不就是0了么。

                  主要是你要提供Urel的边界,如果你知道Urel的值,用fixedValue也可以。SRFVelocity只不过帮你计算了一下

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

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

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

                    @东岳 好的,我先尝试一下

                    1 条回复 最后回复 回复 引用
                    • benqing
                      benqing 最后由 benqing 编辑

                      一年多了还是没有解决这个算例,不知道有没有人做过相类似的,重新描述下这个算例,如图所示一个后台阶流动,但是是在图中左下(0,0,0)位置为原点绕Z轴旋转,主要流动参数,
                      雷诺数$Re_b=\frac{U_b(H-h)}{\nu}=5600$,其中$U_b$是体积平均速度;
                      旋转数$U_b=\frac{\Omega(H-h)}{U_b}$,$\Omega$是角速度;
                      台阶前段槽道流部分壁面剪切雷诺数$Re_\tau=\frac{u_\tau(H-h)}{\nu}=360$,这个与槽道流DNS结果是完全一致的,也就是说台阶前是完全发展的槽道流。目前存在以下问题:
                      1、求解器选择,SRF类还是需要自己修改求解器?
                      2、合理的进口条件,不旋转的相对好给,但是旋转的case,速度进口边界条件和初始值不知道如何给合理。
                      3、无量纲化的流向速度为什么会小于1,如下图所示。以我的经验,槽道流里面最大速度比体积平均速度肯定是大的,也就是$\frac{U_{max}}{U_b}$是大于1的,但是为什么在旋转流动里面就会于1呢?
                      6c444575-46fa-460b-bb57-8d655a396d45-image.png
                      参考文献
                      1、DNS原始文献:Turbulent flow over a backward-facing step. Part 1. Effects of anti-cyclonic system rotation
                      2、图片来自:Development and Application of Rotation and Curvature Correction,是篇博士论文

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