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

courant number explodes with pisoFoam solver

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

Reply
 
LinkBack Thread Tools Display Modes
Old   December 7, 2011, 08:23
Default courant number explodes with pisoFoam solver
  #1
New Member
 
Join Date: Dec 2011
Posts: 2
Rep Power: 0
tol24 is on a distinguished road
Hi,

I'm trying to simulate the unsteady flow arround a cylinder with a velocity inlet of 17 m/s. I'm using a pisoFoam solver with a RAS kEpsilon turbulent model. When I run the simulation the courant numbre explodes.

Time = 5e-05

Courant Number mean: 4.53303e+07 max: 1.0204e+12
DILUPBiCG: Solving for Ux, Initial residual = 0.985101, Final residual = 9.2599e-06, No Iterations 389
DILUPBiCG: Solving for Uy, Initial residual = 0.735666, Final residual = 0.188579, No Iterations 1001
DICPCG: Solving for p, Initial residual = 1, Final residual = 12.4404, No Iterations 1001
time step continuity errors : sum local = 4.28628e+29, global = -1.82086e+27, cumulative = -1.82086e+27
DICPCG: Solving for p, Initial residual = 0.582249, Final residual = 18.1122, No Iterations 1001
time step continuity errors : sum local = 2.2474e+41, global = -6.08311e+32, cumulative = -6.08312e+32
DICPCG: Solving for p, Initial residual = 0.940306, Final residual = 3.70065, No Iterations 1001
time step continuity errors : sum local = 1.08796e+43, global = 1.37468e+36, cumulative = 1.37407e+36
#0 Foam::error:rintStack(Foam::Ostream&) in "/apl/OpenFOAM/OpenFOAM-2.0.0/platforms/linux64Gcc45DPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/apl/OpenFOAM/OpenFOAM-2.0.0/platforms/linux64Gcc45DPOpt/lib/libOpenFOAM.so"
#2 __restore_rt at sigaction.c:0
#3 Foam::LimitedScheme<double, Foam::limitedLinearLimiter<Foam::NVDTVD>, Foam::limitFuncs::magSqr>::limiter(Foam::Geometric Field<double, Foam::fvPatchField, Foam::volMesh> const&) const in "/apl/OpenFOAM/OpenFOAM-2.0.0/platforms/linux64Gcc45DPOpt/lib/libfiniteVolume.so"
#4 Foam::limitedSurfaceInterpolationScheme<double>::w eights(Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) const in "/apl/OpenFOAM/OpenFOAM-2.0.0/platforms/linux64Gcc45DPOpt/lib/libfiniteVolume.so"
#5 Foam::fv::gaussConvectionScheme<double>::fvmDiv(Fo am::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) const in "/apl/OpenFOAM/OpenFOAM-2.0.0/platforms/linux64Gcc45DPOpt/lib/libfiniteVolume.so"
#6 Foam::tmp<Foam::fvMatrix<double> > Foam::fvm::div<double>(Foam::GeometricField<double , Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&, Foam::word const&) in "/apl/OpenFOAM/OpenFOAM-2.0.0/platforms/linux64Gcc45DPOpt/lib/libincompressibleRASModels.so"
#7 Foam::tmp<Foam::fvMatrix<double> > Foam::fvm::div<double>(Foam::GeometricField<double , Foam::fvsPatchField, Foam::surfaceMesh> const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> const&) in "/apl/OpenFOAM/OpenFOAM-2.0.0/platforms/linux64Gcc45DPOpt/lib/libincompressibleRASModels.so"
#8 Foam::incompressible::RASModels::kEpsilon::correct () in "/apl/OpenFOAM/OpenFOAM-2.0.0/platforms/linux64Gcc45DPOpt/lib/libincompressibleRASModels.so"
#9 main in "/apl/OpenFOAM/OpenFOAM-2.0.0/platforms/linux64Gcc45DPOpt/bin/pisoFoam"
#10 __libc_start_main in "/lib64/libc.so.6"
#11 Foam::UOPstream::write(char) in "/apl/OpenFOAM/OpenFOAM-2.0.0/platforms/linux64Gcc45DPOpt/bin/pisoFoam"
Floating exception


I've tried to decrease the deltaT but it doesn't work. I also changed the nNonOrthogonalCorrectors and the ddtSchemes

If I run the simulation without a turbulent model it works.

I tried to modify the boundaryField but i cannot find out the problem. I'm just learning to use openfoam so it may be a basic problem that i cannot see.

I couldn't upload the geometry files since it's too large.

Thanks for your help

Best regards
Attached Images
File Type: jpg cylinder.jpg (95.7 KB, 51 views)
Attached Files
File Type: zip Cylinder.zip (6.5 KB, 21 views)
tol24 is offline   Reply With Quote

Old   December 8, 2011, 04:53
Default
  #2
Senior Member
 
Suresh kumar Kannan
Join Date: Mar 2009
Location: Luxembourg, Luxembourg, Luxembourg
Posts: 129
Rep Power: 8
kumar is on a distinguished road
HI ,
Just a suggestion run a checkMesh and see if your aspect ratio and Non-orthogonality checks are ok.

It may also be caused due to the non-appropriate boundary conditions. So if you give some details about the boundary conditions.

May be some one can help.

regards
K.SUresh kumar
kumar is offline   Reply With Quote

Old   December 8, 2011, 05:32
Default
  #3
Member
 
Rob
Join Date: Sep 2011
Posts: 55
Rep Power: 5
robbirobocop is on a distinguished road
@kumar: He attached his case in a .zip file

@tol24:
By having a quick look at your files I assume that your BC's for the symmetryPlane are wrong. The wallfunctions might not be the problem. But having a look at your "U" file you can see that you assigned 17 m/s on your symmetryPlane. So even if you treat your symmetryPlane like a wall this is not correct!

If the "fix(es)" don't help maybe you should give simpleFoam a try.
robbirobocop is offline   Reply With Quote

Old   December 9, 2011, 07:57
Default
  #4
New Member
 
Join Date: Dec 2011
Posts: 2
Rep Power: 0
tol24 is on a distinguished road
Thanks, for helping!

@kumar: I did check mesh and it was ok

@robbirobocop: yes, I thought it could be that, but actually I want to simulate the cylinder in a freestream and not between two walls. How could I do that?
If I fix the velocity to zero it will generate a boundary layer arround the limits, doesn't it?
tol24 is offline   Reply With Quote

Old   August 19, 2013, 04:19
Default
  #5
Member
 
Xiangyin Meng
Join Date: Jun 2012
Location: London
Posts: 69
Rep Power: 5
mxylondon is on a distinguished road
I have the nearly exactly same trouble now.
mxylondon 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
mach number and selection of the fluent solver turbinesv FLUENT 4 April 24, 2011 00:14
Newbie Question IcoFoam - Courant Number explodes sprobst76 OpenFOAM Running, Solving & CFD 10 November 19, 2010 12:22
Courant number songpen1985 FLUENT 0 April 17, 2009 03:45
LES near wall model & courant number kasim CFX 5 March 16, 2008 19:23
Trimmed cell and embedded refinement mesh conversion issues michele OpenFOAM Other Meshers: ICEM, Star, Ansys, Pointwise, GridPro, Ansa, ... 2 July 15, 2005 04:15


All times are GMT -4. The time now is 23:33.