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

BUG: potentialFoam/cylinder: indexing into non-existing patch in parallel

Browse Source
This commit is contained in:
mattijs
2011-10-20 16:11:45 +01:00
parent f1c7967316
commit fcebb1e6e7
2 changed files with 24 additions and 4 deletions
Download Patch File Download Diff File Expand all files Collapse all files

4
tutorials/basic/potentialFoam/cylinder/0.org/U
View File

@ -37,8 +37,8 @@ boundaryField
left left
{ {
type fixedValue; type uniformFixedValue;
value uniform (1 0 0); uniformValue (1 0 0);
} }
cylinder cylinder

24
tutorials/basic/potentialFoam/cylinder/system/controlDict
View File

@ -62,19 +62,39 @@ functions
Info<< "Reading inlet velocity uInfX\n" << endl; Info<< "Reading inlet velocity uInfX\n" << endl;
scalar ULeft = 0.0;
label leftI = mesh().boundaryMesh().findPatchID("left");
const fvPatchVectorField& fvp = U.boundaryField()[leftI];
if (fvp.size())
{
ULeft = fvp[0].x();
}
reduce(ULeft, maxOp<scalar>());
dimensionedScalar uInfX dimensionedScalar uInfX
( (
"uInfx", "uInfx",
dimensionSet(0, 1, -1, 0, 0), dimensionSet(0, 1, -1, 0, 0),
U.boundaryField()[3][0].x() ULeft
); );
Info << "U at inlet = " << uInfX.value() << " m/s" << endl; Info << "U at inlet = " << uInfX.value() << " m/s" << endl;
scalar magCylinder = 0.0;
label cylI = mesh().boundaryMesh().findPatchID("cylinder");
const fvPatchVectorField& cylFvp = mesh().C().boundaryField()[cylI];
if (cylFvp.size())
{
magCylinder = mag(cylFvp[0]);
}
reduce(magCylinder, maxOp<scalar>());
dimensionedScalar radius dimensionedScalar radius
( (
"radius", "radius",
dimensionSet(0, 1, 0, 0, 0), dimensionSet(0, 1, 0, 0, 0),
mag(U.mesh().boundary()[4].Cf()[0]) magCylinder
); );
Info << "Cylinder radius = " << radius.value() << " m" << endl; Info << "Cylinder radius = " << radius.value() << " m" << endl;