<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=Windows-1252">
<meta name="Generator" content="Microsoft Exchange Server">
<!-- converted from text --><style><!-- .EmailQuote { margin-left: 1pt; padding-left: 4pt; border-left: #800000 2px solid; } --></style>
</head>
<body>
<meta content="text/html; charset=UTF-8">
<style type="text/css" style="">
<!--
p
{margin-top:0;
margin-bottom:0}
-->
</style>
<div dir="ltr">
<div id="x_divtagdefaultwrapper" dir="ltr" style="font-size:12pt; color:#000000; font-family:Calibri,Helvetica,sans-serif">
<p><br>
</p>
<p>Hi Kiran,</p>
<p><br>
</p>
<p>I always like to test code on something I know the exact answer to... so this would be one such case.</p>
<p><br>
</p>
<p>Paul</p>
<p><br>
</p>
</div>
<hr tabindex="-1" style="display:inline-block; width:98%">
<div id="x_divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" color="#000000" style="font-size:11pt"><b>From:</b> nek5000-users-bounces@lists.mcs.anl.gov <nek5000-users-bounces@lists.mcs.anl.gov> on behalf of nek5000-users@lists.mcs.anl.gov <nek5000-users@lists.mcs.anl.gov><br>
<b>Sent:</b> Friday, June 23, 2017 6:06:02 PM<br>
<b>To:</b> nek5000-users@lists.mcs.anl.gov<br>
<b>Subject:</b> Re: [Nek5000-users] Help: compute dissipation</font>
<div> </div>
</div>
</div>
<font size="2"><span style="font-size:10pt;">
<div class="PlainText">Hi Paul,<br>
<br>
No, I haven’t tried with Poiseiulle flow. Any reason as to why I should test it on that flow?<br>
<br>
Thank,<br>
Kiran<br>
<br>
<br>
<br>
> On 23 Jun 2017, at 16:21, nek5000-users-request@lists.mcs.anl.gov wrote:<br>
> <br>
> Send Nek5000-users mailing list submissions to<br>
> nek5000-users@lists.mcs.anl.gov<br>
> <br>
> To subscribe or unsubscribe via the World Wide Web, visit<br>
> <a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
> or, via email, send a message with subject or body 'help' to<br>
> nek5000-users-request@lists.mcs.anl.gov<br>
> <br>
> You can reach the person managing the list at<br>
> nek5000-users-owner@lists.mcs.anl.gov<br>
> <br>
> When replying, please edit your Subject line so it is more specific<br>
> than "Re: Contents of Nek5000-users digest..."<br>
> <br>
> <br>
> Today's Topics:<br>
> <br>
> 1. Re: Help: compute dissipation (nek5000-users@lists.mcs.anl.gov)<br>
> <br>
> <br>
> ----------------------------------------------------------------------<br>
> <br>
> Message: 1<br>
> Date: Fri, 23 Jun 2017 15:17:51 +0000<br>
> From: nek5000-users@lists.mcs.anl.gov<br>
> To: "nek5000-users@lists.mcs.anl.gov"<br>
> <nek5000-users@lists.mcs.anl.gov><br>
> Subject: Re: [Nek5000-users] Help: compute dissipation<br>
> Message-ID:<br>
> <mailman.9487.1498231285.2967.nek5000-users@lists.mcs.anl.gov><br>
> Content-Type: text/plain; charset="windows-1252"<br>
> <br>
> <br>
> Hi,<br>
> <br>
> <br>
> Have you tested your routine on plane Poiseiulle flow,<br>
> <br>
> u=1-y^2 ?<br>
> <br>
> <br>
> Paul<br>
> <br>
> ________________________________<br>
> From: nek5000-users-bounces@lists.mcs.anl.gov <nek5000-users-bounces@lists.mcs.anl.gov> on behalf of nek5000-users@lists.mcs.anl.gov <nek5000-users@lists.mcs.anl.gov><br>
> Sent: Friday, June 23, 2017 8:12:37 AM<br>
> To: nek5000-users@lists.mcs.anl.gov<br>
> Subject: [Nek5000-users] Help: compute dissipation<br>
> <br>
> Hi NEK users,<br>
> <br>
> I am trying to calculate dissipation for fully developed turbulent channel flow. I tried the code below, but result I got from this totally disagrees with Vreman?s 2014 paper. x-axis is streamwise, y-axis is wall-normal and z-axis is spanwise. Any errors
in the code? Any thoughts?<br>
> <br>
> real diss_sqr(lx1,ly1,lz1, lelt)<br>
> real s11<br>
> real s12<br>
> real s13<br>
> real s22<br>
> real s23<br>
> real s33<br>
> <br>
> <br>
> common /plane/ uavg_pl(ly1*lely)<br>
> $ , vavg_pl(ly1*lely)<br>
> $ , wavg_pl(ly1*lely)<br>
> $ , urms_pl(ly1*lely)<br>
> $ , vrms_pl(ly1*lely)<br>
> $ , wrms_pl(ly1*lely)<br>
> $ , uvms_pl(ly1*lely)<br>
> $ , diss_sqr_pl(ly1*lely)<br>
> $ , yy(ly1*lely)<br>
> $ , w1(ly1*lely),w2(ly1*lely)<br>
> $ , ffx_avg, dragx_avg<br>
> <br>
> <br>
> parameter (lr=lx1*ly1*lz1)<br>
> common /scruz/ ur(lr),us(lr),ut(lr)<br>
> $ , vr(lr),vs(lr),vt(lr)<br>
> $ , wr(lr),ws(lr),wt(lr)<br>
> <br>
> <br>
> common /avg/ uavg(lx1,ly1,lz1,lelv)<br>
> & , vavg(lx1,ly1,lz1,lelv)<br>
> & , wavg(lx1,ly1,lz1,lelv)<br>
> & , urms(lx1,ly1,lz1,lelv)<br>
> & , vrms(lx1,ly1,lz1,lelv)<br>
> & , wrms(lx1,ly1,lz1,lelv)<br>
> & , uvms(lx1,ly1,lz1,lelv)<br>
> & , dissavg(lx1,ly1,lz1,lelv)<br>
> <br>
> <br>
> common /scrns/ sij (lx1*ly1*lz1,6,lelv)<br>
> <br>
> if(icalld.eq.0) then<br>
> call rzero(uavg,n)<br>
> call rzero(urms,n)<br>
> call rzero(vrms,n)<br>
> call rzero(wrms,n)<br>
> call rzero(uvms,n)<br>
> call rzero(dissavg,n)<br>
> call rzero(wavg,n)<br>
> call rzero(vavg,n)<br>
> call rzero(diss_sqr,n)<br>
> <br>
> <br>
> call comp_sij(sij,6,vx,vy,vz,ur,us,ut,vr,vs,vt,wr,ws,wt)<br>
> <br>
> do e=1,nelt<br>
> do i = 1,nx1*ny1*nz1<br>
> s11=sij(i,1,e)<br>
> s12=sij(i,4,e)<br>
> s22=sij(i,2,e)<br>
> s13=sij(i,6,e)<br>
> s23=sij(i,5,e)<br>
> s33=sij(i,3,e)<br>
> <br>
> diss_sqr(i,1,1,e) = (s11*s11) +<br>
> $ (s22*s22) +<br>
> $ (s33*s33) +<br>
> $ (s12*s12)*2 + (s23*s23)*2 +<br>
> $ (s13*s13)*2<br>
> <br>
> enddo<br>
> enddo<br>
> <br>
> call avg1(dissavg,diss_sqr,alpha,beta,n,'diss',ifverbose)<br>
> <br>
> call planar_average_s(diss_sqr_pl,dissavg,w1,w2)<br>
> <br>
> do i=1,ny1*nely/2<br>
> uavg_pl(i) = 0.5 * (uavg_pl(i) + uavg_pl(m-i+1))<br>
> vavg_pl(i) = 0.5 * (vavg_pl(i) + vavg_pl(m-i+1))<br>
> wavg_pl(i) = 0.5 * (wavg_pl(i) + wavg_pl(m-i+1))<br>
> diss_sqr_pl(i) = 0.5 * (diss_sqr_pl(i) + diss_sqr_pl(m-i+1))<br>
> urms_pl(i) = 0.5 * (urms_pl(i) + urms_pl(m-i+1))<br>
> vrms_pl(i) = 0.5 * (vrms_pl(i) + vrms_pl(m-i+1))<br>
> wrms_pl(i) = 0.5 * (wrms_pl(i) + wrms_pl(m-i+1))<br>
> <br>
> enddo<br>
> <br>
> <br>
> <br>
> Many thanks.<br>
> <br>
> <br>
> -------------- next part --------------<br>
> An HTML attachment was scrubbed...<br>
> URL: <<a href="http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20170623/c4bae54c/attachment.html">http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20170623/c4bae54c/attachment.html</a>><br>
> <br>
> ------------------------------<br>
> <br>
> _______________________________________________<br>
> Nek5000-users mailing list<br>
> Nek5000-users@lists.mcs.anl.gov<br>
> <a href="https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users">https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users</a><br>
> <br>
> <br>
> End of Nek5000-users Digest, Vol 100, Issue 27<br>
> **********************************************<br>
<br>
_______________________________________________<br>
Nek5000-users mailing list<br>
Nek5000-users@lists.mcs.anl.gov<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></font>
</body>
</html>