[Nek5000-users] Variable density flows

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

Hi Praveen,

That looks right... As always, it's best to run a test
to verify first.

Paul


----- Original Message -----
From: nek5000-users at lists.mcs.anl.gov
To: nek5000-users at lists.mcs.anl.gov
Sent: Saturday, October 12, 2013 11:17:06 AM
Subject: Re: [Nek5000-users] Variable density flows



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 
_______________________________________________
Nek5000-users mailing list
Nek5000-users at lists.mcs.anl.gov
https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users