logSummary.H
Go to the documentation of this file.
1 {
2  const scalar meanp = p.weightedAverage(mesh.V()).value();
3  const scalar meanT = T.weightedAverage(mesh.V()).value();
4  const scalar meanup =
5  (sqrt((2.0/3.0)*turbulence->k()))().weightedAverage(mesh.V()).value();
6 
7  if (Pstream::master())
8  {
9  Info<< "Mean pressure:" << meanp << endl;
10  Info<< "Mean temperature:" << meanT << endl;
11  Info<< "Mean u':" << meanup << endl;
12 
14  << runTime.theta() << tab
15  << meanp << tab
16  << meanT << tab
17  << meanup
18  << endl;
19  }
20 }
static const char tab
Definition: Ostream.H:264
engineTime & runTime
dimensionedScalar sqrt(const dimensionedScalar &ds)
Ostream & endl(Ostream &os)
Add newline and flush stream.
Definition: Ostream.H:256
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::turbulenceModel > turbulence(compressible::turbulenceModel::New(rho, U, phi, thermo))
Definition: createFields.H:94
dynamicFvMesh & mesh
Info<< "Total cylinder mass: "<< fvc::domainIntegrate(rho).value()<< endl;autoPtr< OFstream > logSummaryFile
Definition: startSummary.H:3
const volScalarField & T
messageStream Info
volScalarField & p