Skip to content
Snippets Groups Projects
Commit afa64919 authored by Henry Weller's avatar Henry Weller
Browse files

potentialFoam: simplify the Phi BCs to use only fixedValue and zeroGradient by default

Resolves bug-report http://bugs.openfoam.org/view.php?id=2129
parent 6e10f3a3
Branches
Tags
1 merge request!60Merge foundation
......@@ -40,7 +40,7 @@ word pName("p");
// Update name of the pressure field from the command-line option
args.optionReadIfPresent("pName", pName);
// Infer the pressure BCs from the velocity BCs
// Infer the pressure BCs from the velocity
wordList pBCTypes
(
U.boundaryField().size(),
......@@ -71,6 +71,21 @@ volScalarField p
pBCTypes
);
// Infer the velocity potential BCs from the pressure
wordList PhiBCTypes
(
p.boundaryField().size(),
zeroGradientFvPatchScalarField::typeName
);
forAll(p.boundaryField(), patchi)
{
if (p.boundaryField()[patchi].fixesValue())
{
PhiBCTypes[patchi] = fixedValueFvPatchScalarField::typeName;
}
}
Info<< "Constructing velocity potential field Phi\n" << endl;
volScalarField Phi
(
......@@ -84,7 +99,7 @@ volScalarField Phi
),
mesh,
dimensionedScalar("Phi", dimLength*dimVelocity, 0),
p.boundaryField().types()
PhiBCTypes
);
label PhiRefCell = 0;
......
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment