Analysis of the multi-dimensional semi-discrete Active Flux method using the Fourier transform

Christian Klingenberg; Janina Kern; Lisa Lechner; Wasilij Barsukow

REVIEW 3 major objections 4 minor 2 cited by

Analysis of the multi-dimensional semi-discrete Active Flux method using the Fourier transform

T0 review · 3 major / 4 minor · reviewed 2026-08-11 · deepseek-v4-flash

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

Pith's one-line read For linear acoustics, semi-discrete Active Flux preserves all steady states in 2D and 3D.

desk verdict Solid extension of the Active Flux stationarity-preservation program to the semi-discrete variant; the central claim rests on an unshipped Mathematica rank computation, so the fix is to make that check reproducible. read the letter →

arxiv 2412.03477 v2 pith:HK4GWEWH submitted 2024-12-04 math.NA cs.NA

Wasilij Barsukow , Janina Kern , Christian Klingenberg , Lisa Lechner This is my paper · ORCID

classification math.NAcs.NA MSC 65M2065M7065M0835E15

keywords ActiveFluxstationaritypreservinglinearacousticsFouriertransformsemi-discretemethodnumericaldiffusionCartesiangridsupwindscheme

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

This paper proves that the semi-discrete (generalized) Active Flux method, which replaces the exact evolution operator of classical Active Flux with finite-difference approximations of spatial derivatives, retains stationarity preservation for linear acoustics on Cartesian grids. In two dimensions the evolution matrix has a one-dimensional kernel, and in three dimensions a five-dimensional kernel, whose elements correspond exactly to divergence-free velocity reconstructions with constant pressure. These discrete stationary states coincide with those of the classical Active Flux method in 2D, so the easier-to-implement semi-discrete version loses none of the structure preservation. The analysis also shows that the upwind Jacobian splitting is essential: a Rusanov-type splitting destroys stationarity preservation entirely. Numerical experiments confirm third-order convergence for smooth waves and the long-time preservation of vortices and vortex rings.

What carries the argument

The central object is the evolution matrix $E(k)$ obtained by applying the discrete Fourier transform to the semi-discrete update equations; its kernel is the set of numerical stationary states. The argument that kernel vectors correspond exactly to divergence-free reconstructions rests on the $P^{2,2}$-projected divergence of Definition 5.2 and the unisolvence of the broken space $D^{2,2}_{\mathrm{br}}$ with respect to the eight pointwise degrees of freedom (Lemma 5.1), which lets the paper pass from eight (or 26) pointwise divergence-vanishing conditions to the vanishing of the divergence polynomial itself. The upwind Jacobian splitting matters because it adds the requirement that normal derivatives be continuous at stationary states, and the resulting count of constraints leaves exactly one free parameter in 2D (and five in 3D).

What would settle it

Compute exactly the rank of the 12×12 evolution matrix $E(k)$ of Appendix C (and the 32×32 matrix for 3D) over the field of rational functions in $t_x,t_y,t_z$, for example by evaluating all minors symbolically at generic integer values of the translation factors; the theorem holds only if the generic rank is 11 in 2D and 27 in 3D, and any single wave vector with nullity greater than 1 (2D) or 5 (3D) would disprove stationarity preservation.

Watch

Extended reading notes

Core claim

The central discovery is that the semi-discrete Active Flux method, defined by updating cell averages through exact flux quadrature and point values through one-sided finite differences obtained from a globally continuous biparabolic reconstruction, is stationarity preserving for the linear acoustic equations in two and three spatial dimensions. Concretely, the discrete Fourier transform reduces the method to an evolution matrix $E(k)$, and the paper proves (Theorems 5.3 and 5.4) that for general wave vectors the kernel of $E(k)$ is one-dimensional in 2D and five-dimensional in 3D, consisting precisely of the Fourier modes whose reconstruction has vanishing divergence and constant pressure. The unique kernel vector in 2D is given explicitly, and the 3D basis is provided in an appendix; in both cases the kernels match the analytic stationary states of the acoustic system (the divergence-free velocity fields with zero pressure) plus, in 3D, higher-order terms. The proof uses a unisolvence property of the broken divergence space $D^{2,2}_{\mathrm{br}}$ to show that vanishing of the divergence at all point values forces the polynomial divergence to vanish identically, together with continuity of normal derivatives that the upwind splitting enforces at stationarity.

