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

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

Register Blogs Members List Search Today's Posts Mark Forums Read

Like Tree4Likes
  • 2 Post By ata
  • 2 Post By alberto

Reply
 
LinkBack Thread Tools Display Modes
Old   October 5, 2012, 12:21
Default what is the meaning of this?!??!!?
  #1
Senior Member
 
adambarfi's Avatar
 
Mostafa Mahmoudi
Join Date: Jan 2012
Posts: 314
Rep Power: 6
adambarfi is on a distinguished road
Send a message via Yahoo to adambarfi Send a message via Skype™ to adambarfi
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
adambarfi is offline   Reply With Quote

Old   October 6, 2012, 07:34
Default
  #2
ata
Senior Member
 
ata's Avatar
 
ata kamyabi
Join Date: Aug 2009
Location: Kerman
Posts: 322
Rep Power: 8
ata is on a distinguished road
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 and yesaswi92 like this.
ata is offline   Reply With Quote

Old   October 6, 2012, 12:25
Default
  #3
Senior Member
 
adambarfi's Avatar
 
Mostafa Mahmoudi
Join Date: Jan 2012
Posts: 314
Rep Power: 6
adambarfi is on a distinguished road
Send a message via Yahoo to adambarfi Send a message via Skype™ to adambarfi
Quote:
Originally Posted by ata View Post
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.
adambarfi is offline   Reply With Quote

Old   October 6, 2012, 20:44
Default
  #4
Senior Member
 
Alberto Passalacqua
Join Date: Mar 2009
Location: Ames, Iowa, United States
Posts: 1,894
Rep Power: 26
alberto will become famous soon enoughalberto will become famous soon enough
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 live DVD/USB, hard drive image and virtual image.
OpenQBMM - An open-source implementation of quadrature-based moment methods
alberto is offline   Reply With Quote

Old   October 8, 2012, 02:46
Default
  #5
Senior Member
 
adambarfi's Avatar
 
Mostafa Mahmoudi
Join Date: Jan 2012
Posts: 314
Rep Power: 6
adambarfi is on a distinguished road
Send a message via Yahoo to adambarfi Send a message via Skype™ to adambarfi
Quote:
Originally Posted by alberto View Post
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.
adambarfi is offline   Reply With Quote

Old   October 8, 2012, 05:23
Default
  #6
ata
Senior Member
 
ata's Avatar
 
ata kamyabi
Join Date: Aug 2009
Location: Kerman
Posts: 322
Rep Power: 8
ata is on a distinguished road
Hi
If I want to speak general your schemes and/or methods for solving system of linear equations and/or their adjustments are improper.
ata is offline   Reply With Quote

Old   October 8, 2012, 06:12
Default
  #7
Senior Member
 
adambarfi's Avatar
 
Mostafa Mahmoudi
Join Date: Jan 2012
Posts: 314
Rep Power: 6
adambarfi is on a distinguished road
Send a message via Yahoo to adambarfi Send a message via Skype™ to adambarfi
Quote:
Originally Posted by ata View Post
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.
adambarfi is offline   Reply With Quote

Old   September 12, 2013, 15:24
Default
  #8
New Member
 
Erick Lorenzato
Join Date: Mar 2013
Posts: 8
Rep Power: 4
EricKl. is on a distinguished road
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
EricKl. is offline   Reply With Quote

Old   September 13, 2013, 02:39
Default
  #9
Senior Member
 
Bernhard
Join Date: Sep 2009
Location: Delft
Posts: 790
Rep Power: 12
Bernhard is on a distinguished road
Please check the velocity boundary conditions. What are the velocity boundary conditions in your case?
Bernhard is offline   Reply With Quote

Old   September 13, 2013, 10:53
Default
  #10
New Member
 
Erick Lorenzato
Join Date: Mar 2013
Posts: 8
Rep Power: 4
EricKl. is on a distinguished road
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
EricKl. is offline   Reply With Quote

Reply

Thread Tools
Display Modes

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 On
Pingbacks are On
Refbacks are On


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 CD-adapco 0 July 24, 2008 03:49
scalar meaning usker CD-adapco 3 July 6, 2008 11:26
want to know meaning Sangamesh CD-adapco 0 May 15, 2007 05:15
What's the meaning of "combustion scalar"and.... cfdbeginner CFX 0 November 27, 2003 10:02


All times are GMT -4. The time now is 12:54.