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

Checking courant number

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   November 2, 2014, 09:43
Default Checking courant number
  #1
New Member
 
Jason
Join Date: Oct 2014
Posts: 17
Rep Power: 12
J.Y.Won is on a distinguished road
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
J.Y.Won is offline   Reply With Quote

Old   November 2, 2014, 11:50
Default
  #2
Senior Member
 
anonymous
Join Date: Aug 2014
Posts: 205
Rep Power: 13
ssss is on a distinguished road
Please read the rules and post appropiate information.

Solver used, error log, mesh description, etc.
ssss is offline   Reply With Quote

Old   November 3, 2014, 06:11
Default sorry..
  #3
New Member
 
Jason
Join Date: Oct 2014
Posts: 17
Rep Power: 12
J.Y.Won is on a distinguished road
sorry for that mistakes..

here is my blockMesh, and surfaceMesh and control dictionary
I used interDyMFoam to solve my problem
J.Y.Won is offline   Reply With Quote

Old   November 3, 2014, 06:12
Default
  #4
New Member
 
Jason
Join Date: Oct 2014
Posts: 17
Rep Power: 12
J.Y.Won is on a distinguished road
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:perator++()
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:rintStack(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!
J.Y.Won is offline   Reply With Quote

Old   November 3, 2014, 07:55
Default
  #5
New Member
 
Jason
Join Date: Oct 2014
Posts: 17
Rep Power: 12
J.Y.Won is on a distinguished road
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?
Attached Files
File Type: gz 0.org.tar.gz (1.2 KB, 8 views)
File Type: gz blockMesh_controldict_etc.tar.gz (1.4 KB, 8 views)
J.Y.Won is offline   Reply With Quote

Reply


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
[mesh manipulation] Mesh Refinement Luiz Eduardo Bittencourt Sampaio (Sampaio) OpenFOAM Meshing & Mesh Conversion 42 January 8, 2017 13:55
DecomposePar unequal number of shared faces maka OpenFOAM Pre-Processing 6 August 12, 2010 10:01
critical error during installation of openfoam Fabio88 OpenFOAM Installation 21 June 2, 2010 04:01
IcoFoam parallel woes msrinath80 OpenFOAM Running, Solving & CFD 9 July 22, 2007 03:58
Could anybody help me see this error and give help liugx212 OpenFOAM Running, Solving & CFD 3 January 4, 2006 19:07


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