[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