Ferrite.jl
Ferrite.jl copied to clipboard
How-to: Projection methods for time-dependent Navier-Stokes equation
In the #ferrite-fem channel of Julia Slack, in asking questions about the incompressible Navier-Stokes tutorial, @termi-official stated "we shall also note that the used Navier Stokes time discretization in our examples is not the usual kind of solver used in CFD for this kind of problem." Maybe a how-to guide that builds on the tutorial but instead uses projection methods would be beneficial? Some implementation references/learning resources are provided:
- deal.ii: The step-35 tutorial program
- FEniCSx Tutorial: The Navier-Stokes Equations
- YouTube Video -- FEniCS Tutorial: Navier-Stokes Equation for Lid-Driven Cavity
- Linear Elasticity Weak Form Derivation. Relevant because a theory section in an NS tutorial/how-to could show in an Extra Details sub-section the step-by-step integration by parts one applies. Books and other tutorials seem to gloss over this step for vector fields (i.e., the integration by parts is often shown only for the laplace of a scalar fields e.g., through the heat equation but not vector fields), so I think it'd be nice to show for pedagogical purposes. This step is also shown in Weak-Form for Navier-Stokes with Chorin's projection.
This would be something that I could look into implementing over the next month(s) if it's something that you think would be useful.
I would like to be assigned this issue and will start working on it ASAP. I believe it will be a good first step to get more familiar with Ferrite internals before attempting to contribute to other issues related to functionality (e.g., VEM implementation, etc.).