pith. sign in

arxiv: 2607.00686 · v1 · pith:WI3WYJI5new · submitted 2026-07-01 · 🧮 math.NA · cs.NA

Goal-oriented space-time adaptivity for the Navier--Stokes equations based on the dual weighted residual method

Marius Paul Bruchh\"auser , Nils Margenberg , Markus Bause This is my paper

Pith reviewed 2026-07-02 07:54 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords Navier-Stokes equationsdual weighted residual methodgoal-oriented adaptivityspace-time finite elementsa posteriori error estimationdiscontinuous Galerkingeometric multigridNewton-GMRES
0
0 comments X

The pith

The dual weighted residual method drives goal-oriented space-time mesh adaptivity to control specific target quantities in Navier-Stokes simulations on feasible meshes.

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

This paper develops a goal-oriented a posteriori error estimator based on the dual weighted residual method and pairs it with space-time adaptivity for the Navier-Stokes equations. The discretization uses discontinuous Galerkin in time and inf-sup stable finite element pairs on tensor-product meshes, with nonlinear systems solved by Newton iteration and GMRES preconditioned by slab-wise geometric multigrid. The central claim is that this combination produces reliable error control for chosen quantities of interest while keeping the computational meshes tractable and the algebraic solves robust. A reader would care because many fluid applications require accuracy only in particular outputs rather than everywhere in the domain. Benchmark tests are used to check accuracy, efficiency, and stability of the overall approach.

Core claim

The paper claims that the dual weighted residual method supplies a reliable a posteriori estimator for chosen target functionals on the nonlinear Navier-Stokes system when discretized with the given DG-in-time and inf-sup stable space-time finite element scheme; when this estimator is used to drive space-time mesh adaptivity and the resulting algebraic systems are solved by Newton-GMRES with slab-wise geometric multigrid, the method yields reliable control of the target quantities on computationally feasible meshes together with robust and efficient algebraic solution.

What carries the argument

Dual weighted residual (DWR) method for a posteriori error estimation that drives goal-oriented space-time mesh adaptivity

If this is right

  • Reliable control of target quantities on computationally feasible space-time meshes
  • Robust and efficient solution of the nonlinear algebraic systems via Newton iteration with GMRES and slab-wise geometric multigrid
  • Performance demonstrated through benchmark computations with respect to accuracy, efficiency, and stability
  • MPI-parallel implementation on tensor-product meshes using discontinuous Galerkin time discretization and inf-sup stable elements

Where Pith is reading between the lines

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

  • The same DWR-driven adaptivity could be applied to other nonlinear time-dependent systems where only a few output functionals matter.
  • Because the method separates the choice of target functional from the mesh generation, it may simplify parameter studies or optimization loops that repeatedly solve the flow equations.
  • The reliance on slab-wise multigrid suggests that parallel scalability could remain good when the number of time slabs grows with problem size.

Load-bearing premise

The dual weighted residual method supplies a reliable a posteriori estimator for the chosen target functionals on the nonlinear Navier-Stokes system discretized with the given DG-in-time and inf-sup stable space-time finite element scheme.

What would settle it

A computation in which the DWR estimator fails to bound the true error in the target functional, or in which the adapted meshes require more degrees of freedom than uniform refinement to reach a prescribed tolerance on that functional.

Figures

Figures reproduced from arXiv: 2607.00686 by Marius Paul Bruchh\"auser, Markus Bause, Nils Margenberg.

Figure 5.1
Figure 5.1. Figure 5.1: Algebraic performance over the adaptive loops: nonlinear iterations per [PITH_FULL_IMAGE:figures/full_fig_p019_5_1.png] view at source ↗
Figure 5.2
Figure 5.2. Figure 5.2: Share of time spent on the primal problem, dual problem, and error [PITH_FULL_IMAGE:figures/full_fig_p020_5_2.png] view at source ↗
Figure 5.3
Figure 5.3. Figure 5.3: Distribution of the adaptively determined temporal step size [PITH_FULL_IMAGE:figures/full_fig_p020_5_3.png] view at source ↗
Figure 5.4
Figure 5.4. Figure 5.4: Distribution of the adaptively determined minimal spatial element di [PITH_FULL_IMAGE:figures/full_fig_p021_5_4.png] view at source ↗
Figure 5.5
Figure 5.5. Figure 5.5: Geometry of the test scenario with a parabolic inflow profile Γ [PITH_FULL_IMAGE:figures/full_fig_p021_5_5.png] view at source ↗
Figure 5.6
Figure 5.6. Figure 5.6: Drag over all DWR loops Finally, in [PITH_FULL_IMAGE:figures/full_fig_p023_5_6.png] view at source ↗
Figure 5.7
Figure 5.7. Figure 5.7: Adaptive refinement history and final-loop solution at [PITH_FULL_IMAGE:figures/full_fig_p024_5_7.png] view at source ↗
read the original abstract

