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

Error using simpleFoam : Foam::error::printStack(Foam::Ostream&) at

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   May 16, 2017, 10:59
Default Error using simpleFoam : Foam::error::printStack(Foam::Ostream&) at
  #1
New Member
 
Kim Le
Join Date: May 2016
Location: Lyon, France
Posts: 21
Rep Power: 6
lethu is on a distinguished road
Hello,

I got this error message unsing simpleFoam but I can't understand where the error says it is from...

Can you help me please ?

Thanks !

Code:
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  4.1                                   |
|   \\  /    A nd           | Web:      www.OpenFOAM.org                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 4.1
Exec   : simpleFoam
Date   : May 16 2017
Time   : 14:52:58
Host   : "abc1b7de1187"
PID    : 827
Case   : /home/openfoam/run/22-bassinTestmixerVessel2D
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


SIMPLE: no convergence criteria found. Calculations will run for 500 steps.

Reading field p

Reading field U

Reading/calculating face flux field phi

Selecting incompressible transport model Newtonian
Selecting turbulence model type RAS
Selecting RAS turbulence model kEpsilon
bounding k, min: 0 max: 1 average: 1
bounding epsilon, min: 0 max: 0 average: 0
kEpsilonCoeffs
{
    Cmu             0.09;
    C1              1.44;
    C2              1.92;
    C3              -0.33;
    sigmak          1;
    sigmaEps        1.3;
}

Creating MRF zone list from MRFProperties
    creating MRF zone: MRF1Agitateur1
    creating MRF zone: MRF2Agitateur2
No finite volume options present

#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::divide(Foam::Field<double>&, Foam::UList<double> const&, Foam::UList<double> cons                                                        t&) at ??:?
#4  Foam::tmp<Foam::GeometricField<double, Foam::fvPatchField, Foam::volMesh> > Foam::opera                                                        tor/<Foam::fvPatchField, Foam::volMesh>(Foam::tmp<Foam::GeometricField<double, Foam::fvPatc                                                        hField, Foam::volMesh> > const&, Foam::GeometricField<double, Foam::fvPatchField, Foam::vol                                                        Mesh> const&) at ??:?
#5  Foam::RASModels::kEpsilon<Foam::IncompressibleTurbulenceModel<Foam::transportModel> >::                                                        correctNut() at ??:?
#6  ? at ??:?
#7  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#8  ? at ??:?
lethu is offline   Reply With Quote

Old   May 16, 2017, 12:03
Default
  #2
Member
 
power
Join Date: Jun 2014
Posts: 86
Rep Power: 8
enginpower is on a distinguished road
It seems to be some diverging problem or some division-by-zero problem.
I have the same with my case

Try to review all the boundary condition set, if they are all ok.

I can say just this.
Wait for more expert followers
enginpower is offline   Reply With Quote

Old   May 16, 2017, 14:17
Default
  #3
New Member
 
Rana
Join Date: May 2017
Posts: 10
Rep Power: 5
Rana shaharyar is on a distinguished road
im getting divergence problem what does it mean. during cfd simulation.
Rana shaharyar is offline   Reply With Quote


Old   May 17, 2017, 03:44
Default
  #5
New Member
 
Kim Le
Join Date: May 2016
Location: Lyon, France
Posts: 21
Rep Power: 6
lethu is on a distinguished road
Maybe this could help :

Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  4.x                                   |
|   \\  /    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
{
    agitateur1
    {
        type            noSlip;
    }
	
	agitateur2
    {
        type            noSlip;
    }
	
	bassin_wall
    {
        type            noSlip;
    }

    bassin_bottom
    {
        type            noSlip;
    }

    bassin_top
    {
        type            pressureInletOutletVelocity;
        value           uniform (0 0 0);
    }

}
Code:
    object      p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField   uniform 0;

boundaryField
{
	agitateur1
    {
        type            zeroGradient;
    }
	
	agitateur2
    {
        type            zeroGradient;
    }

    bassin_wall
    {
        type            zeroGradient;
    }
	
	bassin_bottom
    {
        type            zeroGradient;
    }

    bassin_top
    {
        type            totalPressure;
		p0              uniform 0;
    }
}
Code:
    object      k;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField   uniform 1;

