|
[Sponsors] |
March 31, 2014, 14:00 |
Setting of Spalart–Allmaras turbulence model
|
#1 |
Member
Join Date: Nov 2012
Posts: 62
Rep Power: 14 |
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. |
|
|
|
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 |