Hi all,
Have some tried using outflow boundary conditions in a 2D domain? I am trying to simulate a 2D box with outflow boundary conditions along x. I have used Fourier basis along x and given boundary conditions as
# Background conditions
problem.add_equation("dx(u)(x=0) = 0")
problem.add_equation("dx(u)(x=1) = 0")
In the above my intention is to set du/dx = 0 at left and right of the domain, which I belive gives the outflow boundary conditions? Am I doing this correctly?
While I try to run, an error comes up
/home/manohar/miniforge3/envs/dedalus3/lib/python3.12/site-packages/dedalus/core/arithmetic.py:250: RuntimeWarning: invalid value encountered in add
np.add(arg0.data, arg1.data, out=out.data)
I don't suppose this is because of boundary conditions but I am just trying to make sure.
Also if there is a more straightforward way to implement outflow boundary conditions, could someone tell me?