使用OpenFOAM读取网格
使用foamNewApp
命令创建一个新App,其中main函数部分已经自动填充了一些必要的代码.
本文的目的是读取文件夹下constant/polyMesh
的网格信息,并输出单元中心坐标和单元间连接关系. 原代码如下
/*---------------------------------------------------------------------------*\
========= |
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
\\ / O peration | Website: https://openfoam.org
\\ / A nd | Copyright (C) 2022 OpenFOAM Foundation
\\/ M anipulation |
-------------------------------------------------------------------------------
License
This file is part of OpenFOAM.
OpenFOAM is free software: you can redistribute it and/or modify it
under the terms of the GNU General Public License as published by
the Free Software Foundation, either version 3 of the License, or
(at your option) any later version.
OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
for more details.
You should have received a copy of the GNU General Public License
along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
Application
volumeRatioCheck
Description
\*---------------------------------------------------------------------------*/
#include "fvCFD.H"//这个是最基本的头文件,包括了openFOAM的基本库
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
int main(int argc, char *argv[])
{
#include "setRootCase.H"//根据argc和argv[]设置案例,是必须的
#include "createTime.H"//记录时间
#include "createMesh.H"//从constant/polyMesh中读取网格文件
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
const label c=mesh.C().size();//获取网格所有中心坐标的数量
/*label相当于标准C语言中的int
mesh是由createMesh.H文件创立的,记录了网格的所有信息
C()函数返回单元中心坐标,为volVectorField类型,可使用数组的方式调用
size()函数返回C()的长度,为一个整数
*/
Info<< nl << "网格数为" << c << nl << endl;//使用Info可以打印各种各样的数据类型,比较实用
const labelListList &neighbour =mesh.cellCells();
/*
cellCells()返回单元连接关系,是一个二维列表
第i行第j列的含义是:与第i个单元所相邻的第j个单元的序号
*/
forAll(neighbour,cellI)//forAll是openFOAM自定义的循环方式
{
List<label> neighboursOfPrime= neighbour[cellI];
Info<< nl <<"第"<< cellI <<"个单元的体积为" << mesh.V()[cellI] << endl;
Info<< "该单元的中心坐标为" << mesh.C()[cellI] << endl;
Info<< "其相邻的单元序号分别为" << endl;
forAll(neighboursOfPrime,i)
{
label index=neighboursOfPrime[i];
Info<< index << endl;
}
}
//nl表示新建一行,以下返回计算时间
Info<< nl << "ExecutionTime = " << runTime.elapsedCpuTime() << " s"
<< " ClockTime = " << runTime.elapsedClockTime() << " s"
<< nl << endl;
Info<< "End\n" << endl;
return 0;
}
// ************************************************************************* //
使用wmake
命令编译后,在命令行中调用编译好的可执行程序.
本文使用了openFOAM教学案例pitzDaily中的网格,部分计算结果如下
Create time
Create mesh for time = 0
网格数为12225
第0个单元的体积为5.04673e-10
该单元的中心坐标为(-0.0198086 0.00015942 0)
其相邻的单元序号分别为
1
18
第1个单元的体积为4.84509e-10
该单元的中心坐标为(-0.0182574 0.00015942 0)
其相邻的单元序号分别为
0
2
19
第2个单元的体积为4.65152e-10
该单元的中心坐标为(-0.0167681 0.00015942 2.78877e-20)
其相邻的单元序号分别为
1
3
20
第3个单元的体积为4.46567e-10
该单元的中心坐标为(-0.0153384 0.00015942 0)
其相邻的单元序号分别为
2
4
21