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());
-
@王金成
topoSet
命令可以可以提取计算域内的指定区域,执行命令后在constant/polyMesh
会生成cellZoned
的字典文件。
在求解器的主函数中用下面的语句,可以针对这个区域遍历物理场:const cellZone& certainZone = mesh.lookupObject<cellZone>("certainZone"); forAll(certainZone.zoneMesh(), celli) { p[celli] = 0.0; //将压力场的certainZone区域置0 }
这只是我的想法,并且编译通过了,但我没有实际应用过。
-
更新一下,偶然看到的,针对物理场特定网格内的值处理:
cellSet selectedCells(mesh, cellSetName); labelList cells = selectedCells.toc(); forAll(cells, i) { p[i] = 0.0 }
其中
cellSetName
是用topoSet
命令提取的区域名字。
2020年4月14日 11:15
1/8
2020年10月12日 02:14