Efficient high-order Gradient-based Reconstruction for compressible flows

Amareshwara Sainadh Chamarthi

DOI: https://doi.org/10.1016/j.jcp.2023.112119

2023-05-01

Abstract:This paper extends the gradient-based reconstruction approach of Chamarthi \cite{chamarthi2023gradient} to genuine high-order accuracy for inviscid test cases involving smooth flows. A seventh-order accurate scheme is derived using the same stencil as of the explicit fourth-order scheme proposed in Ref. \cite{chamarthi2023gradient}, which also has low dissipation properties. The proposed method is seventh-order accurate under the assumption that the variables at the \textit{cell centres are point values}. A problem-independent discontinuity detector is used to obtain high-order accuracy. Accordingly, primitive or conservative variable reconstruction is performed around regions of discontinuities, whereas smooth solution regions apply flux reconstruction. The proposed approach can still share the derivatives between the inviscid and viscous fluxes, which is the main idea behind the gradient-based reconstruction. Several standard benchmark test cases are presented. The proposed method is more efficient than the seventh-order weighted compact nonlinear scheme (WCNS) for the test cases considered in this paper.

Fluid Dynamics,Numerical Analysis

What problem does this paper attempt to address?

### Problems the paper attempts to solve

The paper aims to achieve true high - order accuracy for smooth flows in inviscid flow test cases by extending the gradient - based reconstruction method proposed by Chamarthi. Specifically, the author develops a seventh - order accurate scheme, using the same stencil as the explicit fourth - order scheme proposed in reference [1] and having low - dissipation characteristics. The keys to this method are:

1. **High - order accuracy**: Under the assumption that the variables at the cell center are point values, this method can achieve seventh - order accuracy.
2. **Discontinuity detector**: Use a problem - independent discontinuity detector to achieve high - order accuracy. Reconstruct the primitive or conserved variables in the discontinuous region, and reconstruct the flux in the smooth region.
3. **Efficiency**: This method can share derivatives when dealing with inviscid and viscous fluxes, which is the main idea of the gradient - based reconstruction method.
4. **Performance advantage**: For the test cases considered in this paper, this method is more efficient than the seventh - order weighted compact nonlinear scheme (WCNS).

### Main contributions

1. **Gradient - based reconstruction of seventh - order accuracy**: A two - parameter gradient - based reconstruction method is proposed to achieve seventh - order accuracy. The fourth - order scheme in reference [1] can be regarded as a subset of the current method.
2. **Discontinuity detection of high - order accuracy**: For nonlinear test cases, a problem - independent discontinuity detector is used to obtain high - order accuracy. Reconstruct the primitive or conserved variables in the discontinuous region, and reconstruct the flux in the smooth region.
3. **Efficient time integration**: Use a third - order TVD Runge - Kutta scheme for time integration to ensure the stability and accuracy of the numerical solution.

### Mathematical models and methods

1. **Navier - Stokes equations**:
\[
\frac{\partial Q}{\partial t}+\frac{\partial F_c}{\partial x}+\frac{\partial G_c}{\partial y}+\frac{\partial F_v}{\partial x}+\frac{\partial G_v}{\partial y} = 0,
\]
where \(Q\) is the vector of conserved variables (\(\rho,\rho u,\rho v,E\)), \(\rho\) is the density, \(u\) and \(v\) are the velocity components in the \(x\) and \(y\) directions, the total energy per unit volume of the fluid is \(E=\rho\left(e+\frac{u^{2}+v^{2}}{2}\right)\), where \(e\) is the specific internal energy. The system is closed by the ideal gas equation of state:
\[
p = (\gamma - 1)\left(E-\rho\frac{u^{2}+v^{2}}{2}\right),
\]
where \(\gamma\) is the specific heat ratio, which is taken as 1.4 in this paper.

2. **Flux definitions**:
\[
F_c=\begin{pmatrix}
\rho u\\
\rho u^{2}+p\\
\rho uv\\
\rho uw\\
u(E + p)
\end{pmatrix},
G_c=\begin{pmatrix}
\rho v\\
\rho uv\\
\rho v^{2}+p\\
\rho vw\\
v(E + p)
\end{pmatrix},
F_v=\begin{pmatrix}
0\\
-\tau_{xx}\\
-\tau_{yx}\\
-\phi_x+q_x
\end{pmatrix},
G_v=\begin{pmatrix}
0\\
-\tau_{xy}\\
-\
\]