pith. sign in

arxiv: 2501.02416 · v1 · submitted 2025-01-05 · ⚛️ physics.flu-dyn · physics.comp-ph

Wave or Physics-Appropriate Multidimensional Upwinding Approach for Compressible Multiphase Flows

Pith reviewed 2026-05-23 06:22 UTC · model grok-4.3

classification ⚛️ physics.flu-dyn physics.comp-ph
keywords multidimensional upwindingmultiphase compressible flowsEuler equationscharacteristic reconstructionTHINC schemenumerical artifactsshock capturingvortical structures
0
0 comments X

The pith

A new multidimensional upwinding method treats acoustic, vorticity, and entropy waves with separate reconstruction schemes to cut numerical artifacts in multiphase flows.

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

The paper introduces multidimensional algorithms that combine the wave structure of the Euler equations in characteristic space with physical-space properties for compressible multiphase flow simulation. Acoustic waves receive upwind reconstruction, vorticity and entropy waves receive central reconstruction, and material interfaces plus contact discontinuities receive THINC reconstruction; in physical space, phasic densities use THINC while tangential velocities use central schemes, with an adaptive switch to primitive variables near shocks identified via the stiffened gas parameter. The approach is shown to eliminate spurious vortices in shear layers, preserve vortical structures in gas-gas and gas-liquid cases, and improve shock-entropy interactions. A sympathetic reader would care because these changes produce visibly cleaner flow features and closer experimental agreement without adding extra dissipation or complexity to the base solver.

Core claim

The central claim is that applying upwind schemes to acoustic waves, central schemes to vorticity and entropy waves, and THINC reconstruction to interfaces and contacts in characteristic space, together with physical-space THINC on densities and central differencing on tangential velocities plus adaptive reconstruction using the stiffened gas parameter, yields higher accuracy, fewer numerical artifacts, and better experimental fidelity than conventional schemes across the tested multiphase cases.

What carries the argument

Characteristic decomposition of the Euler equations that isolates acoustic, vorticity, and entropy waves, each receiving a distinct reconstruction stencil (upwind for acoustic, central for vorticity/entropy, THINC for interfaces), augmented by physical-space THINC on phasic densities and an adaptive primitive/characteristic switch triggered by the stiffened gas parameter.

If this is right

  • Spurious vortices disappear from periodic shear-layer tests.
  • Vortical structures remain sharp in gas-gas and gas-liquid mixing problems.
  • Shock-entropy wave interactions show reduced post-shock oscillations and higher amplitude preservation.
  • Complex test cases match experimental data more closely than traditional one-dimensional upwinding schemes.

Where Pith is reading between the lines

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

  • The separation of wave families could be applied to other hyperbolic systems that admit a similar characteristic structure.
  • The adaptive switch between primitive and characteristic variables might reduce the need for special limiters near material interfaces in existing codes.
  • Extending the same wave-type logic to three-dimensional unstructured grids would test whether the artifact reduction persists when mesh topology changes.
  • The stiffened-gas identification step could be replaced by other phase indicators if the assumption that the parameter remains reliable near strong shocks proves fragile.

Load-bearing premise

The method assumes that the characteristic decomposition of the Euler equations continues to cleanly separate acoustic, vorticity, and entropy waves even in multiphase regions and that the stiffened gas parameter can correctly flag the liquid phase near discontinuities without misclassification or added error.

What would settle it

Running the standard shock-bubble interaction or gas-liquid Riemann problem with both the new scheme and a conventional MUSCL or WENO scheme and measuring whether the new scheme produces measurably lower spurious vorticity or closer agreement with the reference interface position and pressure trace.

Figures

Figures reproduced from arXiv: 2501.02416 by Amareshwara Sainadh Chamarthi.

