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 十分感谢,问题已经解决了
论坛登录问题反馈可联系 li.dy@dyfluid.com