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

icoFoam unstable for ahmed body

Register Blogs Community New Posts Updated Threads Search

 
 
LinkBack Thread Tools Search this Thread Display Modes
Prev Previous Post   Next Post Next
Old   June 5, 2014, 08:37
Default icoFoam unstable for ahmed body
  #1
New Member
 
Join Date: May 2014
Posts: 11
Rep Power: 12
rogibear is on a distinguished road
Hi, I'm trying to run a very simple ahmed body case with low Reynolds number (100 & 1000) with time steps of 0.001 & 0.0001 for each case. The mesh is pretty course so I'm not expecting to get useable results but I am trying to get the simulation to complete. The problem is the courant number keeps increasing >1 by the timestep doesn't seem to have much effect on it.

see below:
Code:
Courant Number mean: 0.000359148 max: 0.793021
DILUPBiCG:  Solving for Ux, Initial residual = 0.923248, Final residual = 3.32352e-09, No Iterations 2
DILUPBiCG:  Solving for Uy, Initial residual = 0.969197, Final residual = 2.0091e-09, No Iterations 2
DILUPBiCG:  Solving for Uz, Initial residual = 0.930231, Final residual = 9.23342e-10, No Iterations 2
DICPCG:  Solving for p, Initial residual = 0.17526, Final residual = 9.82767e-07, No Iterations 146
time step continuity errors : sum local = 6.91405e-09, global = 4.97837e-12, cumulative = -4.69581e-09
DICPCG:  Solving for p, Initial residual = 0.000412786, Final residual = 8.75322e-07, No Iterations 77
time step continuity errors : sum local = 6.42496e-09, global = 1.81337e-10, cumulative = -4.51448e-09
ExecutionTime = 1028.69 s  ClockTime = 1031 s

Time = 0.0351

Courant Number mean: 0.000366532 max: 0.7604
DILUPBiCG:  Solving for Ux, Initial residual = 0.969019, Final residual = 3.0867e-09, No Iterations 2
DILUPBiCG:  Solving for Uy, Initial residual = 0.969175, Final residual = 1.89008e-09, No Iterations 2
DILUPBiCG:  Solving for Uz, Initial residual = 0.935075, Final residual = 1.67906e-09, No Iterations 2
DICPCG:  Solving for p, Initial residual = 0.174915, Final residual = 9.83657e-07, No Iterations 146
time step continuity errors : sum local = 7.24729e-09, global = -5.65332e-12, cumulative = -4.52013e-09
DICPCG:  Solving for p, Initial residual = 0.000413647, Final residual = 9.5122e-07, No Iterations 35
time step continuity errors : sum local = 7.31108e-09, global = -4.48088e-10, cumulative = -4.96822e-09
ExecutionTime = 1032.19 s  ClockTime = 1034 s

Time = 0.0352

Courant Number mean: 0.000373678 max: 0.872445
DILUPBiCG:  Solving for Ux, Initial residual = 0.925273, Final residual = 4.64637e-09, No Iterations 2
DILUPBiCG:  Solving for Uy, Initial residual = 0.968993, Final residual = 2.77035e-09, No Iterations 2
DILUPBiCG:  Solving for Uz, Initial residual = 0.930175, Final residual = 1.26008e-09, No Iterations 2
DICPCG:  Solving for p, Initial residual = 0.174619, Final residual = 9.7952e-07, No Iterations 146
time step continuity errors : sum local = 7.55802e-09, global = 5.9816e-12, cumulative = -4.96224e-09
DICPCG:  Solving for p, Initial residual = 0.000445268, Final residual = 9.10332e-07, No Iterations 77
time step continuity errors : sum local = 7.3276e-09, global = 2.08616e-10, cumulative = -4.75362e-09
ExecutionTime = 1035.99 s  ClockTime = 1038 s

Time = 0.0353

