<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" style="display:none;"><!-- P {margin-top:0;margin-bottom:0;} --></style>
</head>
<body dir="ltr">
<div id="divtagdefaultwrapper" style="font-size:12pt;color:#000000;font-family:Calibri,Helvetica,sans-serif;" dir="ltr">
<p><br>
</p>
<p>Dear Samuel,</p>
<p><br>
</p>
<p>Several of these routines are in the turbchannel example .usr file (and nek5000 source). Please have a look at the .usr file in that example and you'll see the below.</p>
<p><br>
</p>
<p>hth, Paul</p>
<p><br>
</p>
<p><br>
</p>
<p style="margin-right: 0px; margin-left: 0px; font-size: 13px; line-height: normal; font-family: Courier;">
call comp_gije(sij,vx(1,1,1,e),vy(1,1,1,e),vz(1,1,1,e),e)</p>
<p style="margin-right: 0px; margin-left: 0px; font-size: 13px; line-height: normal; font-family: Courier;">
call comp_sije(sij)</p>
<p style="margin-right: 0px; margin-left: 0px; font-size: 13px; line-height: normal; font-family: Courier; min-height: 16px;">
<br>
</p>
<p style="margin-right: 0px; margin-left: 0px; font-size: 13px; line-height: normal; font-family: Courier;">
call mag_tensor_e(snrm(1,e),sij)</p>
<p></p>
<p style="margin-right: 0px; margin-left: 0px; font-size: 13px; line-height: normal; font-family: Courier;">
call cmult(snrm(1,e),2.0,ntot)</p>
<p style="margin-right: 0px; margin-left: 0px; font-size: 13px; line-height: normal; font-family: Courier;">
<br>
</p>
<p style="margin-right: 0px; margin-left: 0px; font-size: 13px; line-height: normal; font-family: Courier;">
<br>
</p>
</div>
<hr style="display:inline-block;width:98%" tabindex="-1">
<div id="divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" style="font-size:11pt" color="#000000"><b>From:</b> Nek5000-users <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, November 10, 2017 8:56:09 AM<br>
<b>To:</b> nek5000-users@lists.mcs.anl.gov<br>
<b>Subject:</b> Re: [Nek5000-users] Dissipation rate</font>
<div> </div>
</div>
<div class="BodyFragment"><font size="2"><span style="font-size:10pt;">
<div class="PlainText">Dear Nek experts,<br>
<br>
I'm still working on the turbulent dissipation rate.<br>
I would like to compute epsilon=2*nu*<Sij Sij> with <br>
Sij=0.5*(grad(U)+grad(U)^T)<br>
<br>
I found in the forum a proposition to calculate the scalar dissipation <br>
rate as follow.<br>
<br>
If I understood well, this code is basically doing "chi := D * |grad <br>
Z|^2" while I would like to compute chi := D * D * |grad Z + grad^T Z|^2<br>
<br>
Does anyone have a suggestion ?<br>
<br>
Many thanks,<br>
Samuel<br>
<br>
c-----------------------------------------------------------------------<br>
subroutine magSqr(a,b1,b2,b3,n)<br>
include 'SIZE'<br>
real a(1)<br>
real b1(1),b2(1),b3(1)<br>
<br>
do i=1,n<br>
a(i) = b1(i)*b1(i) + b2(i)*b2(i)<br>
enddo<br>
<br>
return<br>
end<br>
c-----------------------------------------------------------------------<br>
subroutine scalDisp(chi,Z,D)<br>
c<br>
c compute scalar dissipation rate<br>
c chi := D * |grad Z|^2<br>
c<br>
include 'SIZE'<br>
include 'TOTAL'<br>
<br>
real chi(lx1,ly1,lz1,1)<br>
real Z (lx1,ly1,lz1,1)<br>
real D (lx1,ly1,lz1,1)<br>
<br>
common /scrns/ w1(lx1,ly1,lz1,lelt)<br>
$ ,w2(lx1,ly1,lz1,lelt)<br>
$ ,w3(lx1,ly1,lz1,lelt)<br>
<br>
ntot = nx1*ny1*nz1*nelv<br>
<br>
call opgrad (w1,w2,w3,Z)<br>
call opdssum(w1,w2,w3)<br>
call opcolv (w1,w2,w3,binvm1)<br>
<br>
call magsqr (chi,w1,w2,w3,ntot)<br>
call col2 (chi,D,ntot)<br>
<br>
return<br>
end<br>
c-----------------------------------------------------------------------<br>
<br>
<br>
Le 2017-11-08 10:55, Samuel Ahizi a écrit :<br>
> Dear Neks,<br>
> <br>
> I would like to compute the turbulent dissipation rate, in order to<br>
> calculate the Kolmogorov microscales. Does a subroutine already exist<br>
> or should I compute it from userchk ?<br>
> <br>
> Thank you very much for your answers,<br>
> Samuel<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></div>
</body>
</html>