Go to the source code of this file.
|
Info<< "Reading field Cs"<< endl;areaScalarField Cs(IOobject("Cs", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE), aMesh);dimensionedScalar Cs0("Cs0", dimMass/dimArea, 1.0);const areaVectorField &R=aMesh.areaCentres();Cs=Cs0 *(1.0+R.component(vector::X)/mag(R));dimensionedScalar Ds("Ds", dimViscosity, 1.0);areaVectorField Us(IOobject("Us", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE), aMesh, dimensionedVector(dimVelocity, Zero));dimensioned< scalar > | Uinf ("Uinf", dimVelocity, 1.0) |
|
| forAll (Us, faceI) |
|
edgeScalarField | phis (IOobject("phis", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE), linearEdgeInterpolate(Us) &aMesh.Le()) |
|
|
| Us = aMesh.faceAreaNormals()*(aMesh.faceAreaNormals() & Us) |
|
◆ Uinf()
Info<< "Reading field Cs"<< endl;areaScalarField Cs( IOobject( "Cs", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE), aMesh);dimensionedScalar Cs0("Cs0", dimMass/dimArea, 1.0);const areaVectorField &R=aMesh.areaCentres();Cs=Cs0 *(1.0+R.component(vector::X)/mag(R));dimensionedScalar Ds("Ds", dimViscosity, 1.0);areaVectorField Us( IOobject( "Us", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE), aMesh, dimensionedVector(dimVelocity, Zero));dimensioned< scalar > Uinf |
( |
"Uinf" |
, |
|
|
dimVelocity |
, |
|
|
1. |
0 |
|
) |
| |
◆ forAll()
◆ phis()
edgeScalarField phis |
( |
IOobject("phis", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE) |
, |
|
|
linearEdgeInterpolate(Us) &aMesh.Le() |
|
|
) |
| |
◆ Us
Us = aMesh.faceAreaNormals()*(aMesh.faceAreaNormals() & Us) |