[jackjiang1989@gmail.com]

Hi All,
I'm trying to understand how viscosity is calculated in chtMultiRegionFoam with transport as polynomial.

I saw in the chtMultiRegionFoam - createfluidfield.h, it's using rhoReactionThermo, it's initialized as follows,
thermoFluid.set(i, rhoReactionThermo::New(fluidRegions[i]).ptr())

Then I'm looking into rhoReactionThermo, so it's from
rhoReactionThermo - rhoThermo - fluidThermo - basicThermo

There is a function tmp<volScalarField> mu() defined in fluidThermo
And there is some initial operation defined in basicThermo

But I can't find there is any relationship between "xxxThermo" and "polynomialTransport"
In polynomialTransportl.H, mu() is defined with polynomial function. And one more thing strange is that the return type is a scalar?? instead of a volScalarField?)

Anyone got any idea how/where this polynomialTtransport is used to calculate the viscosity in chtMultiRegionFoam?
I think it should be called like
Foam:olynomialTransport<xxxThermo, somepolynumber> mypolynomialtransport(xxxThermo, mucoeff, kappacoeff)
mypolynomialtransport.mu(p, T)

Thanks!