Merge branch 'develop' of develop.openfoam.com:Development/OpenFOAM-plus into develop

This commit is contained in:
mattijs
2017-06-28 15:07:04 +01:00
3 changed files with 7 additions and 9 deletions

View File

@ -101,7 +101,7 @@ public:
{ {
if (is.good() && !is.eof()) if (is.good() && !is.eof())
{ {
this->ddtPhiCoeff_ = readScalar(is); is >> this->ddtPhiCoeff_;
} }
// Ensure the old-old-time cell volumes are available // Ensure the old-old-time cell volumes are available

View File

@ -171,10 +171,11 @@ tmp<surfaceScalarField> ddtScheme<Type>::fvcDdtPhiCoeff
) )
); );
surfaceScalarField& ddtCouplingCoeff = tddtCouplingCoeff.ref();
if (ddtPhiCoeff_ < 0) if (ddtPhiCoeff_ < 0)
{ {
tddtCouplingCoeff.ref() =- min ddtCouplingCoeff -= min
( (
mag(phiCorr) mag(phiCorr)
/(mag(phi) + dimensionedScalar("small", phi.dimensions(), SMALL)), /(mag(phi) + dimensionedScalar("small", phi.dimensions(), SMALL)),
@ -183,14 +184,11 @@ tmp<surfaceScalarField> ddtScheme<Type>::fvcDdtPhiCoeff
} }
else else
{ {
tddtCouplingCoeff.ref() = ddtCouplingCoeff =
dimensionedScalar("ddtPhiCoeff", dimless, ddtPhiCoeff_); dimensionedScalar("ddtPhiCoeff", dimless, ddtPhiCoeff_);
} }
surfaceScalarField& ddtCouplingCoeff = tddtCouplingCoeff.ref(); surfaceScalarField::Boundary& ccbf = ddtCouplingCoeff.boundaryFieldRef();
surfaceScalarField::Boundary& ccbf =
ddtCouplingCoeff.boundaryFieldRef();
forAll(U.boundaryField(), patchi) forAll(U.boundaryField(), patchi)
{ {

View File

@ -110,12 +110,12 @@ calculateSpeciesRR
if (RR[0] > 0.0) if (RR[0] > 0.0)
{ {
production_[speciei][reactioni] = RR[0]; production_[speciei][reactioni] = RR[0];
productionInt_[speciei][reactioni] =+ dt*RR[0]; productionInt_[speciei][reactioni] += dt*RR[0];
} }
else if (RR[0] < 0.0) else if (RR[0] < 0.0)
{ {
consumption_[speciei][reactioni] = RR[0]; consumption_[speciei][reactioni] = RR[0];
consumptionInt_[speciei][reactioni] =+ dt*RR[0]; consumptionInt_[speciei][reactioni] += dt*RR[0];
} }
else else
{ {