[styleworker]

Hello foamers!

I'm quite familiar with using OpenFoam, but now I want to go further and learn more about programming in OF. At the moment I'm playing around with MRconjugateHeatFoam. I tried to extend it to a two phase flow solver, but I'm not able to call the twoPhaseProperties function.

This part of the boundary function is of interest:

Foam::tmp<Foam::scalarField>
Foam::regionCoupleTemperatureFvPatchScalarField::flux() const
{

    const dimensionedScalar& rho1 = twoPhaseProperties.rho1();

    const fvPatchScalarField& T = *this;

    return rho1.value()*T.snGrad();
}

At the beginning of the file, the header file twoPhaseMixture.H is included. Somehow I'm getting the error: »twoPhaseProperties« is not defined

Did I included the wrong header file?

[sharonyue]

I think you should construct transportProperties like this :

incompressibleTwoPhaseMixture twoPhaseProperties(U, phi);

If I am wrong correct me pls.

[styleworker]

Thank you for your reply. Indeed you're right.

It seems incompressibleTwoPhaseMixture was renamed to twoPhaseMixture in OF. In my case (OF-2.2.0) it should be:

twoPhaseMixture twoPhaseProperties(U, phi,"alpha1");

Certainly the U-, and phi-field has to be defined/read:

const volVectorField& U = db().lookupObject<volVectorField>("U");
const surfaceScalarField& phi = db().lookupObject<surfaceScalarField>("phi");

The complete code-snippet:

Foam::tmp<Foam::scalarField>
Foam::regionCoupleTemperatureFvPatchScalarField::flux() const
{
    const volVectorField& U = db().lookupObject<volVectorField>("U");
    const surfaceScalarField& phi = db().lookupObject<surfaceScalarField>("phi");    
    twoPhaseMixture twoPhaseProperties(U, phi,"alpha1");

    const dimensionedScalar& rho1 = twoPhaseProperties.rho1();
    const fvPatchScalarField& T = *this; 
    return rho1.value()*T.snGrad();
}