Load-bearing premise

The assertion that the evolution matrix has no stationary modes beyond the ones listed is verified by a computer algebra system rather than by a hand-written proof, so the kernel-dimension equality stands or falls with that computation.

Editorial extensions

If this is right

Semi-discrete Active Flux can be applied to nonlinear systems of conservation laws without sacrificing correct representation of steady states of the linearized equations, since the structure preservation is inherited from the reconstruction space, not from the exact evolution operator.
The discrete stationary states in 2D are identical to those of classical Active Flux, so vortices and other divergence-free flows are preserved by both schemes; the semi-discrete method can serve as a drop-in replacement that is easier to extend to nonlinear problems.
The upwind splitting is the structure-preserving choice; switching to a Rusanov-type Jacobian splitting (39) produces a method with no non-trivial stationary states, as $\det E(k) \neq 0$ for generic $k$, so the choice of numerical diffusion terms controls stationarity preservation.
The maximum stable time step for the semi-discrete method on linear acoustics is roughly $\Delta t \approx 0.28$ (with $\Delta x = \Delta y = 1$), about half that of the classical method, so third-order accuracy is retained at the cost of more steps per unit time.
In 3D the 5-dimensional kernel includes the two analytic non-trivial stationary modes plus three higher-order modes, meaning the method captures all stationary states of the PDE at the discrete level.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

The same Fourier-kernel analysis could be applied to other linear symmetrizable hyperbolic systems (e.g., linearized Euler or Maxwell equations) to characterize which splittings and reconstruction spaces are stationarity preserving; the general criterion of Definition 5.1 is system-agnostic.
Because the kernel dimension is determined solely by the reconstruction space and the upwind constraints, replacing the biparabolic reconstruction with higher-degree polynomials would likely enlarge the set of preserved stationary states; testing $P^{3,3}$ reconstructions would show whether the matching of classical and semi-discrete stationary states persists.
A practical implication for code designers: the non-physical eigenmode that limits the CFL number is a separate degree-of-freedom mode (the point values vs. averages); damping or filtering that mode could push the stability limit closer to the classical method's CFL.
The Mathematica rank check could be replaced by a constructive symbolic proof (e.g., computing the Smith normal form of $E(k)$ over the polynomial ring), which would remove the only computational step in the proof and make the theorem fully analytic.

Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

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

Desk Editor's Note

read the letter

The new thing here is the proof that the semi-discrete (generalized) Active Flux method preserves stationary states for linear acoustics in 2D and 3D, with explicit kernel bases and a diffusion/stability analysis. The sufficiency direction is done by hand and is clean: stationarity of point values plus continuity of normal derivatives forces the divergence of the reconstruction to vanish, and the average update becomes redundant. The numerical experiments back this up — the vortex test and the well-prepared Fourier modes decay to machine zero, and the 3D vortex ring behaves as advertised. There is no parameter fitting and no circularity; the reliance on earlier papers in the same program is normal and those results are properly cited.

The soft spot is real but narrow. The upper bound in Definition 5.1 — that the kernel has no extra elements beyond the explicitly constructed basis — is asserted on the basis of a Mathematica check that is not shipped. In Theorem 5.3 the paper says “verified using mathematica due to the excessive length of computations,” and Theorem 5.4 says the same. For a theorem whose force is the sharp equality of kernel dimensions, that is the load-bearing step, and an unreproduced computer algebra check is not the same as a proof. The stress-test worry about denominators vanishing on subvarieties (tx = 1, ty = 1) is also legitimate: the displayed basis vectors have poles there, and the paper does not explicitly analyze those wave vectors. I suspect the statements are still true by continuity or by a separate rank calculation, and the lower bound is constructive, so this is not a fatal flaw. But it is exactly the kind of gap a referee should ask to close: ship the notebook, add a short symbolic rank argument, or state the subvariety cases explicitly.

