1 Info<< "Reading thermophysical properties\n" << endl;
3 autoPtr<hsCombustionThermo> pThermo
5 hsCombustionThermo::New(mesh)
8 hsCombustionThermo& thermo = pThermo();
10 basicMultiComponentMixture& composition = thermo.composition();
25 dimensionedScalar stoicRatio
27 thermo.lookup("stoichiometricAirFuelMassRatio")
30 volScalarField& p = thermo.p();
31 volScalarField& hs = thermo.hs();
33 const volScalarField& psi = thermo.psi();
35 volScalarField& ft = composition.Y("ft");
36 volScalarField& fu = composition.Y("fu");
38 Info<< "Reading field U\n" << endl;
53 #include "compressibleCreatePhi.H"
55 Info<< "Creating turbulence model\n" << endl;
56 autoPtr<compressible::turbulenceModel> turbulence
58 compressible::turbulenceModel::New(rho, U, phi, thermo)
61 IOdictionary combustionProperties
65 "combustionProperties",
73 Info<< "Creating combustion model\n" << endl;
74 autoPtr<combustionModel::combustionModel> combustion
76 combustionModel::combustionModel::New
97 dimensionedScalar("dQ", dimMass/pow3(dimTime)/dimLength, 0.0)
100 Info<< "Creating field dpdt\n" << endl;
101 volScalarField dpdt = fvc::ddt(p);
104 Info<< "Calculating field g.h\n" << endl;
105 volScalarField gh("gh", g & mesh.C());
107 surfaceScalarField ghf("gh", g & mesh.Cf());
113 dimensionedScalar initialMass = fvc::domainIntegrate(rho);
116 multivariateSurfaceInterpolationScheme<scalar>::fieldTable fields;
118 if (composition.contains("ft"))
120 fields.add(composition.Y("ft"));
123 if (composition.contains("fu"))
125 fields.add(composition.Y("fu"));