# how transport polynomial is called?

 Register Blogs Members List Search Today's Posts Mark Forums Read

January 9, 2021, 19:10
how transport polynomial is called?
#1
New Member

jackjiang
Join Date: Feb 2020
Posts: 2
Rep Power: 0
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!
Attached Images
 thermophysicalProperties.PNG (66.5 KB, 3 views) createfluidfield.PNG (164.5 KB, 3 views)