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

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

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