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

pisoFoam - Too high Courant Numbers

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

LinkBack Thread Tools Search this Thread Display Modes
Old   December 30, 2014, 08:47
Default pisoFoam - Too high Courant Numbers
Join Date: Dec 2014
Posts: 50
Rep Power: 10
Harak is on a distinguished road
Hi all,

I want to simulate a 3d case where a flow in a channel with a square cross section passes a barrier. See the sketch below.

It's my very first case after going through the tutorials

So, I decided to use pisoFoam since I believe it's well suitable for my case.
My problem is that I get incredible high Courant numbers so that my calculation crashes.

I've put some info for you. I hope you guys can help me




And that's the error I get:
Time = 0.025

Courant Number mean: 3.03546e+85 max: 1.06297e+89
smoothSolver: Solving for Uy, Initial residual = 1, Final residual = 1.44087e-15, No Iterations 1
#0 Foam::error:rintStack(Foam::Ostream&) at ??:?
#1 Foam::sigFpe::sigHandler(int) at ??:?
#2 in "/lib/x86_64-linux-gnu/"
#3 Foam::GAMGSolver::scale(Foam::Field<double>&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, Foam::Field<double> const&, unsigned char) const at ??:?
#4 Foam::GAMGSolver::Vcycle(Foam::PtrList<Foam::lduMa trix::smoother> const&, Foam::Field<double>&, Foam::Field<double> const&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::Field<double>&, Foam::PtrList<Foam::Field<double> >&, Foam::PtrList<Foam::Field<double> >&, unsigned char) const at ??:?
#5 Foam::GAMGSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#6 Foam::fvMatrix<double>::solveSegregated(Foam::dict ionary const&) at ??:?
#7 Foam::fvMatrix<double>::solve(Foam::dictionary const&) at ??:?
#8 Foam::fvMatrix<double>::solve() at ??:?
at ??:?
#10 __libc_start_main in "/lib/x86_64-linux-gnu/"
at ??:?
Floating point exception (core dumped)
Harak is offline   Reply With Quote

Old   December 30, 2014, 12:32
Senior Member
Alhasan's Avatar
Hasan K.J.
Join Date: Dec 2011
Location: Bristol, United Kingdom
Posts: 200
Rep Power: 14
Alhasan is on a distinguished road

The easiest option is to use pimpleFoam where you can set the maximum courant number.

If you want use pisoFoam, you need to reduce the time step (delta t) in order to reduce the courant number.

Kind Regards,
Hasan K.J
"Real knowledge is to know the extent of one's ignorance." - Confucius
Alhasan is offline   Reply With Quote


core dumped, floating point error, floating point exception, pisofoam

Thread Tools Search this Thread
Search this Thread:

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

Similar Threads
Thread Thread Starter Forum Replies Last Post
Visualizing mesh regions with high courant number? Endel OpenFOAM 2 August 20, 2013 03:57
Does MULES::implicitsolve allow for Courant numbers > 1? Philip OpenFOAM 4 May 31, 2011 06:25
Very high courant number, time step continuity errors and bounding epsilon erncyc OpenFOAM Running, Solving & CFD 1 March 24, 2011 11:00
Energy balance error with high prandtl numbers. David_010 OpenFOAM Running, Solving & CFD 0 January 20, 2011 05:21
high and low reynolds numbers - what does it mean Bo Busk Jensen Main CFD Forum 3 September 5, 2002 05:41

All times are GMT -4. The time now is 12:16.