计算Patch的面积
// 计算某个patch的面积 patchID表示面的边界条件标识符
label patchID = mesh.boundaryMesh().findPatchID("inlet");
// 创建一个polyPatch为了循环
const polyPatch& myPatch = mesh.boundaryMesh()[patchID];
// 初始化patchArea
scalar patchArea = 0.0;
// 循环polyPatch的所有面,添加其表面大小
forAll(myPatch, faceI)
{
patchArea += mesh.magSf().boundaryField()[patchID][faceI]; // 获取指定面的面积
}
获取Patch的名字
Info << "patch name is " << U.boundaryField()[patch].patch().name() << endl;
网格相关量
const scalarField& V = mesh.V(); //网格体积
const surfaceVectorField& Sf = mesh.Sf(); //网格面矢量

最低0.47元/天 解锁文章
8346

被折叠的 条评论
为什么被折叠?



