[Nek5000-users] Uservp and drag_calc routine

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

Are you actually using a variable viscosity ?

I don't know that I've coded up the drag for variable viscosity yet.

Also, I prefer torq_calc() to drag_calc() --- the torq routine is
somewhat newer and drag comes for free.

Paul


On Sat, 14 Aug 2010, nek5000-users at lists.mcs.anl.gov wrote:

> Hello,
>
> I turned on uservp through param(30) in rea file and found that it messes up
> my output from drag_calc(). Shown below are outputs with and w/o uservp [
> Repo : 560 ] .
> Is there any specific flag for calculating integral quantities with uservp ?
>
> With Uservp :
> -----------------
>
> Calculating eddy visosity
>  ediff max :->   1.41399999999999989E-003
>       5  2.2500E-02  3.0084E+00  9.7050E-01  3.0084E+00  9.7050E-01 umx
>        5  2.25000E-02  1.50000E+00  3.46825E-01  6.48743E+00  6.48743E+00
> cdiv
>        5     1377  0.000E+00  6.316E+00  0.000E+00  4.436E-01  6.487E+00
> 1.203E-01divmnmx
>       5  2.2500000E-02  0.0000000E+00  0.0000000E+00  6.5247004E-17
> drgxt*****
>       5  2.2500000E-02  0.0000000E+00  0.0000000E+00  1.0521227E-20
> drgyt*****
>       5  2.2500000E-02  0.0000000E+00  0.0000000E+00 -1.0486441E-20
> drgzt*****
>
> Without Uservp :
> ---------------------
>
>  ediff max :->   1.41399999999999989E-003
>       5  2.2500E-02  3.0107E+00  9.7141E-01  3.0107E+00  9.7141E-01 umx
>        5  2.25000E-02  1.50000E+00  3.46825E-01  6.48743E+00  6.48743E+00
> cdiv
>        5     1377  0.000E+00  6.316E+00  0.000E+00  4.436E-01  6.487E+00
> 1.203E-01divmnmx
>       5  2.2500000E-02  1.3228757E+00  1.0841687E+00  2.3870700E-01  dragx
>   1
>       5  2.2500000E-02 -1.2363700E-04 -1.2361744E-04 -1.9553962E-08  dragy
>   1
>       5  2.2500000E-02  1.6762005E-15  1.5786015E-15  9.7598988E-17  dragz
>   1
> check:  6.28249446E+00  3.43007686E-16  0.00000000E+00        5
>       5  2.2500000E-02  1.3228757E+00  1.0841687E+00  2.3870700E-01  drgxt
>   1
>       5  2.2500000E-02 -1.2363700E-04 -1.2361744E-04 -1.9553962E-08  drgyt
>   1
>       5  2.2500000E-02  1.6762005E-15  1.5786015E-15  9.7598988E-17  drgzt
>   1
>
> Please let me know if rea/usr file is needed. Thanks.
>
> Regards
> Shriram
>