Pith. sign in

REVIEW 3 major objections 4 minor 50 references

Stabilizing only the predicted velocity yields a second-order-accurate, cheaper Navier-Stokes projection scheme.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

T0 review · deepseek-v4-flash

2026-08-05 00:11 UTC pith:PBKC6ULH

load-bearing objection A genuinely new and useful VMS-projection combination with strong benchmarks, carrying a formal error estimate that is more assumed than proved. the 3 major comments →

arxiv 2608.00843 v1 pith:PBKC6ULH submitted 2026-08-01 math.NA cs.NA

A Helmholtz-Leray projection method with variational multiscale stabilization for the Navier-Stokes equations

Biswajit Khara , Suresh Murugaiyan , Makrand Khanwale , Baskar Ganapathysubramanian This is my paper
classification math.NA cs.NA MSC 65M6065M1576D0576M10
keywords projection methodsvariational multiscale stabilizationNavier-Stokes equationspressure-correction schemesfinite element methodserror analysisTaylor-Green vortexequal-order elements
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper tries to establish that the expensive fully coupled form of variational multiscale (VMS) stabilization for incompressible flow can be replaced by a fractional-step projection scheme in which only the intermediate velocity is decomposed into coarse and fine scales. Because the projection step removes the velocity-pressure saddle point, the pressure and the final divergence-free velocity need no fine-scale stabilization; the missing advective stability is supplied by a residual-based model of the predicted velocity's unresolved part, which also feeds a PSPG-like term into the pressure Poisson equation. Under stated regularity and spatial-approximation assumptions, the authors derive a combined velocity error of order dt^2 + h^q for BDF2 time stepping, separating time-discretization, splitting, and spatial errors. Numerical comparisons indicate that omitting the pressure fine scale avoids the monolithic scheme's excess dissipation and drag overprediction at the cost of larger divergence error, while reducing per-step solve time by 1.3x to 2.7x in the 3D Taylor-Green test. If correct, this provides an equal-order, second-order-in-time alternative to monolithic stabilized solvers that is cheaper and less dissipative.

Core claim

The central claim, Proposition 1, is that the BDF2 incremental Helmholtz-Leray projection method with residual-based VMS stabilization on the predicted velocity satisfies a discrete velocity error of order dt^2 + h^q. The multiscale decomposition is applied only to the predicted velocity, whose unresolved part is modeled by the momentum residual; neither pressure nor corrected velocity is decomposed. The modeled fine scale enters all three substeps, producing SUPG-like stabilization in the predictor, a PSPG-like term in the pressure Poisson equation, and no pressure-fine-scale/grad-div term. On equal-order elements the scheme matches reference data for manufactured solutions, lid-driven cavi

What carries the argument

The load-bearing mechanism is the residual-based fine-scale model of the predicted velocity, ~u' = -tau_m * ~r_m, where ~r_m is the strong residual of the momentum predictor and tau_m is an elementwise stabilization time scale. This object supplies advective SUPG-type stabilization in the predictor, a PSPG-type divergence term in the pressure Poisson equation, and a correction contribution in the velocity projection. Around it sits the incremental Helmholtz-Leray splitting: a nonlinear advection-diffusion predictor with the previous pressure, a pressure Poisson equation for the pressure correction, and an L2 projection onto weakly divergence-free velocities. BDF2 gives the O(dt^2) temporal b

Load-bearing premise

The load-bearing premise is that the stabilized spatial discretization of the three fractional subproblems achieves the assumed error bound of order h^q; the proof does not derive this bound, so if the spatial accuracy assumption fails, the mesh part of the headline estimate collapses while the temporal O(dt^2) part remains.

What would settle it

Run the manufactured-solution test with a very fine fixed mesh and successively halve the time step at Re=10^6: the measured velocity error at the final time should fall at about dt^2. If the measured slope drops below 2, or if refining the mesh at fixed dt produces no improvement because the error plateaus, then respectively the temporal or spatial claim in Proposition 1 fails.

Watch this falsifier. Get emailed when new claim-graph text bears on it.

Share X Bluesky LinkedIn Reddit HN

