CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Running, Solving & CFD

Negative value of k causing simulation to stop

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

Reply
 
LinkBack Thread Tools Display Modes
Old   October 16, 2008, 21:00
Default Hi all I am running supersoni
  #1
Member
 
Velan
Join Date: Mar 2009
Location: India
Posts: 50
Rep Power: 8
velan is on a distinguished road
Hi all
I am running supersonic flow along with exit back pressure condition. I am using RNGEpsilon turbulence method. I am getting the good flow development. But suddenly openfoam quits with some error. When i checked the error, only differnce i found was k is negative. This happens after flow was developed fully but not yet converged. Is there any solution to solve this problem ?


Time = 0.00108087

Courant Number mean: 0.00170796 max: 0.0971893
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for Ux, Initial residual = 7.11502e-05, Final residual = 1.35237e-10, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 0.00011828, Final residual = 5.53025e-11, No Iterations 1
DILUPBiCG: Solving for h, Initial residual = 0.000251239, Final residual = 9.49742e-11, No Iterations 1
DILUPBiCG: Solving for p, Initial residual = 0.000376528, Final residual = 1.33934e-14, No Iterations 2
DILUPBiCG: Solving for p, Initial residual = 2.45546e-06, Final residual = 5.77244e-15, No Iterations 2
DILUPBiCG: Solving for p, Initial residual = 4.80871e-08, Final residual = 6.8838e-15, No Iterations 2
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 2.13924e-12, global = 1.25215e-13, cumulative = -2.11148e-12
DILUPBiCG: Solving for p, Initial residual = 1.28753e-06, Final residual = 1.27746e-14, No Iterations 2
DILUPBiCG: Solving for p, Initial residual = 5.51159e-09, Final residual = 5.40315e-15, No Iterations 2
DILUPBiCG: Solving for p, Initial residual = 1.26357e-09, Final residual = 2.10007e-14, No Iterations 2
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.15167e-11, global = -6.40266e-13, cumulative = -2.75174e-12
DILUPBiCG: Solving for epsilon, Initial residual = 6.54999e-05, Final residual = 1.64537e-10, No Iterations 1
DILUPBiCG: Solving for k, Initial residual = 0.000147741, Final residual = 2.10956e-10, No Iterations 1
ExecutionTime = 1.68 s ClockTime = 2 s

Time = 0.00108092

Courant Number mean: 0.00177227 max: 0.0776741
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for Ux, Initial residual = 7.13326e-05, Final residual = 6.69864e-11, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 0.000120215, Final residual = 5.63942e-11, No Iterations 1
DILUPBiCG: Solving for h, Initial residual = 0.000198837, Final residual = 1.50863e-10, No Iterations 1
DILUPBiCG: Solving for p, Initial residual = 0.000254745, Final residual = 1.13091e-14, No Iterations 2
DILUPBiCG: Solving for p, Initial residual = 2.72206e-07, Final residual = 4.31035e-15, No Iterations 2
DILUPBiCG: Solving for p, Initial residual = 7.12951e-09, Final residual = 3.73677e-15, No Iterations 2
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.11735e-12, global = 3.73839e-14, cumulative = -2.71436e-12
DILUPBiCG: Solving for p, Initial residual = 3.15769e-06, Final residual = 8.8416e-15, No Iterations 2
DILUPBiCG: Solving for p, Initial residual = 8.59995e-09, Final residual = 5.96147e-15, No Iterations 2
DILUPBiCG: Solving for p, Initial residual = 1.68437e-09, Final residual = 6.72562e-15, No Iterations 2
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.20992e-12, global = 1.77206e-13, cumulative = -2.53715e-12
DILUPBiCG: Solving for epsilon, Initial residual = 8.10672e-05, Final residual = 1.61564e-10, No Iterations 1
DILUPBiCG: Solving for k, Initial residual = 0.000190269, Final residual = 2.11132e-10, No Iterations 1
bounding k, min: -0.257568 max: 2202.59 average: 159.908
ExecutionTime = 1.76 s ClockTime = 2 s

Time = 0.00108097

Courant Number mean: 0.00178108 max: 0.0838975
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for Ux, Initial residual = 7.17025e-05, Final residual = 7.1539e-11, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 0.000123248, Final residual = 5.75958e-11, No Iterations 1
DILUPBiCG: Solving for h, Initial residual = 0.000208497, Final residual = 1.2543e-10, No Iterations 1
DILUPBiCG: Solving for p, Initial residual = 0.0145355, Final residual = 8.33417e-15, No Iterations 2
DILUPBiCG: Solving for p, Initial residual = 0.000477816, Final residual = 7.66563e-15, No Iterations 2
DILUPBiCG: Solving for p, Initial residual = 9.7461e-05, Final residual = 8.97423e-15, No Iterations 2
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 2.62836e-12, global = 2.72003e-13, cumulative = -2.26515e-12
DILUPBiCG: Solving for p, Initial residual = 0.00632594, Final residual = 8.8338e-15, No Iterations 2
DILUPBiCG: Solving for p, Initial residual = 6.28529e-05, Final residual = 8.58555e-15, No Iterations 2
DILUPBiCG: Solving for p, Initial residual = 1.99815e-07, Final residual = 2.34281e-14, No Iterations 2
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 8.98943e-12, global = 2.80106e-13, cumulative = -1.98504e-12
DILUPBiCG: Solving for epsilon, Initial residual = 0.00343727, Final residual = 1.58363e-10, No Iterations 1
bounding epsilon, min: -12607.9 max: 4.37996e+07 average: 618204
DILUPBiCG: Solving for k, Initial residual = 0.00864829, Final residual = 2.08733e-10, No Iterations 1
ExecutionTime = 1.86 s ClockTime = 2 s

