[Nek5000-users] Interpolation on a finer mesh

nek5000-users at lists.mcs.anl.gov nek5000-users at lists.mcs.anl.gov
Tue Jun 20 08:39:10 CDT 2017


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
> 
> 
> -- 
> Karlsruhe Institute of Technology (KIT)
> Institute of Fluid Mechanics
> 
> M.Sc. Steffen Straub
> Doctoral Researcher
> 
> Kaiserstraße 10
> Building 10.23
> 76131 Karlsruhe, Germany
> 
> Phone: +49 721 608-43027
> E-mail: steffen.straub at kit.edu
> Web: http://www.istm.kit.edu
> 
> KIT – The Research University in the Helmholtz Association
> 
> Since 2010, the KIT has been certified as a family-friendly university.
> 
> _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
> 


More information about the Nek5000-users mailing list