[vasava]

I want to write a new viscosity model where viscosity is also a function of volume fraction of one of the phases in a multiphase flow.

I have played a bit with CrossPowerLaw. I made some minor modifications and everything was nice. I understood that this piece of code does is the main equation for the viscosity model.

Foam::tmp<Foam::volScalarField>
Foam::viscosityModels::NewCrossPowerLaw::calcNu() const
{
    return (nu0_ - nuInf_)/(scalar(1) + pow(m_*strainRate(), n_)) + nuInf_;
}

My question is how do I add information about volume fraction of certain phase into this code?? For example I want to multiply the final viscosity with volume fraction of one of the phases.

My C++ skills are way too basic so you would also have to explain about changes in the header file.

Thanks!!

PS:

1. I would use this with twoPhaseEulerFoam or multiphaseEulerFoam.
2. I have openfoam 2.3.0

[massive_turbulence]

Originally Posted by vasava

I want to write a new viscosity model where viscosity is also a function of volume fraction of one of the phases in a multiphase flow.

I have played a bit with CrossPowerLaw. I made some minor modifications and everything was nice. I understood that this piece of code does is the main equation for the viscosity model.

Foam::tmp<Foam::volScalarField>
Foam::viscosityModels::NewCrossPowerLaw::calcNu() const
{
    return (nu0_ - nuInf_)/(scalar(1) + pow(m_*strainRate(), n_)) + nuInf_;
}

My question is how do I add information about volume fraction of certain phase into this code?? For example I want to multiply the final viscosity with volume fraction of one of the phases.

My C++ skills are way too basic so you would also have to explain about changes in the header file.

Thanks!!

PS:

1. I would use this with twoPhaseEulerFoam or multiphaseEulerFoam.
2. I have openfoam 2.3.0

In the code which variable is the final viscosity?

Foam::tmp<Foam::volScalarField>
Foam::viscosityModels::NewCrossPowerLaw::calcNu() const
{
    return (nu0_ - nuInf_)/(scalar(1) + pow(m_*strainRate(), n_)) + nuInf_;
}

I'm assuming it's

nuInf_

so you would create it (volume fraction)

double vol_frac = function_for_calculation_of_vol_frac();

and then plug it into the code

Foam::tmp<Foam::volScalarField>
Foam::viscosityModels::NewCrossPowerLaw::calcNu() const
{
    double vol_frac = function_for_calculation_of_vol_frac();

    return (nu0_ - nuInf_*vol_frac)/(scalar(1) + pow(m_*strainRate(), n_)) + nuInf_*vol_frac;
}

I'm not sure what the type for nuInf_ is. What file is this in by the way?


For the header file you would need to add the function for

function_for_calculation_of_vol_frac()

if not then you would need to include the volume fraction calculation within the function Foam::viscosityModels::NewCrossPowerLaw::calcNu() and there's no need to change any header file.

[vasava]

Thanks for your response.

This code as I mentioned is from CrossPowerLaw. I have included this code just for illustration. Right now I do not know how to include the information about the volume fraction of a phase.

[massive_turbulence]

Originally Posted by vasava

Thanks for your response.

This code as I mentioned is from CrossPowerLaw. I have included this code just for illustration. Right now I do not know how to include the information about the volume fraction of a phase.

How is the 'volume fraction of the phase' calculated though, surely there's a simple formula (I think I remember something about this from thermo)? You might be able to read data from a dat file using fstream although if time dependence is necessary or some other variable changes then the function call would get more difficult to implement.