Hello Ammar<div><br></div><div>The filter refers to a stabilization technique and is not related to LES. It is explained here</div><div><br></div><div><a href="http://www-unix.mcs.anl.gov/~fischer/pubhtml/filter01.pdf"> "Filter-Based Stabilization of Spectral Element Methods" </a></div>
<div><br></div><div>I do not know the details of the spectral simulation. I took the data from the website I mentioned in my post and some explanation is given in the case description I linked to before.</div><div><br></div>
<div>I used 14^3 and 28^3 elements but within each element the velocity is a Q_9 function and pressure is Q_7 function.</div><div><br></div><div>Thanks for the reference. I have not seen it before and will check it out.</div>
<div><br></div><div>Thanks</div><div>praveen<br><br><div class="gmail_quote">On Sun, Dec 9, 2012 at 2:42 AM,  <span dir="ltr"><<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.gov</a>></span> wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word">Can you elaborate on the 512^3 case? it is impressive you are getting close with 32^3 and some filtering<div>
I was wondering :</div><div>(1) what if you turn off filtering with 32^3  (i.e. set f=0), the question here if filtering is significantly contributing to the dissipation rate. if it does not make a big difference then I would say your LES with 32^3 and f=0.05 is pretty good. </div>
<div>(2) how does the time history in all cases compare to Brachet et al. (see my previous note) at the same Re. </div><div><br></div><div>Ammar</div><div><div class="h5"><div><br><div><div>On Dec 8, 2012, at 3:01 PM, <a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.gov</a> wrote:</div>
<br><blockquote type="cite"><div style="word-wrap:break-word">Your results make sense to me.  Note that you are not really "capturing" peak dissipation rate with filtering. <div>you are only supplementing the resolved dissipation on a given grid by artificially damping energy.   that's why you get higher peak with f=0.05 and N=14^3 relative to  f=0  and N=14^3 </div>
<div>I wonder how much resolution N^3 you need (with f=0) to truly capture the peak dissipation rate?  </div><div><br></div><div>In order to better judge your simulations you should compare DNS (f=0) and LES f=nonzero for gradually increased resolution N^3 </div>
<div><br></div><div>if DNS is expensive , good news DNS is already been done , Check Brachet et al JFM article  <a href="http://journals.cambridge.org/action/displayJournal?jid=FLM" title="Journal of Fluid Mechanics" style="margin:0px;padding:0px;border:0px;outline:0px;font-weight:normal;font-style:normal;font-size:11px;font-family:'Arial Unicode MS',Arial,Helvetica,sans-serif;vertical-align:baseline;color:rgb(4,89,137);text-decoration:none;line-height:14px;font-variant:normal;letter-spacing:normal;text-align:-webkit-auto;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;background-color:rgb(241,241,241)" target="_blank">Journal of Fluid Mechanics</a><span style="text-indent:0px;letter-spacing:normal;font-variant:normal;text-align:-webkit-auto;font-style:normal;background-color:rgb(241,241,241);display:inline!important;font-weight:normal;float:none;line-height:14px;text-transform:none;font-size:11px;white-space:normal;font-family:'Arial Unicode MS',Arial,Helvetica,sans-serif;word-spacing:0px"> / Volume 130 / </span><span style="text-indent:0px;letter-spacing:normal;font-variant:normal;text-align:-webkit-auto;font-style:normal;background-color:rgb(241,241,241);display:inline!important;font-weight:normal;float:none;line-height:14px;text-transform:none;font-size:11px;white-space:normal;font-family:'Arial Unicode MS',Arial,Helvetica,sans-serif;word-spacing:0px">May 1983,    pp 411-452</span></div>
<div>(</div><h2 style="margin:0px;padding:5px 0px;border:0px;outline:0px;font-weight:bold;font-style:normal;font-size:11px;font-family:'Arial Unicode MS',Arial,Helvetica,sans-serif;vertical-align:baseline;color:rgb(4,89,137);font-variant:normal;letter-spacing:normal;line-height:15px;text-align:left;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;background-color:rgb(241,241,241)">
<font color="#336699" size="+1" style="margin:0px;padding:0px;border:0px;outline:0px;font-weight:inherit;font-style:inherit;font-size:11px;font-family:inherit;vertical-align:baseline">Small-scale structure of the Taylor–Green vortex</font></h2>
<h3 style="margin-top:0px;margin-right:0px;margin-bottom:0px;margin-left:0px;padding-top:5px;padding-right:0px;padding-bottom:5px;padding-left:0px;border-top-width:0px;border-right-width:0px;border-bottom-width:0px;border-left-width:0px;border-style:initial;border-color:initial;outline-width:0px;outline-style:initial;outline-color:initial;font-weight:bold;font-style:normal;font-size:11px;font-family:verdana,arial,sans-serif;vertical-align:baseline;color:rgb(98,98,98);font-variant:normal;letter-spacing:normal;line-height:15px;text-align:left;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;background-color:rgb(241,241,241)">
<table border="0" align="top" width="95%" style="margin-top:10px;margin-right:0px;margin-bottom:0px;margin-left:0px;padding-top:0px;padding-right:0px;padding-bottom:0px;padding-left:0px;border-top-width:0px;border-right-width:0px;border-bottom-width:0px;border-left-width:0px;border-style:initial;border-color:initial;outline-width:0px;outline-style:initial;outline-color:initial;font-weight:normal;font-style:normal;font-size:11px;font-family:'Arial Unicode MS',Arial,Helvetica,sans-serif;vertical-align:baseline;border-collapse:separate;color:rgb(98,98,98);font-variant:normal;letter-spacing:normal;line-height:15px;text-align:left;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px;background-color:rgb(241,241,241)">
<tbody style="margin:0px;padding:0px;border:0px;outline:0px;font-weight:inherit;font-style:inherit;font-size:11px;font-family:inherit;vertical-align:baseline"><tr valign="top" style="margin:0px;padding:0px;border:0px;outline:0px;font-weight:inherit;font-style:inherit;font-size:11px;font-family:inherit;vertical-align:baseline">
<td style="margin:0px;padding:0px;border:0px;outline:0px;font-weight:normal;font-style:inherit;font-size:11px;font-family:'Arial Unicode MS',Arial,Helvetica,sans-serif;vertical-align:baseline;text-align:left;color:rgb(98,98,98);line-height:1.4em">
<b><b>Marc E. </b> <b>Brachet </b>^a1<a href="http://journals.cambridge.org/action/displayAbstract;jsessionid=2541E59358A427C052F22128A3D13BCC.journals?fromPage=online&aid=376272#p1" style="margin:0px;padding:0px;border:0px;outline:0px;font-weight:inherit;font-style:inherit;font-size:11px;font-family:inherit;vertical-align:baseline;color:rgb(4,89,137);text-decoration:none" target="_blank">^p1</a></b>, <b><b>Daniel I. </b> <b>Meiron </b>^a1</b>, <b><b>Steven A. </b> <b>Orszag </b>^a1</b>, <b><b>B. G. </b> <b>Nickel </b>^a2</b></td>
</tr></tbody></table></h3><div>)</div><div><br></div><div><br></div><div> and compare the time history of the dissipation rate</div><div>to what you get at the same Reynolds number. </div><div><br><div><br></div><div><br>
</div><div>Ammar</div><div><br></div><div><br></div><div><div><br></div><div><br></div><div><br></div><div><br><div><div>On Dec 8, 2012, at 2:36 AM, <a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank">nek5000-users@lists.mcs.anl.gov</a> wrote:</div>
<br><blockquote type="cite">Hello<div>I am solving the 3-D Taylor-Green problem as described here </div><div><br></div><div><a href="http://www.dlr.de/as/desktopdefault.aspx/tabid-8170/13999_read-35550/" target="_blank">http://www.dlr.de/as/desktopdefault.aspx/tabid-8170/13999_read-35550/</a></div>

