[Alpha001]

Hello Foamers,

This is my first time asking a question here, so please be patient with any mistakes.

System:
OpenFOAM 7
Ubuntu 18.04LTS

Problem-Statement:
I am trying to couple interFoam to a in-house DEM solver.(Using preCICE)
When using pimpleFoam I could just hard set the values for density and viscocity.
I need the rho and nu field for calculating the Bouyant force on particles. But currently just trying to solve the issue with rho.
In interFoam we have two-phase, so different densities through the domain depending on where you see.

I am using the incompressible interFoam solver.


Changes
I have the following changes in the createFields.H

Info<< "Reading transportProperties\n" << endl;
immiscibleIncompressibleTwoPhaseMixture mixture(U, phi);

volScalarField& alpha1(mixture.alpha1());
volScalarField& alpha2(mixture.alpha2());

const dimensionedScalar& rho1 = mixture.rho1();
const dimensionedScalar& rho2 = mixture.rho2();


// Need to store rho for ddt(rho, U)
volScalarField rho
(
    IOobject
    (
        "rho",
        runTime.timeName(),
        mesh,
        IOobject::READ_IF_PRESENT
        IOobject::AUTO_WRITE
    ),
    alpha1*rho1 + alpha2*rho2
);
rho.oldTime();

I get the following error when I compile the solver using wmake command.

./createFields.H: In function ‘int main(int, char**)’:
./createFields.H:110:5: error: expected primary-expression before ‘(’ token
     (
     ^
./createFields.H:115:9: error: expected ‘)’ before ‘IOobject’
         IOobject::AUTO_WRITE
         ^~~~~~~~
In file included from exInterFoam.C:60:0:
createFields.H:110:5: error: expected primary-expression before ‘(’ token
     (
     ^
createFields.H:115:9: error: expected ‘)’ before ‘IOobject’
         IOobject::AUTO_WRITE
         ^~~~~~~~

I am a bit unsure of what I am messing up.
Any tips, advice and help is welcome.


Thank you.

[Wenyuan]

Hi,

It seems that a comma is missing before the newly added line.

[Alpha001]

Thank you.
It works.

[Alpha001]

In the above code in createFields.H, how can I add the nu field?

Info<< "Reading transportProperties\n" << endl;
immiscibleIncompressibleTwoPhaseMixture mixture(U, phi);

volScalarField& alpha1(mixture.alpha1());
volScalarField& alpha2(mixture.alpha2());

const dimensionedScalar& rho1 = mixture.rho1();
const dimensionedScalar& rho2 = mixture.rho2();

const dimensionedScalar& nu1 = mixture.nu1();
const dimensionedScalar& nu2 = mixture.nu2();

// Getting viscocity
volScalarField nu
(
    IOobject
    (
        "nu",
        runTime.timeName(),
        mesh,
        IOobject::NO_READ,
        IOobject::AUTO_WRITE
    ),
    alpha1*nu1 + alpha2*nu2
);

I added this, but I get the following error:

./createFields.H:107:40: error: ‘class Foam::immiscibleIncompressibleTwoPhaseMixture’ has no member named ‘nu1’; did you mean ‘nu’?
 const dimensionedScalar& nu1 = mixture.nu1();
                                        ^~~
                                        nu
./createFields.H:108:40: error: ‘class Foam::immiscibleIncompressibleTwoPhaseMixture’ has no member named ‘nu2’; did you mean ‘nu’?
 const dimensionedScalar& nu2 = mixture.nu2();
                                        ^~~
                                        nu

SO I did some digging and found that in src/transportModels/immiscibleIncompressibleTwoPhaseMixture/immiscibleIncompressibleTwoPhaseMixture.H.

[Wenyuan]

nu1() and nu2() are not defined in incompressibleTwoPhaseMixture, which is the parent class of immiscibleIncompressibleTwoPhaseMixture.

You can call nu() directly, which returns a tmp<volScalarField>.

[Alpha001]

Thank you. The solver is compiling.
I was also able to write out mu field in the similar way since it is present in incompressibleTwoPhaseMixture.H.

//- Return the dynamic laminar viscosity
        tmp<volScalarField> mu() const;

[absrocks]

Hi, I want to use turbulence->nut() and turbulence->nu() in interFoam solver, but I am getting the error while compiling

error: ‘class Foam::incompressibleInterPhaseTransportModel<Foam::immiscibleIncompressibleTwoPhaseMixture>’ has no member named ‘nu’
error: ‘class Foam::incompressibleInterPhaseTransportModel<Foam::immiscibleIncompressibleTwoPhaseMixture>’ has no member named ‘nut’

I am using openfoam2106 version. When I used the same in OpenFOAM-8 it works perfectly. Can anyone have an idea?
If I use mixture->nu() I did not get any error, but how would I access nut?

[Venky_94]

Originally Posted by Alpha001

Thank you. The solver is compiling.
I was also able to write out mu field in the similar way since it is present in incompressibleTwoPhaseMixture.H.

//- Return the dynamic laminar viscosity
        tmp<volScalarField> mu() const;

Hey, did you check the written out mu and nu values to see if they're right? Because I followed the same procedure as you and my mu and nu values are written out but they're constant values and not the calculated values at each time step.

Did you have to do something in addition to what you have mentioned earlier?