assign polynominal function for water viscosity in OpenFOAM-2.2.x?
 

Dear OpenFOAMers,

I am wondering if anyone here has tried using polynominal function for water viscosity in OpenFOAM-2.2.x, specifically for compressibleInterFoam solver.

I adapted the polynominal implementation from one of the tutorial cases, but, got the following error. Do I need to link compressibleInterFoam solver to specific libraries?

Thanks!

Pei-Ying
---------------------------------------
hsieh@hillsdale:~/OpenFOAM/hsieh-2.2.x/run/compressibleInterFoam/newSensorManifold2x> compressibleInterFoam
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.2.x |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 2.2.x-444404f722b4
Exec : compressibleInterFoam
Date : Mar 30 2013
Time : 06:52:48
Host : "hillsdale"
PID : 5711
Case : /home/hsieh/OpenFOAM/hsieh-2.2.x/run/compressibleInterFoam/newSensorManifold2x
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0


Reading g

PIMPLE: Operating solver in PISO mode

Reading field p_rgh

Reading field U

Reading/calculating face flux field phi

Constructing twoPhaseMixtureThermo

Selecting thermodynamics package
{
type heRhoThermo;
mixture pureMixture;
transport polynomial;
thermo hPolynomial;
equationOfState icoPolynomial;
specie specie;
energy sensibleEnthalpy;
}



--> FOAM FATAL IO ERROR:
keyword rhoCoeffs<8> is undefined in dictionary "/home/hsieh/OpenFOAM/hsieh-2.2.x/run/compressibleInterFoam/newSensorManifold2x/constant/thermophysicalPropertieswater.mixture.equationOfSt ate"

file: /home/hsieh/OpenFOAM/hsieh-2.2.x/run/compressibleInterFoam/newSensorManifold2x/constant/thermophysicalPropertieswater.mixture.equationOfSt ate from line 38 to line 38.

From function dictionary::lookupEntry(const word&, bool, bool) const
in file db/dictionary/dictionary.C at line 402.

FOAM exiting

HI,

I just found out that I made a mistake in the water properties dictionary. But, now, I have a new error:

hsieh@hillsdale:~/OpenFOAM/hsieh-2.2.x/run/compressibleInterFoam/newSensorManifold2x> compressibleInterFoam
/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.2.x |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 2.2.x-444404f722b4
Exec : compressibleInterFoam
Date : Mar 30 2013
Time : 07:12:59
Host : "hillsdale"
PID : 6011
Case : /home/hsieh/OpenFOAM/hsieh-2.2.x/run/compressibleInterFoam/newSensorManifold2x
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0


Reading g

PIMPLE: Operating solver in PISO mode

Reading field p_rgh

Reading field U

Reading/calculating face flux field phi

Constructing twoPhaseMixtureThermo

Selecting thermodynamics package
{
type heRhoThermo;
mixture pureMixture;
transport polynomial;
thermo hPolynomial;
equationOfState icoPolynomial;
specie specie;
energy sensibleEnthalpy;
}

Selecting thermodynamics package
{
type heRhoThermo;
mixture pureMixture;
transport const;
thermo hConst;
equationOfState perfectGas;
specie specie;
energy sensibleInternalEnergy;
}



--> FOAM FATAL ERROR:
Supported energy type is ewater, thermodynamics package provides hwater

From function water
in file basicThermo/basicThermo.C at line 332.

FOAM exiting
--------------------------------------
Any suggestions will be highly appreciated.

Pei-Ying

HI,

OK, I changed

energy sensibleEnthalpy;

to

energy sensibleInternalEnthalpy;

and now the case is running. I am not sure if this is the right thing to do. I am not faminar with the thermo lib.

Pei-Ying