Inverse asymptotic treatment: capturing discontinuities in fluid flows via equation modification

Shahab Mirjalili,Søren Taverniers,Henry Collis,Morad Behandish,Ali Mani

2023-09-03

Abstract:A major challenge in developing accurate and robust numerical solutions to multi-physics problems is to correctly model evolving discontinuities in field quantities, which manifest themselves as interfaces between different phases in multi-phase flows, or as shock and contact discontinuities in compressible flows. When a quick response is required to rapidly emerging challenges, the complexity of bespoke discretization schemes impedes a swift transition from problem formulation to computation, which is exacerbated by the need to compose multiple interacting physics. We introduce "inverse asymptotic treatment" (IAT) as a unified framework for capturing discontinuities in fluid flows that enables building directly computable models based on off-the-shelf numerics. By capturing discontinuities through modifications at the level of the governing equations, IAT can seamlessly handle additional physics and thus enable novice end users to quickly obtain numerical results for various multi-physics scenarios. We outline IAT in the context of phase-field modeling of two-phase incompressible flows, and then demonstrate its generality by showing how localized artificial diffusivity (LAD) methods for single-phase compressible flows can be viewed as instances of IAT. Through the real-world example of a laminar hypersonic compression corner, we illustrate IAT's ability to, within just a few months, generate a directly computable model whose wall metrics predictions for sufficiently small corner angles come close to that of NASA's VULKAN-CFD solver. Finally, we propose a novel LAD approach via "reverse-engineered" PDE modifications, inspired by total variation diminishing (TVD) flux limiters, to eliminate the problem-dependent parameter tuning that plagues traditional LAD. We demonstrate that, when combined with second-order central differencing, it can robustly and accurately model compressible flows.

Computational Physics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is how to model the evolving discontinuities in fluid flow accurately and robustly in multi - physics problems. These discontinuities can be manifested as interfaces between different phases in multiphase flow, or as shock waves and contact discontinuities in (single - phase or multiphase) compressible flow. Traditional solutions rely on specialized discretization schemes to capture these discontinuities, but in the face of rapidly emerging new challenges, such as the need to design new hypersonic vehicles, their complexity hinders the rapid transition from problem definition to computation. Moreover, when it is necessary to combine multiple interacting physical processes, these specialized methods may unexpectedly fail due to the need to change and/or add control equations.

To this end, the paper introduces the "Inverse Asymptotic Treatment" (IAT) as a unified framework to capture discontinuities in fluid flow by modifying the control equations rather than relying on specialized discretization schemes. IAT makes it possible to construct directly computable models based on standard off - the - shelf numerical methods. By capturing discontinuities at the control - equation level, IAT can seamlessly handle additional physical processes, enabling novice users to quickly obtain numerical results in a variety of multi - physics scenarios. In addition, IAT also promotes the use of "multi - physics compilers", which can automate the conversion of the modified partial differential equations into numerical source code that can be directly read by popular computing frameworks such as OpenFOAM.

The paper demonstrates the versatility of IAT through examples of phase - field modeling of two - phase incompressible flow and the Localized Artificial Diffusivity (LAD) method in single - phase compressible flow. Specifically, the paper proposes a new LAD method, which is inspired by the Total Variation Diminishing (TVD) flux limiter, aiming to eliminate problem - related parameter adjustments in traditional LAD methods. Through classical numerical tests, the paper proves that this "limiter - inspired" LAD method, when combined with second - order central differences, can simulate compressible flow robustly and accurately.

Overall, the goal of the paper is to provide a systematic method to effectively address the discontinuity problem in fluid dynamics through modification of the control equations, thereby simplifying the numerical simulation of multi - physics problems and improving computational efficiency and accuracy.