How to stop courant number rise more than 1 in icoFoam model

1,149 views
Skip to first unread message

Prasanna Welahetti

unread,
Feb 2, 2016, 1:20:59 PM2/2/16
to OpenFOAM
Haloo,


When I run the openFoam icoFoam in command window i got following error.

How to stop courant number rise more than 1.





meeseeks@meeseeks-Precision-Tower-7810:~/prasa/run/model_22$ icoFoam
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.4.0                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 2.4.0-f0842aea0e77
Exec   : icoFoam
Date   : Feb 02 2016
Time   : 18:32:30
Host   : "meeseeks-Precision-Tower-7810"
PID    : 12783
Case   : /home/meeseeks/prasa/run/model_22
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0

Reading transportProperties

Reading field p

Reading field U

Reading/calculating face flux field phi


Starting time loop

Time = 0.005

Courant Number mean: 0.00614036 max: 1.35776
smoothSolver:  Solving for Ux, Initial residual = 1, Final residual = 9.80075e-06, No Iterations 24
smoothSolver:  Solving for Uy, Initial residual = 0, Final residual = 0, No Iterations 0
smoothSolver:  Solving for Uz, Initial residual = 0, Final residual = 0, No Iterations 0
DICPCG:  Solving for p, Initial residual = 1, Final residual = 8.02393e-07, No Iterations 168
time step continuity errors : sum local = 1.12673e-08, global = -1.65632e-10, cumulative = -1.65632e-10
DICPCG:  Solving for p, Initial residual = 0.0113109, Final residual = 9.14177e-07, No Iterations 139
time step continuity errors : sum local = 2.55703e-06, global = 1.5327e-08, cumulative = 1.51614e-08
ExecutionTime = 0.77 s  ClockTime = 1 s

Time = 0.01

Courant Number mean: 1.52506 max: 9.1826
#0  Foam::error::printStack(Foam::Ostream&) at ??:?
#1  Foam::sigFpe::sigHandler(int) at ??:?
#2  ? in "/lib/x86_64-linux-gnu/libc.so.6"
#3  Foam::symGaussSeidelSmoother::smooth(Foam::word const&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::Field<double> const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, unsigned char, int) at ??:?
#4  Foam::symGaussSeidelSmoother::smooth(Foam::Field<double>&, Foam::Field<double> const&, unsigned char, int) const at ??:?
#5  Foam::smoothSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#6  ? at ??:?
#7  ? at ??:?
#8  ? at ??:?
#9  ? at ??:?
#10  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#11  ? at ??:?
Floating point exception (core dumped)
meeseeks@meeseeks-Precision-Tower-7810:~/prasa/run/model_22$




Best regards,
Prasanna.

ohbuchi

unread,
Feb 2, 2016, 7:30:27 PM2/2/16
to OpenFOAM
Hi,

I can not determine the reason of this phenomena, because your information is insufficient.
At least, following information will be needed.
1. mesh (blockMeshDict)
2. boundary and initial conditions (boundary, 0/U,p)
3. numerical scheme, solver settings (fvSchemes, fvSolution)

Best,
M.Ohbuchi


2016年2月3日水曜日 3時20分59秒 UTC+9 Prasanna Welahetti:

Prasanna Welahetti

unread,
Feb 4, 2016, 6:27:19 AM2/4/16
to open...@googlegroups.com
Hai,

Now I have this issue corant number rising continuesly give error massege, I'm using icoFoam solver. following error massege comes, what can be the reasons ?


