Improvement of Mixing Function for Modified Upwinding Compact Scheme

This paper presents an improved mixing function for a modified upwinding compact scheme that effectively combines the high-order accuracy of compact schemes with the shock-capturing capability of WENO, enabling sharp shock resolution while preserving high resolution in smooth regions for applications like shock-boundary layer and shock-acoustic interactions.

The Gist

Imagine you are trying to paint a picture of the wind blowing around an airplane. You want your painting to be incredibly detailed, showing every tiny swirl and eddy of the air (turbulence), but you also need to capture the sudden, sharp "wall" of air that forms when the plane breaks the sound barrier (a shock wave).

This paper describes a new mathematical tool designed to solve a specific problem in computer simulations: How do you get both extreme detail and sharp edges without the picture getting blurry or shaky?

Here is the breakdown of the problem and the authors' solution, using simple analogies:

The Problem: Two Tools, One Flaw

The researchers were trying to combine two existing "paintbrushes" (mathematical schemes) that each had a major weakness:

1. The Compact Scheme (The Detail Brush): This tool is amazing at painting smooth, flowing air. It captures tiny details and small swirls with high precision. However, it is terrible at painting sharp edges. If you try to use it on a shock wave, the picture gets "wobbly" and oscillates, like a shaky hand trying to draw a straight line. It needs information from both sides of the line to work, which confuses it when a sudden wall appears.
2. The WENO Scheme (The Sharp Edge Brush): This tool is a master at painting sharp, sudden changes (shocks) without wobbling. However, it is too "heavy-handed." It smears out the smooth parts of the picture. If you use it for the tiny swirls of turbulence, it blurs them out, making the air look like thick soup instead of a fine mist.

The Goal: The authors wanted a "Super Brush" that uses the Detail Brush for smooth areas and the Sharp Edge Brush only when it hits a shock, switching between them seamlessly.

The Solution: The "Smart Mixer"

The authors created a Modified Upwinding Compact Scheme (MUCS). Think of this as a smart paintbrush that automatically knows when to switch tools.

1. The Shock Detector (The Eye): The new system has a built-in "eye" that scans the air. It looks for signs of a shock wave (a sudden, violent change). It's very sensitive and can spot even weak shocks.
2. The Control Function (The Hand): This is the most important part of their new work. In previous versions, the switch between the two brushes was like a light switch: ON (shock) or OFF (smooth). This caused jerky transitions.
   • The authors improved this by creating a dimmer switch. Instead of a sudden jump, the system gradually blends the two methods.
   • When the "eye" sees a shock, the "hand" slowly turns up the volume on the Sharp Edge Brush (WENO) and turns down the Detail Brush (Compact).
   • When the air is smooth, it does the opposite.
   • The Analogy: Imagine driving a car. Old methods were like slamming on the brakes when you saw a stop sign. The new method is like gently pressing the brake pedal as you approach the sign, then smoothly accelerating again once you pass it. This prevents the passengers (the data) from being thrown around.

The Results: A Sharper, Clearer Picture

The team tested this new "Smart Brush" on a simulation of air hitting a wall at high speed (2-D Euler equations).

• The Test: They compared their new method against using only the Sharp Edge Brush (WENO).
• The Outcome:
  • Sharpness: The new method captured the shock wave much sharper than the old Sharp Edge Brush. The "wall" of air was distinct and clear.
  • Detail: In the smooth areas, the new method kept the tiny swirls of turbulence visible. The old Sharp Edge Brush had blurred these out.
  • Stability: The new method didn't create the "wobbly" lines that the Detail Brush usually makes near shocks.
  • No "Magic Numbers": A major bonus is that this new system doesn't require the user to guess or tweak specific settings for different scenarios. It works automatically, which makes it much easier to use.

The Bottom Line

The authors successfully built a hybrid system that gets the best of both worlds. It keeps the high-resolution detail needed for studying turbulence while adding the ability to handle the violent, sharp shocks of supersonic flight. They proved it works well on a computer model of 2D airflow, showing that their new "dimmer switch" for mixing math tools is more efficient and accurate than previous methods.

Technical Summary

Technical Summary: Improvement of Mixing Function for Modified Upwinding Compact Scheme

