Navigation

    CFD中文网

    CFD中文网

    • Login
    • Search
    • 最新

    关于在rhoPimpleFoam里面修改EEqn和状态方程的问题

    OpenFOAM
    2
    9
    3649
    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.
    • 1
      100yearsalone last edited by

      大家好,:happy:

      我在使用rhoPimpleFoam来计算空化,在其中植入自己的代码,目前第一步是在改了这个压缩求解器中的状态方程.原来的设置是perfectGas就是理想气体状态方程来计算的,我把它改成了一个这样的方程式:
      在pEqn的头部
      0_1505579116624_2017-09-17 01:24:33屏幕截图.png
      和尾部的地方把rho=thermo.rho()的地方换成了这个.
      0_1505579278636_2017-09-17 01:24:47屏幕截图.png
      my_Eqn_of_State的代码,就是把通过气相和液相求出综合的密度,其中RHOBF和RHOL2F是由自定义数据来实现的,根据T和P的不同决定密度.
      0_1505579537178_2017-09-17 01:32:04屏幕截图.png
      计算后发现是可以计算的,类那边也看过但觉得对自己改来说太复杂了就用这样的方法去改,不知道方法正确与否希望大家给点意见.
      下图是整个主题的代码,加了两行TEqn和ElEqn
      0_1505579382462_2017-09-17 01:28:53屏幕截图.png
      因为一开始计算的是综合的密度,所以得出的也是混合体的压力和速度什么的,由混合物的速度可以求出液体的速度.其中混合体的温度T是通过焓方程求出来的.我老师让我把液体的速度再代入焓方程求液体的温度.
      所以我在pimple的循环之后加了一个ElEqn的来像计算出液体的温度.
      0_1505582464607_2017-09-17 02:20:57屏幕截图.png
      因为现在只是想知道方法,所以直接让混合体速度等于液体速度,这一步根据计算结果来看没有问题.然后我把EEqn文档中的数值符号全部加了l在后面,volScalarField也全部在creatField里面声明了.但计算结果来看,好像有问题,而且也没有输出Tl的值.
      请问一下大家我对于这个求解器的修改有什么不完美的地方吗?
      还望大家赐教.

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

        所以你的问题是没有输出TI?你是说没有写入还是在终端没有输出?没有写入的话试试在求解器后面增加TI.write()函数试试,没有输出的话可能表示你的Info行没有编译进去,
        都不是大问题。目前你可以这样修改,简单点。

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

        1 1 Reply Last reply Reply Quote
        • 1
          100yearsalone @李东岳 last edited by

          @李东岳 谢谢:laughing: ,根据你的指导的确是把Tl成功的输出了。
          0_1505636392910_Screenshot from 2017-09-17 17:13:35.png 上面这个是添加在TlEqn下面的代码。
          但是输出到Tl的时间步文件中的各个格子的值似乎没有被计算,而是直接被输出了。
          0_1505636621632_Screenshot from 2017-09-17 17:21:29.png 我查了下时间步的phi和自己定义的phil输出的也不一样。
          0_1505636661205_Screenshot from 2017-09-17 17:21:49.png下面是phil
          0_1505636665993_Screenshot from 2017-09-17 17:22:00.png 而且输出的文档也有些异样
          0_1505636834325_Screenshot from 2017-09-17 17:26:36.png 在运算是控制台输出时最后一步TlEqn的计算特别慢,好像卡了一下,从这个文档来看,最后的iterations也到了1000啊,好像是有些问题。

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

            各个格子的值似乎没有被计算,而是直接被输出了。

            我觉得你添加的代码没有编译进去。

            你做什么的?听起来好像是做别的方向模拟转过来的?:confused:

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

            1 2 Replies Last reply Reply Quote
            • 1
              100yearsalone @李东岳 last edited by

              @李东岳
              现在是在做空化

              1 Reply Last reply Reply Quote
              • 1
                100yearsalone @李东岳 last edited by

                @李东岳
                0_1506036877347_image.png
                老师,请问下e是在哪个地方可以获取?

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

                  e是什么?能量么?能量方程需要求的啊。EEqn.H。

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

                  1 1 Reply Last reply Reply Quote
                  • 1
                    100yearsalone @李东岳 last edited by

                    @李东岳
                    这几天我一直在试着这个方程,一直都是焓he无法收敛,最后是发散。
                    所以我改了下方程,让他变成了能量的方程。
                    0_1506411167066_UC截图20170926163214.png 其中状态方程是自己定义的,如前面所说
                    替代文字把thermo.rho()换成了这个图里的状态方程,rhoG和rhoL的函数都是自己定义的,这样上述四个式子可以求出内部能量e,EEqn我是这样改了下0_1506411964389_Screenshot from 2017-09-26 16-42-48.png 原来是thermo.he()的,直接使用理想气体状态方程时是没问题的,但一用这种方式就无法收敛以及发散。
                    0_1506412240205_Screenshot from 2017-09-26 16-47-51.png
                    我觉得我是气体状态方程的地方出现了问题,但具体有什么问题我不是很清楚。

                    1 Reply Last reply Reply Quote
                    • 1
                      100yearsalone last edited by

                      自己挖的坑自己填回去。。。。。。。。。。。。。

                      我今天发现是我自己定义的TELF函数出了问题,所以计算出来的结果有问题。
                      但是这种方法加入状态方程我实践过的确是可行的。
                      就是换了以后,连续性误差的地方会变的非常怪,查了continuityErr这个文档,发现里面计算连续性误差还是使用rho-thermo.rho所以会出现误差值十分大,
                      目前还在想着这么去改thermo下的东西,因为自己不能把黏性系数和普朗数设成固定值,
                      怎样把这些thermo文档下设置函数还在研究中。

                      1 Reply Last reply Reply Quote
                      • First post
                        Last post

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