[Nek5000-users] variable properties, uservp

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

uservp() is called if:

IFUSERVP = true
or
param(30) > 0
or
if you have a non-trivial 'Variable Property Data' section in your .rea file with MATYPE=2


Stefan


On May 4, 2010, at 12:53 AM, nek5000-users at lists.mcs.anl.gov wrote:

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