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

diffusivityModel mobility method

Browse source
This commit is contained in:
Yeongdo Park 2019-02-16 10:05:10 -05:00
parent 96011184a3
commit ad9cb13d57
2 changed files with 21 additions and 0 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

18
diffusivityModel/diffusivityModel/diffusivityModel.C
View file

@ -811,6 +811,24 @@ Foam::tmp<Foam::scalarField> Foam::diffusivityModel::kpure(const scalar T, const
}
Foam::tmp<Foam::volScalarField> Foam::diffusivityModel::mu(const label i, const volScalarField& T) const
{
const speciesTable &species_(thermo_.composition().species());
const word namei(species_[i]);
const label zi = thermo_.composition().z(i);
// - Elementary charge (default in [C])
const dimensionedScalar eCharge = constant::electromagnetic::e;
//- Boltzmann constant
const dimensionedScalar kB = constant::physicoChemical::k;
return (this->D(i)/T * (eCharge*zi/kB));
}
// * * * * * * * * * * * * * * Member Operators * * * * * * * * * * * * * * //
void Foam::diffusivityModel::operator=(const diffusivityModel& rhs)

3
diffusivityModel/diffusivityModel/diffusivityModel.H
View file

@ -225,6 +225,9 @@ public:
// Edit
tmp<volScalarField> mu(const label i, const volScalarField& T) const;
// Write