CFD Online Discussion Forums - Thermophysical transport model to include shear rate

Greetings!

I'm currently working on new thermophysical model with its own equation of state and transport properties. I get stuck when I try to implement shear rate dependence of viscosity. Is this feasible? Where would you calculate shear rate and pass it to viscosity function? Goal is to capture shear thinning behaviour of polymers with multiphase flow.

This is where ("/src/thermophysicalModels/specie/transport/crossWLF/crossWLFTransportI.H") viscosity is evaluated. It takes p and T.

Code:

template<class Thermo>

inline Foam::scalar Foam::crossWLFTransport<Thermo>::mu

(

    const scalar p,

    const scalar T

) const

{

    return D1_*::exp(-(A1_*(T-D2_))/(A2_+T-D2_));

}

This is file /src/thermophysicalModels/basic/rhoThermo/heRhoThermo.H, where I suppose mu is called.

Code:

template<class BasicPsiThermo, class MixtureType>

void Foam::heRhoThermo<BasicPsiThermo, MixtureType>::calculate()

{

    const scalarField& hCells = this->he().internalField();

    const scalarField& pCells = this->p_.internalField();

    //const vectorField& UCells = this->U_.internalField(); ofcourse it doesn't work like this



    scalarField& TCells = this->T_.internalField();

    scalarField& psiCells = this->psi_.internalField();

    scalarField& rhoCells = this->rho_.internalField();

    scalarField& muCells = this->mu_.internalField();

    scalarField& alphaCells = this->alpha_.internalField();



    forAll(TCells, celli)

    {

        const typename MixtureType::thermoType& mixture_ =

            this->cellMixture(celli);



        TCells[celli] = mixture_.THE

        (

            hCells[celli],

            pCells[celli],

            TCells[celli]

        );



        psiCells[celli] = mixture_.psi(pCells[celli], TCells[celli]);

        rhoCells[celli] = mixture_.rho(pCells[celli], TCells[celli]);



        muCells[celli] = mixture_.mu(pCells[celli], TCells[celli]); // is this the line where mu is called?

        alphaCells[celli] = mixture_.alphah(pCells[celli], TCells[celli]);

    }

Thank you!