[huangxianbei]

Hi,all:
I implement a model as here http://www.cfd-online.com/Forums/ope...model-les.html

I want to output the model coefficient as it runs, the most simple way may be something like

Info<< "Uncorrected Ubar = " << magUbarStar.value() << tab
            << "pressure gradient = " << gradP.value() << endl;

as in channelfoam
However, when I tried like this, it seems that it doesn't work.
Can anyone point out whether it's possible?
I tried another method, by defining it as a scalarfield to output as the other variables such as U and p,but this requires an initial value, when I set the internal filed of it , the value keeps the same all the time, how can it be?

Xianbei

[huangxianbei]

well, I find that if I define it as a scalarfield, it's almost the case of nuSgs which doesn't output as a field, just a file containing boundary conditions.

So is there any way to do this?

[huangxianbei]

I find a much more easier way!

Defining cV as a volScalarField(don't do the volume average), so cV can be output as a field. Then we can calculate the volumeAverage of cV by using swak4Foam.

However, the question is:
Why the cV field is can output as a field? Here is the model
DVM-Cvout.tar.gz

[huangxianbei]

Here is the initial condition for cV

dimensions      [ 0 0 0 0 0 0 0 ];

internalField   uniform 0.07;

boundaryField
{
    in
    {
        type            calculated;
        value           uniform 0.07;
    }
    out
    {
        type            calculated;
        value           uniform 0.07;
    }
   
    hub
    {
        type            calculated;
        value           uniform 0.07;
    }

    shroud
    {
        type            calculated;
        value           uniform 0.07;
    }
    
    blade
    {
        type            calculated;
        value           uniform 0.07;
    }

    peri1
    {
        type            cyclic;
    }
    peri2
    {
        type            cyclic;
    }
}

I'm considering whether it's right to implement such kind of initial condition

[huangxianbei]

The problem should be in the code that nuSgs is not output correctly, only the boundary values found. How do this come from? Can anyone help?