Figure 1
Figure 1. Figure 1: Numerical solution for shock interface interaction problem in Example 4.1 on a grid size of [PITH_FULL_IMAGE:figures/full_fig_p013_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Numerical solution for shock interface interaction problem in Example 4.1 on a grid size of [PITH_FULL_IMAGE:figures/full_fig_p014_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Numerical solution for isolated contact test case using [PITH_FULL_IMAGE:figures/full_fig_p015_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Numerical solution for Liquid-gas shock tube problem in Example 4.3 on a grid size of [PITH_FULL_IMAGE:figures/full_fig_p016_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: z-vorticity contours on a grid size of 8002 , Example 4.4 [PITH_FULL_IMAGE:figures/full_fig_p017_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: z-vorticity contours of the considered schemes on a grid size of 3202 , Example 4.4. (a) WENO5 (Prim). (b) MP5 (Prim). (c) Wave-MP (Prim) [PITH_FULL_IMAGE:figures/full_fig_p018_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: z-vorticity contours of the considered schemes (using primitive variables) computed on a grid size of 3202 , Example 4.4 [PITH_FULL_IMAGE:figures/full_fig_p018_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Figure is taken from Reference [48], where the simulations are computed on a grid size of 512 [PITH_FULL_IMAGE:figures/full_fig_p018_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Figure shows the z-vorticity contours for the TENO, TENO-THINC and modified TENO-THINC (THINC for entropy wave only) schemes using a grid size of 3202 , Example 4.4. Example 4.5. Compressible triple point This test case examines the multi-species compressible triple-point problem, a two-dimensional Riemann problem involving three states and two distinct materials. The primary objective is to demonstrate th… view at source ↗
Figure 10
Figure 10. Figure 10: Density gradient contours at time t = 5 using various schemes, Example 4.5. 21 [PITH_FULL_IMAGE:figures/full_fig_p021_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Figures 11(a) and 11(b) show sensor location regions. [PITH_FULL_IMAGE:figures/full_fig_p022_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: Density gradient contours at time t = 5 using primitive variable reconstruction schemes, Example 4.5. Example 4.6. 2D shock-entropy wave test This test case demonstrates the impact of reconstructing the entropy wave on the numerical solution. This test case is the two-dimensional shock-entropy wave interaction problem, modified for the multi-species case, proposed in [39] with the following initial condit… view at source ↗
Figure 13
Figure 13. Figure 13: The fine grid solution is shown in Figure 13(a) and coarse grid results computed using Wave-MP are shown in 13(b). [PITH_FULL_IMAGE:figures/full_fig_p023_13.png] view at source ↗
Figure 14
Figure 14. Figure 14: Figure 14(a) shows the local density profiles for the MP5 and Wave-MP schemes, and Figure 14(b) shows the local [PITH_FULL_IMAGE:figures/full_fig_p024_14.png] view at source ↗
Figure 15
Figure 15. Figure 15: Volume fraction contours are shown in 15(a), Example 4.6. Figure 15(b) shows the local volume fraction at [PITH_FULL_IMAGE:figures/full_fig_p025_15.png] view at source ↗
Figure 16
Figure 16. Figure 16: Discontinuity detection locations from various papers using THINC from the literature. Figure 16(a) is reproduced [PITH_FULL_IMAGE:figures/full_fig_p026_16.png] view at source ↗
Figure 17
Figure 17. Figure 17: Cylindrical cavity collapse, Example 4.7, using Wave-MP approach. [PITH_FULL_IMAGE:figures/full_fig_p027_17.png] view at source ↗
Figure 18
Figure 18. Figure 18: Nonlinear function of normalized density gradient magnitude, [PITH_FULL_IMAGE:figures/full_fig_p028_18.png] view at source ↗
Figure 19
Figure 19. Figure 19: Nonlinear function of normalized density gradient magnitude, [PITH_FULL_IMAGE:figures/full_fig_p029_19.png] view at source ↗
Figure 20
Figure 20. Figure 20: Shock-water cylinder interation, Example 4.9, comparison of numerical results with experiment. [PITH_FULL_IMAGE:figures/full_fig_p030_20.png] view at source ↗
Figure 21
Figure 21. Figure 21: Shock-water cylinder interaction, Example 4.9. Numerical and experimental result at [PITH_FULL_IMAGE:figures/full_fig_p030_21.png] view at source ↗
Figure 22
Figure 22. Figure 22: Shock-water cylinder interaction, Example 4.9, on a 3072 [PITH_FULL_IMAGE:figures/full_fig_p031_22.png] view at source ↗
Figure 23
Figure 23. Figure 23: Shock-water cylinder interaction, Example 4.9, on a 3072 [PITH_FULL_IMAGE:figures/full_fig_p032_23.png] view at source ↗
Figure 24
Figure 24. Figure 24: Normalised kinetic energy and enstrophy using various schemes on grid size of 64 [PITH_FULL_IMAGE:figures/full_fig_p034_24.png] view at source ↗
Figure 25
Figure 25. Figure 25: Results from the literature. Figure 25(a) is reproduced from [41] with permission from Elsevier BV 2024, Li [PITH_FULL_IMAGE:figures/full_fig_p034_25.png] view at source ↗
read the original abstract

This paper introduces multidimensional algorithms for simulating multiphase flows, leveraging the wave structure of the Euler equations in characteristic space and the physical properties of variables in physical space. The algorithm applies different reconstruction schemes to acoustic, vorticity, and entropy waves in characteristic space to enhance accuracy and minimize numerical artifacts. In characteristic space, upwind schemes are used for acoustic waves, central schemes for vorticity and entropy waves, and Tangent of Hyper-bola for INterface Capturing (THINC) reconstruction for material interfaces and contact discontinuities (a subset of entropy waves). This approach prevents spurious vortices in periodic shear layers, accurately captures vortical structures in gas-gas and gas-liquid interactions, and improves the accuracy of shock-entropy wave interactions. In physical space, phasic densities are computed using THINC in regions of contact discontinuities and material interfaces, while tangential velocities are calculated with central schemes to improve vortical structures. An adaptive reconstruction technique is also introduced to mitigate oscillations near shocks, which arise from primitive variable reconstruction, by combining primitive and characteristic variable reconstructions with the liquid phase being identified using the stiffened gas parameter. The proposed multidimensional upwinding approach outperforms traditional schemes, demonstrating superior accuracy in capturing physical phenomena, reducing numerical artifacts, and better matching experimental results across complex test cases.

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

1 major / 0 minor

Summary. The paper introduces multidimensional upwinding algorithms for compressible multiphase flows that combine wave-structure decomposition in characteristic space (upwind schemes for acoustic waves, central schemes for vorticity and entropy waves, THINC for material interfaces) with physics-based reconstructions in physical space (THINC for phasic densities at contacts, central schemes for tangential velocities). An adaptive blend of primitive and characteristic variable reconstructions is proposed to suppress oscillations near shocks, with the liquid phase identified via the stiffened-gas parameter. The method is claimed to eliminate spurious vortices, improve vortical structure capture in gas-gas and gas-liquid interactions, and yield superior accuracy and experimental agreement relative to traditional schemes across complex test cases.

Significance. If the central claims hold after verification, the approach could reduce common numerical artifacts (spurious vorticity, interface oscillations) in multiphase compressible flow simulations without introducing additional free parameters. The use of standard wave decomposition and THINC reconstruction is a strength, but the absence of derivations, error analysis, or quantitative metrics in the abstract limits assessment of whether the gains are robust or merely case-specific.

major comments (1)
  1. [Abstract] Abstract (adaptive reconstruction paragraph): the central claim that the primitive-characteristic blend mitigates oscillations near shocks depends on the stiffened-gas parameter reliably flagging the liquid phase for stencil selection. Because this parameter varies continuously across numerically diffused interfaces, any threshold-based switch risks selecting the wrong reconstruction in cells with volume fractions between 0.1 and 0.9; the manuscript provides no validation or analysis of the identification logic in this transition regime, which directly undermines the oscillation-control assertion.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the careful review and the specific comment on the adaptive reconstruction description. We address the point directly below and will revise the manuscript accordingly to strengthen the presentation of the method.

read point-by-point responses
  1. Referee: [Abstract] Abstract (adaptive reconstruction paragraph): the central claim that the primitive-characteristic blend mitigates oscillations near shocks depends on the stiffened-gas parameter reliably flagging the liquid phase for stencil selection. Because this parameter varies continuously across numerically diffused interfaces, any threshold-based switch risks selecting the wrong reconstruction in cells with volume fractions between 0.1 and 0.9; the manuscript provides no validation or analysis of the identification logic in this transition regime, which directly undermines the oscillation-control assertion.

    Authors: We acknowledge that the current manuscript text does not include explicit validation or analysis of the stiffened-gas parameter threshold behavior across numerically diffused interfaces with intermediate volume fractions. The identification logic relies on the parameter value to select between primitive and characteristic reconstructions near shocks, but we agree that demonstrating its robustness in the 0.1–0.9 volume-fraction range would directly support the oscillation-mitigation claim. We will add a short subsection (or appendix) with supporting analysis, including the specific threshold criterion employed and results from representative interface-crossing cells, in the revised manuscript. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation is self-contained from standard wave decomposition

full rationale

The paper constructs its multidimensional upwinding scheme directly from the characteristic decomposition of the Euler equations, applying established upwind/central/THINC reconstructions to acoustic, vorticity, and entropy waves plus an adaptive primitive-characteristic blend flagged by the stiffened-gas parameter. No equation or claim reduces by construction to a fitted parameter renamed as prediction, no self-citation supplies a load-bearing uniqueness theorem, and no ansatz is smuggled via prior work by the same authors. The algorithm is presented as a novel combination of existing techniques whose performance is then demonstrated on test cases; the derivation chain therefore remains independent of its own outputs.

Axiom & Free-Parameter Ledger

0 free parameters · 3 axioms · 0 invented entities

The approach rests on standard domain assumptions from hyperbolic PDE theory and CFD practice; no free parameters or invented entities are explicitly introduced in the abstract.

axioms (3)
  • domain assumption The compressible Euler equations admit a characteristic decomposition separating acoustic, vorticity, and entropy waves.
    Invoked in the description of applying different schemes in characteristic space.
  • domain assumption THINC reconstruction preserves sharp material interfaces and contact discontinuities without excessive diffusion.
    Used for interfaces in both characteristic and physical space reconstructions.
  • domain assumption The stiffened gas equation of state parameter can identify the liquid phase for adaptive reconstruction decisions.
    Cited for mitigating oscillations near shocks via phase identification.

pith-pipeline@v0.9.0 · 5760 in / 1531 out tokens · 54448 ms · 2026-05-23T06:22:14.516471+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

57 extracted references · 57 canonical work pages

  1. [1]

    S. K. Godunov, A difference method for numerical calculation of discontinuous solutions of the equations of hydrodynamics, Matematicheskii Sbornik 89 (3) (1959) 271–306

  2. [2]

    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. 34

  3. [3]

    X. X.-D. Liu, S. Osher, T. Chan, Weighted essentially non-oscillatory schemes (1994)

  4. [4]

    Jiang, C.-W

    G.-S. Jiang, C.-W. Shu, Efficient Implementation of Weighted ENO Schemes, Journal of Computational Physics 126 (126) (1995) 202–228

  5. [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. [6]

    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

  7. [7]

    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

  8. [8]

    Suresh, H

    A. Suresh, H. Huynh, Accurate monotonicity-preserving schemes with runge-kutta time stepping, Jour- nal of Computational Physics 136 (1) (1997) 83–99

  9. [9]

    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

  10. [10]

    P. L. Roe, Discrete models for the numerical analysis of time-dependent multidimensional gas dynamics, Journal of Computational Physics 63 (2) (1986) 458–476

  11. [11]

    Hirsch, Numerical Computation of Internal and External Flows, Volume 2: Computational Methods for Inviscid and Viscous Flows, Wiley, 1990

    C. Hirsch, Numerical Computation of Internal and External Flows, Volume 2: Computational Methods for Inviscid and Viscous Flows, Wiley, 1990

  12. [12]

    Abgrall, How to prevent pressure oscillations in multicomponent flow calculations: a quasi conser- vative approach, Journal of Computational Physics 125 (1) (1996) 150–160

    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

  13. [13]

    Abgrall, S

    R. Abgrall, S. Karni, Computations of compressible multifluids, Journal of computational physics 169 (2) (2001) 594–623

  14. [14]

    Chargy, R

    D. Chargy, R. Abgrall, L. F. Fezoui, B. Larrouturou, Comparisons of several upwind schemes for multi-component one-dimensional inviscid flows, Ph.D. thesis, INRIA (1990)

  15. [15]

    G. K. Batchelor, An introduction to fluid dynamics, Cambridge university press, 1967

  16. [16]

    J. C. Meng, T. Colonius, Numerical simulation of the aerobreakup of a water droplet, Journal of Fluid Mechanics 835 (2018) 1108–1135

  17. [17]

    Coralic, T

    V. Coralic, T. Colonius, Finite-volume weno scheme for viscous compressible multicomponent flows, Journal of Computational Physics 274 (2014) 95–121

  18. [18]

    S. H. Bryngelson, K. Schmidmayer, V. Coralic, J. C. Meng, K. Maeda, T. Colonius, Mfc: An open- source high-order multi-component, multi-phase, and multi-scale compressible flow solver, Computer physics communications 266 (2021) 107396

  19. [19]

    Hoffmann, A

    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

  20. [20]

    van Noordt, S

    W. van Noordt, S. Ganju, C. Brehm, An immersed boundary method for wall-modeled large-eddy simulation of turbulent high-mach-number flows, Journal of Computational Physics 470 (2022) 111583

  21. [21]

    T. R. Hendrickson, A. Kartha, G. V. Candler, An improved ducros sensor for the simulation of com- pressible flows with shocks, in: 2018 Fluid Dynamics Conference, 2018, p. 3710. 35

  22. [22]

    Y. Kuya, K. Totani, S. Kawai, Kinetic energy and entropy preserving schemes for compressible flows by split convective forms, Journal of Computational Physics 375 (2018) 823–853

  23. [23]

    A. S. Chamarthi, N. Hoffmann, S. Bokor, S. H. Frankel, et al., On the role of spectral properties of viscous flux discretization for flow simulations on marginally resolved grids, Computers & Fluids 251 (2023) 105742

  24. [24]

    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

  25. [25]

    A. S. Chamarthi, N. Hoffmann, S. Frankel, A wave appropriate discontinuity sensor approach for compressible flows, Physics of Fluids 35 (6) (2023)

  26. [26]

    A. S. Chamarthi, Efficient high-order gradient-based reconstruction for compressible flows, Journal of Computational Physics 486 (2023) 112119

  27. [27]

    Sainadh Chamarthi, Consistent interface capturing adaptive reconstruction approach for viscous compressible multicomponent flows, arXiv e-prints (2024) arXiv–2410

    A. Sainadh Chamarthi, Consistent interface capturing adaptive reconstruction approach for viscous compressible multicomponent flows, arXiv e-prints (2024) arXiv–2410

  28. [28]

    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

  29. [29]

    Shyue, F

    K.-M. Shyue, F. Xiao, An eulerian interface sharpening algorithm for compressible two-phase flow: The algebraic thinc approach, Journal of Computational Physics 268 (2014) 326–354

  30. [30]

    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

  31. [31]

    X. Deng, B. Xie, R. Loub` ere, Y. Shimizu, F. Xiao, Limiter-free discontinuity-capturing scheme for compressible gas dynamics with reactive fronts, Computers & Fluids 171 (2018) 1–14

  32. [32]

    Allaire, S

    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

  33. [33]

    Toro, Riemann Solvers and Numerical Methods for Fluid Dynamics: A Practical Introduction, Springer Berlin Heidelberg, 2009

    E. Toro, Riemann Solvers and Numerical Methods for Fluid Dynamics: A Practical Introduction, Springer Berlin Heidelberg, 2009

  34. [34]

    Johnsen, T

    E. Johnsen, T. Colonius, Implementation of weno schemes in compressible multicomponent flow prob- lems, Journal of Computational Physics 219 (2) (2006) 715–732

  35. [35]

    Van Leer, Towards the ultimate conservative difference scheme

    B. Van Leer, Towards the ultimate conservative difference scheme. iv. a new approach to numerical convection, Journal of Computational Physics 23 (3) (1977) 276–299

  36. [36]

    van Leer, H

    B. van Leer, H. Nishikawa, Towards the ultimate understanding of muscl: Pitfalls in achieving third- order accuracy, Journal of Computational Physics (2021) 110640

  37. [37]

    Wakimura, S

    H. Wakimura, S. Takagi, F. Xiao, Symmetry-preserving enforcement of low-dissipation method based on boundary variation diminishing principle, Computers & Fluids 233 (2022) 105227

  38. [38]

    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

  39. [39]

    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

  40. [40]

    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. 36

  41. [41]

    Takagi, L

    S. Takagi, L. Fu, H. Wakimura, F. Xiao, A novel high-order low-dissipation teno-thinc scheme for hyperbolic conservation laws, Journal of Computational Physics 452 (2022) 110899

  42. [42]

    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

  43. [43]

    Ducros, V

    F. Ducros, V. Ferrand, F. Nicoud, C. Weber, D. Darracq, C. Gacherieu, T. Poinsot, Large-eddy simu- lation of the shock/turbulence interaction, Journal of Computational Physics 152 (2) (1999) 517–549

  44. [44]

    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

  45. [45]

    M. L. Wong, J. B. Angel, M. F. Barad, C. C. Kiris, A positivity-preserving high-order weighted compact nonlinear scheme for compressible gas-liquid flows, Journal of Computational Physics 444 (2021) 110569

  46. [46]

    D. L. Brown, Performance of under-resolved two-dimensional incompressible flow simulations, Journal of Computational Physics 122 (1) (1995) 165–183

  47. [47]

    M. L. Minion, D. L. Brown, Performance of under-resolved two-dimensional incompressible flow simu- lations, ii, Journal of Computational Physics 138 (2) (1997) 734–765

  48. [48]

    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

  49. [49]

    Cheng, X

    L. Cheng, X. Deng, B. Xie, Y. Jiang, F. Xiao, Low-dissipation bvd schemes for single and multi-phase compressible flows on unstructured grids, Journal of Computational Physics 428 (2021) 110088

  50. [50]

    Huang, X

    H. Huang, X. Li, L. Fu, A new high-order rkdg method based on the teno-thinc scheme for shock- capturing, Journal of Computational Physics 520 (2025) 113459

  51. [51]

    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

  52. [52]

    P. T. Barton, An interface-capturing godunov method for the simulation of compressible solid-fluid problems, Journal of Computational Physics 390 (2019) 25–50

  53. [53]

    Sembian, M

    S. Sembian, M. Liverts, N. Tillmark, N. Apazidis, Plane shock wave interaction with a cylindrical water column, Physics of Fluids 28 (5) (2016)

  54. [54]

    A. N. Barrett, P. K. Subbareddy, G. V. Candler, Development of a low-dissipation diffuse interface method for compressible multiphase flow, in: AIAA SCITECH 2024 Forum, 2024, p. 1756

  55. [55]

    Zhang, N

    W. Zhang, N. Fleischmann, S. Adami, N. A. Adams, A hybrid weno5is-thinc reconstruction scheme for compressible multiphase flows, Journal of Computational Physics 498 (2024) 112672

  56. [56]

    Y. Li, L. Fu, N. A. Adams, A low-dissipation shock-capturing framework with flexible nonlinear dissi- pation control, Journal of Computational Physics 428 (2021) 109960

  57. [57]

    T. Yang, G. Zhao, Q. Zhao, Novel teno schemes with improved accuracy order based on perturbed polynomial reconstruction, Journal of Computational Physics 488 (2023) 112219. 37