[Nek5000-users] Uservp and drag_calc routine
nek5000-users at lists.mcs.anl.gov
nek5000-users at lists.mcs.anl.gov
Sat Aug 14 18:51:54 CDT 2010
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
>
More information about the Nek5000-users
mailing list