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

Realisable K-Epsilon Model in a Channel

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

Reply
 
LinkBack Thread Tools Display Modes
Old   June 28, 2013, 10:49
Default Realisable K-Epsilon Model in a Channel
  #1
New Member
 
Matt Mosquera
Join Date: Jun 2013
Location: Lewisburg, PA
Posts: 18
Rep Power: 4
mosquera is on a distinguished road
I've been struggling with modelling a flow through a channel. I have it set up as a two dimensional flow.

My objective is to model the flow using a realisable k-epsilon model. I've mapped experimental values at the inlet. I'm using simpleFoam for this. Unfortunately, I can't get the model to converge.

Within my model, SymmWall is the boundary at the top "wall". The boundary layer has been tripped in the flow and so I am using a symmetric boundary condition (I am not sure, but I have been assured by my advisor this is the correct method.

Inlet and Outlet refer to the inlet and outlet.

Wall refers to the bottom wall across which the flow is moving.

Anyway, on to my initial boundary conditions:

U
Code:
dimensions      [0 1 -1 0 0 0 0];

internalField   uniform (1 0 0); 

boundaryField
{
    SymmWall
    {
        type            symmetryPlane;
    }
    Wall
    {
        type            fixedValue;
        value           uniform (0 0 0);
    }
    Outlet
    {
        type            zeroGradient;
    }
    Inlet
    {
        type            timeVaryingMappedFixedValue;
       setAverage       0;
        perturb         1e-05;
        offset          (0 0 0);
        value           uniform (0 0 0);
  
    }
    frontAndBackPlanes
    {
        type            empty;
    }
}
k
Code:
dimensions      [0 2 -2 0 0 0 0];

internalField   uniform 1.5e-09;

boundaryField
{
    Wall
    {
        type            kLowReWallFunction;
        Cmu             0.09;
        kappa           0.41;
        E               9.8;
        value           uniform 8.1323e-04;
    }
    Outlet
    {
        type            zeroGradient;
    }
    Inlet
    {
        type            timeVaryingMappedFixedValue;
        setAverage      0;
        perturb         0;
        offset          0;
        value         uniform 8.1323e-04;

    }
     frontAndBackPlanes
     {
        type            empty;
     }
     SymmWall
     {
        type            symmetryPlane;
     }
    
}
epsilon
Code:
dimensions      [0 2 -3 0 0 0 0];

internalField   uniform 0.0008123;

boundaryField
{
    Wall
    {
        type            epsilonWallFunction;
        Cmu             0.09;
        kappa           0.41;
        E               9.8;
        value           uniform 0.0008123;
    }
    Outlet
    {
        type            zeroGradient;
    }
    Inlet
    {
        type            timeVaryingMappedFixedValue;
        setAverage      0;
        perturb         0;
        offset          0;
        value           uniform 0.0008132;
    }
    SymmWall
    {
        type            symmetryPlane;
    }
    frontAndBackPlanes
    {
        type            empty;
    }
}
Also, here is my current fvschemes file. The goal is to have the solution be of second order.

Code:
ddtSchemes
{
    default         steadyState;
}

gradSchemes
{
    default         Gauss linear;
    k            Gauss linear;
    epsilon        Gauss linear;
}

divSchemes
{
    default                          none;
    div(phi,U)                       Gauss linearUpwindV grad(U);
    div(phi,k)                       Gauss linearUpwind grad(k);
    div(phi,epsilon)                   Gauss linearUpwind grad(epsilon);
    div(phi,T)                      Gauss linear grad(T);
    div((nuEff*dev(T(grad(U)))))     Gauss linear;
}

laplacianSchemes
{
    default                    none;
    laplacian(nuEff,U)            Gauss linear corrected;
    laplacian((1|A(U)),p)        Gauss linear corrected;
    laplacian(DkEff,k)                Gauss linear corrected;
    laplacian(DepsilonEff,epsilon) Gauss linear corrected;
}

interpolationSchemes
{
    default         linear;
    interpolate(U)  linear;
    interpolate(T)  linear;
    interpolate(epsilon) linear;
//    interpolate(omega) limitedLinear 0 1000000;
}

snGradSchemes
{
    default         corrected;
}

fluxRequired
{
    default         no;
    p;
}
checkMesh yields that my mesh is OK. The mesh was imported from Gambit.

I could additionally attach the mapped values if anyone feels that neccesary. Are there any values I should be on the lookout for there? I'm already positive there's no zeros anywhere in those files.

Again, any assistance would be appreciated. I've been working at this for a couple of days now.

Thanks,
Matt
mosquera 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
Jump in epsilon values near the wall : low re k-epsilon model malaboss OpenFOAM Verification & Validation 1 February 1, 2013 17:36
Calculation of Epsilon in LES model Yuri Almeida OpenFOAM 2 January 18, 2013 08:43
Behavior of k and epsilon in a channel flow mn1363 Main CFD Forum 1 January 11, 2013 12:49
K - epsilon VS SST turbulence model Maicol Main CFD Forum 0 November 30, 2012 17:25
TCL model for channel flow Hui Main CFD Forum 2 July 25, 2008 23:28


All times are GMT -4. The time now is 11:09.