Published December 15, 2025 | Version Published
Journal Article

Physics appropriate interface capturing reconstruction approach for viscous compressible multicomponent flows

  • 1. ROR icon California Institute of Technology

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. 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 numerical results of the benchmark tests show that the proposed method captured the material interface sharply compared to existing methods. 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 physics appropriate approach also shown to prevent spurious vortices, despite being formally second-order accurate for nonlinear problems, on a coarser mesh than a genuinely high-order accurate method.

Copyright and License

© 2025 Elsevier Ltd. All rights are reserved, including those for text and data mining, AI training, and similar technologies.

Acknowledgement

This work is not funded by anyone. The concerned work was started in April 2019. Some preliminary results, shown in Appendix B, were obtained in the first week of January 2021. Inviscid results were submitted to the APS conference (not attended due to lack of funding) in 2023 [87], and final results are presented in this paper. All the simulations are carried out on the author’s computer, a Mac Mini M1 with 8 GB RAM or the now-defunct MacBook Pro with 16 GB RAM. Some of the results took six months or more to compute. A.S. thanks his wife (sorry you could not even go to your father’s funeral) and his kid for their support despite the horrible situations they have gone through during these six years for the sake of “research and science”.

Data Availability

No data was used for the research described in the article.

Additional details

Dates

Accepted
2025-10-01
Available
2025-10-09
Available online
Available
2025-10-09
Version of record

Caltech Custom Metadata

Caltech groups
Division of Engineering and Applied Science (EAS)
Publication Status
Published