To complement our previous analysis of interactions of large-scale turbulence with strong detonations, the corresponding theory of interactions of small-scale turbulence is presented here. Focusing most directly on the results of greatest interest, the ultimate long-time effects of high-frequency vortical and entropic disturbances on the burnt-gas flow, a normal-mode analysis is selected here, rather than the Laplace-transform approach. The interaction of the planar detonation with a monochromatic pattern of perturbations is addressed first, and then a Fourier superposition for two-dimensional and three-dimensional isotropic turbulent fields is employed to provide integral formulas for the amplification of the kinetic energy, enstrophy, and density fluctuations. Effects of the propagation Mach number and of the chemical heat release and the chemical reaction rate are identified, as well as the similarities and differences from the previous result for the thin-detonation (fast-reaction) limit.
In previous work,1 we analyzed the interaction of a strong detonation with inhomogeneous density fields. The motivation was to improve the descriptions of the influence of compressible turbulence on detonation propagation and, in particular, to determine how passage of a planar detonation modifies the turbulence. In that respect, the work complemented the analyses of Jackson et al.,2,3 who addressed interactions of detonations with constant-density vorticity fields. It is well known4 that, excluding acoustic perturbations, which propagate with respect to the fluid, inhomogeneities that travel with the fluid can be decomposed into vortical and entropic components. The earlier work2,3 had explored the influences of the vortical component, but at high Mach numbers those influences are dominated by influences of the entropic component,5 which therefore was in need of study. The previous investigation1 provided the additional information that was needed for strong detonations that could be treated as discontinuities.
General reasons for interest in interactions between detonations and turbulence have been discussed and referenced earlier,1 and, since the present paper relies on knowledge of the previous paper, those reasons will not be repeated here. One motivation is improvement in the performance of detonative propulsion devices for hypersonic aircraft, where limited times for fuel-air mixing can pose problems that could be attacked through shock-wave-enhanced mixing.6,7 In such situations, detonation thicknesses may be larger than the most important representative turbulence scales. The earlier analysis1 then become inaccurate, and turbulence effects on the evolution of the finite-rate heat release in the interior of the detonation become relevant. The present investigation offers a step towards removing this deficiency.
It is worth observing that numerical studies by Massa et al.8,9 have shown enhanced interactions between turbulence and detonations associated with there being comparable sizes of the inhomogeneities and the unperturbed reaction zones. This increases the interest in pursuing analytical investigations in which turbulence scales are not large compared with reaction-zone thicknesses, for example, for testing scaling hypotheses and to determine whether effects of systematic variations of parameters can be derived. Such investigations become quite difficult when the reaction zones and inhomogeneities are of comparable size, but they can be performed accurately when the dimensions of all disturbances are small compared with reaction-zone thicknesses. This is the limit to be addressed here. Inferences concerning behavior at intermediate scales may be derived from resulting functional dependences and comparisons with previous thin-detonation results. Since neither vortical nor entropic fluctuations have yet been considered in this limit, both will be analyzed in the present work.
The analysis is a perturbation of the steady, planar ZND detonation structure, a shock followed by an inviscid reaction zone,10–12 for an ideal gas with a constant specific heat at constant pressure. The formulation will be entirely in terms of nondimensional variables. For simplicity, we assume that the rate of energy release per unit volume behind the lead shock $\dot{\bar{q}}$
is uniformly distributed over the thickness of the reaction region, which ends once the fuel is depleted at a finite distance from the shock. This thickness $(\bar{l})$
is taken as the unit of length in the analysis, and the sound velocity $\bar{a}_{\scriptscriptstyle {\rm N}}$
at the Neumann state (just behind of the shock) is taken as the unit of velocity, overbars identifying dimensional quantities. All velocities, both streamwise $\bar{u}$
and transverse $\bar{v}$
, as well as sound speeds $\bar{a}$
are nondimensionalized with respect to this same sound velocity, and the conditions at the Neumann state (identified by the subscript N) are also used to normalize the other variables, so that the density $\bar{\rho }$
is normalized by $\bar{\rho }_{\scriptscriptstyle {\rm N}}$
, the pressure $\bar{p}$
with $\bar{\rho }_{\scriptscriptstyle {\rm N}} \bar{a}_{\scriptscriptstyle {\rm N}}^2$
, and the temperature T by $\bar{T}_{\scriptscriptstyle {\rm N}}$
. The subscripts o and b will identify conditions in the fresh mixture and in the burnt gas behind the detonation, respectively, and subscripts l and s will distinguish quantities in the laboratory reference frame in which the fresh mixture is at rest and in the frame moving with the unperturbed detonation, respectively. The three nondimensional parameters that appear in the problem are then the propagation Mach number $M_o=\bar{D}/\bar{a}_o$
of the undisturbed detonation, where $\bar{D}$
is the propagation velocity, the (constant) ratio of specific heats γ, and the total energy released per unit mass of the mixture multiplied by $(\gamma ^2-1)/(2 \bar{a}_o^2)$
, denoted by q. That is, $q=\bar{q}(\gamma ^2-1)/(2 \bar{a}_o^2)$
, where $\bar{q}$
is the heat released per unit mass of mixture, which is related to the heat-release rate by $\dot{\bar{q}}=\bar{\rho }_{\scriptscriptstyle {\rm N}}\bar{u}_{\scriptscriptstyle {\rm N}}\bar{q}/\bar{l}$
, with $\bar{u}_{\scriptscriptstyle {\rm N}}$
representing the Neumann streamwise velocity relative to the shock. The nondimensional streamwise and transverse coordinates are x and y, respectively.
Figure 1, in which D denotes the nondimensional propagation velocity of the detonation, illustrates the particle paths and velocity profiles in the laboratory frame of reference. The upstream conditions can be expressed in units of their corresponding Neumann values by means of the Rankine-Hugoniot equations, which lead to
The Neumann Mach number defined as the nondimensional velocity at which the particles leave the shock front is
For strong detonations, Mo is appreciably greater than its Chapman-Jouget value,
where $\bar{D}_{cj}$
is the Chapman-Jouget detonation propagation velocity. Behind the shock front, the combustion process takes place, and the thermodynamic quantities vary in accordance with the type of heat-deposition involved. To calculate the steady-state profiles of the thermodynamic quantities, it is convenient to write the conservation equations in the shock-fixed reference frame xs = Dt − xl, where the fluid particles travel with velocity us = D − ul. Here $t=\bar{t}\ \bar{l}/\bar{a}_{\scriptscriptstyle {\rm N}}$
. After some straightforward algebra, for a constant heat-release rate, the undisturbed pressure profile as a function of the distance from the shock xs = [0, 1] is found to be
which is related to the velocity profile according to $p=\gamma ^{-1}+M_{\scriptscriptstyle {\rm N}}^2-u_s M_{\scriptscriptstyle {\rm N}}$
and to the density profile, since $\rho =M_{\scriptscriptstyle {\rm N}}/u_s$
. The temperature profile is obtained directly from the perfect-gas equation T = γp/ρ, and the dimensionless local speed of sound is then given by a2 = T. At the end of the reaction zone (xs = 1), there is a weak discontinuity that separates the burning and the burnt flow regions; beyond this point, the unperturbed fluid properties remain constant.
FIG. 1.
Schematic fluid particle paths in the laboratory reference frame.
FIG. 1.
Schematic fluid particle paths in the laboratory reference frame.
Close modal
Complex notation is employed in the normal-mode analysis, with the physical values represented by the real parts. The upstream vortical velocity field is defined in terms of the nondimensional longitudinal $u_o^{\prime }$
and transverse $v_o^{\prime }$
perturbations in the laboratory frame of reference as
where the divergence-free condition $k_x u_o^{\prime }=-k_y v_o^{\prime }$
has been applied. The factor εr represents the amplitude of the upstream velocity disturbances, and it is assumed to be much smaller than unity to remain within the limits of linear theory. Here the wave-number vector $\vec{k}$
, with x and y components kx and ky, is scaled with the inverse of the detonation thickness $1/\bar{l}$
, in accord with our nondimensionalization, so that, in the present analysis, kx and ky are large parameters of expansion. Similarly, the density perturbation field is
where εe is the nondimensional amplitude, assumed to be small, like εr. Since the linear responses are independent, no phase angle is imposed between the rotational (subscript r) and entropic (subscript e) perturbations.
Incoming waves from behind the detonation are excluded. When the lead shock encounters the perturbations in (5) and (6), it develops a time-dependent response proportional to $e^{i \omega _s t}$
, where the nondimensional frequency is defined as $\omega _s= M_{\scriptscriptstyle {\rm N}}k_x/\rho _o$
. The linearized Rankine-Hugoniot relationships then yield for the coefficients of the factor $e^{i(\omega _s t+ k_y y)}$
in the nondimensional perturbations
where ξ
s =
ky[
xl, s(
t) −
Dt] is the amplitude of the perturbation of the shock position. An equation will be derived for
$p_s^{\prime }$
which determines pressure perturbations just behind the shock.
For a given harmonic excitation, such as that considered in (5) and (6), the linear disturbances following the shock front also must be harmonic. At leading order, the perturbations are adiabatic so an acoustic wave form describes the post-shock perturbation field, the acoustic frequency and wave number being related to the shock frequency by the expression $\omega _{a,\scriptscriptstyle{\rm N}}=\omega _s-M_{\scriptscriptstyle {\rm N}}k_{a,\scriptscriptstyle {\rm N}}$
. Evaluating the linear reactive Euler equations at the shock front then yields
where the dimensionless acoustic wavenumbers
$k_{a,\scriptscriptstyle {\rm N}}$
and
ky come from the streamwise and lateral derivatives, respectively, and the dimensionless frequency ω
s is obtained from the material derivative following the shock front. At leading order,
$\dot{q}=0$
, the well-known adiabatic dispersion relationship is recovered,
$\omega _{a,\scriptscriptstyle {\rm N}}^2=k_{a,\scriptscriptstyle {\rm N}}^2+k_y^2$
. The dimensionless acoustic wavenumber
$k_{a,\scriptscriptstyle {\rm N}}$
and the associated frequency are
The effect of the chemical reaction enters through the dimensionless rate of heat release
which is scaled with use made of the detonation thickness $\bar{l}$
and the unperturbed enthalpy flux at the Neumann state $\bar{\rho }_{\scriptscriptstyle {\rm N}}\bar{u}_{\scriptscriptstyle {\rm N}}\bar{a}_{\scriptscriptstyle {\rm N}}^2/(\gamma -1)$
.
There can be stable acoustic radiation right behind the shock wave if $\omega _s>(1- M_{\scriptscriptstyle {\rm N}})^{1/2}k_y$
, but if $\omega _s<(1- M_{\scriptscriptstyle {\rm N}})^{1/2}k_y$
, then these pressure perturbations are evanescent (waves whose amplitudes exponentially decay behind the shock front).13 The latter condition corresponds to the square roots in (9) becoming imaginary, contributing exponential rather than oscillatory behavior, which must decay to satisfy the boundary conditions. The shock oscillation frequency is conveniently described in terms of a dimensionless frequency ζ,
θ being the angle between the incident wavenumber $\vec{k}$
and the direction of propagation of the unperturbed detonation. It is useful to define ζ in this manner here because acoustic radiation can then occur for ζ ⩾ 1 but the acoustics is evanescent for ζ < 1. The wavelength of the downstream disturbances are long for ζ < 1 and short ζ ⩾ 1, and therefore solutions for ζ ⩾ 1, the high-frequency solutions that appear for small angles θ, are called short-wavelength solutions, while those corresponding to the low-frequency, large-θ range ζ < 1 are called long-wavelength solutions. By introducing an expansion up to terms of order $k_y^{-1}$
, Eqs. (7) and (8) can be combined to yield
plus terms of order $\epsilon k_y^{-2}$
, where
Equation (12) is the generalization of the well-known expression for the response of the pressure behind the shock wave, accounting for the exothermic finite-rate chemistry in the heat-release zone of the detonation that follows the shock. The present contribution analyzes the influences of the terms involving $a_{\scriptscriptstyle {\rm R}}$
and $b_{\scriptscriptstyle {\rm R}}$
which describe this effect. It is worthwhile to note that these influences depend only on the local effect caused by the rate of heat release $\dot{\bar{q}}$
right behind the shock. The profile of the heat-release rate downstream becomes irrelevant. These results thus will be general in the sense that they are independent of the particular model that produced (4). All that is required is that there be no regions of rapid heat-release rate. With the length $\bar{l}$
defined in terms of the average rate of heat release, the parameter $k_y^{-1}$
must be small.
Although the formulation leading to (12) is then independent of the heat-release processes occurring downstream, those processes do affect the disturbances that occur within the thick detonation and downstream therefrom. As is well known from the earlier investigations, those disturbances can be considered to be of two general types, namely, acoustic fluctuations that propagate downstream with respect to the fluid locally and fluctuations that remain fixed with respect to the fluid particles. In the present notation, the former are proportional to $e^{i (\omega _s t + k_a x_s + k_y y)}$
, and the latter are proportional to $e^{i [\omega _s t - (\rho /\rho _o) k_x x_s + k_y y]}$
. Here
for the acoustic component, and the corresponding acoustic frequency experienced by a fluid element moving with the local fluid speed is
The results to be given below pertain to the coefficients of these disturbances, in particular, as the fluctuations emerge into the burnt gas. There are contributions to theses modification associated with passage through the heat-release region at order ε as well as order $\epsilon k_y^{-1}$
. The former, however, turn out to be of the same functional form in the burnt gas as those derived in previous studies of thin detonations, but with coefficients corresponding instead to passage through the lead shock only. Unlike the thin detonation, the thick detonation does not modify the amplitudes from those of the lead shock. Aside from this difference, the results are identical in the two limits. Attention therefore will be focused on the fluctuations of order $\epsilon k_y^{-1}$
.
In terms of the expansions in (12), the results are quite different, depending on whether the acoustics is radiating (ζ ⩾ 1) or evanescent (ζ < 1). For radiating conditions, there is no contribution at order $k_y^{-1}$
, and so the first correction that involves the heat release q directly is of order $k_y^{-2}$
, beyond the range of the current analysis. The perturbations to be derived here that explicitly involve q therefore are restricted to the evanescent conditions that occur at the larger angles θ that produce ζ < 1. Otherwise the results from (12) are the same as those for nonreacting shocks.
The aforementioned indirect effect of the heat release in the detonations on the acoustic perturbations at leading order do, however, influence the final radiation condition. As a consequence of the variation of the quantities M and a with the distance behind the shock, the condition for the detonation to be radiating at its downstream boundary is ζ > ζb, where
which is plotted in Fig. 2 as a function of f − 1 (where $f=M_o^2/M_{cj}^2$
is the overdrive factor), showing that it can be greater or less than unity. For sufficiently weak overdrive, ζb is less than unity, and so, when ζ < 1, radiating conditions are encountered at some point within the detonation, but the exponential decay of the evanescent wave from the shock before that will render the amplitude of the radiated wave negligible small for thick detonations, whence it is not considered here. On the other hand, especially at the higher values of the heat release parameter q, above a critical overdrive ζb becomes greater than unity, in which case an evanescent shock will always correspond to an evanescent detonation. The relevant parameter therefore in every case is ζ, and the perturbations calculated here are all for downstream conditions that are effectively evanescent.
FIG. 2.
Frequency ζb for γ = 1.2 as a function of Mo for different heat releases q = 0.1, 1, 10.
FIG. 2.
Frequency ζb for γ = 1.2 as a function of Mo for different heat releases q = 0.1, 1, 10.
Close modal
Because of the dominance of the lead shock for thick detonations, the amplitudes of the perturbation at order ε are the same as those derived previously for shock-vorticity interactions (εe = 0)13,14 and for shock-density interactions (εr = 0).15 Since there are no perturbations at order $\epsilon k_y^{-1}$
for radiating acoustics (ζ ⩾ 1) the exothermicity does not affect the amplitudes of the radiated acoustic fields at this order. The amplitudes of the evanescent fields are, however, affected. To describe such effects, the acoustic pressure fluctuations in the burnt gas at this order can be represented as $\epsilon \dot{q} k_y^{-1} C_p e^{i (\omega _s t + k_{a,b} x_s + k_y y)}$
, where ka, b, having a positive imaginary part in this case, is given by (14) with M = Mb and a = ab. The factor $\dot{q}$
is included here because these corrections are proportional to $\dot{q}$
. Since overall phase angles are irrelevant for subsequent statistical averages, the complex coefficients Cp will not be given here; instead only Ap = ±|Cp| will be given, with the positive sign selected for contributions that increase the amplitude of the perturbation of order ε and the negative sign for perturbations that decrease that amplitude. With this convention, it is found that
where the subscripts r and e identify the rotational and entropic contributions to ε, respectively. The factors Bpr and Bpe refer to the pressure disturbances generated by the inert shock (see the Appendix for details) and are given by
In a similar manner, if the contributions to the density and vorticity fluctuations in the burnt gas at order $\epsilon k_y^{-1}$
are denoted by $\epsilon \dot{q} k_y^{-1} C_{\rho } e^{i [\omega _s t + \frac{\rho _b}{\rho _o} x_s + k_y y]}$
and $\epsilon \dot{q} k_y^{-1} C_{\omega } e^{i [\omega _s t + \frac{\rho _b}{\rho _o} x_s + k_y y]}$
, respectively, then results are given for Aρ = ±|Cρ| and Aω = ±|Cω|, the signs similarly being selected according to the influence on the perturbations of order ε. It is then found that
and, with the Ω’s defined in the Appendix,
where
Bρe,
Bωr, and
Bωe are the order ε contributions for the long-wavelength entropic and vorticity perturbations (see the Appendix for details) and are given by
The velocity field is also modified at order $\epsilon k_y^{-1}$
as a consequence of the associated vorticity perturbations. The corresponding A(u, v)r and A(u, v)e are provided in (A10) and (A12), respectively.
In viewing the corrections proportional to $\epsilon \dot{q} k_y^{-1}$
, it is helpful to plot the results in Eqs. (20) and (21) for selected values of Mo and γ, because the formulas are too complicated for simple conclusions to be readily drawn directly from them. Since the results vary strongly with θ as the critical orientation is approached, their characteristics are seen most clearly on a log scale of 1 − ζ. In Fig. 3, such plots for the coefficients Aρ and Aω therefore are shown for the representative values Mo = 5 and γ = 1.2. The general characteristics of the results are the same for other values of these parameters.
FIG. 3.
Exothermic contribution for the density Aρ and vorticity Aω fields as a function of 1 − ζ for Mo = 5 and γ = 1.2.
FIG. 3.
Exothermic contribution for the density Aρ and vorticity Aω fields as a function of 1 − ζ for Mo = 5 and γ = 1.2.
Close modal
All the curves have certain attributes in common. For example, all of the corrections approach zero as 1 − ζ approaches zero and unity, which, according to (11), correspond, respectively, to the critical angle θ and to the detonation encountering the perturbation field edge-on. It is understandable that an edge-on encounter will not have any effect, and since the corrections are of higher order for the angles of encounter between normal and the critical angle (that is, at the higher frequencies), considerations of continuity suggest that the approach to zero at the critical angle is reasonable. In addition, all of the curves achieve both negative and positive values, although, for purely entropic initial perturbations, the range and extent of the negative segment are very small. The curves, in general, exhibit oscillations, in that they pass through zero once or twice. This complicates the task of drawing general conclusions.
Since the logarithmic scale in Fig. 3 accentuates the range as the critical angle is approached, in determining the dominant overall effect, the scale variation must be considered. With this in mind, it may be inferred from the curves that the corrections are mainly positive for entropic upstream disturbances and mainly negative for rotational upstream disturbances. It was found that the heat release in thin detonations tends to reduce fluctuations for entropic upstream fluctuations.1 Although the curves in the earlier paper2 at first glance suggest otherwise, in fact, with the current scaling, the same is true for rotational fluctuations.2 Hence, the present thick-detonation results are in opposite directions from the thin-detonation results for entropic upstream disturbances but in the same direction for rotational disturbances. This suggests that entropic disturbances, which tend to be dominant at high Mach numbers, may experience enhanced responses at intermediate sizes.
Statistical averages for turbulent flows are generated from the preceding results by integrating the square amplitudes of the perturbations over the spectrum of wave numbers. As in previous work,1,3 the upstream flow is assumed to be homogeneous and isotropic so that the wave-number vector $\vec{k}$
is uniformly distributed over the unit sphere (or around the unit semicircle in two dimensions). Therefore, the analyses consists of a superposition of linear perturbations whose amplitudes are functions only of the wave-number amplitude $|\vec{k}|$
.16 Averaged quantities of interest include the kinetic energy K, the enstrophy W, and the density G.
When considering the interaction with a two-dimensional isotropic vorticity field, the velocity field is represented by divergence-free velocity disturbances $\mathbf {v}_o^{\prime } = \epsilon _k a_o (\sin \theta ,\cos \theta )$
, whose associated wave-number vector is $\vec{k}=k(-\cos \theta ,\sin \theta )$
. The mean square of the longitudinal and transverse velocity perturbations are3,14,17,18
Thanks to isotropy upstream, the spectrum contribution $\int _0^{\infty }\epsilon _k(k)^2k{\rm d}k$
is factored out. It then becomes possible to eliminate this factor in expressing results as ratios of post-detonation intensities to upstream intensities. The three-dimensional problem is conveniently formulated in spherical polar coordinates, so the upstream velocity field is $\mathbf {v}_o^{\prime }=a_o\epsilon _k(\sin \theta ,\cos \theta \cos \varphi ,\cos \theta \sin \varphi )$
and the associated wave-number vector is $\vec{k}=k(-\sin \theta ,\cos \theta \sin \varphi ,\cos \theta \cos \varphi )$
. The upstream mean-square disturbances are expressed as3,14,17,18
in that case. If there are only density disturbances ahead of the detonation wave, then the statistical averages under the same isotropy assumptions are1,15,19
where εe(k) represents the normalized internal-energy spectrum. In all cases results are given as ratios of final to initial normalized intensities, independent of the particular spectral distributions.
Consistent with the monochromatic perturbations to the inert lead-shock results being proportional to $\dot{q} k_y^{-1}$
with respect to the order-ε base-case perturbations, for isotropic turbulence the corresponding relative perturbations to the statistical averages will be of order $\dot{q} |k|^{-1}$
. The interest here therefore focuses on corrections of order $\epsilon \dot{q} |k|^{-1}$
to the order-ε base-case predictions of the changes in the statistical averages. With $K_{\scriptscriptstyle {\rm B}}$
denoting the order-unity ratio of the normalized post-detonation turbulent kinetic energy (normalized by the sound speed at the Neumann state) to the normalized upstream turbulent kinetic energy (normalized by the sound speed in the fresh mixture) and $K_{\scriptscriptstyle {\rm A}}$
denoting the order-unity coefficient of the correction to this ratio of order $\dot{q} |k|^{-1}$
, interest then rests on the fractional correction $\dot{q} |k|^{-1}\Delta _{\scriptscriptstyle {\rm K}}$
, where $\Delta _{\scriptscriptstyle {\rm K}}=K_{\scriptscriptstyle {\rm A}}/K_{\scriptscriptstyle {\rm B}}$
is a quantity of order unity that is independent of $\dot{q}$
and depends only on Mo and γ. While this definition is useful for rotational turbulent fluctuations, it clearly cannot be applied if the upstream fluctuations are purely entropic. In that case, as in earlier work,1 $K_{\scriptscriptstyle {\rm B}}$
is defined as a normalized post-detonation turbulent kinetic energy divided by the mean-square fresh-mixture density fluctuations normalized by the average density of the fresh mixture. However, instead of employing the burnt-gas sound speed for normalizing the post-detonation turbulent kinetic energy, as earlier,1 now the Neumann-state sound speed is used for that purpose, because otherwise and additional dependence on $\dot{q}$
arises. With this revised normalization and $K_{\scriptscriptstyle {\rm B}}$
again denoting the order-unity ratio for the post-detonation turbulent kinetic energy and $K_{\scriptscriptstyle {\rm A}}$
the order-unity coefficient of the correction of order $\dot{q} |k|^{-1}$
, the order-unity measure $\Delta _{\scriptscriptstyle {\rm K}}=K_{\scriptscriptstyle {\rm A}}/K_{\scriptscriptstyle {\rm B}}$
is to be exhibited here. The additional subscripts r and e will distinguish whether rotational or entropic fluctuations are considered, as before.
There are similar breakdowns for enstrophy and density. Here, too, to avoid introducing an additional dependence on $\dot{q}$
, the post-detonation enstrophy and density fluctuations are normalized by properties at the Neumann state, while the upstream perturbations are normalized, as always, by the upstream mean properties. In this way, the associated order-unity corrections are $\Delta _{\scriptscriptstyle {\rm W}}=W_{\scriptscriptstyle {\rm A}}/W_{\scriptscriptstyle {\rm B}}$
and $\Delta _{\scriptscriptstyle {\rm G}}=G_{\scriptscriptstyle {\rm A}}/G_{\scriptscriptstyle {\rm B}}$
. Even though the initial turbulence is isotropic, the interaction with the shock induces anisotropy, and there may be interested separately in longitudinal velocity fluctuations (velocities in the direction of detonation propagation) and transverse velocity fluctuations (velocities transverse to the propagation direction). The longitudinal and transverse contributions to the kinetic-energy ratio for K are denoted by L and T, respectively, and results for $\Delta _{\scriptscriptstyle {\rm L}}=L_{\scriptscriptstyle {\rm A}}/L_{\scriptscriptstyle {\rm B}}$
and for $\Delta _{\scriptscriptstyle {\rm T}}=T_{\scriptscriptstyle {\rm A}}/T_{\scriptscriptstyle {\rm B}}$
also will be given. Specific definitions appear in Subsections VI A–VI C.
For rotational upstream disturbances, following Refs. 3 and 13 and 14, amplification factors for turbulent kinetic energy are defined as
for two-dimensional and three-dimensional disturbances, respectively, where the functions $|u_b^{\prime }|$
and $|v_b^{\prime }|$
refer to the longitudinal and transverse velocity perturbation amplitudes in the burnt gas, and the normalized probability-density distributions are
for the two-dimensional and three-dimensional cases, respectively. The factors $M_{\scriptscriptstyle {\rm N}}^2+\zeta ^2\rho _o^2\left(1-M_{\scriptscriptstyle {\rm N}}^2\right)$
, which appear throughout, arise through Eq. (11) in converting integrations over θ to integrations over ζ. The breakdown $K= K_{\scriptscriptstyle {\rm B}}+ \dot{q} |k|^{-1}K_{\scriptscriptstyle {\rm A}}$
may also be divided into the longitudinal and transverse components, according to K = (1/2)(L + T) in two dimensions and K = (1/3)(L + 2T) in three dimensions. The three-dimensional case can be simplified so as to deal with an equivalent two-dimensional problem by choosing a reference frame that leaves one of the velocity contributions invariant. The base-component terms, which account for the rotational and acoustic contributions at order ε, are provided in the Appendix in (A13) and (A14) and also in the literature.13,14 The corresponding quantities at order $\epsilon \dot{q} |k|^{-1}$
are
for both two and three dimensions, and
for two dimensions and three dimensions, respectively. The factor $\left[1+\zeta ^2\rho _o^2\left(M_{\scriptscriptstyle {\rm N}}^{-2}-1\right)\right]^{1/2}$
, which does not appear for the base-contribution factors, arises from the equality |k|/ky = (sin θ)−1. The amplitudes of the base velocity fluctuations Bur and Bvr that appear here are provided in (A9a), and the corrections Aur and Avr are given in (A10).
As explained above, when the upstream perturbations are purely entropic the ratio K is normalized as in Ref. 1, which results in
and
for the two-dimensional and three-dimensional cases, where the corresponding normalized probability density functions are
Here K is unity when the integrand functions equal unity. Splitting (30) into longitudinal and transverse kinetic-energy contributions with the definition K = L + T for both two-dimensional and three-dimensional cases results at order $\epsilon \dot{q} |k|^{-1}$
in
The amplitudes Bue and Bve are provided in (A11a), and the components Aue and Ave in given in (A12).
Since the normalization of the turbulent kinetic energy is different for upstream vorticity $K_{\scriptscriptstyle {\rm A}r}$
and upstream density $K_{\scriptscriptstyle {\rm A}e}$
disturbances, to facilitate comparisons use is made of the exothermic deviation factor $\Delta _{\scriptscriptstyle {\rm K}}=K_{\scriptscriptstyle {\rm A}}/K_{\scriptscriptstyle {\rm B}}$
defined previously, which holds for both $K_{\scriptscriptstyle {\rm A}r}$
and $K_{\scriptscriptstyle {\rm A}e}$
, as does $\Delta _{\scriptscriptstyle {\rm L}}=L_{\scriptscriptstyle {\rm A}}/L_{\scriptscriptstyle {\rm B}}$
and $\Delta _{\scriptscriptstyle {\rm T}}=T_{\scriptscriptstyle {\rm A}}/T_{\scriptscriptstyle {\rm B}}$
. It can be seen in Fig. 4 that the upstream entropic disturbances (gray curves) produce positive deviations for $\Delta _{\scriptscriptstyle {\rm L}}$
, $\Delta _{\scriptscriptstyle {\rm T}}$
, and $\Delta _{\scriptscriptstyle {\rm K}}$
, while the upstream rotational fluctuations (black curves) generate negative deviations, except for the longitudinal $\Delta _{\scriptscriptstyle {\rm L}}$
at high Mach numbers. Since the effect of the heat-release rate $\dot{q}$
is factored out, the results are universal and valid for any q. They are, however, meaningful only if Mo ⩾ Mcj, and since Mcj increases with q according to (3), the range of relevance decreases with increasing q. To provide an indication of the ranges, the curves are shown dashed in the figure for Mo ≲ 2.4, which corresponds to Mo = Mcj for q = 1. For strong overdrive, the effect of heat-release rate on the turbulent kinetic energy reaches a constant value. As a representative indication of the magnitudes of the effects, it may be observed that a detonation with Mo = 3 will modify the kinetic energy by $\pm \ 10 \dot{q} |k|^{-1}\ \% $
, negative if rotational and positive if entropic. As expected from the results for monochromatic disturbances, the heat release decreases the turbulent kinetic energy of rotational disturbances for both thick and thin detonations, while with purely entropic initial fluctuations, the heat release in the detonation decreases it for thin detonations but increases it for thick detonations
FIG. 4.
Turbulent kinetic-energy deviations as functions of the shock Mach number Mo for two-dimensional (left panel) and three-dimensional (right panel) cases. The black and gray curves correspond to purely rotational (εe = 0) and purely entropic (εr = 0) disturbances upstream, respectively.
FIG. 4.
Turbulent kinetic-energy deviations as functions of the shock Mach number Mo for two-dimensional (left panel) and three-dimensional (right panel) cases. The black and gray curves correspond to purely rotational (εe = 0) and purely entropic (εr = 0) disturbances upstream, respectively.
Close modal
The effect of the inert shock front on the enstrophy levels is given in (A15) and agrees with previous results in shock-turbulence interactions.14,18 The two-dimensional and three-dimensional contributions of order $\epsilon \dot{q} |k|^{-1}$
are, respectively,
for rotational disturbances. For entropic disturbances, there is generation of vorticity, and following the same reasoning as in (30), we normalize the enstrophy downstream as in Refs. 1 and 15, so the contribution for both two and three dimensions become
Figure 5 shows the enstrophy deviation $\Delta _{\scriptscriptstyle {\rm W}}=W_{\scriptscriptstyle {\rm A}}/W_{\scriptscriptstyle {\rm B}}$
for γ = 1.2 as a function of the shock strength Mo, in the same fashion as in Fig. 4. It is seen that $\Delta _{\scriptscriptstyle {\rm W}}$
grows unbounded for very weak shocks when there are no vorticity perturbations upstream. This occurs because, although both $W_{\scriptscriptstyle {\rm A}e}$
and its base contribution $W_{\scriptscriptstyle {\rm B}e}$
(see Fig. 5 in Ref. 1) approach zero when Mo ≅ 1, the base contribution does so faster, but it is unimportant because weak-shock limits are not relevant to the Mach-number domains of common detonations. The general conclusions to be drawn from Fig. 5 are the same as those drawn from Fig. 4.
FIG. 5.
Enstrophy deviation as a function of the shock Mach number Mo for two-dimensional (left panel) and three-dimensional (right panel) cases. The black and gray curves correspond to purely rotational (εe = 0) and purely entropic (εr = 0) disturbances upstream, respectively.
FIG. 5.
Enstrophy deviation as a function of the shock Mach number Mo for two-dimensional (left panel) and three-dimensional (right panel) cases. The black and gray curves correspond to purely rotational (εe = 0) and purely entropic (εr = 0) disturbances upstream, respectively.
Close modal
The average of the downstream square-density perturbations is denoted by $\langle \rho _b^{\prime 2}\rangle$
, and for rotational upstream disturbances, following the same approach as for K, the factor G is defined as
for two-dimensional and three-dimensional cases, respectively, and the breakdown $G= G_{\scriptscriptstyle {\rm B}}+ \dot{q} |k|^{-1}G_{\scriptscriptstyle {\rm A}}$
is introduced. The base contribution from rotational disturbances contains the entropic and acoustic fluctuations shown in (A17), and the contribution at order $\epsilon \dot{q} |k|^{-1}$
is
for both two and three dimensions. For purely entropic fluctuations, the contribution at order εq|k|−1 is
Figure 6 shows the associated average-density deviation factor $\Delta _{\scriptscriptstyle {\rm G}}=G_{\scriptscriptstyle {\rm A}}/G_{\scriptscriptstyle {\rm B}}$
as a function of the shock Mach number for γ = 1.2. It is observed that upstream entropic disturbances generate relative small contributions for the downstream density field at order $\dot{q} |k|^{-1}$
, but the relative contribution of the upstream rotational disturbances is much greater because $G_{\scriptscriptstyle {\rm B}}$
is very small in that case. As in the previous Figs. 4 and 5, the factor $\Delta _{\scriptscriptstyle {\rm G}}$
is in qualitative agreement with the thin-detonation limit for rotational perturbations, but in the opposite direction for entropic perturbations.
FIG. 6.
Density deviation as a function of the shock Mach number Mo for two-dimensional (left panel) and three-dimensional (right panel) cases. The black and gray curves correspond to purely rotational (εe = 0) and purely entropic (εr = 0) disturbances upstream, respectively.
FIG. 6.
Density deviation as a function of the shock Mach number Mo for two-dimensional (left panel) and three-dimensional (right panel) cases. The black and gray curves correspond to purely rotational (εe = 0) and purely entropic (εr = 0) disturbances upstream, respectively.
Close modal
To complement previous analyses of interactions of detonation with turbulence that treated all turbulence scales as being large compared with the thickness of the detonation, in the present work, the turbulence scales were taken to be small compared with the thickness of the reaction zones that follow the lead inert shock in the detonation. Disturbances in the fresh mixture that involve only velocity fluctuations without any fluctuation of thermodynamic properties (rotational disturbances), and disturbances that involve fluctuations of the density of the fresh mixture without any velocity fluctuation (entropic disturbances) were considered separately, although results for cases in which both types of disturbances were present simultaneously are given in the Appendix. It was found that, unlike the previous results for thin detonations, where the jump conditions across the complete detonation determined how the initial fluctuations were modified, for these thick detonations the modifications are affected mainly by the jump conditions across the leading inert shock. The wavelengths, frequencies, and phases of the disturbances vary as they pass through the heat-release zones, but the amplitudes do not, to leading order in the perturbations. Moreover, the effect of the rate of heat release $\dot{q}$
makes itself felt only right behind the lead shock, so that reaction-rate distributions are irrelevant, so long as they do not produce appreciable enthalpy changes over distances comparable with the disturbance wavelengths. Since real detonations often have endothermic pyrolysis regions behind the lead shock, before the heat release sets in, the fluctuations actually may respond as if the detonation were to absorb energy rather than releasing heat.
A finding that was not unexpected is that the influences on monochromatic fluctuations are strongest when the interaction occurs near the critical angle of incidence that separates downstream disturbances behind the lead shock that are supersonic from those that are subsonic. This motivated exhibiting results in a variable stretched about the critical angle, to maximize the efficiency of presentation. The expansion parameter of the analysis $\dot{q} |k|^{-1}$
, essentially the product of the disturbance wavelength $|\bar{k}|^{-1}$
and the rate of heat release right behind the lead shock $\dot{\bar{q}}$
, divided by the unperturbed enthalpy flux at the Neumann state $\bar{\rho }_{\scriptscriptstyle {\rm N}}\bar{u}_{\scriptscriptstyle {\rm N}}\bar{a}_{\scriptscriptstyle {\rm N}}^2/(\gamma -1)$
, then factors out, and by normalizing quantities by their values at the Neumann state, parametric results can be presented with only the propagation Mach number Mo and the ratio of specific heats γ as parameters. A finding that was unexpected is that when the lead shock encounters the monochromatic fluctuations at an angle more head-on than the critical angle, and angle that results in generating propagating rather than evanescent acoustic perturbations, the influence of the heat-release rate on the disturbances is of second order in this small parameter. The leading-order results presented here, then, reflect only the influences of the interactions occurring at the larger angle θ, the encounters being more glancing than the critical angle. This is a peculiar result for thick detonations that does not have a counterpart for thin detonations.
For purely rotational disturbances (those without fluctuations of thermodynamic properties) in the fresh mixture, the heat release in the detonation reduces the relative fluctuation intensities downstream below the levels that would exist behind the inert lead shock, in both thin-detonation and thick-detonation limits. But for purely entropic disturbances in the fresh mixture (those without any velocity fluctuation), the effects of the heat release are the opposite in the two limits, the intensities being decreased for thin detonations but increased for thick detonations. This clearly suggests that when density fluctuations are dominant over velocity fluctuations, a situation thought to arise in turbulent flows at high Mach numbers,5 the influence of the passage of the detonation on the fluctuation field is affected qualitatively by the ratio of the detonation thickness to the range of sizes of the initial fluctuations. Smaller-scale fluctuations may be amplified, while longer-scale fluctuations are damped. The associated enhancement of smaller-scale disturbances may improve mixing rates and thereby promote combustion.
As indicated is Sec. IV, the base perturbations of order ε were obtained in previous work, namely, in Ref. 14 for rotational fluctuations and in Ref. 15 for entropic fluctuations. The formulas are different for short-wavelength (ζ ⩾ 1) and long-wavelength (ζ < 1) cases, and only the latter appears in Sec. IV. In addition, the base pressure perturbations are affected by a second long-wavelength contribution that is orthogonal to the contribution in Sec. IV and that therefore does not appear there. This second long-wavelength contribution will be identified by a prime, and the short-wavelength base perturbations will be identified by a double prime. Both of these influence later results, and therefore both are summarized here.
Addressing first the rotational contributions, it was found that14
The resulting amplitude of the asymptotic vorticity field generated by the upstream rotational mode is then
where Ω1 is the direct amplification contribution generated by compression effects, and Ω2 is the associated shock-curvature effect, here given by
The associated density perturbations are
For the entropic contributions it was found that15
The resulting the vorticity-disturbance coefficients are
where Ω3 is the baroclinic contribution (i.e., it accounts for the vorticity caused by the streamwise pressure jump across the shock and transverse density fluctuations) and is given by
The corresponding amplitudes of the steady density perturbations are
Besides velocity fluctuations associated with the vorticity fluctuations addressed here, there are also acoustic contributions to the velocity perturbations, which will be identified with the superscript a. For rotational disturbances, the associated amplitudes of the x and y components are
while the acoustic contributions are
The relationships in (A9a), written there for the long-wavelength cases, also apply for the short-wavelength cases. The corresponding contributions to the velocity field at order $\epsilon \dot{q} k_y^{-1}$
are
where the factor Aωr is defined in (21). For entropic disturbances the amplitudes are
for the rotational and acoustic contributions, respectively. Similarly, the vortical contributions to the velocity field at order at order $\epsilon \dot{q} k_y^{-1}$
are found to be
where the factors Aωe are also provided in (21).
The longitudinal $L_{\scriptscriptstyle {\rm B}}$
kinetic energy generated by the inert shock when traveling through an isotropic rotational field is
and the transverse $T_{\scriptscriptstyle {\rm B}}$
contributions to kinetic energy in two and three dimensions are, respectively,
and
When the upstream flow consists of isotropic solenoidal perturbations, the amplification of those perturbations in two and three dimensions is, respectively,
and
where the three-dimensional contribution can be split into $W_{\scriptscriptstyle {\rm B}r}=1/3+2/3W_{\perp }$
, W⊥ being the transverse contribution, $W_{\perp }= 1/4 \rho _o^{-2} + 3/4 W_z$
. The behavior of the transverse vorticity W⊥ can be found in Fig. 4 of Ref. 18 and Fig. 26 of Ref. 14. If there is no vorticity perturbations upstream, the enstrophy downstream is redefined as15
and detailed study of the enstrophy generated by shock-entropy interaction may be found in Refs. 15 and 20. The inert base contribution to the square-density perturbations behind the detonation, defined by (35), is1
where $\mathcal {P}$
stands for $\mathcal {P}_e$
or $\mathcal {P}_r$
depending on whether rotational or entropic perturbations are considered.
When there are both rotational and entropic fluctuations upstream, phase shifts introduce interferences that prevent results from being obtained by simply adding the separate contributions. Previous work21,22 has addressed the two cases $\epsilon _e = \epsilon _r M_o^{-1}$
and $\epsilon _e = - \epsilon _r M_o^{-1}$
with the same spectral density distribution for each type of fluctuations. Corresponding results from the present formulation are shown in Fig. 7, where it is seen that, for the moderate-to-strong shocks of interest, all of the deviation factors are negative, irrespective of the sign of the correlation. This may be expected from Figs. 4 to 6, in view of the domination of the rotational fluctuations at high Mo for the selected combination of εr and εe. The negative correlation does show a positive peak at lower Mo, but this is unlikely to lie in the range of realistic detonations. If εe had been chosen to be proportional to εrMo instead of $\epsilon _r M_o^{-1}$
, as it would be for the disturbances suggested by Morkovin,5 then the entropic contribution would be much more important at high Mach numbers. The choice made for Fig. 7 was to facilitate comparison with corresponding results in the literature.21,22
FIG. 7.
Turbulent kinetic-energy, enstrophy, and density deviations as functions of the shock Mach number Mo for two-dimensional (left panel) and three-dimensional (right panel) cases. The black and gray curves correspond to purely rotational (εe = 0) and purely entropic (εr = 0) disturbances upstream, respectively.
FIG. 7.
Turbulent kinetic-energy, enstrophy, and density deviations as functions of the shock Mach number Mo for two-dimensional (left panel) and three-dimensional (right panel) cases. The black and gray curves correspond to purely rotational (εe = 0) and purely entropic (εr = 0) disturbances upstream, respectively.
Close modal