boundaryField
{
	agitateur1
    {
        type            kqRWallFunction;
        value           uniform 0;
    }
	
	agitateur2
    {
        type            kqRWallFunction;
        value           uniform 0;
    }

    bassin_wall
    {
        type            kqRWallFunction;
        value           uniform 0;
    }
	
	bassin_bottom
    {
        type            kqRWallFunction;
        value           uniform 0;
    }

    bassin_top
    {
        type            calculated;
        value           uniform 0;
    }
}
Code:
    object      epsilon;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField   uniform 0;

boundaryField
{	
	agitateur1
    {
        type            epsilonWallFunction;
        value           uniform 0;
    }
	
	agitateur2
    {
        type            epsilonWallFunction;
        value           uniform 0;
    }

    bassin_wall
    {
        type            epsilonWallFunction;
        value           uniform 0;
    }
	
	bassin_bottom
    {
        type            epsilonWallFunction;
        value           uniform 0;
    }

    bassin_top
    {
        type            calculated;
        value           uniform 0;
    }
}
Code:
    object      nut;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField   uniform 0;

boundaryField
{
    agitateur1
    {
        type            nutkWallFunction;
        value           uniform 0;
    }
	
	agitateur2
    {
        type            nutkWallFunction;
        value           uniform 0;
    }

    bassin_wall
    {
        type            nutkWallFunction;
        value           uniform 0;
    }
	
	bassin_bottom
    {
        type            nutkWallFunction;
        value           uniform 0;
    }

    bassin_top
    {
        type            calculated;
		value           uniform 0;
    }
}
lethu is offline   Reply With Quote

Old   May 17, 2017, 08:02
Default
  #6
Member
 
power
Join Date: Jun 2014
Posts: 86
Rep Power: 8
enginpower is on a distinguished road
At a fast look it seems you put everywhere 0 values. So what do you need to calculate if everywhere each variable is 0?

What are your real boundary conditions? And what about the real case you want to simulate? Geometry?
enginpower is offline   Reply With Quote

Old   May 17, 2017, 08:24
Default
  #7
New Member
 
Kim Le
Join Date: May 2016
Location: Lyon, France
Posts: 21
Rep Power: 6
lethu is on a distinguished road
I am trying to simulate a tank full of water with some agitation at the bottom. In MFRProperties, I set two agitators.

Code:
MRF1Agitateur1
{
    cellZone    agitateur1;
    active      yes;

    // Fixed patches (by default they 'move' with the MRF zone)
    nonRotatingPatches ("(bassin_top|bassin_bottom)");

    origin    (1 0 0);
    axis      (0 0 1);
    omega     12.32;
}

MRF2Agitateur2
{
    cellZone    agitateur2;
    active      yes;

    // Fixed patches (by default they 'move' with the MRF zone)
    nonRotatingPatches ("(bassin_top|bassin_bottom)");

    origin    (-1 0 0);
    axis      (0 0 1);
    omega     12.32;
}

I guess that my "bassin_top" BC is not good (?)
Because in reality, it would be a free surface at an atmospheric pressure.

I wanted to do a steady-state and single phase simulation. I supposed that the tank is not filled during that time.

In fine, this is what I would like to observe but for now I just made a simple cylinder.

https://drive.google.com/open?id=0B7...kZaY0FzUjBIaE0

Last edited by lethu; May 17, 2017 at 08:25. Reason: Image not working
lethu is offline   Reply With Quote

Old   May 18, 2017, 07:03
Default
  #8
Member
 
power
Join Date: Jun 2014
Posts: 86
Rep Power: 8
enginpower is on a distinguished road
Sorry,
I can't help more.

Wait for a more expert foamer!
enginpower is offline   Reply With Quote

Reply

Thread Tools Search this Thread
Search this Thread:

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


Similar Threads
Thread Thread Starter Forum Replies Last Post
interFoam vs. simpleFoam channel flow comparison DanM OpenFOAM Running, Solving & CFD 12 January 31, 2020 15:26
MPI error with simpleFoam blaise OpenFOAM Running, Solving & CFD 0 November 7, 2015 14:01
simpleFoam parallel solver & Fluent polyhedral mesh Zlatko OpenFOAM Running, Solving & CFD 3 September 26, 2014 06:53
Laminar simpleFoam and inviscid simpleFoam herenger OpenFOAM Running, Solving & CFD 7 July 11, 2013 06:27
Trying to run a benchmark case with simpleFoam spsb OpenFOAM 3 February 24, 2012 09:07


All times are GMT -4. The time now is 05:52.