Info<< "Constructing pressure field "<< pName<< nl<< endl;volScalarField p(IOobject(pName, runTime.timeName(), mesh, IOobject::READ_IF_PRESENT, IOobject::NO_WRITE, false), mesh, dimensionedScalar(sqr(dimVelocity), Zero), pBCTypes);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(IOobject("Phi", runTime.timeName(), mesh, IOobject::READ_IF_PRESENT, IOobject::NO_WRITE), mesh, dimensionedScalar(dimLength *dimVelocity, Zero), PhiBCTypes);label PhiRefCell=0;scalar PhiRefValue=0;setRefCell(Phi, potentialFlow.dict(), PhiRefCell, PhiRefValue);mesh.setFluxRequired(Phi.name());{ wordHashSet &nonInt=const_cast< wordHashSet & > |
( |
Stencil::New(mesh).nonInterpolatedFields() |
| ) |
|