research/OpenFOAM-2.4.x
1
0
Fork 0
Code Issues Pull requests Projects Releases Packages Wiki Activity Actions

chemistryModel electron reactions are selected with indices and Te update method changed

Browse source
This commit is contained in:
ignis 2016-07-03 14:07:58 +09:00
parent 5a50a80378
commit 765926435b
1 changed files with 5 additions and 9 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

14
src/thermophysicalModels/chemistryModel/chemistryModel/chemistryModel/chemistryModel.C
View file

@ -832,15 +832,15 @@ Foam::scalar Foam::chemistryModel<CompType, ThermoType>::solve
const scalarField& Te = this->mesh().template lookupObject<volScalarField>("Te").internalField();
const scalarField& ne = this->mesh().template lookupObject<volScalarField>("ne").internalField();
PtrList<Reaction<ThermoType> > eReactions_(nReaction_);
label neReactions_ = 0;
labelList eRxnLabels;
for (label i=0; i<nReaction_; i++)
{
if (reactions_[i].species().contains("E-"))
{
// Info << "With E- " << reactions_[i] << endl;
eReactions_.set(neReactions_, const_cast<Reaction<ThermoType>*>(reactions_(i)));
eRxnLabels.append(i);
neReactions_++;
}
else
@ -849,11 +849,6 @@ Foam::scalar Foam::chemistryModel<CompType, ThermoType>::solve
}
}
for (label i=0; i<neReactions_; i++)
{
Info << eReactions_[i] << endl;
}
forAll(rho, celli)
{
const scalar rhoi = rho[celli];
@ -867,8 +862,9 @@ Foam::scalar Foam::chemistryModel<CompType, ThermoType>::solve
for (label i=0; i<neReactions_; i++)
{
// Info << eReactions_[i] << endl;
eReactions_[i].Te() = Tei;
label ei = eRxnLabels[i];
Info << reactions_[ei] << endl;
reactions_[ei].Te(Tei);
}
for (label i=0; i<nSpecie_; i++)