Problem Statement
The paper addresses a fundamental trade-off in computational fluid dynamics (CFD) numerical schemes. The Compact Scheme (CS), specifically the Pade-type schemes introduced by Lele (1992), offers high-order accuracy and high resolution, making it highly effective for Direct Numerical Simulation (DNS) and Large Eddy Simulation (LES) of turbulence. However, standard compact schemes are non-dissipative and rely on both upstream and downstream points, rendering them unstable and incapable of capturing shocks without generating oscillations. Conversely, the Weighted Essentially Non-Oscillatory (WENO) scheme (Jiang & Shu, 1996) is robust for shock capturing but introduces excessive numerical dissipation. This dissipation degrades the resolution of small length scales and turbulence, making standard WENO unsuitable for high-fidelity DNS/LES applications. Previous attempts to hybridize these schemes often relied on complex mixing functions with case-specific adjustable coefficients.

Methodology
The authors propose a Modified Upwinding Compact Scheme (MUCS) that dynamically blends a high-order upwinding compact scheme with the WENO scheme based on local flow smoothness. The methodology involves three primary components:

1. Hybrid Flux Construction: The scheme utilizes a control function, $\Omega$, to mix the fluxes. The final numerical flux is a weighted combination: $(1-\Omega) \times \text{Compact Scheme} + \Omega \times \text{WENO}$.

   • When $\Omega = 0$, the scheme operates as a standard 6th-order compact scheme (or 7th-order upwinding compact).
   • When $\Omega = 1$, the scheme operates as a standard 5th-order WENO scheme.
   • This blending results in a tri-diagonal matrix system, preserving the compact nature of the derivative calculation while allowing for shock adaptation.

2. Shock Detection: The authors employ a shock detector capable of identifying shocks, discontinuities, and weak shocks. This detector calculates a ratio of coarse-grid to fine-grid truncation errors ($MR$) and a local left-right slope ratio ($LR$).

   • $MR$ is derived from the ratio of truncation errors on different grid levels.
   • $LR$ checks the smoothness of the function by comparing left and right slopes.

3. Optimization of the Control Function: A critical contribution of this work is the refinement of the mixing function $\Omega$.

   • The original control function was defined as $\Omega = \min[1.0, 0.8 / MR \times LR]$. The authors noted this sometimes yielded values too small for shock locations, leading to excessive smearing.
   • The new control function introduces a square root operation and a local averaging mechanism to ensure smoothness and appropriate weighting:
     • An intermediate variable $A(i, h)$ is calculated as $\min[1.0, 4.0 / (MR \times LR)]$.
     • The square root is applied to $A(i, h)$ to "recover" the value closer to 1.0 in shock regions, preventing the mixing function from becoming too small.
     • The final weight $\Omega$ is the average of the square-rooted values of the current point and its two immediate neighbors. This averaging reduces misjudgment and ensures the control function remains smooth (avoiding the use of a hard switch function).

Key Contributions

• Development of MUCS: The paper presents a modified upwinding compact scheme that successfully captures shocks sharply while maintaining high-order accuracy and resolution in smooth regions.
• Refined Mixing Function: The authors introduce a new, smoother control function that eliminates the need for case-related adjustable coefficients. This function ensures the scheme remains upwinding near shocks, provides sufficient dissipation before and after the shock, and maintains high accuracy in smooth areas.
• Efficiency and Robustness: The new scheme is designed to be efficient, avoiding the complexity of previous hybrid schemes that required tuning parameters for specific cases.

Computational Results
The scheme was validated using 2-D Euler equations with an incident shock problem (Inflow Mach number 2, attach angle $35.241^\circ$). Comparisons were made against exact solutions and pure 5th-order WENO results on grids ranging from coarse to fine ($129 \times 129$).

• Shock Capturing: The MUCS captured the shock sharper than the pure WENO scheme across all grid resolutions.
• Oscillation Control: While the pure compact scheme exhibited oscillations on fine grids, the MUCS did not show serious oscillations, performing better than WENO after the second shock.
• Resolution: The results demonstrated that MUCS is comparable to the exact solution. Crucially, the pure WENO scheme exhibited significant smearing of the flow, particularly in pressure distributions and at specific contour levels (e.g., $K=30$). The MUCS avoided this smearing, indicating superior capability for resolving small length scales.
• Detector Performance: The shock detector accurately identified shock locations, and the new control function successfully transitioned the scheme from WENO-dominated (near shocks) to UCS-dominated (in smooth regions).

Significance and Claims
The authors claim that the MUCS, equipped with the new shock detector and optimized mixing function, is ready for application to both 2-D and 3-D Euler and Navier-Stokes equations. The primary significance lies in its ability to achieve sharp shock capturing without sacrificing the high resolution required for turbulence simulation. The paper emphasizes that the scheme contains no case-related parameters, making it a robust and generalizable tool for high-fidelity simulations involving shock-boundary layer and shock-acoustic interactions. The work is presented as a continuation of previous efforts to modify the control function, with numerical results confirming its success for 2-D Euler equations.