在OpenFOAM中如何判断一个单元是否为边界单元？

[AppleKiller]

各位道友好！

关于这个问题，我最初的想法是遍历该单元的所有面，然后判断这些面是否在边界上。请问有没有其他的方法呢？（任何方法都行）

[AppleKiller]

@AppleKiller 在 在OpenFOAM中如何判断一个单元是否为边界单元？ 中说：

各位道友好！

关于这个问题，我最初的想法是遍历该单元的所有面，然后判断这些面是否在边界上。请问有没有其他的方法呢？（任何方法都行）

这是这个方法的代码，有没有其他比较快速的方法呢？

forAll (mesh().cells()[cellID], faceI)
{
	const label patchi = mesh().boundaryMesh().which(mesh().cells()[cellID][faceI]);
	
	if
	(
	   patchi > 0
	&& !isA<processorPolyPatch>(mesh().boundaryMesh()[patchi])									
	)
	{
		Info << "Find it" << endl;
	}
}

[学流体的小明]

可以从边界上的face返回这个face属于的cell

const Foam::fvBoundaryMesh &procBoundary = mesh.boundary();
    forAll(procBoundary, patch)
        {
        forAll(procBoundary[patch], facei)
            {
            const label bcell = boundaryMesh[patch].faceCells()[facei];
            }
        }

你可以得到所有的边界单元编号，然后做你其它向做的事就应该比较方便了