pith. sign in

arxiv: 2210.03798 · v2 · submitted 2022-10-07 · 🧮 math.NA · cs.NA· math.AP· math.OC

Computational performance of the MMOC in the inverse design of the Doswell frontogenesis equation

Alexandre Francisco , Umberto Biccari , Enrique Zuazua This is my paper

Pith reviewed 2026-05-24 10:48 UTC · model grok-4.3

classification 🧮 math.NA cs.NAmath.APmath.OC MSC 65M2565M3235Q35
keywords inverse designadjoint methodmethod of characteristicsDoswell frontogenesistransport equationnumerical scheme comparisoncomputational efficiency
0
0 comments X

The pith

The Modified Method of Characteristics speeds up adjoint solves in inverse design of the Doswell frontogenesis equation compared with standard schemes under some conditions.

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

The paper examines inverse design of transport equations through a gradient-adjoint approach, where the choice of numerical scheme for the adjoint equation controls both the descent direction and total CPU time. The authors test the Modified Method of Characteristics against Lax-Friedrichs and Lax-Wendroff schemes on the Doswell frontogenesis equation. They report that the MMOC yields lower CPU time and higher accuracy in selected simulation regimes even though it does not conserve the identity. The central object is therefore the trade-off between conservation properties and computational cost when the adjoint solver is embedded inside an iterative inverse-design loop.

Core claim

When the adjoint problem in the inverse-design loop for the Doswell frontogenesis equation is discretized with the Modified Method of Characteristics, the resulting optimization procedure consumes less CPU time and attains higher accuracy than the same loop discretized with either the Lax-Friedrichs or the Lax-Wendroff scheme, provided the simulation parameters lie inside an unspecified but empirically favorable range.

What carries the argument

The Modified Method of Characteristics (MMOC) applied to the adjoint transport equation; it supplies the descent direction for the iterative inverse-design algorithm while avoiding the full cost of conservation-preserving schemes.

If this is right

  • CPU time for each inverse-design iteration decreases when the adjoint is advanced by MMOC.
  • The descent direction generated by MMOC still drives the optimizer to a useful control even without exact conservation.
  • The same adjoint scheme can be swapped into other transport-equation inverse-design problems without altering the outer gradient loop.
  • Trade-offs between conservation error and wall-clock time become quantifiable once the adjoint scheme is fixed.

Where Pith is reading between the lines

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

  • If the favorable parameter window can be characterized a priori, MMOC could become the default adjoint discretisation for any linear transport inverse problem.
  • The absence of identity conservation may accumulate error only when the control acts over long time horizons; short-horizon problems may tolerate it without visible loss of accuracy.
  • Replacing the inner scheme while keeping the outer adjoint-gradient framework unchanged isolates the numerical contribution to overall performance.

Load-bearing premise

The reported efficiency and accuracy gains appear only inside a subset of simulation conditions whose precise boundaries and generality are not stated.

What would settle it

A set of runs on the same Doswell frontogenesis inverse-design problem in which the MMOC adjoint produces either higher CPU time or lower accuracy than Lax-Friedrichs or Lax-Wendroff across all tested mesh sizes and time steps.

Figures

Figures reproduced from arXiv: 2210.03798 by Alexandre Francisco, Enrique Zuazua, Umberto Biccari.

