[Nek5000-users] Change boundary condition type with time?

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

      include 'TSTEP' ! time, ifield
      include 'INPUT' ! cbc


      if (ifield .eq. 1 .and. time .gt. T) then
         nface = 2*ndim


         do iel=1,nelv


            do ifc=1,nface
               if ("face belongs to part B") then   ! You need to identify
your boundary B here
                  cbc (ifc,iel,ifld) = 'O  ' ! or 'o  ' or 'on ' or 'ON '
depending on what you need there
               endif


            enddo


         enddo
      endif

On Wed, Apr 11, 2018 at 2:08 PM,  <nek5000-users at lists.mcs.anl.gov> wrote:
> Hi Neks,
>
>
> I wonder if it is possible to have boundary conditions that chage type
with
> time?
>
>
> I want to simulate a flow where I make a 'puff' of fluid through a hole
in a
> wall.
>
> After the puff is made it continues to evolve without beeing forced.
>
>
> Mathematically, this is what I want:
>
> Dirichlet boundary conditions (ux,uy,uz) = (0,0,1) for the velocity at a
> part B of the boundary up to a time T.
> Open boundary condition for the velocity at B for time > T.
>
>
> Is this possible to implement?
>
>
> Best,
>
>
> Johan
>
>
>
> _______________________________________________
> Nek5000-users mailing list
> Nek5000-users at lists.mcs.anl.gov
> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20180412/be7cb8dc/attachment.html>