[Nek5000-users] Interpolation on a finer mesh

nek5000-users at lists.mcs.anl.gov nek5000-users at lists.mcs.anl.gov
Wed Jun 21 04:12:19 CDT 2017


Hello Stefan,

thanks for your quick reply.
I have just tested the generic fld reader gfldr as you suggested.

It is super easy and works like a charm.


Cheers,
Steffen


On 06/20/2017 07:00 PM, nek5000-users-request at lists.mcs.anl.gov wrote:
> Message: 2 Date: Tue, 20 Jun 2017 15:39:10 +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 on a finer mesh Message-ID: 
> <mailman.9088.1497966043.2967.nek5000-users at lists.mcs.anl.gov> 
> Content-Type: text/plain; charset=utf-8 You can use our generic fld 
> reader. Just add this to userchk(): if(istep.eq.0) call 
> gfldr('myfldfilename') This will interpolate your field file data to 
> your new mesh. You're all set :) Cheers, Stefan -----Original 
> message-----
>> From:nek5000-users at lists.mcs.anl.gov  <nek5000-users at lists.mcs.anl.gov>
>> Sent: Tuesday 20th June 2017 15:19
>> To:nek5000-users at lists.mcs.anl.gov
>> Subject: [Nek5000-users] Interpolation on a finer mesh
>>
>> Dear Neks,
>>
>> I am simulating turbulent pipe flow with Nek5000. As a first step, I
>> have generated a fully developed turbulent pipe flow within a domain of
>> L_z = 5 D at Re_\tau = 360. Now, I would like to use this flow field as
>> an initial / restart condition for a turbulent flow at a Re_\tau=720
>> with a finer mesh.
>> The geometry of the finer mesh is contained in a .rea file.
>>
>>
>> For doing the interpolation, I have found the "interpolation wrapper for
>> usage in .usr file" in core/intp_usr.f.
>> Can someone give a short example on how to use these subroutines to
>> interpolate results from 3 or 4 restart files onto my new geometry?
>>
>>
>> Kind Regards,
>> Steffen Straub

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20170621/9771ff3b/attachment.html>


More information about the Nek5000-users mailing list