How to implement Neumann BC on non-Poisson PDE

2

Take a look at this thread, for example. There’s no concept of a Neumann boundary condition for an advection problem. You need to consider fluxes of your solution on the boundary.