TY - JOUR

T1 - A true-direction reconstruction of the quiet direct simulation method for inviscid gas flows

AU - Lin, Y. J.

AU - Smith, M. R.

AU - Kuo, F. A.

AU - Cave, H. M.

AU - Huang, J. C.

AU - Wu, Jong-Shinn

PY - 2013/11

Y1 - 2013/11

N2 - In this paper, a true-direction flux reconstruction of the second-order quiet direct simulation (QDS-2N) Smith et al. (2009) [3] as an equivalent Euler equation solver, called QDS-N2, is proposed. Because of the true-directional nature of QDS, where volume-to-volume (true-direction) fluxes are computed, as opposed to fluxes at cell interfaces as employed by traditional finite volume schemes, a volumetric reconstruction is required to reach a totally true-direction scheme. The conserved quantities are permitted to vary (according to a polynomial expression) across all simulated dimensions. Prior to the flux computation, QDS particles are introduced using properties based on weighted moments taken over the polynomial reconstruction of the conserved quantity fields. The resulting flux expressions are shown to exactly reproduce the existing second-order extension for a one-dimensional flow, while providing a means for true multi-dimensional reconstruction. The new reconstruction is demonstrated in several verification studies. These include a shock-bubble interaction problem, an Euler-four-shock interaction problem, and the advection of a vortical disturbance. These results are presented, and the increased computation time and the effect of higher-order extension are discussed in this paper. The results show that the proposed multi-dimensional reconstruction provides a significant increase in the accuracy of the solution. We show that, despite the increase in the computational expense, the computational speed of the proposed QDS-N2 method is several times higher than that of the previously proposed QDS-2N scheme for a fixed degree of numerical accuracy, at least, for the test problem of the advection of vertical disturbances.

AB - In this paper, a true-direction flux reconstruction of the second-order quiet direct simulation (QDS-2N) Smith et al. (2009) [3] as an equivalent Euler equation solver, called QDS-N2, is proposed. Because of the true-directional nature of QDS, where volume-to-volume (true-direction) fluxes are computed, as opposed to fluxes at cell interfaces as employed by traditional finite volume schemes, a volumetric reconstruction is required to reach a totally true-direction scheme. The conserved quantities are permitted to vary (according to a polynomial expression) across all simulated dimensions. Prior to the flux computation, QDS particles are introduced using properties based on weighted moments taken over the polynomial reconstruction of the conserved quantity fields. The resulting flux expressions are shown to exactly reproduce the existing second-order extension for a one-dimensional flow, while providing a means for true multi-dimensional reconstruction. The new reconstruction is demonstrated in several verification studies. These include a shock-bubble interaction problem, an Euler-four-shock interaction problem, and the advection of a vortical disturbance. These results are presented, and the increased computation time and the effect of higher-order extension are discussed in this paper. The results show that the proposed multi-dimensional reconstruction provides a significant increase in the accuracy of the solution. We show that, despite the increase in the computational expense, the computational speed of the proposed QDS-N2 method is several times higher than that of the previously proposed QDS-2N scheme for a fixed degree of numerical accuracy, at least, for the test problem of the advection of vertical disturbances.

KW - CFD

KW - Euler equations

KW - Kinetic theory of gases

KW - QDS

KW - QDSMC

KW - TDEFM

UR - http://www.scopus.com/inward/record.url?scp=84883247780&partnerID=8YFLogxK

U2 - 10.1016/j.cpc.2013.05.007

DO - 10.1016/j.cpc.2013.05.007

M3 - Article

AN - SCOPUS:84883247780

SN - 0010-4655

VL - 184

SP - 2378

EP - 2390

JO - Computer Physics Communications

JF - Computer Physics Communications

IS - 11

ER -