Navigation

    CFD中文网

    CFD中文网

    • Login
    • Search
    • 最新

    OpenFOAM中随时间变化的属性设置

    OpenFOAM
    3
    8
    1596
    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.
    • C
      ChimneyFoam last edited by

      背景:孔隙流求解器中有一个volScalarField permeability类型的量表示孔隙介质的渗透率,一般都是在0文件夹下给定permeability的边界条件和初始条件(渗透率并不需要求解,所以边界条件其实没有什么用,就是个占位符)。如果我想让这个permeability随时间变化(比如k=k0+t*10^-20),怎么实现?

      随时间变化的边界条件好设置,就是不知道如何设置随时间变化的场值?哪位大神有这个经验,求赐教,不胜感激!

      1 Reply Last reply Reply Quote
      • 王
        王金成 last edited by

        请问你解决了吗,我也遇到了相同的问题。

        1 Reply Last reply Reply Quote
        • C
          Cr last edited by

          在createFields.H中声明一个多项式:

          scalar a[8] = {k0, 1e-20, 0, 0, 0, 0, 0, 0};
          Polynomial<8> timeVaryingValue(a);
          

          在while计算循环中的时间步进后更新该场:

          permeability =timeVaryingValue.value(runTime.timeOutputValue());
          
          王 1 Reply Last reply Reply Quote
          • 王
            王金成 @Cr last edited by

            @Cr 你好,非常感谢你的回答,我还有一个问题想请教以下。请问我如果只想使体标量场的一部分值(某一区域的值)随时间变化,其余区域的值还是初始值,请问这种情况能操作吗?

            C 1 Reply Last reply Reply Quote
            • C
              Cr @王金成 last edited by

              @王金成
              topoSet命令可以可以提取计算域内的指定区域,执行命令后在constant/polyMesh会生成cellZoned的字典文件。
              在求解器的主函数中用下面的语句,可以针对这个区域遍历物理场:

              const cellZone& certainZone = mesh.lookupObject<cellZone>("certainZone");
              forAll(certainZone.zoneMesh(), celli)
              {
                  p[celli] = 0.0;  //将压力场的certainZone区域置0
              }
              

              这只是我的想法,并且编译通过了,但我没有实际应用过。

              C 1 Reply Last reply Reply Quote
              • C
                Cr @Cr last edited by

                @Cr 这里面如果cellZone区域和整体mesh区域的cell编号不一致的的话,就得想别的办法来处理了。

                1 Reply Last reply Reply Quote
                • C
                  Cr last edited by

                  更新一下,偶然看到的,针对物理场特定网格内的值处理:

                  cellSet selectedCells(mesh, cellSetName);
                  labelList cells = selectedCells.toc();
                  forAll(cells, i)
                  {
                  p[i] = 0.0
                  }
                  

                  其中cellSetName是用topoSet命令提取的区域名字。

                  王 1 Reply Last reply Reply Quote
                  • 王
                    王金成 @Cr last edited by

                    @Cr 十分感谢,问题已经解决了:xinxin:

                    1 Reply Last reply Reply Quote
                    • First post
                      Last post

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