/*---------------------------------------------------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.4.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
Build : 2.4.0-f0842aea0e77
Exec : icoFoam
Date : Feb 04 2016
Time : 12:11:05
Host : "meeseeks-Precision-Tower-7810"
PID : 3746
Case : /home/meeseeks/prasa/run/model3

nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0

Reading transportProperties

Reading field p

Reading field U

Reading/calculating face flux field phi


Starting time loop

Time = 5e-05

Courant Number mean: 6.71136e-05 max: 0.0518149
smoothSolver: Solving for Ux, Initial residual = 1, Final residual = 1.38584e-06, No Iterations 3
smoothSolver: Solving for Uy, Initial residual = 1, Final residual = 5.94309e-06, No Iterations 3

smoothSolver: Solving for Uz, Initial residual = 0, Final residual = 0, No Iterations 0
DICPCG: Solving for p, Initial residual = 1, Final residual = 7.70462e-07, No Iterations 164
time step continuity errors : sum local = 1.03434e-10, global = 2.72659e-12, cumulative = 2.72659e-12
DICPCG: Solving for p, Initial residual = 0.00784595, Final residual = 9.57557e-07, No Iterations 124
time step continuity errors : sum local = 2.06591e-08, global = 2.60165e-10, cumulative = 2.62891e-10
ExecutionTime = 0.73 s ClockTime = 1 s

Time = 0.0001

Courant Number mean: 0.0151534 max: 0.374485
smoothSolver: Solving for Ux, Initial residual = 0.606285, Final residual = 9.33117e-07, No Iterations 3
smoothSolver: Solving for Uy, Initial residual = 0.13063, Final residual = 7.57316e-06, No Iterations 2
smoothSolver: Solving for Uz, Initial residual = 0.332338, Final residual = 8.93221e-06, No Iterations 2
DICPCG: Solving for p, Initial residual = 0.142013, Final residual = 8.82601e-07, No Iterations 143
time step continuity errors : sum local = 2.16199e-08, global = -6.83841e-11, cumulative = 1.94507e-10
DICPCG: Solving for p, Initial residual = 0.0453243, Final residual = 8.11834e-07, No Iterations 131
time step continuity errors : sum local = 7.20617e-09, global = -2.82476e-11, cumulative = 1.6626e-10
ExecutionTime = 1.08 s ClockTime = 1 s

Time = 0.00015

Courant Number mean: 0.0158302 max: 0.561537
smoothSolver: Solving for Ux, Initial residual = 0.181653, Final residual = 6.75105e-07, No Iterations 3
smoothSolver: Solving for Uy, Initial residual = 0.0507999, Final residual = 2.53766e-06, No Iterations 2
smoothSolver: Solving for Uz, Initial residual = 0.249424, Final residual = 8.27343e-07, No Iterations 3
DICPCG: Solving for p, Initial residual = 0.522183, Final residual = 9.40058e-07, No Iterations 140
time step continuity errors : sum local = 5.53735e-09, global = 9.29897e-12, cumulative = 1.75559e-10
DICPCG: Solving for p, Initial residual = 0.10772, Final residual = 7.87431e-07, No Iterations 128
time step continuity errors : sum local = 1.57164e-09, global = 2.80607e-11, cumulative = 2.03619e-10
ExecutionTime = 1.42 s ClockTime = 1 s

Time = 0.0002

Courant Number mean: 0.0157345 max: 0.557623
smoothSolver: Solving for Ux, Initial residual = 0.0260866, Final residual = 3.24105e-06, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 0.00990692, Final residual = 6.46839e-07, No Iterations 2
smoothSolver: Solving for Uz, Initial residual = 0.0538437, Final residual = 2.93782e-06, No Iterations 2
DICPCG: Solving for p, Initial residual = 0.349896, Final residual = 9.11121e-07, No Iterations 139
time step continuity errors : sum local = 1.33925e-09, global = -1.13943e-11, cumulative = 1.92225e-10
DICPCG: Solving for p, Initial residual = 0.0431763, Final residual = 8.98031e-07, No Iterations 132
time step continuity errors : sum local = 9.07484e-10, global = -2.82289e-11, cumulative = 1.63996e-10
ExecutionTime = 1.76 s ClockTime = 2 s

Time = 0.00025

Courant Number mean: 0.0157019 max: 0.550726
smoothSolver: Solving for Ux, Initial residual = 0.0120567, Final residual = 2.41214e-06, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 0.00481902, Final residual = 6.04422e-07, No Iterations 2
smoothSolver: Solving for Uz, Initial residual = 0.0239811, Final residual = 9.86613e-07, No Iterations 2
DICPCG: Solving for p, Initial residual = 0.190656, Final residual = 8.03269e-07, No Iterations 136
time step continuity errors : sum local = 7.23939e-10, global = -1.28889e-13, cumulative = 1.63867e-10
DICPCG: Solving for p, Initial residual = 0.0155782, Final residual = 9.49934e-07, No Iterations 125
time step continuity errors : sum local = 7.76001e-10, global = -1.30238e-11, cumulative = 1.50843e-10
ExecutionTime = 2.08 s ClockTime = 2 s

Time = 0.0003

Courant Number mean: 0.0156726 max: 0.654244
smoothSolver: Solving for Ux, Initial residual = 0.00879693, Final residual = 1.94545e-06, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 0.00383515, Final residual = 6.77615e-07, No Iterations 2
smoothSolver: Solving for Uz, Initial residual = 0.0171416, Final residual = 6.29931e-07, No Iterations 2
DICPCG: Solving for p, Initial residual = 0.116409, Final residual = 8.701e-07, No Iterations 134
time step continuity errors : sum local = 6.62414e-10, global = 2.8479e-13, cumulative = 1.51128e-10
DICPCG: Solving for p, Initial residual = 0.00560623, Final residual = 9.80388e-07, No Iterations 122
time step continuity errors : sum local = 7.07401e-10, global = 5.86687e-12, cumulative = 1.56995e-10
ExecutionTime = 2.4 s ClockTime = 2 s

Time = 0.00035

Courant Number mean: 0.0156595 max: 0.870647
smoothSolver: Solving for Ux, Initial residual = 0.00727937, Final residual = 1.73164e-06, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 0.00348196, Final residual = 1.55832e-06, No Iterations 2
smoothSolver: Solving for Uz, Initial residual = 0.0140972, Final residual = 7.85586e-07, No Iterations 2
DICPCG: Solving for p, Initial residual = 0.0964285, Final residual = 7.90232e-07, No Iterations 135
time step continuity errors : sum local = 5.6033e-10, global = 5.53276e-12, cumulative = 1.62528e-10
DICPCG: Solving for p, Initial residual = 0.00281709, Final residual = 9.23352e-07, No Iterations 96
time step continuity errors : sum local = 6.50587e-10, global = -3.56866e-11, cumulative = 1.26841e-10
ExecutionTime = 2.69 s ClockTime = 3 s

Time = 0.0004

Courant Number mean: 0.015657 max: 1.19054
smoothSolver: Solving for Ux, Initial residual = 0.00658232, Final residual = 2.7559e-06, No Iterations 2
smoothSolver: Solving for Uy, Initial residual = 0.00352735, Final residual = 3.16808e-06, No Iterations 3
smoothSolver: Solving for Uz, Initial residual = 0.0127143, Final residual = 6.00431e-06, No Iterations 2
DICPCG: Solving for p, Initial residual = 0.111259, Final residual = 9.04817e-07, No Iterations 138
time step continuity errors : sum local = 6.60456e-10, global = 2.80977e-11, cumulative = 1.54939e-10
DICPCG: Solving for p, Initial residual = 0.00427442, Final residual = 7.9171e-07, No Iterations 127
time step continuity errors : sum local = 6.10697e-10, global = -2.55804e-12, cumulative = 1.52381e-10
ExecutionTime = 3.02 s ClockTime = 3 s

Time = 0.00045

Courant Number mean: 0.0156727 max: 1.75486
smoothSolver: Solving for Ux, Initial residual = 0.0066215, Final residual = 8.75173e-06, No Iterations 19
smoothSolver: Solving for Uy, Initial residual = 0.00425823, Final residual = 8.65799e-06, No Iterations 36
smoothSolver: Solving for Uz, Initial residual = 0.0124809, Final residual = 8.67935e-06, No Iterations 32
DICPCG: Solving for p, Initial residual = 0.166059, Final residual = 9.73893e-07, No Iterations 143
time step continuity errors : sum local = 8.89965e-10, global = 4.2187e-11, cumulative = 1.94568e-10
DICPCG: Solving for p, Initial residual = 0.0170893, Final residual = 9.57808e-07, No Iterations 137
time step continuity errors : sum local = 1.06481e-09, global = -2.37156e-11, cumulative = 1.70852e-10
ExecutionTime = 3.5 s ClockTime = 3 s

Time = 0.0005

Courant Number mean: 0.0157411 max: 2.85948
#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::symGaussSeidelSmoother::smooth(Foam::word const&, Foam::Field<double>&, Foam::lduMatrix const&, Foam::Field<double> const&, Foam::FieldField<Foam::Field, double> const&, Foam::UPtrList<Foam::lduInterfaceField const> const&, unsigned char, int) at ??:?
#4 Foam::symGaussSeidelSmoother::smooth(Foam::Field<d ouble>&, Foam::Field<double> const&, unsigned char, int) const at ??:?

#5 Foam::smoothSolver::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const at ??:?
#6 ? at ??:?
#7 ? at ??:?
#8 ? at ??:?
#9 ? at ??:?
#10 __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#11 ? at ??:?
Floating point exception (core dumped)
meeseeks@meeseeks-Precision-Tower-7810:~/prasa/run/model3$ ^C
meeseeks@meeseeks-Precision-Tower-7810:~/prasa/run/model3$



Please find below details,

P

/*--------------------------------*- C++ -*----------------------------------*\

| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.4.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volScalarField;
object p;

}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 2 -2 0 0 0 0];

internalField uniform 0;

boundaryField
{
pressure_outlet
{
type fixedValue;
value uniform 0;

}

velocity_inlet_small
{
type zeroGradient;
}

velocity_inlet_large
{
type zeroGradient;
}
walls
{
type zeroGradient;
}
}


// ************************************************** *********************** //

U dictionary

/*--------------------------------*- C++ -*----------------------------------*\

| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.4.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class volVectorField;
object U;

}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions [0 1 -1 0 0 0 0];

internalField uniform (0 0 0);

boundaryField
{
pressure_outlet
{
type zeroGradient;

}

velocity_inlet_small
{
type fixedValue;
value uniform (0 -15 0);
}

velocity_inlet_large
{
type fixedValue;
value uniform (10 0 0);
}
walls
{
type fixedValue;
value uniform (0 0 0);
}

}

// ************************************************** *********************** //

controldict

/*--------------------------------*- C++ -*----------------------------------*\

| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.4.0 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ M anipulation | |
\*---------------------------------------------------------------------------*/
FoamFile
{
version 2.0;
format ascii;
class dictionary;
location "system";
object controlDict;

}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

