[Nek5000-users] variable properties, uservp

nek5000-users at lists.mcs.anl.gov nek5000-users at lists.mcs.anl.gov
Mon May 3 17:53:27 CDT 2010


Hello all,

I am trying to understand how to set variable fluid properties depending
on spatial location.  I have two cases, for one "uservp" is never
called, for the other it is.  The difference in their *.rea files is as
follows.  I cannot see why in one case "uservp" is called as neither
IFUSERVP=.true. or param(3)>0.  

Cheers,
Frank

[fmuldoo at tetra sim3]# diff
*rea ../sim7/AR-.66-RE-1800-PR-04.0-dt-3.1e-5-time-006.0-grid-1-Intel/delet.rea | more
3c3
<             2 DIMENSIONAL RUN
---
>    3 DIMENSIONAL RUN
6c6
<   -1.00000         VISCOS
---
>    1.00000         VISCOS
12c12
<    1.00000         CONDUCT
---
>    0.25            CONDUCT
15,16c15,16
<   1000.0         NSTEPS
<  -5.e-6            DT
---
>   193548          NSTEPS
>  -3.1e-05            DT
19c19
<   20.00000         IOSTEP
---
>   1290             IOSTEP
25,26c25,26
<   1.000000E-07     DIVERGENCE
<   1.000000E-07     HELMHOLTZ
---
>   1.000000E-05     DIVERGENCE
>   1.000000E-05     HELMHOLTZ
28,29c28,29
<   1.000000E-07     TOLREL
<   1.000000E-07     TOLABS
---
>   1.000000E-05     TOLREL
>   1.000000E-05     TOLABS
39a40,45
>   0.000000E+00     XMAGNET
>   0.000000E+00     NGRIDS
>   0.000000E+00     NORDER2
>   0.000000E+00     NORDER3
>   0.000000E+00     NORDER4
>   0.000000E+00     NORDER5
46,51d51
<   1.000000E+00     p42=1 for std 2-level schwarz
<   1.000000E+00     p43=1 etc.
<   1.000000E+00     p44=1
<   0.000000E+00
<   0.000000E+00
<   0.3              Poisson ratio for mesh
107c107
<    0.05000           p103: filter strength,  < 0 ==> no fil
---
>    0.00000           p103: filter strength,  < 0 ==> no fil
119c119
<   T     IFAXIS
---
>   F     IFAXIS
125,126c125,126
<   T     IFMVBD
<   F     IFCHAR
---
>   F     IFMVBD
>   T     IFCHAR




-- 
Frank Herbert Muldoon, Ph.D. Mechanical Engineering
Technische Universität Wien (Technical University of Vienna)
Inst. f. Strömungsmechanik und Wärmeübertragung (Institute of Fluid
Mechanics and Heat Transfer)
Resselgasse 3
1040 Wien
Tel: +4315880132232
Fax: +4315880132299 
Cell:+436765203470
fmuldoo (skype)
http://tetra.fluid.tuwien.ac.at/fmuldoo/public_html/webpage/frank-muldoon.html




More information about the Nek5000-users mailing list