[Mehrez]

Hi Foamers,

I am trying to build a groovyBC in which I need the average of the velocity (and pressure) field at the last output time.
Do you have an idea how can I read or access the velocity field at the previous output time ?

I know that I can use "oldTime(p)" for the last time step (=/= last output time).

Thank you for your help.

Mhrz

[Mehrez]

I found that using "swak4Foam", I can read a field of the previous time step: oldTime(fieldName) but nothing to read the previous output time.
Mhrz

[Ali Blues]

Hi,
Let's assume we have solved a transient scalar transport equation. Now for each time-step we want to evaluated the following global quantity (sort of an average of the entire field variable, so we end up with a scalar):



Now the issue is how i can read in the initial field , and in a more general context any of the previous time-steps?


Also a side-note that assuming I just have the following to compute:


I was thinking of putting the code like:

Index
{
        type swakExpression;
        valueType internalField;
        variables ( "Tmean = 0.05;" );
         accumulations ( sum );
        expression "mag(T-Tmean)*vol()";
        verbose true;
        outputControlMode timeStep;
       outputInterval  1;
}

but it gives me an error which I think concerns the subtraction of field from scalar:

Parser Error for driver FieldValueExpressionDriver at "1.7-11" :"field Tmean not existing or of wrong type"

So where is the silly mistake that I've made here!

Best
Ali

[wyldckat]

Quick answer: Tested with OpenFOAM 3.0.x, tutorial case "heatTransfer/buoyantSimpleFoam/circuitBoardCooling":

libs ("libthermalBaffleModels.so" "libOpenFOAM.so" "libsimpleFunctionObjects.so" "libsimpleSwakFunctionObjects.so");

functions
{
  Index
  {
          type swakExpression;
          valueType internalField;
          variables ("Tmean=0.05;");
          accumulations ( sum );
          expression "mag(T-Tmean)*vol()";
          verbose true;
          outputControlMode timeStep;
          outputInterval  1;
  }
}

The problem is that you had spaces in the variables definition, namely:

• Good: "Tmean=0.05;"
• Bad: "Tmean = 0.05;"

[Ali Blues]

Thanks a lot Bruno!
Completely overlooked that

I was also wondering if you would know how I can read-in the initial field , and well perhaps more generally previous time-step solutions, as I've indicated in the early part of my post.

Best
Ali

[Ali Blues]

So can anyone suggest on how I can read-in previous time-steps solution into swak4foam?
Thanks
Ali

[wyldckat]

Quick answer: swak4Foam has a function known as "oldTime". Some details and examples are available here: http://www.cfd-online.com/Forums/ope...ormations.html - there you'll also find a few tricks on how to store fields from previous times.

But as demonstrated in that thread, you cannot easily store/retrieve the field for any specific time step. You will need to do some clever gradual storage of values, so that you can continuously update only for the current time step, while accounting only for the previous accumulated values.

[Mehrez]

Thank you for participating in this tread.
The swak4Foam "storedVariables" is a good way to use the old fields averages but there is a bug when using it in parallel.
Regards,
Mhrz

[Mehrez]

Any idea ?
Thanks,
Mhrz