CFD Online Logo CFD Online URL
Home > Forums > OpenFOAM Programming & Development

Modifing Turbulence Production Term

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

LinkBack Thread Tools Display Modes
Old   July 31, 2012, 14:45
Post Modifing Turbulence Production Term
New Member
Join Date: Feb 2012
Posts: 6
Rep Power: 7
twray is on a distinguished road
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

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.

 tmp<volSymmTensorField> SpalartAllmaras::R() const
    return tmp<volSymmTensorField>
        new volSymmTensorField
            ((2.0/3.0)*I)*k() - nut()*twoSymm(fvc::grad(U_))

twray is offline   Reply With Quote


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
momentum source term zwdi FLUENT 14 June 27, 2017 15:40
Turbulence postprocessing Mohsin FLUENT 2 October 3, 2016 14:18
ATTENTION! Reliability problems in CFX 5.7 Joseph CFX 14 April 20, 2010 15:45
bouyancy term in epsilon equation Michael Main CFD Forum 1 June 25, 1999 10:20
Turbulence in Turbomachinery Erich F. Main CFD Forum 15 June 22, 1999 23:37

All times are GMT -4. The time now is 10:30.