1 Info<<
"Reading velocity field U\n" <<
endl;
33 U.correctBoundaryConditions();
48 U.boundaryField().size(),
49 fixedValueFvPatchScalarField::typeName
54 if (
U.boundaryField()[patchi].fixesValue())
56 pBCTypes[patchi] = zeroGradientFvPatchScalarField::typeName;
60 Info<<
"Constructing pressure field " << pName <<
nl <<
endl;
68 IOobject::READ_IF_PRESENT,
79 p.boundaryField().size(),
80 zeroGradientFvPatchScalarField::typeName
83 forAll(
p.boundaryField(), patchi)
85 if (
p.boundaryField()[patchi].fixesValue())
87 PhiBCTypes[patchi] = fixedValueFvPatchScalarField::typeName;
91 Info<<
"Constructing velocity potential field Phi\n" <<
endl;
99 IOobject::READ_IF_PRESENT,
107 label PhiRefCell = 0;
108 scalar PhiRefValue = 0;
116 mesh.setFluxRequired(Phi.name());
118 #include "createMRF.H" forAll(U.boundaryField(), patchi)
tmp< surfaceScalarField > flux(const volVectorField &vvf)
Return the face-flux field obtained from the given volVectorField.
dimensionedSymmTensor sqr(const dimensionedVector &dv)
constexpr char nl
The newline '\n' character (0x0a)
dimensioned< vector > dimensionedVector
Dimensioned vector obtained from generic dimensioned type.
Ostream & endl(Ostream &os)
Add newline and flush stream.
const dictionary & potentialFlow(mesh.solutionDict().subDict("potentialFlow"))
setRefCell(p, pimple.dict(), pRefCell, pRefValue)
T getOrDefault(const word &optName, const T &deflt) const
Get a value from the named option if present, or return default.
GeometricField< vector, fvPatchField, volMesh > volVectorField
GeometricField< scalar, fvPatchField, volMesh > volScalarField
wordList pBCTypes(U.boundaryField().size(), fixedValueFvPatchScalarField::typeName)
volVectorField U(IOobject("U", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE), mesh, dimensionedVector(dimVelocity, Zero))
List< word > wordList
List of word.
const dimensionSet dimLength(0, 1, 0, 0, 0, 0, 0)
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
messageStream Info
Information stream (stdout output on master, null elsewhere)
GeometricField< scalar, fvsPatchField, surfaceMesh > surfaceScalarField
Foam::argList args(argc, argv)
bool found(const word &optName) const
Return true if the named option is found.
static constexpr const zero Zero
Global zero (0)
const dimensionSet dimVelocity