The stability analysis in Section 6 is more heuristic — eigenvalue plots and a CFL estimate — but that is secondary and clearly labeled. The 3D convergence study is honest about the coarse meshes and observed order around 2.5.

Who gets value from this: anyone working on Active Flux, structure-preserving finite volume/finite difference methods, or low-Mach behavior. It deserves a serious referee. My recommendation is to send it to review, with the explicit request that the authors provide the missing computational verification in reproducible form before acceptance.

Referee Report

3 major / 4 minor

Summary. The paper analyzes the semi-discrete (generalized) Active Flux method for the linear acoustic equations on Cartesian grids in two and three spatial dimensions. The method updates cell averages by flux quadrature and point values by finite-difference approximations obtained from a globally continuous biparabolic reconstruction. Using the discrete Fourier transform, the authors derive the evolution matrix E(k) and characterize its kernel. Their central theorems (Theorems 5.3 and 5.4) state that the upwind semi-discrete Active Flux method is stationarity preserving: the numerical stationary states consist exactly of divergence-free velocity reconstructions with constant pressure, with kernel dimensions 1 in 2D and 5 in 3D, matching the non-trivial stationary states of the PDE. The paper also provides an analysis of numerical diffusion and stability and presents numerical experiments for vortex, well-prepared mode, and traveling-wave tests.

Significance. If the results are correct, this is a valuable contribution: it shows that the semi-discrete Active Flux method preserves the same stationary states as the classical fully discrete Active Flux method, extending a desirable structure-preservation property to a method-of-lines formulation that is more easily applied to nonlinear problems. The Fourier-based framework is clean and provides explicit kernel vectors with no fitted parameters. The authors give constructive lower bounds for the kernel, show that well-prepared modes remain stationary to round-off, and demonstrate that the vortex test drives the discrete divergence relations to machine zero. These numerical experiments lend credible support to the main claim. The main weakness is that the sharpness of the kernel dimension (the upper bound in Definition 5.1) is established only by an unreproduced computer algebra computation, which is load-bearing for the stationarity-preservation claim.

major comments (3)

[Section 5.3.1, Theorem 5.3 and Theorem 5.4] The upper bound in Definition 5.1, namely dim ker E(k) ≤ S·Ndof, is the load-bearing part of the stationarity-preservation claim, but the only evidence offered is the statements 'verified using mathematica due to the excessive length of computations' (end of the proof of Theorem 5.3) and 'we confirm using mathematica' (proof of Theorem 5.4). No notebook, script, or algebraic certificate is provided, so the claim is not independently verifiable from the manuscript. Since the lower bound alone shows only that the exhibited modes are stationary and does not rule out additional spurious kernel elements, this gap directly affects the central theorem. Please provide a reproducible symbolic computation (e.g., a notebook or a rank computation over C(tx,ty,tz) with a certificate such as a maximal minor) or a hand proof that the homogeneous system has full rank on the relevant complement.
[Theorems 5.3 and 5.4 and Appendix D] The displayed kernel vectors (93) and Q1,...,Q5 in Appendix D have entries with denominators that vanish on subvarieties of the torus, e.g., tx=1, ty=1, tz=1, 1+tx=0, and 1+ty=0. The paper does not analyze whether the exact equality dim ker E(k)=S holds on these subvarieties, either by continuity, by limits, or by separate computations. Since Definition 5.1 involves min_k over all k, including these special wave vectors, the 'for general k' statement is insufficient as written. Please specify the exceptional set of wave vectors (if any) and verify the dimension claim there or explain why the limit argument applies.
[Section 4.3 and Definition 5.1] The paper defines stationarity preservation through a comparison of kernel dimensions, but the numerical tests in Section 7 (vortex decay to machine zero and well-prepared modes) can only confirm that the constructed vectors lie in the kernel; they cannot certify the absence of extra kernel elements. This does not invalidate the experiments, but it means the numerical evidence does not fill the gap left by the computer algebra step. A statement acknowledging this explicitly would help the reader calibrate the role of the Mathematica verification.

