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
https://ngsolve.org/media/kunena/attachments/1359/H(div)-DG-NS.py
Attachment: Hdiv-DG-NS_2020-02-22.py