Navigation

    CFD中文网

    CFD中文网

    • Login
    • Search
    • 最新

    不同的物性参数材料如何制定?

    OpenFOAM
    4
    12
    5787
    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.
    • piteqiu
      piteqiu last edited by

      计算固体传热问题,区域内有不同物性参数的材料,比如:铜、铝、铁等。
      0_1499739265867_Image 1.png
      上图为剖分截面图,下图为三维剖分图:
      0_1499739441437_Image 2.png
      setFields貌似针对气体或液体质量分数发生变化时设置的吧?适用于固体传热吗?
      我这种情况如何分区域设置不同的物性参数?

      程 1 Reply Last reply Reply Quote
      • piteqiu
        piteqiu last edited by

        用laplacianFoam求解。
        cellZones中包含4个zone。transportProperties文件只设置了DT(导热系数),且只有一个值。
        在哪设置不同的物性参数比如:密度,比热,导热系数?

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

          固体导热为什么要设置密度和比热容?

          laplacianFoam目前不能多区域求解,只能改个求解器,你把DT项设置为AUTO_READ,然后通过setFields设置DT的值。

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

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

            瞬态计算,需要设置密度和比热。
            laplacianFoam中植入的瞬态热传导方程为:
            0_1499751918842_Image 1.png
            不能分区设置密度和比热。
            目前有多区域求解温度场的求解器吗?

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

              setFilelds只能设置矩形区域吗?
              cellZones文件中已经区分好区域了,能根据cellZones的不同elem设置不同的物性参数吗?

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

                不能分区设置密度和比热。

                DT是热到率,你帖的方程并不存在密度。

                目前有多区域求解温度场的求解器吗?

                chtMultiRegionFoam

                setFilelds只能设置矩形区域吗?

                还可以设置别的简单形状,不过从你帖的图片看起来比较简单?

                cellZones文件中已经区分好区域了,能根据cellZones的不同elem设置不同的物性参数吗?

                可以通过fvOptions,不过需要改代码,还不如多用几个矩形区域包围 :cheeky:

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

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

                  1. 上图贴的方程是laplacianFoam中的方程。
                    我的方程是:0_1499767723704_Image 1.png
                    由于是固体传热,忽略速度,简化为:
                    0_1499767752030_Image 2.png
                    有内热源,瞬态问题,不同区域密度,比热,导热系数不一样。

                  2. 上图是我做的简易模型,实际模型几何比较复杂。
                    照您这样说,复杂几何setFields基本不能解决问题?

                  3. chtMultiRegionFoam可以解决我的方程和不同区域材料设置吗?用这个求解器还是需要用fvOptions设置不同区域材料吗?

                  Wayne 1 Reply Last reply Reply Quote
                  • Wayne
                    Wayne @piteqiu last edited by

                    @piteqiu 把这些物性参数都弄成场,然后随便咋赋值。。

                    Blog: http://blog.sina.com.cn/multiphyzks
                    RG:https://www.researchgate.net/profile/Yan_Wang154

                    piteqiu 1 Reply Last reply Reply Quote
                    • piteqiu
                      piteqiu @Wayne last edited by

                      @Wayne 这个如何弄?物性参数当作求解变量?

                      Wayne 1 Reply Last reply Reply Quote
                      • Wayne
                        Wayne @piteqiu last edited by

                        @piteqiu

                            Info<< "Reading field S\n" << endl;
                            volScalarField S  //1.48e-5 kg/s
                            (
                                IOobject
                                (
                                    "S",
                                    runTime.timeName(),
                                    mesh,
                                    IOobject::NO_READ,
                                    IOobject::AUTO_WRITE
                                ),
                                mesh,
                                dimensionedScalar("S",dimensionSet(1 ,-3, -1 ,0 ,0 ,0 ,0),0.0)
                            );
                            forAll(S,i)
                            {
                                S[i]=0;
                            }
                            S[mesh.findCell(point(0,0,0))]=Q.value()/mesh.V()[mesh.findCell(point(0,0,0))];
                        

                        Blog: http://blog.sina.com.cn/multiphyzks
                        RG:https://www.researchgate.net/profile/Yan_Wang154

                        piteqiu 1 Reply Last reply Reply Quote
                        • piteqiu
                          piteqiu @Wayne last edited by

                          @Wayne 好的,多谢,我试试。

                          1 Reply Last reply Reply Quote
                          • 程
                            程迪 @piteqiu last edited by

                            @piteqiu
                            还有就是那个扩散物性用inverse插值比线性插值准确。

                            github: chengdi123000
                            网站:chengdi123000.github.io
                            本作品采用知识共享署名-非商业性使用-相同方式共享 4.0 国际许可协议进行许可。

                            1 Reply Last reply Reply Quote
                            • First post
                              Last post

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