[Nek5000-users] Pn-Pn formulation problem with accl_outflow()

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

Hi Paul,
Thanks for your reply.
I will look into it.

Mohsen

On 1-May-14, at 8:51 PM, <nek5000-users at lists.mcs.anl.gov> wrote:

> Hi Mohsen,
>
> I suggest the turb_outflow() subroutine call.   The routine is in  
> navier5.f ---
>
> You need to declare two arrarys in userchk and then call  
> turb_outflow on each step.
>
>  subroutine turb_outflow(d,m1,rq,uin)
>
> c     . Set div U > 0 in elements with 'O  ' bc.
> c
> c     . rq is nominally the ratio of Qout/Qin and is typically 1.5
> c
> c     . d and m1 are work arrays of size (lx1,ly1,lz1,lelt),  
> assumed persistant
>
> d and m1 should be declared "real"  ... A typical call would be:
>
> c     real m1(lx1,ly1,lz1,lelt),d(lx1,ly1,lz1,lelt)
> c
> c     uin = 1.    ! Default characteristic velocity scale
> c     rq  = 1.5   ! Outflow/Inflow ratio
> c
> c     call turb_outflow(d,m1,rq,uin)
>
> ... (without the "c" in the 1st column)
>
> Paul
>
>
>
> ________________________________________
> From: nek5000-users-bounces at lists.mcs.anl.gov [nek5000-users- 
> bounces at lists.mcs.anl.gov] on behalf of nek5000- 
> users at lists.mcs.anl.gov [nek5000-users at lists.mcs.anl.gov]
> Sent: Thursday, May 01, 2014 1:40 PM
> To: Nek 5000
> Subject: [Nek5000-users] Pn-Pn formulation problem with accl_outflow()
>
> Hello Nek,
> I tried to use "accl_outflow()" subroutine for a turbulent outflow in
> my simulation. I noticed that there is an issue in this subroutine
> when using Pn-Pn formulation. In fact, the user defined divergence is
> added to all nodes beyond those of the outflow elements. Similar
> problem was noticed before as I searched the mailing list (see the
> following link):
>
> https://lists.mcs.anl.gov/mailman/htdig/nek5000-users/2010-July/
> 000862.html
>
> Pn-Pn-2 formulation is not supported in my case in which I use Moab.
> I wondered if anyone came up with a solution to this issue.
>
> Thanks,
> Mohsen
> _______________________________________________
> 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