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

Stockmayer collision integral fit coefficients made constexpr

Browse source
This commit is contained in:
ignis 2018-06-18 16:16:03 +09:00
parent d71978e286
commit 6e4e4a96d6
1 changed files with 4 additions and 4 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

8
diffusivityModel/Stockmayer/StockmayerI.H
View file

@ -39,28 +39,28 @@ License
inline Foam::scalar Foam::Stockmayer::Omega11(const scalar Tstar)
{
const scalar a[6] = { 1.0548, 0.15504, 0.55909, 2.1705, 0.093193, 1.5 };
constexpr scalar a[6] = { 1.0548, 0.15504, 0.55909, 2.1705, 0.093193, 1.5 };
return a[0] * power(Tstar, -a[1]) + power(Tstar + a[2], -a[3]);
}
inline Foam::scalar Foam::Stockmayer::Omega22(const scalar Tstar)
{
const scalar b[6] = { 1.0413, 0.11930, 0.43628, 1.6041, 0.095661, 2.0 };
constexpr scalar b[6] = { 1.0413, 0.11930, 0.43628, 1.6041, 0.095661, 2.0 };
return b[0] * power(Tstar, -b[1]) + power(Tstar + b[2], -b[3]);
}
inline Foam::scalar Foam::Stockmayer::f11(const scalar Tstar)
{
const scalar a[6] = { 1.0548, 0.15504, 0.55909, 2.1705, 0.093193, 1.5 };
constexpr scalar a[6] = { 1.0548, 0.15504, 0.55909, 2.1705, 0.093193, 1.5 };
return 1.0 + sqr(deltaij_)*(expon(a[4]/Tstar)-expon(-a[5]/Tstar))/(2 + 2.5*deltaij_);
}
inline Foam::scalar Foam::Stockmayer::f22(const scalar Tstar)
{
const scalar b[6] = { 1.0413, 0.11930, 0.43628, 1.6041, 0.095661, 2.0 };
constexpr scalar b[6] = { 1.0413, 0.11930, 0.43628, 1.6041, 0.095661, 2.0 };
return 1.0 + sqr(deltaij_)*(expon(b[4]/Tstar)-expon(-b[5]/Tstar))/(2 + 2.5*deltaij_);
}