[Nek5000-users] Torq_calc() position of X0

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

Hi Paul, 

Thanks for the reply.

Yes I am looking for the axial viscous force in the whole domain which I guess is nothing but the viscous drag along the z direction. 

Thanks,

Kamal 


On Feb 3, 2014, at 10:17 PM, nek5000-users at lists.mcs.anl.gov wrote:

> 
> Hi Kamal,
> 
> You definitely want to give an array for x0, else
> you don't know what you'll be referencing in memory.
> 
> You can turn off the torque output and just get the
> drag, which I think is what you want.
> 
> Paul
> 
> 
> On Mon, 3 Feb 2014, nek5000-users at lists.mcs.anl.gov wrote:
> 
>> Hello Neks,
>> 
>> I worked around the set_obj() and torq_calc() to figure it out.
>> 
>> I set my sey_obj() routine to find the wall.
>> 
>> If I am not wrong,  the output of the torque_calc() is Total torque, viscous force torque and pressure force torque taken about a point X0
>> 
>> If I don't define X0, what point does nek5000 automatically take ?
>> 
>> Thank you,
>> Kamal
>> _______________________________________________
>> Nek5000-users mailing list
>> Nek5000-users at lists.mcs.anl.gov
>> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>> 
> _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users