|
[Sponsors] |
December 30, 2014, 08:47 |
pisoFoam - Too high Courant Numbers
|
#1 |
Member
Join Date: Dec 2014
Posts: 50
Rep Power: 11 |
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) |
|
Tags |
core dumped, floating point error, floating point exception, pisofoam |
|
|
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 |