Navigation

    CFD中文网

    CFD中文网

    • Login
    • Search
    • 最新

    能问下关于湍流方程是否可以添加源项吗?

    OpenFOAM
    5
    13
    5614
    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.
    • M
      mohui last edited by

      最近想添加一个湍动能的源项,不知道是否可以直接在k方程添加源项呢?

      1 Reply Last reply Reply Quote
      • W
        wwzhao 教授 last edited by

        可以用 fvm::susp 添加。

        M 1 Reply Last reply Reply Quote
        • M
          mohui @wwzhao last edited by

          @wwzhao 准备怼of湍流模型的调用机制,导致扫了一下K-equation 是可以添加源项的,但是具体的还不是很清楚,比如说如何引入源项(因为源项是通过某个数学模型计算得到的)。

          W 1 Reply Last reply Reply Quote
          • W
            wwzhao 教授 @mohui last edited by

            @mohui 打开源码,找到kEqn,添加源项,完成。

            M 1 Reply Last reply Reply Quote
            • M
              mohui @wwzhao last edited by

              @wwzhao 谢谢,我看到这个方程了,现在我在想如何引进源项。

              W 1 Reply Last reply Reply Quote
              • W
                wwzhao 教授 @mohui last edited by

                @mohui 你可以把源项 IOobject 注册到内存中,再在湍流模型里通过 lookupObject 查找得到源项。

                M 1 Reply Last reply Reply Quote
                • M
                  mohui @wwzhao last edited by

                  @wwzhao 好的,谢谢。现在稍微了解到3.0以后版本的湍流模型的编译较为复杂,需要仔细看看。

                  Y 1 Reply Last reply Reply Quote
                  • zym604
                    zym604 last edited by

                    用pimpleFoam的话可以直接在fvoptions里添加momentum source项。

                    W 1 Reply Last reply Reply Quote
                    • W
                      wwzhao 教授 @zym604 last edited by

                      @zym604 动量方程和湍动能方程是两个不同的概念,fvOptions只能在动量方程中使用。因为只有动量方程的源项中加了fvOptions(U),而湍动能方程没有。

                      zym604 Y 2 Replies Last reply Reply Quote
                      • zym604
                        zym604 @wwzhao last edited by

                        @wwzhao 啊,sorry,没看清题目。

                        1 Reply Last reply Reply Quote
                        • Y
                          yuan_neu @mohui last edited by

                          @mohui 也不复杂吧,比如对于kepsilon,用incompressible turbulence模型,
                          1 在ras目录下直接在kepsilon做更改,
                          2 回到incompressible目录下,编译
                          3 回到求解器,在编译一下
                          ok 还算挺方便的
                          当然也可以添加新的文件 mykepsilon,需要在incompressible目录下修改某子文件夹下的 xxxmodels 添加两句代码,就可以了

                          1 Reply Last reply Reply Quote
                          • Y
                            yuan_neu @wwzhao last edited by

                            @wwzhao 其实很多模型里面已经有ksource项了,当然默认下这项为零,可以直接将这个source给激活就行了。虽然我没用过

                            1 Reply Last reply Reply Quote
                            • Cloud1990
                              Cloud1990 last edited by

                              of5.0里面已经有fvOption了,

                                  tmp<fvScalarMatrix> epsEqn
                                  (
                                      fvm::ddt(alpha, rho, epsilon_)
                                    + fvm::div(alphaRhoPhi, epsilon_)
                                    - fvm::laplacian(alpha*rho*DepsilonEff(), epsilon_)
                                   ==
                                      C1_*alpha()*rho()*G*epsilon_()/k_()
                                    - fvm::SuSp(((2.0/3.0)*C1_ - C3_)*alpha()*rho()*divU, epsilon_)
                                    - fvm::Sp(C2_*alpha()*rho()*epsilon_()/k_(), epsilon_)
                                    + epsilonSource()
                                    + fvOptions(alpha, rho, epsilon_)
                                  );
                              

                              或者复制模仿realizableKE的格式,赋值于epsilonSource() 或者 kSource().

                              1 Reply Last reply Reply Quote
                              • First post
                                Last post

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