Consistent Interface Capturing Adaptive Reconstruction Approach for Viscous Compressible Multicomponent Flows
Pith reviewed 2026-05-23 19:27 UTC · model grok-4.3
The pith
A discretization method applies THINC reconstruction only at detected contact discontinuities and central schemes to tangential velocities to capture material interfaces sharply in viscous multicomponent flows.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The proposed discretization uses a contact discontinuity detector to apply the THINC reconstruction selectively to the contact wave, volume fractions, and phasic densities, while reconstructing tangential velocities with a central scheme wherever the Ducros sensor indicates a material interface; this combination is claimed to capture interfaces sharply without oscillations because tangential velocities are continuous across such interfaces in viscous flows.
What carries the argument
Contact discontinuity detector that triggers THINC reconstruction for entropy waves and phasic densities, paired with Ducros-sensor-driven central reconstruction for tangential velocities.
If this is right
- Dissipation errors near contact discontinuities are reduced by limiting THINC reconstruction to the entropy wave and phasic densities.
- No oscillations appear at material interfaces when tangential velocities are reconstructed centrally using the Ducros sensor.
- Material interfaces are captured more sharply than with standard MP or WENO reconstructions alone in the benchmark tests performed.
- The method remains stable for both inviscid and viscous compressible multicomponent flows in the reported cases.
Where Pith is reading between the lines
- The same selective-reconstruction logic could be tested on three-dimensional unstructured grids to check whether interface sharpness is preserved under mesh irregularity.
- If the Ducros sensor misses weak contacts in certain regimes, replacing it with a dedicated material-interface sensor might further improve robustness.
Load-bearing premise
That applying a central reconstruction scheme to tangential velocities across material interfaces identified by the Ducros sensor will produce no oscillations.
What would settle it
A viscous multicomponent flow simulation in which the central tangential-velocity reconstruction across a detected material interface generates visible oscillations or excessive smearing compared with existing interface-capturing schemes.
Figures
read the original abstract
The paper proposes a physically consistent numerical discretization approach for simulating viscous compressible multicomponent flows. It has two main contributions. First, a contact discontinuity (and material interface) detector is developed. In those regions of contact discontinuities, the THINC (Tangent of Hyperbola for INterface Capturing) approach is used for reconstructing appropriate variables (phasic densities). For other flow regions, the variables are reconstructed using the Monotonicity-preserving (MP) scheme (or Weighted essentially non-oscillatory scheme (WENO)). For reconstruction in the characteristic space, the THINC approach is used only for the contact (or entropy) wave and volume fractions and for the reconstruction of primitive variables, the THINC approach is used for phasic densities and volume fractions only, offering an effective solution for reducing dissipation errors near contact discontinuities. The second contribution is the development of an algorithm that uses a central reconstruction scheme for the tangential velocities, as they are continuous across material interfaces in viscous flows. In this regard, the Ducros sensor (a shock detector that cannot detect material interfaces) is employed to compute the tangential velocities using a central scheme across material interfaces. Using the central scheme does not produce any oscillations at the material interface. The proposed approach is thoroughly validated with several benchmark test cases for compressible multicomponent flows, highlighting its advantages. The numerical results of the benchmark tests show that the proposed method captured the material interface sharply compared to existing methods.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript proposes a numerical discretization for viscous compressible multicomponent flows with two contributions: (1) a contact-discontinuity detector that applies THINC reconstruction to phasic densities and volume fractions (and to the contact wave in characteristic space) while using MP or WENO elsewhere, and (2) a Ducros-sensor-triggered central reconstruction for tangential velocities across material interfaces, asserted to produce no oscillations because the sensor cannot detect interfaces. The method is claimed to be validated on several benchmark tests, with results showing sharper material-interface capture than existing methods.
Significance. If the central claims hold, the approach could reduce numerical dissipation at material interfaces in viscous multicomponent flows without introducing oscillations, which would be useful for high-fidelity simulations in aerodynamics and combustion. The algorithmic construction itself introduces no free parameters or circular definitions.
major comments (2)
- [Abstract] Abstract (second contribution): the assertion that 'Using the central scheme does not produce any oscillations at the material interface' is unsupported by any derivation, eigenvalue analysis, modified-equation argument, or isolated numerical test. No with/without comparison for the tangential-velocity treatment is described, nor is any oscillation metric (e.g., max |Δu_tang| across the interface) reported. If this claim does not hold, the stated advantage of the second contribution is lost.
- [Abstract] Abstract (validation statement): the claim that 'the numerical results of the benchmark tests show that the proposed method captured the material interface sharply compared to existing methods' supplies neither quantitative error norms, convergence rates, nor tabulated comparisons. Without these data the improvement cannot be evaluated and the central claim remains unverified.
minor comments (1)
- [Abstract] The abstract does not name the specific benchmark tests or the existing methods used for comparison, which would help readers assess the scope of the validation.
Simulated Author's Rebuttal
We thank the referee for the constructive and detailed review. We address each major comment point-by-point below, indicating where revisions will be made to strengthen the manuscript.
read point-by-point responses
-
Referee: [Abstract] Abstract (second contribution): the assertion that 'Using the central scheme does not produce any oscillations at the material interface' is unsupported by any derivation, eigenvalue analysis, modified-equation argument, or isolated numerical test. No with/without comparison for the tangential-velocity treatment is described, nor is any oscillation metric (e.g., max |Δu_tang| across the interface) reported. If this claim does not hold, the stated advantage of the second contribution is lost.
Authors: We agree that the abstract statement lacks explicit supporting analysis or tests in its current form. The physical basis is that tangential velocity is continuous across material interfaces for viscous flows, and the Ducros sensor (a shock detector) does not trigger at interfaces, so the central scheme is applied only where the velocity field remains smooth. To address the concern rigorously, we will add an isolated numerical test in the revised manuscript comparing tangential velocity profiles across an interface with and without the central scheme, including a quantitative oscillation metric such as max |Δu_tang|. A brief explanatory paragraph will also be inserted in the methods section. revision: yes
-
Referee: [Abstract] Abstract (validation statement): the claim that 'the numerical results of the benchmark tests show that the proposed method captured the material interface sharply compared to existing methods' supplies neither quantitative error norms, convergence rates, nor tabulated comparisons. Without these data the improvement cannot be evaluated and the central claim remains unverified.
Authors: We concur that quantitative metrics are necessary to substantiate the improvement claims. In the revised manuscript we will augment the validation section with L1 and L2 error norms for key variables (e.g., density and volume fraction), observed convergence rates on successively refined meshes, and tabulated side-by-side comparisons against the reference methods used in the benchmarks. These additions will allow direct, objective evaluation of the interface-capturing gains. revision: yes
Circularity Check
No circularity detected; derivation is an independent algorithmic construction
full rationale
The paper defines an algorithmic discretization that combines a contact-discontinuity detector with THINC reconstruction for phasic densities and volume fractions, MP/WENO elsewhere, and Ducros-triggered central reconstruction for tangential velocities. No equation or central claim reduces by construction to a fitted parameter, self-definition, or self-citation chain; the assertion that the central scheme produces no oscillations is presented as a design choice without supporting derivation, yet this does not create circularity because the result is not forced to equal its inputs. Benchmarks supply external empirical checks, rendering the method self-contained against the enumerated circularity patterns.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption Central reconstruction of tangential velocities across material interfaces identified by the Ducros sensor produces no oscillations
Reference graph
Works this paper leans on
-
[1]
C. Hirsch, Numerical Computation of Internal and External Flows, Volume 2: Computational Methods for Inviscid and Viscous Flows, Wiley, 1990. 37
work page 1990
-
[2]
G. K. Batchelor, An introduction to fluid dynamics, Cambridge university press, 1967
work page 1967
-
[3]
G.-S. Jiang, C.-W. Shu, Efficient Implementation of Weighted ENO Schemes, Journal of Computational Physics 126 (126) (1995) 202–228
work page 1995
-
[4]
V. Coralic, T. Colonius, Finite-volume weno scheme for viscous compressible multicomponent flows, Journal of Computational Physics 274 (2014) 95–121
work page 2014
-
[5]
X. Y. Hu, Q. Wang, N. A. Adams, An adaptive central-upwind weighted essentially non-oscillatory scheme, Journal of Computational Physics 229 (23) (2010) 8952–8965. doi:10.1016/j.jcp.2010.08. 019
- [6]
-
[7]
A. S. Chamarthi, S. H. Frankel, High-order central-upwind shock capturing scheme using a boundary variation diminishing (bvd) algorithm, Journal of Computational Physics 427 (2021) 110067
work page 2021
-
[8]
A. S. Chamarthi, Gradient based reconstruction: Inviscid and viscous flux discretizations, shock cap- turing, and its application to single and multicomponent flows, Computers & Fluids 250 (2023) 105706
work page 2023
- [9]
- [10]
- [11]
-
[12]
X.-D. Liu, S. Osher, T. Chan, Weighted essentially non-oscillatory schemes, Journal of Computational Physics 115 (1) (1994) 200–212
work page 1994
-
[13]
X. Y. Hu, N. A. Adams, Scale separation for implicit large eddy simulation, Journal of Computational Physics 230 (19) (2011) 7240–7249. doi:10.1016/j.jcp.2011.05.023
-
[14]
N. A. Adams, K. Shariff, A high-resolution hybrid compact-eno scheme for shock-turbulence interaction problems, Journal of Computational Physics 127 (1) (1996) 27–51
work page 1996
-
[15]
X. Liu, S. Zhang, H. Zhang, C. W. Shu, A new class of central compact schemes with spectral-like resolution II: Hybrid weighted nonlinear schemes, Journal of Computational Physics 284 (2015) 133– 154
work page 2015
-
[16]
D. S. Balsara, S. Garain, C. W. Shu, An efficient class of WENO schemes with adaptive order, Journal of Computational Physics 326 (2016) 780–804. doi:10.1016/j.jcp.2016.09.009. URL http://dx.doi.org/10.1016/j.jcp.2016.09.009
-
[17]
M. L. Wong, S. K. Lele, High-order localized dissipation weighted compact nonlinear scheme for shock- and interface-capturing in compressible flows, Vol. 339, Elsevier Inc., 2017. doi:10.1016/j.jcp.2017. 03.008
-
[18]
T. Nonomura, K. Fujii, Characteristic finite-difference WENO scheme for multicomponent com- pressible fluid analysis: Overestimated quasi-conservative formulation maintaining equilibriums of velocity, pressure, and temperature, Journal of Computational Physics 340 (2017) 358–388. doi: 10.1016/j.jcp.2017.02.054
-
[19]
T. Nonomura, S. Morizawa, H. Terashima, S. Obayashi, K. Fujii, Numerical (error) issues on com- pressible multicomponent flows using a high-order differencing scheme: Weighted compact nonlinear scheme, Journal of Computational Physics 231 (8) (2012) 3181–3210. 38
work page 2012
-
[20]
D. S. Balsara, C.-W. Shu, Monotonicity Preserving Weighted Essentially Non-oscillatory Schemes with Increasingly High Order of Accuracy, Journal of Computational Physics 160 (2000) 405–452. doi: 10.1006/jcph.2000.6443
-
[21]
R. K. Shukla, C. Pantano, J. B. Freund, An interface capturing method for the simulation of multi-phase compressible flows, Journal of Computational Physics 229 (19) (2010) 7411–7439
work page 2010
-
[22]
A. Chiapolino, R. Saurel, B. Nkonga, Sharpening diffuse interfaces with compressible fluids on unstruc- tured meshes, Journal of Computational Physics 340 (2017) 389–417
work page 2017
-
[23]
Harten, The artificial compression method for computation of shocks and contact discontinuities
A. Harten, The artificial compression method for computation of shocks and contact discontinuities. i-single conservation laws, Communications in Pure Applied Mathematics 30 (1977) 611–638
work page 1977
-
[24]
H. Yang, An artificial compression method for eno schemes: the slope modification method, Journal of Computational Physics 89 (1) (1990) 125–160
work page 1990
-
[25]
Z. He, B. Tian, Y. Zhang, F. Gao, Characteristic-based and interface-sharpening algorithm for high- order simulations of immiscible compressible multi-material flows, Journal of Computational Physics 333 (2017) 247–268
work page 2017
-
[26]
Harten, Eno schemes with subcell resolution, Journal of Computational Physics 83 (1) (1989) 148– 184
A. Harten, Eno schemes with subcell resolution, Journal of Computational Physics 83 (1) (1989) 148– 184
work page 1989
-
[27]
H. T. Huynh, Accurate upwind methods for the euler equations, SIAM Journal on Numerical Analysis 32 (5) (1995) 1565–1619
work page 1995
- [28]
-
[29]
D. P. Garrick, W. A. Hagen, J. D. Regele, An interface capturing scheme for modeling atomization in compressible flows, Journal of Computational Physics 344 (2017) 260–280
work page 2017
- [30]
-
[31]
Z. Sun, S. Inaba, F. Xiao, Boundary variation diminishing (bvd) reconstruction: A new approach to improve godunov schemes, Journal of Computational Physics 322 (2016) 309–325
work page 2016
-
[32]
X. Deng, Y. Shimizu, F. Xiao, A fifth-order shock capturing scheme with two-stage boundary variation diminishing algorithm, Journal of Computational Physics 386 (2019) 323–349
work page 2019
-
[33]
X. Deng, S. Inaba, B. Xie, K.-M. Shyue, F. Xiao, High fidelity discontinuity-resolving reconstruction for compressible multiphase flows with moving interfaces, Journal of Computational Physics 371 (2018) 945–966
work page 2018
- [34]
-
[35]
Q. Li, Y. Lv, L. Fu, A high-order diffuse-interface method with teno-thinc scheme for compressible multiphase flows, International Journal of Multiphase Flow 173 (2024) 104732
work page 2024
-
[36]
A. S. Chamarthi, N. Hoffmann, S. Frankel, A wave appropriate discontinuity sensor approach for compressible flows, Physics of Fluids 35 (6) (2023)
work page 2023
- [37]
-
[38]
N. Hoffmann, A. S. Chamarthi, S. H. Frankel, Centralized gradient-based reconstruction for wall mod- elled large eddy simulations of hypersonic boundary layer transition, Journal of Computational Physics (2024) 113128. 39
work page 2024
-
[39]
G. Allaire, S. Clerc, S. Kokh, A five-equation model for the simulation of interfaces between compressible fluids, Journal of Computational Physics 181 (2) (2002) 577–616
work page 2002
-
[40]
H. Nishikawa, Two ways to extend diffusion schemes to navier-stokes schemes: Gradient formula or upwind flux, 20th AIAA Computational Fluid Dynamics Conference 2011 (2011) 27–30
work page 2011
-
[41]
P. Buchm¨ uller, C. Helzel, Improved accuracy of high-order weno finite volume methods on cartesian grids, Journal of Scientific Computing 61 (2) (2014) 343–368
work page 2014
-
[42]
E. Toro, Riemann Solvers and Numerical Methods for Fluid Dynamics: A Practical Introduction, Springer Berlin Heidelberg, 2009
work page 2009
- [43]
- [44]
-
[45]
A. S. Chamarthi, N. Hoffmann, H. Nishikawa, S. H. Frankel, Implicit gradients based conservative numerical scheme for compressible flows, Journal of Scientific Computing 95 (1) (2023) 17
work page 2023
-
[46]
A. S. Chamarthi, Efficient high-order gradient-based reconstruction for compressible flows, Journal of Computational Physics 486 (2023) 112119
work page 2023
-
[47]
van Leer, Towards the ultimate conservative difference scheme
B. van Leer, Towards the ultimate conservative difference scheme. v. a second-order sequel to godunov’s method, Journal of Computational Physics 32 (1) (1979) 101 – 136
work page 1979
-
[48]
F. Xiao, S. Ii, C. Chen, Revisit to the thinc scheme: a simple algebraic vof algorithm, Journal of Computational Physics 230 (19) (2011) 7086–7092
work page 2011
-
[49]
H. Wakimura, S. Takagi, F. Xiao, Symmetry-preserving enforcement of low-dissipation method based on boundary variation diminishing principle, Computers & Fluids 233 (2022) 105227
work page 2022
-
[50]
D. S. Balsara, Higher-order accurate space-time schemes for computational astrophysics—part i: finite volume methods, Living reviews in computational astrophysics 3 (1) (2017) 2
work page 2017
- [51]
-
[52]
H. Nishikawa, From hyperbolic diffusion scheme to gradient method: Implicit green–gauss gradients for unstructured grids, Journal of Computational Physics 372 (2018) 126–160
work page 2018
-
[53]
B. Van Leer, Upwind and high-resolution methods for compressible flow: From donor cell to residual- distribution schemes, in: 16th AIAA Computational Fluid Dynamics Conference, 2003, p. 3559
work page 2003
-
[54]
R. Abgrall, S. Karni, Computations of compressible multifluids, Journal of computational physics 169 (2) (2001) 594–623
work page 2001
-
[55]
E. Johnsen, T. Colonius, Implementation of weno schemes in compressible multicomponent flow prob- lems, Journal of Computational Physics 219 (2) (2006) 715–732
work page 2006
-
[56]
M. L. Wong, S. K. Lele, Improved Weighted Compact Nonlinear Scheme for Flows with Shocks and Material Interfaces: Algorithm and Assessment, 54th AIAA Aerospace Sciences Meeting (January) (2016) 1807. doi:10.2514/6.2016-1807
-
[57]
C. W. Shu, S. Osher, Efficient implementation of essentially non-oscillatory shock-capturing schemes, Journal of Computational Physics 77 (2) (1988) 439–471. doi:10.1016/0021-9991(88)90177-5
-
[58]
Y. Lv, M. Ihme, Discontinuous galerkin method for multicomponent chemically reacting flows and combustion, Journal of Computational Physics 270 (2014) 105–137. 40
work page 2014
-
[59]
Y. Li, C. Chen, Y.-X. Ren, A class of high-order finite difference schemes with minimized dispersion and adaptive dissipation for solving compressible flows, Journal of Computational Physics 448 (2022) 110770
work page 2022
-
[60]
R. Abgrall, How to prevent pressure oscillations in multicomponent flow calculations: a quasi conser- vative approach, Journal of Computational Physics 125 (1) (1996) 150–160
work page 1996
-
[61]
Fu, A hybrid method with teno based discontinuity indicator for hyperbolic conservation laws, Commun
L. Fu, A hybrid method with teno based discontinuity indicator for hyperbolic conservation laws, Commun. Comput. Phys. 26 (4) (2019) 973–1007
work page 2019
-
[62]
L. Krivodonova, J. Xin, J.-F. Remacle, N. Chevaugeon, J. E. Flaherty, Shock detection and limiting with discontinuous galerkin methods for hyperbolic conservation laws, Applied Numerical Mathematics 48 (3-4) (2004) 323–338
work page 2004
-
[63]
G.-Y. Zhao, M.-B. Sun, S. Pirozzoli, On shock sensors for hybrid compact/weno schemes, Computers & Fluids 199 (2020) 104439
work page 2020
- [64]
-
[65]
D. S. Balsara, C.-W. Shu, Monotonicity preserving weighted essentially non-oscillatory schemes with increasingly high order of accuracy, Journal of Computational Physics 160 (2) (2000) 405–452
work page 2000
-
[66]
H. C. Yee, N. D. Sandham, M. J. Djomehri, Low-dissipative high-order shock-capturing methods using characteristic-based filters, Journal of computational physics 150 (1) (1999) 199–238
work page 1999
-
[67]
Y. Feng, J. Winter, N. A. Adams, F. S. Schranner, A general multi-objective bayesian optimization framework for the design of hybrid schemes towards adaptive complex flow simulations, Journal of Computational Physics 510 (2024) 113088
work page 2024
-
[68]
S. Pan, L. Han, X. Hu, N. A. Adams, A conservative interface-interaction method for compressible multi-material flows, Journal of Computational Physics 371 (2018) 870–895
work page 2018
-
[69]
H. C. Yee, B. Sj¨ ogreen, Simulation of richtmyer–meshkov instability by sixth-order filter methods, Shock Waves 17 (3) (2007) 185–193
work page 2007
-
[70]
J. J. Quirk, S. Karni, On the dynamics of a shock–bubble interaction, Journal of Fluid Mechanics 318 (1996) 129–163
work page 1996
- [71]
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.