[Nek5000-users] problem with applying boundary conditions

nek5000-users at lists.mcs.anl.gov nek5000-users at lists.mcs.anl.gov
Thu Sep 14 11:45:59 CDT 2017


How can i leave the mailing list?

⁣Inviato da BlueMail ​

Il giorno 14 set 2017, 18:41, alle ore 18:41, nek5000-users at lists.mcs.anl.gov ha scritto:
>Hi,
>Not sure if I understand your question well, but I guess what happens 
>during the first time step is that through the continuity equation and 
>the incompressibility constraint you will have velocity about one 
>everywhere. This is just to maintain mass conservation.
>
>Philipp
>
>On 2017-09-14 18:29, nek5000-users at lists.mcs.anl.gov wrote:
>> Dear Neks,
>> 
>> I am trying to run a simulation, replicating the Nek5000 run from
>this 
>> youtube video:
>> https://www.youtube.com/watch?v=RmmN-KFSvAY
>> 
>> The inlet BC is constant velocity (which is set to 1.0), and as can
>be 
>> seen from the video, the velocity magnitude reasonably away from the 
>> inlet is zero at starting times. However, in all of my simulations, I
>
>> see the velocity magnitude as one at all the interior points away
>from 
>> the inlet from the beginning (the inlet BC ux=1.0 I am applying is 
>> somehow propagated to all interior elements). Am I doing something
>wrong 
>> somewhere, or is this the expected behavior? Please see my userbc()
>and 
>> useric() from my .usr file below:
>> 
>>        subroutine userbc (ix,iy,iz,iside,eg)
>>        include 'SIZE'
>>        include 'TOTAL'
>>        include 'NEKUSE'
>>        integer ix,iy,iz,eg,e,i,n
>> 
>> CC "My inlet is at x=0"
>> 
>>        if (x.le.1e-4) then
>>          ux=1.0
>>        endif
>> 
>>        return
>>        end
>> 
>>        subroutine useric (ix,iy,iz,eg)
>>        include 'SIZE'
>>        include 'TOTAL'
>>        include 'NEKUSE'
>> 
>>        ux=0
>>        uy=0
>> 
>>        return
>>        end
>> 
>> I believe userbc() is executed for only 'v' in my case. For the other
>
>> constraints on the elements (like 'W', 'O', 'E'),  I think userbc() 
>> would not run (please correct me if I am wrong, I am learning here). 
>> Since the above functions wouldn't work as expected, I changed
>userchk() 
>> and removed the above if statements in userbc(). My userchk() is as
>follows:
>> 
>>        subroutine userchk
>>        include 'SIZE'
>>        include 'TOTAL'
>>        integer i,n
>> 
>>        n = nx1*ny1*nz1*nelv
>> 
>>        if (istep.eq.0) then
>>          do i=1,n
>>            if (xm1(i,1,1,1).le.1e-4) then
>>              vx(i,1,1,1)=1
>>            endif
>>          enddo
>>        endif
>> 
>>        return
>>        end
>> 
>> This didn't work either. Can somebody point me in the right
>direction?
>> 
>> Thanks,
>> Ajay
>> 
>> 
>> 
>> 
>> _______________________________________________
>> Nek5000-users mailing list
>> Nek5000-users at lists.mcs.anl.gov
>> https://lists.mcs.anl.gov/mailman/listinfo/nek5000-users
>> 
>_______________________________________________
>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/20170914/caf12cce/attachment.html>


More information about the Nek5000-users mailing list