
<html><head><base href="x-msg://300/"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; "><div>Paul: you may want to supplement my answer.</div><div><br></div>To give you an idea of what we're doing here:<div><br></div><div>1. construct a low-pass filter transfer function (to preserve the accuracy the transfer function needs to be smooth). The particular shape of the transfer function depends in the application.</div><div>2. transform the Lagrange polynomials into a hierarchical modal basis (e.g. Legendre polynomials)</div><div>3. apply the filer</div><div>4. apply transpose of transformation (see step 2) <br><div><br></div><div>At the moment we support only two options for LES:</div><div><br></div><div>- dynamic Smagorinsky (only tensor product domains are supported for now)</div><div><br></div><div>- simple filter model</div><div>In this model we simply filter the velocities u after every timestep (e.g. only last 2 modes are filtered using a quadratic transfer function. The last mode is filtered by 5%). The filtering operation does not preserve div(u) however this effect is typically very little. I cannot tell you how much and many modes you need to filter because this model lacks of a physical justification. However it is in the spirit of every LES model. All of them try to remove energy which you cannot present by your grid.</div><div><br></div><div>hth,</div><div>Stefan</div><div><br></div><div><br><div><div>On Oct 16, 2009, at 6:46 PM, 王志成 wrote:</div><br class="Apple-interchange-newline"><blockquote type="cite"><span class="Apple-style-span" style="border-collapse: separate; font-family: Helvetica; font-size: medium; font-style: normal; font-variant: normal; font-weight: normal; letter-spacing: normal; line-height: normal; orphans: 2; text-indent: 0px; text-transform: none; white-space: normal; widows: 2; word-spacing: 0px; -webkit-border-horizontal-spacing: 0px; -webkit-border-vertical-spacing: 0px; -webkit-text-decorations-in-effect: none; -webkit-text-size-adjust: auto; -webkit-text-stroke-width: 0px; "><div style="margin-top: 10px; margin-right: 10px; margin-bottom: 10px; margin-left: 10px; font-family: verdana; font-size: 10pt; "><div><font size="2" face="Verdana">Hi Stefan,</font></div><div> </div><div>I'm puzzled with the subroutine build_1d_filt in<span class="Apple-converted-space"> </span><font size="2" face="Verdana">postpro.f. Could you please give some papers on this problem?</font></div><div> </div><div>Regards,</div><div> </div><div>Zhicheng</div><div> </div><div align="left"><font color="#c0c0c0" size="2" face="Verdana">2009-10-17</font></div><font size="2" face="Verdana"><hr align="left" size="2" style="width: 122px; height: 2px; "><div><font color="#c0c0c0" size="2" face="Verdana"><span>王志成</span></font></div></font>_______________________________________________<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">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br></div></span></blockquote></div><br></div></div></body></html>