[msaravia]

Hi everyone!

I am taking as a template the viscosityModel class to write a class for magnetic materials constitutive laws.

I have some VERY BASIC doubts about the handling of the viscosity in, for example, the Newtonian model. The code is:

Foam::viscosityModels::Newtonian::Newtonian
(
    const word& name,
    const dictionary& viscosityProperties,
    const volVectorField& U,
    const surfaceScalarField& phi
)
:
    viscosityModel(name, viscosityProperties, U, phi),
    nu0_("nu", dimViscosity, viscosityProperties_),
    nu_
    (
        IOobject
        (
            name,
            U_.time().timeName(),
            U_.db(),
            IOobject::NO_READ,
            IOobject::NO_WRITE
        ),
        U_.mesh(),
        nu0_
    )
{}


// * * * * * * * * * * * * * * Member Functions  * * * * * * * * * * * * * * //

bool Foam::viscosityModels::Newtonian::read
(
    const dictionary& viscosityProperties
)
{
    viscosityModel::read(viscosityProperties);

    viscosityProperties_.lookup("nu") >> nu0_;
    nu_ = nu0_;

    return true;
}

My doubts are:

i) Does the operator >> in

viscosityProperties_.lookup("nu") >> nu0_

assigns the found value to mur0_ ?

ii) If nu_ is a private volScalarField and nu0_ is a dimensionedScalar, why the line mur_ = mur0_ works?


Thanks in advance!