<div dir="ltr"><div><div>Dear Nek users, I've been struggling to program an adecuate outlet BC. I simulate a vertical isothermal surface in an open domain, forming a free convection boundary layer. When there are vortex leaving the domain, the code diverges. I've tried turb_outflow subroutine, stabilized outflow BC (Dong 2014), convective BC (sommerfield radiation condition, thanks to Phillip and Eric, from this comunity), etc. But I can't achieve to make the code run propperly, and I can't go on with my research while I have this problem.</div><div><br></div><div>There is a paper, from Desrayaud (2013), "Benchmark solutions for natural convection flows in vertical channels submitted to different open boundary conditions". </div><div><br></div><div>I want to try these BC, which can be stated as: (x: horizontal and y:vertical.)</div><div><br></div><div>At the top face (outlet)</div><div><br></div><div>Ux=0</div><div>dUy/dy=0</div><div><br></div><div>If Uy > 0 (fluid going out)</div><div>dTEMP/dy=0</div><div>pressure=0</div><div><br></div><div>If Uy < 0 (fluid going in because there is backflow)</div><div>TEMP=0</div><div>pressure= -0.5Uy^2 or 0   (these are two options in the paper)</div><div><br></div><div><br></div><div><br></div><div>As far as I know, in Pn-Pn-2 formulation there are not pressure BC, but it is possible to set it using Pn-Pn. Am I right? If so, how can I set BC for pressure in Pn-Pn?</div><div><br></div><div>And for the velocity gradient and temperature gradient equal cero at the top face, how can it be done? Would it be enough to set dirichlet BC and replicate the last row of values to the last one? </div><div><br></div><div>Any help is welcomed, thanks in advance,</div><div><br></div><div>Juan Pablo.</div></div></div>