CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > OpenFOAM Pre-Processing

Cyclic boundary conditions concerning channel with constant flow rate

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

Reply
 
LinkBack Thread Tools Display Modes
Old   November 16, 2012, 13:42
Default Cyclic boundary conditions concerning channel with constant flow rate
  #1
New Member
 
Join Date: Oct 2012
Posts: 16
Rep Power: 4
pascool is on a distinguished road
Hi there,

I'm trying to simulate a ribbed 2D channel in fully developed flow, by applying cyclic boundary conditions in streamwise direction and using simpleFoam. The case shall be run at a given bulk Reynolds number so that the flow rate must be specified.

My blockMeshDict file looks as follows:

Code:
boundary
(
     inlet
     {
          type cyclic;
          neighbourPatch outlet;
          faces
          (
            (0 21 22 1)
            (1 22 23 2)
            (2 23 24 3)
            (3 24 25 4)
          );
     }

     outlet
     {
         type cyclic;
         neighbourPatch inlet;
         faces
         (
            (16 17 38 37)
            (17 18 39 38)
            (18 19 40 39)
            (19 20 41 40)
          );
     }

     wall1
     {
          type wall;
          faces
          (
            (21 0 8 29)
          );
     }
     
     wall2
     {
          type wall;
          faces
          (
            (26 27 6 5)
            (27 28 7 6)
            (28 29 8 7)
          );
     }
     
     wall3
     {
          type wall;
          faces
          (
            (26 5 13 34)
          );
     }

     wall4
     {
          type wall;
          faces
      (
            (13 14 35 34)
            (14 15 36 35)
            (15 16 37 36)
      );
     }

     wall5
     {
          type wall;
          faces
          (
            (4 25 33 12)
            (12 33 41 20)
          );
     }

    frontAndBack
    {
          type empty;
          faces
          (
            (0 1 9 8)
            (1 2 10 9)
            (2 3 11 10)
            (3 4 12 11)
            (5 6 14 13)
            (6 7 15 14)
            (7 8 16 15)
            (8 9 17 16)
            (9 10 18 17)
            (10 11 19 18)
            (11 12 20 19)
            (21 29 30 22)
            (22 30 31 23)
            (23 31 32 24)
            (24 32 33 25)
            (26 34 35 27)
            (27 35 36 28)
            (28 36 37 29)
            (29 37 38 30)
            (30 38 39 31)
            (31 39 40 32)
            (32 40 41 33)
          );
    }
);
I don't know how to specify the initial and BC? My first attempt for p and U was:

Code:
FoamFile
{
    version     2.0;
    format ascii;
    class volScalarField;
    object p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField uniform 0;

boundaryField {
    inlet
    {
        type cyclic;
    }

    outlet
    {
        type cyclic;
    }
    wall1
    {   
            type            zeroGradient;
    }
    wall2
    {
        type            zeroGradient;
    }
    wall3
    {
        type            zeroGradient;
    }
    wall4
    {
        type            zeroGradient;
    }
    wall5
    {
        type            zeroGradient;
    }
    frontAndBack
    {        
   type            empty;
    }
}
Code:
FoamFile
{
    version     2.0;
    format      ascii;
    class       volVectorField;
    object      U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

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

internalField   uniform (2.98 0 0);

boundaryField
{
    inlet
    {
        type            cyclic;
    }
    outlet
    {
        type            cyclic;
    }
    wall1
    {
        type            fixedValue;
        value           uniform (0 0 0);

    }
    wall2
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }
    wall3
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }
    wall4
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }
    wall5
    {
       type            fixedValue;
        value           uniform (0 0 0);
    }
    frontAndBack
    {
        type            empty;
    }
Any help apprechiated!
Best regards,
Pascal
pascool is offline   Reply With Quote

Old   November 16, 2012, 14:04
Default
  #2
Member
 
Meindert de Groot
Join Date: Jun 2012
Location: Netherlands
Posts: 34
Rep Power: 5
meindert is on a distinguished road
Did you have a look at the channelFoam solver?

Regarding the initial conditions, I would suggest using the perturbU utility. It generates a perturbed laminar velocity field. You should be able to find more information about it on the forums.
The boundary conditions look fine, but be aware that imposing periodic boundary conditions will cause problems in the pressure field if you do not treat the pressure field in a special way (the channelFoam solver takes care of this).
meindert 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
Radiation interface hinca CFX 15 January 26, 2014 18:11
An error has occurred in cfx5solve: volo87 CFX 5 June 14, 2013 17:44
CFX13 Post Periodic interface EtaEta CFX 7 December 8, 2011 18:15
Face ordering problem for nonrectangular cyclic boundary conditions cfdmarkus OpenFOAM Other Meshers: ICEM, Star, Ansys, Pointwise, GridPro, Ansa, ... 3 August 17, 2011 15:07
RPM in Wind Turbine Pankaj CFX 9 November 23, 2009 05:05


All times are GMT -4. The time now is 23:13.