This work presents a goal-oriented a posteriori error estimator based on the Dual Weighted Residual (DWR) method together with space-time mesh adaptivity for the Navier--Stokes equations. The resulting nonlinear algebraic systems on the space-time slabs are solved by Newton's method with GMRES, preconditioned by a slab-wise geometric multigrid method. This combination yields reliable control of target quantities on computationally feasible space-time meshes together with a robust and efficient solution of the algebraic systems. The implementation is based on a MPI-parallel programming model in the deal.II library. Further ingredients are a discontinuous Galerkin discretization in time and inf-sup stable finite element pairs with discontinuous pressure on tensor-product meshes. The performance of the approach is investigated in benchmark computations with regard to accuracy, efficiency, and stability.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

0 major / 3 minor

Summary. The manuscript presents a goal-oriented a posteriori error estimator for the Navier-Stokes equations based on the dual weighted residual (DWR) method, combined with space-time mesh adaptivity. It employs a discontinuous Galerkin discretization in time and inf-sup stable finite element pairs with discontinuous pressure on tensor-product meshes. Nonlinear algebraic systems on space-time slabs are solved via Newton's method with GMRES, preconditioned by a slab-wise geometric multigrid method, all implemented in a MPI-parallel deal.II framework. The performance is assessed via benchmark computations with respect to accuracy, efficiency, and stability, with the central assertion that the approach yields reliable control of target quantities on feasible meshes together with robust algebraic solves.

Significance. If the DWR estimator is shown to be reliable for the nonlinear Navier-Stokes system under the stated discretization, the work provides a concrete framework for goal-oriented space-time adaptivity that integrates error control with efficient multigrid solvers. This could be useful for applications where specific output functionals (e.g., drag or lift) must be computed accurately without excessive mesh refinement. The parallel implementation and focus on slab-wise adaptivity address practical computational constraints in time-dependent flows.

minor comments (3)
  1. Abstract, final paragraph: the statement that 'benchmark computations investigate accuracy, efficiency, and stability' is not accompanied by any quantitative indicators (e.g., effectivity indices, error reductions, or CPU-time savings). Adding a one-sentence summary of key numerical outcomes would strengthen the abstract without altering length substantially.
  2. The manuscript refers to 'the chosen target functionals' but does not list them explicitly in the abstract or early introduction. A short enumeration of the functionals used in the benchmarks (e.g., in §4 or §5) would improve readability for readers interested in specific quantities of interest.
  3. Implementation details: while the use of deal.II and MPI is noted, the paper would benefit from a brief statement on how the DWR weights are computed and stored across slabs (e.g., in the section describing the adaptive algorithm) to aid reproducibility.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive summary of the manuscript, recognition of its significance, and recommendation for minor revision. No specific major comments are listed in the report.

Circularity Check

0 steps flagged

No significant circularity detected

full rationale

The paper applies the established dual weighted residual (DWR) method to derive a goal-oriented a posteriori estimator for the Navier-Stokes equations under a discontinuous Galerkin-in-time and inf-sup stable space-time finite element discretization. The central claims concern the combination of this estimator with slab-wise adaptivity and geometric multigrid solvers for the resulting nonlinear systems; these steps rely on standard techniques from the literature rather than any self-referential definition, fitted parameter renamed as prediction, or load-bearing self-citation chain. No equation or claim reduces the target control result to an input by construction, and the numerical benchmarks provide independent verification outside the derivation itself.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the standard assumption that the DWR framework extends reliably to the nonlinear space-time NS problem with the chosen discretization; no free parameters or new entities are mentioned.

axioms (1)
  • domain assumption The dual weighted residual method yields a reliable a posteriori error estimator for target functionals of the Navier-Stokes equations under the given discretization.
    Invoked as the basis for the goal-oriented adaptivity in the first sentence of the abstract.

pith-pipeline@v0.9.1-grok · 5670 in / 1233 out tokens · 28270 ms · 2026-07-02T07:54:23.816803+00:00 · methodology

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