<div><a href="http://www.as.dlr.de/hiocfd/case_c3.5.pdf" target="_blank">C3.5 Direct Numerical Simulation of the Taylor-Green Vortex at Re = 1600</a></div><div><br></div><div>using NEK.</div><div><br></div><div>I use</div>
<div><br></div>
<div>lx1 = 10</div><div>lx2 = lx1-2</div><div>ldx = 15</div><div><br></div><div>No of elements = 14^3 or 28^3</div><div><br></div><div>The evolution of dissipation rate with time is shown in attached pdf. </div><div><br>
</div>
<div><div>f=0.0 indicates no filtering and f=0.05 means filter parameter P103.</div></div><div><br></div><div>Use of filter seems to capture peak dissipation rate better but still some difference is seen. The DG results with similar number of dofs in the problem description give good estimate of peak rate. So I expect NEK results to improve further and would like to ask for any advice. Should I try with larger amount of filtering ?</div>

<div><br></div><div>I have attached my SIZE and rea files.</div><div><br></div><div>Thanks</div><div>praveen</div>
<span><box.rea></span><span><diss.pdf></span><span><SIZE></span>_______________________________________________<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" target="_blank">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br></blockquote></div><br></div></div></div></div>_______________________________________________<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" target="_blank">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
</blockquote></div><br></div></div></div></div><br>_______________________________________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov">Nek5000-users@lists.mcs.anl.gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" target="_blank">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
<br></blockquote></div><br></div>