zhyang-dev/openfoam
1
0
Fork 0
mirror of https://develop.openfoam.com/Development/openfoam.git synced 2025-11-28 03:28:01 +00:00
Code Issues Actions Packages Projects Releases Wiki Activity

Merge branch 'master' of /home/noisy2/OpenFOAM/OpenFOAM-dev

Browse Source
This commit is contained in:
mattijs
2008-07-18 12:59:30 +01:00
parent 537f11310f 3c07d20652
commit d6c87e58c6
3 changed files with 29 additions and 30 deletions
Download Patch File Download Diff File Expand all files Collapse all files

12
applications/utilities/postProcessing/miscellaneous/postChannel/calculateFields.H
View File

@ -1,9 +1,9 @@
/*
volTensorField gradU = fvc::grad(U);
volSymmTensorField D = symm(fvc::grad(U));
volTensorField Dprim = symm(fvc::grad(U - Umean));
volTensorField Dprim = symm(fvc::grad(U - UMean));
volScalarField prod = -((U - Umean)*(U - Umean)) && D;
volScalarField prod = -((U - UMean)*(U - UMean)) && D;
*/
/*
@ -20,7 +20,7 @@
mesh,
dimensionSet(0, 1, -1, 0, 0)
);
txx =sqrt(Txx - (Umeanx*Umeanx));
txx =sqrt(Txx - (UMeanx*UMeanx));
txx.write();
volScalarField tyy
@ -36,7 +36,7 @@
mesh,
dimensionSet(0, 1, -1, 0, 0)
);
tyy = sqrt(Tyy - (Umeany*Umeany));
tyy = sqrt(Tyy - (UMeany*UMeany));
tyy.write();
volScalarField tzz
@ -52,7 +52,7 @@
mesh,
dimensionSet(0, 1, -1, 0, 0)
);
tzz = sqrt(Tzz - (Umeanz*Umeanz));
tzz = sqrt(Tzz - (UMeanz*UMeanz));
tzz.write();
volScalarField txy
@ -68,6 +68,6 @@
mesh,
dimensionSet(0, 2, -2, 0, 0)
);
txy = Txy - (Umeanx*Umeany);
txy = Txy - (UMeanx*UMeany);
txy.write();
*/

32
applications/utilities/postProcessing/miscellaneous/postChannel/collapse.H
View File

@ -20,9 +20,9 @@
scalarField pPrime2MeanValues = channelIndexing.collapse(pPrime2Mean);
/*
scalarField epsilonValues = channelIndexing.collapse(epsilonMean);
/*
scalarField nuMeanValues = channelIndexing.collapse(nuMean);
scalarField nuPrimeValues = channelIndexing.collapse(nuPrime);
@ -34,27 +34,25 @@
scalarField vrmsValues = sqrt(mag(RyyValues));
scalarField wrmsValues = sqrt(mag(RzzValues));
scalarField kValues
= 0.5*(sqr(urmsValues) + sqr(vrmsValues) + sqr(wrmsValues));
scalarField kValues =
0.5*(sqr(urmsValues) + sqr(vrmsValues) + sqr(wrmsValues));
scalarField y = channelIndexing.y(mesh.C());
makeGraph(y, UMeanXvalues, "Uf", R.path(), gFormat);
makeGraph(y, urmsValues, "u", R.path(), gFormat);
makeGraph(y, vrmsValues, "v", R.path(), gFormat);
makeGraph(y, wrmsValues, "w", R.path(), gFormat);
makeGraph(y, RxyValues, "uv", R.path(), gFormat);
makeGraph(y, kValues, "k", R.path(), gFormat);
makeGraph(y, UMeanXvalues, "Uf", UMean.path(), gFormat);
makeGraph(y, urmsValues, "u", UMean.path(), gFormat);
makeGraph(y, vrmsValues, "v", UMean.path(), gFormat);
makeGraph(y, wrmsValues, "w", UMean.path(), gFormat);
makeGraph(y, RxyValues, "uv", UMean.path(), gFormat);
makeGraph(y, kValues, "k", UMean.path(), gFormat);
makeGraph(y, pPrime2MeanValues, "pPrime2Mean", R.path(), gFormat);
makeGraph(y, epsilonValues, "epsilon", R.path(), gFormat);
makeGraph(y, pPrime2MeanValues, "pPrime2Mean", UMean.path(), gFormat);
/*
makeGraph(y, nuMeanValues, "nu", R.path(), gFormat);
makeGraph(y, nuPrimeValues, "nuPrime", R.path(), gFormat);
makeGraph(y, gammaDotMeanValues, "gammaDot", R.path(), gFormat);
makeGraph(y, gammaDotPrimeValues, "gammaDotPrime", R.path(), gFormat);
makeGraph(y, epsilonValues, "epsilon", UMean.path(), gFormat);
makeGraph(y, nuMeanValues, "nu", UMean.path(), gFormat);
makeGraph(y, nuPrimeValues, "nuPrime", UMean.path(), gFormat);
makeGraph(y, gammaDotMeanValues, "gammaDot", UMean.path(), gFormat);
makeGraph(y, gammaDotPrimeValues, "gammaDotPrime", UMean.path(), gFormat);
*/

15
applications/utilities/postProcessing/miscellaneous/postChannel/readFields.H
View File

@ -8,6 +8,7 @@
if (!UMeanHeader.headerOk())
{
Info<< " No UMean field" << endl;
continue;
}
@ -17,21 +18,21 @@
mesh
);
volSymmTensorField R
volSymmTensorField UPrime2Mean
(
IOobject
(
"R",
"UPrime2Mean",
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 Rxx(UPrime2Mean.component(symmTensor::XX));
volScalarField Ryy(UPrime2Mean.component(symmTensor::YY));
volScalarField Rzz(UPrime2Mean.component(symmTensor::ZZ));
volScalarField Rxy(UPrime2Mean.component(symmTensor::XY));
volScalarField pPrime2Mean
(
@ -45,6 +46,7 @@
mesh
);
/*
volScalarField epsilonMean
(
IOobject
@ -57,7 +59,6 @@
mesh
);
/*
volScalarField nuMean
(
IOobject