openfoam寻找与某一边界相邻的一层网格

最新推荐文章于 2023-06-20 21:14:09 发布

zhaopeng162

最新推荐文章于 2023-06-20 21:14:09 发布

阅读量1.2k

点赞数 1

分类专栏： openfoam notes

本文链接：https://blog.csdn.net/zhaopeng162/article/details/100572546

版权

notes 同时被 2 个专栏收录

22 篇文章 5 订阅

订阅专栏

openfoam

13 篇文章 12 订阅

订阅专栏

基本思路

定位到边界的faceID
在owner中根据faceID选择网格cell

volScalarField post
    (
        IOobject
        (
            "post",
            runTime.timeName(),
            mesh,
            IOobject::NO_READ,
            IOobject::AUTO_WRITE
        ),
	mesh,
	dimensionedScalar("0.0",dimensionSet(0,0,0,0,0),0.0)
    );
label patchID = mesh.boundaryMesh().findPatchID("cylinder");
if(patchID == -1){
	FatalErrorIn("can not find patch, please check you patch name")
		<< abort(FatalError);
}
const polyPatch& cPatch = mesh.boundaryMesh()[patchID];

const unallocLabelList& neighbour = mesh.neighbour();
const unallocLabelList& owner = mesh.owner();
forAll(cPatch,faceI){
	label faceId = cPatch.start() + faceI;
	label cellId = owner[faceId];
	post[cellId]=1;
	Info<<"cellId = "<<cellId<<endl;
}
post.write();