1 Pav << runTime.theta() << " " << p.weightedAverage(mesh.V()).value() << endl;
2 Tav << runTime.theta() << " " << T.weightedAverage(mesh.V()).value() << endl;
4 Info << "Max T = " << max(T).value() << " K" << endl;
5 Info << "Min T = " << min(T).value() << " K" << endl;
6 Info << "Max p = " << max(p).value()/1.0e5 << " bar" << endl;
8 massBal << runTime.theta() << " "
9 << fvc::domainIntegrate(rho).value() << endl;
10 volume << runTime.theta() << " " << sum(mesh.V()) << endl;
12 debugT << nl << "Crank angle: " << runTime.theta() << endl;
13 debugT << "Max T = " << max(T).value() << " K" << endl;
14 debugT << "Min T = " << min(T).value() << " K" << endl;
16 kav << runTime.theta() << " "
17 << (turbulence->k())().weightedAverage(mesh.V()).value() << endl;
19 epsav << runTime.theta() << " "
20 << (turbulence->epsilon())().weightedAverage(mesh.V()).value() << endl;
22 debugMinT << runTime.theta() << " " << min(T).value() << endl;
23 debugMaxT << runTime.theta() << " " << max(T).value() << endl;