CFD Online URL
[Sponsors]
Home > Forums > OpenFOAM

problem in running icoFoam

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

Reply
 
LinkBack Thread Tools Display Modes
Old   March 16, 2012, 07:34
Default problem in running icoFoam
  #1
Member
 
valerio
Join Date: Jun 2009
Location: Nancy
Posts: 30
Rep Power: 7
Alucard is on a distinguished road
Good morning I'm dealing with a "simple" problem:
A cylindrical gemoetry where the fluid enters from two lateral openings and flows out from a small central tube.

The reynolds has to be small (i'm doing some parametric study in the range [0:20] with Reynolds defined as:
Re=UR/v v=cinematic viscosity, R= external radius of the geometry and U the inlet horizontal velocity (supposed constant).

Well the boundary conditions are simply an exit (imposed pression value ) at the bottom of the small pipe and an inlet radial velocity. On the top I put a symetry condition
icoFoam doesn"t converge in any case also for a Stokes case ( I removed the convective terms in the momentum equations)...so i think i'm doing something stupid with the mesh.

Can someone help me?
Thank you!


here I copy the blockMeshDict (i copied and modified this file from an ancient discussion: Help could anyone post a simple cylinder mesh)
HTML 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       dictionary;
    object      blockMeshDict;
}


// * * * * * * * * * * * * * * * * * * * * * * * * * *//
//meshGenApp blockMesh;
convertToMeters 0.017677664;
vertices
(
( 1.0 -1.0 -15.0) //0
( 1.0 1.0 -15.0)
( -1.0 1.0 -15.0)
( -1.0 -1.0 -15.0)

( 1.0 -1.0 -2.0) //4
( 1.0 1.0 -2.0)
( -1.0 1.0 -2.0)
( -1.0 -1.0 -2.0)

( 3.0 -3.0 0.0) //8
( 3.0 3.0 0.0)
( -3.0 3.0 0.0)
( -3.0 -3.0 0.0)

/*( 6.0 -6.0 0.0) //12
( 6.0 6.0 0.0)
( -6.0 6.0 0.0)
( -6.0 -6.0 0.0)
*/
( 40.0 -40.0 0.0) //12
( 40.0 40.0 0.0)
( -40.0 40.0 0.0)
( -40.0 -40.0 0.0)

( 3.0 -3.0 20.0) //16
( 3.0 3.0 20.0)
( -3.0 3.0 20.0)
( -3.0 -3.0 20.0)

( 40.0 -40.0 20.0) //20
( 40.0 40.0 20.0)
( -40.0 40.0 20.0)
( -40.0 -40.0 20.0)


/*( 6.0 -6.0 10.0) //20
( 6.0 6.0 10.0)
( -6.0 6.0 10.0)
( -6.0 -6.0 10.0)
*/
);
blocks
(
/*hex (0 1 2 3 4 5 6 7) (20 20 15) simpleGrading (1 1 1)         // tubo
hex (4 5 6 7 8 9 10 11) (20 20 10) simpleGrading (1 1 1)       //raccordo
hex (8 9 10 11 16 17 18 19) (20 20 50) simpleGrading (1 1 1)  //raccordo riportato fino in alto
hex (14 10 9 13 22 18 17 21) (30 20 50) simpleGrading (3 1 1) //
hex (15 11 10 14 23 19 18 22) (30 20 50) simpleGrading (3 1 1)//
hex (12 8 11 15 20 16 19 23) (30 20 50) simpleGrading (3 1 1) // 
hex (13 9 8 12 21 17 16 20) (30 20 50) simpleGrading (3 1 1)  //
*/

hex (0 1 2 3 4 5 6 7) (40 40 80) simpleGrading (1 1 1)         // tubo
hex (4 5 6 7 8 9 10 11) (40 40 80) simpleGrading (1 1 1)       //raccordo
hex (8 9 10 11 16 17 18 19) (40 40 100) simpleGrading (1 1 1)  //raccordo riportato fino in alto
hex (14 10 9 13 22 18 17 21) (60 40 100) simpleGrading (0.2 1 1) //
hex (15 11 10 14 23 19 18 22) (60 40 100) simpleGrading (0.2 1 1)//
hex (12 8 11 15 20 16 19 23) (60 40 100) simpleGrading (0.2 1 1) // 
hex (13 9 8 12 21 17 16 20) (60 40 100) simpleGrading (0.2 1 1)  //

);
edges
(
arc 2 1 (0.0 1.414214 -15.0)
arc 1 0 (1.414214 0.0 -15.0)
arc 0 3 (0.0 -1.414214 -15.0)
arc 3 2 (-1.414214 0.0 -15.0)
arc 6 5 (0.0 1.414214 -2.0)
arc 5 4 (1.414214 0.0 -2.0)
arc 4 7 (0.0 -1.414214 -2.0)
arc 7 6 (-1.414214 0.0 -2.0)
arc 10 9 (0.0 4.242641 0.0)
arc 9 8 (4.242641 0.0 0.0)
arc 8 11 (0.0 -4.242641 0.0)
arc 11 10 (-4.242641 0.0 0.0)
/*arc 14 13 (0.0 8.485281 0.0)
arc 13 12 (8.485281 0.0 0.0)
arc 12 15 (0.0 -8.485281 0.0)
arc 15 14 (-8.485281 0.0 0.0)
*/
arc 14 13 (0.0 56.56856 0.0)
arc 13 12 (56.56856 0.0 0.0)
arc 12 15 (0.0 -56.56856 0.0)
arc 15 14 (-56.56856 0.0 0.0)

arc 18 17 (0.0 4.242641 20.0)
arc 17 16 (4.242641 0.0 20.0)
arc 16 19 (0.0 -4.242641 20.0)
arc 19 18 (-4.242641 0.0 20.0)
/*arc 22 21 (0.0 8.485281 10.0)
arc 21 20 (8.485281 0.0 10.0)
arc 20 23 (0.0 -8.485281 10.0)
arc 23 22 (-8.485281 0.0 10.0)
*/

arc 22 21 (0.0 56.56856 20.0)
arc 21 20 (56.56856 0.0 20.0)
arc 20 23 (0.0 -56.56856 20.0)
arc 23 22 (-56.56856 0.0 20.0)
);

