[Nek5000-users] Average in arbitrary direction (fwd)

Thu May 19 16:15:14 CDT 2011

Hi Josh,

The routine that you want is attached to this email or you can find it here:
http://www.mcs.anl.gov/~jwalker/interp_v.f

You give it a packed list of the points that you want to interpolate

x(1),y(1),z(1),....,x(n),y(n),z(n)

along with the number of points, n, and it gives you back the velocities at
each point

u(1),v(1),w(1),...,u(n),v(n),w(n).

Make sure to change lpart in SIZE to the number of points you want.

Best,
Justin

On Thu, May 19, 2011 at 3:09 PM, Paul Fischer <fischer at mcs.anl.gov> wrote:

>
> Justin,
>
> Can you post or send the latest interp_uvw() routine that
> you used in your hemi case?
>
> I guess we need to get the hemi example into the repo.
>
> Thanks,
>
> Paul
>
>
> ---------- Forwarded message ----------
> Date: Thu, 19 May 2011 15:06:05 -0500 (CDT)
> From: nek5000-users at lists.mcs.anl.gov
> To: nek5000-users at lists.mcs.anl.gov
> Subject: Re: [Nek5000-users] Average in arbitrary direction
>
> Hi Josh,
>
> yes, there are 3 intpts* routines from postpro.f that you can use. I have
> an example but I need to update it for the current repo version of this
> routines -- I'll send it to you once I modify and test it.
>
> Best,
> Aleks
>
>
> ----- Original Message -----
> From: nek5000-users at lists.mcs.anl.gov
> To: nek5000-users at lists.mcs.anl.gov
> Sent: Thursday, May 19, 2011 9:21:03 AM
> Subject: Re: [Nek5000-users] Average in arbitrary direction
>
> Hello Paul,
>
> This sound like exactly what I need.  I didn't find anything like
> interp() while using grep in the repo--are you referring to the
> routines found in postpro.f ??
>
> Thanks,
>
> Josh
>
> On Wed, May 18, 2011 at 3:59 PM,  <nek5000-users at lists.mcs.anl.gov> wrote:
>
>>
>> Hi Josh,
>>
>> You can use the interp() routines --- let me check if we have
>> examples in the repo.
>>
>> These are somewhat expensive, so if you're also integrating in
>> time (even over just 10 steps), it would pay to only interpolate
>> the time-averaged data, rather than interpolating every timestep
>> and then averaging the profiles.   However, I can understand that
>> there are times you'd want the full time signal.
>>
>> The best integration formulae for periodic signals (e.g., in your theta
>> direction and perhaps in z, if you have periodicity)
>> is to simply sum the values, sampled uniformly in space (z or
>> theta), and dived by n, the number of samples.   You do need
>> to sample with sufficient density to capture all the wiggles ---
>> another reason that time averaging has value.
>>
>> Paul
>>
>>
>>
>> On Wed, 18 May 2011, nek5000-users at lists.mcs.anl.gov wrote:
>>
>>  Hello Neks,
>>>
>>> I have looked at the turbChannel example for an example on planar
>>> averaging.
>>>
>>> However, I was wondering if there was any advice  on how to do this if
>>> the elements aren't exactly "nice"... (i.e., the elements may have
>>> arbitrary rotations)
>>>
>>> Specifically, I have a pipe flow, and I am wanting to do an average in
>>> the z (axial) direction.  My initial thoughts were to interpolate my
>>> values to set R, theta components, and then average this along z.  I'm
>>> stuck on the best way to do either one of these operations, however
>>> (especially because I would like to use an integration rule to compute
>>> the average in the z direction).
>>>
>>> If anyone could point me in the right direction, that would be a big
>>> help.  Thanks!
>>>
>>> --
>>> Josh Camp
>>>
>>> "All that is necessary for the triumph of evil is that good men do
>>> nothing" -- Edmund Burke
>>> _______________________________________________
>>> 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
>>
>>
>
>
> --
> Josh Camp
>
> "All that is necessary for the triumph of evil is that good men do
> nothing" -- Edmund Burke
> _______________________________________________
> 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/20110519/8a0781fb/attachment.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: interp_v.f
Type: text/x-fortran
Size: 1446 bytes
Desc: not available
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20110519/8a0781fb/attachment.bin>