logSummary.H
Go to the documentation of this file.
1 {
2  const scalar meanP = p.weightedAverage(mesh.V()).value();
3  const scalar meanT = thermo.T().weightedAverage(mesh.V()).value();
4  const scalar meanUp =
5  (sqrt((2.0/3.0)*turbulence->k()))().weightedAverage(mesh.V()).value();
6  const scalar meanB = b.weightedAverage(mesh.V()).value();
7 
8  Info<< "Mean pressure:" << meanP << nl
9  << "Mean temperature:" << meanT << nl
10  << "Mean u':" << meanUp << endl;
11 
12  if (Pstream::master())
13  {
14  logSummaryFile()
15  << runTime.theta() << tab
16  << meanP << tab
17  << meanT << tab
18  << meanUp << tab
19  << 1 - meanB
20  << endl;
21  }
22 }
Foam::nl
constexpr char nl
The newline &#39;\n&#39; character (0x0a)
Definition: Ostream.H:50
runTime
engineTime & runTime
Definition: createEngineTime.H:13
Foam::sqrt
dimensionedScalar sqrt(const dimensionedScalar &ds)
Definition: dimensionedScalar.C:137
Foam::endl
Ostream & endl(Ostream &os)
Add newline and flush stream.
Definition: Ostream.H:531
Foam::tab
constexpr char tab
The tab &#39;\t&#39; character(0x09)
Definition: Ostream.H:49
thermo
psiReactionThermo & thermo
Definition: createFields.H:28
meanT
const scalar meanT
Definition: logSummary.H:3
mesh
dynamicFvMesh & mesh
Definition: createDynamicFvMesh.H:6
logSummaryFile
Info<< "Total cylinder mass: "<< fvc::domainIntegrate(rho).value()<< endl;autoPtr< OFstream > logSummaryFile
Definition: startSummary.H:3
Foam::constant::physicoChemical::b
const dimensionedScalar b
Wien displacement law constant: default SI units: [m.K].
Definition: createFields.H:27
turbulence
Info<< "Reading field U\"<< endl;volVectorField U(IOobject("U", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE), mesh);volScalarField rho(IOobject("rho", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::AUTO_WRITE), thermo.rho());volVectorField rhoU(IOobject("rhoU", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE), rho *U);volScalarField rhoE(IOobject("rhoE", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE), rho *(e+0.5 *magSqr(U)));surfaceScalarField pos(IOobject("pos", runTime.timeName(), mesh), mesh, dimensionedScalar("pos", dimless, 1.0));surfaceScalarField neg(IOobject("neg", runTime.timeName(), mesh), mesh, dimensionedScalar("neg", dimless, -1.0));surfaceScalarField phi("phi", fvc::flux(rhoU));Info<< "Creating turbulence model\"<< endl;autoPtr< compressible::turbulenceModel > turbulence(compressible::turbulenceModel::New(rho, U, phi, thermo))
Definition: createFields.H:94
meanUp
const scalar meanUp
Definition: logSummary.H:4
meanB
const scalar meanB
Definition: logSummary.H:6
Foam::Info
messageStream Info
Information stream (stdout output on master, null elsewhere)
p
volScalarField & p
Definition: createFieldRefs.H:8