|
[Sponsors] |
November 2, 2014, 09:43 |
Checking courant number
|
#1 |
New Member
Jason
Join Date: Oct 2014
Posts: 17
Rep Power: 12 |
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 |
|
November 2, 2014, 11:50 |
|
#2 |
Senior Member
anonymous
Join Date: Aug 2014
Posts: 205
Rep Power: 13 |
Please read the rules and post appropiate information.
Solver used, error log, mesh description, etc. |
|
November 3, 2014, 06:11 |
sorry..
|
#3 |
New Member
Jason
Join Date: Oct 2014
Posts: 17
Rep Power: 12 |
sorry for that mistakes..
here is my blockMesh, and surfaceMesh and control dictionary I used interDyMFoam to solve my problem |
|
November 3, 2014, 06:12 |
|
#4 |
New Member
Jason
Join Date: Oct 2014
Posts: 17
Rep Power: 12 |
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! |
|
November 3, 2014, 07:55 |
|
#5 |
New Member
Jason
Join Date: Oct 2014
Posts: 17
Rep Power: 12 |
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? |
|
|
|
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 |