[Nek5000-users] Heat transfer

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

I assume that for your isothermal case, you set the boundary conditions for the velocity field with the following piece of code:

do iel=1,nelt
   do ifc=1,2*ndim
     idss = bc(5,ifc,iel,1)
     if (idss.eq.1) cbc(ifc,iel,1)='v ‘
     if (idss.eq.2) cbc(ifc,iel,1)=‘O ‘
     if (idss.eq.3) cbc(ifc,iel,1)=‘W '
  enddo
enddo

The ‘1’ in the arrays cbc and bc is for the velocity field. If you use the same approach for the temperature equation, you will have the replace ‘1’ with ‘2’ which corresponds to the temperature field.

Marco

On Jan 23, 2018, at 5:04 PM, nek5000-users at lists.mcs.anl.gov<mailto:nek5000-users at lists.mcs.anl.gov> wrote:

Hi Nek`s

I generated my mesh in exo2nek. My initial case was a DNS study over pebbles in tandem for iso-thermal flow. I now wish to add heat transfer from the pebbles..

How do i specify the thermal boundary conditions for a mesh that has been generated by exo2nek?

Regards



_______________________________________________
Nek5000-users mailing list
Nek5000-users at lists.mcs.anl.gov<mailto:Nek5000-users at lists.mcs.anl.gov>
https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20180124/785b16b6/attachment.html>