jorgensd / dolfinx-tutorial

A reimplementation of the Springer book: https://github.com/hplgit/fenics-tutorial/, covering new topics as well as transitioning from dolfin to dolfinx
https://jorgensd.github.io/dolfinx-tutorial/
110 stars 63 forks source link

Issue on page /chapter2/ns_code2.html #193

Closed mcastel1 closed 5 months ago

mcastel1 commented 5 months ago

I don't see why the equation Screenshot 2024-05-08 at 18 41 12 containts $p^{n-1/2}$.

In the analogous equation for the tentative velocity $u^\ast$ here, see Screenshot 2024-05-08 at 18 42 51, $p^\ast$ enters the equation, not $p^{n-1/2}$.

jorgensd commented 5 months ago

$p^*$ in the Oasisx tutorial is $p^{n-3/2}$, i.e. the pressure from the last time step. In oasisx I used the notation $u^n$, $p^{n-1/2}$ to denote the solutions at a given time step, while in the tutorial I use a slightly different indexing (add 1 to n everywhere).