[Nek5000-users] Interpolation between two meshes

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

Hi Stefan,


Many thanks for your reply. Please contact me via wzx1989221 at gmail.com.


Best regards,

Tony


________________________________
From: nek5000-users-bounces at lists.mcs.anl.gov <nek5000-users-bounces at lists.mcs.anl.gov> on behalf of nek5000-users-request at lists.mcs.anl.gov <nek5000-users-request at lists.mcs.anl.gov>
Sent: 30 March 2017 13:17
To: nek5000-users at lists.mcs.anl.gov
Subject: Nek5000-users Digest, Vol 97, Issue 37


Today's Topics:

   1. Outflow BC with negative velocity
      (nek5000-users at lists.mcs.anl.gov)
   2. Re: NekNek with PN-PN (nek5000-users at lists.mcs.anl.gov)
   3. Re: NekNek with PN-PN (nek5000-users at lists.mcs.anl.gov)
   4. Re: Interpolation between two meshes
      (nek5000-users at lists.mcs.anl.gov)
   5. Accessing the value of drag and lift coefficients in usrchk
      (nek5000-users at lists.mcs.anl.gov)



Message: 4
Date: Thu, 30 Mar 2017 11:46:40 +0200
From: nek5000-users at lists.mcs.anl.gov
To: nek5000-users at lists.mcs.anl.gov      <nek5000-users at lists.mcs.anl.gov>
Subject: Re: [Nek5000-users] Interpolation between two meshes
Message-ID:
        <mailman.4068.1490867217.2967.nek5000-users at lists.mcs.anl.gov>
Content-Type: text/plain; charset="utf-8"

Please let me know how I can contact you off list.

Stefan?

On 29 Mar 2017, at 14:59, "nek5000-users at lists.mcs.anl.gov <mailto:nek5000-users at lists.mcs.anl.gov> " <nek5000-users at lists.mcs.anl.gov <mailto:nek5000-users at lists.mcs.anl.gov> > wrote:

Hi Stefan,

Thank you very much for your quick reply.

1. Yes, my source field is in double precision since I need to keep the accuracy.

2. I added several print lines in gfldr.f to check where the interpolation hangs (see below). I also put 'call nekgsync()' before and after findpts as you suggested. I can confirm that findpts is exactly where the interpolation hangs. The last output in the logfile is 'CHECK interp findpts before' and it never gets updated again.

3. I've just tested gfldr() with turbChannel example. Everything worked fine and the?source field was interpolated onto the turbChannel mesh (512 elements). So I moved on increasing the number of elements of turbChannel case.

First, I used 234300 elements?in the turbChannel example (half of my current simulation size) and the interpolation was done successfully. Then, I used 468600 elements for?the turbChannel example (similar to my current simulation size), the interpolation hangs again at findpts.

I guess it might be something to do with the memory. What do you think? Is there any solution to this?

Best regards,
Tony

??????? if(nid==0) print*,'CHECK interp findpts before'
??????? call nekgsync()
??????? call findpts(inth_gfldr,
???? &?????????????? rcode,1,
???? &?????????????? proc,1,
???? &?????????????? elid,1,
???? &?????????????? rst,ndim,
???? &?????????????? dist,1,
???? &?????????????? xm1,1,
???? &?????????????? ym1,1,
???? &?????????????? zm1,1,ntot)
??????? call nekgsync()
??????? if(nid==0) print*,'CHECK interp findpts after'
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20170330/af4210eb/attachment-0001.html>