CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Running, Solving & CFD

Supercritical flow in a duct_pressure blowup at initial time step

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

LinkBack Thread Tools Display Modes
Old   March 7, 2017, 12:49
Default Supercritical flow in a duct_pressure blowup at initial time step
New Member
Mahdi Nabil
Join Date: Sep 2015
Posts: 7
Rep Power: 3
Mhmnabil is on a distinguished road
Hi everyone,

I am trying to simulate supercritical CO2 turbulent flow and heat transfer in a duct using buoyantPimpleFOAM. So far, I have been using const fluid properties and it was working fine.
Currently, I am using a new EOS in which rho and psi (compressiblity) are both functions of pressure and temperature (since we are in a supercritical region. In other words,
psi=drho/dp = psi(p,T)
And, I am sure about the psi formulation and implementation in the code which is the exact analytical differentiation of rho wrt p. However, when I use psi(p,T) it makes the pressure to be unstable at the very first time step and the code blows up. I've printed the output:

diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
PIMPLE: iteration 1
DILUPBiCG: Solving for Ux, Initial residual = 1, Final residual = 5.15665e-08, No Iterations 2
DILUPBiCG: Solving for Uy, Initial residual = 1, Final residual = 4.17258e-07, No Iterations 3
DILUPBiCG: Solving for Uz, Initial residual = 1, Final residual = 3.51138e-07, No Iterations 3
DILUPBiCG: Solving for h, Initial residual = 0.999991, Final residual = 1.05686e-07, No Iterations 4
GAMG: Solving for p_rgh, Initial residual = 0.999569, Final residual = 4.66604e+56, No Iterations 1000
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 123.314, global = -15.475, cumulative = -15.475
************************************************** ********
**** Max Pressure = 1.50284e+61 **** Min Pressure = -8.4671e+60
**** Max Temperature = 313.004 **** Min Temperature = 312.996
************************************************** ********

I have attached my p, T, and U boundary conditions in the 0 folder. Does anyone have any idea if I am making any mistake in the boundary conditions? These boundary conditions work perfectly for a constant psi value, but blows up the pressure for psi(p,T).

Thanks a lot in advance.
Attached Files
File Type: zip (2.2 KB, 0 views)
Mhmnabil is offline   Reply With Quote


blowing up, boundary condition, buoyantpimplefoam, pressure, turbulent flow

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
Star cd es-ice solver error ernarasimman STAR-CD 2 September 12, 2014 00:01
Help for the small implementation in turbulence model shipman OpenFOAM Programming & Development 25 March 19, 2014 11:08
AMI interDyMFoam for mixer nu problem danny123 OpenFOAM Programming & Development 8 September 6, 2013 02:34
pisoFoam with k-epsilon turb blows up - Some questions Heroic OpenFOAM Running, Solving & CFD 26 December 17, 2012 04:34
IcoFoam parallel woes msrinath80 OpenFOAM Running, Solving & CFD 9 July 22, 2007 02:58

All times are GMT -4. The time now is 20:41.