what is the meaning of this?!??!!?

adambarfi · October 5, 2012, 12:21

Hi everybody

I'm modeling a rotating square between water and air. when I run my model, the following messages appeared in my terminal:

Code:

solidBodyMotionFunctions::rotatingMotion::transformation(): Time = 0.0827406 transformation: ((0 0 0) (-0.856565 (0 0 0.51604)))
Execution time for mesh.update() = 0.05 s
time step continuity errors : sum local = 5.19359e-14, global = 1.91199e-15, cumulative = 1.45354e-06
GAMGPCG:  Solving for pcorr, Initial residual = 1, Final residual = 3.28882e-06, No Iterations 9
time step continuity errors : sum local = 3.27857e-15, global = -4.81172e-17, cumulative = 1.45354e-06
MULES: Solving for alpha1
Liquid phase volume fraction = 0.129117  Min(alpha1) = -3.3351e-36  Max(alpha1) = 0.999977
MULES: Solving for alpha1
Liquid phase volume fraction = 0.129169  Min(alpha1) = -3.33127e-36  Max(alpha1) = 0.999977
MULES: Solving for alpha1
Liquid phase volume fraction = 0.12922  Min(alpha1) = -3.32675e-36  Max(alpha1) = 0.999977
GAMG:  Solving for p_rgh, Initial residual = 0.0226012, Final residual = 0.000126703, No Iterations 3
time step continuity errors : sum local = 1.35142e-09, global = 7.04088e-12, cumulative = 1.45355e-06
GAMGPCG:  Solving for p_rgh, Initial residual = 0.00713841, Final residual = 5.44584e-09, No Iterations 6
time step continuity errors : sum local = 5.81134e-14, global = 2.97304e-16, cumulative = 1.45355e-06
smoothSolver:  Solving for epsilon, Initial residual = 0.00084478, Final residual = 7.54014e-08, No Iterations 3
smoothSolver:  Solving for k, Initial residual = 0.000782825, Final residual = 2.30447e-07, No Iterations 3
ExecutionTime = 305.2 s  ClockTime = 306 s

Interface Courant Number mean: 0.0379477 max: 0.490161
Courant Number mean: 0.248671 max: 0.499417
deltaT = 7.19143e-05
Time = 0.0828125

solidBodyMotionFunctions::rotatingMotion::transformation(): Time = 0.0828125 transformation: ((0 0 0) (-0.857728 (0 0 0.514103)))
Execution time for mesh.update() = 0.04 s
time step continuity errors : sum local = 5.86878e-14, global = 2.49004e-16, cumulative = 1.45355e-06

can anybody help me?
what is the meaning of this message:

Code:

time step continuity errors : sum local = 5.86878e-14, global = 2.49004e-16, cumulative = 1.45355e-06

thank you

ata · October 6, 2012, 07:34

Hi
I think
sum local is absolute summation of continuty error in all CVs at present time step.
global is summation of continuty error in all CVs at present time stem and
cumulative is summation of global continuity errors from the first time step.

adambarfi · October 6, 2012, 12:25

Quote:

Originally Posted by ata

Hi
I think
sum local is absolute summation of continuty error in all CVs at present time step.
global is summation of continuty error in all CVs at present time stem and
cumulative is summation of global continuity errors from the first time step.

thank you Ata for your reply.
but does this error has any influence on my results? I think the answer is: No, It doesn't.

alberto · October 6, 2012, 20:44

Yes, it does

The solver reports the error for informative purposes, as a measure of how good your solution is. You want this error to be as close to machine precision as possible.

adambarfi · October 8, 2012, 02:46

Quote:

Originally Posted by alberto

Yes, it does

The solver reports the error for informative purposes, as a measure of how good your solution is. You want this error to be as close to machine precision as possible.

thank you Alberto,

but what do you think about the source of below error:

Code:

