Journal of Engineering Physics and Thermophysics, Vol. 86, No. 2, March, 2013
HEAT AND MASS TRANSFER IN COMBUSTION PROCESSES TALBOT METHOD OF INVESTIGATION OF THE CONCENTRATION DISTRIBUTION OF METHANE IN A TURBULENT AXISYMMETRIC JET P. P. Khramtsov, O. G. Penyazkov, I. N. Shatan, and I. A. Shikh
Fig. 1. Model of the experimental stand for investigation of turbulent flows.
Fig. 2. Scheme of the Talbot grating. was installed on a moving table for accurate positioning of the volume observed in the flow field of the jet. Commercial-purity methane (with a volume fraction of >99.9%) was used for the experiments. The initial velocity of outflow of methane at exit from the nozzle was u = 14 m/s (Re = uD/ν = 3900), which corresponds to a developed turbulent flow. An experimental setup for investigation of turbulent flows using the Talbot effect is shown in Fig. 1. The radiation of a gas laser with a wavelength of 514.5 nm was used as the probing radiation. The laser radiation was converted, with a collimator, to a homogeneous wide-aperture light beam of diameter 0.2 m with an approximately plane wave front. The Talbot matrix represented a reflecting grating with a system of equidistant holes of identical diameter (Fig. 2). It was arranged in the region of the laser beam immediately behind the collimator normally to the vector of propagation of the laser radiation and possessed the following parameters: optical density 4, period of the grating p = 1000 μm, and hole diameter 200 μm. For coherent radiation of the selected wavelength, the distance to the first self-reproduction plane was L = p2/λ = 1.944 m. The specimen was manufactured from a quartz plate coated with chromium. The repetition frequency of the period was ±2 μm. The investigated turbulent jet of methane was at a certain distance from the Talbot grating, which did not exceed the distance to the first self-reproduction plane. A light-sensitive matrix of a Nikon-D700 camera with a physical dimension of 24 × 36 mm and a number of active pixels of 12.1 Mp was placed in the corresponding zone of reproduction depending on the required sensitivity of the method. In an effort to match the apertures of the image of the Talbot matrix and the light-sensitive matrix of the camera, an aperture converter was installed immediately ahead of the camera. The probing radiation was recorded for the reference Talbot image and for the Talbot image of the methane jet. In the absence of the object under study, the field of radiation intensity in the self-reproduction plane represents a regular structure 260
Fig. 3. Differential Talbot image of the turbulent axisymmetric methane jet. of the set of light spots that corresponds to the structure of holes in the Talbot grating. The reference Talbot image recorded in this case determines the initial intensity field. Changes in the refractive index in the flow field of methane induce distortions of the wave front which lead to a displacement of the light spots in the Talbot image of the jet. Measurement of the value of displacement of each light spot enables us to determine the angles of deflection of light beams. Figure 3 gives the differential Talbot image of the methane jet under study, obtained for an exposure time of 10 s. Computational Methods for Determination of the Characteristics of the Axisymmetric Flow. Let us guide the symmetry axis of the turbulent jet vertically upward along the coordinate axis y and perpendicularly to the optional axis of the setup z (Fig. 4). For optical inhomogeneities of axial symmetry, the relationship between the deflection of the light beams and the distribution of the refractive index is expressed by the equations [14, 15]
2 εx = n0
x r2 − x2
2 dr , ε y = n0
r r2 − x2
where εx and εy are the projections of the angle of deflection of the light beams along the corresponding coordinate axes. It is assumed that the gradient angles of the refractive index is equal to zero for a radius r exceeding a certain finite radius R = R(y). Equations (1) represent a special form of the Abelian integral equation of the first kind. As a result of the replacement of the variables r = R 2 − t and x = R 2 − s , they take the form s
ϕ( s ) =
f (t )
dt . s−t
1 ∂n nε For the first equation of (1), the function f(t) corresponds to the function f(r) = , and the function φ(s) – φ(x) = 0 x . r ∂r x For Eq. (2), we know of the analytical solution
1 d f (t ) = π dt
ϕ( s )ds , t−s
⎞ ⎟. x 2 − r 2 ⎠⎟
which, for the first equation of (1), is of the form R n0 d ⎛⎜ ∂n = ∂r π dr ⎜ ∫r ⎝
ε x dx
Fig. 4. Coordinate system for mathematical modeling of the process of diagnostics of the jet. One usually does not use directly the analytical solution of the Abelian equation in integrating experimental data because of the presence of the derivative of the experimental function [16, 17]. Integration of the radial component for the gradient of the refractive index yields [14, 15] r
n( r ) − n( R ) =
∂n dr = ∂r
⎛n ⎜ 0 ⎜π ⎝
⎞ ⎟ dr = − n0 π x 2 − r 2 ⎟⎠ ε x dx
ε x dx x2 − r2
The given equation, in contrast to (4), contains no derivative of the experimental function. The singularity at the point x = r is easily integrated numerically, if we take into account that near the point x = r = 0, the function of experimental εx data must take on a value of zero. In diagnosing by the method of Talbot images, experimental data represent a discrete set of the parameters εx(x) ≡ ε kx dependent on the period of the Talbot grating. Numerical integration of Eq. (5) can be based on interpolation of the experimental data by the Lagrange polynomial L2(x) of the second degree on each segment [xk, xk+2] . In the present work, evaluation of the integral of (5) is based on interpolation of the εx data by a cubic spline. In this case the integral of (5) will be equal to the sum of the integrals on each segment of the corresponding polynomial of the third degree. The distribution (restored from formula (5)) of the averaged refractive index of the medium in the case of the isothermal jet may be used for calculation of the averaged concentration of methane in the flow field. The possibility of the concentration being determined further derives from the additivity property of the refractions of the components forming the mixture. Since the refractive indices of gaseous substances differ little from unity, the initial system of equations can be written as  n −1 =
1 ( N1a1 + N 2 a 2 ) , 2
N1 + N 2 = N ,
where the first equation follows from the Lorentz−Lorenz relation for a gaseous mixture, whereas the second represents the sum of particles in a unit volume. The concentration of the components Xi in this case is proportional to the change in the medium’s refractive index and can be expressed, for methane, by the equation X CH 4 =
n − n0 . nCH 4 − n0
Fig. 5. Three-dimensional distributions of the averaged deflection angles in the flow field of the axisymmetric methane jet: a) 〈εx〉; b) 〈εy〉.
Fig. 6. Distributions of the averaged deflection angles 〈εx〉 (a) and 〈εy〉 (b) in different cross sections of the axisymmetric methane jet: 1) y = 10D, 2) 15D, and 3) 20D. Results of Experimental Investigation of the Processes of Turbulent Mixing in the Axisymmetric Methane Jet. Figure 5 gives three-dimensional distributions (resulting from the processing of experimental Talbot images) of the averaged deflection angles 〈εx〉 and 〈εy〉 in the axisymmetric methane jet flowing out into an air space. Figure 6 shows the distributions of the deflection angles for cross sections of the jet. Different signs of deflections point to the fact that the displacement of disturbed individual elements of the Talbot image from the initial undisturbed ones occurred both right and left along the x axis and up and down along the y axis. The maximum values of the deflection angles along the x axis do not exceed 3 ⋅10 −4 rad ( 60′′ ) in magnitude; along the y axis, they are no higher than 2.5 ⋅10 −5 rad ( 5′′ ). On the symmetry axis, the values of the two projections of the deflection angles are equal to zero. Figure 7 gives the distribution of the averaged refractive index in the axisymmetric methane jet, that has been obtained from the distribution of the deflection angles. Results for the averaged methane concentration with a maximum value on the jet axis, which has been calculated from these data, are presented in Fig. 8. A reduction in the methane concentration downstream of the jet is due to the turbulent mixing of methane and air. Figure 9 shows the reciprocal concentration of methane along the jet axis as a function of the dimensionless distance to the nozzle exit section y/D. The straight line drawn through the experimental data is a result of approximation by a linear function. To construct the approximation we have used the data whose abscissa satisfied the condition y/D ≥ 7. The functional dependence of the concentration along the jet axis on the distance to the nozzle exit section is expressed by the formula [2, 5]
( X CH 4 )0 ( X CH 4 )C
C1 ( y − y01 ) , D
Fig. 7. Three-dimensional radial distribution of the refractive index throughout the flow field of the axisymmetric jet (a) and for the jet cross sections (b): 1) y = 10D, 2) 15D, and 3) 20D.
Fig. 8. Distributions of the averaged concentration of methane in the axisymmetric jet throughout the flow field (a) and for the jet cross sections (b): 1) y = 10D, 2) 15D, and 3) 20D.
where ( X CH 4 )0 = 1 is the concentration of methane at exit from the nozzle, C1 is the constant determining the slope of the straight line, and y01 is the position of the virtual jet source. For the data given in the present work, the indicated parameters are C1 = 0.16 and y01 = 1.5D. The opening angle of the jet is characterized with the aid of r1/2, i.e., of the radial distance at which the concentration is halved compared to the value on the jet axis. At a sufficient distance from the nozzle’s edge, the value of r1/2 can be approximated by a function proportional to the distance from the virtual source y02. This dependence is expressed by the formula  r1 2 = C2 ( y − y02 ) ,
where C2 is the parameter characterizing the opening angle of the jet, and y02 is the position of the virtual source corresponding to the opening angle of the jet. Figure 10 gives r1/2 values which have been calculated from the experimental data obtained in the present work and as a result of their approximation by the function (10) in the domain y/D ≥ 7. A comparison of the results of measurement of the
Fig. 9. Inverse concentration of methane along the jet axis vs. distance to the nozzle exit section. Fig. 10. Dependence of r1/2 on the distance to the nozzle exit section.
Fig. 11. Radial distributions of the averaged concentration of methane: 1) data of the present work for the cross section y = 22.7D; 2) data of . methane concentration on the basis of the method of averaged Talbot images to the results obtained in other investigations was based on the approximations of (9) and (10). We selected  and  in which methane-air flows have also been investigated. In , a study was made of the outflow of a methane jet into a slowly moving concurrent air flow, whereas in , the methane jet flowed out into a quiescent air space. The value of the Reynolds number for the jet investigated in the present work does not differ, in practice, from the value given in  but noticeably differs from its value in . The method employed in  for measurements of the concentration of the components of a binary mixture is based on the Rayleigh scattering of laser radiation. In , the methane jet was investigated using the method based on Raman scattering. Table 1 gives the value of the parameters in Eqs. (9) and (10) for the data measured in the present work and published in  and . In comparing the values of C1 given in the table, we see that the disagreement between these data amounts to as much as 10%. It should be taken into account that this coefficient is sensitive to the selection of the approximation domain y/D. Therefore, it is likely that on selection of the corresponding region, the disagreement between the data will decrease. The values of the parameter C2 for the approximation of (9) are in good agreement with the results of other authors. From Table 1, it is clear that there is a strong similarity of the results, despite the fact that the Reynolds number in  is four times larger, whereas in , the jet flows out into a concurrent flow. The differences in virtual sources point to dissimilar initial conditions of jet outflow in all the three cases. Figure 11 shows the radial concentration distribution of methane for the cross section 22.7D, obtained in the present work, and the data presented in Fig. 10 in . The cross sections are selected so that the distance to the virtual source y02 is the same. As is clear from the figure, the differences in the data are minor despite the absence of a concurrent air flow. 265
TABLE 1. Coefficients of Eqs. (9) and (10) for Experimental Data Source Experimental results
Data from 
Data from 
C2 0.107 0.108 0.105 0.097
y02 1.95D 0 −4.6D 0
On the basis of the good agreement between the characteristics of the concentration field that have been measured in the present work and have been published in  and  and of the absence of experimental proofs pointing to the significant difference in the values of the parameters measured, we can infer that the presented data provide a reliable description of the concentration field in the axisymmetric jet. Conclusions. In this work, we have presented the results of application of the method of Talbot images to determination of the concentration field in the binary mixture of an isothermal turbulent flow. It has been shown that the procedure ensures measurements of the field of deflection angles of probing radiation with a high spatial resolution. From the time-averaged values of the deflection angles, we have restored the distribution of the refractive index in an axisymmetric turbulent jet. Taking into account that the change in the refractive index in the binary isothermal jet is directly proportional to the change in the quantitative composition of the mixture, we have calculated the averaged concentration of methane along the jet flow. A comparison of the presented data to the results of concentration measurements, obtained by other methods, has been made. The possibility of measuring the spatial distribution of the averaged concentration of the mixture components makes this method highly promising for investigation of the processes of mixing. The use of high-speed image recording together with the presented procedure will enable one to investigate turbulent processes with a high temporal resolution. In the future, it is planned to use the method of Talbot images for measurement of the dispersion dependence of the refractive index in the flow field of a reactive methane jet.
NOTATION a, molecular polarizability, m3; D, diameter of the nozzle outlet cross section, m; F, focal length of the objective, mm; L, self-reproduction distance, m; n, refractive index of the medium; N, number of particles in a unit volume, m −3 ; p, period of the grating, m; r, radius in a cylindrical coordinate system, m; R, radius of an optical inhomogeneity, m; u, velocity of flow of the methane jet out of the nozzle, m/s; X, relative concentration of the component of the mixture; x, y, and z, Cartesian coordinate system, m; ε, angle of deflection of light, angular seconds; λ, probing-radiation wavelength, nm; ν, kinematic viscosity, m2/s. Subscripts: 0, undisturbed medium; k, reproduction-zone number.
REFERENCES 1. K. Kohse-Höinghaus, Laser techniques for the quantitative detection of reactive intermediates in combustion systems, Prog. Energy Combust. Sci., 20, No. 3, 203–279 (1994). 2. W. M. Pitts and T. Kashiwagi, The application of laser-induced Rayleigh light scattering to the study of turbulent mixing, J. Fluid Mech., 141, 391–429 (1984). 3. R. W. Schefer and R. W. Dibble, Rayleigh scattering measurements of mixture fraction in a turbulent nonreacting propane jet, AIAA 24th Aerospace Sci. Meeting, Paper No. AIAA-86-0278 (1986). 4. R. W. Dibble, W. Kollmann, and R. W. Schefer, Conserved scalar fluxes measured in a turbulent nonpremixed flame by combined laser Doppler velocimetry and laser Raman scattering. Combust. Flame, 55, 307–321 (1984). 5. A. D. Birch, D. R. Brown, M. G. Dodson, and J. R. Thomas, The turbulent concentration ﬁeld of a methane jet, J. Fluid Mech., 88, 431–449 (1978). 6. H. F. Talbot, Facts Relating to Optical Science. Phil. Mag., 9, No. IV, 401–407 (1836). 7. A. S. Koryakovskii and V. M. Marchenko, Talbot-effect-based interferometry of optical inhomogeneities of active laser media, Kvantovaya Élektron., 7, No. 5, 1048–1057 (1980). 8. A. S. Koryakovskii and V. M. Marchenko, The Talbot-effect-based wave front transducers, Zh. Tekh. Fiz., 51, No. 7, 1432–1438 (1981). 266
9. M. V. Doroshko, O. G. Penyazkov, G. Rankin, K. L. Sevruk, P. P. Khramtsov, and I. A. Shikh, Talbot interferometry in measurements of the parameters of an axisymmetric turbulent jet, Inzh.-Fiz. Zh., 79, No. 5, 94–99 (2006). 10. M. V. Doroshko, O. G. Penyazkov, G. Rankin, K. L. Sevruk, P. P. Khramtsov, and I. A. Shikh, Investigation of twodimensional turbulent jet flow around a horizontal cylinder by the Talbot interferometry method, Coll. of Sci. Papers "Heat and Mass Transfer–2005," Minsk (2005), pp. 210−214. 11. M. V. Doroshko, P. P. Khramtsov, O. G. Penyazkov, and I. A. Shikh, Measurements of admixture concentration fluctuations in a turbulent shear flow using an averaged Talbot image. Exp. Fluids, 44, 461–468 (2008). 12. M. V. Doroshko, O. G. Penyazkov, P. P. Khramtsov, and I. A. Shikh, Measurement of concentration pulsations in a shear turbulent flow by the method of average Talbot images, Inzh.-Fiz. Zh., 81, No. 1, 49–54 (2008). 13. M. V. Doroshko, O. G. Penyazkov, P. P. Khramtsov, and I. A. Shikh, Statistical characteristics of density pulsations in the boundary layer behind the shock wave front, Coll. of Sci. Papers "Heat and Mass Transfer–2007," Minsk (2007), pp. 213−219. 14. L. A. Vasil’ev, Shadow Methods [in Russian], Nauka, Moscow (1968). 15. M. M. Skotnikov, Shadow Quantitative Methods in Gas Dynamics [in Russian], Nauka, Moscow (1976). 16. G. N. Minerbo and M. E. Levy, Inversion of Abel’s integral equation by means of orthogonal polynomials, SIAM J. Numer. Anal., 598–616 (1969). 17. V. K. Singh, R. K. Pandey, and O. P. Singh, New stable numerical solution of singular integral equations of Abel type by using normalized Bernstein polynomials, Appl. Mat. Sci., 3, No. 5, 241–255 (2009).