<html dir="ltr">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=Windows-1252">
<style>
<!--
@font-face
{font-family:Calibri}
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0in;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri","sans-serif"}
a:link, span.MsoHyperlink
{color:blue;
text-decoration:underline}
a:visited, span.MsoHyperlinkFollowed
{color:purple;
text-decoration:underline}
span.EmailStyle17
{font-family:"Calibri","sans-serif";
color:windowtext}
.MsoChpDefault
{font-family:"Calibri","sans-serif"}
@page WordSection1
{margin:1.0in 1.0in 1.0in 1.0in}
-->
</style><style type="text/css" id="owaParaStyle"></style>
</head>
<body lang="EN-US" link="blue" vlink="purple" fpstyle="1" ocsi="0">
<div style="direction: ltr;font-family: Tahoma;color: #000000;font-size: 10pt;">
<div><br>
</div>
Hi Padyumna,
<div><br>
</div>
<div>I just checked that gradm1() operates on all elements (1:nelt) </div>
<div><br>
</div>
<div>If you want a field that is zero in the solid elements, here is one approach:</div>
<div><br>
</div>
<div><br>
</div>
<div> real tm1(lx1,ly1,lz1,lelt,ldimt) [ note, use lelt here ]</div>
<div><br>
</div>
<div> nv = nx1*ny1*nz1*nelv</div>
<div> nt = nx1*ny1*nz1*nelt</div>
<div><br>
</div>
<div> do i=1,ldimt</div>
<div> call rzero(tm1(1,1,1,1,i),nt)</div>
<div> call copy (tm1(1,1,1,1,i),t(1,1,1,1,i),nv)</div>
<div> enddo</div>
<div><br>
</div>
<div>Note that if you're trying to compute (say) grad V_x * T, you can do so via:</div>
<div><br>
</div>
<div> call col2(vxt,vx,t,nt)</div>
<div> call gradm1(vxtx,vxty,vxtz,vxt)</div>
<div><br>
</div>
<div>Note that vxt will be zero in the solid elements because vx is zero there.</div>
<div><br>
</div>
<div>Best, Paul</div>
<div><br>
<div style="font-family: Times New Roman; color: #000000; font-size: 16px">
<hr tabindex="-1">
<div id="divRpF522941" style="direction: ltr;"><font face="Tahoma" size="2" color="#000000"><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> Thursday, May 19, 2016 5:27 AM<br>
<b>To:</b> nek5000-users@lists.mcs.anl.gov<br>
<b>Subject:</b> [Nek5000-users] Temperature field values in mesh 1 only<br>
</font><br>
</div>
<div></div>
<div>
<div class="WordSection1">
<p class="MsoNormal">Hi Neks!</p>
<p class="MsoNormal"> </p>
<p class="MsoNormal">I am running a Conjugate Heat Transfer simulation with a fluid box (velocity mesh) and two solid boxes (all boxes- temperature mesh). I am interested in computing the heat flux budget for this simulation. While doing this, on multiple occasions,
I need to multiply the temperature or passive scalars throughout the fluid domain field with the velocity component values at the respective locations. I also need to calculate the gradients of these variables throughout the fluid domain. For this I use the
sub-routine “<b>gradm1</b>” present in navier5.f . But I believe gradm1 calculates the gradients for field variables in the velocity mesh only. For all these functions, I think it’s best if I could isolate an array with the values of the
<b>temperature and multiple passive scalars in the fluid domain </b>(or velocity mesh) only, which should be a subset of the whole temperature array.</p>
<p class="MsoNormal"><br>
Thus my question is, is there any way to get this? I am using multiple processors which complicates things. I want to create an array that contains the values of the temperature and the passive scalars at all the velocity mesh elements and this I have to do
across all processors. In other words, from the standard NEK array <b>t(lx2,ly2,lz2,lelt,ldimt)</b> which contains the temperature and the passive scalars, I have to create an array like, say,
<b>tm1(lx1,ly1,lz1,lelv,ldimt)</b> which would have these values. I am not sure if I have been able to explain my query clearly enough, but I would greatly appreciate if anybody could help me with this or if you could think of an easier way to calculate the
gradient of the temperature and passive scalars throughout the velocity field. Let me know if you want more clarity with the question, I will try to explain better.</p>
<p class="MsoNormal"> </p>
<p class="MsoNormal">Thanks in advance and cheers!</p>
<p class="MsoNormal"> </p>
<p class="MsoNormal">Pradyumna</p>
</div>
<br>
<br>
</div>
</div>
</div>
</div>
</body>
</html>