[vkrastev]

Hi all, I'm trying to modify the chtMultiregionFoam solver of OpenFOAM 1.6 adding some of the characteristics of the 1.7.1 version. Apart from the details, I need the 1.6 solver to recognize p_rgh as the new pressure variable for which the solution procedure is performed. I've modified the pEqn.H file in order to reach the goal, but after a seemingly correct compilation when I try to run a test case this is the error message coming out:

Solving for fluid region fluid1
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 2.5150392e-08, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 1, Final residual = 9.5553996e-09, No Iterations 2
DILUPBiCG: Solving for Uz, Initial residual = 1, Final residual = 6.502438e-09, No Iterations 2
DILUPBiCG: Solving for h, Initial residual = 1, Final residual = 3.847947e-09, No Iterations 2
Min/max T:300 363


keyword p is undefined in dictionary "/work/alh/vesselin/myMultiRegionCase_4_liquid/system/fluid1/fvSolution::solvers"

file: /work/alh/vesselin/myMultiRegionCase_4_liquid/system/fluid1/fvSolution::solvers from line 21 to line 124.

From function dictionary::subDict(const word& keyword) const
in file db/dictionary/dictionary.C at line 449.

FOAM exiting

Any suggestion to fix the problem? Is it possible to redefine the keyword for the pressure variable?

Thank you in advance

V.

PS-Please, do not answer: "Switch to OF-1.7.1 instead of breaking your head on how to modify the 1.6 version...". I know this could be the simplest solution bur unfortunately I must do it the above way....

[alberto]

Hi Vesselin,

the keyword is based on the database entry of the field name, so that's what you have to change, not the name of the variable in the code. In other words the field name between " " in the field definition.

Best,

[vkrastev]

Originally Posted by alberto

Hi Vesselin,

the keyword is based on the database entry of the field name, so that's what you have to change, not the name of the variable in the code. In other words the field name between " " in the field definition.

Best,

Hi Alberto, and thanks a lot for your answer. Sorry If I'll pose you another (maybe) obvious question, but where can I find this database? I'm trying to search for it with the grep command inside the OpenFOAM folders, but after a couple of hours of searching still I cannot find anything that resembles a database with the fields' definition...

Thanks once again

V.

[alberto]

The database with the list of objects used in your run is created in memory at runtime by the code. Take a look at the objectRegistry.

[vkrastev]

Originally Posted by alberto

The database with the list of objects used in your run is created in memory at runtime by the code. Take a look at the objectRegistry.

Well, it seems that after a slight modification of the pEqn.H file now the p_rgh variable is recognized properly...However, after a first iteration loop on p_rgh, when It comes to solve the epsilon equation something strange happens and the code returns a floating point exception:

Solving for fluid region fluid1
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 2.5150392e-08, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 1, Final residual = 9.5553996e-09, No Iterations 2
DILUPBiCG: Solving for Uz, Initial residual = 1, Final residual = 6.502438e-09, No Iterations 2
DILUPBiCG: Solving for h, Initial residual = 1, Final residual = 3.847947e-09, No Iterations 2
Min/max T:300 363
GAMG: Solving for p_rgh, Initial residual = 0.81016096, Final residual = 0.0049955858, No Iterations 18
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors (fluid1): sum local = 3.2337637e-05, global = 1.1811838e-07, cumulative = 1.1811838e-07
GAMG: Solving for p_rgh, Initial residual = 0.0019349904, Final residual = 1.4132342e-05, No Iterations 12
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors (fluid1): sum local = 5.6645344e-06, global = -7.9588068e-07, cumulative = -6.7776229e-07
GAMG: Solving for p_rgh, Initial residual = 9.3519277e-05, Final residual = 4.8533869e-08, No Iterations 16
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors (fluid1): sum local = 1.8940969e-08, global = -1.140037e-09, cumulative = -6.7890233e-07
#0 Foam::error:rintStack(Foam::Ostream&) in "/work/alh/vesselin/OpenFOAM///OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#1 Foam::sigFpe::sigFpeHandler(int) in "/work/alh/vesselin/OpenFOAM///OpenFOAM-1.6/lib/linux64GccDPOpt/libOpenFOAM.so"
#2 ?? in "/lib64/libc.so.6"
#3 Foam::compressible::RASModels::epsilonWallFunction FvPatchScalarField::updateCoeffs() in "/work/alh/vesselin/OpenFOAM///OpenFOAM-1.6/lib/linux64GccDPOpt/libcompressibleRASModels.so"
#4 Foam::compressible::RASModels::kEpsilon::correct() in "/work/alh/vesselin/OpenFOAM///OpenFOAM-1.6/lib/linux64GccDPOpt/libcompressibleRASModels.so"
#5 main in "/work/alh/vesselin/OpenFOAM/OpenFOAM-1.6/applications/bin/linux64GccDPOpt/myChtMultiRegionFoam"
#6 __libc_start_main in "/lib64/libc.so.6"
#7 Foam::regIOobject::writeObject(Foam::IOstream::str eamFormat, Foam::IOstream::versionNumber, Foam::IOstream::compressionType) const in "/work/alh/vesselin/OpenFOAM/OpenFOAM-1.6/applications/bin/linux64GccDPOpt/myChtMultiRegionFoam"
Floating exception

Any ideas on what is going on? BCs and initial conditions are identical to those that work perfectly with chtMRF in OpenFOAM 1.7.1 (this is the updated version with the capability to handle liquid fluids governed by an incompressible polynomial equation of state)...

Thanks

V.