44 extracted references

  1. [1]

    Pasquale C. Africa, Daniel Arndt, Wolfgang Bangerth, Bruno Blais, Marc Fehling, Rene Gassm¨ oller, Timo Heister, Luca Heltai, Sebastian Kinnewig, Martin Kronbichler, Matthias Maier, Peter Munch, Magdalena Schreter- Fleischhacker, Jan P. Thiele, Bruno Turcksin, David Wells, and Vladimir Yushutin. The deal.II library, version 9.6.Journal of Numerical Mathem...

    2024

  2. [2]

    Mathias Anselmann and Markus Bause. A geometric multigrid method for space- time finite element discretizations of the navier–stokes equations and its appli- cation to 3d flow simulation.ACM Transactions on Mathematical Software, 49(1):1–25, 2023

    2023

  3. [3]

    Bangerth, M

    W. Bangerth, M. Geiger, and R. Rannacher. Adaptive galerkin finite element methods for the wave equation.Computational Methods in Applied Mathematics, 10(1):3–48, 2010

    2010

  4. [4]

    Bangerth and R

    W. Bangerth and R. Rannacher.Adaptive Finite Element Methods for Differen- tial Equations. Birkh¨ auser Verlag, Basel, 2003

    2003

  5. [5]

    Bause, M

    M. Bause, M. P. Bruchh¨ auser, B. Endtmayer, N. Margenberg, I. Toulopoulos, and T. Wick. Anisotropic space-time goal-oriented error control and mesh adap- tivity for convection-diffusion-reaction equations.J. Numer. Math, 2025. 25

    2025

  6. [6]

    Becker and R

    R. Becker and R. Rannacher. An optimal control approach to a posteriori error estimation in finite element methods.Acta Numer., 10:1–102, 2001

    2001

  7. [7]

    Besier and R

    M. Besier and R. Rannacher. Goal-oriented space-time adaptivity in the finite element Galerkin method for the computation of nonstationary incompressible flow.Int. J. Num. Methods Fluids, 70(9):1139–1166, 2012

    2012

  8. [8]

    Michael Besier and Winnifried Wollner. On the pressure approximation in nonstationary incompressible flow simulations on dynamically varying spatial meshes.International Journal for Numerical Methods in Fluids, 69(6):1045– 1064, 2012

    2012

  9. [9]

    On the pressure approximation in nonstationary incompressible flow simulations on dynamically varying spa- tial meshes: PRESSURE APPROXIMATION ON DYNAMIC MESHES

    Michael Besier and Winnifried Wollner. On the pressure approximation in nonstationary incompressible flow simulations on dynamically varying spa- tial meshes: PRESSURE APPROXIMATION ON DYNAMIC MESHES. 69(6):1045–1064, 2012

    2012

  10. [10]

    Adaptive computation of reactive flows with local mesh refinement and model adaptation

    Malte Braack and Alexandre Ern. Adaptive computation of reactive flows with local mesh refinement and model adaptation. pages 159–168. Springer, 2004

    2004

  11. [11]

    A priori error analysis for finite element approximations of the stokes problem on dynamic meshes.IMA Journal of Numerical Analysis, 34(1):123–146, 05 2013

    Andreas Brenner, Eberhard B¨ ansch, and Markus Bause. A priori error analysis for finite element approximations of the stokes problem on dynamic meshes.IMA Journal of Numerical Analysis, 34(1):123–146, 05 2013

    2013

  12. [12]

    M. P. Bruchh¨ auser.Goal-oriented space-time adaptivity for a multirate approach to coupled flow and transport. PhD thesis, Helmut-Schmidt- University/University of the German Federal Armed Forces Hamburg, 2022

    2022

  13. [13]

    Numerical study of approximation techniques for the temporal weights to the dwr method

    Marius Paul Bruchh¨ auser and Markus Bause. Numerical study of approximation techniques for the temporal weights to the dwr method. In Ad´ elia Sequeira, Ana Silvestre, Svilen S. Valtchev, and Jo˜ ao Janela, editors,Numerical Mathematics and Advanced Applications ENUMATH 2023, Volume 1, pages 177–187, Cham. Springer Nature Switzerland

    2023

  14. [14]

    On the implementa- tion of an adaptive multirate framework for coupled transport and flow.Journal of Scientific Computing, 93(59)

    Marius Paul Bruchh¨ auser, Uwe K¨ ocher, and Markus Bause. On the implementa- tion of an adaptive multirate framework for coupled transport and flow.Journal of Scientific Computing, 93(59)

  15. [15]

    Carey and J.T

    G.F. Carey and J.T. Oden.Finite elements, computational aspects, Vol. III. The Texas finite element series. Prentice-Hall, Englewood Cliffs, New Jersey, 1984

    1984

  16. [16]

    Alexandre Ern, Iain Smears, and Martin Vohral´ ık. Guaranteed, locally space- time efficient, and polynomial-degree robust a posteriori error estimates for high- order discretizations of parabolic problems.SIAM Journal on Numerical Anal- ysis, 55(6):2811–2834, 2017

    2017

  17. [17]

    R. D. Falgout, S. Friedhoff, Tz. V. Kolev, S. P. MacLachlan, J. B. Schroder, and S. Vandewalle. Multigrid methods with space–time concurrency.Computing and Visualization in Science, 18(4-5):123–143, 8 2017

    2017

  18. [18]

    Towards Textbook Efficiency for Parallel Multigrid.Numerical Math- ematics Theory Methods and Applications, 8(1):22–46, 2 2015

    Bj¨ orn Gmeiner, Ulrich R¨ ude, Holger Stengel, Christian Waluga, and Barbara Wohlmuth. Towards Textbook Efficiency for Parallel Multigrid.Numerical Math- ematics Theory Methods and Applications, 8(1):22–46, 2 2015. 26

    2015

  19. [19]

    Springer-Verlag, Berlin Heidelberg New York Tokyo, 1985

    Wolfgang Hackbusch.Multi-Grid Methods and Applications, volume 4 ofSpringer Series in Computational Mathematics. Springer-Verlag, Berlin Heidelberg New York Tokyo, 1985

    1985

  20. [20]

    Heywood and Rolf Rannacher

    John G. Heywood and Rolf Rannacher. Finite element approximation of the nonstationary navier–stokes problem. i. regularity of solutions and second-order error estimates for spatial discretization.SIAM Journal on Numerical Analysis, 19(2):275–311, 1982

    1982

  21. [21]

    Artificial boundaries and flux and pressure conditions for the incompressible navier–stokes equations.In- ternational Journal for numerical methods in fluids, 22(5):325–352, 1996

    John G Heywood, Rolf Rannacher, and Stefan Turek. Artificial boundaries and flux and pressure conditions for the incompressible navier–stokes equations.In- ternational Journal for numerical methods in fluids, 22(5):325–352, 1996

    1996

  22. [22]

    S. Hussain. A Note on Accurate and Efficient Higher Order Galerkin Time Stepping Schemes for the Nonstationary Stokes Equations.The Open Numerical Methods Journal, 4(1):35–45, 1 2012

    2012

  23. [23]

    Springer Cham, Berlin; Heidelberg, 1

    Volker John.Finite Element Methods for Incompressible Flow Problems, vol- ume 51 ofSpringer series in computational mathematics. Springer Cham, Berlin; Heidelberg, 1. edition edition, 2016

    2016

  24. [24]

    Higher-order finite element discretizations in a benchmark problem for incompressible flows.International Journal for Numerical Methods in Fluids, 37(8):885–903, 2001

    Volker John and Gunar Matthies. Higher-order finite element discretizations in a benchmark problem for incompressible flows.International Journal for Numerical Methods in Fluids, 37(8):885–903, 2001

    2001

  25. [25]

    SIAM, 1995

    Carl T Kelley.Iterative methods for linear and nonlinear equations. SIAM, 1995

    1995

  26. [26]

    Adaptive space–time finite element meth- ods for parabolic optimal control problems.Journal of Numerical Mathematics, 30(4):247–266, 11 2021

    Ulrich Langer and Andreas Schafelner. Adaptive space–time finite element meth- ods for parabolic optimal control problems.Journal of Numerical Mathematics, 30(4):247–266, 11 2021

    2021

  27. [27]

    De Gruyter, Berlin, Boston, 2019

    Ulrich Langer and Olaf Steinbach, editors.Space-Time Methods. De Gruyter, Berlin, Boston, 2019

    2019

  28. [28]

    Anhpmultigrid approach for tensor-product space-time finite element discretizations of the stokes equations

    Nils Margenberg, Markus Bause, and Peter Munch. Anhpmultigrid approach for tensor-product space-time finite element discretizations of the stokes equations. SIAM Journal on Scientific Computing, 47(6):B1503–B1529, 2025

    2025

  29. [29]

    Matthies and L

    G. Matthies and L. Tobiska. The inf-sup condition for the mapped q k- p k-1 disc element in arbitrary space dimensions.Computing, 69(2):119–139, 2002

    2002

  30. [30]

    McCormick.Multigrid Methods

    Stephen F. McCormick.Multigrid Methods. Society for Industrial and Applied Mathematics, 1987

    1987

  31. [31]

    Time-parallel iterative solvers for parabolic evolution equations.SIAM Journal on Scientific Computing, 41(1):C28–C51, 2019

    Martin Neum¨ uller and Iain Smears. Time-parallel iterative solvers for parabolic evolution equations.SIAM Journal on Scientific Computing, 41(1):C28–C51, 2019

    2019

  32. [32]

    ¨Uber ein variationsprinzip zur l¨ osung von dirichlet-problemen bei verwendung von teilr¨ aumen, die keinen randbedingungen unterworfen sind

    Joachim Nitsche. ¨Uber ein variationsprinzip zur l¨ osung von dirichlet-problemen bei verwendung von teilr¨ aumen, die keinen randbedingungen unterworfen sind. InAbhandlungen aus dem mathematischen Seminar der Universit¨ at Hamburg, volume 36, pages 9–15. Springer, 1971

    1971

  33. [33]

    Springer, 2017

    Thomas Richter.Fluid-structure interactions: models, analysis and finite ele- ments, volume 118. Springer, 2017. 27

    2017

  34. [34]

    Tensor- Product Space-Time Goal-Oriented Error Control and Adaptivity With Partition-of-Unity Dual-Weighted Residuals for Nonstationary Flow Problems

    Julian Roth, Jan Philipp Thiele, Uwe K¨ ocher, and Thomas Wick. Tensor- Product Space-Time Goal-Oriented Error Control and Adaptivity With Partition-of-Unity Dual-Weighted Residuals for Nonstationary Flow Problems. Comput. Methods Appl. Math., 24(1):185–214, 2024

    2024

  35. [35]

    SIAM, 2003

    Yousef Saad.Iterative methods for sparse linear systems. SIAM, 2003

    2003

  36. [36]

    Benchmark computations of laminar flow around a cylinder

    Michael Sch¨ afer, Stefan Turek, Franz Durst, Egon Krause, and Rolf Rannacher. Benchmark computations of laminar flow around a cylinder. InFlow simulation with high-performance computers II: DFG priority research programme results 1993–1995, pages 547–566. Springer, 1996

    1993

  37. [37]

    Schmich and B

    M. Schmich and B. Vexler. Adaptivity with dynamic meshes for space-time finite element discretizations of parabolic equations.SIAM J. Sci. Comput., 30(1):369–393, 2008

    2008

  38. [38]

    Olaf Steinbach and Huidong Yang.An Algebraic Multigrid Method for an Adap- tive Space–Time Finite Element Discretization. 1 2018

    2018

  39. [39]

    Space-time finite element methods for parabolic evolution equations: discretization, a posteriori error estimation, adap- tivity and solution

    Olaf Steinbach and Huidong Yang. Space-time finite element methods for parabolic evolution equations: discretization, a posteriori error estimation, adap- tivity and solution. In H. Taghavy et al., editors,Space-time methods: Applica- tions to partial differential equations, De Gruyter, pages 207–248. De Gruyter, Berlin, 2019

    2019

  40. [40]

    North-Holland, Amsterdam, 1984

    Roger Temam.Navier–Stokes Equations. North-Holland, Amsterdam, 1984

    1984

  41. [41]

    Springer-Verlag, Berlin, Hei- delberg, 2006

    Vidar Thom´ ee.Galerkin Finite Element Methods for Parabolic Problems (Springer Series in Computational Mathematics). Springer-Verlag, Berlin, Hei- delberg, 2006

    2006

  42. [42]

    Block-implicit multigrid solution of navier-stokes equations in primitive variables.Journal of Computational Physics, 65(1):138–158, 1986

    S Pratap Vanka. Block-implicit multigrid solution of navier-stokes equations in primitive variables.Journal of Computational Physics, 65(1):138–158, 1986

    1986

  43. [43]

    Optimization with nonstationary, nonlin- ear monolithic fluid-structure interaction.International Journal for Numerical Methods in Engineering, 122(19):5430–5449, 2021

    Thomas Wick and Winnifried Wollner. Optimization with nonstationary, nonlin- ear monolithic fluid-structure interaction.International Journal for Numerical Methods in Engineering, 122(19):5430–5449, 2021

    2021

  44. [44]

    A space-time discontinuous galerkin discretization for the linear transport equation.Computers & Mathematics with Applications, 152:294– 307, 2023

    Christian Wieners. A space-time discontinuous galerkin discretization for the linear transport equation.Computers & Mathematics with Applications, 152:294– 307, 2023. 28

    2023