CFD Online Discussion Forums

CFD Online Discussion Forums (
-   OpenFOAM (
-   -   Dissabling turbulent dissipation (

lindner July 14, 2010 04:10

Dissabling turbulent dissipation
Hi all,
I have a solver that calculates the local mean age of air, it's the simpleFoam solver with a scalarTransportFoam with source term = 1.
The code of the scalar transport is this

fvScalarMatrix ageEqn
+ fvm::div(phi, age)
- fvm::laplacian(turbulence->nuEff(),age)
dimensionedScalar("ageSource", age.dimensions()*dimensionSet(0,0,-1,0,0), 1)


Now i want to test the local mean age solver but turning off the turbulent dissipation (turbulence->nuEff = 1).
I tried just removing the term from the equation but it didnt work. Is there a way to simple make the turbulent dissipation = 1 or add a simple constant?


santos July 14, 2010 10:16

Hi there,

What if you replace it with:
- fvm::laplacian(Dage,age)
where Dage is the diffusivity of age (dont know whether it make sense!)

Dage will have to be declared in createFields.H, and then defined in constant/transportProperties.

Jose Santos

lindner July 15, 2010 04:48

it worked!
Thanks for the help Jose, the code was implemented successfully.


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