createFields.H

Go to the documentation of this file.

if (mesh.nCells() != 1)
{
    FatalErrorInFunction
        << "Solver only applicable to single cell cases"
        << exit(FatalError);
}
 
#include "createInitialFields.H"
 
Info<< "Reading thermophysical properties\n" << endl;
autoPtr<fluidMulticomponentThermo> pThermo
(
    fluidMulticomponentThermo::New(mesh)
);
fluidMulticomponentThermo& thermo = pThermo();
thermo.validate(args.executable(), "h");
 
autoPtr<chemistryModel> pChemistry(chemistryModel::New(thermo));
 
volScalarField rho
(
    IOobject
    (
        "rho",
        runTime.name(),
        mesh,
        IOobject::NO_READ,
        IOobject::AUTO_WRITE
    ),
    thermo.renameRho()
);
 
volScalarField Rspecific
(
    IOobject
    (
        "Rspecific",
        runTime.name(),
        mesh,
        IOobject::NO_READ,
        IOobject::AUTO_WRITE
    ),
    mesh,
    dimensionedScalar
    (
        dimensionSet(dimEnergy/dimMass/dimTemperature),
        0
    )
);
 
volVectorField U
(
    IOobject
    (
        "U",
        runTime.name(),
        mesh,
        IOobject::NO_READ,
        IOobject::NO_WRITE
    ),
    mesh,
    dimensionedVector(dimVelocity, Zero)
);
 
#include "createPhi.H"
 
OFstream post(args.path()/"chemFoam.out");
post<< "# Time" << token::TAB << "Temperature [K]" << token::TAB
    << "Pressure [Pa]" << endl;