patches
(
patch outlet
(
(0 3 2 1)
)
patch symmetry
(
(16 17 18 19)
(17 21 22 18)
(18 22 23 19)
(20 16 19 23)
(20 21 17 16)
)

patch inlet1
(
(14 15 23 22)
)
patch inlet2
(
(13 12 20 21)
)

patch walls
(
(13 14 22 21) //lat walls
(15 12 20 23) //lat walls

(12 13 9 8) //bottom walls
(13 14 10 9) //bottom walls
(14 15 11 10) //bottom walls
(15 12 8 11) //bottom walls

(4 5 9 8)
(5 6 10 9)
(6 7 11 10)
(7 4 8 11)

(0 1 5 4)
(1 2 6 5)
(2 3 7 6)
(3 0 4 7)

)

);
mergePatchPairs
(
);
And for U and p:
HTML Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.1.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
{
    inlet1
    {
        type            fixedValue;
        value           uniform (1. 0 0);
    }
    inlet2
    {
        type            fixedValue;
        value           uniform (-1. 0 0);
    }
    walls
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }
    symmetry
    {
        type            zeroGradient;
//        type            fixedValue;
//        value           uniform (0 0 0);

    }

    outlet
    {
        type            zeroGradient;
    }

   
}

// ************************************************************************* //
HTML Code:
/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  2.1.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
{
    inlet1
    {
        type            zeroGradient;
    }
    inlet2
    {
        type            zeroGradient;
    }
    symmetry
    {
        type            zeroGradient;
    }

    walls
    {
        type            zeroGradient;
    }


    outlet
    {
        type            fixedValue;
        value           uniform 0;
    }

}

