112 lines
2.0 KiB
C
112 lines
2.0 KiB
C
IOobject UmeanHeader
|
|
(
|
|
"Umean",
|
|
runTime.times()[i].name(),
|
|
mesh,
|
|
IOobject::MUST_READ
|
|
);
|
|
|
|
if (!UmeanHeader.headerOk())
|
|
{
|
|
continue;
|
|
}
|
|
|
|
volVectorField Umean
|
|
(
|
|
UmeanHeader,
|
|
mesh
|
|
);
|
|
|
|
volSymmTensorField R
|
|
(
|
|
IOobject
|
|
(
|
|
"R",
|
|
runTime.times()[i].name(),
|
|
mesh,
|
|
IOobject::MUST_READ
|
|
),
|
|
mesh
|
|
);
|
|
volScalarField Rxx(R.component(symmTensor::XX));
|
|
volScalarField Ryy(R.component(symmTensor::YY));
|
|
volScalarField Rzz(R.component(symmTensor::ZZ));
|
|
volScalarField Rxy(R.component(symmTensor::XY));
|
|
|
|
volScalarField pPrime2Mean
|
|
(
|
|
IOobject
|
|
(
|
|
"pPrime2Mean",
|
|
runTime.times()[i].name(),
|
|
mesh,
|
|
IOobject::MUST_READ
|
|
),
|
|
mesh
|
|
);
|
|
|
|
volScalarField epsilonMean
|
|
(
|
|
IOobject
|
|
(
|
|
"epsilonMean",
|
|
runTime.times()[i].name(),
|
|
mesh,
|
|
IOobject::MUST_READ
|
|
),
|
|
mesh
|
|
);
|
|
|
|
/*
|
|
volScalarField nuMean
|
|
(
|
|
IOobject
|
|
(
|
|
"nuMean",
|
|
runTime.times()[i].name(),
|
|
mesh,
|
|
IOobject::MUST_READ
|
|
),
|
|
mesh
|
|
);
|
|
|
|
volScalarField gammaDotMean
|
|
(
|
|
IOobject
|
|
(
|
|
"gammaDotMean",
|
|
runTime.times()[i].name(),
|
|
mesh,
|
|
IOobject::MUST_READ
|
|
),
|
|
mesh
|
|
);
|
|
|
|
volScalarField nuPrime2
|
|
(
|
|
IOobject
|
|
(
|
|
"nuPrime",
|
|
runTime.times()[i].name(),
|
|
mesh,
|
|
IOobject::MUST_READ
|
|
),
|
|
mesh
|
|
);
|
|
volScalarField nuPrime = sqrt(mag(nuPrime2 - sqr(nuMean)));
|
|
|
|
|
|
volScalarField gammaDotPrime2
|
|
(
|
|
IOobject
|
|
(
|
|
"gammaDotPrime",
|
|
runTime.times()[i].name(),
|
|
mesh,
|
|
IOobject::MUST_READ
|
|
),
|
|
mesh
|
|
);
|
|
volScalarField gammaDotPrime = sqrt(mag(gammaDotPrime2 -sqr(gammaDotMean)));
|
|
*/
|