Time = 0.00108102

Courant Number mean: 0.0018193 max: 0.872636
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for Ux, Initial residual = 0.000284262, Final residual = 6.7951e-11, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 0.0015673, Final residual = 6.46087e-11, No Iterations 1
DILUPBiCG: Solving for h, Initial residual = 0.00829178, Final residual = 4.08737e-10, No Iterations 1
DILUPBiCG: Solving for p, Initial residual = 0.0170355, Final residual = 7.43462e-15, No Iterations 3
DILUPBiCG: Solving for p, Initial residual = 0.000160206, Final residual = 6.32461e-15, No Iterations 2
DILUPBiCG: Solving for p, Initial residual = 2.24195e-08, Final residual = 5.34237e-15, No Iterations 2
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.9573e-12, global = -2.59168e-13, cumulative = -2.24421e-12
DILUPBiCG: Solving for p, Initial residual = 6.45607e-06, Final residual = 2.06521e-14, No Iterations 2
DILUPBiCG: Solving for p, Initial residual = 8.2883e-09, Final residual = 6.84259e-14, No Iterations 2
DILUPBiCG: Solving for p, Initial residual = 1.18195e-09, Final residual = 5.45102e-15, No Iterations 2
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
time step continuity errors : sum local = 1.88121e-12, global = 1.1234e-13, cumulative = -2.13187e-12
DILUPBiCG: Solving for epsilon, Initial residual = 0.000681664, Final residual = 1.63971e-10, No Iterations 1
bounding epsilon, min: -399912 max: 4.37998e+07 average: 619573
DILUPBiCG: Solving for k, Initial residual = 0.00185613, Final residual = 2.23178e-10, No Iterations 1
bounding k, min: -198.053 max: 4258.97 average: 161.253
ExecutionTime = 1.94 s ClockTime = 2 s

Time = 0.00108107

Courant Number mean: 0.00184301 max: 0.170472
diagonal: Solving for rho, Initial residual = 0, Final residual = 0, No Iterations 0
DILUPBiCG: Solving for Ux, Initial residual = 0.000111599, Final residual = 6.88509e-11, No Iterations 1
DILUPBiCG: Solving for Uy, Initial residual = 0.000439034, Final residual = 6.90532e-11, No Iterations 1
DILUPBiCG: Solving for h, Initial residual = 0.00458797, Final residual = 5.4244e-10, No Iterations 1


Maximum number of iterations exceeded#0 Foam::error::printStack(Foam:stream&) in "/home/vc/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libOpenFOAM.so"
#1 Foam::error::abort() in "/home/vc/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libOpenFOAM.so"
#2 Foam::hThermo<foam::puremixture<foam::consttranspo rt<foam::speciethermo<foam::hc onstthermo<foam::perfectgas> > > > >::calculate() in "/home/vc/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libbasicThermophysicalModels.s o"
#3 Foam::hThermo<foam::puremixture<foam::consttranspo rt<foam::speciethermo<foam::hc onstthermo<foam::perfectgas> > > > >::correct() in "/home/vc/OpenFOAM/OpenFOAM-1.5/lib/linuxGccDPOpt/libbasicThermophysicalModels.s o"
#4 main in "/home/vc/OpenFOAM/OpenFOAM-1.5/applications/bin/linuxGccDPOpt/sonicTurbFoam"
#5 __libc_start_main in "/lib/libc.so.6"
#6 Foam::regIOobject::readIfModified() in "/home/vc/OpenFOAM/OpenFOAM-1.5/applications/bin/linuxGccDPOpt/sonicTurbFoam"


From function specieThermo<thermo>::T(scalar f, scalar T0, scalar (specieThermo<thermo>::*F)(const scalar) const, scalar (specieThermo<thermo>::*dFdT)(const scalar) const) const
in file /home/vc/OpenFOAM/OpenFOAM-1.5/src/thermophysicalModels/specie/lnInclude/specieT hermoI.H at line 78.

FOAM aborting

Abort
velan is offline   Reply With Quote

Old   October 17, 2008, 05:36
Default Your simulation has blown up a
  #2
Senior Member
 
Hrvoje Jasak
Join Date: Mar 2009
Location: London, England
Posts: 1,763
Rep Power: 21
hjasak will become famous soon enough
Your simulation has blown up a while before - have a look at your turbulence levels.

What you need to do is to adjust the discretisation on the k and epsilon. FOr starters, use upwind on convection terms, and add some Gauss linear limited 0.5 on diffusion for convection and diffusion.

This should fix it - good luck!

Hrv
__________________
Hrvoje Jasak
Providing commercial FOAM/OpenFOAM and CFD Consulting: http://wikki.co.uk
hjasak 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
Can Fluent start a simulation from a negative time Ted Chu FLUENT 10 November 26, 2008 22:06
Expert parameter to stop the fluid flow simulation KK CFX 1 February 25, 2008 17:29
External Aerodynamic Simulation+negative densities Ramkumar CD-adapco 2 February 19, 2008 05:18
How stop simulation if convergence is not reached? Andreas FLUENT 0 October 16, 2007 12:35
Outflow conditions causing ERROR #064 Coriolius CD-adapco 3 October 10, 2004 10:57


All times are GMT -4. The time now is 09:42.