MultiComponentPhaseModel: Fixed missing density time-derivative

This commit is contained in:
Will Bainbridge
2018-05-25 12:30:38 +01:00
parent 34d2bf7af8
commit 93c1b5644f

View File

@ -149,7 +149,7 @@ Foam::MultiComponentPhaseModel<BasePhaseModel>::YiEqn(volScalarField& Yi)
{ {
const volScalarField& alpha = *this; const volScalarField& alpha = *this;
const surfaceScalarField alphaRhoPhi(this->alphaRhoPhi()); const surfaceScalarField alphaRhoPhi(this->alphaRhoPhi());
const volScalarField rho(this->rho()); const volScalarField& rho = this->thermo().rho();
return return
( (