[Nek5000-users] Setting viscosity in usr file

nek5000-users at lists.mcs.anl.gov nek5000-users at lists.mcs.anl.gov
Thu Apr 24 10:00:19 CDT 2014


On Thu, Apr 24, 2014 at 12:07 PM, <nek5000-users at lists.mcs.anl.gov> wrote:

> Hello Praveen,
>
>
>
> Try to set parameter 'udiff' for field 1 ('ifield.eq.1', i.e. velocity
> field) in subroutine uservp.
>
> Param(30) in .rea file has to be equal to 1..
>
>
>
> Regards,
>
> Donato
>

Thanks Donato. I know this method. My viscosity is a constant, so I want to
set it just once. I have done it in usrdat2 and it seems to work, but
wanted to confirm that it is the right place to set a constant viscosity. I
do the following in usrdat2

param(2) = 1.0/1000.0
cpfld(1,1) = param(2)

Best
praveen
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20140424/259495fd/attachment-0001.html>


More information about the Nek5000-users mailing list