DICPCG:  Solving for p, Initial residual = 1, Final residual = 6.72375e-08, No Iterations 62
time step continuity errors : sum local = 4.77002e+25, global = -9.93552e+15, cumulative = -9.93552e+15
Model mode 1
#0  Foam::error::printStack(Foam::Ostream&) in "/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#1  Foam::sigFpe::sigHandler(int) in "/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#2  Uninterpreted: 
#3  Foam::PBiCG::solve(Foam::Field<double>&, Foam::Field<double> const&, unsigned char) const in "/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/platforms/linuxGccDPOpt/lib/libOpenFOAM.so"
#4  Foam::fvMatrix<Foam::SymmTensor<double> >::solve(Foam::dictionary const&) in "/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/platforms/linuxGccDPOpt/lib/libviscoelasticTransportModels.so"
#5  Foam::fvMatrix<Foam::SymmTensor<double> >::solve() in "/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/platforms/linuxGccDPOpt/lib/libviscoelasticTransportModels.so"
#6  Foam::Giesekus::correct() in "/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/platforms/linuxGccDPOpt/lib/libviscoelasticTransportModels.so"
#7  Foam::multiMode::correct() in "/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/platforms/linuxGccDPOpt/lib/libviscoelasticTransportModels.so"
#8  Foam::viscoelasticModel::correct() in "/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/platforms/linuxGccDPOpt/lib/libviscoelasticTransportModels.so"
#9  
 in "/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/platforms/linuxGccDPOpt/bin/BuoyantBoussinesqViscoelasticFluidFoam"
#10  __libc_start_main in "/lib/i386-linux-gnu/libc.so.6"
#11  
 in "/home/mostafa/OpenFOAM/OpenFOAM-2.1.0/platforms/linuxGccDPOpt/bin/BuoyantBoussinesqViscoelasticFluidFoam"
Floating point exception

I'm so confused. please help me.
my solver is the viscoelasticFluidFoam that I made some changes to it for solving natural convection. when I ran it, after some iterate the above error appeared.

anyone is here can help me????

thanks.

ata · October 8, 2012, 05:23

Hi
If I want to speak general your schemes and/or methods for solving system of linear equations and/or their adjustments are improper.

adambarfi · October 8, 2012, 06:12

Quote:

Originally Posted by ata

Hi
If I want to speak general your schemes and/or methods for solving system of linear equations and/or their adjustments are improper.

hi ata,.
can you tell me the details? I'm about two weeks that staying behind this error.

thank you.

EricKl. · September 12, 2013, 15:24

Hello, I am trying to simulate a runoff, but I'm having the following problem:
What can I do?
Thanks

Code:

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM Extend Project: Open source CFD        |
|  \\    /   O peration     | Version:  1.6-ext                               |
|   \\  /    A nd           | Web:      www.extend-project.de                 |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 1.6-ext-e61a8f0a3194
Exec   : viscoelasticInterFoam
Date   : Sep 12 2013
Time   : 16:23:40
Host   : administrador-desktop
PID    : 10189
Case   : /home/administrador/tutorials/viscoelastic/VIFErick/Escoamento
nProcs : 1
SigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).

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

Create mesh for time = 0

// using new solver syntax:
U
{
    solver          PBiCG;
    preconditioner
    {
        preconditioner  DILU;
    }
    minIter         0;
    maxIter         1000;
    tolerance       1e-06;
    relTol          0;
}

// using new solver syntax:
taufirst
{
    solver          BiCGStab;
    preconditioner
    {
        preconditioner  Cholesky;
    }
    minIter         0;
    maxIter         1000;
    tolerance       1e-06;
    relTol          0;
}

// using new solver syntax:
tau
{
    solver          BiCGStab;
    preconditioner
    {
        preconditioner  Cholesky;
    }
    minIter         0;
    maxIter         1000;
    tolerance       1e-06;
    relTol          0;
}


Reading g

Reading transportProperties

Reading field pd

Reading field alpha

Reading field U

Reading/calculating face flux field phi

Calculating field g.h

Selecting viscoelastic model Giesekus



time step continuity errors : sum local = 6.31081e-07, global = -6.31081e-07, cumulative = -6.31081e-07


--> FOAM FATAL ERROR: 
Continuity error cannot be removed by adjusting the outflow.
Please check the velocity boundary conditions and/or run potentialFoam to initialise the outflow.
Specified mass inflow   : 7.65133e-08
Specified mass outflow  : 0
Difference              : 7.65133e-08
Adjustable mass outflow : 0


    From function adjustPhi
(
    surfaceScalarField& phi,
    const volVectorField& U,
    const volScalarField& p
)
    in file cfdTools/general/adjustPhi/adjustPhi.C at line 141.

FOAM exiting

Bernhard · September 13, 2013, 02:39

Please check the velocity boundary conditions. What are the velocity boundary conditions in your case?

EricKl. · September 13, 2013, 10:53

Hi,
Thanks for reply

My code for velocity is:

