Skip to content
Permalink
Browse files

chtMultiRegionFoam: Removed Qdot from top-level solver

Qdot is only relevant for reacting cases, and even then it is only
written out for post-processing purposes. It has been removed from
chtMultiRegionFoam as this solvers' typical usage is for non-reacting
simulations.

If it is desired to obtain Qdot from a chtMultiRegionFoam simulation,
then a better way would be to implement a function object to look up the
reaction model and write it out.
  • Loading branch information...
Will Bainbridge
Will Bainbridge committed Oct 30, 2018
1 parent b7c6499 commit b59c71a15e7cd790e83da6e1b84be3c094bba300
@@ -19,7 +19,7 @@
==
rho*(U&g)
+ rad.Sh(thermo, he)
+ Qdot
+ reaction.Qdot()
+ fvOptions(rho, he)
);

@@ -16,13 +16,21 @@ if (Y.size())

{
reaction.correct();
Qdot = reaction.Qdot();
volScalarField Yt
(
IOobject("Yt", runTime.timeName(), mesh),
mesh,
dimensionedScalar("Yt", dimless, 0)
);

tmp<volScalarField> Yt(nullptr);

if (Y.size())
{
Yt = tmp<volScalarField>
(
new volScalarField
(
IOobject("Yt", runTime.timeName(), mesh),
mesh,
dimensionedScalar("Yt", dimless, 0)
)
);
}

forAll(Y, i)
{
@@ -49,7 +57,7 @@ if (Y.size())
fvOptions.correct(Yi);

Yi.max(0.0);
Yt += Yi;
Yt.ref() += Yi;
}
}

@@ -15,7 +15,6 @@ PtrList<volScalarField> KFluid(fluidRegions.size());
PtrList<volScalarField> dpdtFluid(fluidRegions.size());
PtrList<multivariateSurfaceInterpolationScheme<scalar>::fieldTable>
fieldsFluid(fluidRegions.size());
PtrList<volScalarField> QdotFluid(fluidRegions.size());

List<scalar> initialMassFluid(fluidRegions.size());

@@ -245,25 +244,6 @@ forAll(fluidRegions, i)
}
fieldsFluid[i].add(thermoFluid[i].he());

Info<< " Adding to QdotFluid\n" << endl;
QdotFluid.set
(
i,
new volScalarField
(
IOobject
(
"Qdot",
runTime.timeName(),
fluidRegions[i],
IOobject::READ_IF_PRESENT,
IOobject::AUTO_WRITE
),
fluidRegions[i],
dimensionedScalar("Qdot", dimEnergy/dimVolume/dimTime, 0.0)
)
);

Info<< " Adding MRF\n" << endl;
MRFfluid.set
(
@@ -43,8 +43,6 @@
multivariateSurfaceInterpolationScheme<scalar>::fieldTable& fields =
fieldsFluid[i];

volScalarField& Qdot = QdotFluid[i];

radiation::radiationModel& rad = radiation[i];

IOMRFZoneList& MRF = MRFfluid[i];

0 comments on commit b59c71a

Please sign in to comment.
You can’t perform that action at this time.