<html><head></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; ">Hi Praveen<div><br><div>Filtering removes energy at the small scales (small length scales/high frequency components)</div><div>If you estimate the turbulent kinetic energy removed by the filtering and compare it to the resolved kinetic energy or the total (resolved+subgrid scale)</div><div>then your simulation :</div><div><br></div><div>(1) May be considered DNS ; if the above ratio is small</div><div>(2) LES if it is relatively big</div><div><br></div><div>in (2) the hope is that subgrid scale motion does not affect the scales of interest in a given problem. </div><div><br></div><div>Filtering stabilizes simulations since it provides a way for removing energy at the small unresolved scales. I ran turbulent wake simulations in the past</div><div>where filtering was absolutely necessary and the argument was that it did not affect the scales of interest.</div><div><br></div><div>So the bottomline is you need to at least be aware of what scales are impacted by filtering and if indeed you are interested in such scales. That is why</div><div>you see some comparisons in the paper for a given case with different values of the filtering parameter. </div><div><br></div><div>so let me know what happens when 32^3 is run with f=0 and also the comparison to Brachet et al. </div><div><br></div><div>Ammar</div></div><div><br></div><div><br></div><div><br></div><div><br></div><div><br></div><div><div>On Dec 8, 2012, at 10:31 PM, <a href="mailto:nek5000-users@lists.mcs.anl.gov">nek5000-users@lists.mcs.anl.gov</a> wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite">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><sup>a1</sup><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"><sup>p1</sup></a></b>, <b><b>Daniel I. </b> <b>Meiron </b><sup>a1</sup></b>, <b><b>Steven A. </b> <b>Orszag </b><sup>a1</sup></b>, <b><b>B. G. </b> <b>Nickel </b><sup>a2</sup></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>
_______________________________________________<br>Nek5000-users mailing list<br><a href="mailto:Nek5000-users@lists.mcs.anl.gov">Nek5000-users@lists.mcs.anl.gov</a><br>https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users<br></blockquote></div><br></body></html>