<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>