Code:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM Extend Project: Open Source CFD        |
|  \\    /   O peration     | Version:  1.6-ext                               |
|   \\  /    A nd           | Web:      www.extend-project.de                 |
|    \\/     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
{
    inlet
    {
        type            fixedValue;
        value           uniform (0.01868 0 0);
    }
    fixedWalls
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }
    frontAndBack
    {
        type            empty;
    }
    outlet
    {
        type            zeroGradient;
    }
    simetry
    {
        type            symmetryPlane;
    }
    atmosphere
    {
        type            pressureInletOutletVelocity;
        phi             phi;
        value           uniform (0 0 0);
    }
   
}

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

When I run potentialFoam, the program gives me

Code:

/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM Extend Project: Open source CFD        |
|  \\    /   O peration     | Version:  1.6-ext                               |
|   \\  /    A nd           | Web:      www.extend-project.de                 |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 1.6-ext-e61a8f0a3194
Exec   : potentialFoam
Date   : Sep 13 2013
Time   : 11:55:07
Host   : administrador-desktop
PID    : 16452
Case   : /home/administrador/tutorials/viscoelastic/VIFErick/Escoamento
nProcs : 1
SigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).

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

Create mesh for time = 0

// using new solver syntax:
U
{
    solver          PBiCG;
    preconditioner
    {
        preconditioner  DILU;
    }
    minIter         0;
    maxIter         1000;
    tolerance       1e-06;
    relTol          0;
}

// using new solver syntax:
taufirst
{
    solver          BiCGStab;
    preconditioner
    {
        preconditioner  Cholesky;
    }
    minIter         0;
    maxIter         1000;
    tolerance       1e-06;
    relTol          0;
}

// using new solver syntax:
tau
{
    solver          BiCGStab;
    preconditioner
    {
        preconditioner  Cholesky;
    }
    minIter         0;
    maxIter         1000;
    tolerance       1e-06;
    relTol          0;
}

Reading field p

Reading field U



--> FOAM FATAL IO ERROR: 
keyword SIMPLE is undefined in dictionary "/home/administrador/tutorials/viscoelastic/VIFErick/Escoamento/system/fvSolution"

file: /home/administrador/tutorials/viscoelastic/VIFErick/Escoamento/system/fvSolution at line 21.

    From function dictionary::subDict(const word& keyword) const
    in file db/dictionary/dictionary.C at line 457.

FOAM exiting

I tried move and ajust velocity, but the errors continues
I'm really lost

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
What's meaning of UDF FUNCTION	zhaoxinyu	Fluent UDF and Scheme Programming	0	March 31, 2010 08:04
scalar meaning	usker	Siemens	0	July 24, 2008 03:49
scalar meaning	usker	Siemens	3	July 6, 2008 11:26
want to know meaning	Sangamesh	Siemens	0	May 15, 2007 05:15
What's the meaning of "combustion scalar"and....	cfdbeginner	CFX	0	November 27, 2003 09:02

October 6, 2012, 07:34		#2
ata Senior Member ata kamyabi Join Date: Aug 2009 Location: Kerman Posts: 323 Rep Power: 17	Hi I think sum local is absolute summation of continuty error in all CVs at present time step. global is summation of continuty error in all CVs at present time stem and cumulative is summation of global continuity errors from the first time step. adambarfi, yesaswi92 and granzer like this.

October 6, 2012, 20:44		#4
alberto Senior Member Alberto Passalacqua Join Date: Mar 2009 Location: Ames, Iowa, United States Posts: 1,912 Rep Power: 36	Yes, it does The solver reports the error for informative purposes, as a measure of how good your solution is. You want this error to be as close to machine precision as possible. adambarfi and yesaswi92 like this. __________________ Alberto Passalacqua GeekoCFD - A free distribution based on openSUSE 64 bit with CFD tools, including OpenFOAM. Available as in both physical and virtual formats (current status: http://albertopassalacqua.com/?p=1541) OpenQBMM - An open-source implementation of quadrature-based moment methods. To obtain more accurate answers, please specify the version of OpenFOAM you are using.

October 8, 2012, 05:23		#6
ata Senior Member ata kamyabi Join Date: Aug 2009 Location: Kerman Posts: 323 Rep Power: 17	Hi If I want to speak general your schemes and/or methods for solving system of linear equations and/or their adjustments are improper.

September 13, 2013, 02:39		#9
Bernhard Senior Member Bernhard Join Date: Sep 2009 Location: Delft Posts: 790 Rep Power: 21	Please check the velocity boundary conditions. What are the velocity boundary conditions in your case?