[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