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

refactoring thermal conductivity calculation

Browse source
This commit is contained in:
Yeongdo Park 2018-12-17 01:27:55 -05:00
parent 36163004be
commit f8cdf0c875
3 changed files with 52 additions and 14 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

26
diffusivityModel/Particle/ParticleI.H
View file

@ -83,15 +83,14 @@ inline Foam::scalar Foam::Particle::Zrot()
inline Foam::scalar Foam::Particle::zeta(const scalar T, const scalar wellDepth)
const
{
const scalar a3 = sqrt(pow3(pi));
const scalar a1 = a3/2.0;
const scalar pi32 = sqrt(pow3(pi));
const scalar a2 = sqr(pi) / 4.0 + 2.0;
const scalar t1 = sqrt(wellDepth/T);
const scalar t2 = wellDepth/T;
const scalar t3 = t1 * t2;
const scalar rSqrtTstar = 1.0/sqrt(T/wellDepth);
const scalar rTstar = wellDepth/T;
return 1.0 + a1 * t1 + a2 * t2 + a3 * t3;
return 1.0 + pi32 * rSqrtTstar * (0.5 + rTstar) + a2 * rTstar;
}
inline Foam::scalar Foam::Particle::R()
@ -125,7 +124,7 @@ inline Foam::scalar Foam::Particle::CvVib(const scalar p, const scalar T)
inline Foam::scalar Foam::Particle::Cv(const scalar p, const scalar T)
const
{
return thermo_.Cv(p,T) * W();
return thermo_.cv(p,T);
}
@ -147,15 +146,14 @@ inline Foam::scalar Foam::Particle::fTrans(const scalar rSc)
inline Foam::scalar Foam::Particle::lambda(const scalar mu, const scalar D, const scalar Zrot, const GasState &state)
const
{
const scalar rSc = state.rho() * D / mu;
const scalar fVib = state.rho() * D / mu;
const scalar A = 5./2. - rSc;
const scalar B = Zrot + (2./pi) * ((5./3.)*CvRot()/R() + rSc);
const scalar AB = (2./pi) * (A/B);
const scalar A = 5./2. - fVib;
const scalar B = Zrot + (2./pi) * ((5./3.)*CvRot()/R() + fVib);
const scalar c1 = (2./pi) * (A/B);
const scalar fTrans = (5./2.) * (1.0 - AB * CvRot() / CvTrans());
const scalar fRot = rSc * (1.0 + AB);
const scalar fVib = rSc;
const scalar fTrans = (5./2.) * (1.0 - c1 * CvRot() / CvTrans());
const scalar fRot = fVib * (1.0 + c1);
return (mu/W()) * (fTrans*CvTrans() + fRot*CvRot() + fVib*CvVib(state.p(), state.T()));
}

38
diffusivityModel/diffusivityModel/diffusivityModel.C
View file

@ -714,6 +714,44 @@ void Foam::diffusivityModel::correct()
}
Foam::tmp<Foam::scalarField> Foam::diffusivityModel::kpure(const scalar T, const scalar p, const UList<scalar> &Y)
{
const speciesTable &species_(thermo_.composition().species());
scalarSymmetricSquareMatrix Dij(species_.size());
scalarField localY(Y);
scalarField Di(species_.size());
scalarField Dii(species_.size());
scalarField muI(species_.size());
scalarField *pkI = new scalarField(species_.size());
scalarField &kI = *pkI;
{
const scalar pi = p;
const scalar Ti = T;
GasState state (pi, Ti, localY);
const scalarField &localX = state.X();
calculateMuD (muI, Dij, state);
forAll (Dii, i)
{
Dii[i] = Dij(i,i);
}
// Pure Thermal conductivities
calculateK ( kI, muI, Dii, state);
}
return tmp<scalarField>(pkI);
}
// * * * * * * * * * * * * * * Member Operators * * * * * * * * * * * * * * //
void Foam::diffusivityModel::operator=(const diffusivityModel& rhs)

2
diffusivityModel/diffusivityModel/diffusivityModel.H
View file

@ -217,6 +217,8 @@ public:
// Check
tmp<scalarField> kpure(const scalar T, const scalar p, const UList<scalar> &Y);
// Edit
// Write