[Nek5000-users] gradient of temperature in conj_ht

[nek5000-users at lists.mcs.anl.gov]

Hello all,
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.
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.

I added the following in the userchk subroutine
      parameter (lt=lx1*ly1*lz1*lelt)
      common /mystuff/ dtdx(lt),dtdy(lt),dtdz(lt)
      common /myvort/ cgtx(lt),cgty(lt),cgtz3(lt)

       call gradm1(dtdx,dtdy,dtdz,t)
       call outpost(dtdx,dtdy,dtdz,pr,t,'gdt') ! write to gdtblah.f000...
       call comp_vort3(cgtx,cgty,cgtz,dtdx,dtdy,dtdz)


I attach usr, rea and box files for reference.
Thank you in advance.

Yours sincerely,
Sandeep
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20170105/f07f4179/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: conj_ht.box
Type: application/octet-stream
Size: 640 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20170105/f07f4179/attachment-0003.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: conj_ht.usr
Type: application/octet-stream
Size: 4067 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20170105/f07f4179/attachment-0004.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: conj_ht.rea
Type: application/octet-stream
Size: 130535 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20170105/f07f4179/attachment-0005.obj>