[Nek5000-users] New Recycling BC based on interpolation
nek5000-users at lists.mcs.anl.gov
nek5000-users at lists.mcs.anl.gov
Fri Jul 23 11:26:00 CDT 2010
Hi Aleks and Stephan,
I have a question directly related to this. I would like to use the hpts() routine for my conjugate case, however like Markus mentioned, there is no support for a conjugate case at the moment. Were ya'll able to fix this?
Thanks,
Michael
----- Original Message -----
From: nek5000-users at lists.mcs.anl.gov
To: nek5000-users at lists.mcs.anl.gov
Sent: Monday, July 19, 2010 8:38:05 AM GMT -06:00 US/Canada Central
Subject: Re: [Nek5000-users] New Recycling BC based on interpolation
Hi Aleks,
what was your plan to fix it?
In fact I think findpts_eval() should work as is but there is no-err handling if you try to interpolate on a point in the solid using an input field based on the v-mesh.
Ok we don't need an err-handling at the moment. You'll get a ZERO value back although it is an invalid point for an input field based on the v-mesh.
Why is that?
At the moment, all field arrays have the same layout and size (lx1,ly1,lz1,lelt) independent of the mesh type (v or T). We simply set all values of a v-variable to ZERO for iel>nelv and we don't touch these segments later on.
Stefan
On Jul 19, 2010, at 3:25 PM, nek5000-users at lists.mcs.anl.gov wrote:
> Hi Stefan,
>
> Yes, we have noticed you changed the interface -- and the latest version of Recycling Bcs w/ interpolation works for the latest repo version (548).
>
> Eventually we can get it working for CHT too.
>
> Best,
> Aleks
>
>
>
> ----- Original Message -----
> From: nek5000-users at lists.mcs.anl.gov
> To: nek5000-users at lists.mcs.anl.gov
> Sent: Monday, July 19, 2010 7:05:21 AM GMT -06:00 US/Canada Central
> Subject: Re: [Nek5000-users] New Recycling BC based on interpolation
>
> Hi Guys,
>
> I worked on the interpolation wrappers over the last months and I am not sure if the new recycling BC is compatible to the latest repo-version?
>
> Paul/Aleks: can you comment on that?
>
>
> Now, in the case of CHT (nelgt>nelgv) findpts_eval() won't work! I guess the size of the input field has to match the dimensions of XM1,YM1,ZM1 specified in inpts_setup().
>
> I will talk to James if we can come up with an easy fix for that.
>
>
> Stefan
>
>
>
> -----Original Message-----
> From: nek5000-users-bounces at lists.mcs.anl.gov on behalf of nek5000-users at lists.mcs.anl.gov
> Sent: Sun 7/18/2010 13:53
> To: nek5000-users at lists.mcs.anl.gov
> Subject: Re: [Nek5000-users] New Recycling BC based on interpolation
>
>
>
> Hi Markus,
>
> this stuff is all brand new and you are in the leading edge
> of the testing (since none of us are exercising it in our
> own applications), so there will be a bit of development as
> we go forward. We can readily fix interp_v (but I advise
> using it now in order to identify other issues).
>
> interpv is _much_ more expensive than gs_op.
>
> However, if one takes 1 part in a million and the other 1
> part in 10,000 of the run time, then both are negligible.
> So the question is really relative to the wall-clock time.
>
> Paul
>
>
>
> On Sun, 18 Jul 2010, nek5000-users at lists.mcs.anl.gov wrote:
>
>> Hi,
>>
>> no, I haven't done any speed comparison. So I take it there is no principal
>> difference between intpts and gs_op in terms of computational cost?
>>
>> I was referring to line 172f of the new recycling BC code:
>> "
>> if (nelgt.ne.nelgv) call exitti
>> $ ('ABORT: interp_v() nelgt.ne.nelgv not yet supported!$',nelgv)
>> "
>> and was wondering if that is due to something in intpts and if so, if the
>> same constraint would apply to findpts.
>>
>> Markus
>>
>>
>> nek5000-users at lists.mcs.anl.gov wrote:
>>>
>>> Hi Markus,
>>>
>>> I don't know why this wouldn't work for nelgt > nelgv.
>>>
>>> I would say we don't know the performance hit yet - I think
>>> someone is looking into this question. Performance will of course vary
>>> significantly from problem to problem and platform
>>> to platform. I suspect that it should not be a major hit ---
>>> have you observed otherwise?
>>>
>>> Paul
>>>
>>>
>>> On Sat, 17 Jul 2010, nek5000-users at lists.mcs.anl.gov wrote:
>>>
>>>> Hi Aleks,
>>>>
>>>> thanks for the code, I looked over it and have two questions:
>>>> -at some point, it is mentioned that this method does not work if
>>>> nelgt>nelgv. Why is that? I am asking because this will be applied to a
>>>> conjugate HT problem (although the recycling part is purely concerning
>>>> flow+convection). Will this problem also occur when only findpts is used
>>>> instead of intpts?
>>>> -is there a performance difference between calling intpts as opposed to
>>>> gs_op (as is done in the turbJet example) every time step for large scale
>>>> production runs? The latter is what I used in my version.
>>>>
>>>> Thanks,
>>>> Markus
>>>>
>>>>
>>>> nek5000-users at lists.mcs.anl.gov wrote:
>>>>> Hi Michael and Markus,
>>>>>
>>>>> Here is a test case for an unsteady pipe flow with Michael's
>>>>> GambitPipe.rea file which sets inlet Recycling BCs based on interpolation
>>>>> from a cross-section downstream that you specify by a multiple of the
>>>>> inlet normal:
>>>>>
>>>>> http://www.mcs.anl.gov/~obabko/z4.tgz
>>>>>
>>>>>
>>>>> Let me know if you have any questions.
>>>>> Best,
>>>>> Aleks
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> ----- Original Message -----
>>>>> From: nek5000-users at lists.mcs.anl.gov
>>>>> To: nek5000-users at lists.mcs.anl.gov
>>>>> Sent: Thursday, July 8, 2010 5:06:48 PM GMT -06:00 US/Canada Central
>>>>> Subject: Re: [Nek5000-users] New Recycling BC method?
>>>>>
>>>>>
>>>>>
>>>>>
>>>>>
>>>>> No worries, just checkin in to see how it was going, thanks for the
>>>>> update! Look forward to trying it out!
>>>>>
>>>>> - Michael _______________________________________________
>>>>> 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
>>>>
>>> _______________________________________________
>>> 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
>>
> _______________________________________________
> 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
> _______________________________________________
> 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/20100723/58885daf/attachment.html>
More information about the Nek5000-users
mailing list