[Nek5000-users] Impose inflow conditions from another simulation

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

Hello,

As a first step, here's a code snippet that would give you the xyz
coordinates of the inlet (assuming there is only one inlet):

do iel=1,nelv
  do ifc=1,2*ndim
    if (cbc(ifc,iel,1).eq.'v  ') then
      call facind(kx1,kx2,ky1,ky2,kz1,kz2,nx1,ny1,nz1,ifc)
       do iz=kz1,kz2
       do iy=ky1,ky2
       do ix=kx1,kx2
         x = xm1(ix,iy,iz,iel)
          y = ym1(ix,iy,iz,iel)
          z = zm1(ix,iy,iz,iel)
         !  Do whatever printing you want to do here
        enddo
       enddo
        enddo
    endif
   enddo
enddo

I don't think nek has any built-in routines for interpolating face data,
but one could certainly build one using some of the functions from
speclib.f.

Josh


On Wed, Feb 27, 2013 at 5:52 AM, <nek5000-users at lists.mcs.anl.gov> wrote:

> Hello nek5000-users,****
>
> ** **
>
> I have a data file obtained through a RANS simulation that I would like to
> use as the inflow condition for my subsequent nek5000 simulation. ****
>
> ** **
>
> May I know whether there is an efficient subroutine which allows me to
> interpolate the RANS data onto the spectral-element mesh or, at least,
> allows me to dump out the nodal points at the inlet so that I can perform
> the interpolation and assignment? ****
>
> ** **
>
> I will be performing a parallel simulation involving 300000 elements.
> Hence I will be grateful if someone could provide me with an efficient
> method.****
>
> ** **
>
> Thanks in advance.****
>
> ** **
>
> *Ng Jee Hann (Mr.)*****
>
> *PhD Student.*
>
> Department of Mechanical Engineering,****
>
> National University of Singapore,****
>
> 9 Engineering Drive 1, ****
>
> Singapore 117576.****
>
> ** **
>
> _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>
>


-- 
Josh Camp

"All that is necessary for the triumph of evil is that good men do nothing"
-- Edmund Burke
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20130308/b48c4131/attachment.html>