zhyang-dev/OpenFOAM-12
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

thermophysicalTransportModels::FickianEddyDiffusivity: Updated diffusivity coefficients to Function2

Browse Source 
The specie and Soret diffusivity coefficients are now functions of pressure and
temperature utilising the new run-time selectable Function2.
This commit is contained in:
Henry Weller
2020-11-12 18:10:30 +00:00
parent 2cd197a536
commit 922c172352
2 changed files with 16 additions and 8 deletions
Download Patch File Download Diff File Expand all files Collapse all files

18
src/ThermophysicalTransportModels/turbulence/FickianEddyDiffusivity/FickianEddyDiffusivity.C
View File

@ -29,7 +29,7 @@ License
#include "fvcSnGrad.H"
#include "fvmSup.H"
#include "surfaceInterpolate.H"
#include "Function1Evaluate.H"
#include "Function2Evaluate.H"
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
@ -97,7 +97,7 @@ FickianEddyDiffusivity<TurbulenceThermophysicalTransportModel>::read()
forAll(species, i)
{
D_.set(i, Function1<scalar>::New(species[i], Ddict).ptr());
D_.set(i, Function2<scalar>::New(species[i], Ddict).ptr());
}
if (this->coeffDict_.found("DT"))
@ -106,7 +106,7 @@ FickianEddyDiffusivity<TurbulenceThermophysicalTransportModel>::read()
forAll(species, i)
{
DT_.set(i, Function1<scalar>::New(species[i], DTdict).ptr());
DT_.set(i, Function2<scalar>::New(species[i], DTdict).ptr());
}
}
@ -133,7 +133,13 @@ FickianEddyDiffusivity<TurbulenceThermophysicalTransportModel>::DEff
(
"DEff",
this->momentumTransport().rho()
*evaluate(D_[composition.index(Yi)], dimViscosity, this->thermo().T())
*evaluate
(
D_[composition.index(Yi)],
dimViscosity,
this->thermo().p(),
this->thermo().T()
)
+ (this->Prt_/Sct_)*this->alphat()
);
}
@ -154,6 +160,7 @@ FickianEddyDiffusivity<TurbulenceThermophysicalTransportModel>::DEff
this->momentumTransport().rho().boundaryField()[patchi]
*D_[composition.index(Yi)].value
(
this->thermo().p().boundaryField()[patchi],
this->thermo().T().boundaryField()[patchi]
)
+ this->Prt_.value()/Sct_.value()*this->alphat(patchi);
@ -170,6 +177,7 @@ FickianEddyDiffusivity<TurbulenceThermophysicalTransportModel>::j
if (DT_.size())
{
const basicSpecieMixture& composition = this->thermo().composition();
const volScalarField& p = this->thermo().T();
const volScalarField& T = this->thermo().T();
return
@ -177,7 +185,7 @@ FickianEddyDiffusivity<TurbulenceThermophysicalTransportModel>::j
j(Yi)
- fvc::interpolate
(
evaluate(DT_[composition.index(Yi)], dimDynamicViscosity, T)
evaluate(DT_[composition.index(Yi)], dimDynamicViscosity, p, T)
)
*fvc::snGrad(T)/fvc::interpolate(T);
}

6
src/ThermophysicalTransportModels/turbulence/FickianEddyDiffusivity/FickianEddyDiffusivity.H
View File

@ -54,7 +54,7 @@ SourceFiles
\*---------------------------------------------------------------------------*/
#include "unityLewisEddyDiffusivity.H"
#include "Function1.H"
#include "Function2.H"
#ifndef FickianEddyDiffusivity_H
#define FickianEddyDiffusivity_H
@ -87,11 +87,11 @@ protected:
//- List of specie mass diffusion coefficient functions
// w.r.t the mixture [m^2/s]
PtrList<Function1<scalar>> D_;
PtrList<Function2<scalar>> D_;
//- List of specie Soret thermal diffusion coefficient
// functions [kg/m/s]
PtrList<Function1<scalar>> DT_;
PtrList<Function2<scalar>> DT_;
public: