[vivek05]

Hi,
I want to access "turbulent kinetic energy" and "nut" value in each cell. I am using OF-5.0 version. So I have implemented in following manner.

PHP Code:

volScalarField turb_nut = turbulence->nut();
     volScalarField turb_k = turbulence->k();
     forAll(mesh.cells(),celli) {
     if (condition) {
     Source[celli] = (...*(turb_k[celli]/(a+turb_nut[celli]))*...).value();
     } 
     else {
     Source[celli] = scalar(0.0);
      }
     } 


where Source is volScalarField
a is dimensionedScalar
... represents combination of constant and volScalarField
The way I am accessing "k" and "nut" is it correct??
What does -> operator do in turbulence->k()
Thanks,
Vivek S

[Santiago]

Originally Posted by vivek05

Hi,
I want to access "turbulent kinetic energy" and "nut" value in each cell. I am using OF-5.0 version. So I have implemented in following manner.

PHP Code:

volScalarField turb_nut = turbulence->nut();     volScalarField turb_k = turbulence->k();     forAll(mesh.cells(),celli) {     if (condition) {     Source[celli] = (...*(turb_k[celli]/(a+turb_nut[celli]))*...).value();     }      else {     Source[celli] = scalar(0.0);      }     } 


where Source is volScalarField
a is dimensionedScalar
... represents combination of constant and volScalarField
The way I am accessing "k" and "nut" is it correct??
What does -> operator do in turbulence->k()
Thanks,
Vivek S

turbulence->k() is the RESIDUAL turbulent kinetic energy. the total tke would be the sum of the resolved tke and said residual tke.

[vivek05]

Originally Posted by Santiago

turbulence->k() is the RESIDUAL turbulent kinetic energy. the total tke would be the sum of the resolved tke and said residual tke.

Thanks Santiago for your reply.

In LES simulation k() gives us subgrid kinetic energy. But I want to know how can I access ksgs in each cell? First I am storing turbulence-k() in variable called turb_k which is volumetric scalar field.

volScalarField turb_k = turbulence->k();
volSclarField turb_nut = turbulence->nut();

Then I am using forAll loop for spatial marching over cells

PHP Code:

forAll(sourceElsa , celli) {
If(condition)
sourceElsa = (...*turb_k[celli]/turb_nut[celli]*...);
} 


Is it correct way to access turbulent kinetic energy?
It has also taking lot of time to calculate sourceElsa. How can I program it efficiently?


Thanks
Vivek S