<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>
</head>
<body dir="ltr">
<div id="divtagdefaultwrapper" style="font-size:12pt;color:#000000;font-family:Calibri,Helvetica,sans-serif;" dir="ltr">
<p style="margin-top:0;margin-bottom:0"><br>
</p>
<p style="margin-top:0;margin-bottom:0">Usually I just build a sequence of files with a different name such that the first one has the geometry. This is easily done with a unix script, e.g.:</p>
<p style="margin-top:0;margin-bottom:0"><br>
</p>
<p style="margin-top:0;margin-bottom:0">ln case0.f00001 view0.f00001</p>
<p style="margin-top:0;margin-bottom:0">ln avgcase0.f00001 view0.f00002</p>
<p style="margin-top:0;margin-bottom:0">ln avgcase0.f00002 view0.f00003</p>
<p style="margin-top:0;margin-bottom:0"><br>
</p>
<p style="margin-top:0;margin-bottom:0">etc.</p>
<p style="margin-top:0;margin-bottom:0"><br>
</p>
<p style="margin-top:0;margin-bottom:0">Then use visit to look at case "view"</p>
<p style="margin-top:0;margin-bottom:0"><br>
</p>
<p style="margin-top:0;margin-bottom:0">Paul</p>
<p style="margin-top:0;margin-bottom:0"><br>
</p>
</div>
<hr style="display:inline-block;width:98%" tabindex="-1">
<div id="divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt" color="#000000"><b>From:</b> Nek5000-users <nek5000-users-bounces@lists.mcs.anl.gov> on behalf of nek5000-users@lists.mcs.anl.gov <nek5000-users@lists.mcs.anl.gov><br>
<b>Sent:</b> Sunday, February 4, 2018 2:08:08 PM<br>
<b>To:</b> nek5000-users@lists.mcs.anl.gov<br>
<b>Subject:</b> Re: [Nek5000-users] Questions regarding avg_all</font>
<div> </div>
</div>
<div>
<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="x_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="x_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="x_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>
</div>
</body>
</html>