Consistent Interface Capturing Adaptive Reconstruction Approach for Viscous Compressible Multicomponent Flows

TL;DR

This paper introduces a physically consistent numerical method for simulating viscous compressible multicomponent flows, improving interface sharpness and reducing dissipation errors.

Contribution

It develops a contact discontinuity detector and a hybrid reconstruction scheme combining THINC and MP/WENO, along with a central scheme for tangential velocities across interfaces.

Findings

Captured material interfaces sharply compared to existing methods.

Reduced dissipation errors near contact discontinuities.

Validated with several benchmark test cases showing improved performance.

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.