application icoFoam;

startFrom startTime;

startTime 0;

stopAt endTime;

endTime 1;

deltaT 0.00005;

writeControl adjustableRunTime;

writeInterval 0.01;

purgeWrite 0;

writeFormat ascii;

writePrecision 6;

writeCompression off;

timeFormat general;

timePrecision 6;

runTimeModifiable true;

adjustTimeStep yes;

maxCo 1;
maxAlphaCo 1;
maxDeltaT 1;



// ************************************************** *********************** //


Regards,
Prasanna.

--
このメールは Google グループのグループ「OpenFOAM」のトピックを登録しているユーザーに送られています。
このトピックの登録を解除するには https://groups.google.com/d/topic/openfoam/0fjCu1ah2jA/unsubscribe にアクセスしてください。
このグループを退会し、グループのすべてのトピックの登録を解除するには openfoam+u...@googlegroups.com にメールを送信してください。
このグループに投稿するには open...@googlegroups.com にメールを送信してください。
https://groups.google.com/group/openfoam からこのグループにアクセスしてください。
その他のオプションについては https://groups.google.com/d/optout にアクセスしてください。

大渕真志

unread,
Feb 4, 2016, 2:45:03 PM2/4/16
to OpenFOAM
Hi, Prasanna,

Did you check your mesh with checkMesh?
And icoFoam does not support adjustable time step feature.
I recommend using pisoFoam.

Best Regards,
M.Ohbuchi


--
このメールは Google グループのグループ「OpenFOAM」に登録しているユーザーに送られています。
このグループから退会し、グループからのメールの配信を停止するには openfoam+u...@googlegroups.com にメールを送信してください。
Reply all
Reply to author
Forward
0 new messages