Problems with simpleFoam

samiam1000 · October 4, 2011, 08:52

Dear all,

when I try to run simpleFoam (to simulate wind over terrain - let's say that I am using the turbine siting tutorial, changing the stl fil!) I get this error:

Code:

sammy@nash:~/Desktop/mesh_test/cases/suisse$ simpleFoam 
/*---------------------------------------------------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.0.1                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
Build  : 2.0.1-cce6c678443c
Exec   : simpleFoam
Date   : Oct 04 2011
Time   : 14:48:43
Host   : nash
PID    : 20118
Case   : /home/sammy/Desktop/mesh_test/cases/suisse
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster
allowSystemOperations : Disallowing user-supplied system call operations

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

Create mesh for time = 0

Reading field p

Reading field U

#0  Foam::error::printStack(Foam::Ostream&) in "/opt/openfoam201/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#1  Foam::sigFpe::sigHandler(int) in "/opt/openfoam201/platforms/linux64GccDPOpt/lib/libOpenFOAM.so"
#2   in "/lib/x86_64-linux-gnu/libc.so.6"
#3   in "/lib/x86_64-linux-gnu/libm.so.6"
#4  log in "/lib/x86_64-linux-gnu/libm.so.6"
#5  Foam::incompressible::atmBoundaryLayerInletVelocityFvPatchVectorField::updateCoeffs() in "/opt/openfoam201/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
#6  Foam::incompressible::atmBoundaryLayerInletVelocityFvPatchVectorField::atmBoundaryLayerInletVelocityFvPatchVectorField(Foam::fvPatch const&, Foam::DimensionedField<Foam::Vector<double>, Foam::volMesh> const&, Foam::dictionary const&) in "/opt/openfoam201/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
#7  Foam::fvPatchField<Foam::Vector<double> >::adddictionaryConstructorToTable<Foam::incompressible::atmBoundaryLayerInletVelocityFvPatchVectorField>::New(Foam::fvPatch const&, Foam::DimensionedField<Foam::Vector<double>, Foam::volMesh> const&, Foam::dictionary const&) in "/opt/openfoam201/platforms/linux64GccDPOpt/lib/libincompressibleRASModels.so"
#8  Foam::fvPatchField<Foam::Vector<double> >::New(Foam::fvPatch const&, Foam::DimensionedField<Foam::Vector<double>, Foam::volMesh> const&, Foam::dictionary const&) in "/opt/openfoam201/platforms/linux64GccDPOpt/bin/simpleFoam"
#9  Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>::GeometricBoundaryField::GeometricBoundaryField(Foam::fvBoundaryMesh const&, Foam::DimensionedField<Foam::Vector<double>, Foam::volMesh> const&, Foam::dictionary const&) in "/opt/openfoam201/platforms/linux64GccDPOpt/bin/simpleFoam"
#10  Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>::readField(Foam::dictionary const&) in "/opt/openfoam201/platforms/linux64GccDPOpt/bin/simpleFoam"
#11  Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>::readField(Foam::Istream&) in "/opt/openfoam201/platforms/linux64GccDPOpt/bin/simpleFoam"
#12  Foam::GeometricField<Foam::Vector<double>, Foam::fvPatchField, Foam::volMesh>::GeometricField(Foam::IOobject const&, Foam::fvMesh const&) in "/opt/openfoam201/platforms/linux64GccDPOpt/bin/simpleFoam"
#13  
 in "/opt/openfoam201/platforms/linux64GccDPOpt/bin/simpleFoam"
#14  __libc_start_main in "/lib/x86_64-linux-gnu/libc.so.6"
#15  
 in "/opt/openfoam201/platforms/linux64GccDPOpt/bin/simpleFoam"
Floating point exception

Could anyone help?

Thanks a lot,

Samuele.

PS: the mesh is ok (the command chackMesh says this!)..

philipp. · October 4, 2011, 11:21

Hello Samuele,

there seems to be a problem with the "Reading/calculating face flux field phi". Maybe you should check your initial and boundary conditions again. For more help you probably need to post more details.

Regards, Philipp

samiam1000 · October 4, 2011, 11:22

Hi Philipp,

I can share my case, if you want. Maybe it is better to copy in a message just few files. Please, tell me which files I should post in order to allow you to help!

Thanks,

Samuele

philipp. · October 4, 2011, 13:32

Hello Samuele,

I just checked the tutorials - you are using OF2.0.1, right? There is a turbineSiting tutorial for the windSimpleFoam solver. Maybe you should use that solver...
If that does not work, which files you edited and which are like the ones of the tutorial.

Philipp

samiam1000 · October 4, 2011, 13:57

Dear Philipp,

it is curios that I get the same errore using the solver windSimpleFoam.

Also, if I use the file terrain.stl (the one that you can find in the turbineSitting tutorial) the command simpleFoam in my case works.

So, let me copy here all my files.. That's very strange, imho..

Here are the files:

1. epsilon

Code:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  1.7.1                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    location    "0";
    object      epsilon;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

#include        "include/initialConditions"

internalField   uniform $turbulentEpsilon;

boundaryField
{
    #include "include/ABLConditions"

    "terrain_.*"
    {
        type            epsilonWallFunction;
        Cmu             0.09;
        kappa           0.4;
        E               9.8;
        value           $internalField;
    }

    outlet
    {
        type            zeroGradient;
    }

    inlet
    {
        type            atmBoundaryLayerInletEpsilon;
        Ustar           $Ustar;
        z               $zDirection;
        z0              $z0;
        value           $internalField;
        zGround         $zGround;
    }

    ground
    {
        type            zeroGradient;
    }

    #include "include/sideAndTopPatches"
}


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

2. p

Code:

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

#include        "include/initialConditions"

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

internalField   uniform $pressure;

boundaryField
{
    inlet
    {
        type            zeroGradient;
    }

    outlet
    {
        type            fixedValue;
        value           $internalField;
    }

    "terrain_.*"
    {
        type            zeroGradient;
    }

    ground
    {
        type            zeroGradient;
    }

    #include "include/sideAndTopPatches"
}

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

3. nut

Code:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  1.7.1                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    location    "0";
    object      nut;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField   uniform 0;

boundaryField
{
    inlet
    {
        type            calculated;
        value           uniform 0;
    }

    outlet
    {
        type            calculated;
        value           uniform 0;
    }

    "terrain_.*"
    {
        type            nutkRoughWallFunction;
        Ks              0.2; //Ks = 20 Z0
        Cs              0.5;
        value           uniform 0.0;
    }

    ground
    {
        type            calculated;
        value           uniform 0;
    }

    #include "include/sideAndTopPatches"
}


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

4. k

Code:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  1.7.1                                 |
|   \\  /    A nd           | Web:      http://www.OpenFOAM.com               |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    object      k;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

#include        "include/initialConditions"

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

internalField   uniform $turbulentKE;

boundaryField
{
    #include "include/ABLConditions"

    outlet
    {
        type            inletOutlet;
        inletValue      uniform 0.0;
        value           $internalField;
    }

    inlet
    {
        type            fixedValue;
        value           $internalField;
    }

    "terrain_.*"
    {
        type            kqRWallFunction;
        value           uniform 0.0;
    }

    ground
    {
        type            zeroGradient;
    }

    #include "include/sideAndTopPatches"
}


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

5. U

Code:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  1.7.1                                 |
|   \\  /    A nd           | Web:      http://www.OpenFOAM.com               |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volVectorField;
    location    "0";
    object      U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

#include        "include/initialConditions"

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

internalField   uniform $flowVelocity;

boundaryField
{
    #include "include/ABLConditions"

    outlet
    {
        type            inletOutlet;
        inletValue      uniform (0 0 0);
        value           $internalField;
    }

    inlet
    {
        type            atmBoundaryLayerInletVelocity;
        Uref            $Uref;
        Href            $Href;
        n               $windDirection;
        z               $zDirection;
        z0              $z0;
        value           $internalField;
        zGround         $zGround;
    }

    "terrain_.*"
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }

    ground
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }

    #include "include/sideAndTopPatches"
}


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

6. ABLConditions

Code:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  1.7.1                                 |
|   \\  /    A nd           | Web:      http://www.OpenFOAM.com               |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/

Ustar                0.82;
Uref                 10.0;
Href                 20;
z0                   0.1;
turbulentKE          1.3;
windDirection        (1 0 0);
zDirection           (0 0 1);
zGround              935.0;

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

7. initialConditions

Code:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  1.7.1                                 |
|   \\  /    A nd           | Web:      http://www.OpenFOAM.com               |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/

flowVelocity         (0 0 0);
pressure             0;
turbulentKE          1.3;
turbulentEpsilon     0.01;

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

8. sideAndTopPatches

Code:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  1.7.1                                 |
|   \\  /    A nd           | Web:      http://www.OpenFOAM.com               |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/

top
{
    type slip;
}

sides
{
    type slip;
}

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

philipp. · October 5, 2011, 05:30

Hi Samuele,
if it works with the stl-file of the tutorial there has to be a problem with your file respectly with your mesh. Even if checkMesh says that the mesh is ok there still can be problems. Maybe you can use a simpler geometriem at first.
Otherwise you can try different initial and boundary conditions. Maybe that will help.
Philipp

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
Help! problems in calculating forces with SimpleFoam	DLC	OpenFOAM	3	April 25, 2022 23:43
[ICEM] Problems with coedge curves and surfaces	tommymoose	ANSYS Meshing & Geometry	6	December 1, 2020 11:12
Convergence Problems SimpleFOAM	Kutti	OpenFOAM	16	June 14, 2010 08:12
Convergence problems with simpleFoam on human airway	CedricVH	OpenFOAM Running, Solving & CFD	8	June 3, 2010 09:05
convergence problems using simpleFoam	sebbi	OpenFOAM Running, Solving & CFD	3	November 25, 2009 15:27

October 4, 2011, 11:21		#2
philipp. New Member Philipp Hofemeier Join Date: Sep 2011 Location: Freiberg, Germany Posts: 7 Rep Power: 14	Hello Samuele, there seems to be a problem with the "Reading/calculating face flux field phi". Maybe you should check your initial and boundary conditions again. For more help you probably need to post more details. Regards, Philipp

October 4, 2011, 11:22		#3
samiam1000 Senior Member Samuele Z Join Date: Oct 2009 Location: Mozzate - Co - Italy Posts: 520 Rep Power: 18	Hi Philipp, I can share my case, if you want. Maybe it is better to copy in a message just few files. Please, tell me which files I should post in order to allow you to help! Thanks, Samuele

October 4, 2011, 13:32		#4
philipp. New Member Philipp Hofemeier Join Date: Sep 2011 Location: Freiberg, Germany Posts: 7 Rep Power: 14	Hello Samuele, I just checked the tutorials - you are using OF2.0.1, right? There is a turbineSiting tutorial for the windSimpleFoam solver. Maybe you should use that solver... If that does not work, which files you edited and which are like the ones of the tutorial. Philipp

October 5, 2011, 05:30		#6
philipp. New Member Philipp Hofemeier Join Date: Sep 2011 Location: Freiberg, Germany Posts: 7 Rep Power: 14	Hi Samuele, if it works with the stl-file of the tutorial there has to be a problem with your file respectly with your mesh. Even if checkMesh says that the mesh is ok there still can be problems. Maybe you can use a simpler geometriem at first. Otherwise you can try different initial and boundary conditions. Maybe that will help. Philipp