zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

TUT: add coded example for sampling cellZoneId (#1383)

Browse Source
This commit is contained in:
Mark Olesen
2019-08-21 12:13:51 +02:00
committed by Andrew Heather
parent c0fce5c762
commit 2ca4792734
2 changed files with 128 additions and 0 deletions
Download Patch File Download Diff File Expand all files Collapse all files

5
tutorials/compressible/rhoPorousSimpleFoam/angledDuct/implicit/system/controlDict
View File

@ -47,5 +47,10 @@ graphFormat raw;
runTimeModifiable true;
functions
{
#include "sampling"
}
// ************************************************************************* //

123
tutorials/compressible/rhoPorousSimpleFoam/angledDuct/implicit/system/sampling Normal file
View File

@ -0,0 +1,123 @@
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: v1906 |
| \\ / A nd | Web: www.OpenFOAM.com |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
cellZoneID
{
type coded;
libs ("utilityFunctionObjects");
name cellZoneID;
executeControl writeTime;
localCode
#{
static const word fieldName("cellZoneID");
// Create and populate a "cellZoneID" volScalarField for cellZones.
// Return nullptr if there are not any such zones.
//
// Could be improved for mesh motion etc
// but is fairly low overhead anyhow
volScalarField* getZoneField(const fvMesh& mesh)
{
volScalarField* volZonePtr =
mesh.getObjectPtr<volScalarField>(fieldName);
const cellZoneMesh& zones = mesh.cellZones();
if (!zones.empty())
{
if (!volZonePtr)
{
volZonePtr = new volScalarField
(
IOobject
(
fieldName,
mesh.time().timeName(),
mesh,
IOobject::NO_READ,
IOobject::NO_WRITE,
true // Register
),
mesh,
dimless,
zeroGradientFvPatchField<scalar>::typeName
);
volZonePtr->store();
Info<< "Creating " << fieldName
<< " field for postProcessing" << nl;
}
}
else
{
Info<< "No cellZones - not creating " << fieldName
<< " field for postProcessing" << nl;
}
if (volZonePtr)
{
// Copy zone id as scalar.
// For consistency, do this whenever the volField exists,
// even if there are no zones.
auto& field = *volZonePtr;
auto& internal = field.primitiveFieldRef();
internal = scalar(-1);
for (const cellZone& zn : zones)
{
UIndirectList<scalar>(internal, zn) = scalar(zn.index());
}
field.correctBoundaryConditions();
}
return volZonePtr;
}
#};
codeExecute
#{
// Don't need the return value, just the contents on the registry.
(void) getZoneField(mesh());
#};
}
plane
{
type surfaces;
libs ("sampling");
writeControl writeTime;
surfaceFormat vtk;
fields ( cellZoneID );
surfaces
(
zNormal
{
type cuttingPlane;
planeType pointAndNormal;
pointAndNormalDict
{
point (0 0 0);
normal (0 0 1);
}
interpolate false;
}
);
}
// ************************************************************************* //