[lindstroem]

Hi everybody,

I want to extend a solver from single to multiphase flows. At one point I am stuck because of the reading of the transportPropterties.

The "old" code says:

dimensionedScalar nu(transportProperties.lookup("nu"));

If I get it right, I will have a nu1 and nu2 for multiphase flows. Can I only access them via createFields or can it be read by transportProperties...?

Thanks for your help!
Lindstroem

[marupio]

createFields reads from transportProperties. If you want nu1 and nu2 defined, you need to modify your solver - and why not change createFields?

[lindstroem]

Hi David,

thanks for your suggestion. Actually I need them in a function where I calculate some forces (outside the solver). Up to now, I did not use createFields for this function, only for the solver. I hoped there was an easy way to read them just like for a single phase flow.

Thanks..

[Bernhard]

You can have a peek in the createFields.H file of interFoam to see how it is done.

[lindstroem]

Originally Posted by Bernhard

You can have a peek in the createFields.H file of interFoam to see how it is done.

So in a similar manner like

 twoPhaseMixture twoPhaseProperties(U, phi, "alpha1");
    const dimensionedScalar& rho1 = twoPhaseProperties.rho1();
    const dimensionedScalar& rho2 = twoPhaseProperties.rho2();

Ok, then I need to extend the functions as expected. I hoped to have sth like

dimensionedScalar nu1(transportProperties.lookup("nu1"));

Thanks!

[Bernhard]

That second option is of course also possible. Did you try it already?

[lindstroem]

Yes, but... it compiles, but I get the error

keyword nu1 is undefined in dictionary "[..]/constant/transportProperties"

file: [..]constant/transportProperties from line 19 to line 32.

[lindstroem]

twoPhaseMixture.H provides a volScalarField nu so thats the way i will get access to it:

twoPhaseMixture twoPhaseProperties(U, phi, "alpha1");
const volScalarField& nu = twoPhaseProperties.nu();

Thanks for your help!

[Bernhard]

Originally Posted by lindstroem

Yes, but... it compiles, but I get the error

keyword nu1 is undefined in dictionary "[..]/constant/transportProperties"

file: [..]constant/transportProperties from line 19 to line 32.

Which means that you have to provide a nu1 in the transportProperties file, but you already found your solution.

[lindstroem]

When I provided a nu1 it was complaining that it did not find a nu (which i did not call any longer from my force-calculation-function). Yes, I'll try it with the scalarField.

Thanks!