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::ddt(age)
+ fvm::div(phi, age)
- fvm::laplacian(turbulence->nuEff(),age)
==
dimensionedScalar("ageSource", age.dimensions()*dimensionSet(0,0,-1,0,0), 1)
);

ageEqn.solve();

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?

Thanks

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.

Regards,
Jose Santos

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

Guilherme

Set the term `alphaDt` to zero in these code blocks to turn off the turbulent diffusion contributions into the local mean age of air: https://www.cfd-online.com/Forums/op...tml#post795525