[Nek5000-users] hpts for vorticity field

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

Dear Sandeep,


Please take a look at hemi.usr in the nekexamples repo.


There is a routine interp_v()  that returns the interpolated

velocity field at a user-prescribed set of points.  You could

easily modify this routine to pass in vorticity and interpolate

that (or any other scalar or vector field).


Please note that you should not duplicate interrogation points

on multiple processors (i.e., each processor should have a

unique set of points, xyz();  one possible configuration is that

MPI rank 0 has a list of points of length n and the other ranks

have a list of length 0;  this approach is valid if the length of the

list is not too large).


hth,


Paul


________________________________
From: Nek5000-users <nek5000-users-bounces at lists.mcs.anl.gov> on behalf of nek5000-users at lists.mcs.anl.gov <nek5000-users at lists.mcs.anl.gov>
Sent: Tuesday, December 5, 2017 8:57:20 AM
To: nek5000-users at lists.mcs.anl.gov
Subject: Re: [Nek5000-users] hpts for vorticity field

poor man`s solution: Store your vorticity as an additional scalar field.

-----Original message-----
> From:nek5000-users at lists.mcs.anl.gov <nek5000-users at lists.mcs.anl.gov>
> Sent: Tuesday 5th December 2017 15:55
> To: nek5000-users at lists.mcs.anl.gov
> Subject: [Nek5000-users] hpts for vorticity field
>
> Hello Neks,
> Does anyone knows how to output the vorticity probes using call hpts()
> When I use the function "call hpts()", it outputs only the probes for native arrays like velocity, prssure (and scalar fields, if they exist).
> I compute the vorticity using the function comp_vort3(vort,w1,w2,vx,vy,vz). How can I output the values of the vorticity probes in the flow field using hpts?
> i will be grateful for your answers.
>
> Thank you in advance.
>
> Cheers,
> Sandeep
> _______________________________________________
> 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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20171205/2c8bc7e3/attachment.html>