what to do with the upwind term for Navier-stokes

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())

Hi! Guosheng. Since I just started to contact this software, I have a few problems. Please help me to have a look . Thank you!

Wang Shasha


Attachment: Hdiv-DG-NS_2020-02-22.py