[Nek5000-users] variable properties, uservp
nek5000-users at lists.mcs.anl.gov
nek5000-users at lists.mcs.anl.gov
Tue May 4 06:01:53 CDT 2010
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
More information about the Nek5000-users
mailing list