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!