[Thangam]

I have a domain inside which a plate is vibrating as a sinusoidal function. I have managed to impart the mesh motion. The issue is now obtaining the time averaged fields with fieldAverage. Because, what I really need is an ensemble average with respect to different phase shift/positions of the plate. If I turn on the fieldAverage from the controlDict, this is going to average over all the time steps (in-turn all the phases of impingement plate motion) which is going to result in an erroneous time averaged flow field? So, to avoid this, I decided to keep the fieldAverage option off and just write the instantaneous fields at different phase locations of the impingement plate motion. At this point I have about 10 time directories at each phase location of the plate motion(ie 0, pi/2, pi, 3pi/2, 2pi. ).

Now I decided to use execFlowFunctionObjects to get the time average of the fields during post-processing. I am looking to average all the 10 time directories at a particular phase location(for example at pi/2). I am not sure if execFlowFunctionObjects does what I need. Is there a cleaner/correct way of doing this in runtime?If so, the fieldAverage function has to be modified I guess. Also now I am confused how Umean and the fluctuating components and in-turn the average over time is calculated. Is calculating the time averaged fields based on this mean velocity locked at a particular phase shift correct? I'd be really pleased to get some pointers!



Thanks a lot.


Thangam.