Navigation

    CFD中文网

    CFD中文网

    • Login
    • Search
    • 最新

    dnsFoam初始化问题

    OpenFOAM
    5
    9
    4903
    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.
    • W
      wzchen last edited by

      刚刚接触openfoam,问个简单的问题:对于tut里面的dnsFoam的例子,初始化时内部速度场的数据是怎么来呢

      C 1 Reply Last reply Reply Quote
      • C
        CFD中文网 @wzchen last edited by

        @wzchen

        你好,你可以把速度设置为uniform后,使用boxTurb程序来生成人工随机量。算例自带的例子只不过帮你完成了这个过程。:laughing:

        CFD中国标准用户测试帐号
        目前由徐笑笑登录

        1 Reply Last reply Reply Quote
        • W
          wzchen last edited by

          :happy: 谢谢,在想问一下
          0_1461811236659_upload-a8c2479c-7dcc-49a1-bc7f-e32db3586c1d
          Ea k0的具体含义是什么,如何设置它们的值?有关于openfoam初始化方面的资料推荐一下吗:happy:

          李东岳 T 2 Replies Last reply Reply Quote
          • 李东岳
            李东岳 管理员 @wzchen last edited by

            @wzchen

            我半个月前在写dnsFoam的资料,不过后来被别的任务打断了。dnsFoam模拟的是一盒子湍流,就是不带壁面的DNS模拟。因为没有壁面所以能量会被消耗掉导致无湍。dnsFoam在源项中添加了人工噪声,噪声通过快速傅里叶变换来生成力来保持湍流不被消耗掉。目前我只进行到这里。并没有对公式进行推导。最近事情非常非常的多。做一个类似的dnsFoam解析可能要2,3个月之后了。如果你感兴趣自己推导最好了!:cheeky:

            另外,dsFoam的资料非常少,目前仅有一个查尔姆斯大学的tutorial,我建议你看你专业的DNS文献,里面提及了如何在源项添加力来保证有湍流。

            线上CFD课程 7月1日报名截止 http://dyfluid.com/class.html
            CFD高性能服务器 http://dyfluid.com/servers.html

            1 Reply Last reply Reply Quote
            • W
              wzchen last edited by

              好的,谢谢:happy:

              1 Reply Last reply Reply Quote
              • T
                TangShangyu @wzchen last edited by

                @wzchen 请问您知道如何设置Ea和k0的值了吗?方便的话,请指教,谢谢。

                1 Reply Last reply Reply Quote
                • 李东岳
                  李东岳 管理员 last edited by

                  他们的计算在这个文件中Ek.H,目前如何指定的理论背景不详。

                  namespace Foam
                  {
                  
                  inline tmp<scalarField> Ek
                  (
                      const scalar Ea,
                      const scalar k0,
                      const scalarField& k
                  )
                  {
                      tmp<scalarField> tEk = Ea*pow(k/k0, 4.0)*exp(-2.0*sqr(k/k0));
                  
                      /*
                      scalarField& Ekf = tEk();
                  
                      label i;
                      forAll(Ekf, i)
                      {
                          if (k[i] < 2 || k[i] > 10)
                          {
                              Ekf[i] = 0.0;
                          }
                      }
                      */
                  
                      return tEk;
                  }
                  

                  线上CFD课程 7月1日报名截止 http://dyfluid.com/class.html
                  CFD高性能服务器 http://dyfluid.com/servers.html

                  W 1 Reply Last reply Reply Quote
                  • W
                    Will_Wang @李东岳 last edited by

                    @东岳 他这个是Von-Karman风速谱傅里叶变换的经验公式,李老师我想问一下这个模型算转捩的话扰动应该怎么给会好一些呢

                    李东岳 1 Reply Last reply Reply Quote
                    • 李东岳
                      李东岳 管理员 @Will_Wang last edited by

                      不太了解这个Von-Karman风速谱傅里叶变换公式。你哪有参考文献么?

                      线上CFD课程 7月1日报名截止 http://dyfluid.com/class.html
                      CFD高性能服务器 http://dyfluid.com/servers.html

                      1 Reply Last reply Reply Quote
                      • First post
                        Last post

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