changed molecular weight

This commit is contained in:
ekinaci
2017-02-24 09:41:21 +01:00
parent 44d111345d
commit cf3a6a5cbe
2 changed files with 6 additions and 5 deletions

View File

@ -50,8 +50,11 @@ tmp<fv::convectionScheme<scalar> > mvConvection
}
}
// calculate total molar number in volume
N = Yt*rhoeps/W;
Info << "molecular weight of mixture" << W << endl;
Y[inertIndex] = scalar(1) - Yt;
Y[inertIndex].max(0.0);
}

View File

@ -18,9 +18,8 @@
basicSpecieMixture& composition = thermo.composition();
PtrList<volScalarField>& Y = composition.Y();
// The molecular weight [kj/mol] ====== EKI
//const volScalarField W(composition.W());
const volScalarField& W = composition.W();
// The molecular weight [kj/mol]
const volScalarField W(composition.W());
const word inertSpecie(thermo.lookup("inertSpecie"));
@ -42,7 +41,7 @@
}
fields.add(thermo.he());
// define total number of moles in volume ====== EKI
// Total Mole in Volume
volScalarField N
(
IOobject
@ -55,7 +54,6 @@
),
mesh
);
// =============================================
volScalarField rho
(