简介
通常在流体计算中需要对某个特定区域进行处理（比如添加源项,可参考这篇文章OpenFOAM编程：VOF法与多孔介质模型相结合)，这是就需要用到cellZone.
通常有两种产生cellZone的方式：
（1）从其他划分网格的工具中导入时，自然会对不同实体用cellZone区分。可参考这篇文章将fluentMeshing网格转换为openFoam网格
（2）在已有的网格上使用topoSet命令生成cellZone
本文只关心第二种方式
设置cellZone
本文使用icoFoam的cavity案例
为了告诉openFoam哪个区域要设置成cellZone，需要在system目录下编写topoSetDict文件
/*--------------------------------*- C++ -*----------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Version: 8
\\/ M anipulation |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
object topoSetDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
actions
(
{
name g; // 单元集合命名为g
type cellSet;
action new;
source cylinderToCell; // 圆柱区域
sourceInfo
{
p1 (0.0 0.0 0); // 圆柱体上表面的中心点
p2 (0.05 0 0); // 圆柱体下表面的中心点
radius 0.03; // 圆柱半径
}
}
{
name grain; // 区域命名grain
type cellZoneSet;
action new;
source setToCellZone;
sourceInfo
{
set g; // 给定source信息,来自单元集合g
}
}
);
// ************************************************************************* //
首先生成一个单元集合g，g由一个圆柱体确定。然后根据g再生成名称为grain的cellZone
接下来运行topoSet命令即可获得名称为grain的cellZone
使用paraFoam命令检查一下网格
下图是划分前的网格
下图是划分后的网格
蓝色区域就是名称为grain的cellZone
使用cellZone
cellZone本质上是一个列表，每个元素记录了cell的序号。所以可以通过访问cellZone来设置给定区域的物理量。比如为了将grain处的压力p设置为30，可以将如下代码插入到icoFoam求解器当中
labelUList cellListOfGrain=mesh.cellZones()["grain"];
forAll(cellListOfGrain,i)
{
p[cellListOfGrain[i]]=30;
}
p.write();
运行编译后的icoFoam，可以得到如下结果
