[Nek5000-users] Uservp and drag_calc routine

nek5000-users at lists.mcs.anl.gov nek5000-users at lists.mcs.anl.gov
Sun Aug 15 20:43:55 CDT 2010


Shriram,

I'm not having this difficulty either with T IFUSERVP
nor with p30=1 in the .rea file.

If you'd like to tar up .rea/.usr/SIZE and send them to
me off-list I can take a look.

Paul


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

>
> Interesting... I'll try this.
>
> Paul
>
>
> On Sat, 14 Aug 2010, nek5000-users at lists.mcs.anl.gov wrote:
>
>> Paul,
>> 
>> Thanks for the reply.
>> 
>> Yes. I set my wall as an obj through the same routine. It prints the value
>> if param(30) is 0 but doesn't if param(30) is 1.
>> 
>> I looked at the drag_calc() routine and found that it works for constant
>> viscosity since param(2) is initialized to all the points while  torq_calc
>> allows variable viscosity through visc array and so thought it might work
>> with param(30)=1.
>> 
>> Regards
>> Shriram
>> 
> _______________________________________________
> 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