Checking courant number
 

Hi, Foamers.:)

I have some problem using OpenFOAM.
My calculation blow up everytime I run application.

I think courant number is main reason....:(
block Mesh domain seems OK. because calculation doesn't blow up with domain block itself.
my surface mesh file in the domain seems have a big courant number...

Is there any method for checking courant number before running application?

I use PointWise to make surface mesh...

I wish someone could help me;)

Please read the rules and post appropiate information.

Solver used, error log, mesh description, etc.

sorry..
 

sorry for that mistakes..:(

here is my blockMesh, and surfaceMesh and control dictionary
I used interDyMFoam to solve my problem

Interface Courant Number mean: 9545.13 max: 4.27743e+08
Courant Number mean: 4.15919e+09 max: 3.45743e+14
deltaT = 1.21214e-87
--> FOAM Warning :
From function Time::operator++()
in file db/Time/Time.C at line 1055
Increased the timePrecision from 13 to 14 to distinguish between timeNames at time 0.00080893
Time = 0.00080893041356983


Restraint translationDamper: force (-0 -0 -5.17534e+35)
Restraint rotationDamper: moment (-0 -2.8171e+34 -0)
Centre of mass: (2.92954 0 0.0226606)
Linear velocity: (0 0 -3.38829e+66)
Angular velocity: (0 -1.45672e+66 0)
Execution time for mesh.update() = 0.42 s
GAMG: Solving for pcorr, Initial residual = 1, Final residual = 0.277579, No Iterations 1000
time step continuity errors : sum local = 0.000112955, global = 1.85283e-07, cumulative = 2.7206e+32
smoothSolver: Solving for alpha.water, Initial residual = 5.61573e-08, Final residual = 1.31736e-14, No Iterations 1
Phase-1 volume fraction = -2.83525e+27 Min(alpha1) = -3.6172e+46 Max(alpha1) = 3.18976e+46
Applying the previous iteration compression flux
MULES: Correcting alpha.water
MULES: Correcting alpha.water
MULES: Correcting alpha.water
MULES: Correcting alpha.water
#0 Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam230/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1 Foam::sigFpe::sigHandler(int) in "/opt/openfoam230/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2 in "/lib/x86_64-linux-gnu/libc.so.6"
#3 Foam::tmp<Foam::GeometricField<double, Foam::fvsPatchField, Foam::surfaceMesh> > Foam::mag<Foam::Vector<double>, Foam::fvsPatchField, Foam::surfaceMesh>(Foam::GeometricField<Foam::Vect or<double>, Foam::fvsPatchField, Foam::surfaceMesh> const&) in "/opt/openfoam230/platforms/linux64GccDPOpt/lib/libinterfaceProperties.so"
#4 Foam::interfaceProperties::calculateK() in "/opt/openfoam230/platforms/linux64GccDPOpt/lib/libinterfaceProperties.so"
#5
in "/opt/openfoam230/platforms/linux64GccDPOpt/bin/interDyMFoam"
#6 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#7
in "/opt/openfoam230/platforms/linux64GccDPOpt/bin/interDyMFoam"
Floating point exception (core dumped)



Here is my error log!

oops, here is my control dicts, blockMesh
if other conditions are OK, then how should I correct my surface Mesh?
make it sparse? or more separated?