<div dir="ltr"><div><div>Hi Phillip, <br><br></div>Thanks for the reply, that was very clear. However, I'm now facing a problem while trying to load the avvgblah.f000* files from my case. Currently I'm running visnek on these averaged files, generating a .nek5000 file and trying to read that using visit or paraview. This procedure results in an error. The reader complains that the first step of a nek file must contain a mesh which the files generated by avg_all do not seem to have. What is the procedure for reading this data then?<br><br><br></div>Sincerely,<br></div><br><div class="gmail_quote"><div dir="ltr">On Sun, Feb 4, 2018 at 5:39 AM <<a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Dear Amitvikram,<br>
<br>
As far as I understand, the basic version of avg_all perform only a time<br>
average, i.e. the statistics are computed in 3D in every point of the<br>
domain. In case your flow has additional directions of homogeneity, then<br>
one can do an average over those directions as well; either by taking<br>
the output of avg_all and doing the average there, or direction during<br>
runtime (for the latter there is an implementation described here:<br>
<a href="http://www.ipd.anl.gov/anlpubs/2017/03/133610.pdf" rel="noreferrer" target="_blank">http://www.ipd.anl.gov/anlpubs/2017/03/133610.pdf</a>)<br>
<br>
Best regards,<br>
PHilipp<br>
<br>
On 2018-02-03 17:36, <a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.gov</a> wrote:<br>
> Hi all,<br>
><br>
> I understand that a call to avg_all computes the running averages of u,<br>
> u^2 etc and then dumps the results after a fixed timestep as given by<br>
> param68.<br>
><br>
> What I do not understand is where in the domain does this computation<br>
> take place? Does it average all the points in the domain? I understand<br>
> how this would by useful in the turbChannel case where fully developed<br>
> turbulence covers essentially the entire domain. However in cases where<br>
> the flow develops turbulence within the domain and thus there are strong<br>
> regions of turbulent anisotropy, it would be imperative to calculate<br>
> these running averages along a specified plane or line. I did see some<br>
> code defining these in the turbChannel example<br>
> viz. common /plane/ uavg_pl(ly1*lely)<br>
> $ , urms_pl(ly1*lely)<br>
> $ , vrms_pl(ly1*lely)<br>
> $ , wrms_pl(ly1*lely)<br>
> $ , uvms_pl(ly1*lely)<br>
> $ , yy(ly1*lely)<br>
> $ , w1(ly1*lely),w2(ly1*lely)<br>
> $ , ffx_avg, dragx_avg<br>
><br>
><br>
> but I'm not exactly sure how this is defining a particular section of<br>
> the geometry. In my particular case I'm trying to simulate a rectangular<br>
> jet exhausting into a cuboidal domain. How can I set avg_all to compute<br>
> the said average along a streamwise plane of the jet for example?<br>
><br>
> Sincerely,<br>
> --<br>
><br>
> *Amitvikram Dutta*<br>
><br>
> Graduate Research Assistant<br>
><br>
> Fluid Mechanics Research Lab<br>
><br>
> Multi-Physics Interaction Lab<br>
><br>
> University of Waterloo<br>
><br>
><br>
><br>
> _______________________________________________<br>
> Nek5000-users mailing list<br>
> <a href="mailto:Nek5000-users@lists.mcs.anl.gov" target="_blank">Nek5000-users@lists.mcs.anl.gov</a><br>
> <a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" rel="noreferrer" target="_blank">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
><br>
_______________________________________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov" target="_blank">Nek5000-users@lists.mcs.anl.gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" rel="noreferrer" target="_blank">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
</blockquote></div>-- <br><div dir="ltr" class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><p dir="ltr"><b><font color="#000033">Amitvikram Dutta</font></b></p>
<p dir="ltr">Graduate Research Assistant</p>
<p dir="ltr">Fluid Mechanics Research Lab</p><p>Multi-Physics Interaction Lab</p>
<p dir="ltr">University of Waterloo</p>
</div></div>