[cryple]

Dear all,

I thought I would be standing in front of just a little problem, when I wanted to add some run-tim-calculated fields to the rhoCentralFoam solver.

I would like to calculate e.g. mach number on-the-fly for some easier sampling and so added following to the createFields.H (testet on wedge15Ma5 tutorial):

volScalarField Ma
(
    IOobject
    (
        "Ma",
        runTime.timeName(),
        mesh,
        IOobject::NO_READ,
        IOobject::AUTO_WRITE
    ),
    mag(U)/sqrt(1.4/psi)
);

The result is a file "Ma" in the run time folders, which has a uniform dimensionless internal field of 5.000001 in it, whereas the boundary "obstacle" is comprising of a nonuniform list of 80 scalars.

I want that nonuniform list of scalars for my internal field. What am I doing wrong?

Thanks in advance!
Best
Daniel

[gschaider]

Originally Posted by cryple

Dear all,

I thought I would be standing in front of just a little problem, when I wanted to add some run-tim-calculated fields to the rhoCentralFoam solver.

I would like to calculate e.g. mach number on-the-fly for some easier sampling and so added following to the createFields.H (testet on wedge15Ma5 tutorial):

volScalarField Ma
(
    IOobject
    (
        "Ma",
        runTime.timeName(),
        mesh,
        IOobject::NO_READ,
        IOobject::AUTO_WRITE
    ),
    mag(U)/sqrt(1.4/psi)
);

The result is a file "Ma" in the run time folders, which has a uniform dimensionless internal field of 5.000001 in it, whereas the boundary "obstacle" is comprising of a nonuniform list of 80 scalars.

I want that nonuniform list of scalars for my internal field. What am I doing wrong?

Thanks in advance!
Best
Daniel

You've calculated the field for the initial conditions. To have the current state written out you've got to update it at the end of the timestep (something like "Ma=mag(U)/sqrt(1.4/psi);"

More elegant (without the modification of the solver) would be a functionObject

[cryple]

Thanks a lot for your fast reply, gschaider.

I realised that actualy what I did, was just create an IOobject at the beginning of the solver run, without taking care of performing any calculation for Ma !inside! the time loop at all...

I also thought about using functionObjects, which would be much better, I think but I can't find any documentation on the possibilities and options I have there.

I will post the solution when I got more into it.

Thanks again!
Daniel