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 meanc = 1 - b.weightedAverage(mesh.V()).value();
7 
8  if (Pstream::master())
9  {
10  Info<< "Mean pressure:" << meanp << endl;
11  Info<< "Mean temperature:" << meanT << endl;
12  Info<< "Mean u':" << meanup << endl;
13 
14  logSummaryFile()
15  << runTime.theta() << tab
16  << meanp << tab
17  << meanT << tab
18  << meanup << tab
19  << meanc
20  << endl;
21  }
22 }
Foam::tab
static const char tab
Definition: Ostream.H:259
runTime
engineTime & runTime
Definition: createEngineTime.H:13
Foam::sqrt
dimensionedScalar sqrt(const dimensionedScalar &ds)
Definition: dimensionedScalar.C:142
Foam::endl
Ostream & endl(Ostream &os)
Add newline and flush stream.
Definition: Ostream.H:251
thermo
rhoReactionThermo & thermo
Definition: createFields.H:28
mesh
dynamicFvMesh & mesh
Definition: createDynamicFvMesh.H:18
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(dimless, 1.0));surfaceScalarField neg(IOobject("neg", runTime.timeName(), mesh), mesh, dimensionedScalar(dimless, -1.0));surfaceScalarField phi("phi", fvc::flux(rhoU));Info<< "Creating turbulence model\"<< endl;autoPtr< compressible::momentumTransportModel > turbulence(compressible::momentumTransportModel::New(rho, U, phi, thermo))
Definition: createFields.H:94
Foam::Info
messageStream Info
p
volScalarField & p
Definition: createFieldRefs.H:12