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

pisoFoam - Too high Courant Numbers

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

Reply
 
LinkBack Thread Tools Display Modes
Old   December 30, 2014, 09:47
Default pisoFoam - Too high Courant Numbers
  #1
Member
 
Join Date: Dec 2014
Posts: 33
Rep Power: 2
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

fvSolution:
https://dl.dropboxusercontent.com/u/59780963/fvSolution

fvSchemes:
https://dl.dropboxusercontent.com/u/59780963/fvSchemes

checkMesh:
https://dl.dropboxusercontent.com/u/59780963/checkMesh

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/libc.so.6"
#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 ??:?
#9
at ??:?
#10 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#11
at ??:?
Floating point exception (core dumped)
Harak is offline   Reply With Quote

Old   December 30, 2014, 13:32
Default
  #2
Senior Member
 
Alhasan's Avatar
 
Hasan K.J.
Join Date: Dec 2011
Location: Bristol, United Kingdom
Posts: 181
Rep Power: 5
Alhasan is on a distinguished road
Hey,

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

Reply

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

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
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 12:00
Energy balance error with high prandtl numbers. David_010 OpenFOAM Running, Solving & CFD 0 January 20, 2011 06: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 19:53.