If this is right

  • On equal-order elements, the scheme should show second-order temporal velocity convergence without any inf-sup-stable element pair, because pressure stability comes from the PSPG-like term rather than a pressure fine scale.
  • Monolithic VMS's pressure fine scale / grad-div term is the actionable source of excess modeled dissipation and cylinder-drag overprediction; removing it lets the projection scheme track DNS kinetic-energy decay on coarser meshes than the monolithic scheme.
  • The tradeoff is explicit: suppressing the pressure fine scale raises velocity-divergence error to about the level of monolithic VMS with p'=0, so users accepting this method accept larger pointwise divergence in exchange for better energy behavior and lower cost.
  • Because the coupled nonlinear saddle-point system becomes three sequential subproblems, the method reduces per-step wall-clock solve time by factors 1.3-2.7x in the 3D Taylor-Green tests under identical solver settings.
  • The total error estimate O(dt^2 + h^q) means time step and mesh should be refined together; refining only one leaves the other term dominating.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • If the grad-div term is indeed the sole cause of the monolithic scheme's excess dissipation, then a one-parameter family interpolating between full VMS and p'=0 should trace a smooth tradeoff between divergence error and over-dissipation; this is directly testable in the monolithic solver and would sharpen the attribution beyond the binary ablation.
  • The same splitting logic should transfer to rotational pressure correction, giving improved pressure accuracy while keeping the VMS machinery unchanged; the paper notes the rotational variant but does not combine it with the stabilized projection.
  • Because all three substeps are elliptic or nonlinear advection-diffusion and avoid saddle-point linear algebra, the measured cost advantage with the chosen Krylov solver is likely a lower bound; multigrid or GPU solvers on the substeps could widen the speedup without changing the formulation.
  • The formal h^q component rests entirely on Assumption 3.1, so proving that a standard SUPG/PSPG stability estimate delivers the assumed spatial consistency would turn Proposition 1 into a fully derived a priori bound rather than a conditional one.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

3 major / 4 minor

Summary. The paper proposes a VMS-stabilized incremental projection method for the incompressible Navier-Stokes equations. The multiscale decomposition is applied only to the predicted (intermediate) velocity; the pressure and the corrected velocity are not decomposed. The resulting scheme replaces the monolithic saddle-point system by a nonlinear momentum predictor, a pressure Poisson equation, and a velocity projection, with SUPG-like and PSPG-like terms arising from the modeled fine scales. A formal error analysis decomposes the total velocity error into BDF2 time-discretization, projection-splitting, and spatial-VMS contributions and states a combined O(dt^2 + h^q) estimate under regularity, stability, and spatial-approximation assumptions. Numerical experiments cover manufactured solutions, lid-driven cavity, flow past a cylinder, and the Taylor-Green vortex, and include comparisons with monolithic VMS and reference data, as well as wall-clock cost measurements.

Significance. If the central claims hold, the method is practically attractive: it enables equal-order velocity-pressure elements without the pressure fine scale, avoids the monolithic saddle point, and still supplies advective and pressure stabilization. The numerical verification is extensive and careful: the MMS tests confirm the expected second-order temporal rate across a wide Reynolds range, the cavity and cylinder benchmarks agree with established references, and the Taylor-Green study reports both accuracy and computational speedups. The temporal part of the error analysis is on solid ground, being based on the well-known incremental pressure-correction splitting. However, the spatial half of the headline estimate is not actually proved; it is an assumption that already contains the conclusion. The paper is transparent about its assumptions, but the formal result as stated is weaker than it appears.

