19
19
python diagnostics are used to form the sigma term that is assoicated with a p0 dg mesh (ie. element wise) and the interstitial velocity
23
21
to get a darcy equation the time, stress and advection terms are not included in a bousinessq formulation (to avoid needing to specify a density). Ideally it would be better to actually have a darcy momentum option to simplify the input.
25
23
the absorption term is included in the pressure correction (being a necessity as it is the only term in the matrix equation)
29
27
NOTE: Because the velocity and pressure are time independent the simulation is broken into two.
30
28
First the velocity and pressure are solved and checkpointed.