Compare commits
19 commits
| Author | SHA1 | Date | |
|---|---|---|---|
|
ed25887ecc | ||
|
|
2c655680a6 | ||
|
|
cf0d9f904f | ||
|
|
b510f145d2 | ||
|
|
7873f46038 | ||
|
|
e5fd6a4caf | ||
|
|
e2268053c1 | ||
|
|
ccb67c70d2 | ||
|
|
c0264dcdfa | ||
|
|
211cb97540 | ||
|
|
ba0c6c6063 | ||
|
|
df0957a855 | ||
|
|
2e799dbc77 | ||
|
|
786d9fba2e | ||
|
|
cb57fe5a57 | ||
|
|
f8b83b579e | ||
|
|
367133f440 | ||
|
|
7f45daf526 | ||
|
|
3b7f67cd7b |
248 changed files with 332626 additions and 2091515 deletions
2
.gitignore
vendored
2
.gitignore
vendored
|
|
@ -7,6 +7,8 @@
|
|||
*.bak
|
||||
*.bak[0-9][0-9]
|
||||
\#*\#
|
||||
.*.swp
|
||||
.*.swo
|
||||
|
||||
# File-browser settings - anywhere
|
||||
.directory
|
||||
|
|
|
|||
|
|
@ -1,5 +0,0 @@
|
|||
NOxFoam.C
|
||||
|
||||
EXE = $(FOAM_APPBIN)/NOxFoam
|
||||
|
||||
|
||||
|
|
@ -1,10 +0,0 @@
|
|||
EXE_INC = \
|
||||
-I$(LIB_SRC)/finiteVolume/lnInclude \
|
||||
-I${LIB_SRC}/meshTools/lnInclude \
|
||||
-I$(LIB_SRC)/sampling/lnInclude \
|
||||
|
||||
EXE_LIBS = \
|
||||
-lfiniteVolume \
|
||||
-lmeshTools \
|
||||
-lfvOptions \
|
||||
-lsampling
|
||||
|
|
@ -1,17 +0,0 @@
|
|||
// NH3 transport equation
|
||||
while (simple.correctNonOrthogonal())
|
||||
{
|
||||
fvScalarMatrix NH3Eqn
|
||||
(
|
||||
fvm::ddt(rho, NH3)
|
||||
+ fvm::div(phi, NH3)
|
||||
- fvm::laplacian(mut, NH3)
|
||||
==
|
||||
Snh3
|
||||
);
|
||||
|
||||
NH3Eqn.relax();
|
||||
fvOptions.constrain(NH3Eqn);
|
||||
NH3Eqn.solve();
|
||||
fvOptions.correct(NH3);
|
||||
}
|
||||
|
|
@ -1,17 +0,0 @@
|
|||
// NO transport equation
|
||||
while (simple.correctNonOrthogonal())
|
||||
{
|
||||
fvScalarMatrix NOEqn
|
||||
(
|
||||
fvm::ddt(rho, NO)
|
||||
+ fvm::div(phi, NO)
|
||||
- fvm::laplacian(mut, NO)
|
||||
==
|
||||
Sno
|
||||
);
|
||||
|
||||
NOEqn.relax();
|
||||
fvOptions.constrain(NOEqn);
|
||||
NOEqn.solve();
|
||||
fvOptions.correct(NO);
|
||||
}
|
||||
|
|
@ -1,42 +0,0 @@
|
|||
Xno = rho*NO/Wno; // Molar concentration of NO [mol/m^3]
|
||||
|
||||
// Calculation of source term. (Thermal NO)
|
||||
forAll (mesh.cells(),celli)
|
||||
{
|
||||
SthermNO[celli] = Wno.value()*2*kf1[celli]*Xo[celli]*Xn2[celli]*
|
||||
((1-(kr1[celli]*kr2[celli]*Xno[celli]*Xno[celli])/
|
||||
(kf1[celli]*Xn2[celli]*kf2[celli]*Xo2[celli]))/
|
||||
(1+(kr1[celli]*Xno[celli])/
|
||||
(kf2[celli]*Xo2[celli]+kf3[celli]*Xoh[celli])))
|
||||
/runTime.time().deltaTValue();
|
||||
}
|
||||
|
||||
Sno = SthermNO;
|
||||
|
||||
if (calculateFuelNOx)
|
||||
{
|
||||
Xnh3 = rho*NH3/Wnh3; // Molar concentration of nh3[mol/m^3]
|
||||
|
||||
forAll (mesh.cells(),celli)
|
||||
{
|
||||
// Conversion rate of NH3
|
||||
R1[celli] = A1*Xnh3[celli]*Foam::pow(Xo2[celli],a[celli])*Foam::exp(-E1/(RR*T[celli])); // [1/s]
|
||||
R2[celli] = A2*Xnh3[celli]*Xno[celli]*Foam::exp(-E2/(RR*T[celli])); // [1/s]
|
||||
}
|
||||
|
||||
// NH3 sinks
|
||||
Snh3_1 = -R1*Wnh3*p/(Rgas*T); // NH3 consumption -> NO (oxidation)
|
||||
Snh3_2 = -R2*Wnh3*p/(Rgas*T); // NH3 consumption -> N2 (reduction)
|
||||
|
||||
// net NH3 source (NH3 production, consumption 1, and consumption 2)
|
||||
Snh3 = Snh3_p + Snh3_1 + Snh3_2;
|
||||
|
||||
// Fuel NO
|
||||
// Calculation of source term. (Fuel NO)
|
||||
SfuelNO_1 = R1*Wno*p/(Rgas*T); // NH3 + O2 -> NO (source)
|
||||
SfuelNO_2 = -R2*Wno*p/(Rgas*T); // NH3 + NO -> N2 (sink)
|
||||
SfuelNO = SfuelNO_1 + SfuelNO_2; // Sum of fuel NO source terms
|
||||
|
||||
// Sum of NO Sources (source of thermal NO and fuel NO)
|
||||
Sno += SfuelNO;
|
||||
}
|
||||
|
|
@ -1,660 +0,0 @@
|
|||
|
||||
IOdictionary NOxProperties
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"NOxProperties",
|
||||
runTime.constant(),
|
||||
mesh,
|
||||
IOobject::MUST_READ_IF_MODIFIED,
|
||||
IOobject::NO_WRITE
|
||||
)
|
||||
);
|
||||
|
||||
wordList fuelSourceNames(NOxProperties.lookup("fuelSourceNames"));
|
||||
|
||||
PtrList<volScalarField::Internal> fuelSources(fuelSourceNames.size());
|
||||
|
||||
forAll (fuelSources, si)
|
||||
{
|
||||
fuelSources.set
|
||||
(
|
||||
si,
|
||||
new volScalarField::Internal
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
fuelSourceNames[si],
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::MUST_READ,
|
||||
IOobject::AUTO_WRITE
|
||||
),
|
||||
mesh
|
||||
)
|
||||
);
|
||||
}
|
||||
|
||||
scalarField nitrogenMassInFuels(NOxProperties.lookup("nitrogenMassInFuels"));
|
||||
|
||||
bool instantaneousRadicals
|
||||
(
|
||||
NOxProperties.lookupOrDefault("instantaneousRadicals", true)
|
||||
);
|
||||
|
||||
bool calculateFuelNOx
|
||||
(
|
||||
NOxProperties.lookupOrDefault("calculateFuelNOx", false)
|
||||
);
|
||||
|
||||
volScalarField T
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"T",
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::MUST_READ,
|
||||
IOobject::AUTO_WRITE
|
||||
),
|
||||
mesh
|
||||
);
|
||||
|
||||
volScalarField rho
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"rho",
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::MUST_READ,
|
||||
IOobject::AUTO_WRITE
|
||||
),
|
||||
mesh
|
||||
);
|
||||
|
||||
volScalarField p
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"p",
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::MUST_READ,
|
||||
IOobject::AUTO_WRITE
|
||||
),
|
||||
mesh
|
||||
);
|
||||
|
||||
|
||||
Info<< "\nReading field U\n" << endl;
|
||||
volVectorField U
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"U",
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::MUST_READ,
|
||||
IOobject::AUTO_WRITE
|
||||
),
|
||||
mesh
|
||||
);
|
||||
|
||||
Info<< "Reading/calculating face flux field phi\n" << endl;
|
||||
surfaceScalarField phi
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"phi",
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::READ_IF_PRESENT,
|
||||
IOobject::AUTO_WRITE
|
||||
),
|
||||
linearInterpolate(rho*U) & mesh.Sf()
|
||||
);
|
||||
|
||||
volScalarField nut
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"nut",
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::MUST_READ,
|
||||
IOobject::AUTO_WRITE
|
||||
),
|
||||
mesh
|
||||
);
|
||||
|
||||
volScalarField mut
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"mut",
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::NO_READ,
|
||||
IOobject::AUTO_WRITE
|
||||
),
|
||||
rho * nut
|
||||
);
|
||||
|
||||
volScalarField Xno
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"Xno",
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::READ_IF_PRESENT,
|
||||
IOobject::NO_WRITE
|
||||
),
|
||||
mesh,
|
||||
dimensionedScalar("Xno",dimMoles/dimVolume,0.0)
|
||||
);
|
||||
|
||||
volScalarField Xn2
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"Xn2",
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::NO_READ,
|
||||
IOobject::NO_WRITE
|
||||
),
|
||||
mesh,
|
||||
dimensionedScalar("Xn2",dimMoles/dimVolume,0.0)
|
||||
);
|
||||
|
||||
volScalarField Xo2
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"Xo2",
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::NO_READ,
|
||||
IOobject::NO_WRITE
|
||||
),
|
||||
mesh,
|
||||
dimensionedScalar("Xo2",dimMoles/dimVolume,0.0)
|
||||
);
|
||||
|
||||
volScalarField Xh2o
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"Xh2o",
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::NO_READ,
|
||||
IOobject::NO_WRITE
|
||||
),
|
||||
mesh,
|
||||
dimensionedScalar("Xh2o",dimMoles/dimVolume,0.0)
|
||||
);
|
||||
|
||||
volScalarField Xnh3
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"Xnh3",
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::NO_READ,
|
||||
IOobject::NO_WRITE
|
||||
),
|
||||
mesh,
|
||||
dimensionedScalar("Xnh3",dimMoles/dimVolume,0.0)
|
||||
);
|
||||
|
||||
volScalarField Xo
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"Xo",
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::NO_READ,
|
||||
IOobject::NO_WRITE
|
||||
),
|
||||
mesh,
|
||||
dimensionedScalar("Xo",dimMoles/dimVolume,0.0)
|
||||
);
|
||||
|
||||
volScalarField Xoh
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"Xoh",
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::NO_READ,
|
||||
IOobject::NO_WRITE
|
||||
),
|
||||
mesh,
|
||||
dimensionedScalar("Xoh",dimMoles/dimVolume,0.0)
|
||||
);
|
||||
|
||||
volScalarField NO
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"NO",
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::MUST_READ,
|
||||
IOobject::AUTO_WRITE
|
||||
),
|
||||
mesh
|
||||
);
|
||||
|
||||
volScalarField N2
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"N2",
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::MUST_READ,
|
||||
IOobject::AUTO_WRITE
|
||||
),
|
||||
mesh
|
||||
);
|
||||
|
||||
volScalarField O2
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"O2",
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::MUST_READ,
|
||||
IOobject::AUTO_WRITE
|
||||
),
|
||||
mesh
|
||||
);
|
||||
|
||||
volScalarField H2O
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"H2O",
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::MUST_READ,
|
||||
IOobject::AUTO_WRITE
|
||||
),
|
||||
mesh
|
||||
);
|
||||
|
||||
volScalarField O
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"O",
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::MUST_READ,
|
||||
IOobject::AUTO_WRITE
|
||||
),
|
||||
mesh
|
||||
);
|
||||
|
||||
volScalarField OH
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"OH",
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::MUST_READ,
|
||||
IOobject::AUTO_WRITE
|
||||
),
|
||||
mesh
|
||||
);
|
||||
|
||||
volScalarField NH3
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"NH3",
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::READ_IF_PRESENT,
|
||||
IOobject::NO_WRITE
|
||||
),
|
||||
mesh,
|
||||
dimensionedScalar("NH3",dimless,0.0)
|
||||
);
|
||||
|
||||
if (calculateFuelNOx)
|
||||
{
|
||||
NH3.writeOpt() = IOobject::AUTO_WRITE;
|
||||
}
|
||||
|
||||
volScalarField Sno
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"Sno",
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::NO_READ,
|
||||
IOobject::AUTO_WRITE
|
||||
),
|
||||
mesh,
|
||||
dimensionedScalar("Sno",dimDensity/dimTime,0.0)
|
||||
);
|
||||
|
||||
volScalarField SthermNO
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"SthermNO",
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::NO_READ,
|
||||
IOobject::AUTO_WRITE
|
||||
),
|
||||
mesh,
|
||||
dimensionedScalar("SthermNO",dimDensity/dimTime,0.0)
|
||||
);
|
||||
|
||||
volScalarField SfuelNO
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"SfuelNO",
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::NO_READ,
|
||||
IOobject::AUTO_WRITE
|
||||
),
|
||||
mesh,
|
||||
dimensionedScalar("SfuelNO",dimDensity/dimTime,0.0)
|
||||
);
|
||||
|
||||
volScalarField SfuelNO_1
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"SfuelNO_1",
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::NO_READ,
|
||||
IOobject::AUTO_WRITE
|
||||
),
|
||||
mesh,
|
||||
dimensionedScalar("SfuelNO_1",dimDensity/dimTime,0.0)
|
||||
);
|
||||
|
||||
volScalarField SfuelNO_2
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"SfuelNO_2",
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::NO_READ,
|
||||
IOobject::AUTO_WRITE
|
||||
),
|
||||
mesh,
|
||||
dimensionedScalar("SfuelNO_2",dimDensity/dimTime,0.0)
|
||||
);
|
||||
|
||||
volScalarField a
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"a",
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::NO_READ,
|
||||
IOobject::AUTO_WRITE
|
||||
),
|
||||
mesh,
|
||||
dimensionedScalar("a",dimless,0.0)
|
||||
);
|
||||
|
||||
volScalarField R1
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"R1",
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::NO_READ,
|
||||
IOobject::AUTO_WRITE
|
||||
),
|
||||
mesh,
|
||||
dimensionedScalar("R1",dimless/dimTime,0.0)
|
||||
);
|
||||
|
||||
volScalarField R2
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"R2",
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::NO_READ,
|
||||
IOobject::AUTO_WRITE
|
||||
),
|
||||
mesh,
|
||||
dimensionedScalar("R1",dimless/dimTime,0.0)
|
||||
);
|
||||
|
||||
volScalarField Snh3_p
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"Snh3_p",
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::NO_READ,
|
||||
IOobject::AUTO_WRITE
|
||||
),
|
||||
mesh,
|
||||
dimensionedScalar("Snh3_p",dimDensity/dimTime,0.0)
|
||||
);
|
||||
|
||||
volScalarField Snh3_1
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"Snh3_1",
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::NO_READ,
|
||||
IOobject::AUTO_WRITE
|
||||
),
|
||||
mesh,
|
||||
dimensionedScalar("Snh3_1",dimDensity/dimTime,0.0)
|
||||
);
|
||||
|
||||
volScalarField Snh3_2
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"Snh3_2",
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::NO_READ,
|
||||
IOobject::AUTO_WRITE
|
||||
),
|
||||
mesh,
|
||||
dimensionedScalar("Snh3_2",dimDensity/dimTime,0.0)
|
||||
);
|
||||
|
||||
volScalarField Snh3
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"Snh3",
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::NO_READ,
|
||||
IOobject::AUTO_WRITE
|
||||
),
|
||||
mesh,
|
||||
dimensionedScalar("Snh3",dimDensity/dimTime,0.0)
|
||||
);
|
||||
|
||||
volScalarField kf1
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"kf1",
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::NO_READ,
|
||||
IOobject::AUTO_WRITE
|
||||
),
|
||||
mesh,
|
||||
dimensionedScalar("kf1",dimVolume/dimTime/dimMoles,0.0)
|
||||
);
|
||||
|
||||
volScalarField kf2
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"kf2",
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::NO_READ,
|
||||
IOobject::AUTO_WRITE
|
||||
),
|
||||
mesh,
|
||||
dimensionedScalar("kf2",dimVolume/dimTime/dimMoles,0.0)
|
||||
);
|
||||
|
||||
volScalarField kf3
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"kf3",
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::NO_READ,
|
||||
IOobject::AUTO_WRITE
|
||||
),
|
||||
mesh,
|
||||
dimensionedScalar("kf3",dimVolume/dimTime/dimMoles,0.0)
|
||||
);
|
||||
|
||||
volScalarField kr1
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"kr1",
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::NO_READ,
|
||||
IOobject::AUTO_WRITE
|
||||
),
|
||||
mesh,
|
||||
dimensionedScalar("kr1",dimVolume/dimTime/dimMoles,0.0)
|
||||
);
|
||||
|
||||
volScalarField kr2
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"kr2",
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::NO_READ,
|
||||
IOobject::AUTO_WRITE
|
||||
),
|
||||
mesh,
|
||||
dimensionedScalar("kr2",dimVolume/dimTime/dimMoles,0.0)
|
||||
);
|
||||
|
||||
const dimensionedScalar Wno ("Wno", dimMass/dimMoles, 0.030006); // Molecular weight of NO [kg/mol]
|
||||
const dimensionedScalar Wn2 ("Wn2", dimMass/dimMoles, 0.028013); // Molecular weight of N2 [kg/mol]
|
||||
const dimensionedScalar Wo2 ("Wo2", dimMass/dimMoles, 0.031998); // Molecular weight of O2 [kg/mol]
|
||||
const dimensionedScalar Wh2o ("Wh2o", dimMass/dimMoles, 0.018015); // Molecular weight of H2O[kg/mol]
|
||||
const dimensionedScalar Wo ("Wo", dimMass/dimMoles, 0.015999); // Molecular weight of O [kg/mol]
|
||||
const dimensionedScalar Woh ("Woh", dimMass/dimMoles, 0.017008); // Molecular weight of OH [kg/mol]
|
||||
const dimensionedScalar Wnh3 ("Wnh3", dimMass/dimMoles, 0.017031); // Molecular weight of NH3[kg/mol]
|
||||
const dimensionedScalar Wn ("Wn", dimMass/dimMoles, 0.014007); // Molecular weight of N [kg/mol]
|
||||
|
||||
const dimensionedScalar Rgas (Foam::constant::physicoChemical::R);
|
||||
const scalar RR = Rgas.value(); // Universal gas constant [J/(mol*K)]
|
||||
const scalar A1 = 4.0E+06; // [1/s]
|
||||
const scalar A2 = 1.8E+08; // [1/s]
|
||||
const scalar E1 = 133947.2; // [J/mol]
|
||||
const scalar E2 = 113017.95; // [J/mol]
|
||||
|
||||
// molar concentrations
|
||||
Xn2 = rho*N2/Wn2; // Molar concentration of N2 [mol/m^3]
|
||||
Xo2 = rho*O2/Wo2; // Molar concentration of O2 [mol/m^3]
|
||||
Xh2o = rho*H2O/Wh2o; // Molar concentration of H2O[mol/m^3]
|
||||
|
||||
// O,OH molar concentrations
|
||||
if(instantaneousRadicals)
|
||||
{
|
||||
Xo = rho*O/Wo; // Molar concentration of O [mol/m^3]
|
||||
Xoh = rho*OH/Woh; // Molar concentration of OH[mol/m^3]
|
||||
}
|
||||
else
|
||||
{
|
||||
forAll (mesh.cells(),celli)
|
||||
{
|
||||
// Molar concentration of O [mol/m^3], Partial Equilibrium Approach
|
||||
Xo[celli] = 36.64 * Foam::pow(T[celli],0.5)
|
||||
* Foam::pow(Xo2[celli],0.5) * Foam::exp(-27123.0/T[celli]);
|
||||
// Molar concentration of OH [mol/m^3], Partial Equilibrium Approach
|
||||
Xoh[celli] = 2.129E+02 * Foam::pow(T[celli],-0.57) * Foam::exp(-4595.0/T[celli])
|
||||
* Foam::pow(Xo[celli],0.5) * Foam::pow(Xh2o[celli],0.5);
|
||||
}
|
||||
}
|
||||
|
||||
forAll (mesh.cells(),celli)
|
||||
{
|
||||
// Thermal NO
|
||||
// Reaction rate constants [m^3/(mol*s)]
|
||||
kf1[celli] = 1.8E+08*Foam::exp(-38370/T[celli]);
|
||||
kf2[celli] = 1.8E+04*T[celli]*Foam::exp(-4680/T[celli]);
|
||||
kf3[celli] = 7.1E+07*Foam::exp(-450/T[celli]);
|
||||
kr1[celli] = 3.8E+07*Foam::exp(-425/T[celli]);
|
||||
kr2[celli] = 3.81E+03*T[celli]*Foam::exp(-20820/T[celli]);
|
||||
}
|
||||
|
||||
// a, Order of Reaction "NH3 + O2 -> NO"
|
||||
forAll (mesh.cells(),celli)
|
||||
{
|
||||
if (Xo2[celli] <= 4.1E-03)
|
||||
{
|
||||
a[celli] = 1.0;
|
||||
}
|
||||
else if ( Xo2[celli] > 4.1E-03 and Xo2[celli] <= 1.11E-02 )
|
||||
{
|
||||
a[celli] = -3.95 - 0.9*Foam::log(Xo2[celli]);
|
||||
}
|
||||
else if ( Xo2[celli] > 1.11E-02 and Xo2[celli] <= 0.03 )
|
||||
{
|
||||
a[celli] = -0.35 - 0.1*Foam::log(Xo2[celli]);
|
||||
}
|
||||
else // if ( Xo2[celli] > 0.03 )
|
||||
{
|
||||
a[celli] = 0.0;
|
||||
}
|
||||
}
|
||||
|
||||
// NH3 source
|
||||
forAll (mesh.cells(),celli)
|
||||
{
|
||||
scalar nh3Source = 0.0;
|
||||
forAll (fuelSources, si)
|
||||
{
|
||||
nh3Source += fuelSources[si][celli] * nitrogenMassInFuels[si];
|
||||
}
|
||||
|
||||
Snh3_p[celli] = nh3Source*Wnh3.value()/Wn.value()/mesh.V()[celli]; // NH3 production
|
||||
}
|
||||
|
||||
49
applications/solvers/combustion/SLFMFoam/AMC/AMC.C
Normal file
49
applications/solvers/combustion/SLFMFoam/AMC/AMC.C
Normal file
|
|
@ -0,0 +1,49 @@
|
|||
/*---------------------------------------------------------------------------*\
|
||||
========= |
|
||||
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
||||
\\ / O peration |
|
||||
\\ / A nd | Copyright (C) 2016 OpenFOAM Foundation
|
||||
\\/ M anipulation |
|
||||
-------------------------------------------------------------------------------
|
||||
License
|
||||
This file is part of OpenFOAM.
|
||||
|
||||
OpenFOAM is free software: you can redistribute it and/or modify it
|
||||
under the terms of the GNU General Public License as published by
|
||||
the Free Software Foundation, either version 3 of the License, or
|
||||
(at your option) any later version.
|
||||
|
||||
OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
|
||||
ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
|
||||
FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
|
||||
for more details.
|
||||
|
||||
You should have received a copy of the GNU General Public License
|
||||
along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
|
||||
|
||||
\*---------------------------------------------------------------------------*/
|
||||
|
||||
#include "AMC.H"
|
||||
|
||||
#include "error.H"
|
||||
|
||||
// * * * * * * * * * * * * * Functions * * * * * * * * * * * * //
|
||||
|
||||
//Amplitude Mapping Closure (from KIVA)
|
||||
Foam::tmp<Foam::scalarField> Foam::CMC::AMC(const scalarField& eta)
|
||||
{
|
||||
tmp<scalarField> tRes(new scalarField(eta.size(), 0.0));
|
||||
|
||||
scalarField& Res = tRes.ref();
|
||||
|
||||
for(label i=0 ; i<eta.size() ; i++)
|
||||
{
|
||||
Res[i] = AMC(eta[i]);
|
||||
}
|
||||
|
||||
return tRes;
|
||||
}
|
||||
|
||||
|
||||
|
||||
// ************************************************************************* //
|
||||
74
applications/solvers/combustion/SLFMFoam/AMC/AMC.H
Normal file
74
applications/solvers/combustion/SLFMFoam/AMC/AMC.H
Normal file
|
|
@ -0,0 +1,74 @@
|
|||
/*---------------------------------------------------------------------------*\
|
||||
========= |
|
||||
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
||||
\\ / O peration |
|
||||
\\ / A nd | Copyright (C) 2016 OpenFOAM Foundation
|
||||
\\/ M anipulation |
|
||||
-------------------------------------------------------------------------------
|
||||
License
|
||||
This file is part of OpenFOAM.
|
||||
|
||||
OpenFOAM is free software: you can redistribute it and/or modify it
|
||||
under the terms of the GNU General Public License as published by
|
||||
the Free Software Foundation, either version 3 of the License, or
|
||||
(at your option) any later version.
|
||||
|
||||
OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
|
||||
ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
|
||||
FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
|
||||
for more details.
|
||||
|
||||
You should have received a copy of the GNU General Public License
|
||||
along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
|
||||
|
||||
Function
|
||||
Foam::CMC::AMC
|
||||
|
||||
Description
|
||||
|
||||
SourceFiles
|
||||
AMC.C
|
||||
|
||||
\*---------------------------------------------------------------------------*/
|
||||
|
||||
#ifndef AMC_H
|
||||
#define AMC_H
|
||||
|
||||
#include "tmp.H"
|
||||
#include "scalarField.H"
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
namespace Foam
|
||||
{
|
||||
namespace CMC
|
||||
{
|
||||
|
||||
/*---------------------------------------------------------------------------*\
|
||||
Function AMC Declaration
|
||||
\*---------------------------------------------------------------------------*/
|
||||
|
||||
const scalar pi = 3.141592;
|
||||
|
||||
const scalar spi = Foam::sqrt(pi);
|
||||
|
||||
//Amplitude Mapping Closure (from KIVA)
|
||||
tmp<scalarField> AMC(const scalarField& eta);
|
||||
|
||||
|
||||
//Amplitude Mapping Closure (from KIVA)
|
||||
//Define exp(-2*(erf^-1(2*eta - 1))^2)
|
||||
inline scalar AMC(const scalar eta);
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
} // End namespace CMC
|
||||
} // End namespace Foam
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
#include "AMCI.H"
|
||||
|
||||
#endif
|
||||
|
||||
// ************************************************************************* //
|
||||
54
applications/solvers/combustion/SLFMFoam/AMC/AMCI.H
Normal file
54
applications/solvers/combustion/SLFMFoam/AMC/AMCI.H
Normal file
|
|
@ -0,0 +1,54 @@
|
|||
/*---------------------------------------------------------------------------*\
|
||||
========= |
|
||||
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
||||
\\ / O peration |
|
||||
\\ / A nd | Copyright (C) 2016 OpenFOAM Foundation
|
||||
\\/ M anipulation |
|
||||
-------------------------------------------------------------------------------
|
||||
License
|
||||
This file is part of OpenFOAM.
|
||||
|
||||
OpenFOAM is free software: you can redistribute it and/or modify it
|
||||
under the terms of the GNU General Public License as published by
|
||||
the Free Software Foundation, either version 3 of the License, or
|
||||
(at your option) any later version.
|
||||
|
||||
OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
|
||||
ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
|
||||
FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
|
||||
for more details.
|
||||
|
||||
You should have received a copy of the GNU General Public License
|
||||
along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
|
||||
|
||||
\*---------------------------------------------------------------------------*/
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
// * * * * * * * * * * * * * Private Member Functions * * * * * * * * * * * //
|
||||
|
||||
|
||||
//Amplitude Mapping Closure (from KIVA)
|
||||
//Define exp(-2*(erf^-1(2*eta - 1))^2)
|
||||
inline Foam::scalar Foam::CMC::AMC(const scalar eta)
|
||||
{
|
||||
const scalar a0 = (2.0*eta - 1.0);
|
||||
scalar a1 = 0.5;
|
||||
scalar slope;
|
||||
scalar da = GREAT;
|
||||
|
||||
while(mag(da) > SMALL)
|
||||
{
|
||||
slope = (2.0/spi)*Foam::exp(-1.0*Foam::pow(a1,2.0));
|
||||
da = (a0 - Foam::erf(a1))/slope;
|
||||
a1 = a1 + da;
|
||||
}
|
||||
|
||||
return Foam::exp(-2.0*Foam::pow(a1,2.0));
|
||||
}
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
|
||||
// ************************************************************************* //
|
||||
|
|
@ -0,0 +1,120 @@
|
|||
/*---------------------------------------------------------------------------*\
|
||||
========= |
|
||||
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
||||
\\ / O peration |
|
||||
\\ / A nd | Copyright (C) 2016 OpenFOAM Foundation
|
||||
\\/ M anipulation |
|
||||
-------------------------------------------------------------------------------
|
||||
License
|
||||
This file is part of OpenFOAM.
|
||||
|
||||
OpenFOAM is free software: you can redistribute it and/or modify it
|
||||
under the terms of the GNU General Public License as published by
|
||||
the Free Software Foundation, either version 3 of the License, or
|
||||
(at your option) any later version.
|
||||
|
||||
OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
|
||||
ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
|
||||
FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
|
||||
for more details.
|
||||
|
||||
You should have received a copy of the GNU General Public License
|
||||
along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
|
||||
|
||||
\*---------------------------------------------------------------------------*/
|
||||
|
||||
#include "BetaFunction.H"
|
||||
|
||||
#include "error.H"
|
||||
|
||||
// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
|
||||
|
||||
// const scalar Foam::CMC::BetaFunction::staticData();
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * Static Member Functions * * * * * * * * * * * * //
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * Private Member Functions * * * * * * * * * * * //
|
||||
|
||||
void Foam::CMC::BetaFunction::_setParameters()
|
||||
{
|
||||
scalar gamma = mf_*(1.0 - mf_)/(mfVar_ + SMALL) - 1.0;
|
||||
|
||||
if (gamma < SMALL || (mfVar_)/(mf_*(1 - mf_) + SMALL) < 0.001)
|
||||
{
|
||||
fdelta_ = true;
|
||||
}
|
||||
else
|
||||
{
|
||||
a_ = max(mf_*gamma, 0.0);
|
||||
b_ = max((1.0 - mf_)*gamma, 0.0);
|
||||
|
||||
if(a_ < 1.0)
|
||||
{
|
||||
delta0_ = true;
|
||||
}
|
||||
if(b_ < 1.0)
|
||||
{
|
||||
delta1_ = true;
|
||||
}
|
||||
|
||||
_limitAB();
|
||||
}
|
||||
}
|
||||
|
||||
|
||||
void Foam::CMC::BetaFunction::_limitAB()
|
||||
{
|
||||
scalar fmax = 1.0/(1.0 + (b_ - 1.0)/(a_ - 1.0));
|
||||
|
||||
if(a_ > 500.0)
|
||||
{
|
||||
a_ = 500.0;
|
||||
b_ = (a_ - 1.0 - fmax*(a_ - 2.0))/fmax;
|
||||
}
|
||||
else if(b_ > 500.0)
|
||||
{
|
||||
b_ = 500.0;
|
||||
a_ = (1.0 + fmax*(b_ - 2.0))/(1.0 - fmax);
|
||||
}
|
||||
}
|
||||
|
||||
// * * * * * * * * * * * * Protected Member Functions * * * * * * * * * * * //
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
|
||||
|
||||
Foam::CMC::BetaFunction::BetaFunction(const scalar mf, const scalar mfVar)
|
||||
:
|
||||
mf_(mf),
|
||||
mfVar_(mfVar),
|
||||
a_(0.0),
|
||||
b_(0.0),
|
||||
fdelta_(false),
|
||||
delta0_(false),
|
||||
delta1_(false)
|
||||
{
|
||||
_setParameters();
|
||||
}
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * * Destructor * * * * * * * * * * * * * * * //
|
||||
|
||||
Foam::CMC::BetaFunction::~BetaFunction()
|
||||
{}
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * * //
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * Member Operators * * * * * * * * * * * * * * //
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * Friend Functions * * * * * * * * * * * * * * //
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * Friend Operators * * * * * * * * * * * * * * //
|
||||
|
||||
|
||||
// ************************************************************************* //
|
||||
|
|
@ -0,0 +1,199 @@
|
|||
/*---------------------------------------------------------------------------*\
|
||||
========= |
|
||||
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
||||
\\ / O peration |
|
||||
\\ / A nd | Copyright (C) 2016 OpenFOAM Foundation
|
||||
\\/ M anipulation |
|
||||
-------------------------------------------------------------------------------
|
||||
License
|
||||
This file is part of OpenFOAM.
|
||||
|
||||
OpenFOAM is free software: you can redistribute it and/or modify it
|
||||
under the terms of the GNU General Public License as published by
|
||||
the Free Software Foundation, either version 3 of the License, or
|
||||
(at your option) any later version.
|
||||
|
||||
OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
|
||||
ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
|
||||
FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
|
||||
for more details.
|
||||
|
||||
You should have received a copy of the GNU General Public License
|
||||
along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
|
||||
|
||||
Class
|
||||
Foam::BetaFunction
|
||||
|
||||
Description
|
||||
|
||||
SourceFiles
|
||||
BetaFunctionI.H
|
||||
BetaFunction.C
|
||||
BetaFunctionIO.C
|
||||
|
||||
\*---------------------------------------------------------------------------*/
|
||||
|
||||
#ifndef BetaFunction_H
|
||||
#define BetaFunction_H
|
||||
|
||||
#include "scalar.H"
|
||||
#include "Switch.H"
|
||||
#include "scalarField.H"
|
||||
#include "tmp.H"
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
namespace Foam
|
||||
{
|
||||
|
||||
// Forward declaration of classes
|
||||
class Istream;
|
||||
class Ostream;
|
||||
|
||||
// Forward declaration of friend functions and operators
|
||||
namespace CMC
|
||||
{
|
||||
class BetaFunction;
|
||||
}
|
||||
Istream& operator>>(Istream&, CMC::BetaFunction&);
|
||||
Ostream& operator<<(Ostream&, const CMC::BetaFunction&);
|
||||
|
||||
|
||||
/*---------------------------------------------------------------------------*\
|
||||
Class BetaFunction Declaration
|
||||
\*---------------------------------------------------------------------------*/
|
||||
|
||||
namespace CMC
|
||||
{
|
||||
|
||||
class BetaFunction
|
||||
{
|
||||
// Private data
|
||||
|
||||
//- Mean mixture fraction
|
||||
scalar mf_;
|
||||
|
||||
//- Mixture fraction variance
|
||||
scalar mfVar_;
|
||||
|
||||
//- Beta pdf parameter alpha
|
||||
scalar a_;
|
||||
|
||||
//- Beta pdf parameter beta
|
||||
scalar b_;
|
||||
|
||||
//- Flag beta pdf special case: forced-delta function
|
||||
Switch fdelta_;
|
||||
|
||||
//- delta function at pure oxidizer stream
|
||||
Switch delta0_;
|
||||
|
||||
//- delta function at pure fuel stream
|
||||
Switch delta1_;
|
||||
|
||||
|
||||
// Private Member Functions
|
||||
|
||||
//- Calculate alpha and beta and set special case flags
|
||||
void _setParameters();
|
||||
|
||||
//- Limit too large value of beta pdf parameters
|
||||
void _limitAB();
|
||||
|
||||
//- Disallow default bitwise copy construct
|
||||
// BetaFunction(const BetaFunction&);
|
||||
|
||||
//- Disallow default bitwise assignment
|
||||
// void operator=(const BetaFunction&);
|
||||
|
||||
|
||||
public:
|
||||
|
||||
// Static data members
|
||||
|
||||
//- Static data staticData
|
||||
// static const scalar staticData;
|
||||
|
||||
|
||||
// Constructors
|
||||
|
||||
//- Construct from components
|
||||
BetaFunction(const scalar mf, const scalar mfVar);
|
||||
|
||||
//- Construct from Istream
|
||||
BetaFunction(Istream&);
|
||||
|
||||
//- Construct as copy
|
||||
// BetaFunction(const BetaFunction&);
|
||||
|
||||
|
||||
//- Destructor
|
||||
~BetaFunction();
|
||||
|
||||
|
||||
// Member Functions
|
||||
|
||||
tmp<scalarField> etaFunc(const UList<scalar>& eta) const
|
||||
{
|
||||
return pow(eta, a_-1.0)*pow(1.0-eta, b_-1.0);
|
||||
}
|
||||
|
||||
scalar etaFunc(const scalar eta) const
|
||||
{
|
||||
return Foam::pow(eta, a_-1.0)*Foam::pow(1.0-eta, b_-1.0);
|
||||
}
|
||||
|
||||
|
||||
// Access
|
||||
|
||||
Switch fdelta() const {return fdelta_;}
|
||||
|
||||
Switch delta0() const {return delta0_;}
|
||||
|
||||
Switch delta1() const {return delta1_;}
|
||||
|
||||
scalar mf() const {return mf_;}
|
||||
|
||||
scalar mfVar() const {return mfVar_;}
|
||||
|
||||
scalar alpha() const {return a_;}
|
||||
|
||||
scalar beta() const {return b_;}
|
||||
|
||||
// Check
|
||||
|
||||
// Edit
|
||||
|
||||
// Write
|
||||
|
||||
|
||||
// Member Operators
|
||||
|
||||
// void operator=(const BetaFunction&);
|
||||
|
||||
|
||||
// Friend Functions
|
||||
|
||||
// Friend Operators
|
||||
|
||||
// IOstream Operators
|
||||
|
||||
friend Istream& ::Foam::operator>>(Istream&, BetaFunction&);
|
||||
friend Ostream& ::Foam::operator<<(Ostream&, const BetaFunction&);
|
||||
};
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
} // End namespace CMC
|
||||
} // End namespace Foam
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
#include "BetaFunctionI.H"
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
#endif
|
||||
|
||||
// ************************************************************************* //
|
||||
|
|
@ -0,0 +1,58 @@
|
|||
/*---------------------------------------------------------------------------*\
|
||||
========= |
|
||||
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
||||
\\ / O peration |
|
||||
\\ / A nd | Copyright (C) 2016 OpenFOAM Foundation
|
||||
\\/ M anipulation |
|
||||
-------------------------------------------------------------------------------
|
||||
License
|
||||
This file is part of OpenFOAM.
|
||||
|
||||
OpenFOAM is free software: you can redistribute it and/or modify it
|
||||
under the terms of the GNU General Public License as published by
|
||||
the Free Software Foundation, either version 3 of the License, or
|
||||
(at your option) any later version.
|
||||
|
||||
OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
|
||||
ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
|
||||
FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
|
||||
for more details.
|
||||
|
||||
You should have received a copy of the GNU General Public License
|
||||
along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
|
||||
|
||||
\*---------------------------------------------------------------------------*/
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
// * * * * * * * * * * * * * Private Member Functions * * * * * * * * * * * //
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * * Destructor * * * * * * * * * * * * * * * //
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * Member Operators * * * * * * * * * * * * * //
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * Friend Functions * * * * * * * * * * * * * //
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * Friend Operators * * * * * * * * * * * * * //
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * IOstream Operators * * * * * * * * * * * * * //
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * Ostream Operator * * * * * * * * * * * * * //
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
|
||||
// ************************************************************************* //
|
||||
|
|
@ -0,0 +1,65 @@
|
|||
/*---------------------------------------------------------------------------*\
|
||||
========= |
|
||||
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
||||
\\ / O peration |
|
||||
\\ / A nd | Copyright (C) 2016 OpenFOAM Foundation
|
||||
\\/ M anipulation |
|
||||
-------------------------------------------------------------------------------
|
||||
License
|
||||
This file is part of OpenFOAM.
|
||||
|
||||
OpenFOAM is free software: you can redistribute it and/or modify it
|
||||
under the terms of the GNU General Public License as published by
|
||||
the Free Software Foundation, either version 3 of the License, or
|
||||
(at your option) any later version.
|
||||
|
||||
OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
|
||||
ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
|
||||
FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
|
||||
for more details.
|
||||
|
||||
You should have received a copy of the GNU General Public License
|
||||
along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
|
||||
|
||||
\*---------------------------------------------------------------------------*/
|
||||
|
||||
#include "BetaFunction.H"
|
||||
#include "IOstreams.H"
|
||||
|
||||
// * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
|
||||
|
||||
Foam::CMC::BetaFunction::BetaFunction(Istream& is)
|
||||
{
|
||||
// Check state of Istream
|
||||
is.check("Foam::CMC::BetaFunction::BetaFunction(Foam::Istream&)");
|
||||
}
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * IOstream Operators * * * * * * * * * * * * //
|
||||
|
||||
Foam::Istream& Foam::operator>>(Istream& is, CMC::BetaFunction&)
|
||||
{
|
||||
// Check state of Istream
|
||||
is.check
|
||||
(
|
||||
"Foam::Istream& Foam::operator>>(Foam::Istream&, Foam::CMC::BetaFunction&)"
|
||||
);
|
||||
|
||||
return is;
|
||||
}
|
||||
|
||||
|
||||
Foam::Ostream& Foam::operator<<(Ostream& os, const CMC::BetaFunction&)
|
||||
{
|
||||
// Check state of Ostream
|
||||
os.check
|
||||
(
|
||||
"Foam::Ostream& Foam::operator<<(Foam::Ostream&, "
|
||||
"const Foam::CMC::BetaFunction&)"
|
||||
);
|
||||
|
||||
return os;
|
||||
}
|
||||
|
||||
|
||||
// ************************************************************************* //
|
||||
126
applications/solvers/combustion/SLFMFoam/BetaGrid/BetaGrid.C
Normal file
126
applications/solvers/combustion/SLFMFoam/BetaGrid/BetaGrid.C
Normal file
|
|
@ -0,0 +1,126 @@
|
|||
/*---------------------------------------------------------------------------*\
|
||||
========= |
|
||||
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
||||
\\ / O peration |
|
||||
\\ / A nd | Copyright (C) 2016 OpenFOAM Foundation
|
||||
\\/ M anipulation |
|
||||
-------------------------------------------------------------------------------
|
||||
License
|
||||
This file is part of OpenFOAM.
|
||||
|
||||
OpenFOAM is free software: you can redistribute it and/or modify it
|
||||
under the terms of the GNU General Public License as published by
|
||||
the Free Software Foundation, either version 3 of the License, or
|
||||
(at your option) any later version.
|
||||
|
||||
OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
|
||||
ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
|
||||
FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
|
||||
for more details.
|
||||
|
||||
You should have received a copy of the GNU General Public License
|
||||
along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
|
||||
|
||||
\*---------------------------------------------------------------------------*/
|
||||
|
||||
#include "BetaGrid.H"
|
||||
|
||||
#include "error.H"
|
||||
#include "dictionary.H"
|
||||
#include "SubList.H"
|
||||
|
||||
// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
|
||||
|
||||
// const scalar Foam::CMC::BetaGrid::staticData();
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * Static Member Functions * * * * * * * * * * * * //
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * Private Member Functions * * * * * * * * * * * //
|
||||
|
||||
|
||||
// * * * * * * * * * * * * Protected Member Functions * * * * * * * * * * * //
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
|
||||
|
||||
Foam::CMC::BetaGrid::BetaGrid(const dictionary& dict)
|
||||
:
|
||||
etaCut_(dict.lookup("detailedEta")),
|
||||
nSpacings_(dict.lookup("detailedN")),
|
||||
etaSpace_(sum(nSpacings_)+3)
|
||||
{
|
||||
// Check for input list size
|
||||
if (etaCut_.size() != nSpacings_.size() + 1)
|
||||
{
|
||||
FatalErrorInFunction
|
||||
<< "Number of grid interval points does not match number of grid sizes"
|
||||
<< abort(FatalError);
|
||||
}
|
||||
|
||||
// Check for evenness of number of spacings in each interval
|
||||
forAll (nSpacings_, i)
|
||||
{
|
||||
if ((nSpacings_[i] % 2) == 1)
|
||||
{
|
||||
FatalErrorInFunction
|
||||
<< "Number of spacings in intervals should be even"
|
||||
<< abort(FatalError);
|
||||
}
|
||||
}
|
||||
|
||||
forAll(nSpacings_, i)
|
||||
{
|
||||
const label ns = nSpacings_[i];
|
||||
const labelList::subList prev(nSpacings_, i);
|
||||
const label baseI = sum(prev) + 1;
|
||||
|
||||
const scalar low = etaCut_[i];
|
||||
const scalar upp = etaCut_[i+1];
|
||||
const scalar delta = (upp - low) / scalar(ns);
|
||||
|
||||
for (label j = 0; j < ns; j++)
|
||||
{
|
||||
etaSpace_[baseI+j] = low + delta * j;
|
||||
}
|
||||
}
|
||||
etaSpace_[0] = 0.0;
|
||||
etaSpace_[etaSpace_.size()-2] = etaCut_.last();
|
||||
etaSpace_[etaSpace_.size()-1] = 1.0;
|
||||
}
|
||||
|
||||
|
||||
Foam::CMC::BetaGrid::BetaGrid(const BetaGrid&)
|
||||
{}
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * * Destructor * * * * * * * * * * * * * * * //
|
||||
|
||||
Foam::CMC::BetaGrid::~BetaGrid()
|
||||
{}
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * * //
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * Member Operators * * * * * * * * * * * * * * //
|
||||
|
||||
void Foam::CMC::BetaGrid::operator=(const BetaGrid& rhs)
|
||||
{
|
||||
// Check for assignment to self
|
||||
if (this == &rhs)
|
||||
{
|
||||
FatalErrorInFunction
|
||||
<< "Attempted assignment to self"
|
||||
<< abort(FatalError);
|
||||
}
|
||||
}
|
||||
|
||||
// * * * * * * * * * * * * * * Friend Functions * * * * * * * * * * * * * * //
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * Friend Operators * * * * * * * * * * * * * * //
|
||||
|
||||
|
||||
// ************************************************************************* //
|
||||
161
applications/solvers/combustion/SLFMFoam/BetaGrid/BetaGrid.H
Normal file
161
applications/solvers/combustion/SLFMFoam/BetaGrid/BetaGrid.H
Normal file
|
|
@ -0,0 +1,161 @@
|
|||
/*---------------------------------------------------------------------------*\
|
||||
========= |
|
||||
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
||||
\\ / O peration |
|
||||
\\ / A nd | Copyright (C) 2016 OpenFOAM Foundation
|
||||
\\/ M anipulation |
|
||||
-------------------------------------------------------------------------------
|
||||
License
|
||||
This file is part of OpenFOAM.
|
||||
|
||||
OpenFOAM is free software: you can redistribute it and/or modify it
|
||||
under the terms of the GNU General Public License as published by
|
||||
the Free Software Foundation, either version 3 of the License, or
|
||||
(at your option) any later version.
|
||||
|
||||
OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
|
||||
ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
|
||||
FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
|
||||
for more details.
|
||||
|
||||
You should have received a copy of the GNU General Public License
|
||||
along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
|
||||
|
||||
Class
|
||||
Foam::BetaGrid
|
||||
|
||||
Description
|
||||
|
||||
SourceFiles
|
||||
BetaGridI.H
|
||||
BetaGrid.C
|
||||
BetaGridIO.C
|
||||
|
||||
\*---------------------------------------------------------------------------*/
|
||||
|
||||
#ifndef BetaGrid_H
|
||||
#define BetaGrid_H
|
||||
|
||||
#include "scalarField.H"
|
||||
#include "labelList.H"
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
namespace Foam
|
||||
{
|
||||
|
||||
// Forward declaration of classes
|
||||
class Istream;
|
||||
class Ostream;
|
||||
|
||||
// Forward declaration of friend functions and operators
|
||||
namespace CMC
|
||||
{
|
||||
class BetaGrid;
|
||||
}
|
||||
class dictionary;
|
||||
|
||||
Istream& operator>>(Istream&, CMC::BetaGrid&);
|
||||
Ostream& operator<<(Ostream&, const CMC::BetaGrid&);
|
||||
|
||||
|
||||
/*---------------------------------------------------------------------------*\
|
||||
Class BetaGrid Declaration
|
||||
\*---------------------------------------------------------------------------*/
|
||||
|
||||
namespace CMC
|
||||
{
|
||||
|
||||
class BetaGrid
|
||||
{
|
||||
// Private data
|
||||
|
||||
//- Piecewise uniform interval points
|
||||
scalarField etaCut_;
|
||||
|
||||
//- Sizes of piecewise uniform intervals
|
||||
labelList nSpacings_;
|
||||
|
||||
//- Mixture fraction eta grid points
|
||||
scalarField etaSpace_;
|
||||
|
||||
|
||||
// Private Member Functions
|
||||
|
||||
//- Disallow default bitwise copy construct
|
||||
BetaGrid(const BetaGrid&);
|
||||
|
||||
//- Disallow default bitwise assignment
|
||||
void operator=(const BetaGrid&);
|
||||
|
||||
|
||||
public:
|
||||
|
||||
// Static data members
|
||||
|
||||
//- Static data staticData
|
||||
// static const scalar staticData;
|
||||
|
||||
|
||||
// Constructors
|
||||
|
||||
//- Construct from components
|
||||
BetaGrid(const dictionary& dict);
|
||||
|
||||
//- Construct from Istream
|
||||
BetaGrid(Istream&);
|
||||
|
||||
//- Construct as copy
|
||||
// BetaGrid(const BetaGrid&);
|
||||
|
||||
|
||||
//- Destructor
|
||||
~BetaGrid();
|
||||
|
||||
|
||||
// Member Functions
|
||||
|
||||
// Access
|
||||
const labelList& N() const {return nSpacings_;}
|
||||
|
||||
const scalarField& etaCut() const {return etaCut_;}
|
||||
|
||||
const scalarField& etaSpace() const {return etaSpace_;}
|
||||
|
||||
// Check
|
||||
|
||||
// Edit
|
||||
|
||||
// Write
|
||||
|
||||
|
||||
// Member Operators
|
||||
|
||||
// void operator=(const BetaGrid&);
|
||||
|
||||
|
||||
// Friend Functions
|
||||
|
||||
// Friend Operators
|
||||
|
||||
// IOstream Operators
|
||||
|
||||
friend Istream& ::Foam::operator>>(Istream&, BetaGrid&);
|
||||
friend Ostream& ::Foam::operator<<(Ostream&, const BetaGrid&);
|
||||
};
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
} // End namespace CMC
|
||||
} // End namespace Foam
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
#include "BetaGridI.H"
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
#endif
|
||||
|
||||
// ************************************************************************* //
|
||||
|
|
@ -0,0 +1,58 @@
|
|||
/*---------------------------------------------------------------------------*\
|
||||
========= |
|
||||
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
||||
\\ / O peration |
|
||||
\\ / A nd | Copyright (C) 2016 OpenFOAM Foundation
|
||||
\\/ M anipulation |
|
||||
-------------------------------------------------------------------------------
|
||||
License
|
||||
This file is part of OpenFOAM.
|
||||
|
||||
OpenFOAM is free software: you can redistribute it and/or modify it
|
||||
under the terms of the GNU General Public License as published by
|
||||
the Free Software Foundation, either version 3 of the License, or
|
||||
(at your option) any later version.
|
||||
|
||||
OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
|
||||
ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
|
||||
FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
|
||||
for more details.
|
||||
|
||||
You should have received a copy of the GNU General Public License
|
||||
along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
|
||||
|
||||
\*---------------------------------------------------------------------------*/
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
// * * * * * * * * * * * * * Private Member Functions * * * * * * * * * * * //
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * * Destructor * * * * * * * * * * * * * * * //
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * Member Operators * * * * * * * * * * * * * //
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * Friend Functions * * * * * * * * * * * * * //
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * Friend Operators * * * * * * * * * * * * * //
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * IOstream Operators * * * * * * * * * * * * * //
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * Ostream Operator * * * * * * * * * * * * * //
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
|
||||
// ************************************************************************* //
|
||||
|
|
@ -0,0 +1,65 @@
|
|||
/*---------------------------------------------------------------------------*\
|
||||
========= |
|
||||
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
||||
\\ / O peration |
|
||||
\\ / A nd | Copyright (C) 2016 OpenFOAM Foundation
|
||||
\\/ M anipulation |
|
||||
-------------------------------------------------------------------------------
|
||||
License
|
||||
This file is part of OpenFOAM.
|
||||
|
||||
OpenFOAM is free software: you can redistribute it and/or modify it
|
||||
under the terms of the GNU General Public License as published by
|
||||
the Free Software Foundation, either version 3 of the License, or
|
||||
(at your option) any later version.
|
||||
|
||||
OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
|
||||
ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
|
||||
FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
|
||||
for more details.
|
||||
|
||||
You should have received a copy of the GNU General Public License
|
||||
along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
|
||||
|
||||
\*---------------------------------------------------------------------------*/
|
||||
|
||||
#include "BetaGrid.H"
|
||||
#include "IOstreams.H"
|
||||
|
||||
// * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
|
||||
|
||||
Foam::CMC::BetaGrid::BetaGrid(Istream& is)
|
||||
{
|
||||
// Check state of Istream
|
||||
is.check("Foam::CMC::BetaGrid::BetaGrid(Foam::Istream&)");
|
||||
}
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * IOstream Operators * * * * * * * * * * * * //
|
||||
|
||||
Foam::Istream& Foam::operator>>(Istream& is, CMC::BetaGrid&)
|
||||
{
|
||||
// Check state of Istream
|
||||
is.check
|
||||
(
|
||||
"Foam::Istream& Foam::operator>>(Foam::Istream&, Foam::CMC::BetaGrid&)"
|
||||
);
|
||||
|
||||
return is;
|
||||
}
|
||||
|
||||
|
||||
Foam::Ostream& Foam::operator<<(Ostream& os, const CMC::BetaGrid&)
|
||||
{
|
||||
// Check state of Ostream
|
||||
os.check
|
||||
(
|
||||
"Foam::Ostream& Foam::operator<<(Foam::Ostream&, "
|
||||
"const Foam::CMC::BetaGrid&)"
|
||||
);
|
||||
|
||||
return os;
|
||||
}
|
||||
|
||||
|
||||
// ************************************************************************* //
|
||||
|
|
@ -0,0 +1,180 @@
|
|||
/*---------------------------------------------------------------------------*\
|
||||
========= |
|
||||
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
||||
\\ / O peration |
|
||||
\\ / A nd | Copyright (C) 2016 OpenFOAM Foundation
|
||||
\\/ M anipulation |
|
||||
-------------------------------------------------------------------------------
|
||||
License
|
||||
This file is part of OpenFOAM.
|
||||
|
||||
OpenFOAM is free software: you can redistribute it and/or modify it
|
||||
under the terms of the GNU General Public License as published by
|
||||
the Free Software Foundation, either version 3 of the License, or
|
||||
(at your option) any later version.
|
||||
|
||||
OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
|
||||
ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
|
||||
FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
|
||||
for more details.
|
||||
|
||||
You should have received a copy of the GNU General Public License
|
||||
along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
|
||||
|
||||
\*---------------------------------------------------------------------------*/
|
||||
|
||||
#include "BetaIntegrator.H"
|
||||
|
||||
|
||||
#include "error.H"
|
||||
#include "SubList.H"
|
||||
#include "SubField.H"
|
||||
#include "interpolateXY.H"
|
||||
|
||||
// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
|
||||
|
||||
// const scalar Foam::CMC::BetaIntegrator::staticData();
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * Static Member Functions * * * * * * * * * * * * //
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * Private Member Functions * * * * * * * * * * * //
|
||||
|
||||
//extended Simpson's rule (Numerical recipes, 2nd Ed. p.128)
|
||||
//for equally spaced and even N intervals (or odd N+1 points)
|
||||
Foam::scalar Foam::CMC::BetaIntegrator::simps
|
||||
(const scalar xl, const scalar xh, const label n, const UList<scalar>& fx)
|
||||
const
|
||||
{
|
||||
scalar evensum(0.0), oddsum(0.0), sum(0.0);
|
||||
|
||||
scalar h = (xh - xl)/scalar(n);
|
||||
|
||||
for(label i=0 ; i<fx.size() ; i++)
|
||||
{
|
||||
if(i%2 == 0)
|
||||
{
|
||||
evensum += fx[i];
|
||||
}
|
||||
else
|
||||
{
|
||||
oddsum += fx[i];
|
||||
}
|
||||
}
|
||||
|
||||
sum = -1.0*fx.first() + 2.0*evensum + 4.0*oddsum - 1.0*fx.last();
|
||||
|
||||
return sum*(h/3.0);
|
||||
}
|
||||
|
||||
|
||||
Foam::scalar Foam::CMC::BetaIntegrator::sumSimps(const UList<scalar>& f) const
|
||||
{
|
||||
scalar total = 0.0;
|
||||
|
||||
const labelList &N(bg_.N());
|
||||
const scalarField &etaCut(bg_.etaCut());
|
||||
|
||||
forAll(N, i)
|
||||
{
|
||||
const labelList::subList prev(N, i);
|
||||
const label baseI = sum(prev) + 1;
|
||||
|
||||
const scalarField::subField subInterval(f, N[i]+1, baseI);
|
||||
|
||||
total += simps(etaCut[i], etaCut[i+1], N[i], subInterval);
|
||||
}
|
||||
|
||||
return total;
|
||||
}
|
||||
|
||||
// * * * * * * * * * * * * Protected Member Functions * * * * * * * * * * * //
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
|
||||
|
||||
Foam::CMC::BetaIntegrator::BetaIntegrator(const BetaFunction& bf, const BetaGrid &bg)
|
||||
:
|
||||
etaSpace_(bg.etaSpace()),
|
||||
pdfNum_(bg.etaSpace().size(), 0.0),
|
||||
pdfDen_(0.0),
|
||||
bf_(bf),
|
||||
bg_(bg)
|
||||
{
|
||||
const scalarField::subField domain(etaSpace_, etaSpace_.size()-2, 1);
|
||||
scalarField::subField range(pdfNum_, etaSpace_.size()-2, 1);
|
||||
|
||||
if (bf_.fdelta())
|
||||
{
|
||||
pdfDen_ = 1.0;
|
||||
}
|
||||
else
|
||||
{
|
||||
range = bf_.etaFunc(domain);
|
||||
|
||||
if (!bf_.delta0())
|
||||
{
|
||||
pdfNum_.first() = bf_.etaFunc(etaSpace_.first());
|
||||
}
|
||||
|
||||
if (!bf_.delta1())
|
||||
{
|
||||
pdfNum_.last() = bf_.etaFunc(etaSpace_.last());
|
||||
}
|
||||
|
||||
pdfDen_ = sumSimps(pdfNum_)
|
||||
+ Foam::pow(etaSpace_[1], bf_.alpha())/(bf_.alpha() + SMALL)
|
||||
+ Foam::pow(etaSpace_[1], bf_.beta())/(bf_.beta() + SMALL);
|
||||
}
|
||||
}
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * * Destructor * * * * * * * * * * * * * * * //
|
||||
|
||||
Foam::CMC::BetaIntegrator::~BetaIntegrator()
|
||||
{}
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * * //
|
||||
|
||||
// beta-pdf weighted integration for given mf, mfVar and f
|
||||
Foam::scalar Foam::CMC::BetaIntegrator::betaIntegrate(const scalarField& f) const
|
||||
{
|
||||
scalar result = 0.0;
|
||||
|
||||
if(bf_.fdelta())
|
||||
{
|
||||
result = interpolateXY(bf_.mf(), etaSpace_, f);
|
||||
}
|
||||
else
|
||||
{
|
||||
scalar num = sumSimps (f * pdfNum_)
|
||||
+ f.first()*Foam::pow(etaSpace_[1], bf_.alpha())/(bf_.alpha() + SMALL)
|
||||
+ f.last()*Foam::pow(etaSpace_[1], bf_.beta())/(bf_.beta() + SMALL);
|
||||
result = num/pdfDen_;
|
||||
}
|
||||
|
||||
return result;
|
||||
}
|
||||
|
||||
// * * * * * * * * * * * * * * Member Operators * * * * * * * * * * * * * * //
|
||||
|
||||
void Foam::CMC::BetaIntegrator::operator=(const BetaIntegrator& rhs)
|
||||
{
|
||||
// Check for assignment to self
|
||||
if (this == &rhs)
|
||||
{
|
||||
FatalErrorInFunction
|
||||
<< "Attempted assignment to self"
|
||||
<< abort(FatalError);
|
||||
}
|
||||
}
|
||||
|
||||
// * * * * * * * * * * * * * * Friend Functions * * * * * * * * * * * * * * //
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * Friend Operators * * * * * * * * * * * * * * //
|
||||
|
||||
|
||||
// ************************************************************************* //
|
||||
|
|
@ -0,0 +1,168 @@
|
|||
/*---------------------------------------------------------------------------*\
|
||||
========= |
|
||||
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
||||
\\ / O peration |
|
||||
\\ / A nd | Copyright (C) 2016 OpenFOAM Foundation
|
||||
\\/ M anipulation |
|
||||
-------------------------------------------------------------------------------
|
||||
License
|
||||
This file is part of OpenFOAM.
|
||||
|
||||
OpenFOAM is free software: you can redistribute it and/or modify it
|
||||
under the terms of the GNU General Public License as published by
|
||||
the Free Software Foundation, either version 3 of the License, or
|
||||
(at your option) any later version.
|
||||
|
||||
OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
|
||||
ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
|
||||
FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
|
||||
for more details.
|
||||
|
||||
You should have received a copy of the GNU General Public License
|
||||
along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
|
||||
|
||||
Class
|
||||
Foam::CMC::BetaIntegrator
|
||||
|
||||
Description
|
||||
|
||||
SourceFiles
|
||||
BetaIntegratorI.H
|
||||
BetaIntegrator.C
|
||||
BetaIntegratorIO.C
|
||||
|
||||
\*---------------------------------------------------------------------------*/
|
||||
|
||||
#ifndef BetaIntegrator_H
|
||||
#define BetaIntegrator_H
|
||||
|
||||
#include "scalarField.H"
|
||||
#include "BetaFunction.H"
|
||||
#include "BetaGrid.H"
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
namespace Foam
|
||||
{
|
||||
|
||||
// Forward declaration of classes
|
||||
class Istream;
|
||||
class Ostream;
|
||||
|
||||
// Forward declaration of friend functions and operators
|
||||
namespace CMC
|
||||
{
|
||||
class BetaIntegrator;
|
||||
}
|
||||
Istream& operator>>(Istream&, CMC::BetaIntegrator&);
|
||||
Ostream& operator<<(Ostream&, const CMC::BetaIntegrator&);
|
||||
|
||||
|
||||
/*---------------------------------------------------------------------------*\
|
||||
Class BetaIntegrator Declaration
|
||||
\*---------------------------------------------------------------------------*/
|
||||
|
||||
namespace CMC
|
||||
{
|
||||
|
||||
class BetaIntegrator
|
||||
{
|
||||
// Private data
|
||||
|
||||
//- Quadrature points
|
||||
const scalarField &etaSpace_;
|
||||
|
||||
//- Beta pdf numerators
|
||||
scalarField pdfNum_;
|
||||
|
||||
//- Beta pdf denominator
|
||||
scalar pdfDen_;
|
||||
|
||||
//- Description of data_
|
||||
BetaFunction bf_;
|
||||
|
||||
//- Description of data_
|
||||
const BetaGrid &bg_;
|
||||
|
||||
|
||||
// Private Member Functions
|
||||
|
||||
//- Calculate integral using Simpson's rule
|
||||
scalar simps(const scalar xl, const scalar xh, const label n, const UList<scalar>& fx) const;
|
||||
|
||||
//- Calculate integral using Simpson's rule
|
||||
scalar sumSimps(const UList<scalar>& f) const;
|
||||
|
||||
//- Disallow default bitwise copy construct
|
||||
BetaIntegrator(const BetaIntegrator&);
|
||||
|
||||
//- Disallow default bitwise assignment
|
||||
void operator=(const BetaIntegrator&);
|
||||
|
||||
|
||||
public:
|
||||
|
||||
// Static data members
|
||||
|
||||
//- Static data staticData
|
||||
// static const scalar staticData;
|
||||
|
||||
|
||||
// Constructors
|
||||
|
||||
//- Construct from components
|
||||
BetaIntegrator(const BetaFunction& bf, const BetaGrid &bg);
|
||||
|
||||
//- Construct as copy
|
||||
// BetaIntegrator(const BetaIntegrator&);
|
||||
|
||||
|
||||
//- Destructor
|
||||
~BetaIntegrator();
|
||||
|
||||
|
||||
// Member Functions
|
||||
|
||||
scalar betaIntegrate(const scalarField& f) const;
|
||||
|
||||
// Access
|
||||
|
||||
scalar pdfDen () const {return pdfDen_;}
|
||||
|
||||
// Check
|
||||
|
||||
// Edit
|
||||
|
||||
// Write
|
||||
|
||||
|
||||
// Member Operators
|
||||
|
||||
// void operator=(const BetaIntegrator&);
|
||||
|
||||
|
||||
// Friend Functions
|
||||
|
||||
// Friend Operators
|
||||
|
||||
// IOstream Operators
|
||||
|
||||
friend Istream& ::Foam::operator>>(Istream&, BetaIntegrator&);
|
||||
friend Ostream& ::Foam::operator<<(Ostream&, const BetaIntegrator&);
|
||||
};
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
} // End namespace CMC
|
||||
} // End namespace Foam
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
#include "BetaIntegratorI.H"
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
#endif
|
||||
|
||||
// ************************************************************************* //
|
||||
|
|
@ -0,0 +1,58 @@
|
|||
/*---------------------------------------------------------------------------*\
|
||||
========= |
|
||||
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
||||
\\ / O peration |
|
||||
\\ / A nd | Copyright (C) 2016 OpenFOAM Foundation
|
||||
\\/ M anipulation |
|
||||
-------------------------------------------------------------------------------
|
||||
License
|
||||
This file is part of OpenFOAM.
|
||||
|
||||
OpenFOAM is free software: you can redistribute it and/or modify it
|
||||
under the terms of the GNU General Public License as published by
|
||||
the Free Software Foundation, either version 3 of the License, or
|
||||
(at your option) any later version.
|
||||
|
||||
OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
|
||||
ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
|
||||
FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
|
||||
for more details.
|
||||
|
||||
You should have received a copy of the GNU General Public License
|
||||
along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
|
||||
|
||||
\*---------------------------------------------------------------------------*/
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
// * * * * * * * * * * * * * Private Member Functions * * * * * * * * * * * //
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * * Destructor * * * * * * * * * * * * * * * //
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * Member Operators * * * * * * * * * * * * * //
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * Friend Functions * * * * * * * * * * * * * //
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * Friend Operators * * * * * * * * * * * * * //
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * IOstream Operators * * * * * * * * * * * * * //
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * Ostream Operator * * * * * * * * * * * * * //
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
|
||||
// ************************************************************************* //
|
||||
|
|
@ -0,0 +1,56 @@
|
|||
/*---------------------------------------------------------------------------*\
|
||||
========= |
|
||||
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
||||
\\ / O peration |
|
||||
\\ / A nd | Copyright (C) 2016 OpenFOAM Foundation
|
||||
\\/ M anipulation |
|
||||
-------------------------------------------------------------------------------
|
||||
License
|
||||
This file is part of OpenFOAM.
|
||||
|
||||
OpenFOAM is free software: you can redistribute it and/or modify it
|
||||
under the terms of the GNU General Public License as published by
|
||||
the Free Software Foundation, either version 3 of the License, or
|
||||
(at your option) any later version.
|
||||
|
||||
OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
|
||||
ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
|
||||
FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
|
||||
for more details.
|
||||
|
||||
You should have received a copy of the GNU General Public License
|
||||
along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
|
||||
|
||||
\*---------------------------------------------------------------------------*/
|
||||
|
||||
#include "BetaIntegrator.H"
|
||||
#include "IOstreams.H"
|
||||
|
||||
// * * * * * * * * * * * * * * * IOstream Operators * * * * * * * * * * * * //
|
||||
|
||||
Foam::Istream& Foam::operator>>(Istream& is, CMC::BetaIntegrator&)
|
||||
{
|
||||
// Check state of Istream
|
||||
is.check
|
||||
(
|
||||
"Foam::Istream& Foam::operator>>(Foam::Istream&, Foam::CMC::BetaIntegrator&)"
|
||||
);
|
||||
|
||||
return is;
|
||||
}
|
||||
|
||||
|
||||
Foam::Ostream& Foam::operator<<(Ostream& os, const CMC::BetaIntegrator&)
|
||||
{
|
||||
// Check state of Ostream
|
||||
os.check
|
||||
(
|
||||
"Foam::Ostream& Foam::operator<<(Foam::Ostream&, "
|
||||
"const Foam::CMC::BetaIntegrator&)"
|
||||
);
|
||||
|
||||
return os;
|
||||
}
|
||||
|
||||
|
||||
// ************************************************************************* //
|
||||
|
|
@ -0,0 +1,206 @@
|
|||
/*---------------------------------------------------------------------------*\
|
||||
========= |
|
||||
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
||||
\\ / O peration |
|
||||
\\ / A nd | Copyright (C) 2016 OpenFOAM Foundation
|
||||
\\/ M anipulation |
|
||||
-------------------------------------------------------------------------------
|
||||
License
|
||||
This file is part of OpenFOAM.
|
||||
|
||||
OpenFOAM is free software: you can redistribute it and/or modify it
|
||||
under the terms of the GNU General Public License as published by
|
||||
the Free Software Foundation, either version 3 of the License, or
|
||||
(at your option) any later version.
|
||||
|
||||
OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
|
||||
ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
|
||||
FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
|
||||
for more details.
|
||||
|
||||
You should have received a copy of the GNU General Public License
|
||||
along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
|
||||
|
||||
\*---------------------------------------------------------------------------*/
|
||||
|
||||
#include "FlameStructure.H"
|
||||
|
||||
#include "fileName.H"
|
||||
#include "IFstream.H"
|
||||
#include "hashedWordList.H"
|
||||
#include "interpolateXY.H"
|
||||
|
||||
// * * * * * * * * * * * * * Private Member Functions * * * * * * * * * * * //
|
||||
|
||||
|
||||
// * * * * * * * * * * * * Protected Member Functions * * * * * * * * * * * //
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
|
||||
|
||||
Foam::SLFM::FlameStructure::FlameStructure(const string& NstDir, scalarList NstList, const scalarField &etaGrid, const hashedWordList &spTable)
|
||||
:
|
||||
NstDir_(NstDir),
|
||||
NstList_(NstList),
|
||||
eta_(etaGrid),
|
||||
species_(spTable),
|
||||
Y_(spTable.size(), scalarFieldArray1d(NstList_.size(), scalarField(eta_.size(), 0.0))),
|
||||
W_(spTable.size(), scalarFieldArray1d(NstList_.size(), scalarField(eta_.size(), 0.0))),
|
||||
T_(NstList_.size(), scalarField(eta_.size(), 0.0)),
|
||||
Q_(NstList_.size(), scalarField(eta_.size(), 0.0)),
|
||||
h_(NstList_.size(), scalarField(eta_.size(), 0.0)),
|
||||
rho_(NstList_.size(), scalarField(eta_.size(), 0.0)),
|
||||
Rgas_(NstList_.size(), scalarField(eta_.size(), 0.0))
|
||||
{
|
||||
for (label n = 0; n < NstList_.size(); n++)
|
||||
{
|
||||
scalar N = NstList_[n];
|
||||
|
||||
fileName fname = "sdr"+Foam::name(N)+".inp";
|
||||
IFstream fin(NstDir_/fname);
|
||||
|
||||
string gbg;
|
||||
label etamax_SLFM, NoSpecies;
|
||||
|
||||
//INPUT FILE...
|
||||
fin.getLine(gbg);
|
||||
|
||||
//No.GRID...
|
||||
fin.getLine(gbg);
|
||||
|
||||
//Number of eta-point and species
|
||||
fin >> etamax_SLFM >> NoSpecies;
|
||||
|
||||
//blank line
|
||||
fin.getLine(gbg);
|
||||
|
||||
//PRESSURE...
|
||||
fin.getLine(gbg);
|
||||
|
||||
//1.00...
|
||||
fin.getLine(gbg);
|
||||
|
||||
//MIXTURE...
|
||||
fin.getLine(gbg);
|
||||
|
||||
//read eta space (from SLFM library)
|
||||
scalarField etaValue_SLFM(etamax_SLFM, 0.0);
|
||||
for(label j = 0 ; j < etamax_SLFM ; j++)
|
||||
{
|
||||
fin>>etaValue_SLFM[j];
|
||||
}
|
||||
|
||||
//blank line
|
||||
fin.getLine(gbg);
|
||||
|
||||
//INITIAL...
|
||||
fin.getLine(gbg);
|
||||
|
||||
//species loop
|
||||
scalarField temp(etamax_SLFM, 0.0);
|
||||
for(label i = 0 ; i < NoSpecies ; i++)
|
||||
{
|
||||
word spName;
|
||||
//H2... (species name)
|
||||
fin >> spName;
|
||||
const label spI (species_[spName]);
|
||||
|
||||
// Info << spI << spName << endl;
|
||||
|
||||
//read species mass fraction (from SLFM library)
|
||||
for(label j=0 ; j<etamax_SLFM ; j++)
|
||||
{
|
||||
scalar Yj = 0.0;
|
||||
fin >> Yj;
|
||||
temp[j] = max(0.0, Yj);
|
||||
}
|
||||
|
||||
Y_[spI][n] = interpolateXY(etaGrid, etaValue_SLFM, temp);
|
||||
|
||||
//read reaction rate (from SLFM library)
|
||||
for(label j = 0 ; j < etamax_SLFM ; j++)
|
||||
{
|
||||
fin>>temp[j];
|
||||
}
|
||||
|
||||
W_[spI][n] = interpolateXY(etaGrid, etaValue_SLFM, temp);
|
||||
|
||||
//blank line
|
||||
fin.getLine(gbg);
|
||||
}
|
||||
|
||||
//read temperature
|
||||
{
|
||||
//H2... (species name)
|
||||
fin.getLine(gbg);
|
||||
|
||||
//read temperature
|
||||
for(label j = 0 ; j < etamax_SLFM ; j++)
|
||||
{
|
||||
scalar Yj = 0.0;
|
||||
fin >> Yj;
|
||||
temp[j] = max(0.0, Yj);
|
||||
}
|
||||
|
||||
T_[n] = interpolateXY(etaGrid, etaValue_SLFM, temp);
|
||||
|
||||
//read heat source
|
||||
for(label j = 0 ; j < etamax_SLFM ; j++)
|
||||
{
|
||||
fin>>temp[j];
|
||||
}
|
||||
|
||||
Q_[n] = interpolateXY(etaGrid, etaValue_SLFM, temp);
|
||||
|
||||
//blank line
|
||||
fin.getLine(gbg);
|
||||
}
|
||||
|
||||
//DENSITY...
|
||||
fin.getLine(gbg);
|
||||
|
||||
//read density (from SLFM library)
|
||||
for(label j = 0 ; j < etamax_SLFM ; j++)
|
||||
{
|
||||
fin>>temp[j];
|
||||
}
|
||||
//Info<<rho_temp<<endl;
|
||||
|
||||
rho_[n] = interpolateXY(etaGrid, etaValue_SLFM, temp);
|
||||
|
||||
//blank line
|
||||
fin.getLine(gbg);
|
||||
|
||||
//ENTHALPY...
|
||||
fin.getLine(gbg);
|
||||
|
||||
//read enthalpy (from SLFM library)
|
||||
for(label j = 0 ; j < etamax_SLFM ; j++)
|
||||
{
|
||||
fin>>temp[j];
|
||||
}
|
||||
|
||||
h_[n] = interpolateXY(etaGrid, etaValue_SLFM, temp);
|
||||
}
|
||||
}
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * * Destructor * * * * * * * * * * * * * * * //
|
||||
|
||||
Foam::SLFM::FlameStructure::~FlameStructure()
|
||||
{}
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * * //
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * Member Operators * * * * * * * * * * * * * * //
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * Friend Functions * * * * * * * * * * * * * * //
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * Friend Operators * * * * * * * * * * * * * * //
|
||||
|
||||
|
||||
// ************************************************************************* //
|
||||
|
|
@ -0,0 +1,215 @@
|
|||
/*---------------------------------------------------------------------------*\
|
||||
========= |
|
||||
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
||||
\\ / O peration |
|
||||
\\ / A nd | Copyright (C) 2016 OpenFOAM Foundation
|
||||
\\/ M anipulation |
|
||||
-------------------------------------------------------------------------------
|
||||
License
|
||||
This file is part of OpenFOAM.
|
||||
|
||||
OpenFOAM is free software: you can redistribute it and/or modify it
|
||||
under the terms of the GNU General Public License as published by
|
||||
the Free Software Foundation, either version 3 of the License, or
|
||||
(at your option) any later version.
|
||||
|
||||
OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
|
||||
ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
|
||||
FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
|
||||
for more details.
|
||||
|
||||
You should have received a copy of the GNU General Public License
|
||||
along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
|
||||
|
||||
Class
|
||||
Foam::SLFM::FlameStructure
|
||||
|
||||
Description
|
||||
|
||||
SourceFiles
|
||||
FlameStructureI.H
|
||||
FlameStructure.C
|
||||
FlameStructureIO.C
|
||||
|
||||
\*---------------------------------------------------------------------------*/
|
||||
|
||||
#ifndef FlameStructure_H
|
||||
#define FlameStructure_H
|
||||
|
||||
#include "error.H"
|
||||
#include "scalarField.H"
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
namespace Foam
|
||||
{
|
||||
|
||||
// Forward declaration of classes
|
||||
class Istream;
|
||||
class Ostream;
|
||||
|
||||
class hashedWordList;
|
||||
|
||||
namespace SLFM
|
||||
{
|
||||
class FlameStructure;
|
||||
}
|
||||
|
||||
// Forward declaration of friend functions and operators
|
||||
Istream& operator>>(Istream&, SLFM::FlameStructure&);
|
||||
Ostream& operator<<(Ostream&, const SLFM::FlameStructure&);
|
||||
|
||||
namespace SLFM
|
||||
{
|
||||
|
||||
/*---------------------------------------------------------------------------*\
|
||||
Class FlameStructure Declaration
|
||||
\*---------------------------------------------------------------------------*/
|
||||
|
||||
class FlameStructure
|
||||
{
|
||||
// Private data
|
||||
typedef List<scalarField> scalarFieldArray1d;
|
||||
|
||||
typedef List<scalarFieldArray1d> scalarFieldArray2d;
|
||||
|
||||
typedef List<scalarFieldArray2d> scalarFieldArray3d;
|
||||
|
||||
//- Species mass fraction
|
||||
word NstDir_;
|
||||
|
||||
//- Species mass fraction
|
||||
const scalarList NstList_;
|
||||
|
||||
//- Species mass fraction
|
||||
const scalarField eta_;
|
||||
|
||||
//- Species mass fraction
|
||||
const hashedWordList &species_;
|
||||
|
||||
//- Species mass fraction
|
||||
scalarFieldArray2d Y_;
|
||||
|
||||
//- Species production rate
|
||||
scalarFieldArray2d W_;
|
||||
|
||||
//- Temperature
|
||||
scalarFieldArray1d T_;
|
||||
|
||||
//- Heat source
|
||||
scalarFieldArray1d Q_;
|
||||
|
||||
//- Enthalpy
|
||||
scalarFieldArray1d h_;
|
||||
|
||||
//- Density
|
||||
scalarFieldArray1d rho_;
|
||||
|
||||
//- Specific gas constant for mean W
|
||||
scalarFieldArray1d Rgas_;
|
||||
|
||||
|
||||
// Private Member Functions
|
||||
|
||||
//- Disallow Construct null
|
||||
FlameStructure();
|
||||
|
||||
//- Disallow default bitwise copy construct
|
||||
// FlameStructure(const FlameStructure&);
|
||||
|
||||
//- Disallow default bitwise assignment
|
||||
// void operator=(const FlameStructure&);
|
||||
|
||||
|
||||
public:
|
||||
|
||||
// Static data members
|
||||
|
||||
//- Static data staticData
|
||||
// static const dataType staticData;
|
||||
|
||||
|
||||
// Constructors
|
||||
|
||||
//- Construct from components
|
||||
FlameStructure(const string& NstDir, scalarList NstList, const scalarField &etaGrid, const hashedWordList& spTable);
|
||||
|
||||
//- Construct from Istream
|
||||
// FlameStructure(Istream&);
|
||||
|
||||
//- Construct as copy
|
||||
// FlameStructure(const FlameStructure&);
|
||||
|
||||
|
||||
// Selectors
|
||||
|
||||
//- Select null constructed
|
||||
// static autoPtr<FlameStructure> New();
|
||||
|
||||
|
||||
//- Destructor
|
||||
~FlameStructure();
|
||||
|
||||
|
||||
// Member Functions
|
||||
|
||||
// Access
|
||||
|
||||
//- Select null constructed
|
||||
scalarFieldArray2d& Y() {return Y_;}
|
||||
|
||||
//- Select null constructed
|
||||
scalarFieldArray2d& W() {return W_;}
|
||||
|
||||
//- Select null constructed
|
||||
scalarFieldArray1d& T() {return T_;}
|
||||
|
||||
//- Select null constructed
|
||||
scalarFieldArray1d& Q() {return Q_;}
|
||||
|
||||
//- Select null constructed
|
||||
scalarFieldArray1d& rho() {return rho_;}
|
||||
|
||||
//- Select null constructed
|
||||
scalarFieldArray1d& h() {return h_;}
|
||||
|
||||
//- Select null constructed
|
||||
scalarFieldArray1d& Rgas() {return Rgas_;}
|
||||
|
||||
// Check
|
||||
|
||||
// Edit
|
||||
|
||||
// Write
|
||||
|
||||
|
||||
// Member Operators
|
||||
|
||||
// void operator=(const FlameStructure&);
|
||||
|
||||
|
||||
// Friend Functions
|
||||
|
||||
// Friend Operators
|
||||
|
||||
// IOstream Operators
|
||||
|
||||
friend Istream& operator>>(Istream&, FlameStructure&);
|
||||
friend Ostream& operator<<(Ostream&, const FlameStructure&);
|
||||
};
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
} // End namespace SLFM
|
||||
} // End namespace Foam
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
#include "FlameStructureI.H"
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
#endif
|
||||
|
||||
// ************************************************************************* //
|
||||
|
|
@ -0,0 +1,58 @@
|
|||
/*---------------------------------------------------------------------------*\
|
||||
========= |
|
||||
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
||||
\\ / O peration |
|
||||
\\ / A nd | Copyright (C) 2016 OpenFOAM Foundation
|
||||
\\/ M anipulation |
|
||||
-------------------------------------------------------------------------------
|
||||
License
|
||||
This file is part of OpenFOAM.
|
||||
|
||||
OpenFOAM is free software: you can redistribute it and/or modify it
|
||||
under the terms of the GNU General Public License as published by
|
||||
the Free Software Foundation, either version 3 of the License, or
|
||||
(at your option) any later version.
|
||||
|
||||
OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
|
||||
ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
|
||||
FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
|
||||
for more details.
|
||||
|
||||
You should have received a copy of the GNU General Public License
|
||||
along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
|
||||
|
||||
\*---------------------------------------------------------------------------*/
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
// * * * * * * * * * * * * * Private Member Functions * * * * * * * * * * * //
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * * Destructor * * * * * * * * * * * * * * * //
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * Member Operators * * * * * * * * * * * * * //
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * Friend Functions * * * * * * * * * * * * * //
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * Friend Operators * * * * * * * * * * * * * //
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * IOstream Operators * * * * * * * * * * * * * //
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * Ostream Operator * * * * * * * * * * * * * //
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
|
||||
// ************************************************************************* //
|
||||
|
|
@ -0,0 +1,72 @@
|
|||
/*---------------------------------------------------------------------------*\
|
||||
========= |
|
||||
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
||||
\\ / O peration |
|
||||
\\ / A nd | Copyright (C) 2016 OpenFOAM Foundation
|
||||
\\/ M anipulation |
|
||||
-------------------------------------------------------------------------------
|
||||
License
|
||||
This file is part of OpenFOAM.
|
||||
|
||||
OpenFOAM is free software: you can redistribute it and/or modify it
|
||||
under the terms of the GNU General Public License as published by
|
||||
the Free Software Foundation, either version 3 of the License, or
|
||||
(at your option) any later version.
|
||||
|
||||
OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
|
||||
ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
|
||||
FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
|
||||
for more details.
|
||||
|
||||
You should have received a copy of the GNU General Public License
|
||||
along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
|
||||
|
||||
\*---------------------------------------------------------------------------*/
|
||||
|
||||
#include "FlameStructure.H"
|
||||
#include "IOstreams.H"
|
||||
|
||||
// * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
|
||||
|
||||
/*
|
||||
Foam::SLFM::FlameStructure::FlameStructure(Istream& is)
|
||||
:
|
||||
base1(is),
|
||||
base2(is),
|
||||
member1(is),
|
||||
member2(is)
|
||||
{
|
||||
// Check state of Istream
|
||||
is.check("Foam::SLFM::FlameStructure::FlameStructure(Foam::Istream&)");
|
||||
}
|
||||
*/
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * IOstream Operators * * * * * * * * * * * * //
|
||||
|
||||
Foam::Istream& Foam::operator>>(Istream& is, SLFM::FlameStructure&)
|
||||
{
|
||||
// Check state of Istream
|
||||
is.check
|
||||
(
|
||||
"Foam::Istream& Foam::operator>>(Foam::Istream&, Foam::SLFM::FlameStructure&)"
|
||||
);
|
||||
|
||||
return is;
|
||||
}
|
||||
|
||||
|
||||
Foam::Ostream& Foam::operator<<(Ostream& os, const SLFM::FlameStructure&)
|
||||
{
|
||||
// Check state of Ostream
|
||||
os.check
|
||||
(
|
||||
"Foam::Ostream& Foam::operator<<(Foam::Ostream&, "
|
||||
"const Foam::SLFM::FlameStructure&)"
|
||||
);
|
||||
|
||||
return os;
|
||||
}
|
||||
|
||||
|
||||
// ************************************************************************* //
|
||||
17
applications/solvers/combustion/SLFMFoam/Make/files
Normal file
17
applications/solvers/combustion/SLFMFoam/Make/files
Normal file
|
|
@ -0,0 +1,17 @@
|
|||
FlameStructure/FlameStructure.C
|
||||
FlameStructure/FlameStructureIO.C
|
||||
|
||||
./AMC/AMC.C
|
||||
|
||||
./BetaFunction/BetaFunction.C
|
||||
./BetaFunction/BetaFunctionIO.C
|
||||
|
||||
./BetaGrid/BetaGrid.C
|
||||
./BetaGrid/BetaGridIO.C
|
||||
|
||||
./BetaIntegrator/BetaIntegrator.C
|
||||
./BetaIntegrator/BetaIntegratorIO.C
|
||||
|
||||
SLFMFoam.C
|
||||
|
||||
EXE = $(FOAM_APPBIN)/SLFMFoam
|
||||
34
applications/solvers/combustion/SLFMFoam/Make/options
Normal file
34
applications/solvers/combustion/SLFMFoam/Make/options
Normal file
|
|
@ -0,0 +1,34 @@
|
|||
EXE_INC = \
|
||||
-I. \
|
||||
-I./AMC \
|
||||
-I./BetaFunction \
|
||||
-I./BetaGrid \
|
||||
-I./BetaIntegrator \
|
||||
-I$(LIB_SRC)/finiteVolume/lnInclude \
|
||||
-I$(LIB_SRC)/meshTools/lnInclude \
|
||||
-I$(LIB_SRC)/TurbulenceModels/turbulenceModels/lnInclude \
|
||||
-I$(LIB_SRC)/TurbulenceModels/compressible/lnInclude \
|
||||
-I$(LIB_SRC)/thermophysicalModels/specie/lnInclude \
|
||||
-I$(LIB_SRC)/transportModels/compressible/lnInclude \
|
||||
-I$(LIB_SRC)/thermophysicalModels/basic/lnInclude \
|
||||
-I$(LIB_SRC)/thermophysicalModels/reactionThermo/lnInclude \
|
||||
-I$(LIB_SRC)/thermophysicalModels/chemistryModel/lnInclude \
|
||||
-I$(LIB_SRC)/ODE/lnInclude \
|
||||
-I$(LIB_SRC)/combustionModels/lnInclude \
|
||||
-I$(LIB_SRC)/sampling/lnInclude
|
||||
|
||||
|
||||
EXE_LIBS = \
|
||||
-lfiniteVolume \
|
||||
-lmeshTools \
|
||||
-lturbulenceModels \
|
||||
-lcompressibleTurbulenceModels \
|
||||
-lspecie \
|
||||
-lcompressibleTransportModels \
|
||||
-lfluidThermophysicalModels \
|
||||
-lreactionThermophysicalModels \
|
||||
-lchemistryModel \
|
||||
-lODE \
|
||||
-lcombustionModels \
|
||||
-lfvOptions \
|
||||
-lsampling
|
||||
10
applications/solvers/combustion/SLFMFoam/Mixturefraction.H
Normal file
10
applications/solvers/combustion/SLFMFoam/Mixturefraction.H
Normal file
|
|
@ -0,0 +1,10 @@
|
|||
// mixture fraction equation
|
||||
fvScalarMatrix mfEqn
|
||||
(
|
||||
fvm::div(phi, mf)
|
||||
- fvm::laplacian((1/Sc)*turbulence->mut(), mf) //let 1/Sc = 1.47
|
||||
);
|
||||
|
||||
mfEqn.relax();
|
||||
|
||||
mfEqn.solve();
|
||||
|
|
@ -0,0 +1,17 @@
|
|||
//mixture fraction variance equation
|
||||
//cf) mfVar should not be ZERO
|
||||
SDR = turbulence->epsilon() * mfVar / turbulence->k();
|
||||
volVectorField Gradmf = fvc::grad(mf);
|
||||
|
||||
fvScalarMatrix mfVarEqn
|
||||
(
|
||||
fvm::div(phi, mfVar)
|
||||
- fvm::laplacian((1/Sc) * turbulence->mut(), mfVar)
|
||||
==
|
||||
2 * ((1/Sc) * turbulence->mut()) * (Gradmf & Gradmf)
|
||||
- 2 * rho * SDR
|
||||
);
|
||||
|
||||
mfVarEqn.relax();
|
||||
|
||||
mfVarEqn.solve();
|
||||
118
applications/solvers/combustion/SLFMFoam/SLFMFoam.C
Normal file
118
applications/solvers/combustion/SLFMFoam/SLFMFoam.C
Normal file
|
|
@ -0,0 +1,118 @@
|
|||
/*---------------------------------------------------------------------------*\
|
||||
========= |
|
||||
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
||||
\\ / O peration |
|
||||
\\ / A nd | Copyright (C) 2013-2016 OpenFOAM Foundation
|
||||
\\/ M anipulation |
|
||||
-------------------------------------------------------------------------------
|
||||
License
|
||||
This file is part of OpenFOAM.
|
||||
|
||||
OpenFOAM is free software: you can redistribute it and/or modify it
|
||||
under the terms of the GNU General Public License as published by
|
||||
the Free Software Foundation, either version 3 of the License, or
|
||||
(at your option) any later version.
|
||||
|
||||
OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
|
||||
ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
|
||||
FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
|
||||
for more details.
|
||||
|
||||
You should have received a copy of the GNU General Public License
|
||||
along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
|
||||
|
||||
Application
|
||||
SLFMFoam
|
||||
|
||||
Description
|
||||
Steady Laminar Flamelet Model(SLFM) solver for turbulent combustion. SLFM
|
||||
assumes a turbulent flame composed of thin stretched laminar flamelets.
|
||||
Each flamelet library is generated by external program and imported by
|
||||
the solver.
|
||||
|
||||
References
|
||||
A.Y. Klimenko, R.W. Bilger, Progress in Energy and Combustion Science 25 (1999) 595-687
|
||||
N. Peters, Turbulent Combustion, Cambridge University Press (2000)
|
||||
|
||||
Contact
|
||||
POSTECH combustion lab.
|
||||
huh@postech.ac.kr
|
||||
|
||||
\*---------------------------------------------------------------------------*/
|
||||
|
||||
#include "fvCFD.H"
|
||||
#include "turbulentFluidThermoModel.H"
|
||||
#include "rhoCombustionModel.H"
|
||||
#include "fvOptions.H"
|
||||
#include "simpleControl.H"
|
||||
#include "FlameStructure/FlameStructure.H"
|
||||
#include "interpolateXY.H"
|
||||
|
||||
#include "AMC.H"
|
||||
#include "BetaFunction.H"
|
||||
#include "BetaGrid.H"
|
||||
#include "BetaIntegrator.H"
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
int main(int argc, char *argv[])
|
||||
{
|
||||
#include "postProcess.H"
|
||||
|
||||
#include "setRootCase.H"
|
||||
#include "createTime.H"
|
||||
#include "createMesh.H"
|
||||
#include "createControl.H"
|
||||
#include "createFields.H"
|
||||
#include "createFieldRefs.H"
|
||||
#include "createFvOptions.H"
|
||||
#include "initContinuityErrs.H"
|
||||
|
||||
turbulence->validate();
|
||||
|
||||
//SLFM
|
||||
Info<<"Read SLFM library"<<endl;
|
||||
#include "readSLFMlib.H" //read SLFM library
|
||||
|
||||
Info<<"Construct Beta-PDF"<<endl;
|
||||
#include "preIntegration.H" //AMC C1 coefficient pdf integration
|
||||
|
||||
Info<< "\nStarting time loop\n" << endl;
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
Info<< "\nStarting time loop\n" << endl;
|
||||
|
||||
while (simple.loop())
|
||||
{
|
||||
Info<< "Time = " << runTime.timeName() << nl << endl;
|
||||
|
||||
// --- Pressure-velocity SIMPLE corrector loop
|
||||
{
|
||||
#include "UEqn.H"
|
||||
#include "Mixturefraction.H" //mean mixture fraction
|
||||
#include "MixturefractionVar.H" //mean mixture fraction variance
|
||||
#include "setSDR.H" //calculate scalar dissipation rate
|
||||
#include "updateT_RHO.H" //update temperature and density
|
||||
#include "pEqn.H"
|
||||
}
|
||||
|
||||
turbulence->correct();
|
||||
|
||||
if (runTime.write())
|
||||
{
|
||||
#include "updateYi.H" //update species
|
||||
}
|
||||
|
||||
Info<< "ExecutionTime = " << runTime.elapsedCpuTime() << " s"
|
||||
<< " ClockTime = " << runTime.elapsedClockTime() << " s"
|
||||
<< nl << endl;
|
||||
}
|
||||
|
||||
Info<< "End\n" << endl;
|
||||
|
||||
return 0;
|
||||
}
|
||||
|
||||
|
||||
// ************************************************************************* //
|
||||
20
applications/solvers/combustion/SLFMFoam/UEqn.H
Normal file
20
applications/solvers/combustion/SLFMFoam/UEqn.H
Normal file
|
|
@ -0,0 +1,20 @@
|
|||
MRF.correctBoundaryVelocity(U);
|
||||
|
||||
tmp<fvVectorMatrix> tUEqn
|
||||
(
|
||||
fvm::div(phi, U)
|
||||
+ MRF.DDt(rho, U)
|
||||
+ turbulence->divDevRhoReff(U)
|
||||
==
|
||||
rho()*g
|
||||
+ fvOptions(rho, U)
|
||||
);
|
||||
fvVectorMatrix& UEqn = tUEqn.ref();
|
||||
|
||||
UEqn.relax();
|
||||
|
||||
fvOptions.constrain(UEqn);
|
||||
|
||||
solve(UEqn == -fvc::grad(p));
|
||||
|
||||
fvOptions.correct(U);
|
||||
|
|
@ -0,0 +1,2 @@
|
|||
const volScalarField& psi = thermo.psi();
|
||||
//const volScalarField& T = thermo.T();
|
||||
374
applications/solvers/combustion/SLFMFoam/createFields.H
Normal file
374
applications/solvers/combustion/SLFMFoam/createFields.H
Normal file
|
|
@ -0,0 +1,374 @@
|
|||
#include "readGravitationalAcceleration.H"
|
||||
|
||||
Info<< "Creating combustion model\n" << endl;
|
||||
|
||||
autoPtr<combustionModels::rhoCombustionModel> combustion
|
||||
(
|
||||
combustionModels::rhoCombustionModel::New(mesh)
|
||||
);
|
||||
|
||||
rhoReactionThermo& thermo = combustion->thermo();
|
||||
thermo.validate(args.executable(), "h", "e");
|
||||
|
||||
basicSpecieMixture& composition = thermo.composition();
|
||||
PtrList<volScalarField>& Y = composition.Y();
|
||||
|
||||
const word inertSpecie(thermo.lookup("inertSpecie"));
|
||||
|
||||
if (!composition.contains(inertSpecie))
|
||||
{
|
||||
FatalErrorInFunction
|
||||
<< "Specified inert specie '" << inertSpecie << "' not found in "
|
||||
<< "species list. Available species:" << composition.species()
|
||||
<< exit(FatalError);
|
||||
}
|
||||
|
||||
volScalarField& p = thermo.p();
|
||||
volScalarField& T = thermo.T();
|
||||
|
||||
wordList rhoBoundaryTypes
|
||||
(
|
||||
T.boundaryField().size(),
|
||||
zeroGradientFvPatchScalarField::typeName
|
||||
);
|
||||
|
||||
forAll(composition.Y(inertSpecie).boundaryField(), patchi)
|
||||
{
|
||||
rhoBoundaryTypes[patchi]
|
||||
= composition.Y(inertSpecie).boundaryField()[patchi].type();
|
||||
}
|
||||
|
||||
|
||||
volScalarField rho
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"rho",
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::NO_READ,
|
||||
IOobject::AUTO_WRITE
|
||||
),
|
||||
thermo.rho(),
|
||||
rhoBoundaryTypes
|
||||
);
|
||||
|
||||
Info<< "\nReading field U\n" << endl;
|
||||
volVectorField U
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"U",
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::MUST_READ,
|
||||
IOobject::AUTO_WRITE
|
||||
),
|
||||
mesh
|
||||
);
|
||||
|
||||
#include "compressibleCreatePhi.H"
|
||||
|
||||
mesh.setFluxRequired(p.name());
|
||||
|
||||
dimensionedScalar rhoMax
|
||||
(
|
||||
dimensionedScalar::lookupOrDefault
|
||||
(
|
||||
"rhoMax",
|
||||
simple.dict(),
|
||||
dimDensity,
|
||||
GREAT
|
||||
)
|
||||
);
|
||||
|
||||
dimensionedScalar rhoMin
|
||||
(
|
||||
dimensionedScalar::lookupOrDefault
|
||||
(
|
||||
"rhoMin",
|
||||
simple.dict(),
|
||||
dimDensity,
|
||||
0
|
||||
)
|
||||
);
|
||||
|
||||
Info<< "Creating turbulence model\n" << endl;
|
||||
autoPtr<compressible::turbulenceModel> turbulence
|
||||
(
|
||||
compressible::turbulenceModel::New
|
||||
(
|
||||
rho,
|
||||
U,
|
||||
phi,
|
||||
thermo
|
||||
)
|
||||
);
|
||||
|
||||
// Set the turbulence into the combustion model
|
||||
combustion->setTurbulence(turbulence());
|
||||
|
||||
Info<< "Creating multi-variate interpolation scheme\n" << endl;
|
||||
multivariateSurfaceInterpolationScheme<scalar>::fieldTable fields;
|
||||
|
||||
forAll(Y, i)
|
||||
{
|
||||
fields.add(Y[i]);
|
||||
}
|
||||
fields.add(thermo.he());
|
||||
|
||||
|
||||
// SLFM-related fields
|
||||
|
||||
Info<<"Creating field mf, mfVar, SDR\n"<<endl;
|
||||
|
||||
volScalarField mf //mixture fraction
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"mf",
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::MUST_READ,
|
||||
IOobject::AUTO_WRITE
|
||||
),
|
||||
mesh
|
||||
);
|
||||
|
||||
volScalarField mfVar //mixture fraction variance
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"mfVar",
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::MUST_READ,
|
||||
IOobject::AUTO_WRITE
|
||||
),
|
||||
mesh
|
||||
);
|
||||
|
||||
volScalarField SDR //scalar dissipation rate
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"SDR",
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::NO_READ,
|
||||
IOobject::AUTO_WRITE
|
||||
),
|
||||
mesh,
|
||||
dimensionedScalar("SDR", dimEnergy/dimTime/dimEnergy, 0.0)
|
||||
);
|
||||
|
||||
volScalarField jlc //eta-index, lower
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"jlc",
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::NO_READ,
|
||||
IOobject::NO_WRITE
|
||||
),
|
||||
mesh,
|
||||
dimensionedScalar("jlc", dimless, 0.0)
|
||||
);
|
||||
|
||||
volScalarField jhc //eta-index, higher
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"jhc",
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::NO_READ,
|
||||
IOobject::NO_WRITE
|
||||
),
|
||||
mesh,
|
||||
dimensionedScalar("jhc", dimless, 0.0)
|
||||
);
|
||||
|
||||
volScalarField jfc //eta-factor
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"jfc",
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::NO_READ,
|
||||
IOobject::NO_WRITE
|
||||
),
|
||||
mesh,
|
||||
dimensionedScalar("jfc", dimless, 0.0)
|
||||
);
|
||||
|
||||
volScalarField vlc //var-index, lower
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"vlc",
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::NO_READ,
|
||||
IOobject::NO_WRITE
|
||||
),
|
||||
mesh,
|
||||
dimensionedScalar("vlc", dimless, 0.0)
|
||||
);
|
||||
|
||||
volScalarField vhc //var-index, higher
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"vhc",
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::NO_READ,
|
||||
IOobject::NO_WRITE
|
||||
),
|
||||
mesh,
|
||||
dimensionedScalar("vhc", dimless, 0.0)
|
||||
);
|
||||
|
||||
volScalarField vfc //var-factor
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"vfc",
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::NO_READ,
|
||||
IOobject::NO_WRITE
|
||||
),
|
||||
mesh,
|
||||
dimensionedScalar("vfc", dimless, 0.0)
|
||||
);
|
||||
|
||||
volScalarField nlc //sdr-index, lower
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"nlc",
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::NO_READ,
|
||||
IOobject::NO_WRITE
|
||||
),
|
||||
mesh,
|
||||
dimensionedScalar("nlc", dimless, 0.0)
|
||||
);
|
||||
|
||||
volScalarField nhc //sdr-index, higher
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"nhc",
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::NO_READ,
|
||||
IOobject::NO_WRITE
|
||||
),
|
||||
mesh,
|
||||
dimensionedScalar("nhc", dimless, 0.0)
|
||||
);
|
||||
|
||||
volScalarField nfc //sdr-factor
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"nfc",
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::NO_READ,
|
||||
IOobject::NO_WRITE
|
||||
),
|
||||
mesh,
|
||||
dimensionedScalar("nfc", dimless, 0.0)
|
||||
);
|
||||
|
||||
//SLFM property read
|
||||
|
||||
IOdictionary SLFMdict //SLFM properties dictionary
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"SLFMdict",
|
||||
runTime.constant(),
|
||||
mesh,
|
||||
IOobject::MUST_READ,
|
||||
IOobject::NO_WRITE
|
||||
)
|
||||
);
|
||||
|
||||
const wordList postSpecieNames(SLFMdict.lookup("postSpecieNames"));
|
||||
labelList postSpecieIndices(postSpecieNames.size(),-1);
|
||||
forAll(postSpecieNames, i)
|
||||
{
|
||||
const word &spName (postSpecieNames[i]);
|
||||
const label spI (composition.species()[spName]);
|
||||
postSpecieIndices[i] = spI;
|
||||
}
|
||||
scalar stoiMF(readScalar(SLFMdict.lookup("stoiMF")));
|
||||
scalar stretchFac(readScalar(SLFMdict.lookup("stretchFactor")));
|
||||
scalar lowerN(readScalar(SLFMdict.lookup("lowerN")));
|
||||
scalar upperN(readScalar(SLFMdict.lookup("upperN")));
|
||||
scalar Sc(readScalar(SLFMdict.lookup("Sc")));
|
||||
const label etamax(lowerN+upperN+1);
|
||||
word outletName(SLFMdict.lookup("outletName"));
|
||||
|
||||
labelList nIntervals (SLFMdict.lookup("detailedN"));
|
||||
|
||||
CMC::BetaGrid bg(SLFMdict);
|
||||
const scalarField &etaGrid(bg.etaSpace());
|
||||
const scalarField amcGridValue(CMC::AMC(etaGrid));
|
||||
|
||||
Info << etaGrid << endl;
|
||||
|
||||
scalarField etaValue(etamax,0);
|
||||
scalarField etaIndex(etamax,0);
|
||||
scalarField deta_h(etamax,0);
|
||||
scalarField deta_l(etamax,0);
|
||||
scalarField deta_m(etamax,0);
|
||||
|
||||
#include "setEtaSpace.H"
|
||||
|
||||
Info<<"\nNumber of eta-point = "<<etamax<<endl;
|
||||
|
||||
|
||||
scalar NVar(readScalar(SLFMdict.lookup("NVar")));
|
||||
scalarField varValue(NVar+1,0);
|
||||
scalarField varIndex(NVar+1,0);
|
||||
|
||||
#include "setVarSpace.H"
|
||||
|
||||
|
||||
|
||||
PtrList<volScalarField> Neta(etamax);
|
||||
|
||||
for(label j=0 ; j<etamax ; j++)
|
||||
{
|
||||
|
||||
|
||||
word Netai = "Neta_"+Foam::name(j);
|
||||
Neta.set
|
||||
(
|
||||
j,
|
||||
new volScalarField
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
Netai,
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::NO_READ,
|
||||
IOobject::NO_WRITE
|
||||
),
|
||||
mesh,
|
||||
dimensionedScalar(Netai, dimless, 0.0)
|
||||
)
|
||||
);
|
||||
}
|
||||
|
||||
#include "createMRF.H"
|
||||
55
applications/solvers/combustion/SLFMFoam/pEqn.H
Normal file
55
applications/solvers/combustion/SLFMFoam/pEqn.H
Normal file
|
|
@ -0,0 +1,55 @@
|
|||
{
|
||||
// Thermodynamic density needs to be updated by psi*d(p) after the
|
||||
// pressure solution - done in 2 parts. Part 1:
|
||||
thermo.rho() -= psi*p;
|
||||
|
||||
volScalarField rAU(1.0/UEqn.A());
|
||||
surfaceScalarField rhorAUf("rhorAUf", fvc::interpolate(rho*rAU));
|
||||
volVectorField HbyA(constrainHbyA(rAU*UEqn.H(), U, p));
|
||||
tUEqn.clear();
|
||||
surfaceScalarField phiHbyA
|
||||
(
|
||||
"phiHbyA",
|
||||
fvc::interpolate(rho)*fvc::flux(HbyA)
|
||||
);
|
||||
|
||||
MRF.makeRelative(fvc::interpolate(rho), phiHbyA);
|
||||
|
||||
// Update the pressure BCs to ensure flux consistency
|
||||
constrainPressure(p, rho, U, phiHbyA, rhorAUf, MRF);
|
||||
|
||||
while (simple.correctNonOrthogonal())
|
||||
{
|
||||
fvScalarMatrix pEqn
|
||||
(
|
||||
fvc::div(phiHbyA)
|
||||
- fvm::laplacian(rhorAUf, p)
|
||||
==
|
||||
fvOptions(psi, p, rho.name())
|
||||
);
|
||||
|
||||
pEqn.solve();
|
||||
|
||||
if (simple.finalNonOrthogonalIter())
|
||||
{
|
||||
phi = phiHbyA + pEqn.flux();
|
||||
}
|
||||
}
|
||||
|
||||
p.relax();
|
||||
|
||||
// Second part of thermodynamic density update
|
||||
thermo.rho() += psi*p;
|
||||
|
||||
#include "compressibleContinuityErrs.H"
|
||||
|
||||
U = HbyA - rAU*fvc::grad(p);
|
||||
U.correctBoundaryConditions();
|
||||
fvOptions.correct(U);
|
||||
|
||||
rho = max(rho, rhoMin);
|
||||
rho = min(rho, rhoMax);
|
||||
rho.relax();
|
||||
|
||||
Info<< "p min/max = " << min(p).value() << ", " << max(p).value() << endl;
|
||||
}
|
||||
78
applications/solvers/combustion/SLFMFoam/preIntegration.H
Normal file
78
applications/solvers/combustion/SLFMFoam/preIntegration.H
Normal file
|
|
@ -0,0 +1,78 @@
|
|||
scalarFieldArray1d C1table //AMC C1 coefficient table
|
||||
(
|
||||
etamax,
|
||||
scalarField(NVar+1, 0)
|
||||
);
|
||||
|
||||
scalarFieldArray2d Ttable
|
||||
(
|
||||
etamax,
|
||||
scalarFieldArray1d
|
||||
(
|
||||
NVar+1,
|
||||
scalarField(NstList.size(), 0)
|
||||
)
|
||||
);
|
||||
|
||||
scalarFieldArray2d Rtable
|
||||
(
|
||||
etamax,
|
||||
scalarFieldArray1d
|
||||
(
|
||||
NVar+1,
|
||||
scalarField(NstList.size(), 0)
|
||||
)
|
||||
);
|
||||
|
||||
scalarFieldArray3d Ytable
|
||||
(
|
||||
Y.size(),
|
||||
scalarFieldArray2d
|
||||
(
|
||||
etamax,
|
||||
scalarFieldArray1d
|
||||
(
|
||||
NVar+1,
|
||||
scalarField(NstList.size(), 0)
|
||||
)
|
||||
)
|
||||
);
|
||||
|
||||
Info<<"Construct C1 coefficient table"<<endl;
|
||||
for(label j = 1; j < etamax - 1; j++)
|
||||
{
|
||||
for(label v = 1; v < varValue.size() - 1; v++)
|
||||
{
|
||||
const scalar mf = etaValue[j];
|
||||
const scalar mfVar = varValue[v]*mf*(1.0-mf);
|
||||
|
||||
CMC::BetaFunction bf(mf, mfVar);
|
||||
CMC::BetaIntegrator bi(bf, bg);
|
||||
|
||||
C1table[j][v] = 1.0/bi.betaIntegrate(amcGridValue);
|
||||
}
|
||||
}
|
||||
|
||||
Info<<"Construct Temp and Rgas table"<<endl;
|
||||
for(label j=0 ; j<etamax ; j++)
|
||||
{
|
||||
for(label v=0 ; v<=NVar ; v++)
|
||||
{
|
||||
const scalar mf = etaValue[j];
|
||||
const scalar mfVar = varValue[v]*mf*(1.0-mf);
|
||||
|
||||
CMC::BetaFunction bf(mf, mfVar);
|
||||
CMC::BetaIntegrator bi(bf, bg);
|
||||
|
||||
for(label n=0 ; n<NstList.size() ; n++)
|
||||
{
|
||||
Ttable[j][v][n] = bi.betaIntegrate(T_SLFM[n]);
|
||||
Rtable[j][v][n] = bi.betaIntegrate(Rgas_SLFM[n]);
|
||||
forAll(Y, yi)
|
||||
{
|
||||
Ytable[yi][j][v][n] = bi.betaIntegrate(Y_SLFM[yi][n]);
|
||||
}
|
||||
}
|
||||
}
|
||||
Info<<scalar(j+1)/etamax*100<<"%"<<endl;
|
||||
}
|
||||
33
applications/solvers/combustion/SLFMFoam/readSLFMlib.H
Normal file
33
applications/solvers/combustion/SLFMFoam/readSLFMlib.H
Normal file
|
|
@ -0,0 +1,33 @@
|
|||
Switch Equilibrium = SLFMdict.lookup("Equilibrium");
|
||||
word NstFolder = SLFMdict.lookup("NstFolder");
|
||||
scalarField NstList = SLFMdict.lookup("NstList");
|
||||
scalarField NstIndex(NstList.size(),0);
|
||||
|
||||
forAll(NstList, n)
|
||||
{
|
||||
NstIndex[n] = n;
|
||||
}
|
||||
|
||||
typedef List<scalarField> scalarFieldArray1d;
|
||||
typedef List<scalarFieldArray1d> scalarFieldArray2d;
|
||||
typedef List<scalarFieldArray2d> scalarFieldArray3d;
|
||||
|
||||
SLFM::FlameStructure slfmLibrary(runTime.constant()/NstFolder, NstList, etaGrid, composition.species());
|
||||
|
||||
scalarFieldArray2d& Y_SLFM (slfmLibrary.Y());
|
||||
scalarFieldArray2d& W_SLFM (slfmLibrary.W());
|
||||
scalarFieldArray1d& T_SLFM (slfmLibrary.T());
|
||||
scalarFieldArray1d& Q_SLFM (slfmLibrary.Q());
|
||||
scalarFieldArray1d& rho_SLFM (slfmLibrary.rho());
|
||||
scalarFieldArray1d& h_SLFM (slfmLibrary.h());
|
||||
scalarFieldArray1d& Rgas_SLFM (slfmLibrary.Rgas());
|
||||
|
||||
for(label n=0 ; n<NstList.size() ; n++)
|
||||
{
|
||||
for(label i=0 ; i<Y.size() ; i++)
|
||||
{
|
||||
Rgas_SLFM[n] += Y_SLFM[i][n]/composition.W(i); // = 1/meanMW
|
||||
}
|
||||
Rgas_SLFM[n] = 8314.46*Rgas_SLFM[n]; //universal gas constant = 8314.46 [J/K.Kmole]
|
||||
|
||||
}
|
||||
51
applications/solvers/combustion/SLFMFoam/setEtaSpace.H
Normal file
51
applications/solvers/combustion/SLFMFoam/setEtaSpace.H
Normal file
|
|
@ -0,0 +1,51 @@
|
|||
//set eta-space
|
||||
scalar upperSum(0), upperDelta(0);
|
||||
scalar lowerSum(0), lowerDelta(0);
|
||||
|
||||
//summation of geometric progression (eta > stoi)
|
||||
upperSum = Foam::pow(stretchFac,upperN)-1.0;
|
||||
upperSum = upperSum/(stretchFac-1.0);
|
||||
|
||||
//summation of geometric progression (eta < stoi)
|
||||
lowerSum = Foam::pow(stretchFac,lowerN)-1.0;
|
||||
lowerSum = lowerSum/(stretchFac-1.0);
|
||||
|
||||
//set reference delta eta
|
||||
upperDelta = (1.0-stoiMF)/upperSum;
|
||||
lowerDelta = stoiMF/lowerSum;
|
||||
|
||||
//set boundary and stoichiometric value of eta
|
||||
etaValue[0] = 0.0;
|
||||
etaValue[lowerN] = stoiMF;
|
||||
etaValue.last() = 1.0;
|
||||
|
||||
//set internal value of eta (eta > stoi)
|
||||
for(label j=0 ; j<upperN-1 ; j++)
|
||||
{
|
||||
etaValue[lowerN+j+1] = etaValue[lowerN+j] + upperDelta*pow(stretchFac,j);
|
||||
}
|
||||
|
||||
//set internal value of eta (eta < stoi)
|
||||
for(label j=0 ; j<lowerN-1 ; j++)
|
||||
{
|
||||
etaValue[lowerN-j-1] = etaValue[lowerN-j] - lowerDelta*pow(stretchFac,j);
|
||||
}
|
||||
|
||||
Info<<etaValue<<endl;
|
||||
|
||||
//calculate delta eta
|
||||
for(label j=1 ; j<etamax-1 ; j++)
|
||||
{
|
||||
deta_h[j] = etaValue[j+1] - etaValue[j];
|
||||
deta_l[j] = etaValue[j] - etaValue[j-1];
|
||||
deta_m[j] = 0.5*(etaValue[j+1] - etaValue[j-1]);
|
||||
}
|
||||
|
||||
//Info<<deta_h<<endl;
|
||||
//Info<<deta_l<<endl;
|
||||
//Info<<deta_m<<endl;
|
||||
|
||||
for(label j=0 ; j<etamax ; j++)
|
||||
{
|
||||
etaIndex[j] = j;
|
||||
}
|
||||
79
applications/solvers/combustion/SLFMFoam/setSDR.H
Normal file
79
applications/solvers/combustion/SLFMFoam/setSDR.H
Normal file
|
|
@ -0,0 +1,79 @@
|
|||
scalarField AMCcoeff(CMC::AMC(etaValue));
|
||||
|
||||
//cell loop
|
||||
|
||||
forAll(mf, cellI)
|
||||
{
|
||||
|
||||
|
||||
scalar jl(0), jh(0), vl(0), vh(0);
|
||||
scalar jfac(0), vfac(0);
|
||||
scalar C1coeff(0);
|
||||
|
||||
//find eta-index, factor
|
||||
if(mf[cellI] < etaValue[1])
|
||||
{
|
||||
jl = 0;
|
||||
jh = 1;
|
||||
}
|
||||
else if(mf[cellI] > etaValue[etamax-2])
|
||||
{
|
||||
jl = etamax-2;
|
||||
jh = etamax-1;
|
||||
}
|
||||
else
|
||||
{
|
||||
jl = label( interpolateXY(mf[cellI], etaValue, etaIndex) );
|
||||
jh = jl+1;
|
||||
}
|
||||
jfac = max(0.0, (mf[cellI]-etaValue[jl])/(etaValue[jh]-etaValue[jl]));
|
||||
|
||||
//find var-index, factor
|
||||
scalar scaledVar = min(0.99999, mfVar[cellI]/(mf[cellI]*(1.0-mf[cellI])+SMALL));
|
||||
if(scaledVar < varValue[1])
|
||||
{
|
||||
vl = 0;
|
||||
vh = 1;
|
||||
}
|
||||
else if(scaledVar > varValue[NVar-1])
|
||||
{
|
||||
vl = NVar-1;
|
||||
vh = NVar;
|
||||
}
|
||||
else
|
||||
{
|
||||
vl = label( interpolateXY(scaledVar, varValue, varIndex) );
|
||||
vh = vl+1;
|
||||
}
|
||||
vfac = max(0.0, (scaledVar-varValue[vl])/(varValue[vh]-varValue[vl]));
|
||||
|
||||
|
||||
//Bi-linear interpolation on j and v
|
||||
//Numerical recipes, 2nd Ed. p.117
|
||||
C1coeff
|
||||
= C1table[jl][vl]*(1-jfac)*(1-vfac)
|
||||
+ C1table[jh][vl]*(jfac)*(1-vfac)
|
||||
+ C1table[jh][vh]*(jfac)*(vfac)
|
||||
+ C1table[jl][vh]*(1-jfac)*(vfac);
|
||||
|
||||
jlc[cellI] = jl;
|
||||
jhc[cellI] = jh;
|
||||
jfc[cellI] = jfac;
|
||||
vlc[cellI] = vl;
|
||||
vhc[cellI] = vh;
|
||||
vfc[cellI] = vfac;
|
||||
|
||||
|
||||
scalar Coeff = AMCcoeff[lowerN];
|
||||
|
||||
if (Equilibrium == true)
|
||||
{
|
||||
Neta[lowerN][cellI] = 0; //Equilibrium
|
||||
}
|
||||
else
|
||||
{
|
||||
Neta[lowerN][cellI] = SDR[cellI]*Coeff*C1coeff; //CSDR at stoichiometic m.f.
|
||||
}
|
||||
}
|
||||
|
||||
|
||||
23
applications/solvers/combustion/SLFMFoam/setVarSpace.H
Normal file
23
applications/solvers/combustion/SLFMFoam/setVarSpace.H
Normal file
|
|
@ -0,0 +1,23 @@
|
|||
//set normalized variance space
|
||||
{
|
||||
scalar sum(0), del(0), fac(1.075);
|
||||
|
||||
sum = Foam::pow(fac,NVar) - 1.0;
|
||||
sum = sum/(fac-1.0);
|
||||
|
||||
del = 1.0/sum;
|
||||
|
||||
varValue[0] = 0.0;
|
||||
varValue[NVar] = 1.0;
|
||||
|
||||
for(label v=0 ; v<NVar-1 ; v++)
|
||||
{
|
||||
varValue[v+1] = varValue[v] + del*Foam::pow(fac,v);
|
||||
}
|
||||
|
||||
for(label v=0 ; v<=NVar ; v++)
|
||||
{
|
||||
varIndex[v] = v;
|
||||
}
|
||||
Info<<varValue<<endl;
|
||||
}
|
||||
60
applications/solvers/combustion/SLFMFoam/updateT_RHO.H
Normal file
60
applications/solvers/combustion/SLFMFoam/updateT_RHO.H
Normal file
|
|
@ -0,0 +1,60 @@
|
|||
//force the saving of the old-time values
|
||||
//important!!
|
||||
|
||||
rho.oldTime();
|
||||
|
||||
forAll(rho, cellI)
|
||||
{
|
||||
scalar& Nst = Neta[lowerN][cellI];
|
||||
|
||||
//find Nst index
|
||||
if(Nst < NstList.first())
|
||||
{
|
||||
nlc[cellI] = 0;
|
||||
nfc[cellI] = 0;
|
||||
}
|
||||
else if(Nst > NstList[NstList.size()-2])
|
||||
{
|
||||
nhc[cellI] = NstList.size()-2; //just before extinction
|
||||
nfc[cellI] = 1;
|
||||
}
|
||||
else
|
||||
{
|
||||
nlc[cellI] = label( interpolateXY(Nst, NstList, NstIndex) );
|
||||
nhc[cellI] = nlc[cellI]+1;
|
||||
nfc[cellI] = (Nst-NstList[nlc[cellI]])
|
||||
/(NstList[nhc[cellI]]-NstList[nlc[cellI]]);
|
||||
}
|
||||
|
||||
//update Favre mean temp. and rho(=1/RT)
|
||||
scalar jl = jlc[cellI], jh = jhc[cellI];
|
||||
scalar vl = vlc[cellI], vh = vhc[cellI];
|
||||
scalar nl = nlc[cellI], nh = nhc[cellI];
|
||||
scalar jfac = jfc[cellI], vfac = vfc[cellI], nfac = nfc[cellI];
|
||||
|
||||
scalar T00 = Ttable[jl][vl][nl]*(1-jfac)+Ttable[jh][vl][nl]*jfac;
|
||||
scalar T01 = Ttable[jl][vl][nh]*(1-jfac)+Ttable[jh][vl][nh]*jfac;
|
||||
scalar T10 = Ttable[jl][vh][nl]*(1-jfac)+Ttable[jh][vh][nl]*jfac;
|
||||
scalar T11 = Ttable[jl][vh][nh]*(1-jfac)+Ttable[jh][vh][nh]*jfac;
|
||||
|
||||
scalar T0 = T00*(1-vfac)+T10*vfac;
|
||||
scalar T1 = T01*(1-vfac)+T11*vfac;
|
||||
|
||||
T[cellI] = T0*(1-nfac)+T1*nfac;
|
||||
|
||||
scalar R00 = Rtable[jl][vl][nl]*(1-jfac)+Rtable[jh][vl][nl]*jfac;
|
||||
scalar R01 = Rtable[jl][vl][nh]*(1-jfac)+Rtable[jh][vl][nh]*jfac;
|
||||
scalar R10 = Rtable[jl][vh][nl]*(1-jfac)+Rtable[jh][vh][nl]*jfac;
|
||||
scalar R11 = Rtable[jl][vh][nh]*(1-jfac)+Rtable[jh][vh][nh]*jfac;
|
||||
|
||||
scalar R0 = R00*(1-vfac)+R10*vfac;
|
||||
scalar R1 = R01*(1-vfac)+R11*vfac;
|
||||
|
||||
scalar Rgas = R0*(1-nfac)+R1*nfac;
|
||||
rho[cellI] = p[cellI]/(Rgas*T[cellI]);
|
||||
}
|
||||
|
||||
//correct processor boundary value of T and rho
|
||||
T.correctBoundaryConditions();
|
||||
rho.correctBoundaryConditions();
|
||||
|
||||
56
applications/solvers/combustion/SLFMFoam/updateYi.H
Normal file
56
applications/solvers/combustion/SLFMFoam/updateYi.H
Normal file
|
|
@ -0,0 +1,56 @@
|
|||
//force the saving of the old-time values
|
||||
//important!!
|
||||
|
||||
forAll(rho, cellI)
|
||||
{
|
||||
scalar& Nst = Neta[lowerN][cellI];
|
||||
|
||||
//find Nst index
|
||||
if(Nst < NstList.first())
|
||||
{
|
||||
nlc[cellI] = 0;
|
||||
nfc[cellI] = 0;
|
||||
}
|
||||
else if(Nst > NstList[NstList.size()-2])
|
||||
{
|
||||
nhc[cellI] = NstList.size()-2; //just before extinction
|
||||
nfc[cellI] = 1;
|
||||
}
|
||||
else
|
||||
{
|
||||
nlc[cellI] = label( interpolateXY(Nst, NstList, NstIndex) );
|
||||
nhc[cellI] = nlc[cellI]+1;
|
||||
nfc[cellI] = (Nst-NstList[nlc[cellI]])
|
||||
/(NstList[nhc[cellI]]-NstList[nlc[cellI]]);
|
||||
}
|
||||
|
||||
//update Favre mean temp. and rho(=1/RT)
|
||||
scalar jl = jlc[cellI], jh = jhc[cellI];
|
||||
scalar vl = vlc[cellI], vh = vhc[cellI];
|
||||
scalar nl = nlc[cellI], nh = nhc[cellI];
|
||||
scalar jfac = jfc[cellI], vfac = vfc[cellI], nfac = nfc[cellI];
|
||||
|
||||
|
||||
forAll(Y, yi)
|
||||
{
|
||||
scalar y00 = Ytable[yi][jl][vl][nl]*(1-jfac)+Ytable[yi][jh][vl][nl]*jfac;
|
||||
scalar y01 = Ytable[yi][jl][vl][nh]*(1-jfac)+Ytable[yi][jh][vl][nh]*jfac;
|
||||
scalar y10 = Ytable[yi][jl][vh][nl]*(1-jfac)+Ytable[yi][jh][vh][nl]*jfac;
|
||||
scalar y11 = Ytable[yi][jl][vh][nh]*(1-jfac)+Ytable[yi][jh][vh][nh]*jfac;
|
||||
|
||||
scalar y0 = y00*(1-vfac)+y10*vfac;
|
||||
scalar y1 = y01*(1-vfac)+y11*vfac;
|
||||
|
||||
scalar ygas = y0*(1-nfac)+y1*nfac;
|
||||
|
||||
Y[yi][cellI] = ygas;
|
||||
}
|
||||
}
|
||||
|
||||
//correct processor boundary value of T and rho
|
||||
forAll(Y, yi)
|
||||
{
|
||||
Y[yi].correctBoundaryConditions();
|
||||
}
|
||||
|
||||
|
||||
|
|
@ -1,4 +0,0 @@
|
|||
SOxFoam.C
|
||||
|
||||
EXE = $(FOAM_USER_APPBIN)/SOxFoam
|
||||
|
||||
|
|
@ -1,11 +0,0 @@
|
|||
EXE_INC = \
|
||||
-I$(LIB_SRC)/finiteVolume/lnInclude \
|
||||
-I${LIB_SRC}/meshTools/lnInclude \
|
||||
-I$(LIB_SRC)/sampling/lnInclude \
|
||||
|
||||
|
||||
EXE_LIBS = \
|
||||
-lfiniteVolume \
|
||||
-lmeshTools \
|
||||
-lfvOptions \
|
||||
-lsampling
|
||||
|
|
@ -1,14 +0,0 @@
|
|||
const scalar Wso2 = 0.064066; // Molecular weight of SO2[kg/mol]
|
||||
const scalar Ws = 0.032065; // Molecular weight of S [kg/mol]
|
||||
|
||||
forAll (mesh.cells(),celli)
|
||||
{
|
||||
scalar sulfurSource = 0.0;
|
||||
|
||||
forAll (fuelSources, si)
|
||||
{
|
||||
sulfurSource += fuelSources[si][celli] * sulfurMassInFuels[si];
|
||||
}
|
||||
|
||||
Sso2[celli] = sulfurSource*Wso2/Ws/mesh.V()[celli];
|
||||
}
|
||||
|
|
@ -1,130 +0,0 @@
|
|||
|
||||
volScalarField rho
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"rho",
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::MUST_READ,
|
||||
IOobject::AUTO_WRITE
|
||||
),
|
||||
mesh
|
||||
);
|
||||
|
||||
volVectorField U
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"U",
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::MUST_READ,
|
||||
IOobject::AUTO_WRITE
|
||||
),
|
||||
mesh
|
||||
);
|
||||
|
||||
surfaceScalarField phi
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"phi",
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::READ_IF_PRESENT,
|
||||
IOobject::AUTO_WRITE
|
||||
),
|
||||
linearInterpolate(rho*U) & mesh.Sf()
|
||||
);
|
||||
|
||||
volScalarField nut
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"nut",
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::MUST_READ,
|
||||
IOobject::AUTO_WRITE
|
||||
),
|
||||
mesh
|
||||
);
|
||||
|
||||
volScalarField mut
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"mut",
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::NO_READ,
|
||||
IOobject::AUTO_WRITE
|
||||
),
|
||||
rho * nut
|
||||
);
|
||||
|
||||
volScalarField SO2
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"SO2",
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::MUST_READ,
|
||||
IOobject::AUTO_WRITE
|
||||
),
|
||||
mesh
|
||||
);
|
||||
|
||||
volScalarField Sso2
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"Sso2",
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::MUST_READ,
|
||||
IOobject::AUTO_WRITE
|
||||
),
|
||||
mesh,
|
||||
dimensionedScalar("Sso2",dimDensity/dimTime,0.0)
|
||||
);
|
||||
|
||||
IOdictionary SOxProperties
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"SOxProperties",
|
||||
runTime.constant(),
|
||||
mesh,
|
||||
IOobject::MUST_READ_IF_MODIFIED,
|
||||
IOobject::NO_WRITE
|
||||
)
|
||||
);
|
||||
|
||||
wordList fuelSourceNames(SOxProperties.lookup("fuelSourceNames"));
|
||||
|
||||
PtrList<volScalarField::Internal> fuelSources(fuelSourceNames.size());
|
||||
|
||||
forAll (fuelSources, si)
|
||||
{
|
||||
fuelSources.set
|
||||
(
|
||||
si,
|
||||
new volScalarField::Internal
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
fuelSourceNames[si],
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::MUST_READ,
|
||||
IOobject::AUTO_WRITE
|
||||
),
|
||||
mesh
|
||||
)
|
||||
);
|
||||
}
|
||||
|
||||
scalarField sulfurMassInFuels(SOxProperties.lookup("sulfurMassInFuels"));
|
||||
40
applications/solvers/lagrangian/ThermoMPPICFoam/EEqn.H
Normal file
40
applications/solvers/lagrangian/ThermoMPPICFoam/EEqn.H
Normal file
|
|
@ -0,0 +1,40 @@
|
|||
{
|
||||
volScalarField& he = thermo.he();
|
||||
|
||||
fvScalarMatrix EEqn
|
||||
(
|
||||
fvm::ddt(alphac,rho, he) + mvConvection->fvmDiv(alphacPhi, he)
|
||||
+ fvc::ddt(alphac,rho, K) + fvc::div(alphacPhi, K)
|
||||
+ (
|
||||
he.name() == "e"
|
||||
? fvc::div
|
||||
(
|
||||
fvc::absolute(phi/fvc::interpolate(rho), U),
|
||||
p,
|
||||
"div(phiv,p)"
|
||||
)
|
||||
: -dpdt
|
||||
)
|
||||
- fvm::laplacian(alphac*turbulence->alphaEff(), he)
|
||||
==
|
||||
alphac*rho*(U&g)
|
||||
+ parcels.Sh(he)
|
||||
+ radiation->Sh(thermo)
|
||||
+ combustion->Sh()
|
||||
+ fvOptions(alphac,rho, he)
|
||||
);
|
||||
|
||||
EEqn.relax();
|
||||
|
||||
fvOptions.constrain(EEqn);
|
||||
|
||||
EEqn.solve();
|
||||
|
||||
fvOptions.correct(he);
|
||||
|
||||
thermo.correct();
|
||||
radiation->correct();
|
||||
|
||||
Info<< "T gas min/max = " << min(T).value() << ", "
|
||||
<< max(T).value() << endl;
|
||||
}
|
||||
|
|
@ -0,0 +1,3 @@
|
|||
ThermoMPPICFoam.C
|
||||
|
||||
EXE = $(FOAM_APPBIN)/ThermoMPPICFoam
|
||||
56
applications/solvers/lagrangian/ThermoMPPICFoam/Make/options
Normal file
56
applications/solvers/lagrangian/ThermoMPPICFoam/Make/options
Normal file
|
|
@ -0,0 +1,56 @@
|
|||
EXE_INC = \
|
||||
-I. \
|
||||
-I$(LIB_SRC)/finiteVolume/lnInclude \
|
||||
-I${LIB_SRC}/meshTools/lnInclude \
|
||||
-I$(LIB_SRC)/TurbulenceModels/turbulenceModels/lnInclude \
|
||||
-I$(LIB_SRC)/TurbulenceModels/compressible/lnInclude \
|
||||
-I$(LIB_SRC)/lagrangian/basic/lnInclude \
|
||||
-I$(LIB_SRC)/lagrangian/intermediate/lnInclude \
|
||||
-I$(LIB_SRC)/lagrangian/coalCombustion/lnInclude \
|
||||
-I$(LIB_SRC)/lagrangian/distributionModels/lnInclude \
|
||||
-I$(LIB_SRC)/thermophysicalModels/specie/lnInclude \
|
||||
-I$(LIB_SRC)/transportModels/compressible/lnInclude \
|
||||
-I$(LIB_SRC)/thermophysicalModels/basic/lnInclude \
|
||||
-I$(LIB_SRC)/thermophysicalModels/properties/liquidProperties/lnInclude \
|
||||
-I$(LIB_SRC)/thermophysicalModels/properties/liquidMixtureProperties/lnInclude \
|
||||
-I$(LIB_SRC)/thermophysicalModels/properties/solidProperties/lnInclude \
|
||||
-I$(LIB_SRC)/thermophysicalModels/properties/solidMixtureProperties/lnInclude \
|
||||
-I$(LIB_SRC)/thermophysicalModels/thermophysicalFunctions/lnInclude \
|
||||
-I$(LIB_SRC)/thermophysicalModels/reactionThermo/lnInclude \
|
||||
-I$(LIB_SRC)/thermophysicalModels/SLGThermo/lnInclude \
|
||||
-I$(LIB_SRC)/thermophysicalModels/chemistryModel/lnInclude \
|
||||
-I$(LIB_SRC)/thermophysicalModels/radiation/lnInclude \
|
||||
-I$(LIB_SRC)/ODE/lnInclude \
|
||||
-I$(LIB_SRC)/regionModels/regionModel/lnInclude \
|
||||
-I$(LIB_SRC)/regionModels/surfaceFilmModels/lnInclude \
|
||||
-I$(LIB_SRC)/combustionModels/lnInclude \
|
||||
-I$(LIB_SRC)/sampling/lnInclude \
|
||||
-I$(FOAM_SOLVERS)/combustion/reactingFoam
|
||||
|
||||
|
||||
EXE_LIBS = \
|
||||
-lfiniteVolume \
|
||||
-lmeshTools \
|
||||
-lturbulenceModels \
|
||||
-lcompressibleTurbulenceModels \
|
||||
-llagrangian \
|
||||
-llagrangianIntermediate \
|
||||
-llagrangianTurbulence \
|
||||
-lspecie \
|
||||
-lcompressibleTransportModels \
|
||||
-lfluidThermophysicalModels \
|
||||
-lliquidProperties \
|
||||
-lliquidMixtureProperties \
|
||||
-lsolidProperties \
|
||||
-lsolidMixtureProperties \
|
||||
-lthermophysicalFunctions \
|
||||
-lreactionThermophysicalModels \
|
||||
-lSLGThermo \
|
||||
-lchemistryModel \
|
||||
-lradiationModels \
|
||||
-lODE \
|
||||
-lregionModels \
|
||||
-lsurfaceFilmModels \
|
||||
-lcombustionModels \
|
||||
-lfvOptions \
|
||||
-lsampling
|
||||
|
|
@ -0,0 +1,147 @@
|
|||
/*---------------------------------------------------------------------------*\
|
||||
========= |
|
||||
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
||||
\\ / O peration |
|
||||
\\ / A nd | Copyright (C) 2011-2016 OpenFOAM Foundation
|
||||
\\/ M anipulation |
|
||||
-------------------------------------------------------------------------------
|
||||
License
|
||||
This file is part of OpenFOAM.
|
||||
|
||||
OpenFOAM is free software: you can redistribute it and/or modify it
|
||||
under the terms of the GNU General Public License as published by
|
||||
the Free Software Foundation, either version 3 of the License, or
|
||||
(at your option) any later version.
|
||||
|
||||
OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
|
||||
ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
|
||||
FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
|
||||
for more details.
|
||||
|
||||
You should have received a copy of the GNU General Public License
|
||||
along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
|
||||
|
||||
Application
|
||||
ThermoMPPICFoam
|
||||
|
||||
Description
|
||||
Transient solver for compressible, turbulent flow with a reacting,
|
||||
multiphase particle cloud, and optional sources/constraints.
|
||||
|
||||
\*---------------------------------------------------------------------------*/
|
||||
|
||||
#include "fvCFD.H"
|
||||
#include "turbulentFluidThermoModel.H"
|
||||
#include "basicThermoKinematicMPPICCloud.H"
|
||||
#include "rhoCombustionModel.H"
|
||||
#include "radiationModel.H"
|
||||
#include "fvOptions.H"
|
||||
#include "SLGThermo.H"
|
||||
#include "pimpleControl.H"
|
||||
#include "localEulerDdtScheme.H"
|
||||
#include "fvcSmooth.H"
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
int main(int argc, char *argv[])
|
||||
{
|
||||
#include "postProcess.H"
|
||||
|
||||
#include "setRootCase.H"
|
||||
#include "createTime.H"
|
||||
#include "createMesh.H"
|
||||
#include "createControl.H"
|
||||
#include "createTimeControls.H"
|
||||
#include "createRDeltaT.H"
|
||||
#include "createFields.H"
|
||||
#include "createFieldRefs.H"
|
||||
#include "createFvOptions.H"
|
||||
#include "initContinuityErrs.H"
|
||||
|
||||
turbulence->validate();
|
||||
|
||||
#include "compressibleCourantNo.H"
|
||||
#include "setInitialDeltaT.H"
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
Info<< "\nStarting time loop\n" << endl;
|
||||
|
||||
while (runTime.run())
|
||||
{
|
||||
#include "readTimeControls.H"
|
||||
#include "compressibleCourantNo.H"
|
||||
#include "setDeltaT.H"
|
||||
|
||||
runTime++;
|
||||
|
||||
Info<< "Time = " << runTime.timeName() << nl << endl;
|
||||
|
||||
parcels.evolve();
|
||||
|
||||
// Update continuous phase volume fraction field
|
||||
alphac = max(1.0 - parcels.theta(), alphacMin);
|
||||
alphac.correctBoundaryConditions();
|
||||
alphacf = fvc::interpolate(alphac);
|
||||
alphacPhi = alphacf*phi;
|
||||
|
||||
fvVectorMatrix cloudSU(parcels.SU(U));
|
||||
volVectorField cloudVolSUSu
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"cloudVolSUSu",
|
||||
runTime.timeName(),
|
||||
mesh
|
||||
),
|
||||
mesh,
|
||||
dimensionedVector
|
||||
(
|
||||
"0",
|
||||
cloudSU.dimensions()/dimVolume,
|
||||
vector::zero
|
||||
),
|
||||
zeroGradientFvPatchVectorField::typeName
|
||||
);
|
||||
|
||||
cloudVolSUSu.primitiveFieldRef() = -cloudSU.source()/mesh.V();
|
||||
cloudVolSUSu.correctBoundaryConditions();
|
||||
cloudSU.source() = vector::zero;
|
||||
|
||||
#include "rhoEqn.H"
|
||||
|
||||
// --- Pressure-velocity PIMPLE corrector loop
|
||||
while (pimple.loop())
|
||||
{
|
||||
#include "UEqn.H"
|
||||
#include "YEqn.H"
|
||||
#include "EEqn.H"
|
||||
|
||||
// --- Pressure corrector loop
|
||||
while (pimple.correct())
|
||||
{
|
||||
#include "pEqn.H"
|
||||
}
|
||||
|
||||
if (pimple.turbCorr())
|
||||
{
|
||||
turbulence->correct();
|
||||
}
|
||||
}
|
||||
|
||||
rho = thermo.rho();
|
||||
|
||||
runTime.write();
|
||||
|
||||
Info<< "ExecutionTime = " << runTime.elapsedCpuTime() << " s"
|
||||
<< " ClockTime = " << runTime.elapsedClockTime() << " s"
|
||||
<< nl << endl;
|
||||
}
|
||||
|
||||
Info<< "End\n" << endl;
|
||||
|
||||
return 0;
|
||||
}
|
||||
|
||||
|
||||
// ************************************************************************* //
|
||||
43
applications/solvers/lagrangian/ThermoMPPICFoam/UEqn.H
Normal file
43
applications/solvers/lagrangian/ThermoMPPICFoam/UEqn.H
Normal file
|
|
@ -0,0 +1,43 @@
|
|||
fvVectorMatrix UEqn
|
||||
(
|
||||
fvm::ddt(alphac, rho, U)
|
||||
+ fvm::div(alphacPhi, U)
|
||||
- fvm::laplacian(alphac*turbulence->muEff(),U)
|
||||
- fvc::div(alphac*turbulence->muEff()*dev2(fvc::grad(U)().T()))
|
||||
- fvm::Sp(fvc::ddt(alphac,rho) + fvc::div(alphacPhi), U)
|
||||
==
|
||||
cloudSU
|
||||
+ fvOptions(alphac,rho, U)
|
||||
);
|
||||
|
||||
UEqn.relax();
|
||||
|
||||
volScalarField rAU(1.0/UEqn.A());
|
||||
surfaceScalarField rAUf("rAUf", fvc::interpolate(rAU));
|
||||
|
||||
surfaceScalarField rhorAUf("rhorAUf", fvc::interpolate(rho*rAU));
|
||||
|
||||
surfaceScalarField phicForces
|
||||
(
|
||||
(fvc::interpolate(rAU*cloudVolSUSu) & mesh.Sf())
|
||||
+ alphacf*rhorAUf*(g & mesh.Sf())
|
||||
);
|
||||
|
||||
fvOptions.constrain(UEqn);
|
||||
|
||||
if (pimple.momentumPredictor())
|
||||
{
|
||||
solve
|
||||
(
|
||||
UEqn
|
||||
==
|
||||
fvc::reconstruct
|
||||
(
|
||||
phicForces/rAUf
|
||||
)
|
||||
- fvc::grad(p)
|
||||
);
|
||||
|
||||
fvOptions.correct(U);
|
||||
K = 0.5*magSqr(U);
|
||||
}
|
||||
53
applications/solvers/lagrangian/ThermoMPPICFoam/YEqn.H
Normal file
53
applications/solvers/lagrangian/ThermoMPPICFoam/YEqn.H
Normal file
|
|
@ -0,0 +1,53 @@
|
|||
tmp<fv::convectionScheme<scalar>> mvConvection
|
||||
(
|
||||
fv::convectionScheme<scalar>::New
|
||||
(
|
||||
mesh,
|
||||
fields,
|
||||
phi,
|
||||
mesh.divScheme("div(phi,Yi_h)")
|
||||
)
|
||||
);
|
||||
|
||||
{
|
||||
combustion->correct();
|
||||
dQ = combustion->dQ();
|
||||
label inertIndex = -1;
|
||||
volScalarField Yt(0.0*Y[0]);
|
||||
|
||||
forAll(Y, i)
|
||||
{
|
||||
if (Y[i].name() != inertSpecie)
|
||||
{
|
||||
volScalarField& Yi = Y[i];
|
||||
|
||||
fvScalarMatrix YEqn
|
||||
(
|
||||
fvm::ddt(alphac,rho, Yi)
|
||||
+ mvConvection->fvmDiv(alphacPhi, Yi)
|
||||
- fvm::laplacian(alphac*turbulence->muEff(), Yi)
|
||||
==
|
||||
combustion->R(Yi)
|
||||
+ fvOptions(alphac,rho, Yi)
|
||||
);
|
||||
|
||||
YEqn.relax();
|
||||
|
||||
fvOptions.constrain(YEqn);
|
||||
|
||||
YEqn.solve(mesh.solver("Yi"));
|
||||
|
||||
fvOptions.correct(Yi);
|
||||
|
||||
Yi.max(0.0);
|
||||
Yt += Yi;
|
||||
}
|
||||
else
|
||||
{
|
||||
inertIndex = i;
|
||||
}
|
||||
}
|
||||
|
||||
Y[inertIndex] = scalar(1) - Yt;
|
||||
Y[inertIndex].max(0.0);
|
||||
}
|
||||
|
|
@ -0,0 +1,26 @@
|
|||
Info<< "\nConstructing reacting cloud" << endl;
|
||||
basicThermoKinematicMPPICCloud parcels
|
||||
(
|
||||
"reactingCloud",
|
||||
rho,
|
||||
U,
|
||||
slgThermo.thermo().mu(),
|
||||
g
|
||||
);
|
||||
|
||||
|
||||
// Continuouse phase volume fraction lower limit
|
||||
scalar alphacMin
|
||||
(
|
||||
1.0
|
||||
- readScalar
|
||||
(
|
||||
parcels.particleProperties().subDict("constantProperties").lookup("alphaMax")
|
||||
)
|
||||
);
|
||||
|
||||
// Update alphac from the particle locations
|
||||
alphac = max(1.0 - parcels.theta(), alphacMin);
|
||||
alphac.correctBoundaryConditions();
|
||||
alphacf = fvc::interpolate(alphac);
|
||||
alphacPhi = alphacf*phi;
|
||||
|
|
@ -0,0 +1,2 @@
|
|||
const volScalarField& T = thermo.T();
|
||||
const volScalarField& psi = thermo.psi();
|
||||
158
applications/solvers/lagrangian/ThermoMPPICFoam/createFields.H
Normal file
158
applications/solvers/lagrangian/ThermoMPPICFoam/createFields.H
Normal file
|
|
@ -0,0 +1,158 @@
|
|||
#include "readGravitationalAcceleration.H"
|
||||
|
||||
Info<< "Creating combustion model\n" << endl;
|
||||
|
||||
autoPtr<combustionModels::rhoCombustionModel> combustion
|
||||
(
|
||||
combustionModels::rhoCombustionModel::New(mesh)
|
||||
);
|
||||
|
||||
rhoReactionThermo& thermo = combustion->thermo();
|
||||
thermo.validate(args.executable(), "h", "e");
|
||||
|
||||
SLGThermo slgThermo(mesh, thermo);
|
||||
|
||||
basicSpecieMixture& composition = thermo.composition();
|
||||
PtrList<volScalarField>& Y = composition.Y();
|
||||
|
||||
const word inertSpecie(thermo.lookup("inertSpecie"));
|
||||
|
||||
if (!composition.contains(inertSpecie))
|
||||
{
|
||||
FatalErrorInFunction
|
||||
<< "Specified inert specie '" << inertSpecie << "' not found in "
|
||||
<< "species list. Available species:" << composition.species()
|
||||
<< exit(FatalError);
|
||||
}
|
||||
|
||||
volScalarField& p = thermo.p();
|
||||
|
||||
volScalarField rho
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"rho",
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::NO_READ,
|
||||
IOobject::AUTO_WRITE
|
||||
),
|
||||
thermo.rho()
|
||||
);
|
||||
|
||||
Info<< "\nReading field U\n" << endl;
|
||||
volVectorField U
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"U",
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::MUST_READ,
|
||||
IOobject::AUTO_WRITE
|
||||
),
|
||||
mesh
|
||||
);
|
||||
|
||||
#include "compressibleCreatePhi.H"
|
||||
|
||||
mesh.setFluxRequired(p.name());
|
||||
|
||||
dimensionedScalar rhoMax
|
||||
(
|
||||
dimensionedScalar::lookupOrDefault
|
||||
(
|
||||
"rhoMax",
|
||||
pimple.dict(),
|
||||
dimDensity,
|
||||
GREAT
|
||||
)
|
||||
);
|
||||
|
||||
dimensionedScalar rhoMin
|
||||
(
|
||||
dimensionedScalar::lookupOrDefault
|
||||
(
|
||||
"rhoMin",
|
||||
pimple.dict(),
|
||||
dimDensity,
|
||||
0
|
||||
)
|
||||
);
|
||||
|
||||
Info<< "Creating turbulence model\n" << endl;
|
||||
autoPtr<compressible::turbulenceModel> turbulence
|
||||
(
|
||||
compressible::turbulenceModel::New
|
||||
(
|
||||
rho,
|
||||
U,
|
||||
phi,
|
||||
thermo
|
||||
)
|
||||
);
|
||||
|
||||
// Set the turbulence into the combustion model
|
||||
combustion->setTurbulence(turbulence());
|
||||
|
||||
Info<< "Creating field dpdt\n" << endl;
|
||||
volScalarField dpdt
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"dpdt",
|
||||
runTime.timeName(),
|
||||
mesh
|
||||
),
|
||||
mesh,
|
||||
dimensionedScalar("dpdt", p.dimensions()/dimTime, 0)
|
||||
);
|
||||
|
||||
Info<< "Creating field kinetic energy K\n" << endl;
|
||||
volScalarField K("K", 0.5*magSqr(U));
|
||||
|
||||
Info<< "Creating multi-variate interpolation scheme\n" << endl;
|
||||
multivariateSurfaceInterpolationScheme<scalar>::fieldTable fields;
|
||||
|
||||
forAll(Y, i)
|
||||
{
|
||||
fields.add(Y[i]);
|
||||
}
|
||||
fields.add(thermo.he());
|
||||
|
||||
volScalarField dQ
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"dQ",
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::NO_READ,
|
||||
IOobject::AUTO_WRITE
|
||||
),
|
||||
mesh,
|
||||
dimensionedScalar("dQ", dimEnergy/dimTime, 0.0)
|
||||
);
|
||||
|
||||
volScalarField alphac
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"alphac",
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::READ_IF_PRESENT,
|
||||
IOobject::AUTO_WRITE
|
||||
),
|
||||
mesh,
|
||||
dimensionedScalar("0", dimless, 1)
|
||||
);
|
||||
|
||||
surfaceScalarField alphacf("alphacf", fvc::interpolate(alphac));
|
||||
surfaceScalarField alphacPhi("alphacPhi", alphacf*phi);
|
||||
|
||||
volScalarField alphacrho("alphacrho", alphac*rho);
|
||||
|
||||
#include "createMRF.H"
|
||||
#include "createRadiationModel.H"
|
||||
#include "createClouds.H"
|
||||
81
applications/solvers/lagrangian/ThermoMPPICFoam/pEqn.H
Normal file
81
applications/solvers/lagrangian/ThermoMPPICFoam/pEqn.H
Normal file
|
|
@ -0,0 +1,81 @@
|
|||
{
|
||||
rho = thermo.rho();
|
||||
|
||||
// Thermodynamic density needs to be updated by psi*d(p) after the
|
||||
// pressure solution - done in 2 parts. Part 1:
|
||||
thermo.rho() -= psi*p;
|
||||
|
||||
rhorAUf= fvc::interpolate(rho*rAU);
|
||||
surfaceScalarField rhof("rhof", fvc::interpolate(rho));
|
||||
|
||||
volVectorField HbyA("HbyA", U);
|
||||
HbyA = rAU*UEqn.H();
|
||||
|
||||
surfaceScalarField phiHbyA
|
||||
(
|
||||
"phiHbyA",
|
||||
(
|
||||
(fvc::interpolate(rho*HbyA) & mesh.Sf())
|
||||
+ rhorAUf*fvc::ddtCorr(rho, U, phi)
|
||||
+ phicForces*rhof
|
||||
|
||||
)
|
||||
);
|
||||
|
||||
MRF.makeRelative(fvc::interpolate(rho), phiHbyA);
|
||||
|
||||
// Update the pressure BCs to ensure flux consistency
|
||||
constrainPressure(p, rho, U, phiHbyA, rhorAUf, MRF);
|
||||
|
||||
fvScalarMatrix pDDtEqn
|
||||
(
|
||||
fvc::ddt(alphac,rho) + alphac*psi*correction(fvm::ddt(p))
|
||||
+ fvc::div(alphacf*phiHbyA)
|
||||
==
|
||||
fvOptions(alphac*psi, p, rho.name())
|
||||
);
|
||||
|
||||
while (pimple.correctNonOrthogonal())
|
||||
{
|
||||
fvScalarMatrix pEqn
|
||||
(
|
||||
pDDtEqn
|
||||
- fvm::laplacian(alphacf*rhorAUf, p)
|
||||
);
|
||||
|
||||
pEqn.solve(mesh.solver(p.select(pimple.finalInnerIter())));
|
||||
|
||||
if (pimple.finalNonOrthogonalIter())
|
||||
{
|
||||
|
||||
phi = phiHbyA + pEqn.flux()/alphacf;
|
||||
|
||||
alphacPhi = alphacf*phi;
|
||||
|
||||
p.relax();
|
||||
|
||||
thermo.rho() += psi*p;
|
||||
|
||||
#include "rhoEqn.H"
|
||||
#include "compressibleContinuityErrs.H"
|
||||
|
||||
U = HbyA + rAU*fvc::reconstruct(phicForces/rAUf) - rAU*fvc::grad(p);
|
||||
|
||||
}
|
||||
}
|
||||
|
||||
U.correctBoundaryConditions();
|
||||
fvOptions.correct(U);
|
||||
K = 0.5*magSqr(U);
|
||||
|
||||
if (thermo.dpdt())
|
||||
{
|
||||
dpdt = fvc::ddt(p);
|
||||
}
|
||||
|
||||
rho = thermo.rho();
|
||||
rho = max(rho, rhoMin);
|
||||
rho = min(rho, rhoMax);
|
||||
|
||||
Info<< "p min/max = " << min(p).value() << ", " << max(p).value() << endl;
|
||||
}
|
||||
52
applications/solvers/lagrangian/ThermoMPPICFoam/rhoEqn.H
Normal file
52
applications/solvers/lagrangian/ThermoMPPICFoam/rhoEqn.H
Normal file
|
|
@ -0,0 +1,52 @@
|
|||
/*---------------------------------------------------------------------------*\
|
||||
========= |
|
||||
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
||||
\\ / O peration |
|
||||
\\ / A nd | Copyright (C) 2011-2015 OpenFOAM Foundation
|
||||
\\/ M anipulation |
|
||||
-------------------------------------------------------------------------------
|
||||
License
|
||||
This file is part of OpenFOAM.
|
||||
|
||||
OpenFOAM is free software: you can redistribute it and/or modify it
|
||||
under the terms of the GNU General Public License as published by
|
||||
the Free Software Foundation, either version 3 of the License, or
|
||||
(at your option) any later version.
|
||||
|
||||
OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
|
||||
ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
|
||||
FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
|
||||
for more details.
|
||||
|
||||
You should have received a copy of the GNU General Public License
|
||||
along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
|
||||
|
||||
Global
|
||||
rhoEqn
|
||||
|
||||
Description
|
||||
Solve the continuity for density.
|
||||
|
||||
\*---------------------------------------------------------------------------*/
|
||||
|
||||
{
|
||||
fvScalarMatrix rhoEqn
|
||||
(
|
||||
fvm::ddt(alphac,rho)
|
||||
+ fvc::div(alphacPhi)
|
||||
==
|
||||
fvOptions(alphac,rho)
|
||||
);
|
||||
|
||||
rhoEqn.solve();
|
||||
|
||||
fvOptions.correct(rho);
|
||||
|
||||
Info<< "rho min/max = " << min(rho).value() << ", " << max(rho).value()
|
||||
<< endl;
|
||||
|
||||
rho = max(rho, rhoMin);
|
||||
rho = min(rho, rhoMax);
|
||||
}
|
||||
|
||||
// ************************************************************************* //
|
||||
134
applications/solvers/lagrangian/ThermoMPPICFoam/setRDeltaT.H
Normal file
134
applications/solvers/lagrangian/ThermoMPPICFoam/setRDeltaT.H
Normal file
|
|
@ -0,0 +1,134 @@
|
|||
/*---------------------------------------------------------------------------*\
|
||||
========= |
|
||||
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
||||
\\ / O peration |
|
||||
\\ / A nd | Copyright (C) 2011-2016 OpenFOAM Foundation
|
||||
\\/ M anipulation |
|
||||
-------------------------------------------------------------------------------
|
||||
License
|
||||
This file is part of OpenFOAM.
|
||||
|
||||
OpenFOAM is free software: you can redistribute it and/or modify it
|
||||
under the terms of the GNU General Public License as published by
|
||||
the Free Software Foundation, either version 3 of the License, or
|
||||
(at your option) any later version.
|
||||
|
||||
OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
|
||||
ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
|
||||
FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
|
||||
for more details.
|
||||
|
||||
You should have received a copy of the GNU General Public License
|
||||
along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
|
||||
|
||||
\*---------------------------------------------------------------------------*/
|
||||
|
||||
{
|
||||
volScalarField& rDeltaT = trDeltaT.ref();
|
||||
|
||||
const dictionary& pimpleDict = pimple.dict();
|
||||
|
||||
// Maximum flow Courant number
|
||||
scalar maxCo(readScalar(pimpleDict.lookup("maxCo")));
|
||||
|
||||
// Maximum time scale
|
||||
scalar maxDeltaT(pimpleDict.lookupOrDefault<scalar>("maxDeltaT", GREAT));
|
||||
|
||||
// Smoothing parameter (0-1) when smoothing iterations > 0
|
||||
scalar rDeltaTSmoothingCoeff
|
||||
(
|
||||
pimpleDict.lookupOrDefault<scalar>("rDeltaTSmoothingCoeff", 0.1)
|
||||
);
|
||||
|
||||
// Damping coefficient (1-0)
|
||||
scalar rDeltaTDampingCoeff
|
||||
(
|
||||
pimpleDict.lookupOrDefault<scalar>("rDeltaTDampingCoeff", 0.2)
|
||||
);
|
||||
|
||||
// Maximum change in cell temperature per iteration
|
||||
// (relative to previous value)
|
||||
scalar alphaTemp(pimpleDict.lookupOrDefault("alphaTemp", 0.05));
|
||||
|
||||
|
||||
Info<< "Time scales min/max:" << endl;
|
||||
|
||||
// Cache old reciprocal time scale field
|
||||
volScalarField rDeltaT0("rDeltaT0", rDeltaT);
|
||||
|
||||
// Flow time scale
|
||||
{
|
||||
rDeltaT.ref() =
|
||||
(
|
||||
fvc::surfaceSum(mag(phi))()()
|
||||
/((2*maxCo)*mesh.V()*rho())
|
||||
);
|
||||
|
||||
// Limit the largest time scale
|
||||
rDeltaT.max(1/maxDeltaT);
|
||||
|
||||
Info<< " Flow = "
|
||||
<< gMin(1/rDeltaT.primitiveField()) << ", "
|
||||
<< gMax(1/rDeltaT.primitiveField()) << endl;
|
||||
}
|
||||
|
||||
// Reaction source time scale
|
||||
{
|
||||
volScalarField::Internal rDeltaTT
|
||||
(
|
||||
mag
|
||||
(
|
||||
parcels.hsTrans()/(mesh.V()*runTime.deltaT())
|
||||
+ combustion->Sh()()
|
||||
)
|
||||
/(
|
||||
alphaTemp
|
||||
*rho()
|
||||
*thermo.Cp()()()
|
||||
*T()
|
||||
)
|
||||
);
|
||||
|
||||
Info<< " Temperature = "
|
||||
<< gMin(1/(rDeltaTT.field() + VSMALL)) << ", "
|
||||
<< gMax(1/(rDeltaTT.field() + VSMALL)) << endl;
|
||||
|
||||
rDeltaT.ref() = max
|
||||
(
|
||||
rDeltaT(),
|
||||
rDeltaTT
|
||||
);
|
||||
}
|
||||
|
||||
// Update tho boundary values of the reciprocal time-step
|
||||
rDeltaT.correctBoundaryConditions();
|
||||
|
||||
// Spatially smooth the time scale field
|
||||
if (rDeltaTSmoothingCoeff < 1.0)
|
||||
{
|
||||
fvc::smooth(rDeltaT, rDeltaTSmoothingCoeff);
|
||||
}
|
||||
|
||||
// Limit rate of change of time scale
|
||||
// - reduce as much as required
|
||||
// - only increase at a fraction of old time scale
|
||||
if
|
||||
(
|
||||
rDeltaTDampingCoeff < 1.0
|
||||
&& runTime.timeIndex() > runTime.startTimeIndex() + 1
|
||||
)
|
||||
{
|
||||
rDeltaT = max
|
||||
(
|
||||
rDeltaT,
|
||||
(scalar(1.0) - rDeltaTDampingCoeff)*rDeltaT0
|
||||
);
|
||||
}
|
||||
|
||||
Info<< " Overall = "
|
||||
<< gMin(1/rDeltaT.primitiveField())
|
||||
<< ", " << gMax(1/rDeltaT.primitiveField()) << endl;
|
||||
}
|
||||
|
||||
|
||||
// ************************************************************************* //
|
||||
|
|
@ -0,0 +1,32 @@
|
|||
{
|
||||
volScalarField& he = thermo.he();
|
||||
|
||||
fvScalarMatrix EEqn
|
||||
(
|
||||
mvConvection->fvmDiv(phi, he)
|
||||
+ (
|
||||
he.name() == "e"
|
||||
? fvc::div(phi, volScalarField("Ekp", 0.5*magSqr(U) + p/rho))
|
||||
: fvc::div(phi, volScalarField("K", 0.5*magSqr(U)))
|
||||
)
|
||||
- fvm::laplacian(turbulence->alphaEff(), he)
|
||||
==
|
||||
parcels.Sh(he)
|
||||
+ radiation->Sh(thermo)
|
||||
+ combustion->Sh()
|
||||
+ fvOptions(rho, he)
|
||||
);
|
||||
|
||||
EEqn.relax();
|
||||
|
||||
fvOptions.constrain(EEqn);
|
||||
|
||||
EEqn.solve();
|
||||
|
||||
fvOptions.correct(he);
|
||||
thermo.correct();
|
||||
radiation->correct();
|
||||
|
||||
Info<< "T gas min/max = " << min(T).value() << ", "
|
||||
<< max(T).value() << endl;
|
||||
}
|
||||
|
|
@ -0,0 +1,3 @@
|
|||
simpleCoalcombustionFoam.C
|
||||
|
||||
EXE = $(FOAM_APPBIN)/simpleCoalcombustionFoam
|
||||
|
|
@ -0,0 +1,57 @@
|
|||
EXE_INC = \
|
||||
-I. \
|
||||
-I$(LIB_SRC)/finiteVolume/lnInclude \
|
||||
-I${LIB_SRC}/meshTools/lnInclude \
|
||||
-I$(LIB_SRC)/TurbulenceModels/turbulenceModels/lnInclude \
|
||||
-I$(LIB_SRC)/TurbulenceModels/compressible/lnInclude \
|
||||
-I$(LIB_SRC)/lagrangian/basic/lnInclude \
|
||||
-I$(LIB_SRC)/lagrangian/intermediate/lnInclude \
|
||||
-I$(LIB_SRC)/lagrangian/coalCombustion/lnInclude \
|
||||
-I$(LIB_SRC)/lagrangian/distributionModels/lnInclude \
|
||||
-I$(LIB_SRC)/thermophysicalModels/specie/lnInclude \
|
||||
-I$(LIB_SRC)/transportModels/compressible/lnInclude \
|
||||
-I$(LIB_SRC)/thermophysicalModels/basic/lnInclude \
|
||||
-I$(LIB_SRC)/thermophysicalModels/properties/liquidProperties/lnInclude \
|
||||
-I$(LIB_SRC)/thermophysicalModels/properties/liquidMixtureProperties/lnInclude \
|
||||
-I$(LIB_SRC)/thermophysicalModels/properties/solidProperties/lnInclude \
|
||||
-I$(LIB_SRC)/thermophysicalModels/properties/solidMixtureProperties/lnInclude \
|
||||
-I$(LIB_SRC)/thermophysicalModels/thermophysicalFunctions/lnInclude \
|
||||
-I$(LIB_SRC)/thermophysicalModels/reactionThermo/lnInclude \
|
||||
-I$(LIB_SRC)/thermophysicalModels/SLGThermo/lnInclude \
|
||||
-I$(LIB_SRC)/thermophysicalModels/chemistryModel/lnInclude \
|
||||
-I$(LIB_SRC)/thermophysicalModels/radiation/lnInclude \
|
||||
-I$(LIB_SRC)/ODE/lnInclude \
|
||||
-I$(LIB_SRC)/regionModels/regionModel/lnInclude \
|
||||
-I$(LIB_SRC)/regionModels/surfaceFilmModels/lnInclude \
|
||||
-I$(LIB_SRC)/combustionModels/lnInclude \
|
||||
-I$(LIB_SRC)/sampling/lnInclude \
|
||||
-I$(FOAM_SOLVERS)/combustion/reactingFoam
|
||||
|
||||
|
||||
EXE_LIBS = \
|
||||
-lfiniteVolume \
|
||||
-lmeshTools \
|
||||
-lturbulenceModels \
|
||||
-lcompressibleTurbulenceModels \
|
||||
-llagrangian \
|
||||
-llagrangianIntermediate \
|
||||
-lcoalCombustion \
|
||||
-llagrangianTurbulence \
|
||||
-lspecie \
|
||||
-lcompressibleTransportModels \
|
||||
-lfluidThermophysicalModels \
|
||||
-lliquidProperties \
|
||||
-lliquidMixtureProperties \
|
||||
-lsolidProperties \
|
||||
-lsolidMixtureProperties \
|
||||
-lthermophysicalFunctions \
|
||||
-lreactionThermophysicalModels \
|
||||
-lSLGThermo \
|
||||
-lchemistryModel \
|
||||
-lradiationModels \
|
||||
-lODE \
|
||||
-lregionModels \
|
||||
-lsurfaceFilmModels \
|
||||
-lcombustionModels \
|
||||
-lfvOptions \
|
||||
-lsampling
|
||||
|
|
@ -0,0 +1,21 @@
|
|||
MRF.correctBoundaryVelocity(U);
|
||||
|
||||
tmp<fvVectorMatrix> tUEqn
|
||||
(
|
||||
fvm::div(phi, U)
|
||||
+ MRF.DDt(rho, U)
|
||||
+ turbulence->divDevRhoReff(U)
|
||||
==
|
||||
rho()*g
|
||||
+ parcels.SU(U)
|
||||
+ fvOptions(rho, U)
|
||||
);
|
||||
fvVectorMatrix& UEqn = tUEqn.ref();
|
||||
|
||||
UEqn.relax();
|
||||
|
||||
fvOptions.constrain(UEqn);
|
||||
|
||||
solve(UEqn == -fvc::grad(p));
|
||||
|
||||
fvOptions.correct(U);
|
||||
|
|
@ -0,0 +1,53 @@
|
|||
tmp<fv::convectionScheme<scalar>> mvConvection
|
||||
(
|
||||
fv::convectionScheme<scalar>::New
|
||||
(
|
||||
mesh,
|
||||
fields,
|
||||
phi,
|
||||
mesh.divScheme("div(phi,Yi_h)")
|
||||
)
|
||||
);
|
||||
|
||||
{
|
||||
combustion->correct();
|
||||
dQ = combustion->dQ();
|
||||
label inertIndex = -1;
|
||||
volScalarField Yt(0.0*Y[0]);
|
||||
|
||||
forAll(Y, i)
|
||||
{
|
||||
if (Y[i].name() != inertSpecie)
|
||||
{
|
||||
volScalarField& Yi = Y[i];
|
||||
|
||||
fvScalarMatrix YEqn
|
||||
(
|
||||
mvConvection->fvmDiv(phi, Yi)
|
||||
- fvm::laplacian(turbulence->muEff(), Yi)
|
||||
==
|
||||
parcels.Srho(i)
|
||||
+ combustion->R(Yi)
|
||||
+ fvOptions(rho, Yi)
|
||||
);
|
||||
|
||||
YEqn.relax();
|
||||
|
||||
fvOptions.constrain(YEqn);
|
||||
|
||||
YEqn.solve(mesh.solver("Yi"));
|
||||
|
||||
fvOptions.correct(Yi);
|
||||
|
||||
Yi.max(0.0);
|
||||
Yt += Yi;
|
||||
}
|
||||
else
|
||||
{
|
||||
inertIndex = i;
|
||||
}
|
||||
}
|
||||
|
||||
Y[inertIndex] = scalar(1) - Yt;
|
||||
Y[inertIndex].max(0.0);
|
||||
}
|
||||
|
|
@ -0,0 +1,8 @@
|
|||
Info<< "\nConstructing reacting cloud" << endl;
|
||||
coalCloudList parcels
|
||||
(
|
||||
rho,
|
||||
U,
|
||||
g,
|
||||
slgThermo
|
||||
);
|
||||
|
|
@ -0,0 +1,2 @@
|
|||
const volScalarField& T = thermo.T();
|
||||
const volScalarField& psi = thermo.psi();
|
||||
|
|
@ -0,0 +1,123 @@
|
|||
#include "readGravitationalAcceleration.H"
|
||||
|
||||
Info<< "Creating combustion model\n" << endl;
|
||||
|
||||
autoPtr<combustionModels::rhoCombustionModel> combustion
|
||||
(
|
||||
combustionModels::rhoCombustionModel::New(mesh)
|
||||
);
|
||||
|
||||
rhoReactionThermo& thermo = combustion->thermo();
|
||||
thermo.validate(args.executable(), "h", "e");
|
||||
|
||||
SLGThermo slgThermo(mesh, thermo);
|
||||
|
||||
basicSpecieMixture& composition = thermo.composition();
|
||||
PtrList<volScalarField>& Y = composition.Y();
|
||||
|
||||
const word inertSpecie(thermo.lookup("inertSpecie"));
|
||||
|
||||
if (!composition.contains(inertSpecie))
|
||||
{
|
||||
FatalErrorInFunction
|
||||
<< "Specified inert specie '" << inertSpecie << "' not found in "
|
||||
<< "species list. Available species:" << composition.species()
|
||||
<< exit(FatalError);
|
||||
}
|
||||
|
||||
volScalarField& p = thermo.p();
|
||||
|
||||
volScalarField rho
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"rho",
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::NO_READ,
|
||||
IOobject::AUTO_WRITE
|
||||
),
|
||||
thermo.rho()
|
||||
);
|
||||
|
||||
Info<< "\nReading field U\n" << endl;
|
||||
volVectorField U
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"U",
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::MUST_READ,
|
||||
IOobject::AUTO_WRITE
|
||||
),
|
||||
mesh
|
||||
);
|
||||
|
||||
#include "compressibleCreatePhi.H"
|
||||
|
||||
mesh.setFluxRequired(p.name());
|
||||
|
||||
dimensionedScalar rhoMax
|
||||
(
|
||||
dimensionedScalar::lookupOrDefault
|
||||
(
|
||||
"rhoMax",
|
||||
simple.dict(),
|
||||
dimDensity,
|
||||
GREAT
|
||||
)
|
||||
);
|
||||
|
||||
dimensionedScalar rhoMin
|
||||
(
|
||||
dimensionedScalar::lookupOrDefault
|
||||
(
|
||||
"rhoMin",
|
||||
simple.dict(),
|
||||
dimDensity,
|
||||
0
|
||||
)
|
||||
);
|
||||
|
||||
Info<< "Creating turbulence model\n" << endl;
|
||||
autoPtr<compressible::turbulenceModel> turbulence
|
||||
(
|
||||
compressible::turbulenceModel::New
|
||||
(
|
||||
rho,
|
||||
U,
|
||||
phi,
|
||||
thermo
|
||||
)
|
||||
);
|
||||
|
||||
// Set the turbulence into the combustion model
|
||||
combustion->setTurbulence(turbulence());
|
||||
|
||||
Info<< "Creating multi-variate interpolation scheme\n" << endl;
|
||||
multivariateSurfaceInterpolationScheme<scalar>::fieldTable fields;
|
||||
|
||||
forAll(Y, i)
|
||||
{
|
||||
fields.add(Y[i]);
|
||||
}
|
||||
fields.add(thermo.he());
|
||||
|
||||
volScalarField dQ
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"dQ",
|
||||
runTime.timeName(),
|
||||
mesh,
|
||||
IOobject::NO_READ,
|
||||
IOobject::AUTO_WRITE
|
||||
),
|
||||
mesh,
|
||||
dimensionedScalar("dQ", dimEnergy/dimTime, 0.0)
|
||||
);
|
||||
|
||||
#include "createMRF.H"
|
||||
#include "createRadiationModel.H"
|
||||
#include "createClouds.H"
|
||||
|
|
@ -0,0 +1,57 @@
|
|||
{
|
||||
// Thermodynamic density needs to be updated by psi*d(p) after the
|
||||
// pressure solution - done in 2 parts. Part 1:
|
||||
thermo.rho() -= psi*p;
|
||||
|
||||
volScalarField rAU(1.0/UEqn.A());
|
||||
surfaceScalarField rhorAUf("rhorAUf", fvc::interpolate(rho*rAU));
|
||||
volVectorField HbyA(constrainHbyA(rAU*UEqn.H(), U, p));
|
||||
tUEqn.clear();
|
||||
surfaceScalarField phiHbyA
|
||||
(
|
||||
"phiHbyA",
|
||||
fvc::interpolate(rho)*fvc::flux(HbyA)
|
||||
);
|
||||
|
||||
MRF.makeRelative(fvc::interpolate(rho), phiHbyA);
|
||||
|
||||
// Update the pressure BCs to ensure flux consistency
|
||||
constrainPressure(p, rho, U, phiHbyA, rhorAUf, MRF);
|
||||
|
||||
while (simple.correctNonOrthogonal())
|
||||
{
|
||||
fvScalarMatrix pEqn
|
||||
(
|
||||
fvc::div(phiHbyA)
|
||||
- fvm::laplacian(rhorAUf, p)
|
||||
==
|
||||
parcels.Srho()
|
||||
+ fvOptions(psi, p, rho.name())
|
||||
);
|
||||
|
||||
pEqn.solve();
|
||||
|
||||
if (simple.finalNonOrthogonalIter())
|
||||
{
|
||||
phi = phiHbyA + pEqn.flux();
|
||||
}
|
||||
}
|
||||
|
||||
p.relax();
|
||||
|
||||
// Second part of thermodynamic density update
|
||||
thermo.rho() += psi*p;
|
||||
|
||||
#include "compressibleContinuityErrs.H"
|
||||
|
||||
U = HbyA - rAU*fvc::grad(p);
|
||||
U.correctBoundaryConditions();
|
||||
fvOptions.correct(U);
|
||||
|
||||
rho = thermo.rho();
|
||||
rho = max(rho, rhoMin);
|
||||
rho = min(rho, rhoMax);
|
||||
rho.relax();
|
||||
|
||||
Info<< "p min/max = " << min(p).value() << ", " << max(p).value() << endl;
|
||||
}
|
||||
|
|
@ -22,7 +22,7 @@ License
|
|||
along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
|
||||
|
||||
Application
|
||||
simpleReactingParcelFoam
|
||||
simpleCoalcombustionFoam
|
||||
|
||||
Description
|
||||
Steady state solver for compressible, turbulent flow with reacting,
|
||||
|
|
@ -31,19 +31,31 @@ Description
|
|||
\*---------------------------------------------------------------------------*/
|
||||
|
||||
#include "fvCFD.H"
|
||||
#include "turbulentFluidThermoModel.H"
|
||||
#include "coalCloudList.H"
|
||||
#include "rhoCombustionModel.H"
|
||||
#include "radiationModel.H"
|
||||
#include "IOporosityModelList.H"
|
||||
#include "fvOptions.H"
|
||||
#include "SLGThermo.H"
|
||||
#include "simpleControl.H"
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
int main(int argc, char *argv[])
|
||||
{
|
||||
#include "postProcess.H"
|
||||
|
||||
#include "setRootCase.H"
|
||||
#include "createTime.H"
|
||||
#include "createMesh.H"
|
||||
#include "createControl.H"
|
||||
#include "createFields.H"
|
||||
#include "createFieldRefs.H"
|
||||
#include "createFvOptions.H"
|
||||
#include "initContinuityErrs.H"
|
||||
|
||||
turbulence->validate();
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
|
|
@ -53,25 +65,18 @@ int main(int argc, char *argv[])
|
|||
{
|
||||
Info<< "Time = " << runTime.timeName() << nl << endl;
|
||||
|
||||
#include "calculateSso2.H"
|
||||
parcels.evolve();
|
||||
|
||||
while (simple.correctNonOrthogonal())
|
||||
// --- Pressure-velocity SIMPLE corrector loop
|
||||
{
|
||||
fvScalarMatrix SO2Eqn
|
||||
(
|
||||
fvm::ddt(rho, SO2)
|
||||
+ fvm::div(phi, SO2)
|
||||
- fvm::laplacian(mut, SO2)
|
||||
==
|
||||
Sso2
|
||||
);
|
||||
|
||||
SO2Eqn.relax();
|
||||
fvOptions.constrain(SO2Eqn);
|
||||
SO2Eqn.solve();
|
||||
fvOptions.correct(SO2);
|
||||
#include "UEqn.H"
|
||||
#include "YEqn.H"
|
||||
#include "EEqn.H"
|
||||
#include "pEqn.H"
|
||||
}
|
||||
|
||||
turbulence->correct();
|
||||
|
||||
runTime.write();
|
||||
|
||||
Info<< "ExecutionTime = " << runTime.elapsedCpuTime() << " s"
|
||||
227
src/combustionModels/EDM/EDM.C
Normal file
227
src/combustionModels/EDM/EDM.C
Normal file
|
|
@ -0,0 +1,227 @@
|
|||
/*---------------------------------------------------------------------------*\
|
||||
========= |
|
||||
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
||||
\\ / O peration |
|
||||
\\ / A nd | Copyright (C) 2013-2016 OpenFOAM Foundation
|
||||
\\/ M anipulation |
|
||||
-------------------------------------------------------------------------------
|
||||
License
|
||||
This file is part of OpenFOAM.
|
||||
|
||||
OpenFOAM is free software: you can redistribute it and/or modify it
|
||||
under the terms of the GNU General Public License as published by
|
||||
the Free Software Foundation, either version 3 of the License, or
|
||||
(at your option) any later version.
|
||||
|
||||
OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
|
||||
ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
|
||||
FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
|
||||
for more details.
|
||||
|
||||
You should have received a copy of the GNU General Public License
|
||||
along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
|
||||
|
||||
\*---------------------------------------------------------------------------*/
|
||||
|
||||
#include "EDM.H"
|
||||
#include "fvmSup.H"
|
||||
|
||||
// * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
|
||||
|
||||
|
||||
template<class Type>
|
||||
Foam::combustionModels::EDM<Type>::EDM
|
||||
(
|
||||
const word& modelType,
|
||||
const fvMesh& mesh,
|
||||
const word& phaseName
|
||||
)
|
||||
:
|
||||
Type(modelType, mesh, phaseName),
|
||||
finiteRate_ (this->coeffs().lookupOrDefault("finiteRate", false)),
|
||||
A_ (this->coeffs().lookupOrDefault("A", 4.0)),
|
||||
B_ (this->coeffs().lookupOrDefault("B", 0.5))
|
||||
{
|
||||
if (finiteRate_)
|
||||
{
|
||||
Info<< " using Finite-rate/Eddy Dissipation Model" << endl;
|
||||
}
|
||||
else
|
||||
{
|
||||
Info<< " using Eddy Dissipation Model" << endl;
|
||||
}
|
||||
Info<< " A = " << A_ << endl;
|
||||
Info<< " B = " << B_ << endl;
|
||||
}
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * * Destructor * * * * * * * * * * * * * * * //
|
||||
|
||||
template<class Type>
|
||||
Foam::combustionModels::EDM<Type>::~EDM()
|
||||
{}
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * * //
|
||||
|
||||
template<class Type>
|
||||
Foam::tmp<Foam::volScalarField>
|
||||
Foam::combustionModels::EDM<Type>::tc() const
|
||||
{
|
||||
return this->chemistryPtr_->tc();
|
||||
}
|
||||
|
||||
|
||||
template<class Type>
|
||||
void Foam::combustionModels::EDM<Type>::correct()
|
||||
{
|
||||
if (this->active())
|
||||
{
|
||||
tmp<volScalarField> tk(this->turbulence().k());
|
||||
const volScalarField& k = tk();
|
||||
|
||||
tmp<volScalarField> tepsilon(this->turbulence().epsilon());
|
||||
const volScalarField& epsilon = tepsilon();
|
||||
|
||||
this->chemistryPtr_->edm(epsilon/k, A_, B_);
|
||||
|
||||
if (finiteRate_)
|
||||
{
|
||||
// Copy edm rate stored in chemistryPtr_->RR(i)
|
||||
PtrList<DimensionedField<scalar, volMesh>> edmRates
|
||||
(
|
||||
this->thermo().composition().Y().size()
|
||||
);
|
||||
forAll(edmRates, i)
|
||||
{
|
||||
edmRates.set
|
||||
(
|
||||
i,
|
||||
new DimensionedField<scalar, volMesh>
|
||||
(
|
||||
this->chemistryPtr_->RR(i)
|
||||
)
|
||||
);
|
||||
}
|
||||
|
||||
this->chemistryPtr_->calculate();
|
||||
|
||||
forAll(edmRates, i)
|
||||
{
|
||||
DimensionedField<scalar, volMesh> &finiteRatei
|
||||
(
|
||||
this->chemistryPtr_->RR(i)
|
||||
);
|
||||
|
||||
forAll(finiteRatei, celli)
|
||||
{
|
||||
finiteRatei[celli] = min(finiteRatei[celli], edmRates[i][celli]);
|
||||
finiteRatei[celli] =
|
||||
finiteRatei[celli]*edmRates[i][celli]
|
||||
/ (finiteRatei[celli] + edmRates[i][celli] + SMALL);
|
||||
}
|
||||
|
||||
}
|
||||
}
|
||||
}
|
||||
}
|
||||
|
||||
|
||||
template<class Type>
|
||||
Foam::tmp<Foam::fvScalarMatrix>
|
||||
Foam::combustionModels::EDM<Type>::R(volScalarField& Y) const
|
||||
{
|
||||
tmp<fvScalarMatrix> tSu(new fvScalarMatrix(Y, dimMass/dimTime));
|
||||
|
||||
fvScalarMatrix& Su = tSu.ref();
|
||||
|
||||
if (this->active())
|
||||
{
|
||||
const label specieI =
|
||||
this->thermo().composition().species()[Y.member()];
|
||||
|
||||
Su += this->chemistryPtr_->RR(specieI);
|
||||
}
|
||||
|
||||
return tSu;
|
||||
}
|
||||
|
||||
|
||||
template<class Type>
|
||||
Foam::tmp<Foam::volScalarField>
|
||||
Foam::combustionModels::EDM<Type>::dQ() const
|
||||
{
|
||||
tmp<volScalarField> tdQ
|
||||
(
|
||||
new volScalarField
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
IOobject::groupName(typeName + ":dQ", this->phaseName_),
|
||||
this->mesh().time().timeName(),
|
||||
this->mesh(),
|
||||
IOobject::NO_READ,
|
||||
IOobject::NO_WRITE,
|
||||
false
|
||||
),
|
||||
this->mesh(),
|
||||
dimensionedScalar("dQ", dimEnergy/dimTime, 0.0)
|
||||
)
|
||||
);
|
||||
|
||||
if (this->active())
|
||||
{
|
||||
tdQ.ref() = this->chemistryPtr_->dQ();
|
||||
}
|
||||
|
||||
return tdQ;
|
||||
}
|
||||
|
||||
|
||||
template<class Type>
|
||||
Foam::tmp<Foam::volScalarField>
|
||||
Foam::combustionModels::EDM<Type>::Sh() const
|
||||
{
|
||||
tmp<volScalarField> tSh
|
||||
(
|
||||
new volScalarField
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
IOobject::groupName(typeName + ":Sh", this->phaseName_),
|
||||
this->mesh().time().timeName(),
|
||||
this->mesh(),
|
||||
IOobject::NO_READ,
|
||||
IOobject::NO_WRITE,
|
||||
false
|
||||
),
|
||||
this->mesh(),
|
||||
dimensionedScalar("zero", dimEnergy/dimTime/dimVolume, 0.0)
|
||||
)
|
||||
);
|
||||
|
||||
if (this->active())
|
||||
{
|
||||
tSh.ref() = this->chemistryPtr_->Sh();
|
||||
}
|
||||
|
||||
return tSh;
|
||||
}
|
||||
|
||||
|
||||
template<class Type>
|
||||
bool Foam::combustionModels::EDM<Type>::read()
|
||||
{
|
||||
if (Type::read())
|
||||
{
|
||||
this->coeffs().lookup("finiteRate") >> finiteRate_;
|
||||
return true;
|
||||
}
|
||||
else
|
||||
{
|
||||
return false;
|
||||
}
|
||||
}
|
||||
|
||||
|
||||
// ************************************************************************* //
|
||||
147
src/combustionModels/EDM/EDM.H
Normal file
147
src/combustionModels/EDM/EDM.H
Normal file
|
|
@ -0,0 +1,147 @@
|
|||
/*---------------------------------------------------------------------------*\
|
||||
========= |
|
||||
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
||||
\\ / O peration |
|
||||
\\ / A nd | Copyright (C) 2013-2016 OpenFOAM Foundation
|
||||
\\/ M anipulation |
|
||||
-------------------------------------------------------------------------------
|
||||
License
|
||||
This file is part of OpenFOAM.
|
||||
|
||||
OpenFOAM is free software: you can redistribute it and/or modify it
|
||||
under the terms of the GNU General Public License as published by
|
||||
the Free Software Foundation, either version 3 of the License, or
|
||||
(at your option) any later version.
|
||||
|
||||
OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
|
||||
ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
|
||||
FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
|
||||
for more details.
|
||||
|
||||
You should have received a copy of the GNU General Public License
|
||||
along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
|
||||
|
||||
Class
|
||||
Foam::combustionModels::EDM
|
||||
|
||||
Description
|
||||
Eddy dissipation model with finite-rate chemistry.
|
||||
28.Feb.2017
|
||||
Combustion Lab. POSTECH
|
||||
Karam Han
|
||||
(Version upgraded by Jinwoo Park)
|
||||
|
||||
SourceFiles
|
||||
EDM.C
|
||||
|
||||
\*---------------------------------------------------------------------------*/
|
||||
|
||||
#ifndef EDM_H
|
||||
#define EDM_H
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
namespace Foam
|
||||
{
|
||||
namespace combustionModels
|
||||
{
|
||||
|
||||
/*---------------------------------------------------------------------------*\
|
||||
Class EDM Declaration
|
||||
\*---------------------------------------------------------------------------*/
|
||||
|
||||
template<class Type>
|
||||
class EDM
|
||||
:
|
||||
public Type
|
||||
{
|
||||
// Private data
|
||||
|
||||
//- finite-rate/EDM using Arrhenius formula
|
||||
bool finiteRate_;
|
||||
|
||||
//- EDM Parameter A
|
||||
scalar A_;
|
||||
|
||||
//- EDM Parameter B
|
||||
scalar B_;
|
||||
|
||||
protected:
|
||||
|
||||
// Protected Member Functions
|
||||
|
||||
//- Return the chemical time scale
|
||||
tmp<volScalarField> tc() const;
|
||||
|
||||
private:
|
||||
|
||||
// Private Member Functions
|
||||
|
||||
//- Disallow copy construct
|
||||
EDM(const EDM&);
|
||||
|
||||
//- Disallow default bitwise assignment
|
||||
void operator=(const EDM&);
|
||||
|
||||
|
||||
public:
|
||||
|
||||
//- Runtime type information
|
||||
TypeName("EDM");
|
||||
|
||||
|
||||
// Constructors
|
||||
|
||||
//- Construct from components
|
||||
EDM
|
||||
(
|
||||
const word& modelType,
|
||||
const fvMesh& mesh,
|
||||
const word& phaseName
|
||||
);
|
||||
|
||||
|
||||
//- Destructor
|
||||
virtual ~EDM();
|
||||
|
||||
|
||||
// Member Functions
|
||||
|
||||
// Evolution
|
||||
|
||||
//- Correct combustion rate
|
||||
virtual void correct();
|
||||
|
||||
//- Fuel consumption rate matrix.
|
||||
virtual tmp<fvScalarMatrix> R(volScalarField& Y) const;
|
||||
|
||||
//- Heat release rate calculated from fuel consumption rate matrix
|
||||
virtual tmp<volScalarField> dQ() const;
|
||||
|
||||
//- Return source for enthalpy equation [kg/m/s3]
|
||||
virtual tmp<volScalarField> Sh() const;
|
||||
|
||||
|
||||
// IO
|
||||
|
||||
//- Update properties from given dictionary
|
||||
virtual bool read();
|
||||
};
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
} // End namespace combustionModels
|
||||
} // End namespace Foam
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
#ifdef NoRepository
|
||||
#include "EDM.C"
|
||||
#endif
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
#endif
|
||||
|
||||
// ************************************************************************* //
|
||||
38
src/combustionModels/EDM/EDMs.C
Normal file
38
src/combustionModels/EDM/EDMs.C
Normal file
|
|
@ -0,0 +1,38 @@
|
|||
/*---------------------------------------------------------------------------*\
|
||||
========= |
|
||||
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
||||
\\ / O peration |
|
||||
\\ / A nd | Copyright (C) 2013-2015 OpenFOAM Foundation
|
||||
\\/ M anipulation |
|
||||
-------------------------------------------------------------------------------
|
||||
License
|
||||
This file is part of OpenFOAM.
|
||||
|
||||
OpenFOAM is free software: you can redistribute it and/or modify it
|
||||
under the terms of the GNU General Public License as published by
|
||||
the Free Software Foundation, either version 3 of the License, or
|
||||
(at your option) any later version.
|
||||
|
||||
OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
|
||||
ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
|
||||
FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
|
||||
for more details.
|
||||
|
||||
You should have received a copy of the GNU General Public License
|
||||
along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
|
||||
|
||||
\*---------------------------------------------------------------------------*/
|
||||
|
||||
#include "makeCombustionTypes.H"
|
||||
|
||||
#include "psiChemistryCombustion.H"
|
||||
#include "rhoChemistryCombustion.H"
|
||||
#include "EDM.H"
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
makeCombustionTypes(EDM, psiChemistryCombustion, psiCombustionModel);
|
||||
makeCombustionTypes(EDM, rhoChemistryCombustion, rhoCombustionModel);
|
||||
|
||||
|
||||
// ************************************************************************* //
|
||||
|
|
@ -13,6 +13,8 @@ rhoCombustionModel/rhoChemistryCombustion/rhoChemistryCombustion.C
|
|||
diffusion/diffusions.C
|
||||
infinitelyFastChemistry/infinitelyFastChemistrys.C
|
||||
|
||||
EDM/EDMs.C
|
||||
|
||||
PaSR/PaSRs.C
|
||||
|
||||
laminar/laminars.C
|
||||
|
|
|
|||
|
|
@ -34,6 +34,7 @@ License
|
|||
#include "COxidationKineticDiffusionLimitedRate.H"
|
||||
#include "COxidationHurtMitchell.H"
|
||||
#include "COxidationMurphyShaddix.H"
|
||||
#include "COxidationKD_CO.H"
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
|
|
@ -47,6 +48,7 @@ License
|
|||
); \
|
||||
makeSurfaceReactionModelType(COxidationIntrinsicRate, CloudType); \
|
||||
makeSurfaceReactionModelType(COxidationHurtMitchell, CloudType); \
|
||||
makeSurfaceReactionModelType(COxidationKD_CO, CloudType); \
|
||||
makeSurfaceReactionModelType(COxidationMurphyShaddix, CloudType);
|
||||
|
||||
|
||||
|
|
|
|||
|
|
@ -0,0 +1,178 @@
|
|||
/*---------------------------------------------------------------------------*\
|
||||
========= |
|
||||
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
||||
\\ / O peration |
|
||||
\\ / A nd | Copyright (C) 2011 OpenFOAM Foundation
|
||||
\\/ M anipulation |
|
||||
-------------------------------------------------------------------------------
|
||||
License
|
||||
This file is part of OpenFOAM.
|
||||
|
||||
OpenFOAM is free software: you can redistribute it and/or modify it
|
||||
under the terms of the GNU General Public License as published by
|
||||
the Free Software Foundation, either version 3 of the License, or
|
||||
(at your option) any later version.
|
||||
|
||||
OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
|
||||
ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
|
||||
FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
|
||||
for more details.
|
||||
|
||||
You should have received a copy of the GNU General Public License
|
||||
along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
|
||||
|
||||
\*---------------------------------------------------------------------------*/
|
||||
|
||||
#include "COxidationKD_CO.H"
|
||||
#include "mathematicalConstants.H"
|
||||
#include "thermodynamicConstants.H"
|
||||
|
||||
// * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
|
||||
|
||||
template<class CloudType>
|
||||
Foam::COxidationKD_CO<CloudType>::
|
||||
COxidationKD_CO
|
||||
(
|
||||
const dictionary& dict,
|
||||
CloudType& owner
|
||||
)
|
||||
:
|
||||
SurfaceReactionModel<CloudType>(dict, owner, typeName),
|
||||
Sb_(readScalar(this->coeffDict().lookup("Sb"))),
|
||||
C1_(readScalar(this->coeffDict().lookup("C1"))),
|
||||
C2_(readScalar(this->coeffDict().lookup("C2"))),
|
||||
E_(readScalar(this->coeffDict().lookup("E"))),
|
||||
CsLocalId_(-1),
|
||||
O2GlobalId_(owner.composition().carrierId("O2")),
|
||||
COGlobalId_(owner.composition().carrierId("CO")),
|
||||
WC_(0.0),
|
||||
WO2_(0.0),
|
||||
HcCO_(0.0)
|
||||
{
|
||||
// Determine Cs ids
|
||||
label idSolid = owner.composition().idSolid();
|
||||
CsLocalId_ = owner.composition().localId(idSolid, "C");
|
||||
|
||||
// Set local copies of thermo properties
|
||||
WO2_ = owner.thermo().carrier().W(O2GlobalId_);
|
||||
const scalar WCO = owner.thermo().carrier().W(COGlobalId_);
|
||||
WC_ = WCO - 0.5*WO2_;
|
||||
|
||||
HcCO_ = owner.thermo().carrier().Hc(COGlobalId_);
|
||||
|
||||
const scalar YCloc = owner.composition().Y0(idSolid)[CsLocalId_];
|
||||
const scalar YSolidTot = owner.composition().YMixture0()[idSolid];
|
||||
Info<< " C(s): particle mass fraction = " << YCloc*YSolidTot << endl;
|
||||
}
|
||||
|
||||
|
||||
template<class CloudType>
|
||||
Foam::COxidationKD_CO<CloudType>::
|
||||
COxidationKD_CO
|
||||
(
|
||||
const COxidationKD_CO<CloudType>& srm
|
||||
)
|
||||
:
|
||||
SurfaceReactionModel<CloudType>(srm),
|
||||
Sb_(srm.Sb_),
|
||||
C1_(srm.C1_),
|
||||
C2_(srm.C2_),
|
||||
E_(srm.E_),
|
||||
CsLocalId_(srm.CsLocalId_),
|
||||
O2GlobalId_(srm.O2GlobalId_),
|
||||
COGlobalId_(srm.COGlobalId_),
|
||||
WC_(srm.WC_),
|
||||
WO2_(srm.WO2_),
|
||||
HcCO_(srm.HcCO_)
|
||||
{}
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * * Destructor * * * * * * * * * * * * * * * //
|
||||
|
||||
template<class CloudType>
|
||||
Foam::COxidationKD_CO<CloudType>::
|
||||
~COxidationKD_CO()
|
||||
{}
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
|
||||
|
||||
template<class CloudType>
|
||||
Foam::scalar Foam::COxidationKD_CO<CloudType>::calculate
|
||||
(
|
||||
const scalar dt,
|
||||
const label cellI,
|
||||
const scalar d,
|
||||
const scalar T,
|
||||
const scalar Tc,
|
||||
const scalar pc,
|
||||
const scalar rhoc,
|
||||
const scalar mass,
|
||||
const scalarField& YGas,
|
||||
const scalarField& YLiquid,
|
||||
const scalarField& YSolid,
|
||||
const scalarField& YMixture,
|
||||
const scalar N,
|
||||
scalarField& dMassGas,
|
||||
scalarField& dMassLiquid,
|
||||
scalarField& dMassSolid,
|
||||
scalarField& dMassSRCarrier
|
||||
) const
|
||||
{
|
||||
// Fraction of remaining combustible material
|
||||
const label idSolid = CloudType::parcelType::SLD;
|
||||
const scalar fComb = YMixture[idSolid]*YSolid[CsLocalId_];
|
||||
|
||||
// Surface combustion active combustible fraction is consumed
|
||||
if (fComb < SMALL)
|
||||
{
|
||||
return 0.0;
|
||||
}
|
||||
|
||||
const SLGThermo& thermo = this->owner().thermo();
|
||||
|
||||
// Local mass fraction of O2 in the carrier phase
|
||||
const scalar YO2 = thermo.carrier().Y(O2GlobalId_)[cellI];
|
||||
|
||||
// Diffusion rate coefficient
|
||||
const scalar D0 = C1_/d*pow(0.5*(T + Tc), 0.75);
|
||||
|
||||
// Kinetic rate
|
||||
const scalar Rk = C2_*exp(-E_/(constant::thermodynamic::RR*T)); //
|
||||
|
||||
// Particle surface area
|
||||
const scalar Ap = constant::mathematical::pi*sqr(d);
|
||||
|
||||
// Change in C mass [kg]
|
||||
scalar dmC = Ap*rhoc*constant::thermodynamic::RR*Tc*YO2/WO2_*D0*Rk/(D0 + Rk)*dt;
|
||||
|
||||
// Limit mass transfer by availability of C
|
||||
dmC = min(mass*fComb, dmC);
|
||||
|
||||
// Molar consumption
|
||||
const scalar dOmega = dmC/WC_;
|
||||
|
||||
// Change in O2 mass [kg]
|
||||
const scalar dmO2 = dOmega*Sb_*WO2_;
|
||||
|
||||
// Mass of newly created CO [kg]
|
||||
const scalar dmCO = dOmega*(WC_ + Sb_*WO2_);
|
||||
|
||||
// Update local particle C mass
|
||||
dMassSolid[CsLocalId_] += dOmega*WC_;
|
||||
|
||||
// Update carrier O2 and CO mass
|
||||
dMassSRCarrier[O2GlobalId_] -= dmO2;
|
||||
dMassSRCarrier[COGlobalId_] += dmCO;
|
||||
|
||||
const scalar HsC = thermo.solids().properties()[CsLocalId_].Hs(T);
|
||||
|
||||
// carrier sensible enthalpy exchange handled via change in mass
|
||||
|
||||
// Heat of reaction [J]
|
||||
//Info<<"Heat of reaction(Char)"<<tab<<(dmC*HsC - dmCO*HcCO_)/dmC<<endl;
|
||||
return dmC*HsC - dmCO*HcCO_;
|
||||
}
|
||||
|
||||
|
||||
// ************************************************************************* //
|
||||
|
|
@ -0,0 +1,176 @@
|
|||
/*---------------------------------------------------------------------------*\
|
||||
========= |
|
||||
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
||||
\\ / O peration |
|
||||
\\ / A nd | Copyright (C) 2011 OpenFOAM Foundation
|
||||
\\/ M anipulation |
|
||||
-------------------------------------------------------------------------------
|
||||
License
|
||||
This file is part of OpenFOAM.
|
||||
|
||||
OpenFOAM is free software: you can redistribute it and/or modify it
|
||||
under the terms of the GNU General Public License as published by
|
||||
the Free Software Foundation, either version 3 of the License, or
|
||||
(at your option) any later version.
|
||||
|
||||
OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
|
||||
ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
|
||||
FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
|
||||
for more details.
|
||||
|
||||
You should have received a copy of the GNU General Public License
|
||||
along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
|
||||
|
||||
Class
|
||||
COxidationKD_CO
|
||||
|
||||
Description
|
||||
Kinetic/diffusion limited rate surface reaction model for coal parcels.
|
||||
Limited to:
|
||||
|
||||
C(s) + Sb*O2 -> CO
|
||||
|
||||
where Sb is the stoichiometry of the reaction
|
||||
|
||||
\*---------------------------------------------------------------------------*/
|
||||
|
||||
#ifndef COxidationKD_CO_H
|
||||
#define COxidationKD_CO_H
|
||||
|
||||
#include "SurfaceReactionModel.H"
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
namespace Foam
|
||||
{
|
||||
|
||||
// Forward class declarations
|
||||
template<class CloudType>
|
||||
class COxidationKD_CO;
|
||||
|
||||
/*---------------------------------------------------------------------------*\
|
||||
Class COxidationKD_CO Declaration
|
||||
\*---------------------------------------------------------------------------*/
|
||||
|
||||
template<class CloudType>
|
||||
class COxidationKD_CO
|
||||
:
|
||||
public SurfaceReactionModel<CloudType>
|
||||
{
|
||||
// Private data
|
||||
|
||||
// Model constants
|
||||
|
||||
//- Stoichiometry of reaction
|
||||
const scalar Sb_;
|
||||
|
||||
//- Mass diffusion limited rate constant, C1
|
||||
const scalar C1_;
|
||||
|
||||
//- Kinetics limited rate pre-exponential constant, C2
|
||||
const scalar C2_;
|
||||
|
||||
//- Kinetics limited rate activation energy
|
||||
const scalar E_;
|
||||
|
||||
|
||||
// Addressing
|
||||
|
||||
//- Cs positions in global/local lists
|
||||
label CsLocalId_;
|
||||
|
||||
//- O2 position in global list
|
||||
label O2GlobalId_;
|
||||
|
||||
//- CO2 positions in global list
|
||||
label COGlobalId_;
|
||||
|
||||
|
||||
// Local copies of thermo properties
|
||||
|
||||
//- Molecular weight of C [kg/kmol]
|
||||
scalar WC_;
|
||||
|
||||
//- Molecular weight of O2 [kg/kmol]
|
||||
scalar WO2_;
|
||||
|
||||
//- Formation enthalpy for CO2 [J/kg]
|
||||
scalar HcCO_;
|
||||
|
||||
|
||||
public:
|
||||
|
||||
//- Runtime type information
|
||||
TypeName("COxidationKD_CO");
|
||||
|
||||
|
||||
// Constructors
|
||||
|
||||
//- Construct from dictionary
|
||||
COxidationKD_CO
|
||||
(
|
||||
const dictionary& dict,
|
||||
CloudType& owner
|
||||
);
|
||||
|
||||
//- Construct copy
|
||||
COxidationKD_CO
|
||||
(
|
||||
const COxidationKD_CO<CloudType>& srm
|
||||
);
|
||||
|
||||
//- Construct and return a clone
|
||||
virtual autoPtr<SurfaceReactionModel<CloudType> > clone() const
|
||||
{
|
||||
return autoPtr<SurfaceReactionModel<CloudType> >
|
||||
(
|
||||
new COxidationKD_CO<CloudType>(*this)
|
||||
);
|
||||
}
|
||||
|
||||
|
||||
//- Destructor
|
||||
virtual ~COxidationKD_CO();
|
||||
|
||||
|
||||
// Member Functions
|
||||
|
||||
//- Update surface reactions
|
||||
virtual scalar calculate
|
||||
(
|
||||
const scalar dt,
|
||||
const label cellI,
|
||||
const scalar d,
|
||||
const scalar T,
|
||||
const scalar Tc,
|
||||
const scalar pc,
|
||||
const scalar rhoc,
|
||||
const scalar mass,
|
||||
const scalarField& YGas,
|
||||
const scalarField& YLiquid,
|
||||
const scalarField& YSolid,
|
||||
const scalarField& YMixture,
|
||||
const scalar N,
|
||||
scalarField& dMassGas,
|
||||
scalarField& dMassLiquid,
|
||||
scalarField& dMassSolid,
|
||||
scalarField& dMassSRCarrier
|
||||
) const;
|
||||
};
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
} // End namespace Foam
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
#ifdef NoRepository
|
||||
#include "COxidationKD_CO.C"
|
||||
#endif
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
#endif
|
||||
|
||||
// ************************************************************************* //
|
||||
|
|
@ -50,6 +50,12 @@ $(KINEMATICMPPICPARCEL)/defineBasicKinematicMPPICParcel.C
|
|||
$(KINEMATICMPPICPARCEL)/makeBasicKinematicMPPICParcelSubmodels.C
|
||||
|
||||
|
||||
/* thermo kinematic MPPIC parcel sub-models */
|
||||
THERMOMPPICPARCEL=$(DERIVEDPARCELS)/basicThermoKinematicMPPICParcel
|
||||
$(THERMOMPPICPARCEL)/defineBasicThermoKinematicMPPICParcel.C
|
||||
$(THERMOMPPICPARCEL)/makeBasicThermoKinematicMPPICParcelSubmodels.C
|
||||
|
||||
|
||||
/* bolt-on models */
|
||||
RADIATION=submodels/addOns/radiation
|
||||
$(RADIATION)/absorptionEmission/cloudAbsorptionEmission/cloudAbsorptionEmission.C
|
||||
|
|
|
|||
|
|
@ -213,6 +213,88 @@ Foam::ThermoCloud<CloudType>::ThermoCloud
|
|||
}
|
||||
|
||||
|
||||
template<class CloudType>
|
||||
Foam::ThermoCloud<CloudType>::ThermoCloud
|
||||
(
|
||||
const word& cloudName,
|
||||
const volScalarField& rho,
|
||||
const volVectorField& U,
|
||||
const volScalarField& mu,
|
||||
const dimensionedVector& g,
|
||||
bool readFields
|
||||
)
|
||||
:
|
||||
CloudType
|
||||
(
|
||||
cloudName,
|
||||
rho,
|
||||
U,
|
||||
mu,
|
||||
g,
|
||||
false
|
||||
),
|
||||
thermoCloud(),
|
||||
cloudCopyPtr_(NULL),
|
||||
constProps_(this->particleProperties()),
|
||||
thermo_(this->mesh().template lookupObject<SLGThermo>("SLGThermo")),
|
||||
T_(thermo_.thermo().T()),
|
||||
p_(thermo_.thermo().p()),
|
||||
heatTransferModel_(NULL),
|
||||
TIntegrator_(NULL),
|
||||
radiation_(false),
|
||||
radAreaP_(NULL),
|
||||
radT4_(NULL),
|
||||
radAreaPT4_(NULL),
|
||||
hsTrans_
|
||||
(
|
||||
new DimensionedField<scalar, volMesh>
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
this->name() + ":hsTrans",
|
||||
this->db().time().timeName(),
|
||||
this->db(),
|
||||
IOobject::READ_IF_PRESENT,
|
||||
IOobject::AUTO_WRITE
|
||||
),
|
||||
this->mesh(),
|
||||
dimensionedScalar("zero", dimEnergy, 0.0)
|
||||
)
|
||||
),
|
||||
hsCoeff_
|
||||
(
|
||||
new DimensionedField<scalar, volMesh>
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
this->name() + ":hsCoeff",
|
||||
this->db().time().timeName(),
|
||||
this->db(),
|
||||
IOobject::READ_IF_PRESENT,
|
||||
IOobject::AUTO_WRITE
|
||||
),
|
||||
this->mesh(),
|
||||
dimensionedScalar("zero", dimEnergy/dimTemperature, 0.0)
|
||||
)
|
||||
)
|
||||
{
|
||||
if (this->solution().active())
|
||||
{
|
||||
setModels();
|
||||
|
||||
if (readFields)
|
||||
{
|
||||
parcelType::readFields(*this);
|
||||
}
|
||||
}
|
||||
|
||||
if (this->solution().resetSourcesOnStartup())
|
||||
{
|
||||
resetSourceTerms();
|
||||
}
|
||||
}
|
||||
|
||||
|
||||
template<class CloudType>
|
||||
Foam::ThermoCloud<CloudType>::ThermoCloud
|
||||
(
|
||||
|
|
|
|||
|
|
@ -180,6 +180,17 @@ public:
|
|||
bool readFields = true
|
||||
);
|
||||
|
||||
//- Construct given carrier gas fields
|
||||
ThermoCloud
|
||||
(
|
||||
const word& cloudName,
|
||||
const volScalarField& rho,
|
||||
const volVectorField& U,
|
||||
const volScalarField& mu,
|
||||
const dimensionedVector& g,
|
||||
bool readFields = true
|
||||
);
|
||||
|
||||
//- Copy constructor with new name
|
||||
ThermoCloud(ThermoCloud<CloudType>& c, const word& name);
|
||||
|
||||
|
|
|
|||
|
|
@ -0,0 +1,66 @@
|
|||
/*---------------------------------------------------------------------------*\
|
||||
========= |
|
||||
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
||||
\\ / O peration |
|
||||
\\ / A nd | Copyright (C) 2013 OpenFOAM Foundation
|
||||
\\/ M anipulation |
|
||||
-------------------------------------------------------------------------------
|
||||
License
|
||||
This file is part of OpenFOAM.
|
||||
|
||||
OpenFOAM is free software: you can redistribute it and/or modify it
|
||||
under the terms of the GNU General Public License as published by
|
||||
the Free Software Foundation, either version 3 of the License, or
|
||||
(at your option) any later version.
|
||||
|
||||
OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
|
||||
ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
|
||||
FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
|
||||
for more details.
|
||||
|
||||
You should have received a copy of the GNU General Public License
|
||||
along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
|
||||
|
||||
Class
|
||||
Foam::basicThermoKinematicMPPICCloud
|
||||
|
||||
Description
|
||||
Cloud class to introduce Thermo kinematic MPPIC parcels
|
||||
|
||||
\*---------------------------------------------------------------------------*/
|
||||
|
||||
#ifndef basicThermoKinematicMPPICCloud_H
|
||||
#define basicThermoKinematicMPPICCloud_H
|
||||
|
||||
#include "Cloud.H"
|
||||
#include "KinematicCloud.H"
|
||||
#include "MPPICCloud.H"
|
||||
#include "basicThermoKinematicMPPICParcel.H"
|
||||
|
||||
#include "ThermoCloud.H"
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
namespace Foam
|
||||
{
|
||||
typedef MPPICCloud
|
||||
<
|
||||
ThermoCloud
|
||||
<
|
||||
KinematicCloud
|
||||
<
|
||||
Cloud
|
||||
<
|
||||
basicThermoKinematicMPPICParcel
|
||||
>
|
||||
>
|
||||
>
|
||||
>
|
||||
basicThermoKinematicMPPICCloud;
|
||||
}
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
#endif
|
||||
|
||||
// ************************************************************************* //
|
||||
|
|
@ -21,60 +21,43 @@ License
|
|||
You should have received a copy of the GNU General Public License
|
||||
along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
|
||||
|
||||
Application
|
||||
NOxFoam
|
||||
Class
|
||||
Foam::basicThermoKinematicMPPICParcel
|
||||
|
||||
Description
|
||||
Definition of basic Thermo kinematic MPPIC parcel
|
||||
|
||||
SourceFiles
|
||||
basicThermoKinematicMPPICParcel.H
|
||||
|
||||
\*---------------------------------------------------------------------------*/
|
||||
|
||||
#include "fvCFD.H"
|
||||
#include "fvOptions.H"
|
||||
#include "simpleControl.H"
|
||||
#ifndef basicThermoKinematicMPPICParcel_H
|
||||
#define basicThermoKinematicMPPICParcel_H
|
||||
|
||||
#include "contiguous.H"
|
||||
#include "particle.H"
|
||||
#include "KinematicParcel.H"
|
||||
#include "MPPICParcel.H"
|
||||
|
||||
#include "ThermoParcel.H"
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
int main(int argc, char *argv[])
|
||||
namespace Foam
|
||||
{
|
||||
#include "setRootCase.H"
|
||||
#include "createTime.H"
|
||||
#include "createMesh.H"
|
||||
#include "createControl.H"
|
||||
#include "createFields.H"
|
||||
#include "createFvOptions.H"
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
typedef MPPICParcel < ThermoParcel < KinematicParcel < particle > > > basicThermoKinematicMPPICParcel;
|
||||
|
||||
Info<< "\nStarting time loop\n" << endl;
|
||||
|
||||
if(instantaneousRadicals)
|
||||
template<>
|
||||
inline bool contiguous<basicThermoKinematicMPPICParcel>()
|
||||
{
|
||||
Info<< "Instantaneous mass fraction field will be used for O, OH radicals" << endl;
|
||||
return true;
|
||||
}
|
||||
else
|
||||
{
|
||||
Info<< "Partial Equilibrium Approach is selected for O, OH radicals" << endl;
|
||||
}
|
||||
|
||||
while (simple.loop())
|
||||
{
|
||||
Info<< "Time = " << runTime.timeName() << nl << endl;
|
||||
|
||||
#include "SourceCalc.H"
|
||||
#include "NH3Eqn.H"
|
||||
#include "NOEqn.H"
|
||||
|
||||
runTime.write();
|
||||
|
||||
Info<< "ExecutionTime = " << runTime.elapsedCpuTime() << " s"
|
||||
<< " ClockTime = " << runTime.elapsedClockTime() << " s"
|
||||
<< nl << endl;
|
||||
}
|
||||
|
||||
Info<< "End\n" << endl;
|
||||
|
||||
return 0;
|
||||
}
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
#endif
|
||||
|
||||
// ************************************************************************* //
|
||||
|
|
@ -0,0 +1,38 @@
|
|||
/*---------------------------------------------------------------------------*\
|
||||
========= |
|
||||
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
||||
\\ / O peration |
|
||||
\\ / A nd | Copyright (C) 2013 OpenFOAM Foundation
|
||||
\\/ M anipulation |
|
||||
-------------------------------------------------------------------------------
|
||||
License
|
||||
This file is part of OpenFOAM.
|
||||
|
||||
OpenFOAM is free software: you can redistribute it and/or modify it
|
||||
under the terms of the GNU General Public License as published by
|
||||
the Free Software Foundation, either version 3 of the License, or
|
||||
(at your option) any later version.
|
||||
|
||||
OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
|
||||
ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
|
||||
FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
|
||||
for more details.
|
||||
|
||||
You should have received a copy of the GNU General Public License
|
||||
along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
|
||||
|
||||
\*---------------------------------------------------------------------------*/
|
||||
|
||||
#include "basicThermoKinematicMPPICParcel.H"
|
||||
#include "Cloud.H"
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
namespace Foam
|
||||
{
|
||||
defineTemplateTypeNameAndDebug(basicThermoKinematicMPPICParcel, 0);
|
||||
defineTemplateTypeNameAndDebug(Cloud<basicThermoKinematicMPPICParcel>, 0);
|
||||
}
|
||||
|
||||
|
||||
// ************************************************************************* //
|
||||
|
|
@ -0,0 +1,67 @@
|
|||
/*---------------------------------------------------------------------------*\
|
||||
========= |
|
||||
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
||||
\\ / O peration |
|
||||
\\ / A nd | Copyright (C) 2013-2015 OpenFOAM Foundation
|
||||
\\/ M anipulation |
|
||||
-------------------------------------------------------------------------------
|
||||
License
|
||||
This file is part of OpenFOAM.
|
||||
|
||||
OpenFOAM is free software: you can redistribute it and/or modify it
|
||||
under the terms of the GNU General Public License as published by
|
||||
the Free Software Foundation, either version 3 of the License, or
|
||||
(at your option) any later version.
|
||||
|
||||
OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
|
||||
ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
|
||||
FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
|
||||
for more details.
|
||||
|
||||
You should have received a copy of the GNU General Public License
|
||||
along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
|
||||
|
||||
\*---------------------------------------------------------------------------*/
|
||||
|
||||
#include "basicThermoKinematicMPPICCloud.H"
|
||||
|
||||
#include "makeParcelCloudFunctionObjects.H"
|
||||
|
||||
// Kinematic sub-models
|
||||
#include "makeParcelForces.H"
|
||||
#include "makeParcelDispersionModels.H"
|
||||
#include "makeParcelInjectionModels.H"
|
||||
#include "makeParcelPatchInteractionModels.H"
|
||||
#include "makeParcelStochasticCollisionModels.H"
|
||||
#include "makeParcelSurfaceFilmModels.H"
|
||||
|
||||
// MPPIC sub-models
|
||||
#include "makeMPPICParcelDampingModels.H"
|
||||
#include "makeMPPICParcelIsotropyModels.H"
|
||||
#include "makeMPPICParcelPackingModels.H"
|
||||
|
||||
// Thermo sub-models
|
||||
#include "makeParcelHeatTransferModels.H"
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
makeParcelCloudFunctionObjects(basicThermoKinematicMPPICCloud);
|
||||
|
||||
// Kinematic sub-models
|
||||
makeParcelForces(basicThermoKinematicMPPICCloud);
|
||||
|
||||
makeParcelDispersionModels(basicThermoKinematicMPPICCloud);
|
||||
makeParcelInjectionModels(basicThermoKinematicMPPICCloud);
|
||||
makeParcelPatchInteractionModels(basicThermoKinematicMPPICCloud);
|
||||
makeParcelStochasticCollisionModels(basicThermoKinematicMPPICCloud);
|
||||
makeParcelSurfaceFilmModels(basicThermoKinematicMPPICCloud);
|
||||
|
||||
// MPPIC sub-models
|
||||
makeMPPICParcelDampingModels(basicThermoKinematicMPPICCloud);
|
||||
makeMPPICParcelIsotropyModels(basicThermoKinematicMPPICCloud);
|
||||
makeMPPICParcelPackingModels(basicThermoKinematicMPPICCloud);
|
||||
|
||||
// Thermo sub-models
|
||||
makeParcelHeatTransferModels(basicThermoKinematicMPPICCloud);
|
||||
|
||||
// ************************************************************************* //
|
||||
|
|
@ -86,4 +86,17 @@ Foam::basicChemistryModel::~basicChemistryModel()
|
|||
{}
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
|
||||
|
||||
void Foam::basicChemistryModel::edm
|
||||
(
|
||||
const scalarField& mixing,
|
||||
const scalar A,
|
||||
const scalar B
|
||||
)
|
||||
{
|
||||
NotImplemented;
|
||||
}
|
||||
|
||||
|
||||
// ************************************************************************* //
|
||||
|
|
|
|||
|
|
@ -159,6 +159,14 @@ public:
|
|||
//- Calculates the reaction rates
|
||||
virtual void calculate() = 0;
|
||||
|
||||
//- Calculates the eddy dissipation reaction rates
|
||||
virtual void edm
|
||||
(
|
||||
const scalarField& mixing,
|
||||
const scalar A,
|
||||
const scalar B
|
||||
);
|
||||
|
||||
//- Solve the reaction system for the given time step
|
||||
// and return the characteristic time
|
||||
virtual scalar solve(const scalar deltaT) = 0;
|
||||
|
|
|
|||
|
|
@ -756,6 +756,104 @@ void Foam::chemistryModel<CompType, ThermoType>::calculate()
|
|||
}
|
||||
|
||||
|
||||
template<class CompType, class ThermoType>
|
||||
void Foam::chemistryModel<CompType, ThermoType>::edm
|
||||
(
|
||||
const scalarField& mixing,
|
||||
const scalar A,
|
||||
const scalar B
|
||||
)
|
||||
{
|
||||
if (!this->chemistry_)
|
||||
{
|
||||
return;
|
||||
}
|
||||
|
||||
const volScalarField rho
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"rho",
|
||||
this->time().timeName(),
|
||||
this->mesh(),
|
||||
IOobject::NO_READ,
|
||||
IOobject::NO_WRITE,
|
||||
false
|
||||
),
|
||||
this->thermo().rho()
|
||||
);
|
||||
|
||||
forAll(rho, celli)
|
||||
{
|
||||
const scalar rhoi = rho[celli];
|
||||
const scalar mixingi = mixing[celli];
|
||||
|
||||
scalarField Yi(nSpecie_, 0.0);
|
||||
for (label k=0; k<nSpecie_; k++)
|
||||
{
|
||||
Yi[k] = Y_[k][celli];
|
||||
}
|
||||
|
||||
scalarField om(nSpecie_, 0.0);
|
||||
|
||||
scalar minYR = 1.0;
|
||||
scalar YP = 0.0;
|
||||
forAll(reactions_, m)
|
||||
{
|
||||
const Reaction<ThermoType>& R = reactions_[m];
|
||||
|
||||
scalar num = 0.0;
|
||||
scalar den = 0.0;
|
||||
|
||||
forAll(R.lhs(), s)
|
||||
{
|
||||
const label k = R.lhs()[s].index;
|
||||
const scalar vk = R.lhs()[s].stoichCoeff + SMALL;
|
||||
const scalar Wk = specieThermo_[k].W();
|
||||
|
||||
minYR = min(minYR, Yi[k]/Wk/vk);
|
||||
}
|
||||
|
||||
forAll(R.rhs(), s)
|
||||
{
|
||||
const label k = R.rhs()[s].index;
|
||||
const scalar vk = R.rhs()[s].stoichCoeff;
|
||||
const scalar Wk = specieThermo_[k].W();
|
||||
|
||||
num += Yi[k];
|
||||
den += vk*Wk;
|
||||
}
|
||||
YP = num/den;
|
||||
|
||||
scalar omegam = min (A*rhoi*mixingi*minYR, A*B*rhoi*mixingi*YP);
|
||||
|
||||
forAll(R.lhs(), s)
|
||||
{
|
||||
const label k = R.lhs()[s].index;
|
||||
const scalar vk = R.lhs()[s].stoichCoeff;
|
||||
const scalar Wk = specieThermo_[k].W();
|
||||
|
||||
om[k] -= vk*Wk*omegam;
|
||||
}
|
||||
|
||||
forAll(R.rhs(), s)
|
||||
{
|
||||
const label k = R.rhs()[s].index;
|
||||
const scalar vk = R.rhs()[s].stoichCoeff;
|
||||
const scalar Wk = specieThermo_[k].W();
|
||||
|
||||
om[k] += vk*Wk*omegam;
|
||||
}
|
||||
}
|
||||
|
||||
for (label k=0; k<nSpecie_; k++)
|
||||
{
|
||||
RR_[k][celli] = om[k];
|
||||
}
|
||||
}
|
||||
}
|
||||
|
||||
|
||||
template<class CompType, class ThermoType>
|
||||
template<class DeltaTType>
|
||||
Foam::scalar Foam::chemistryModel<CompType, ThermoType>::solve
|
||||
|
|
|
|||
|
|
@ -191,6 +191,14 @@ public:
|
|||
//- Calculates the reaction rates
|
||||
virtual void calculate();
|
||||
|
||||
//- Calculates the eddy dissipation reaction rates
|
||||
virtual void edm
|
||||
(
|
||||
const scalarField& mixing,
|
||||
const scalar A,
|
||||
const scalar B
|
||||
);
|
||||
|
||||
|
||||
// Chemistry model functions (overriding abstract functions in
|
||||
// basicChemistryModel.H)
|
||||
|
|
|
|||
|
|
@ -25,6 +25,7 @@ submodels/absorptionEmissionModel/binaryAbsorptionEmission/binaryAbsorptionEmiss
|
|||
submodels/absorptionEmissionModel/greyMeanAbsorptionEmission/greyMeanAbsorptionEmission.C
|
||||
submodels/absorptionEmissionModel/wideBandAbsorptionEmission/wideBandAbsorptionEmission.C
|
||||
submodels/absorptionEmissionModel/greyMeanSolidAbsorptionEmission/greyMeanSolidAbsorptionEmission.C
|
||||
submodels/absorptionEmissionModel/wsggmAbsorptionEmission/wsggmAbsorptionEmission.C
|
||||
|
||||
/* Soot model */
|
||||
submodels/sootModel/sootModel/sootModel.C
|
||||
|
|
|
|||
|
|
@ -0,0 +1,206 @@
|
|||
/*---------------------------------------------------------------------------*\
|
||||
========= |
|
||||
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
||||
\\ / O peration |
|
||||
\\ / A nd | Copyright (C) 2011-2013 OpenFOAM Foundation
|
||||
\\/ M anipulation |
|
||||
-------------------------------------------------------------------------------
|
||||
License
|
||||
This file is part of OpenFOAM.
|
||||
|
||||
OpenFOAM is free software: you can redistribute it and/or modify it
|
||||
under the terms of the GNU General Public License as published by
|
||||
the Free Software Foundation, either version 3 of the License, or
|
||||
(at your option) any later version.
|
||||
|
||||
OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
|
||||
ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
|
||||
FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
|
||||
for more details.
|
||||
|
||||
You should have received a copy of the GNU General Public License
|
||||
along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
|
||||
|
||||
\*---------------------------------------------------------------------------*/
|
||||
|
||||
#include "wsggmAbsorptionEmission.H"
|
||||
#include "addToRunTimeSelectionTable.H"
|
||||
#include "unitConversion.H"
|
||||
#include "zeroGradientFvPatchFields.H"
|
||||
#include "basicSpecieMixture.H"
|
||||
|
||||
// * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * //
|
||||
|
||||
namespace Foam
|
||||
{
|
||||
namespace radiation
|
||||
{
|
||||
defineTypeNameAndDebug(wsggmAbsorptionEmission, 0);
|
||||
|
||||
addToRunTimeSelectionTable
|
||||
(
|
||||
absorptionEmissionModel,
|
||||
wsggmAbsorptionEmission,
|
||||
dictionary
|
||||
);
|
||||
}
|
||||
}
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * * Constructors * * * * * * * * * * * * * * //
|
||||
|
||||
Foam::radiation::wsggmAbsorptionEmission::wsggmAbsorptionEmission
|
||||
(
|
||||
const dictionary& dict,
|
||||
const fvMesh& mesh
|
||||
)
|
||||
:
|
||||
absorptionEmissionModel(dict, mesh),
|
||||
coeffsDict_((dict.subDict(typeName + "Coeffs"))),
|
||||
thermo_(mesh.lookupObject<fluidThermo>("thermophysicalProperties")),
|
||||
emissivityCoeffs_(coeffsDict_.lookup("emissivityCoeffs")),
|
||||
fittingFactors_(coeffsDict_.lookup("fittingFactors")),
|
||||
pathLength_(coeffsDict_.lookup("pathLength"))
|
||||
{
|
||||
if (!isA<basicMultiComponentMixture>(thermo_))
|
||||
{
|
||||
FatalErrorIn
|
||||
(
|
||||
"radiation::wsggmAbsorptionEmission::wsggmAbsorptionEmission"
|
||||
"("
|
||||
"const dictionary&, "
|
||||
"const fvMesh&"
|
||||
")"
|
||||
) << "Model requires a multi-component thermo package"
|
||||
<< abort(FatalError);
|
||||
}
|
||||
}
|
||||
|
||||
// * * * * * * * * * * * * * * * * Destructor * * * * * * * * * * * * * * * //
|
||||
|
||||
Foam::radiation::wsggmAbsorptionEmission::~wsggmAbsorptionEmission()
|
||||
{}
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * //
|
||||
|
||||
Foam::tmp<Foam::volScalarField>
|
||||
Foam::radiation::wsggmAbsorptionEmission::aCont(const label bandI) const
|
||||
{
|
||||
|
||||
const basicSpecieMixture& mixture =
|
||||
dynamic_cast<const basicSpecieMixture&>(thermo_);
|
||||
|
||||
const volScalarField& T = thermo_.T();
|
||||
const volScalarField& p = thermo_.p();
|
||||
|
||||
tmp<volScalarField> ta
|
||||
(
|
||||
new volScalarField
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"aCont" + name(bandI),
|
||||
mesh().time().timeName(),
|
||||
mesh(),
|
||||
IOobject::NO_READ,
|
||||
IOobject::NO_WRITE
|
||||
),
|
||||
mesh(),
|
||||
dimensionedScalar("a", dimless/dimLength, 0.0),
|
||||
zeroGradientFvPatchVectorField::typeName
|
||||
)
|
||||
);
|
||||
|
||||
scalarField& a = ta.ref().primitiveFieldRef();
|
||||
|
||||
label indexCO2, indexH2O;
|
||||
scalar wCO2, wH2O;
|
||||
forAll(mixture.Y(), specieI)
|
||||
{
|
||||
if(mixture.Y(specieI).name() == "CO2")
|
||||
{
|
||||
indexCO2 = specieI;
|
||||
wCO2 = mixture.W(indexCO2);
|
||||
}
|
||||
else if(mixture.Y(specieI).name() == "H2O")
|
||||
{
|
||||
indexH2O = specieI;
|
||||
wH2O = mixture.W(indexH2O);
|
||||
}
|
||||
}
|
||||
|
||||
forAll(a,celli)
|
||||
{
|
||||
scalar invWt = 0.0;
|
||||
forAll(mixture.Y(), specieI)
|
||||
{
|
||||
invWt += mixture.Y(specieI)[celli]/mixture.W(specieI);
|
||||
}
|
||||
scalar meanWt = 1.0/invWt; //Mean molecular weight [kg/kmole]
|
||||
|
||||
// Sum of partial pressures of all absorbing gases [atm]
|
||||
scalar Pcell_atm = paToAtm(meanWt*p[celli]*(mixture.Y(indexCO2)[celli]/wCO2 + mixture.Y(indexH2O)[celli]/wH2O));
|
||||
|
||||
// P*S at Eq.(1) of Ref. [atm*m]
|
||||
scalar presPathLength = Pcell_atm*pathLength_.value();
|
||||
|
||||
// Limit cell temperature [K]
|
||||
scalar Tcell = min(T[celli],2400.0);
|
||||
|
||||
scalar emissivity = 0.0;
|
||||
forAll(emissivityCoeffs_, i) //k_i at Eq.(1) of Ref.
|
||||
{
|
||||
scalar weightingFactor = 0.0;
|
||||
for(label j=0 ; j<4 ; j++)
|
||||
{
|
||||
//a_epsilon_i at Eq.(3) of Ref.
|
||||
weightingFactor += fittingFactors_[i][j]*pow(Tcell, j);
|
||||
}
|
||||
//epsilon at Eq.(1) of Ref.
|
||||
emissivity += weightingFactor * (1.0 - exp(-1.0*emissivityCoeffs_[i] * presPathLength));
|
||||
}
|
||||
emissivity = min(emissivity, 0.99999);
|
||||
|
||||
//Fluent theory guide(ver. 13) Eq.(5-106) [1/m]
|
||||
a[celli] = max(1e-2, (-1.0)*log(1.0-emissivity)/pathLength_.value());
|
||||
}
|
||||
|
||||
ta.ref().correctBoundaryConditions();
|
||||
|
||||
return ta;
|
||||
}
|
||||
|
||||
|
||||
Foam::tmp<Foam::volScalarField>
|
||||
Foam::radiation::wsggmAbsorptionEmission::eCont(const label bandI) const
|
||||
{
|
||||
return aCont(bandI);
|
||||
}
|
||||
|
||||
|
||||
Foam::tmp<Foam::volScalarField>
|
||||
Foam::radiation::wsggmAbsorptionEmission::ECont(const label bandI) const
|
||||
{
|
||||
tmp<volScalarField> E
|
||||
(
|
||||
new volScalarField
|
||||
(
|
||||
IOobject
|
||||
(
|
||||
"ECont" + name(bandI),
|
||||
mesh_.time().timeName(),
|
||||
mesh_,
|
||||
IOobject::NO_READ,
|
||||
IOobject::NO_WRITE
|
||||
),
|
||||
mesh_,
|
||||
dimensionedScalar("E", dimMass/dimLength/pow3(dimTime), 0.0)
|
||||
)
|
||||
);
|
||||
|
||||
return E;
|
||||
}
|
||||
|
||||
|
||||
// ************************************************************************* //
|
||||
|
|
@ -0,0 +1,163 @@
|
|||
/*---------------------------------------------------------------------------*\
|
||||
========= |
|
||||
\\ / F ield | OpenFOAM: The Open Source CFD Toolbox
|
||||
\\ / O peration |
|
||||
\\ / A nd | Copyright (C) 2011-2012 OpenFOAM Foundation
|
||||
\\/ M anipulation |
|
||||
-------------------------------------------------------------------------------
|
||||
License
|
||||
This file is part of OpenFOAM.
|
||||
|
||||
OpenFOAM is free software: you can redistribute it and/or modify it
|
||||
under the terms of the GNU General Public License as published by
|
||||
the Free Software Foundation, either version 3 of the License, or
|
||||
(at your option) any later version.
|
||||
|
||||
OpenFOAM is distributed in the hope that it will be useful, but WITHOUT
|
||||
ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
|
||||
FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
|
||||
for more details.
|
||||
|
||||
You should have received a copy of the GNU General Public License
|
||||
along with OpenFOAM. If not, see <http://www.gnu.org/licenses/>.
|
||||
|
||||
Class
|
||||
Foam::radiation::wsggmAbsorptionEmission
|
||||
|
||||
Description
|
||||
wsggmAbsorptionEmission radiation absorption and emission coefficients
|
||||
for continuous phase
|
||||
|
||||
Ref. T.F Smith, Z.F. Shen and J.N. Friedman
|
||||
Evaluation of Coefficients for the Weighted Sum of Gray Gases Model,
|
||||
Transactions of the ASME, Vol. 104, pp. 602-608
|
||||
|
||||
Valid when Ptotal = 1atm, 0.001 < P*S < 10.0atm.m, 600 < T < 2400K
|
||||
(where P*S = sum of partial pressures of absorbing gases
|
||||
* mean beam length)
|
||||
|
||||
i.e. dictionary
|
||||
\verbatim
|
||||
wsggmAbsorptionEmissionCoeffs
|
||||
{
|
||||
pathLength pathLength [ 0 1 0 0 0 0 0] 0.251; //3.6V/S
|
||||
// Pw/Pc = 2.0
|
||||
emissivityCoeffs 3(0.4201 6.516 131.9);
|
||||
fittingFactors
|
||||
3
|
||||
(
|
||||
4(6.508e-1 -5.551e-4 3.029e-7 -5.353e-11)
|
||||
4(-0.2504e-1 6.112e-4 -3.882e-7 6.528e-11)
|
||||
4(2.718e-1 -3.118e-4 1.221e-7 -1.612e-11)
|
||||
);
|
||||
}
|
||||
\endverbatim
|
||||
|
||||
|
||||
SourceFiles
|
||||
wsggmAbsorptionEmission.C
|
||||
|
||||
\*---------------------------------------------------------------------------*/
|
||||
|
||||
#ifndef wsggmAbsorptionEmission_H
|
||||
#define wsggmAbsorptionEmission_H
|
||||
|
||||
#include "interpolationLookUpTable.H"
|
||||
#include "absorptionEmissionModel.H"
|
||||
#include "HashTable.H"
|
||||
#include "absorptionCoeffs.H"
|
||||
#include "fluidThermo.H"
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
namespace Foam
|
||||
{
|
||||
namespace radiation
|
||||
{
|
||||
|
||||
/*---------------------------------------------------------------------------*\
|
||||
Class wsggmAbsorptionEmission Declaration
|
||||
\*---------------------------------------------------------------------------*/
|
||||
|
||||
class wsggmAbsorptionEmission
|
||||
:
|
||||
public absorptionEmissionModel
|
||||
{
|
||||
public:
|
||||
|
||||
private:
|
||||
|
||||
// Private data
|
||||
|
||||
//- Absorption model dictionary
|
||||
dictionary coeffsDict_;
|
||||
|
||||
//- SLG thermo package
|
||||
const fluidThermo& thermo_;
|
||||
|
||||
//- absorption coefficitent for individual grey gases
|
||||
scalarList emissivityCoeffs_;
|
||||
|
||||
//- fitting Factors
|
||||
scalarListList fittingFactors_;
|
||||
|
||||
//characteristic path length [m] from dict
|
||||
dimensionedScalar pathLength_;
|
||||
|
||||
public:
|
||||
|
||||
//- Runtime type information
|
||||
TypeName("wsggmAbsorptionEmission");
|
||||
|
||||
|
||||
// Constructors
|
||||
|
||||
//- Construct from components
|
||||
wsggmAbsorptionEmission(const dictionary& dict, const fvMesh& mesh);
|
||||
|
||||
|
||||
//- Destructor
|
||||
virtual ~wsggmAbsorptionEmission();
|
||||
|
||||
|
||||
// Member Functions
|
||||
|
||||
// Access
|
||||
|
||||
// Absorption coefficient
|
||||
|
||||
//- Absorption coefficient for continuous phase
|
||||
tmp<volScalarField> aCont(const label bandI = 0) const;
|
||||
|
||||
|
||||
// Emission coefficient
|
||||
|
||||
//- Emission coefficient for continuous phase
|
||||
tmp<volScalarField> eCont(const label bandI = 0) const;
|
||||
|
||||
|
||||
// Emission contribution
|
||||
|
||||
//- Emission contribution for continuous phase
|
||||
tmp<volScalarField> ECont(const label bandI = 0) const;
|
||||
|
||||
|
||||
// Member Functions
|
||||
|
||||
inline bool isGrey() const
|
||||
{
|
||||
return true;
|
||||
}
|
||||
};
|
||||
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
} // End namespace radiation
|
||||
} // End namespace Foam
|
||||
|
||||
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
|
||||
|
||||
#endif
|
||||
|
||||
// ************************************************************************* //
|
||||
File diff suppressed because it is too large
Load diff
File diff suppressed because it is too large
Load diff
File diff suppressed because it is too large
Load diff
File diff suppressed because it is too large
Load diff
File diff suppressed because it is too large
Load diff
File diff suppressed because it is too large
Load diff
File diff suppressed because it is too large
Load diff
File diff suppressed because it is too large
Load diff
File diff suppressed because it is too large
Load diff
File diff suppressed because it is too large
Load diff
|
|
@ -1,25 +0,0 @@
|
|||
/*--------------------------------*- C++ -*----------------------------------*\
|
||||
| ========= | |
|
||||
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
|
||||
| \\ / O peration | Version: 4.x |
|
||||
| \\ / A nd | Web: www.OpenFOAM.org |
|
||||
| \\/ M anipulation | |
|
||||
\*---------------------------------------------------------------------------*/
|
||||
FoamFile
|
||||
{
|
||||
version 2;
|
||||
format ascii;
|
||||
class dictionary;
|
||||
location "constant";
|
||||
object SOxProperties;
|
||||
}
|
||||
|
||||
fuelSourceNames ( reactingCloud1:rhoTrans_C7H16 );
|
||||
|
||||
nitrogenMassInFuels ( 0.01 );
|
||||
|
||||
instantaneousRadicals false;
|
||||
|
||||
calculateFuelNOx false;
|
||||
|
||||
// ************************************************************************* //
|
||||
Some files were not shown because too many files have changed in this diff Show more
Loading…
Add table
Reference in a new issue