// ************************************************************************* //
Alucard is offline   Reply With Quote

Old   March 17, 2012, 10:50
Default
  #2
Senior Member
 
niaz's Avatar
 
A_R
Join Date: Jun 2009
Posts: 116
Rep Power: 7
niaz is on a distinguished road
plz copy your error firstly.
niaz is offline   Reply With Quote

Old   March 17, 2012, 11:01
Default
  #3
Member
 
valerio
Join Date: Jun 2009
Location: Nancy
Posts: 30
Rep Power: 7
Alucard is on a distinguished road
Quote:
Originally Posted by niaz View Post
plz copy your error firstly.
Hello,
Im still working on that and now everything it's fine (I did everything again from 0) , thanks anyway
Just to be "useful" for someone else I'll copy the mesh file I've done (and that now is working). The mesh size has to be "calibrated between the parts but everything it's fine now).

HTML 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       dictionary;
    object      blockMeshDict;
}

/*-----------------------------------------------------*\
| ========= | |
| \ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \ / O peration | Version: 1.3 |
| \ / A nd | Web: http://www.openfoam.org |
| \/ M anipulation | |
\*-----------------------------------------------------*/

FoamFile
{
version 2.0;
format ascii;

root "";
case "";
instance "";
local "";

class dictionary;
object blockMeshDict;
}

// * * * * * * * * * * * * * * * * * * * * * * * * * *//
//meshGenApp blockMesh;
convertToMeters 0.01;
vertices
(



( 5.0 0 0) //0     
( 0 5.0 0)
( -5.0 0 0)
( 0 -5.0 0)

( 5.0 0 50.) //4     
( 0 5.0 50.)
( -5.0 0 50.)
( 0 -5.0 50.)


( 10. 0 0) //8   
( 0 10. 0)
( -10. 0 0)
( 0 -10. 0)

( 10. 0 50.) //12     
( 0 10. 50.)
( -10. 0 50.)
( 0 -10. 50.)

( 5.0 0 -50) //16     
( 0 5.0 -50)
( -5.0 0 -50)
( 0 -5.0 -50)

( 10. 0 -50) //20     
( 0 10. -50)
( -10. 0 -50)
( 0 -10. -50)

( 100. 0 0) //24     
( 0 100. 0)
( -100. 0 0)
( 0 -100. 0)

( 100. 0 50) //28     
( 0 100. 50)
( -100. 0 50)
( 0 -100. 50)

);
blocks
(
// pipe
hex (0 1 2 3 4 5 6 7) (30 30 60) simpleGrading (1 1 2)
hex (0 8 9 1 4 12 13 5) (30 30 60) simpleGrading (1 1 2)   //1st quater
hex (1 9 10 2 5 13 14 6) (30 30 60) simpleGrading (1 1 2)  //2nd quater
hex (2 10 11 3 6 14 15 7) (30 30 60) simpleGrading (1 1 2) //3rd quater
hex (0 3 11 8 4 7 15 12) (30 30 60) simpleGrading (1 1 2)   //4th quater
//internal part of the domain
hex (16 17 18 19 0 1 2 3) (30 30 60) simpleGrading (1 1 1)
hex (16 20 21 17 0 8 9 1) (30 30 60) simpleGrading (1 1 1)   //1st quater
hex (17 21 22 18 1 9 10 2) (30 30 60) simpleGrading (1 1 1)  //2nd quater
hex (18 22 23 19 2 10 11 3) (30 30 60) simpleGrading (1 1 1) //3rd quater
hex (16 19 23 20 0 3 11 8) (30 30 60) simpleGrading (1 1 1)   //4th quater*/
// external part of the domain
hex (8 24 25 9 12 28 29 13) (60 30 60) simpleGrading (7 1 2)   //1st quater
hex (9 25 26 10 13 29 30 14) (60 30 60) simpleGrading (7 1 2)  //2nd quater
hex (10 26 27 11 14 30 31 15) (60 30 60) simpleGrading (7 1 2) //3rd quater
hex (11 27 24 8 15 31 28 12) (60 30 60) simpleGrading (7 1 2)   //4th quater


);
edges
(
 
arc 8 9    (7.071067812 7.071067812 0.0)   //1st quater
arc 12 13  (7.071067812 7.071067812 50.0) //1st quater
arc 20 21  (7.071067812 7.071067812 -50.0) //1st quater

arc 9 10   (-7.071067812 7.071067812 0.0)  //2nd quater
arc 13 14  (-7.071067812 7.071067812 50.0) //2nd quater
arc 21 22  (-7.071067812 7.071067812 -50.0) //2nd quater

arc 10 11  (-7.071067812 -7.071067812 0.0) //3rd quater
arc 14 15  (-7.071067812 -7.071067812 50.0) //3rd quater
arc 22 23  (-7.071067812 -7.071067812 -50.0) //3rd quater

arc 11 8   (7.071067812 -7.071067812 0.0)  //4th quater
arc 15 12  (7.071067812 -7.071067812 50.0) //4th quater
arc 23 20  (7.071067812 -7.071067812 -50.0) //4th quater

//******** EXTERNAL PART *********************************
arc 24 25  (70.71067812 70.71067812 0) //1st quater
arc 28 29  (70.71067812 70.71067812 50.0) //1st quater

arc 25 26  (-70.71067812 70.71067812 0.0)  //2nd quater
arc 29 30  (-70.71067812 70.71067812 50.0) //2nd quater

arc 26 27  (-70.71067812 -70.71067812 0.0) //3rd quater
arc 30 31  (-70.71067812 -70.71067812 50.0) //3rd quater

arc 27 24  (70.71067812 -70.71067812 0.0)  //4th quater
arc 31 28  (70.71067812 -70.71067812 50.0) //4th quater



);


