[Nek5000-users] Variable density flows

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

Thanks to Paul and Ammar. I will try the variable density approach.

I am first trying a case in which viscosity depends on temperature and a
passive scalar. My uservp looks like this. Can I use "temp" and "PS" to get
temperature and passive scalar as below ?

      subroutine uservp (ix,iy,iz,ieg)
      include 'SIZE'
      include 'TOTAL'
      include 'NEKUSE'
      common /mygeom/ xmin,xmax
      common /myparam/ r_k,Rey_no,Pr_no,Sc_no,visc_rat

      if(ifield.eq.1)then   ! momentum equation
         utrans= 1.0
         udiff = exp(PS(1)*log(visc_rat) - temp)/Rey_no
      else if(ifield.eq.2)then ! temperature eqn
         utrans= 1.0
         udiff = PS(1)*(r_k - 1.0) + 1.0
         udiff = udiff/(Rey_no*Pr_no)
      else if(ifield.eq.3)then  ! passive scalar
         utrans= 1.0
         udiff = 1.0/(Rey_no*Sc_no)
      else
         print*,'Only 3 fields'
         call exitt()
      endif

      return
      end

Thanks
praveen
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20131012/85a9d277/attachment.html>