On numerical stabilization in modeling double-diffusive viscous fingering
Pith reviewed 2026-05-25 00:40 UTC · model grok-4.3
The pith
Stabilized finite element methods for double-diffusive viscous fingering violate non-negativity and the maximum principle and can suppress physical instabilities.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The SUPG formulations and SOLD modifications do not respect the non-negative constraint and the maximum principle for the concentration field. Representative numerical results illustrate that these violations affect viscous finger development, and that the stabilizations often used to suppress numerical instabilities may also suppress physical instabilities such as viscous fingering.
What carries the argument
The SUPG (Streamline Upwind Petrov-Galerkin) formulations and SOLD modifications, which add stabilization to convection-dominated transport equations but fail to preserve physical bounds on the concentration field.
If this is right
- Existing simulations using these methods can produce unphysical negative concentration values.
- Violations of the maximum principle alter the predicted development of viscous fingers.
- Stabilization terms intended to control numerical oscillations can also eliminate physical instabilities.
- Formulations that respect both numerical stability and physical constraints are required for reliable modeling.
Where Pith is reading between the lines
- Similar bound violations may appear when the same stabilized methods are applied to other convection-dominated transport problems.
- Quantifying the difference in finger growth rates between stabilized and constraint-preserving schemes on benchmark cases would measure the practical impact.
- Extending the tests to three-dimensional domains or different diffusivity ratios could show whether the suppression effect persists.
Load-bearing premise
That the representative numerical results are sufficient to demonstrate general limitations of the formulations and that the observed violations meaningfully affect viscous finger development.
What would settle it
A direct comparison on the same double-diffusive fingering problem between a stabilized formulation and a bound-preserving reference method, checking whether negative concentrations appear or whether finger patterns differ.
read the original abstract
A firm understanding and control of viscous fingering (VF) and miscible displacement will be vital to a wide range of industrial, environmental, and pharmaceutical applications, such as geological carbon-dioxide sequestration, enhanced oil recovery, and drug delivery. We restrict our study to VF, a well-known hydrodynamic instability, in miscible fluid systems but consider double-diffusive (DD) effects---the combined effect of compositional changes because of solute transport and temperature. One often uses numerical formulations to study VF with DD effects. The primary aim of the current study is to show that popular formulations have limitations to study VF with DD effect. These limitations include exhibiting node-to-node spurious oscillations, violating physical constraints such as the non-negativity of the concentration field or mathematical principles such as the maximum principle, and suppressing physical instabilities. We will use several popular stabilized finite element formulations---the SUPG formulations and three modifications based on the SOLD approach---in our study. Using representative numerical results, we will illustrate two critical limitations. First, we document that these formulations do not respect the non-negative constraint and the maximum principle for the concentration field. We will also show the impact of these violations on how viscous fingers develop. Second, we show that these stabilized formulations, often used to suppress numerical instabilities, may also suppress physical instabilities, such as viscous fingering. Our study will be valuable to practitioners who use existing numerical formulations and to computational mathematicians who develop new formulations.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript presents an empirical numerical study of stabilized finite-element formulations (SUPG and three SOLD variants) for double-diffusive viscous fingering. It claims that these methods violate the non-negativity constraint and maximum principle on the concentration field and can suppress the physical instability of viscous fingering, with the evidence consisting of representative numerical examples that illustrate node-to-node oscillations, negative concentrations, and damped finger development.
Significance. The identification of practical limitations in widely used stabilization techniques for convection-dominated double-diffusive problems would be useful to practitioners in carbon sequestration and enhanced oil recovery if the violations prove general. The work supplies concrete counter-examples rather than parameter-free derivations or machine-checked proofs, so its impact hinges on whether the reported failures are formulation-inherent.
major comments (2)
- [Abstract] Abstract and numerical-results section: the central claim that the formulations 'do not respect the non-negative constraint and the maximum principle' rests solely on representative examples for selected Rayleigh numbers, mesh aspect ratios, and SOLD parameters. No a-priori bound, positivity proof, or exhaustive parameter sweep is supplied to show that violations must occur for arbitrary discretizations or stabilization coefficients.
- [Numerical results] Numerical-results section: the reported suppression of viscous fingering is demonstrated on specific test cases, yet no comparison against a provably positivity-preserving scheme (e.g., a monotone finite-volume or algebraic-flux-correction method) on identical problems is provided. Consequently it remains unclear whether the damping is caused by the stabilization itself or by the particular choice of time-step, mesh, or physical parameters.
minor comments (1)
- [Abstract] Abstract: repeated use of future tense ('We will use', 'we will illustrate', 'Our study will be valuable') is inconsistent with a completed manuscript.
Simulated Author's Rebuttal
We thank the referee for the detailed and constructive report. We address each major comment below, clarifying the empirical scope of the study and indicating where revisions will be made.
read point-by-point responses
-
Referee: [Abstract] Abstract and numerical-results section: the central claim that the formulations 'do not respect the non-negative constraint and the maximum principle' rests solely on representative examples for selected Rayleigh numbers, mesh aspect ratios, and SOLD parameters. No a-priori bound, positivity proof, or exhaustive parameter sweep is supplied to show that violations must occur for arbitrary discretizations or stabilization coefficients.
Authors: The manuscript presents an empirical investigation whose goal is to document concrete limitations of SUPG and SOLD methods when applied to double-diffusive viscous fingering. The selected Rayleigh numbers, mesh aspect ratios, and stabilization parameters are representative of those encountered in applications such as carbon sequestration. While we do not supply a general positivity proof or exhaustive sweep, the reported examples demonstrate that violations occur and affect the computed finger patterns. We will revise the abstract and the opening of the numerical-results section to state explicitly that the violations are observed for the tested, application-relevant discretizations rather than asserted for all possible choices. revision: partial
-
Referee: [Numerical results] Numerical-results section: the reported suppression of viscous fingering is demonstrated on specific test cases, yet no comparison against a provably positivity-preserving scheme (e.g., a monotone finite-volume or algebraic-flux-correction method) on identical problems is provided. Consequently it remains unclear whether the damping is caused by the stabilization itself or by the particular choice of time-step, mesh, or physical parameters.
Authors: Within each stabilized formulation we vary the stabilization parameter while keeping the mesh, time step, and physical parameters fixed; the resulting increase in damping with stronger stabilization indicates that the added artificial diffusion is responsible. In addition, the appearance of non-physical negative concentrations coincides with the observed suppression, providing internal evidence that the stabilization is altering the solution in a non-physical way. A side-by-side comparison with an algebraic-flux-correction or monotone finite-volume scheme on the identical double-diffusive problem would be informative but lies outside the scope of the present study, which focuses on the behavior of commonly used stabilized finite-element methods. We will add a short discussion of this limitation and of the internal evidence already present in the revised manuscript. revision: partial
Circularity Check
No circularity: empirical numerical demonstration with no derivation chain
full rationale
The paper is an empirical study that uses representative numerical experiments on specific test cases to illustrate limitations (non-negativity violations, maximum principle breaches, and suppression of physical instabilities) in existing SUPG and SOLD stabilized formulations. No first-principles derivation, parameter fitting, or predictive step is claimed or performed; the central claims rest directly on the observed simulation outputs rather than reducing to any self-referential definition, fitted input renamed as prediction, or self-citation load-bearing theorem. Self-citations, if present, are not invoked to justify uniqueness or ansatz choices that close a loop. This is the standard structure of a numerical methods critique paper and remains self-contained against external benchmarks.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption The continuous problem satisfies the maximum principle and non-negativity for the concentration field
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.