[villet]

Hi Gijs,

check that "createMyFields.H" is included after "createFields.H". Also check the equation for dissipation rate. It doesn't look like the right one...

Anyway, I have another solution. You can call "sgsModel->epsilon()" that computes TKE dissipation estimation (ce*ksgs*sqrt(ksgs)/delta for eddy viscosity models).

Hope this helps!
Ville

[gwierink]

Hi Ville,

How are things? Thanks for your reaction. I did put "#include createMyFields.H" after createFields.H. The equation is just a guess, because I also tried:



Info<< "Calculating the dissipation.\n" << endl;
volScalarField dissipation
(
IOobject
(
"dissipation",
runTime.timeName(),
mesh,
IOobject::MUST_READ,
IOobject::AUTO_WRITE
),
sgsModel->epsilon()
);



and that didn't work. Well, it does write out a file called "dissipation" every timestep, but it only contains the same info as the initial conditions (i.e. zerGradient for most things and some fixed values for the two inlets). Any idea how to fix that?

Thanks, Gijs

[villet]

Uups, you need to update epsilon somewhere in the code:

dissipation = sgsModel->epsilon();

[gwierink]

So you mean I should put "dissipation = sgsModel->epsilon();" in myoodles.C and have it also in createMyFields.H ?

[gwierink]

Hi Ville,

Thanks for your hint. To recap for others, this is what worked (at least for a few time steps, still need to test it for longer times):

Put in myoodles.C the lines "#include createMyFields.H" and "dissipation = sgsModel->epsilon()". Then, create a file called "createMyFields.H" (or whatever you like of course) with the code:

Info<< "Calculating the dissipation.\n" << endl;
volScalarField dissipation
(
IOobject
(
"dissipation",
runTime.timeName(),
mesh,
IOobject::MUST_READ,
IOobject::AUTO_WRITE
),
mesh
);

Now every writing step a file called "dissipation" is written that contains a field of subgrid epsilon.

Many thanks for your help Ville ( and also many thanks to Ville Vuorinen for helping me out)!

Cheers,

Gijs

[villet]

That's good to hear you got it working!