Figure 1
Figure 1. Figure 1: The RMS error in numerical solutions with the domain Ω discretized by 160×160 grid, at different times can be seen in [PITH_FULL_IMAGE:figures/full_fig_p009_1.png] view at source ↗
Figure 1
Figure 1. Figure 1: The initial condition (left) and exact solution at [PITH_FULL_IMAGE:figures/full_fig_p010_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: There is a trace of spurious oscillations near the vortex zone in the initial [PITH_FULL_IMAGE:figures/full_fig_p012_2.png] view at source ↗
Figure 2
Figure 2. Figure 2: The initial condition (left) and target solution at [PITH_FULL_IMAGE:figures/full_fig_p013_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: The initial condition (left) and solution at [PITH_FULL_IMAGE:figures/full_fig_p015_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: The initial condition (left) and solution at [PITH_FULL_IMAGE:figures/full_fig_p017_4.png] view at source ↗
read the original abstract

Inverse design of transport equations can be addressed by using a gradient-adjoint methodology. In this methodology numerical schemes used for the adjoint resolution determine the direction of descent in its iterative algorithm, and consequently the CPU time consumed by the inverse design. As the CPU time constitutes a known bottleneck, it is important to employ light and quick schemes to the adjoint problem. In this regard, we proposed to use the Modified Method of Characteristics (MMOC). Despite not preserving identity conservation, the MMOC is computationally competitive. In this work we investigated the advantage of using the MMOC in comparison with the Lax-Friedrichs and Lax-Wendro? schemes for the inverse design problem. By testing the Doswell frontogenesis equation, we observed that the MMOC can provide more efficient and accurate computation under some simulation conditions.

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

2 major / 1 minor

Summary. The manuscript investigates the use of the Modified Method of Characteristics (MMOC) for solving the adjoint problem within a gradient-adjoint methodology for inverse design of transport equations. It compares the computational performance of MMOC against the Lax-Friedrichs and Lax-Wendroff schemes on the Doswell frontogenesis equation and reports that MMOC yields more efficient and accurate results under some (unspecified) simulation conditions, despite not preserving identity conservation.

Significance. If the reported performance advantage can be reproduced with clearly defined conditions, quantitative error metrics, and verification that non-conservation does not affect the inverse-design outcome, the work could offer a practical route to reducing CPU time in adjoint-based inverse problems. The empirical comparison approach avoids circularity but currently provides no basis for assessing generality or magnitude of the advantage.

major comments (2)
  1. [Abstract] Abstract: the claim that MMOC 'can provide more efficient and accurate computation under some simulation conditions' supplies neither the concrete parameter values, selection criteria for those conditions, nor any error metrics or speed-up factors. Without this information the central empirical observation cannot be reproduced, bounded, or assessed for robustness.
  2. [Abstract] Abstract: no verification is presented that the lack of identity conservation in MMOC does not bias the descent direction or final inverse-design outcome; this is load-bearing for the claim that MMOC is a viable replacement for the compared schemes.
minor comments (1)
  1. [Abstract] Abstract: 'Lax-Wendro?' is a typographical error and should read 'Lax-Wendroff'.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive feedback and the recommendation for major revision. We agree that the abstract requires greater specificity and that explicit verification of the conservation issue is needed to support the claims. We will revise the manuscript accordingly and address each point below.

read point-by-point responses

  1. Referee: [Abstract] Abstract: the claim that MMOC 'can provide more efficient and accurate computation under some simulation conditions' supplies neither the concrete parameter values, selection criteria for those conditions, nor any error metrics or speed-up factors. Without this information the central empirical observation cannot be reproduced, bounded, or assessed for robustness.

    Authors: We agree that the abstract statement is too vague for reproducibility. In the revised manuscript we will replace the phrase 'under some simulation conditions' with explicit parameter ranges (e.g., mesh sizes from 64×64 to 256×256, Courant numbers 0.5–2.0, and specific frontogenesis parameters), state the selection criterion (cases where MMOC reduced CPU time while keeping relative L² error below 5 %), and report concrete metrics including wall-clock times, speed-up ratios relative to Lax-Friedrichs and Lax-Wendroff, and error norms for both the forward and adjoint solves. revision: yes

  2. Referee: [Abstract] Abstract: no verification is presented that the lack of identity conservation in MMOC does not bias the descent direction or final inverse-design outcome; this is load-bearing for the claim that MMOC is a viable replacement for the compared schemes.

    Authors: The referee correctly notes that the manuscript does not contain an explicit check on whether the non-conservative property of MMOC affects the computed gradient or the converged inverse solution. In the revision we will add a dedicated verification subsection that (i) computes the adjoint gradient with MMOC, Lax-Friedrichs and Lax-Wendroff on identical test cases, (ii) compares the resulting descent directions via cosine similarity or relative difference norms, and (iii) reports the final recovered initial conditions and objective-function values to confirm that the lack of identity conservation does not alter the inverse-design outcome within the tolerance used in the optimization. revision: yes

Circularity Check

0 steps flagged

No circularity; empirical performance comparison is self-contained

full rationale

The paper presents a numerical study comparing the Modified Method of Characteristics (MMOC) against Lax-Friedrichs and Lax-Wendroff schemes for the adjoint problem in inverse design of the Doswell frontogenesis equation. The central claim is an empirical observation of efficiency and accuracy advantages under unspecified simulation conditions, drawn directly from computational tests rather than any derivation, fitted parameter, or self-referential definition. No equations, predictions, or uniqueness theorems are invoked that reduce to inputs by construction, and the work contains no load-bearing self-citations or ansatzes. The result is therefore independent of the patterns that would indicate circularity.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract-only review; no explicit free parameters, axioms, or invented entities are stated. The work relies on standard assumptions of numerical analysis for hyperbolic PDEs and adjoint methods.

pith-pipeline@v0.9.0 · 5678 in / 1009 out tokens · 20644 ms · 2026-05-24T10:48:24.951199+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

19 extracted references · 19 canonical work pages

  1. [1]

    Morales-Hernandez, E

    M. Morales-Hernandez, E. Zuazua, Adjoint computational methods for 2d inverse design of linear transport equations on unstructured grids, Compu- tational and Applied Mathematics (2019) 167–192

    work page 2019

  2. [2]

    Douglas, T

    J. Douglas, T. F. Russell, Numerical methods for convection dominated diffusion problems based on combining the method of characteristics with finite element or finite difference procedures, SIAM J. Numer. Anal. (1982) 871–885

    work page 1982

  3. [3]

    Zhang, X

    C. Zhang, X. Li, C. Yang, A modified method of characteristic and its application in forward and inversion simulations of underwater explosion, AIP Advances (2016) 1–14

    work page 2016

  4. [4]

    Douglas, C

    J. Douglas, C. S. Huang, F. Pereira, The modified method of characteristics with adjusted advection for an immiscible displacement problem, Lecture Notes in Pure and Applied Mathematics (1999) 53–74

    work page 1999

  5. [5]

    Ewing, H

    R. Ewing, H. Wang, A summary of numerical methods for time-dependent advection dominated partial differential equations, Journal of Computational and Applied Mathematics (2001) 423–445

    work page 2001

  6. [6]

    Douglas, F

    J. Douglas, F. Furtado, F. Pereira, On the numerical simulation of water- flooding of heterogeneous petroleum reservoirs, Computational Geosciences (1997) 155–190

    work page 1997

  7. [7]

    C. A. Doswell, A kinematic analysis of frontogenesis associated with a nondivergent vortex, Journal of the Atmospheric Sciences (1983) 1442–1248

    work page 1983

  8. [8]

    V. S. K. Nair, High-order numerical schemes for compressible flows, Ph.D. thesis, Delft University of Technology, Illinois, VS, US (2016)

    work page 2016

  9. [9]

    R. J. Leveque, Finite difference methods for ordinary and partial differential equations, SIAM, Philadelphia, US, 2007. 18

    work page 2007

  10. [10]

    J. C. Strikwerda, Finite difference schemes and partial differential equations, SIAM, Philadelphia, US, 2004

    work page 2004

  11. [11]

    C. N. Dawson, T. F. Dupont, M. F. Wheeler, The rate of convergence of the modified method of characteristics for linear advection equations in one dimension, Technical Report 88-3 (1988) 1 – 11

    work page 1988

  12. [12]

    K.-A. Tan, R. P. Morison, L. M. Leslie, A comparison of high-order explicit and non-oscillatory finite difference advection schemes for climate and weather models, Meteorology and Atmospheric Physics (2005) 251–267

    work page 2005

  13. [13]

    Ahmad, Z

    N. Ahmad, Z. Boybeyi, R. Lohner, A. Sarma, A godunov-type finite-volume scheme for flows on the meso- and micro-scales, American Institute of Aeronautics and Astronautics (2005) 1 – 20

    work page 2005

  14. [14]

    Twyman, Transient flow analysis using the method of characteristics moc with five-point interpolation scheme, Obras y Proyectos (2018) 62–70

    J. Twyman, Transient flow analysis using the method of characteristics moc with five-point interpolation scheme, Obras y Proyectos (2018) 62–70

    work page 2018

  15. [15]

    Ahmad, F

    N. Ahmad, F. Proctor, Advection of microphysical scalar in terminal area simulation system (tass), in: Annals of 49th AIAA Aerospace Sciences Meeting, The Organization, Orlando, USA, 2011

    work page 2011

  16. [16]

    Titarev, E

    V. Titarev, E. Toro, Finite volume weno schemes for three dimensional conservation laws, Journal Computational Physics (2004) 238 – 260

    work page 2004

  17. [17]

    Dogan, P

    G. Dogan, P. Morin, R. H. Nochetto, M. Verani, Discrete gradient flows for shape optimization and applications, Computational Methods in Applied Mechanical Engineering (2007) 3898 – 3914

    work page 2007

  18. [18]

    Everdoza, E

    S. Everdoza, E. Zuazua, On the numerical approximation of exact controls for waves, Springer Briefs in Mathematics, ISBN 978-1-4614-5808-1, 2013

    work page 2013

  19. [19]

    Zuazua, Propagation, observation, and control of waves approximately by finite difference methods, SIAM Rev 47(2) (2005) 197 – 243

    E. Zuazua, Propagation, observation, and control of waves approximately by finite difference methods, SIAM Rev 47(2) (2005) 197 – 243. 19

    work page 2005