Hostname: page-component-cd9895bd7-lnqnp Total loading time: 0 Render date: 2024-12-22T21:10:55.943Z Has data issue: false hasContentIssue false

Numerical modelling of supersonic boundary-layer receptivity to solid particulates – CORRIGENDUM

Published online by Cambridge University Press:  10 May 2024

Abstract

Type
Corrigendum
Copyright
© The Author(s), 2024. Published by Cambridge University Press

In our paper Chuvakhov et al. (Reference Chuvakhov, Fedorov and Obraz2019), the energy equation (2.6) is written for the gas temperature T as

(1) \begin{align} \rho \left[ {\dfrac{{\partial T}}{{\partial t}} + {u_j}\dfrac{{\partial T}}{{\partial {r_j}}}} \right] & = \dfrac{1}{{PrRe}}\dfrac{\partial }{{\partial {r_j}}}\left( {\mu \dfrac{{\partial T}}{{\partial {r_j}}}} \right) + (\gamma - 1)M_\infty ^2\left( {\dfrac{{\partial p}}{{\partial t}} + {u_j}\dfrac{{\partial p}}{{\partial {r_j}}}} \right) + \dfrac{{(\gamma - 1)M_\infty ^2}}{{Re}}\varPhi \nonumber\\ & \quad + \dfrac{{{R_p}}}{{PrRe}}{{\bar{Q}}_p}\delta (\boldsymbol{r} - {\boldsymbol{r}_p}) + \underline {R_p^2(\gamma - 1)M_\infty ^2({u_{pj}} - {u_j}){{\bar{F}}_{pj}}\delta (\boldsymbol{r} - {\boldsymbol{r}_p})} , \end{align}

where underlined is the particle-induced source term associated with the drag force. This form is convenient for analyses of receptivity and stability problems (Fedorov Reference Fedorov2013). However, in our numerical simulations the energy equation is used in the conservative form

(2) \begin{align} \begin{aligned} \dfrac{{\partial \rho e}}{{\partial t}} + \dfrac{{\partial \rho {u_j}H}}{{\partial {r_j}}} & = \dfrac{1}{{Re}}\dfrac{\partial }{{\partial {r_j}}}\left( {\dfrac{\mu }{{(\gamma - 1)M_\infty^2Pr}}\dfrac{{\partial T}}{{\partial {r_j}}} + {\tau_{ij}}{u_j}} \right)\\ & \quad + \dfrac{{{R_p}}}{{(\gamma - 1)M_\infty ^2PrRe}}{{\bar{Q}}_p}\delta (\boldsymbol{r} - {\boldsymbol{r}_p}) + \underline {R_p^2{u_{pj}}{{\bar{F}}_{pj}}\delta (\boldsymbol{r} - {\boldsymbol{r}_p})} , \end{aligned} \end{align}

where $e = p/(\rho (\gamma - 1)) + 0.5{u_i}{u_i}$ is gas total energy per unit mass, and $H = e + p/\rho $ is total enthalpy per unit mass, and the underlined term represents the power per unit volume which is produced by the particle passing through the gas. In our computations the underlined term of (2) was erroneously replaced by the underlined term of (1). After correcting this error, our computations showed that the disturbance amplitude reduced by 2.7 times; see figure 1 that is the correct version of figure 10 of Chuvakhov et al. (Reference Chuvakhov, Fedorov and Obraz2019). Note that the theoretical distributions were not changed.

Figure 1. Theoretical (lines with white symbols) and numerical (black symbols) distributions of the hump amplitude ${p^{\prime}_{w,max}}({x_{max}})$ for the collision points ${x_c} = 0.067$ (circles) and ${x_c} = 0.134$ (stars); ${F_s}$ – frequency parameter of the dominant wave at the observation point ${x_{max}}$.

Funding

The work completed to find these corrected findings was carried out in Moscow Institute of Physics and Technology under the support of Russian Science Foundation (project no. 19-79-10132).

References

Chuvakhov, P.V., Fedorov, A.V. & Obraz, A.O. 2019 Numerical modeling of supersonic boundary-layer receptivity to solid particulates. J. Fluid Mech. 859, 949971.CrossRefGoogle Scholar
Fedorov, A.V. 2013 Receptivity of a supersonic boundary layer to solid particulates. J. Fluid Mech. 737, 105131.CrossRefGoogle Scholar
Figure 0

Figure 1. Theoretical (lines with white symbols) and numerical (black symbols) distributions of the hump amplitude ${p^{\prime}_{w,max}}({x_{max}})$ for the collision points ${x_c} = 0.067$ (circles) and ${x_c} = 0.134$ (stars); ${F_s}$ – frequency parameter of the dominant wave at the observation point ${x_{max}}$.