Courant Number mean: 0.000381645 max: 0.818073
DILUPBiCG:  Solving for Ux, Initial residual = 0.96842, Final residual = 4.16886e-09, No Iterations 2
DILUPBiCG:  Solving for Uy, Initial residual = 0.968974, Final residual = 2.53332e-09, No Iterations 2
DILUPBiCG:  Solving for Uz, Initial residual = 0.935703, Final residual = 2.32248e-09, No Iterations 2
DICPCG:  Solving for p, Initial residual = 0.174293, Final residual = 9.80834e-07, No Iterations 146
time step continuity errors : sum local = 7.92459e-09, global = -6.76192e-12, cumulative = -4.76038e-09
DICPCG:  Solving for p, Initial residual = 0.000447096, Final residual = 9.42898e-07, No Iterations 35
time step continuity errors : sum local = 7.94622e-09, global = -3.75958e-10, cumulative = -5.13634e-09
ExecutionTime = 1039.07 s  ClockTime = 1041 s

Time = 0.0354

Courant Number mean: 0.000389596 max: 0.960098
DILUPBiCG:  Solving for Ux, Initial residual = 0.927154, Final residual = 6.49988e-09, No Iterations 2
DILUPBiCG:  Solving for Uy, Initial residual = 0.968825, Final residual = 3.79038e-09, No Iterations 2
DILUPBiCG:  Solving for Uz, Initial residual = 0.930027, Final residual = 2.14034e-09, No Iterations 2
DICPCG:  Solving for p, Initial residual = 0.174019, Final residual = 9.76577e-07, No Iterations 146
time step continuity errors : sum local = 8.26213e-09, global = 7.12199e-12, cumulative = -5.12922e-09
DICPCG:  Solving for p, Initial residual = 0.000480689, Final residual = 9.48658e-07, No Iterations 77
time step continuity errors : sum local = 8.37163e-09, global = 2.40243e-10, cumulative = -4.88898e-09
ExecutionTime = 1042.65 s  ClockTime = 1045 s

Time = 0.0355

Courant Number mean: 0.000398198 max: 0.877074
DILUPBiCG:  Solving for Ux, Initial residual = 0.967843, Final residual = 5.58373e-09, No Iterations 2
DILUPBiCG:  Solving for Uy, Initial residual = 0.968759, Final residual = 3.34602e-09, No Iterations 2
DILUPBiCG:  Solving for Uz, Initial residual = 0.936349, Final residual = 3.19647e-09, No Iterations 2
DICPCG:  Solving for p, Initial residual = 0.17371, Final residual = 9.78324e-07, No Iterations 146
time step continuity errors : sum local = 8.66548e-09, global = -8.01879e-12, cumulative = -4.897e-09
DICPCG:  Solving for p, Initial residual = 0.000483549, Final residual = 9.38136e-07, No Iterations 35
time step continuity errors : sum local = 8.66638e-09, global = -2.75065e-10, cumulative = -5.17206e-09
ExecutionTime = 1045.74 s  ClockTime = 1048 s

Time = 0.0356

Courant Number mean: 0.000407007 max: 1.06032
DILUPBiCG:  Solving for Ux, Initial residual = 0.928911, Final residual = 9.04348e-09, No Iterations 2
DILUPBiCG:  Solving for Uy, Initial residual = 0.968677, Final residual = 5.06989e-09, No Iterations 2
DILUPBiCG:  Solving for Uz, Initial residual = 0.929783, Final residual = 6.0862e-09, No Iterations 2
DICPCG:  Solving for p, Initial residual = 0.173457, Final residual = 9.73925e-07, No Iterations 146
time step continuity errors : sum local = 9.032e-09, global = 8.40815e-12, cumulative = -5.16365e-09
DICPCG:  Solving for p, Initial residual = 0.000519323, Final residual = 9.90508e-07, No Iterations 77
time step continuity errors : sum local = 9.58036e-09, global = 2.7687e-10, cumulative = -4.88678e-09
ExecutionTime = 1049.41 s  ClockTime = 1051 s

If you need me to upload the full case I'll put a dropbox link up.

Thanks,

Adam
rogibear is offline   Reply With Quote

 


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 Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
6dof model - floating body flinde FLUENT 13 March 6, 2021 14:32
How to add a source term (body force) to icoFoam? chris Liu OpenFOAM Programming & Development 4 October 27, 2014 15:38
Beginner questions - modelling lift and drag on a towed body ABF Main CFD Forum 15 October 19, 2014 23:27
[Need some advise] Am I doing right?? Simulating a rotational mixer setasena STAR-CCM+ 4 March 10, 2013 10:32
Kubuntu uses dash breaks All scripts in tutorials platopus OpenFOAM Bugs 8 April 15, 2008 07:52


All times are GMT -4. The time now is 02:10.