[Nek5000-users] gradient of temperature in conj_ht

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

Hi Sandeep,

My guess is that prepost/outpost is not writing the vector field in the solid elements.

Using outpost in this manner, which in fact is what I would suggest, the code does
not realize that you are writing a field that should be nonzero in the solid.

I'll try to think of a work around.   One quick one that might work:

      nelf = nelv  ! Save number of fluid elements
      nelfld(1) = nelt
      nelv = nelt
      call outpost(tx,ty,tz,pr,t,'gdt')
      nelv = nelf
      nelfld(1) = nelf


... this might be sufficient to trigger all the flags that prepost may look for.

Paul

________________________________
From: nek5000-users-bounces at lists.mcs.anl.gov [nek5000-users-bounces at lists.mcs.anl.gov] on behalf of nek5000-users at lists.mcs.anl.gov [nek5000-users at lists.mcs.anl.gov]
Sent: Thursday, January 05, 2017 10:40 AM
To: nek5000-users at lists.mcs.anl.gov
Subject: [Nek5000-users] gradient of temperature in conj_ht

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/4ca179aa/attachment.html>