minor comments (4)

[Equation (96)] There is an unbalanced parenthesis in the formula for det E(k): '... + 2c2(∆x2 − ∆x∆y + ∆y2)))' appears to have an extra closing parenthesis; please correct.
[Section 6.2] The stability analysis reports that the instability appears 'between ∆t = 0.28 and ∆t = 0.3' for the 2D acoustics case. Since ∆x = ∆y = 1 in that figure, it may be helpful to state explicitly that ∆t equals the CFL number there, and to clarify whether this bound is for the given fixed wave-number parametrization.
[Eq. (114) and Section 7.2.2] In the expression (2kz(ky − i), −2kxkz, 2ikx, 0), the symbol i is used for the imaginary unit, but this is not defined at that point; consider writing 'i' as 'i' with a comment that this is the imaginary unit, or use ı to avoid confusion with an index.
[General notation] The notation '1 + tx(4 + tx)' is unambiguous but could be typeset more clearly as '1 + 4tx + tx^2' for readability; similarly for other kernel vectors. This is only a presentation issue.

Circularity Check

0 steps flagged · score 1.0 of 10

No significant circularity: the stationarity-preservation analysis is derived from the explicit evolution matrix, with only a Mathematica rank check delegated rather than proven by hand.

full rationale

The central claims (Theorems 5.3 and 5.4) are obtained by writing the semi-discrete Active Flux scheme as d/dt qhat + E(k) qhat = 0 and analyzing the kernel of E(k). The lower bound on dim ker E is proved constructively: the paper exhibits explicit kernel vectors, proves that vanishing projected divergence and continuous normal derivatives imply stationarity, and invokes Lemma 5.1 (unisolvence of the broken divergence space) with a proof by explicit computation. No parameter is fitted to the target result, and the scheme constants are fixed by the reconstruction and quadrature rules before the analysis begins. The upper bound on dim ker E, i.e. the absence of additional kernel elements, is delegated to Mathematica ('That there are no other linearly independent elements in the kernel has been verified using mathematica due to the excessive length of computations'). This is a computational verification rather than a self-citation or a fitted-input-as-prediction, so it is not circular; it is a reproducibility and rigor caveat. The paper's self-citations to [Bar19], [BHKR19], and [AB23a] supply definitions, the classical Active Flux method, and the semi-discrete construction under study, but the stationarity-preservation proof does not assume the target theorem from those works. Corollary 5.2 compares the derived stationary states with those of classical Active Flux, and the numerical experiments in Section 7 are illustrative rather than load-bearing for the kernel-dimension claim. Accordingly, no step reduces to its own input by construction, and the residual concern is limited to the unshipped Mathematica verification of the exact kernel dimension.

Assumptions & free parameters 0 free parameters · 6 assumptions · 0 invented entities

No free parameters are fitted; the only numerical parameters are grid spacings, wave vectors, and CFL numbers, chosen for tests. The analysis relies on two domain assumptions: linear acoustics with constant coefficients on an infinite or periodic Cartesian grid, and the upwind Jacobian splitting. The paper also depends on CAS-verified kernel computations for the necessity parts. No new physical entities are introduced.

assumptions (6)

standard math The acoustic system is hyperbolic, so J·k is diagonalizable for all k.
Used throughout Section 2 to define Fourier modes and stationary states.
domain assumption Infinite or periodic equidistant Cartesian grid, so the discrete Fourier transform diagonalizes the scheme.
Section 4 assumes translation-invariant lattices with spacings Δx, Δy, Δz.
domain assumption The method is linear with constant coefficients, so superposition of Fourier modes is valid.
Analysis restricted to linear acoustics (4).
standard math Unisolvence of the P^{2,2} reconstruction space and D^{2,2}_br with respect to the chosen degrees of freedom (Lemma 5.1).
Proven by explicit computation in the paper; used in Theorems 5.2 and 5.3.
domain assumption Upwind Jacobian splitting (36) is used for the point value updates; Rusanov-type splitting is not stationarity preserving.
Section 3.3 and Section 5.3.1; the alternative (39) is shown to fail.
domain assumption Stationarity preservation is defined via the kernel dimension inequality of Definition 5.1.
The paper adopts the criterion from Bar19 to compare discrete and continuous stationary states.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Analysis of the multi-dimensional semi-discrete Active Flux method using the Fourier transform." pith.science (2026). https://pith.science/paper/HK4GWEWH

