Add Stokes/Poiseuille 1D FD playground: temporal convergence for index-1 semi-explicit DAE with constraint lifting, differentiated-constraint remedy, coupled block-solve exploration, and monolithic RADAU-IIA solver - #641

Draft

pancetta with Copilot wants to merge 9 commits into

copilot/test-temporal-order-convergence-again

Commits on Mar 26, 2026

Redesign: use genericImplicitConstrained sweeper with MeshDAE state

Copilot
and
pancetta
authored
a92be2d
Implement Stokes/Poiseuille DAE playground with SemiImplicitDAE sweeper

Copilot
and
pancetta
authored
18f3a88
Add constraint-lifting variant stokes_poiseuille_1d_fd_lift; update run_convergence.py to compare both

Copilot
and
pancetta
authored
d6aa801
Use nu=0.1, Tend=1.0 for clean M+1=4 convergence orders; document U-formulation order limit

Copilot
and
pancetta
authored
289e581
Add FullyImplicitDAE-compatible Stokes classes; update run_convergence.py with 4-variant comparison

Copilot
and
pancetta
authored
6d21c46
Add differentiated-constraint classes achieving full 2M-1 RADAU order for velocity and pressure

Copilot
and
pancetta
authored
870f342
Fix diffconstr order claim: M+2 (not 2M-1); for M=3 M+2=5=2M-1 is coincidental

Copilot
and
pancetta
authored
6e79a2c
Add stokes_poiseuille_1d_fd_coupled: (N+1)x(N+1) block solve with optional projection

Copilot
and
pancetta
authored
0aad670
Add run_monolithic.py: monolithic RADAU-IIA collocation (algebraic and differentiated-constraint variants)

Copilot
and
pancetta
authored
2016d90