makeGraphs.H
Go to the documentation of this file.
2 (
3  IOobject
4  (
5  "R",
6  runTime.timeName(),
7  mesh,
8  IOobject::NO_READ,
9  IOobject::AUTO_WRITE
10  ),
11  turbulence->sigma()
12 );
13 
14 runTime.write();
15 
16 setWriter::New(runTime.graphFormat())->write
17 (
18  runTime.globalPath()
19  /functionObjects::writeFile::outputPrefix
20  /args.executable()
21  /runTime.timeName(),
22 
23  args.executable(),
24 
25  coordSet(true, "y", y),
26 
27  "Uf", (flowDirection & U)().primitiveField(),
28 
29  "nu", turbulence->nu()().primitiveField(),
30  "k", turbulence->nu()().primitiveField(),
31  "epsilon", turbulence->nu()().primitiveField(),
32 
33  "Rff", (flowDirection & R & flowDirection)().primitiveField(),
34  "Rww", (wallNormal & R & wallNormal)().primitiveField(),
35  "Rfw", (flowDirection & R & wallNormal)().primitiveField(),
36 
37  "u", sqrt(mag(R.component(symmTensor::XX)))().primitiveField(),
38  "v", sqrt(mag(R.component(symmTensor::YY)))().primitiveField(),
39  "w", sqrt(mag(R.component(symmTensor::ZZ)))().primitiveField(),
40 
41  "uv", R.component(symmTensor::XY)().primitiveField(),
42 
43  "gammaDot", mag(fvc::grad(U))().primitiveField()
44 );
volSymmTensorField R(IOobject("R", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::AUTO_WRITE), turbulence->sigma())
GeometricField< symmTensor, fvPatchField, volMesh > volSymmTensorField
Definition: volFieldsFwd.H:62
tmp< GeometricField< typename outerProduct< vector, Type >::type, fvPatchField, volMesh >> grad(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
Definition: fvcGrad.C:52
const word & executable() const
Name of executable without the path.
Definition: argListI.H:36
autoPtr< CompressibleMomentumTransportModel > New(const volScalarField &rho, const volVectorField &U, const surfaceScalarField &phi, const viscosity &viscosity)
U
Definition: pEqn.H:72
vector flowDirection
Definition: createFields.H:41
dimensionedScalar sqrt(const dimensionedScalar &ds)
fvMesh & mesh
scalar y
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
dimensioned< scalar > mag(const dimensioned< Type > &)
vector wallNormal(Zero)
Foam::argList args(argc, argv)
runTime write()