@misc{pith2026241203477,
  author       = {Pith},
  title        = {Pith review of: Analysis of the multi-dimensional semi-discrete Active Flux method using the Fourier transform},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/HK4GWEWH}},
  note         = {Machine review of arXiv:2412.03477}
}

read the original abstract

The degrees of freedom of Active Flux are cell averages and point values along the cell boundaries. These latter are shared between neighbouring cells, which gives rise to a globally continuous reconstruction. The semi-discrete Active Flux method uses its degrees of freedom to obtain Finite Difference approxi\-mations to the spatial derivatives which are used in the point value update. The averages are updated using a quadrature of the flux and making use of the point values as quadrature points. The integration in time employs standard Runge-Kutta methods. We show that this generalization of the Active Flux method in two and three spatial dimensions is stationarity preserving for linear acoustics on Cartesian grids, and present an analysis of numerical diffusion and stability.

Figures

Figures reproduced from arXiv: 2412.03477 by the authors.

**Figure 2.** Shape functions B7 ( top left), B8 ( top right), B0 ( bottom left) in two spatial dimensions and a reconstruction devised by using these shape functions ( bottom right). The same applies to the formulas Dy [PITH_FULL_IMAGE:figures/full_fig_p009_2.png] view at source ↗

**Figure 3.** Illustration of the coefficients occurring in the difference formulas for node values and [PITH_FULL_IMAGE:figures/full_fig_p011_3.png] view at source ↗

Figures from the paper (16 more)

**Figure 4.** Figure 4: Numerical diffusion analysis of Active Flux for 1-d linear advection. The absolute value [PITH_FULL_IMAGE:figures/full_fig_p024_4.png]

**Figure 5.** Figure 5: Numerical diffusion analysis of Active Flux for 2-d linear acoustics. The absolute value [PITH_FULL_IMAGE:figures/full_fig_p026_5.png]

**Figure 6.** Figure 6: The setup of a spherical Gaussian in the pressure. [PITH_FULL_IMAGE:figures/full_fig_p027_6.png]

**Figure 7.** Figure 7: The error of the numerical solution for the spherical Gaussian wave at time [PITH_FULL_IMAGE:figures/full_fig_p027_7.png]

