scalarField、surfaceScalarField与slicedSurfaceScalarField

Lidong Cheng

已于 2023-04-28 16:45:10 修改

阅读量785

点赞数

分类专栏： OpenFOAM 文章标签：经验分享

于 2022-04-27 21:04:11 首次发布

本文链接：https://blog.csdn.net/Jerry_CLD/article/details/124460318

版权

OpenFOAM 专栏收录该内容

11 篇文章 3 订阅

订阅专栏

看到一段与这三个类相关的经典代码，记录备查。

来源：https://www.cfd-online.com/Forums/openfoam

surfaceScalarField D(fvc::interpolate(D_cell));
scalarField allLambda(mesh.nFaces(), 1.0);

slicedSurfaceScalarField lambda
(
    IOobject
    (
        "lambda",
        mesh.time().timeName(),
        mesh,
        IOobject::NO_READ,
        IOobject::NO_WRITE,
        false
    ),
    mesh,
    dimless,
    allLambda,
    false   // Use slices for the couples
);


scalarField& lambdaIf = lambda;
surfaceScalarField::Boundary& lambdaBf = lambda.boundaryFieldRef();

forAll(lambdaIf, facei)
{
    if(P[facei] < epsilon) {lambdaIf[facei] = 0;}
}

forAll(lambdaBf, patchi)
{
    scalarField& lambdaPf = lambdaBf[patchi];
    forAll(lambdaPf, pFacei)
    {
        if(P[pFacei] < epsilon) {lambdaPf[pFacei] = 0;}
    }
}

syncTools::syncFaceList(mesh,allLambda,minEqOp<scalar>());

D*=lambda;

scalarField 是标量场；
surfaceScalarField 为几何场量，包含 internalField 和 boundaryField；
slicedSurfaceScalarField 生成 scalarField 的引用，其数据结构与surfaceScalarField相同，主要用于自动处理 coupled 边界。