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



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

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



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



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

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

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

    permeability =timeVaryingValue.value(runTime.timeOutputValue());
    


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



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

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

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



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



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

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

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



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


Log in to reply
 

CFD中文网 2016 - 2020 | 京ICP备15017992号-2