CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Programming & Development

p_rgh in OpenFOAM-1.6

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Display Modes
Old   March 29, 2011, 04:59
Default p_rgh in OpenFOAM-1.6
  #1
Senior Member
 
Vesselin Krastev
Join Date: Jan 2010
Location: University of Tor Vergata, Rome
Posts: 361
Rep Power: 10
vkrastev is on a distinguished road
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....
vkrastev is offline   Reply With Quote

Old   March 30, 2011, 04:31
Default
  #2
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,894
Rep Power: 26
alberto will become famous soon enoughalberto will become famous soon enough
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,
__________________
Alberto Passalacqua

GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as live DVD/USB, hard drive image and virtual image.
OpenQBMM - An open-source implementation of quadrature-based moment methods
alberto is offline   Reply With Quote

Old   March 30, 2011, 06:29
Default
  #3
Senior Member
 
Vesselin Krastev
Join Date: Jan 2010
Location: University of Tor Vergata, Rome
Posts: 361
Rep Power: 10
vkrastev is on a distinguished road
Quote:
Originally Posted by alberto View Post
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.
vkrastev is offline   Reply With Quote

Old   March 30, 2011, 15:00
Default
  #4
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,894
Rep Power: 26
alberto will become famous soon enoughalberto will become famous soon enough
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.
__________________
Alberto Passalacqua

GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as live DVD/USB, hard drive image and virtual image.
OpenQBMM - An open-source implementation of quadrature-based moment methods
alberto is offline   Reply With Quote

Old   March 31, 2011, 09:35
Default
  #5
Senior Member
 
Vesselin Krastev
Join Date: Jan 2010
Location: University of Tor Vergata, Rome
Posts: 361
Rep Power: 10
vkrastev is on a distinguished road
Quote:
Originally Posted by alberto View Post
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.
vkrastev is offline   Reply With Quote

Reply

Thread Tools
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are On
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
OpenFOAM 1.6 and 1.7 with interFoam, groovyBC give different strange results Arnoldinho OpenFOAM 7 December 9, 2010 17:29
kOmegaSST in openfoam 1.6 Gearb0x OpenFOAM 2 March 3, 2010 07:02
Cross-compiling OpenFOAM 1.6 on Linux for Windows 32 and 64bits with Mingw-w64 wyldckat OpenFOAM Announcements from Other Sources 7 January 19, 2010 16:39
OpenFOAM 1.6 package - CentOS 5.3 x86_64 linnemann OpenFOAM Installation 36 September 2, 2009 12:46
OpenFOAM Version 1.6 Released opencfd OpenFOAM Announcements from ESI-OpenCFD 0 July 27, 2009 17:55


All times are GMT -4. The time now is 08:26.