We have moved the forum to https://forum.ngsolve.org. This is an archived version of the topics until 05/05/23. All the topics were moved to the new forum and conversations can be continued there. This forum is just kept as legacy to not invalidate old links. If you want to continue a conversation just look for the topic in the new forum.
Hello everyone! I'm debugging the H (div)-DG code for steady-steady Navier-Stokes now, but I don't know what to do with the upwind term. Which command code should I use?.
You need two things in your DG formulation:
(1) fes = FESpace([V,Q], dgjumps=True)
(2) For upwinding, simply use the IfPos function:
IfPos(u*n, u, u.Other())