-
-
Notifications
You must be signed in to change notification settings - Fork 182
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Add divergence conforming discontinuous Galerkin demo for the Navier-Stokes equations #2390
Conversation
If anyone feels this would be better kept outside the main repository, it could be added to the FEniCSx tutorial instead if @jorgensd is happy |
Co-authored-by: Nate <[email protected]>
python/demo/demo_navier_stokes.py
Outdated
# $\Omega \subset \mathbb{R}^d$, $d \in \{2, 3\}$, and time interval | ||
# $(0, \infty)$, given by | ||
# $$ | ||
# \partial_t u - \nu \Delta u + (u \cdot \nabla)u + \nabla p = f |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Does \nu
actually get used in the discretisation?
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
The discretisation uses the dimensionless form, where \nu
is incorporated into the Reynolds number. I'll make this clearer
Co-authored-by: Nate <[email protected]>
…inx into jpdean/div_con_navier_stokes
UFL doesn't properly support DG-type operators with complex.
…inx into jpdean/div_con_navier_stokes
Need to figure out how to include extra latex packages
This PR adds a divergence conforming discontinuous Galerkin demo for the Navier-Stokes equations.
It demonstrates:
As far as I'm aware, none of the other demos demonstrate the above functionality.
I still need to so some tidying. Any feedback would be much appreciated.