|
[Sponsors] |
July 31, 2012, 15:45 |
Modifing Turbulence Production Term
|
#1 |
New Member
Join Date: Feb 2012
Posts: 6
Rep Power: 14 |
I wish to add something along the lines of f1()*P() to the production term of the Spalart Allmaras turbulence model. Where P() is defined as follows with relations based on the Boussinesq approximation:
I thought I could accomplish this using R() as defined in SpallartAllmaras.C along with Code:
const volScalarField P = fvm::Sp(fvc::grad(U_),R()); but this has not worked. Am I using Sp correctly? Or does my error come from a tensor/scalar field mismatch? SpalartAllmaras::R() is shown below for your reference. Any insight is appreciated. Code:
tmp<volSymmTensorField> SpalartAllmaras::R() const { return tmp<volSymmTensorField> ( new volSymmTensorField ( IOobject ( "R", runTime_.timeName(), mesh_, IOobject::NO_READ, IOobject::NO_WRITE ), ((2.0/3.0)*I)*k() - nut()*twoSymm(fvc::grad(U_)) ) ); } |
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
momentum source term | zwdi | FLUENT | 14 | June 27, 2017 16:40 |
Turbulence postprocessing | Mohsin | FLUENT | 2 | October 3, 2016 15:18 |
ATTENTION! Reliability problems in CFX 5.7 | Joseph | CFX | 14 | April 20, 2010 16:45 |
bouyancy term in epsilon equation | Michael | Main CFD Forum | 1 | June 25, 1999 11:20 |
Turbulence in Turbomachinery | Erich F. | Main CFD Forum | 15 | June 23, 1999 00:37 |