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

[ChimneyFoam]

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

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

[王金成]

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

[Cr]

在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 你好，非常感谢你的回答，我还有一个问题想请教以下。请问我如果只想使体标量场的一部分值（某一区域的值）随时间变化，其余区域的值还是初始值，请问这种情况能操作吗？

[Cr]

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

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

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

[Cr]

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

[Cr]

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

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

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

[王金成]

@Cr 十分感谢，问题已经解决了