<html dir="ltr">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<style type="text/css" id="owaParaStyle"></style>
</head>
<body fpstyle="1" ocsi="0">
<div style="direction: ltr;font-family: Tahoma;color: #000000;font-size: 10pt;">
<div>Hi Sandeep,</div>
<div><br>
</div>
My guess is that prepost/outpost is not writing the vector field in the solid elements.
<div><br>
</div>
<div>Using outpost in this manner, which in fact is what I would suggest, the code does</div>
<div>not realize that you are writing a field that should be nonzero in the solid.</div>
<div><br>
</div>
<div>I'll try to think of a work around.   One quick one that might work:</div>
<div><br>
</div>
<div>      nelf = nelv  ! Save number of fluid elements</div>
<div>      nelfld(1) = nelt</div>
<div>      nelv = nelt</div>
<div>      call outpost(tx,ty,tz,pr,t,'gdt')</div>
<div>      nelv = nelf</div>
<div>      nelfld(1) = nelf</div>
<div><br>
</div>
<div><br>
</div>
<div>... this might be sufficient to trigger all the flags that prepost may look for.</div>
<div><br>
</div>
<div>Paul</div>
<div><br>
<div style="font-family: Times New Roman; color: #000000; font-size: 16px">
<hr tabindex="-1">
<div id="divRpF286236" 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, January 05, 2017 10:40 AM<br>
<b>To:</b> nek5000-users@lists.mcs.anl.gov<br>
<b>Subject:</b> [Nek5000-users] gradient of temperature in conj_ht<br>
</font><br>
</div>
<div></div>
<div>
<div dir="ltr">
<div class="gmail_default" style="font-size:small">Hello all,</div>
<div class="gmail_default" style="font-size:small">I computed the gradient of temperature in conjugate heat transfer case. The gradient of temperature is being computed correctly for the fluid domain (I validated analytically). But, in the solid domain grad
 T is showing to be zero, which is incorrect. Can anyone please tell me how to compute gradient of a scalar in solid domain for conj_ht case. </div>
<div class="gmail_default" style="font-size:small">Is it possible to compute the gradient of a scalar in this case? The scalar is C0 and not differentiable at the interface between solid and fluid domain.</div>
<div class="gmail_default" style="font-size:small"><br>
</div>
<div class="gmail_default" style="font-size:small">I added the following in the userchk subroutine</div>
<div class="gmail_default" style="font-size:small">
<div class="gmail_default">      parameter (lt=lx1*ly1*lz1*lelt)</div>
<div class="gmail_default">      common /mystuff/ dtdx(lt),dtdy(lt),dtdz(lt)</div>
<div class="gmail_default">
<div class="gmail_default">      common /myvort/ cgtx(lt),cgty(lt),cgtz3(lt)</div>
<div><br>
</div>
<div>       call gradm1(dtdx,dtdy,dtdz,t)<br>
</div>
</div>
</div>
<div class="gmail_default" style="font-size:small">
<div class="gmail_default">
<div class="gmail_default">       call outpost(dtdx,dtdy,dtdz,pr,t,'gdt') ! write to gdtblah.f000...</div>
<div class="gmail_default">
<div class="gmail_default">       call comp_vort3(cgtx,cgty,cgtz,dtdx,dtdy,dtdz)</div>
<div><br>
</div>
</div>
<div><br>
</div>
</div>
<div class="gmail_default">I attach usr, rea and box files for reference.</div>
<div class="gmail_default">Thank you in advance.</div>
<div class="gmail_default"><br>
</div>
<div class="gmail_default">Yours sincerely,</div>
<div class="gmail_default">Sandeep</div>
</div>
</div>
</div>
</div>
</div>
</div>
</body>
</html>