major comments (3)
  1. [§3.1.3, Assumption 3.1 and Eq. (60)-(61)] The spatial part of the headline estimate is an assumption, not a derived result. Assumption 3.1 (Eq. (59)) assumes exactly the combined bound on the predicted-velocity error and the pressure-correction gradient error. Subtracting the continuous and discrete velocity-update equations and applying the triangle inequality (Eq. (60)) then yields the O(h^q) spatial error in Proposition 1. Thus the h^q term in (61) is Assumption 3.1 restated, not a proved consequence of the VMS stabilization or the fractional-step structure. The temporal O(dt^2) part is independently grounded, but the spatial rate remains a modeling hypothesis. Please either supply a proof (or a convincing derivation) of (59) for the nonlinear predictor and pressure Poisson equation under explicit regularity and SUPG/VMS hypotheses, or explicitly reframe Proposition 1 as a conditional statement in which the spatial rate is as
  2. [§2.3 and §2.4.2, Algorithm 1, Step 1] The derivation of the SUPG-type fine-scale advection term uses the assumption that the total velocity is divergence-free: the text after Eq. (23) states that the term (div(u_h+u'))u' is dropped because 'we have assumed that the total velocity u = u_h + u' is divergence free.' This assumption is not valid for the intermediate velocity of the projection scheme. The paper itself notes (before Eq. (29)) that the intermediate velocity is not generally solenoidal. Nevertheless, the same expansion is used in (43a) and in Algorithm 1, and the corresponding term involving div(tilde u_h + tilde u') is absent. Assumption 3.1 does not state or justify that this omitted term is O(h^q). This is a consistency gap in the derivation of the predictor and in its propagation to the pressure equation. The authors should either include the omitted term in the fine-scale model, prove that it is of higher order
  3. [§4.4 and Reference [49]] The Taylor-Green vortex DNS reference is cited as 'Wang et al., 2013 (512^3)' in Figures 12-15 and in the text. The cited reference [49] is 'High-order CFD methods: current status and perspective,' which does not appear to be a pseudo-spectral DNS dataset of the Taylor-Green vortex. As cited, this does not support the claim of agreement with a 512^3 pseudo-spectral DNS. If the data are taken from another source, please cite that source; if [49] does contain the DNS results used, please clarify the provenance in the text. This is important because the Taylor-Green comparison is a central validation claim for the method's accuracy.
minor comments (4)
  1. [§3.1.2 and §3.2] Lemma 3.1 is referred to as 'Theorem 3.1' in the paragraph following Eq. (58) and in the proof of Proposition 1; Remark 3.2 refers to 'Theorem 3.2'. Please correct the cross-references.
  2. [§3.1.2, Remark 3.2] The explanation of the splitting defect is a bit terse. The statement that the pressure increment phi^n is O(dt) for the incremental scheme should be justified or accompanied by a citation; currently it is presented as self-evident.
  3. [Algorithm 1, Step 1] The boundary integral in the penultimate line of the momentum predictor is said to vanish when fine scales are neglected on the boundary (Remark 2.1). This is consistent, but it would help to explicitly mark the term as '(vanishing by assumption)' in Algorithm 1, since the equation as written appears to retain a boundary contribution.
  4. [§4.3, Figures 10] The legend entry 'VMS/monolithic, p' = 0' uses the prime symbol ambiguously; it should be clarified that this is the monolithic solver with the pressure fine scale set to zero.

Circularity Check

2 steps flagged

Spatial half of Prop. 1 is Assumption 3.1 restated; predictor inherits invalid divergence-free fine-scale manipulation.

specific steps
  1. self definitional [Section 3.1.3 (Assumption 3.1, Eq. (59)) and Section 3.2 (Proposition 1, Eq. (61) and proof)]
    "We assume that the VMS discretizations of the three fractional subproblems are stable and consistent and that the total predicted velocity and pressure correction satisfy ||˜u^n − (˜u_h + ˜u′)|| + (1/σ)||∇(φ_n − φ_h)|| ≤ C h^q. ... Subtracting the continuous and discrete velocity-update equations and applying the triangle inequality together with Assumption 3.1 gives ||ˆu^n − ˆu_h|| ≤ C h^q. ... The spatial VMS term ||ˆu^n − ˆu_h|| is O(h^q) by (60)."

    The h^q component of Proposition 1 is not derived. Eq. (60) is exactly Assumption 3.1 after subtracting the velocity-update equations, and the proof of Proposition 1 invokes (60) to conclude the O(h^q) term in (61). Thus the spatial part of the headline error estimate is assumed in the same form in which it is output. Only the temporal O(dt^2) part is independent, coming from BDF2 and Shen/van Kan splitting. This is a partial circularity: the formal error estimate for the spatial rate is an assumption restated as a result.

  2. other [Section 2.3 (fine-scale advection expansion, before Remark 2.1) and Section 2.4.2 / Algorithm 1 (VMS predictor)]
    "where in the second step, we have assumed that the total velocity u = u_h + u′ is divergence free ... Note that the simplifications discussed in Remark 2.1 are also applied to ˜u_h in (43a)."

    The monolithic derivation drops the term ((∇·(u_h+u'))u', v_h) using div(u_h+u')=0. The projection method says the same simplifications apply to the predictor, but the intermediate velocity is not divergence-free: div(˜u_h+˜u') is the very source term in the pressure Poisson equation. The dropped term is therefore not zero, and no equivalent term appears in Algorithm 1. This is not a logical circle in itself, but it is the mechanism that lets the spatial error claim rest on Assumption 3.1: the omitted term is never shown to be O(h^q) and is silently absorbed into the assumed 'effects of the modeled velocity fine scale'.

full rationale

The paper is genuinely self-contained on the temporal side: Eq. (54) is standard BDF2 consistency, Lemma 3.1 is Shen's independent splitting estimate, and the manufactured-solution tests confirm second-order temporal convergence. The numerical benchmarks (Ghia, Wang, Schäfer/Kanaris, Taylor-Green) are external and not fitted. The only self-citation, [33], is a literature pointer and not load-bearing. The circularity is concentrated in the spatial error estimate: Assumption 3.1 (Eq. 59) already contains the total predicted-velocity and pressure-correction approximation at rate h^q; Eq. (60) is just that assumption applied to the velocity update; Proposition 1 then outputs the same h^q. No spatial convergence proof is supplied for the VMS closure, the nonlinear predictor, or the pressure Poisson step. In addition, the momentum predictor imports the monolithic 'total velocity divergence-free' manipulation although the intermediate velocity is not solenoidal; the missing term is folded into Assumption 3.1 rather than analyzed. These two points make the spatial half of the central claim conditional/assumed rather than derived, while the temporal half remains independent. Hence score 6, not higher.

Axiom & Free-Parameter Ledger

1 free parameters · 5 axioms · 0 invented entities

The central estimate rests on two nonstandard assumptions: Assumption 3.1 (the spatial error bound to be proven) and the divergence-free total velocity used to derive the SUPG term for a non-solenoidal predictor. The temporal machinery is external and standard. C_I is an unspecified input to the stabilization parameter.

free parameters (1)
  • C_I (dimensionless constant in tau_m) = not stated (literature VMS constant)
    tau_m in Eq. (21a)/(78) depends on C_I; the paper gives no chosen value, so the implemented stabilization parameter is not fully specified.
axioms (5)
  • ad hoc to paper Spatial approximation bound of Assumption 3.1: combined predicted-velocity and pressure-correction error is O(h^q) for the three VMS-discretized subproblems
    Eq. (59) states the spatial error that Proposition 1 then reports; it is the crux of the fully discrete estimate and is not proved.
  • ad hoc to paper Total velocity u_h+u' is divergence-free when deriving the SUPG-type fine-scale advection term
    Used in Section 2.3 before Eq. (24), then applied to the predictor ~u_h in (43a), where the total predicted velocity is not divergence-free; the mismatch is not discussed.
  • domain assumption Fine-scale velocity vanishes on the boundary and fine-scale time derivative and viscous contributions are neglected
    Remark 2.1, following Bazilevs et al. [4]; standard for algebraic VMS closures but an approximation.
  • domain assumption Sufficient smoothness of the exact NSE solution and a stable time-marching regime so per-step errors do not amplify
    Used in Section 3 before Lemma 3.1 and in the discrete Gronwall step of Proposition 1.
  • standard math BDF2 consistency and Shen's splitting error theorem for the non-rotational incremental pressure-correction scheme
    External results [11, 41, 24] provide the O(dt^2) velocity and O(dt) pressure temporal rates.

pith-pipeline@v1.3.0-alltime-deepseek · 25581 in / 14805 out tokens · 163014 ms · 2026-08-05T00:11:31.366914+00:00 · methodology

0 comments
Cite this review

Pith. "Pith review of A Helmholtz-Leray projection method with variational multiscale stabilization for the Navier-Stokes equations." pith.science (2026). https://pith.science/paper/PBKC6ULH

@misc{pith2026260800843,
  author       = {Pith},
  title        = {Pith review of: A Helmholtz-Leray projection method with variational multiscale stabilization for the Navier-Stokes equations},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/PBKC6ULH}},
  note         = {Machine review of arXiv:2608.00843}
}
Share X Bluesky LinkedIn Reddit HN
read the original abstract

The Galerkin finite element formulation of the incompressible Navier-Stokes equations presents two principal challenges: maintaining stable velocity-pressure coupling and controlling instability in advection-dominated regimes. Moreover, the monolithic formulation produces a coupled nonlinear saddle-point problem. In this work, we present a residual-based variational multiscale (VMS) stabilization of an incremental Helmholtz-Leray projection method that replaces this coupled saddle-point problem with a nonlinear velocity predictor, a pressure Poisson equation, and a velocity projection. The multiscale decomposition is applied only to the predicted velocity; neither the pressure nor the corrected, weakly divergence-free velocity is decomposed into coarse and fine scales. The modeled velocity fine scale contributes consistently to all three subproblems, introducing SUPG-like stabilization in the momentum predictor and a PSPG-like residual contribution in the pressure Poisson equation. We provide a formal error decomposition that separates the BDF2 time-discretization, projection-splitting, and spatial-VMS errors and, under stated stability and spatial-approximation assumptions, yields a combined velocity error estimate with second-order temporal accuracy. Numerical results for manufactured solutions, lid-driven cavity flow, flow past a cylinder, and the Taylor-Green vortex agree closely with established reference data. Comparisons with monolithic VMS indicate that omitting the pressure fine scale reduces drag overprediction and excess modeled dissipation, at the cost of increased divergence error. In the Taylor-Green tests, the projection formulation also reduces the average solution time per step by factors ranging from approximately $1.3\times$ to $2.7\times$ under identical solver settings.

Figures

Figures reproduced from arXiv: 2608.00843 by Baskar Ganapathysubramanian, Biswajit Khara, Makrand Khanwale, Suresh Murugaiyan.

Figure 1
Figure 1. Figure 1: Sketch of a general domain with different conditions applied at different boundaries. Neu [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: A schematic diagram of the path to the VMS-stabilized projection method ( [PITH_FULL_IMAGE:figures/full_fig_p014_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Convergence of velocity and pressure errors, [PITH_FULL_IMAGE:figures/full_fig_p018_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Boundary conditions for the lid-driven cavity problem. The blue dashed lines are the [PITH_FULL_IMAGE:figures/full_fig_p019_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Midline velocity profiles for the lid-driven cavity problem. Left column shows [PITH_FULL_IMAGE:figures/full_fig_p019_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: An example of divergence fields at steady state for the lid-driven cavity at [PITH_FULL_IMAGE:figures/full_fig_p020_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Velocity correction applied by the projection step for the lid-driven cavity at [PITH_FULL_IMAGE:figures/full_fig_p020_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Schematic diagram of the flow past a circular cylinder in 2D. The inlet velocity profile [PITH_FULL_IMAGE:figures/full_fig_p021_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Instantaneous vorticity contours for two-dimensional flow past a circular cylinder at Reynolds [PITH_FULL_IMAGE:figures/full_fig_p021_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Results for flow past a circular cylinder in 2D: drag coefficient [PITH_FULL_IMAGE:figures/full_fig_p023_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Taylor-Green vortex at Re = 1600 from the projection-based RBVMS scheme on the 1283 mesh, visualized using Q-criterion iso-surfaces at Q = 0.3, 0.5, 0.8, coloured by velocity magnitude. (Left) Initial condition at t = 0, showing the organized cellular structure. (Centre) Developed turbu￾lent state at t = 15, after the breakdown of the initial vortex sheets. (Right) Decaying turbulent state at t = 20, show… view at source ↗
Figure 12
Figure 12. Figure 12: Taylor–Green vortex at Re = 1600 on three mesh resolutions (643 , 1283 , 2563 ), compared against the pseudo-spectral DNS of Wang et al. [49] on a 5123 grid. Top row: monolithic VMS scheme. Bottom row: VMS-stabilized projection scheme. Left column: volume-averaged kinetic energy Ek(t). Right column: energy-based dissipation rate εk(t) [PITH_FULL_IMAGE:figures/full_fig_p025_12.png] view at source ↗
Figure 13
Figure 13. Figure 13: Three-dimensional kinetic-energy spectrum [PITH_FULL_IMAGE:figures/full_fig_p026_13.png] view at source ↗
Figure 14
Figure 14. Figure 14: Energy-based (εk, solid) and enstrophy-based (εΩ, dashed) dissipation rates for the Taylor– Green vortex at Re = 1600. Top row: monolithic VMS solver. Bottom row: VMS/projection solver. Columns correspond to the 643 , 1283 , and 2563 meshes. The shaded region indicates the gap εk − εΩ, representing the contribution of unresolved scales (model and numerical dissipation) to the total kinetic-energy decay. O… view at source ↗
Figure 15
Figure 15. Figure 15: Effect of the continuity (grad–div) stabilization on the Taylor–Green vortex at [PITH_FULL_IMAGE:figures/full_fig_p027_15.png] view at source ↗
Figure 16
Figure 16. Figure 16: Average wall-clock solve time per time step for the Taylor-Green vortex at [PITH_FULL_IMAGE:figures/full_fig_p028_16.png] view at source ↗

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

50 extracted references · 41 canonical work pages

  1. [1]

    A stable finite element for the stokes equations.Calcolo, 21(4):337–344, 1984

    Douglas N Arnold, Franco Brezzi, and Michel Fortin. A stable finite element for the stokes equations.Calcolo, 21(4):337–344, 1984

    work page 1984

  2. [2]

    The finite element method with lagrangian multipliers.Numerische Mathematik, 20(3):179–192, 1973

    Ivo Babuška. The finite element method with lagrangian multipliers.Numerische Mathematik, 20(3):179–192, 1973

    work page 1973

  3. [3]

    Claudio Baiocchi, Franco Brezzi, and Leopoldo P. Franca. Virtual bubbles and galerkin-least- squares type methods.Computer Methods in Applied Mechanics and Engineering, 105(1):125–141, 1993

    work page 1993

  4. [4]

    Y Bazilevs, VM Calo, JA Cottrell, TJR Hughes, A Reali, and G Scovazzi. Variational multiscale residual-based turbulence modeling for large eddy simulation of incompressible flows.Computer methods in applied mechanics and engineering, 197(1-4):173–201, 2007

    work page 2007

  5. [5]

    Local projection stabilization for the Oseen problem and its interpretation as a variational multiscale method.SIAM Journal on Numerical Analysis, 43(6): 2544–2566, 2006

    Malte Braack and Erik Burman. Local projection stabilization for the Oseen problem and its interpretation as a variational multiscale method.SIAM Journal on Numerical Analysis, 43(6): 2544–2566, 2006. doi: 10.1137/050631227

    doi:10.1137/050631227 2006

  6. [6]

    On the existence, uniqueness and approximation of saddle-point problems arising from lagrangian multipliers.Publications mathématiques et informatique de Rennes, S4:1–26, 1974

    Franco Brezzi. On the existence, uniqueness and approximation of saddle-point problems arising from lagrangian multipliers.Publications mathématiques et informatique de Rennes, S4:1–26, 1974

    work page 1974

  7. [7]

    Springer Science & Business Media, 2012

    Franco Brezzi and Michel Fortin.Mixed and hybrid finite element methods, volume 15. Springer Science & Business Media, 2012

    work page 2012

  8. [8]

    Choosing bubbles for advection-diffusion problems.Math- ematical Models and Methods in Applied Sciences, 4(4):571–587, 1994

    Franco Brezzi and Alessandro Russo. Choosing bubbles for advection-diffusion problems.Math- ematical Models and Methods in Applied Sciences, 4(4):571–587, 1994

    work page 1994

  9. [9]

    Alexander N Brooks and Thomas JR Hughes. Streamline upwind/petrov-galerkin formulations for convection dominated flows with particular emphasis on the incompressible navier-stokes equations.Computer methods in applied mechanics and engineering, 32(1-3):199–259, 1982

    work page 1982

  10. [10]

    Edge stabilization for Galerkin approximations of convection– diffusion–reaction problems.Computer Methods in Applied Mechanics and Engineering, 193(15– 16):1437–1453, 2004

    Erik Burman and Peter Hansbo. Edge stabilization for Galerkin approximations of convection– diffusion–reaction problems.Computer Methods in Applied Mechanics and Engineering, 193(15– 16):1437–1453, 2004. doi: 10.1016/j.cma.2004.01.030

    work page doi:10.1016/j.cma.2004.01.030 2004

  11. [11]

    John Wiley & Sons, 2016

    John Charles Butcher.Numerical methods for ordinary differential equations. John Wiley & Sons, 2016

    2016

  12. [12]

    Alexandre J. Chorin. Numerical solution of the navier–stokes equations.Mathematics of Com- putation, 22(104):745–762, 1968. 29

    work page 1968

  13. [13]

    Runge–Kutta discontinuous Galerkin methods for convection-dominated problems.Journal of Scientific Computing, 16(3):173–261, 2001

    Bernardo Cockburn and Chi-Wang Shu. Runge–Kutta discontinuous Galerkin methods for convection-dominated problems.Journal of Scientific Computing, 16(3):173–261, 2001. doi: 10.1023/A:1012873910884

    doi:10.1023/a:1012873910884 2001

  14. [14]

    Comparison of some finite element methods for solving the diffusion-convection- reaction equation.Computer methods in applied mechanics and engineering, 156(1-4):185–210, 1998

    Ramon Codina. Comparison of some finite element methods for solving the diffusion-convection- reaction equation.Computer methods in applied mechanics and engineering, 156(1-4):185–210, 1998

    work page 1998

  15. [15]

    Ramon Codina. Stabilization of incompressibility and convection through orthogonal sub-scales in finite element methods.Computer methods in applied mechanics and engineering, 190(13-14): 1579–1599, 2000

    work page 2000

  16. [16]

    Ramon Codina. Stabilized finite element approximation of transient incompressible flows using orthogonal subscales.Computer methods in applied mechanics and engineering, 191(39-40):4295– 4321, 2002

    work page 2002

  17. [17]

    Conforming and nonconforming finite element methods for solving the stationary stokes equations i.Revue française d’automatique informatique recherche opérationnelle

    Michel Crouzeix and P-A Raviart. Conforming and nonconforming finite element methods for solving the stationary stokes equations i.Revue française d’automatique informatique recherche opérationnelle. Mathématique, 7(R3):33–75, 1973

    work page 1973

  18. [18]

    John Wiley & Sons, 2003

    Jean Donea and Antonio Huerta.Finite element methods for flow problems. John Wiley & Sons, 2003

    2003

  19. [19]

    Interior penalty procedures for elliptic and parabolic Galerkin methods

    Jim Douglas and Todd Dupont. Interior penalty procedures for elliptic and parabolic Galerkin methods. In R. Glowinski and J.-L. Lions, editors,Computing Methods in Applied Sciences, volume 58 ofLecture Notes in Physics, pages 207–216. Springer, Berlin, Heidelberg, 1976. doi: 10.1007/BFb0120591

    work page doi:10.1007/bfb0120591 1976

  20. [20]

    An absolutely stabilized finite element method for the stokes problem.Mathematics of computation, 52(186):495–508, 1989

    Jim Douglas and Jun Ping Wang. An absolutely stabilized finite element method for the stokes problem.Mathematics of computation, 52(186):495–508, 1989

    work page 1989

  21. [21]

    Cambridge University Press, 2001

    Ciprian Foias, Oscar Manley, Ricardo Rosa, and Roger Temam.Navier–Stokes Equations and Turbulence. Cambridge University Press, 2001

    work page 2001

  22. [22]

    High-re solutions for incompressible flow using the navier-stokes equations and a multigrid method.Journal of computational physics, 48(3):387–411, 1982

    UKNG Ghia, Kirti N Ghia, and CT Shin. High-re solutions for incompressible flow using the navier-stokes equations and a multigrid method.Journal of computational physics, 48(3):387–411, 1982

    1982

  23. [23]

    A multistep technique with implicit difference schemes for calculating two- or three-dimensional cavity flows.Journal of Computational Physics, 30(1):76–95, 1979

    Katuhiko Goda. A multistep technique with implicit difference schemes for calculating two- or three-dimensional cavity flows.Journal of Computational Physics, 30(1):76–95, 1979

    work page 1979

  24. [24]

    An overview of projection methods for incom- pressible flows.Computer Methods in Applied Mechanics and Engineering, 195(44–47):6011–6045, 2006

    Jean-Luc Guermond, Peter Minev, and Jie Shen. An overview of projection methods for incom- pressible flows.Computer Methods in Applied Mechanics and Engineering, 195(44–47):6011–6045, 2006

    work page 2006

  25. [25]

    Heywood and Rolf Rannacher

    John G. Heywood and Rolf Rannacher. Finite element approximation of the nonstationary Navier–Stokes problem. IV. error analysis for second-order time discretization.SIAM Journal on Numerical Analysis, 27(2):353–384, 1990. doi: 10.1137/0727022

    doi:10.1137/0727022 1990

  26. [26]

    Thomas JR Hughes. Multiscale phenomena: Green’s functions, the dirichlet-to-neumann formu- lation, subgrid scale models, bubbles and the origins of stabilized methods.Computer methods in applied mechanics and engineering, 127(1-4):387–401, 1995

    work page 1995

  27. [27]

    A new finite element formulation for computational fluid dynamics: V

    Thomas JR Hughes, Leopoldo P Franca, and Marc Balestra. A new finite element formulation for computational fluid dynamics: V. circumventing the babuška-brezzi condition: A stable petrov- galerkin formulation of the stokes problem accommodating equal-order interpolations.Computer Methods in Applied Mechanics and Engineering, 59(1):85–99, 1986. 30

    work page 1986

  28. [28]

    A new finite element formu- lation for computational fluid dynamics: Viii

    Thomas JR Hughes, Leopoldo P Franca, and Gregory M Hulbert. A new finite element formu- lation for computational fluid dynamics: Viii. the galerkin/least-squares method for advective- diffusive equations.Computer methods in applied mechanics and engineering, 73(2):173–189, 1989

    work page 1989

  29. [29]

    Large eddy simulation and the varia- tional multiscale method.Computing and visualization in science, 3(1):47–59, 2000

    Thomas JR Hughes, Luca Mazzei, and Kenneth E Jansen. Large eddy simulation and the varia- tional multiscale method.Computing and visualization in science, 3(1):47–59, 2000

    work page 2000

  30. [30]

    Springer, 2016

    Volker John.Finite element methods for incompressible flow problems. Springer, 2016

    work page 2016

  31. [31]

    Afiniteelementvariationalmultiscalemethodforthenavier–stokes equations.SIAM Journal on Scientific Computing, 26(5):1485–1503, 2005

    VolkerJohnandSongulKaya. Afiniteelementvariationalmultiscalemethodforthenavier–stokes equations.SIAM Journal on Scientific Computing, 26(5):1485–1503, 2005

    work page 2005

  32. [32]

    Three dimensional flow around a circular cylinder confined in a plane channel.Physics of fluids, 23(6), 2011

    Nicolas Kanaris, Dimokratis Grigoriadis, and Stavros Kassinos. Three dimensional flow around a circular cylinder confined in a plane channel.Physics of fluids, 23(6), 2011

    work page 2011

  33. [33]

    A semi-implicit variational multiscale formulation for the incompressible navier-stokes equations via exact adjoint linearization.arXiv preprint arXiv:2512.21773, 2025

    Biswajit Khara, Suresh Murugaiyan, Suriya Dhakshinamoorthy, Makrand Khanwale, Ming-Chen Hsu, and Baskar Ganapathysubramanian. A semi-implicit variational multiscale formulation for the incompressible navier-stokes equations via exact adjoint linearization.arXiv preprint arXiv:2512.21773, 2025

    arXiv 2025

  34. [34]

    Flux correction tools for finite elements.Journal of Compu- tational Physics, 175(2):525–558, 2002

    Dmitri Kuzmin and Stefan Turek. Flux correction tools for finite elements.Journal of Compu- tational Physics, 175(2):525–558, 2002. doi: 10.1006/jcph.2001.6955

    arXiv 2002

  35. [35]

    Blair Perot

    J. Blair Perot. An analysis of the fractional step method.Journal of Computational Physics, 108 (1):51–58, 1993. doi: 10.1006/jcph.1993.1162

    work page arXiv 1993

  36. [36]

    Pope.Turbulent Flows

    Stephen B. Pope.Turbulent Flows. Cambridge University Press, Cambridge, 2000. ISBN 9780521598866. doi: 10.1017/CBO9780511840531

    doi:10.1017/cbo9780511840531 2000

  37. [37]

    Factorization methods for the nu- merical approximation of Navier–Stokes equations.Computer Methods in Applied Mechanics and Engineering, 188(1–3):505–526, 2000

    Alfio Quarteroni, Fausto Saleri, and Alessandro Veneziani. Factorization methods for the nu- merical approximation of Navier–Stokes equations.Computer Methods in Applied Mechanics and Engineering, 188(1–3):505–526, 2000. doi: 10.1016/S0045-7825(99)00192-9

    work page doi:10.1016/s0045-7825(99)00192-9 2000

  38. [38]

    Springer, 1996

    Michael Schäfer, Stefan Turek, Franz Durst, Egon Krause, and Rolf Rannacher.Benchmark computations of laminar flow around a cylinder. Springer, 1996

    work page 1996

  39. [39]

    Ridgway Scott and Michael Vogelius

    L. Ridgway Scott and Michael Vogelius. Norm estimates for a maximal right inverse of the divergence operator in spaces of piecewise polynomials.ESAIM: Mathematical Modelling and Numerical Analysis, 19(1):111–143, 1985. URLhttp://eudml.org/doc/193439

    work page 1985

  40. [40]

    A new finite element formulation for computational fluid dynamics: X

    Farzin Shakib, Thomas JR Hughes, and Zdeněk Johan. A new finite element formulation for computational fluid dynamics: X. the compressible euler and navier-stokes equations.Computer Methods in Applied Mechanics and Engineering, 89(1-3):141–219, 1991

    work page 1991

  41. [41]

    On error estimates of the projection methods for the navier-stokes equations: second- order schemes.Mathematics of computation, 65(215):1039–1065, 1996

    Jie Shen. On error estimates of the projection methods for the navier-stokes equations: second- order schemes.Mathematics of computation, 65(215):1039–1065, 1996

    work page 1996

  42. [42]

    Flow past a circular cylinder between parallel walls at low reynolds numbers.Ocean Engineering, 37(8-9):757–769, 2010

    Sintu Singha and KP Sinhamahapatra. Flow past a circular cylinder between parallel walls at low reynolds numbers.Ocean Engineering, 37(8-9):757–769, 2010

    work page 2010

  43. [43]

    A numerical solution of the navier-stokes equations using the finite element technique.Computers & Fluids, 1(1):73–100, 1973

    Cedric Taylor and Paul Hood. A numerical solution of the navier-stokes equations using the finite element technique.Computers & Fluids, 1(1):73–100, 1973

    work page 1973

  44. [44]

    Sur l’approximation de la solution des équations de navier–stokes par la méthode des pas fractionnaires.Archive for Rational Mechanics and Analysis, 32:135–153, 1969

    Roger Temam. Sur l’approximation de la solution des équations de navier–stokes par la méthode des pas fractionnaires.Archive for Rational Mechanics and Analysis, 32:135–153, 1969

    work page 1969

  45. [45]

    Tezduyar and Yong J

    Tayfun E. Tezduyar and Yong J. Park. Discontinuity-capturing finite element formulations for nonlinear convection-diffusion-reaction equations.Computer Methods in Applied Mechanics and Engineering, 59(3):307–325, 1986. doi: 10.1016/0045-7825(86)90003-4. 31

    work page doi:10.1016/0045-7825(86)90003-4 1986

  46. [46]

    Tayfun E Tezduyar, Sanjay Mittal, SE Ray, and R Shih. Incompressible flow computations with stabilized bilinear and linear equal-order-interpolation velocity-pressure elements.Computer Methods in Applied Mechanics and Engineering, 95(2):221–242, 1992

    work page 1992

  47. [47]

    L. J. P. Timmermans, P. D. Minev, and F. N. van de Vosse. An approximate projection scheme for incompressible flow using spectral elements.International Journal for Numerical Methods in Fluids, 22(7):673–688, 1996

    work page 1996

  48. [48]

    J. van Kan. A second-order accurate pressure-correction scheme for viscous incompressible flow. SIAM Journal on Scientific and Statistical Computing, 7(3):870–891, 1986

    work page 1986

  49. [49]

    High-order cfd methods: current status and perspective.International Journal for Numerical Methods in Fluids, 72(8):811–845, 2013

    Zhijian J Wang, Krzysztof Fidkowski, Rémi Abgrall, Francesco Bassi, Doru Caraeni, Andrew Cary, Herman Deconinck, Ralf Hartmann, Koen Hillewaert, Hung T Huynh, et al. High-order cfd methods: current status and perspective.International Journal for Numerical Methods in Fluids, 72(8):811–845, 2013

    work page 2013

  50. [50]

    Flow about a circular cylinder between parallel walls

    Luigino Zovatto and Gianni Pedrizzetti. Flow about a circular cylinder between parallel walls. Journal of Fluid Mechanics, 440:1–25, 2001. 32 A Element metric quantities for the VMS parameters Letx=x(ξ)denote the map from the reference element, with coordinatesξ, to a physical element, with coordinatesx. Define the inverse-Jacobian components by Kij = ∂ξi...

    work page 2001