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

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

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/645db45c/attachment-0001.html>