[Nek5000-users] Open Boundary condition in NEK user Documentation

nek5000-users at lists.mcs.anl.gov nek5000-users at lists.mcs.anl.gov
Thu Oct 6 08:12:15 CDT 2016


Thank you,


this makes sense to me.


/Johan


________________________________
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, October 6, 2016 2:22 PM
To: nek5000-users at lists.mcs.anl.gov
Subject: Re: [Nek5000-users] Open Boundary condition in NEK user Documentation


Dear Johan,

Actually - I was mistaken.

What we use is:

     d/dn (u_i) = 0

for i=1,...,3

We could however (and sometimes do) require Neumann only on the normal velocity component --

grad (phi) . nhat = 0

where phi = U . nhat

The tangential components can be prescribed.   That is the capability provided by the "on "
BC.

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, October 06, 2016 6:58 AM
To: nek5000-users at lists.mcs.anl.gov
Subject: Re: [Nek5000-users] Open Boundary condition in NEK user Documentation


Thank you for a good answer!


Just to clarify, is it in LaTeX written as this:


(\hat{n} \cdot \nabla)(u \cdot \hat{n}) = 0 ?


If the wall normal is in the z-direction, then the open boundary condition for the velocity is:


d/dz u_z = 0?


Best regards,


Johan




________________________________
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, October 6, 2016 12:40 PM
To: nek5000-users at lists.mcs.anl.gov
Subject: Re: [Nek5000-users] Open Boundary condition in NEK user Documentation


Dear Johan,

The outflow boundary condition, 'O  ', is indeed the same as 'I  ' for the thermal problem
since they both enforce grad T . nhat = 0.

However, both also admit thermal flux by advection for cases where U . nhat is > 0, and that is
the intent of outflow --- to allow thermal flux to be carried out.

For velocity, the situation is a bit more delicate.  To leading order, however, 'O   ' imposes
d/dn (U.nhat) = 0 and p=0.    That is, Neumann conditions for the viscous fluxes and Dirichlet
for pressure.

hth,

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, October 06, 2016 4:59 AM
To: nek5000-users at lists.mcs.anl.gov
Subject: [Nek5000-users] Open Boundary condition in NEK user Documentation


Dear Neks,


I am trying to understand what is actually meant by an "open boundary" in Nek.

Un the user documentation, eqn (4.5) sais that, for the temperature, it is


grad(T) dot n = 0, where n is the normal to the boundary.


Since the heat flux is q = - grad(T), I would rather interpret this as a zero flux boundary condition.


And indeed, the insulated boundary described by equation (4.6) is identical to the open boundary condition.


Is equation (4.5) a type-O or is the open boundary for the temperature actually the same as an insulated boundary?


Also: Should not the open boundary condition for the "no-stress forumlation" also be given by the same expression as for the "stress-formulation" eqn (4.2)?


Best Regards,


Johan
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.mcs.anl.gov/pipermail/nek5000-users/attachments/20161006/393177f0/attachment-0001.html>


More information about the Nek5000-users mailing list