<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
</head>
<body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">
Hi Tanmoy,
<div class=""><br class="">
</div>
<div class="">thank you for your reply. I did find the correlation in the reference you provided.</div>
<div class=""><br class="">
</div>
<div class="">I have one more question regarding the implementation of the gradient model:</div>
<div class="">the gradient model consists of computing \tau defined as proportional to the tensor product \grad u \dot \grad u. Then, once \tau is computed, \grad \tau is added to the momentum equation as a source term. Should I add \grad \tau to the momentum
equation using userf function in the user file? </div>
<div class=""><br class="">
</div>
<div class="">Thanks,</div>
<div class="">Marco</div>
<div class=""><br class="">
</div>
<div class=""><br class="">
<div>
<blockquote type="cite" class="">
<div class="">On Mar 21, 2017, at 9:07 PM, <a href="mailto:nek5000-users@lists.mcs.anl.gov" class="">
nek5000-users@lists.mcs.anl.gov</a> wrote:</div>
<br class="Apple-interchange-newline">
<div class="">
<div dir="ltr" class="">
<div class="">
<div class="">Hi Marco,<br class="">
<br class="">
</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 class="">
</div>
<div class="">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 class="">
</div>
</div>
<div class="gmail_extra"><br clear="all" class="">
<div class="">
<div class="gmail_signature" data-smartmail="gmail_signature">
<div dir="ltr" class="">
<div class="">Best Regards,<br class="">
</div>
Tanmoy<br class="">
</div>
</div>
</div>
<br class="">
<div class="gmail_quote">On Tue, Mar 21, 2017 at 12:25 PM, <span dir="ltr" class="">
<<a href="mailto:nek5000-users@lists.mcs.anl.gov" target="_blank" class="">nek5000-users@lists.mcs.anl.gov</a>></span> wrote:<br class="">
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hello,<br class="">
<br class="">
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 class="">
<a href="https://arxiv.org/pdf/math/0205264.pdf" rel="noreferrer" target="_blank" class="">https://arxiv.org/pdf/math/<wbr class="">0205264.pdf</a><br class="">
<br class="">
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 class="">
<br class="">
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 class="">
<br class="">
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 class="">
<br class="">
What is the relationship between Re_tau and the Reynolds number (param(2)) set in the rea file?<br class="">
<br class="">
Any help will be appreciated.<br class="">
<br class="">
Thanks,<br class="">
Marco<br class="">
______________________________<wbr class="">_________________<br class="">
Nek5000-users mailing list<br class="">
<a href="mailto:Nek5000-users@lists.mcs.anl.gov" class="">Nek5000-users@lists.mcs.anl.<wbr class="">gov</a><br class="">
<a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users" rel="noreferrer" target="_blank" class="">https://lists.mcs.anl.gov/<wbr class="">mailman/listinfo/nek5000-users</a><br class="">
<br class="">
</blockquote>
</div>
<br class="">
</div>
_______________________________________________<br class="">
Nek5000-users mailing list<br class="">
<a href="mailto:Nek5000-users@lists.mcs.anl.gov" class="">Nek5000-users@lists.mcs.anl.gov</a><br class="">
https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users<br class="">
</div>
</blockquote>
</div>
<br class="">
</div>
</body>
</html>