**Figure 8.** Figure 8: Left: Initial data of the stationary vortex test. The magnitude of the velocity is colorcoded. Right: Error of the numerical solution as a function of time. One observes the stationarization of the setup (the velocity components are on top of each other) [PITH_FULL_…

**Figure 9.** Figure 9: Stationary vortex test case. Decay of the 7 discrete divergences [PITH_FULL_IMAGE:figures/full_fig_p028_9.png]

**Figure 10.** Figure 10: Rusanov-type Jacobian split (39) does not lead to a stationarity preserving method. Left: Solution of the stationary vortex test at t = 1000, the magnitude of the velocity is colorcoded. The vortex has become two overlapping shear flows. Right: Error of the numerical…

**Figure 11.** Figure 11: Stationary mode for 2-d acoustics. Top. Initial data in u, nodal values are shown. Left: Well-prepared data (107). Right: Not well-prepared data directly from (110). There is no visible difference. Bottom: Error of the numerical solution is shown as a function of time…

**Figure 12.** Figure 12: Convergence study for the spherical Gaussian in the pressure. Computations on finer [PITH_FULL_IMAGE:figures/full_fig_p032_12.png]

**Figure 13.** Figure 13: The L 1 distance between the initial data of a Fourier mode in the kernel of the evolution matrix E for Active Flux in three spatial dimensions shown as function of time. Left: Non-wellprepared initialization according to (114). The mode is not discretely stationary,…

**Figure 14.** Figure 14: Coordinate setup for the vortex ring test case. [PITH_FULL_IMAGE:figures/full_fig_p033_14.png]

**Figure 15.** Figure 15: View of the vortex ring setup, showing the toroidal boundary of the support of the [PITH_FULL_IMAGE:figures/full_fig_p034_15.png]

**Figure 16.** Figure 16: Vortex ring test case. Top: Cut of the initial data through the plane z = 0 ( left) and y = 0 ( right). Color coded is the magnitude of the velocity, and the arrows (normalized in length) indicate the direction. Bottom: The same for the numerical solution at time t = …

**Figure 17.** Figure 17: Vortex ring test case. L 1 error of the numerical solution on grids of 203 and 353 cells, computed on the nodal point values. One observes the stationarization of the setup, as it converges towards the numerical stationary state. 7.3 Comparison to the fully discrete A…

**Figure 18.** Figure 18: Spherical Riemann problem. Numerical results for the pressure [PITH_FULL_IMAGE:figures/full_fig_p037_18.png]

**Figure 19.** Figure 19: Spherical Riemann problem. Numerical results for the pressure [PITH_FULL_IMAGE:figures/full_fig_p037_19.png]

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 2 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score. Full citation record

An Active Flux method for the Euler equations based on the exact acoustic evolution operator
math.NA 2025-06 conditional novelty 7.0 of 10

A new Active Flux scheme for the 2D Euler equations combines an exactly solved locally linearized acoustic operator with a third-order advective evolution operator inside an additive splitting, using primitive variabl...
Semi-discrete Active Flux as a Petrov-Galerkin method: the case of one-dimensional and Cartesian grids
math.NA 2025-08 unverdicted novelty 6.0 of 10

Semi-discrete Active Flux methods are shown to be Petrov-Galerkin schemes with discontinuous biorthogonal test functions, with explicit constructions for arbitrary order in one dimension and third order on Cartesian grids.

Reference graph

Works this paper leans on

15 extracted references · 4 canonical work pages · cited by 2 Pith papers

[1]

Extensions of A ctive F lux to arbitrary order of accuracy

Remi Abgrall and Wasilij Barsukow. Extensions of A ctive F lux to arbitrary order of accuracy. ESAIM: Mathematical Modelling and Numerical Analysis , 57(2):991--1027, 2023

2023
[2]

A hybrid finite element--finite volume method for conservation laws

R \'e mi Abgrall and Wasilij Barsukow. A hybrid finite element--finite volume method for conservation laws. Applied Mathematics and Computation , 447:127846, 2023

work page 2023
[3]

R \'e mi Abgrall. A combination of R esidual D istribution and the A ctive F lux formulations or a new class of schemes that can combine several writings of the same hyperbolic problem: application to the 1d E uler equations. Communications on Applied Mathematics and Computation , pages 1--33, 2022

2022
[4]

A semi-discrete A ctive F lux method for the E uler equations on C artesian grids

Rémi Abgrall, Wasilij Barsukow, and Christian Klingenberg. A semi-discrete A ctive F lux method for the E uler equations on C artesian grids. accepted 2024 in J. Sci. Comp., arXiv preprint arXiv:2310.00683 , 2024

arXiv 2024
[5]

Stationarity preserving schemes for multi-dimensional linear systems

Wasilij Barsukow. Stationarity preserving schemes for multi-dimensional linear systems. Mathematics of Computation , 88(318):1621--1645, 2019

work page 2019
[6]

Stationarity preservation properties of the active flux scheme on cartesian grids

Wasilij Barsukow. Stationarity preservation properties of the active flux scheme on cartesian grids. Communications on Applied Mathematics and Computation , pages 1--15, 2020

work page 2020
[7]

The active flux scheme for nonlinear problems

Wasilij Barsukow. The active flux scheme for nonlinear problems. Journal of Scientific Computing , 86(1):1--34, 2021

2021
[8]

The active flux scheme on C artesian grids and its low M ach number limit

Wasilij Barsukow, Jonathan Hohm, Christian Klingenberg, and Philip L Roe. The active flux scheme on C artesian grids and its low M ach number limit. Journal of Scientific Computing , 81(1):594--622, 2019

2019

Show all 15 references

[9]

Exact solution and a truly multidimensional G odunov scheme for the acoustic equations

Wasilij Barsukow and Christian Klingenberg. Exact solution and a truly multidimensional G odunov scheme for the acoustic equations. ESAIM: M2AN , 56(1), 2022

2022
[10]

The C artesian G rid A ctive F lux method: L inear stability and bound preserving limiting

Erik Chudzik, Christiane Helzel, and David Kerkmann. The C artesian G rid A ctive F lux method: L inear stability and bound preserving limiting. Applied Mathematics and Computation , 393:125501, 2021

2021
[11]

Multidimensional active flux schemes

Timothy A Eymann and Philip L Roe. Multidimensional active flux schemes. In 21st AIAA computational fluid dynamics conference , 2013

2013
[12]

A new ADER method inspired by the active flux method

Christiane Helzel, David Kerkmann, and Leonardo Scandurra. A new ADER method inspired by the active flux method. Journal of Scientific Computing , 80(3):1463--1497, 2019

2019
[13]

Designing CFD methods for bandwidth—a physical approach

Philip Roe. Designing CFD methods for bandwidth—a physical approach. Computers & Fluids , 214:104774, 2021

2021
[14]

Towards the ultimate conservative difference scheme

Bram van Leer. Towards the ultimate conservative difference scheme. IV . A new approach to numerical convection. Journal of computational physics , 23(3):276--299, 1977

1977
[15]

Linear hybrid-variable methods for advection equations

Xianyi Zeng. Linear hybrid-variable methods for advection equations. Advances in Computational Mathematics , 45(2):929--980, 2019

2019

Pith tools

Reviewed August 11, 2026 · model on record in the stance chip above.

[1] [1]

Extensions of A ctive F lux to arbitrary order of accuracy

Remi Abgrall and Wasilij Barsukow. Extensions of A ctive F lux to arbitrary order of accuracy. ESAIM: Mathematical Modelling and Numerical Analysis , 57(2):991--1027, 2023

2023

[2] [2]

A hybrid finite element--finite volume method for conservation laws

R \'e mi Abgrall and Wasilij Barsukow. A hybrid finite element--finite volume method for conservation laws. Applied Mathematics and Computation , 447:127846, 2023

work page 2023

[3] [3]

R \'e mi Abgrall. A combination of R esidual D istribution and the A ctive F lux formulations or a new class of schemes that can combine several writings of the same hyperbolic problem: application to the 1d E uler equations. Communications on Applied Mathematics and Computation , pages 1--33, 2022

2022

[4] [4]

A semi-discrete A ctive F lux method for the E uler equations on C artesian grids

Rémi Abgrall, Wasilij Barsukow, and Christian Klingenberg. A semi-discrete A ctive F lux method for the E uler equations on C artesian grids. accepted 2024 in J. Sci. Comp., arXiv preprint arXiv:2310.00683 , 2024

arXiv 2024

[5] [5]

Stationarity preserving schemes for multi-dimensional linear systems

Wasilij Barsukow. Stationarity preserving schemes for multi-dimensional linear systems. Mathematics of Computation , 88(318):1621--1645, 2019

work page 2019

[6] [6]

Stationarity preservation properties of the active flux scheme on cartesian grids

Wasilij Barsukow. Stationarity preservation properties of the active flux scheme on cartesian grids. Communications on Applied Mathematics and Computation , pages 1--15, 2020

work page 2020

[7] [7]

The active flux scheme for nonlinear problems

Wasilij Barsukow. The active flux scheme for nonlinear problems. Journal of Scientific Computing , 86(1):1--34, 2021

2021

[8] [8]

The active flux scheme on C artesian grids and its low M ach number limit

Wasilij Barsukow, Jonathan Hohm, Christian Klingenberg, and Philip L Roe. The active flux scheme on C artesian grids and its low M ach number limit. Journal of Scientific Computing , 81(1):594--622, 2019

2019