<div dir="ltr"><div><div>Hi Marco,<br><br></div>For the gradient model, you can use gradm1() subroutine to calculate the gradient of vx, vy, and vz individually and then the implementation is straightforward.<br></div><div>For the turbChannel example check for the grid (genbox) and the Reynolds number used, such that it matches the paper., param(2) in negative is basically the inverse of dynamic viscosity. So, a properly defined scale, to generate the Reynolds number relationship for you. Eg, if your bulk flow streamwise velocity is 1, and the channel half width is 1, then Re_Delta = param(2), (Reynolds number based on half channel width). The relationship between Re_tau and Re_delta, can be found from Dean's correlation, that is provided in the seminal paper of Kim and Moin (1987) for DNS of Re_tau = 180.<br></div></div><div class="gmail_extra"><br clear="all"><div><div class="gmail_signature" data-smartmail="gmail_signature"><div dir="ltr"><div>Best Regards,<br></div>Tanmoy<br></div></div></div>
<br><div class="gmail_quote">On Tue, Mar 21, 2017 at 12:25 PM,  <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">Hello,<br>
<br>
I am trying to reproduce some of the turbChannel results that were published with Nek5000 in the paper that can be found in the following link:<br>
<a href="https://arxiv.org/pdf/math/0205264.pdf" rel="noreferrer" target="_blank">https://arxiv.org/pdf/math/<wbr>0205264.pdf</a><br>
<br>
In this paper, two turbulent models are tested: the gradient model and the static Smagorinsky model. The paper also shows results from a coarse DNS for comparison.<br>
<br>
The implementation of the static Smagorinsky model is easily done when starting from the turbChannel example but I am not sure how to implement the gradient model since it involves a tensor product \grad u \dot \grad u. Can anyone point me to an example that has a similar implementation or explain how to implement the gradient model in Nek5000?<br>
<br>
Is it correct to assume that the current setup of the turbChannel example in Nek5000 should allow to reproduce the coarse DNS results presented in the paper?<br>
<br>
What is the relationship between Re_tau and the Reynolds number (param(2)) set in the rea file?<br>
<br>
Any help will be appreciated.<br>
<br>
Thanks,<br>
Marco<br>______________________________<wbr>_________________<br>
Nek5000-users mailing list<br>
<a href="mailto:Nek5000-users@lists.mcs.anl.gov">Nek5000-users@lists.mcs.anl.<wbr>gov</a><br>
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" rel="noreferrer" target="_blank">https://lists.mcs.anl.gov/<wbr>mailman/listinfo/nek5000-users</a><br>
<br></blockquote></div><br></div>