High-order implicit shock tracking boundary conditions for flows with parametrized shocks

Tianci Huang,Charles Naudet,Matthew J. Zahr

2023-04-22

Abstract:High-order implicit shock tracking (fitting) is a class of high-order, optimization-based numerical methods to approximate solutions of conservation laws with non-smooth features by aligning elements of the computational mesh with non-smooth features. This ensures the non-smooth features are perfectly represented by inter-element jumps and high-order basis functions approximate smooth regions of the solution without nonlinear stabilization, which leads to accurate approximations on traditionally coarse meshes. In this work, we introduce a robust implicit shock tracking framework specialized for problems with parameter-dependent lead shocks (i.e., shocks separating a farfield condition from the downstream flow), which commonly arise in high-speed aerodynamics and astrophysics applications. After a shock-aligned mesh is produced at one parameter configuration, all elements upstream of the lead shock are removed and the nodes on the lead shock are positioned for new parameter configurations using the implicit shock tracking solver. The proposed framework can be used for most many-query applications involving parametrized lead shocks such as optimization, uncertainty quantification, parameter sweeps, "what-if" scenarios, or parameter-based continuation. We demonstrate the robustness and flexibility of the framework using a one-dimensional space-time Riemann problem, and two- and three-dimensional supersonic and hypersonic benchmark problems.

Numerical Analysis,Optimization and Control

What problem does this paper attempt to address?

### What problem does this paper attempt to solve?

This paper aims to solve the high - precision numerical simulation problems involving parameterized leading shock waves (i.e., shock waves that separate far - field conditions and downstream flows), especially the challenges encountered in multi - query applications. Specifically, the paper proposes a new method based on the Implicit Shock Tracking (IST) framework to handle problems containing parameterized leading shock waves efficiently and robustly.

#### Main problem background

1. **Shock capture and grid alignment**
- In computational fluid dynamics (CFD), accurate shock capture is crucial for predicting key physical quantities such as aerodynamic heating.
- Traditional methods require the generation of very fine grids, and these grids must be aligned with the geometric shape of the shock wave, which results in huge computational costs and user input.

2. **Multi - query analysis**
- Multi - query applications (such as optimization, uncertainty quantification, parameter scanning, etc.) require a large number of simulations under different parameter configurations.
- When the shock wave position changes with parameters, traditional grid generation and simulation methods become very complex and time - consuming.

#### Solutions proposed in the paper

1. **Implicit Shock Tracking (IST) framework**
- Use a high - order implicit shock - tracking method to ensure that shock - aligned grids can be automatically generated and optimized.
- By solving an optimization problem, automatically adjust the positions of grid nodes so that the grid boundaries are aligned with the shock wave, thus avoiding the need for explicit identification of shock wave positions.

2. **Special - purpose framework for parameterized leading shock waves**
- Propose a new framework specifically for handling parameterized leading shock wave problems.
- After the initial grid generation, remove all elements upstream of the shock wave and only keep the grid downstream of the shock wave.
- For new parameter configurations, only optimize the node positions of the shock wave boundary, while other nodes are determined by smoothing partial differential equations (PDEs).

3. **Efficient multi - query analysis**
- Significantly improve computational efficiency by reducing degrees of freedom and optimizing initial guesses.
- Propose a method for using partially converged solutions to further accelerate the solution process in the parameter continuation method.

#### Application examples

- **One - dimensional space - time Riemann problem**: Demonstrates the robustness and flexibility of the framework.
- **Two - and three - dimensional supersonic and hypersonic benchmark problems**: Verifies the effectiveness of the method in actual complex flows.

Through these improvements, the paper provides an efficient and accurate method for handling flow problems containing parameterized leading shock waves, especially suitable for multi - query application scenarios.