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

Setting of Spalart–Allmaras turbulence model

Register Blogs Community New Posts Updated Threads Search

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   March 31, 2014, 14:00
Default Setting of Spalart–Allmaras turbulence model
  #1
Member
 
Join Date: Nov 2012
Posts: 62
Rep Power: 14
Naruto is on a distinguished road
Hello.

Currently I am simulating flow over a cylinder. For that purpose I am using pimpleFoam solver. I am planning to use Spalart–Allmaras turbulence model.
I am using Dirichlet boundary condition. I am having a little difficulty in setting up the turbulence model boundary condition. I have used the formula
nutilda = root(1.5) * U * I * l

Here,
U = inflow velocity = 5.07, I = Turbulence intensity = 0.1% and l = turbulent length scale = 0.02
Now my settings of NuTilda looks like this:

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

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

internalField uniform 0.00124;

boundaryField
{
Rotor-symmetry
{
type symmetryPlane;
}
Rotor-AMI
{
type cyclicAMI;
value uniform 0;
}
cylinder
{
type nutUSpaldingWallFunction;
value uniform 0;
}
Stator-symmetry
{
type symmetryPlane;
}
Top
{
type nutUSpaldingWallFunction;
value uniform 0;
}
Stator-AMI
{
type cyclicAMI;
value uniform 0;
}
inlet
{
type fixedValue;
value uniform 0.00124;
}
sides
{
type nutUSpaldingWallFunction;
value uniform 0;
}
outlet
{
type inletOutlet;
inletValue uniform 0.00124;
value uniform 0.00124;
}
}


// ************************************************** *********************** //
and Nu looks like this:
/*--------------------------------*- C++ -*----------------------------------*\
| ========= | |
| \\ / F ield | OpenFOAM: The Open Source CFD Toolbox |
| \\ / O peration | Version: 2.2.2 |
| \\ / A nd | Web: www.OpenFOAM.org |
| \\/ 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
{
Rotor-symmetry
{
type symmetryPlane;
}
Rotor-AMI
{
type cyclicAMI;
value uniform 0;
}
cylinder
{
type nutUSpaldingWallFunction;
value uniform 0;
}
Stator-symmetry
{
type symmetryPlane;
}
Top
{
type nutUSpaldingWallFunction;
value uniform 0;
}
Stator-AMI
{
type cyclicAMI;
value uniform 0;
}
inlet
{
type calculated;
value uniform 0;
}
sides
{
type nutUSpaldingWallFunction;
value uniform 0;
}
outlet
{
type calculated;
value uniform 0;
}
}


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

Are my methods and boundary conditions ok?

Thank you.
__________________
Happy Foaming

Last edited by Naruto; April 1, 2014 at 02:41.
Naruto is offline   Reply With Quote

Reply


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
SimpleFoam case with SpalartAllmaras turbulence model implemented nedved OpenFOAM Running, Solving & CFD 2 November 30, 2014 23:43
[snappyHexMesh] determining displacement for added points CFDnewbie147 OpenFOAM Meshing & Mesh Conversion 1 October 22, 2013 10:53
Wrong calculation of nut in the kOmegaSST turbulence model FelixL OpenFOAM Bugs 27 March 27, 2012 10:02
Low Reynolds k-epsilon model YJZ ANSYS 1 August 20, 2010 14:57
Cells with t below lower limit Purushothama Siemens 2 May 31, 2010 22:58


All times are GMT -4. The time now is 16:04.