boundary
(
    inlet1
    {
        type patch;
        faces
        (
            (24 25 29 28)
            
        );
    }

    inlet2
    {
        type patch;
        faces
        (
         (26 27 31 30) 
            
        );
    }

    outlet
    {
        type patch;
        faces
        (

      //  exit of the pipe
         (16 17 18 19)
         (16 20 21 17)
         (17 21 22 18)
         (18 22 23 19) 
         (19 23 20 16) 
      

        );
    }

    symmetry
    {
        type patch;
        faces
        (

               // top faces
         (12 28 29 13)
         (13 29 30 14)
         (14 30 31 15)
         (12 15 31 28) 

         (4 12 13 5) 
         (5 13 14 6)
         (6 14 15 7)
         (4 7 15 12)
         (4 5 6 7) 

        );
    }



    fixedWalls
    {
        type wall;
        faces
        ( 
         (25 26 30 29) //lat walls
         (27 24 28 31)

  
         (8 24 25 9)
         (9 25 26 10)
         (10 26 27 11)
         (11 27 24 8)

         // lateral draining pipe walls
         (9 8 20 21)
         (10 9 21 22)
         (11 10 22 23)
         (8 11 23 20)
        
        );

    }

);


mergePatchPairs
(
);


Alucard 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
IcoFoam Divergent Duct Simulation Problem Technoyoungman OpenFOAM Running, Solving & CFD 1 November 9, 2011 23:22
Problem running a parallel fluent job on local machine via mpd highhopes FLUENT 0 March 3, 2011 06:07
Problem running parallel. Hernán STAR-CCM+ 1 December 23, 2009 13:04
problem when running icoFoam on a complex shape flow field wendywu OpenFOAM 1 May 21, 2009 00:40
Problem with IcoFoam in parallel skabilan OpenFOAM Running, Solving & CFD 12 April 1, 2008 06:55


All times are GMT -4. The time now is 19:51.