Quantitative Schlieren Image-Noise Reduction Using Inverse Process and Multi-Path Integration

Ahmad Zaid Nazari^{1},
Yojiro Ishino^{1}^{*},
Fumiya Ito^{1},
Harumi Kondo^{1},
Ryoya Yamada^{1},
Takanori Motohiro^{1},
Yu Saiki^{1},
Yoshiaki Miyazato^{2},
Shinichiro Nakao^{2}

Show more

1. Introduction

Schlieren imaging technique is a common tool in science and technology to visualize density gradients and investigate phenomena with non-uniform density flows in transparent media. Recently, this technique is reviewed in [1] and is shown that schlieren photography has been developed from a qualitative visualization method into a quantitative measurement method. In previous works [2] - [7], by employing and developing a non-scanning three-dimensional computerized tomography (3D-CT) technique using a delicate multi (20)-directional quantitative schlieren optical system with a flashlight source, measurement of the instantaneous density distributions of several types of flames in laminar and turbulent flows [2] [3] [4] [5] [6] and supersonic micro-jets [7] have been successfully obtained.

In the previous work [2], the target high-speed turbulent burner flames are investigated and reported. Figure 1 depicts the target flame, details dimensions of burner nozzle, coordinate system and 3D bird’s-eye view of the CT-reconstructed density data. For more information related to the flow conditions, measured parameters of target flames and burner details refer to the previous work [2].

In the present report, two new techniques are introduced based on a novel concept of complex brightness gradient in the quantitative schlieren images, “inverse process” and “multi-path integration” for image-noise reduction. The new techniques partially are presented in international conferences [8] [9], in detail and entirely will be discussed here. Furthermore, it will be shown an important ability of complex schlieren brightness gradient, which is independence on the path and starting point in the integration process.

The schlieren images brightness change gradually (brightness gradient), while

(a) (b) (c)

Figure 1. High-speed turbulent flame, burner, coordinate system and bird’s-eye view of 3D model. (a) Target flame & schlieren image; (b) Burner tip details; (c) 3D model view.

image noise is a random variation of brightness and is generally considered undesirable. Noise can deduct the quality of the image and quantitative analysis [10]. Noise in the schlieren images can arise from density gradients anywhere between the source slit and the knife-edge in the wind-tunnel research, and that is why a multiple-source schlieren system is introduced for noise reduction in [11]. Quantitative schlieren imaging requires low noise level images accordingly high data accuracy. Noise in the schlieren imaging is unavoidable. There are very few noise reduction techniques studies in the works of literature. Therefore, the lack of proper noise reduction techniques is clear. Some image noise generating sources are as follows: dusty photography surroundings, dusty and dirty lenses, defective lens surface, digital cameras, and lenses. In the present study, a new approach called “multi-path integration” novel technique against former technique “single-path integration” is introduced for noise reduction in projections (density thickness) images of computerized tomography (CT).The new multi-path integration technique is required schlieren brightness gradients in both horizontal (*x*-directional) and vertical (*z*-directional) directions, or in other words, in two orthogonal directions (e.g. horizontal and vertical or two perpendicular diagonal directions). The 20-directional quantitative schlieren optical system gives only images of schlieren brightness in the horizontal gradient (*x*-directional) and another 20-directional optical system seems necessary to obtain vertical (*z*-directional) schlieren brightness gradient, simultaneously. The present investigation represents the new “inverse process” novel technique for obtaining vertical and diagonal schlieren brightness gradient from horizontal experimental data (with vertical knife-edge orientation) without the necessity of a new optical system. The main advantages of the new approach employed are its relative simplicity, the ease of application, without extra cost, and the satisfactory accuracy of results and can be used for obtaining any optional direction of schlieren brightness gradients. Recently [8], we presented the “inverse process” technique but in a different way than this time. In the previous work, the “inverse process” was performed after the CT-reconstruction, however, now is done before the CT-reconstruction procedure (Figure 2). Indeed, in the previous work phantom (virtual) data of schlieren images are obtained by the “inverse process” and used to evaluate noise reduction technique while in the new “inverse process” technique actual experimental data is used for reproducing other directions of schlieren brightness gradient.

By using both white-light and monochromatic sources, color schlieren images have been produced from combining horizontal and vertical gradient brightness of common schlieren images by [12], but for this purpose, two independent knife-edges along two different beam paths are employed. About 70 years ago, a more complex arrangement is presented for obtaining two spatially separated images sensitive to orthogonal path deviations [13]. A good example of the same numerical supersonic jet-exit flow field is rendered as shadowgraph, bright-field schlieren (circular cutoff), vertical and horizontal knife-edge schlieren in [14].

Figure 2. Flow diagram of the “inverse process” new technique for the reproduction of schlieren brightness gradient in horizontal, vertical, diagonal directions and conversion to the new projections (density thickness) of CT by new “multi-path integration” technique.

In the current investigation, the horizontal schlieren brightness gradient data of high-speed turbulent flame [2] with average flow velocity 10 m/s is used for reproducing the 20-directional images of schlieren brightness in horizontal, vertical and diagonal gradients.

First, as shown in Figure 2 (top), the initial density thickness (*Dt*) images are obtained from images of schlieren brightness in the horizontal gradient (*x*-directional). Next, using the inverse process from density thickness (*Dt*) images, the horizontal (*x*-directional), vertical (*z*-directional) and two diagonals (45° and 135˚) schlieren brightness gradient data of target flame have been successfully obtained. Finally, by applying the multi-path integration technique, the new set of projections images of density thickness are successfully obtained (Figure 2 (bottom)). Then, by using the new set of density thickness images “with decreased noise” as projections for CT-reconstruction the 3D density distributions of target flame have been successfully reconstructed.

2. Experimental Apparatus Setup and Methods

2.1. Multi-Directional Quantitative Schlieren Optical System

Figures 3(a)-(c) illustrates the concept of a 20-directional schlieren camera. In the camera system, the target flame/non-uniform density field is observed from 180˚ direction using numerous schlieren optical systems simultaneously from *θ* = –85.5˚ to +85.5˚ at an interval of 9˚. Here, angle *θ* is defined as the horizontal angle from the *x*-axis. For pre-investigation and for time series observation (high-speed Schlieren movies) of the target flow, a high-speed camera (HSC) is used simultaneously with 20-direction Schlieren photographing apparatus (between cameras No. 13 and No. 14, Figure 3(b)).

(a) (b) (c) (d)

Figure 3. Multi-direction schlieren optical system and coordinate system [2] - [9]. (a) Schlieren camera system; (b) Coordinate system; (c) Schematic of 20-direction schlieren camera; (d) Single unit of quantitative schlieren optical system.

The diagram of a single instance of the multi-directional quantitative schlieren camera system is depicted in Figure 3(d). This system is composed of two convex achromatic lenses of 50 mm diameter and 300 mm focal length, a flashlight source unit, a vertical knife-edge (for obtaining horizontal schlieren brightness gradient images), and digital camera.

The neutral-density (ND) filter (Fuji 1.5, exposure adjustment multiple 3.16) and stepped neutral-density filter is used for calibrating the cameras. The light unit is a xenon flashtube that emits full-spectrum white light with a uniform luminance rectangular area of 1 mm × 2 mm and a 35 μs duration.

2.2. Image Processing

Figure 4 shows a conceptual framework for the formation of a projection image (density thickness) in the proposed schlieren 3D-CT system. Detailed descriptions of the various sections of Figures 4(a)-(k), along with the supporting formula, are given in the previous reports [2] - [9] and many other references cited therein. Here, for brevity, the main points of Figure 4 are explained as follows. The first goal is obtaining the density thickness of deviation density (*Dt**(*X*(*θ*))) which shown in Figure 4(d) with a unit of (kg/m^{3})(m). The density thickness of deviation density *Dt**(*X*(*θ*)) is obtained automatically from schlieren observation using spatial integration of deviation density *Δρ**(*X*(*θ*),*Y*) along the line of sight. In practice for obtaining density thickness of deviation density (Figure 4(i)), the image processing activity starts from Figure 4(e) by obtaining two sets of images, “with target” and “without target” (without any disturbance in the test section) with a horizontal brightness gradient in the *x*-direction. Two sets of images present by schlieren observation as *B*(*X*) and *B** _{nf}*(

Figure 4. Processes of formation of schlieren brightness and conversion to projections (density thickness (*Dt*)) of CT.

image in no-flame condition). To obtain the density thickness of deviation density *Dt*(*X*(*θ*)) from *B*(*X*) and *B** _{nf}*(

$\Delta B=B\left(X\right)-{B}_{nf}\left(X\right)$ (1)

is scaled to d(*Dt*)/d*X* by

$\text{d}\left(Dt\right)/\text{d}X=-\left(1/K\right)\left(\Delta s/f\right)\left[\Delta B\left(X\right)/{B}_{nf}\left(X\right)\right]$ (2)

where *K* is Gladstone-Dale constant for air (*K* = 2.26 × 10^{−4} m^{3}/kg), *Δs* (= 1 mm × 2 mm (Hor*.** *× Ver*.*)) is the transparent width of the light source image on schlieren stop location and *f** *(= 300 mm) is the focal length of the convergent lens. Deviation density thickness *Dt*(*X*(*θ*)) is, therefore, reproduced by transverse integration of d(*Dt*)/d*X* from schlieren images, as shown in Figure 4(h). Finally, the “density thickness” *Dt**'*(*X*(*θ*)) is obtained by adding the thickness of (
${\rho}_{ref}\u2013{\rho}_{a}^{\ast}$ ) to the deviation density thickness *Dt*(*X*(*θ*)) on peripheral of the observed range of radius R (Figure 4(i)), as expressed by

$Dt\prime \left(X\left(\theta \right)\right)=Dt\left(X\left(\theta \right)\right)+2{\left({R}^{2}-{X}^{2}\left(\theta \right)\right)}^{1/2}\left({\rho}_{ref}-{\rho}_{a}^{\ast}\right)$ (3)

In the present study, ${\rho}_{ref}={\rho}_{a}^{\ast}$, accordingly, $Dt\prime \left(X\left(\theta \right)\right)=Dt\left(X\left(\theta \right)\right)$.

2.3. CT-Reconstruction

Density thickness images are used for CT-reconstruction by maximum likelihood-expectation maximization (ML-EM) [15] an appropriate CT algorithm to obtain the 3D reconstruction of density distribution, the ML-EM method [2] - [9] is employed for CT-reconstruction. The CT procedure is carried out in each horizontal plane of *z*-axis for the reconstruction of deviation density distribution *Δρ*(*x*,*y*) from a linear data set of density thickness* Dt'*(*X*(*θ*)) (Figure 4(j)). Finally, 2D density distribution (Figure 4(k)) is obtained as follow.

$\rho \left(x,y\right)={\rho}_{a}^{\ast}-\Delta \rho \left(x,y\right)$ (4)

The reconstruction was performed cross-section by cross-section and then the cross-sections were stacked to form a three-dimensional density distribution. Therefore, a 2D distribution *ρ*(*x*,*y*) is accumulated in layers to form the 3D-CT distribution *ρ*(*x*,*y*,*z*). In the present study, the “density thickness” projections images of 400 (*H*) × 500 (*V*) pixel (32 mm × 40 mm) are used for CT-reconstruction to produce 3D data 400 (*x*) × 400 (*y*) × 500 (*z*) pixel (32 mm × 32 mm × 40 mm). The voxel size is 0.08 mm in each direction.

3. New Techniques and Treatments

3.1. Inverse Process

As mentioned before, in the present study, using “inverse process” new approach from density thickness (*Dt*) images (Figure 2), the horizontal (*x*-directional), vertical (*z*-directional) and two diagonals (45˚ and 135˚) schlieren brightness gradient data of target flame have been successfully obtained. The relation between deviation brightness ∆*B* and density thickness *Dt* is expressed by Equation (2). Equation (2) can be rewritten in the following simple forms.

$\text{d}\left(Dt\right)/\text{d}X\propto \Delta B\left(X\right)$ (5)

and

$\text{d}\left(Dt\right)/\text{d}X=\alpha \times \Delta B\left(X\right)$ (6)

where *α* is a coefficient with parameters that explained in Equation (2). The deviation density thickness can be expressed as,

$\text{d}\left(Dt\left(x,z\right)\right)/\text{d}X=Dt\left(x,z\right)-Dt\left(x-1,z\right)$ (7)

Finally, in the inverse process, by combining the Equation (6) and (7), the deviation brightness of the quantitative schlieren image in each direction will be given by:

$\Delta B{\left(x,z\right)}_{Horizontal,+x}=\left(1/\alpha \right)\times \left[Dt\left(x,z\right)-Dt\left(x-1,z\right)\right]$ (8)

$\Delta B{\left(x,z\right)}_{Vertical,+z}=\left(1/\alpha \right)\times \left[Dt\left(x,z\right)-Dt\left(x,z-1\right)\right]$ (9)

$\Delta B{\left(x,z\right)}_{Diagonal,{45}^{\circ}}=(1/\alpha )\times \left[Dt\left(x,z\right)-Dt\left(x-1,z-1\right)\right]$ (10)

$\Delta B{\left(x,z\right)}_{Diagonal,{135}^{\circ}}=(1/\alpha )\times \left[Dt\left(x,z\right)-Dt\left(x+1,z-1\right)\right]$ (11)

The direction of the schlieren brightness gradient can be changed in the opposite direction (Figures 5(c2)-(f2)) only by applying one minus sign in the corresponding equation. Therefore, by using density thickness *Dt* in the inverse process, the schlieren brightness in horizontal, vertical and diagonal directions can be calculated (Equation (8)-(11)) without the necessity of a new optical system.

One of the benefits of using this approach is that now we can apply the “multi-path integration” noise reduction technique. The sample results for only camera number 10 are depicted for all mentioned directions in Figure 5. The corresponding result for numerical simulation (refer to section 3.4) of the flame is depicted in the top right side and the corresponding knife-edge orientation is shown underneath of each schlieren image in Figure 5.

A key element in the schlieren setup is the knife-edge. In traditional schlieren imaging, gradients in the refractive index normal to the imaging plane are viewed by focusing the deflected light rays onto a knife-edge. The rays that are deflected above the knife-edge appear as bright regions in the image while rays that are deflected onto the knife-edge appear as dark regions in the image. The knife-edge serves as a cut-off filter for light intensity (brightness). The knife-edge orientation will reveal directional density gradients (brightness gradients). For example, vertical knife-edge orientation present horizontal

Figure 5. Inverse process for reproduction of schlieren brightness in horizontal, vertical and diagonal gradient.

schlieren brightness gradient data and horizontal knife-edge orientation present vertical schlieren brightness gradient data. In the 20-directional quantitative schlieren optical system, a vertical knife-edge (Figure 5(c1)) is used which for adjusting moves from left to right in the positive *x*-direction.

3.2. Complex Schlieren Brightness Gradient

Figure 6 shows an example of contour lines display of density deviation thickness (density thickness). At this figure, the brightness value of the Schlieren image is given by the gradient (slope) of the curved surface. For example, when a knife-edge parallel to the *z*-axis is used as a schlieren stop (vertical knife-edge) and the light beam is cut in the +*x* direction, the gradient (slope) of the red curve can be obtained by measuring the brightness of the schlieren image. On the other hand, if the knife-edge is inclined 90°, the slope of the blue curve can be obtained.

Therefore, by simultaneously considering the schlieren brightness gradient information in the horizontal direction (red, real part) and vertical direction (blue, imaginary part), the complex gradient of the density thickness distribution at any position can be obtained. The spatial integration to obtain the density thickness distributions can now be performed on any arbitrary path, and noise reduction now becomes possible. In this study, a schlieren image

Figure 6. Complex schlieren image explanation.

with information on the gradient of density thickness in two orthogonal directions is called a “complex schlieren image”.

3.3. Cauchy Integration Theorem

The Cauchy integration theorem in complex analysis is an important statement about line integrals for regular functions in the complex plane. Cauchy’s integration theorem states that if *f*(*z*) be a regular function in the simply connected domain *D*, its contour integral on a path does not depend on the integration path and depends only on the beginning and ending points of the path. Therefore, if let *C* be a piecewise continuously differentiable path in *D* with start point *a* and end point *b*, and *F* is a complex antiderivative of *f*, then,

${\int}_{C}f\left(z\right)}={\displaystyle {\int}_{a}^{b}f\left(z\right)}=F\left(b\right)-F\left(a\right)$ (12)

If *a* = *b*, then *C* is a closed curve (loop) within *D*, accordingly,

${\oint}_{C}f\left(z\right)}=0$ (13)

We applied the above theorem to develop the noise reduction technique. As shown in Figure 7 and Figure 10, for obtaining density thickness images “with decreased noise” the average of transverse integration from multi-path is employed.

It is remarkable, as a new finding, a new important concept in the complex schlieren image could be expressed as follow; The complex schlieren image for obtaining density thickness at an arbitrary position (ending point of the integration path) has one more potential power not depends on the starting point of the integration path. Therefore, the calculation of density thickness can be performed just by starting integration from any point on any path (Figure 10).

3.4. Numerical Simulation

First, as a preliminary investigation and for a deep understanding of complex schlieren brightness, a numerical simulation for schlieren brightness of a flame is performed. The results of the numerical simulation are depicted in Figure 7. Schlieren images in horizontal, vertical and diagonal directions of brightness gradients are shown in Figures 7(a1)-(h1). One hundred noises with various brightness in different positions are added to the schlieren images randomly.

Figure 7. (a1)-(h1) Numerical simulation of schlieren brightness gradient in the horizontal, vertical and diagonal directions of flame, (a2)-(h2) Noisy density thickness (*Dt*) brightness using "single-path integration" underneath of the corresponding each schlieren image and (i) Density thickness brightness using "multi-path integration" with reduced noise effects.

The noisy density thickness images obtained by single-path integration (former technique) are depicted underneath of corresponding each schlieren images (Figures 7(a2)-(h2)). Figure 7(i) illustrates the density thickness image by employing the “multi-path integration” new technique. Figure 7(i2) depicts the result of the density thickness image in the absence of noise in the schlieren image.

Figure 8 depicts a schlieren brightness gradient (Figure 8 (left)) at the yellow dashed line position for the numerical simulation result (Figure 8 (right)).The brightness of schlieren image changes gradually (brightness gradient), while image noise is the variation of brightness suddenly (Figure 8) and is generally considered undesirable. The existence of some noises in the schlieren image affects the density thickness image (Figure 9). Noise effect like an infectious disease spreads from a noisy pixel to the next pixel in the direction of single-path integration former technique (Figure 9(c)). One spot noise in the schlieren image appears in a line shape in the density thickness image, which starts from noise position and continues as a line in the direction of single-path integration (Figure 9(c)). As seen in Figure 7(i1) and Figure 9(d), using the new technique the noise effects are reduced to an acceptable level and slightly present around

Figure 8. Schlieren brightness gradient at the yellow dashed line position with noise effects.

Figure 9. Noise effects on the density thickness (*Dt*) brightness gradient (simulation result).

initial noise spots (this may be a disadvantage of the proposed technique).

Noise in schlieren imaging is unavoidable and can detract from the quality of the image and subsequent visual and quantitative analysis. Therefore, the need for proper noise reduction techniques is clear.

3.5. Multi-Path Integration Technique

As mentioned before, a new approach called “multi-path integration” novel technique against the former technique “single-path integration” is introduced for noise reduction in the projections (density thickness) images of CT. The new multi-path integration technique is required schlieren brightness gradients in both horizontal (*x*-directional) and vertical (*z*-directional) directions, or in other words, in two orthogonal directions (e.g. horizontal and vertical or two perpendicular diagonal directions) (Figure 2).

In the usual procedure, density thickness image reproduces by transverse integration of the schlieren image brightness in the single-path along schlieren brightness gradient (path No. 4 in Figure 10 (left)).

In the present study, for obtaining density thickness images “with decreased noise” the average of transverse integration from multi-path is employed. These paths are in four main horizontal and vertical directions (totally 28 paths, 7 paths in positive and 7 paths in negative of *x* and *z* directions) (Figure 10), as well as, the paths can be used from two perpendicular diagonal directions. The new procedure “multi-path integration” is one comprehensive technique that including the usual former “single-path integration” technique as well. In the former technique “single-path integration” (path No. 4, Figure 10), the density thickness for pixel (*x*, *z*) can be expressed from Equation (7) as follow.

$Dt\left(x,z\right)=Dt\left(x-1,z\right)+\text{d}\left(Dt\left(x,z\right)\right)/\text{d}x$ (14)

On the right side of Equation (14), the first term is the density thickness data from prior adjacent pixel (the density thickness for the first-pixel column (1, *z*) is assumed zero) and the second term comes from the schlieren image brightness (Equation (2)). This second term for *z* (vertical) direction can be obtained by applying the inverse process technique. Calculation of two sample paths (No. 6 and No. 16) in the new technique “multi-path integration” is shown in the following equations.

$Dt\left(x,z\right)=Dt\left(x-1,z+2\right)+\frac{\text{d}\left(Dt\left(x,z+2\right)\right)}{\text{d}x}-\frac{\text{d}\left(Dt\left(x,z+1\right)\right)}{\text{d}z}-\frac{\text{d}\left(Dt\left(x,z\right)\right)}{\text{d}z}$ (15)

$Dt\left(x,z\right)=Dt\left(x-2,z-1\right)+\frac{\text{d}\left(Dt\left(x-2,z\right)\right)}{\text{d}z}+\frac{\text{d}\left(Dt\left(x-1,z\right)\right)}{\text{d}x}+\frac{\text{d}\left(Dt\left(x,z\right)\right)}{\text{d}x}$ (16)

First, the density thickness for pixel (*x*,*y*) is obtained from the multi-path (28 paths). Then, the average of all these paths is considered as the density thickness for the target pixel.

In the schlieren imaging system, brightness gradually increases (or decreases) from ambient condition (with constant brightness) to boundary condition and vice versa. This fact is effective whether schlieren photography is containing a

Figure 10. Integration paths (28 paths) for calculating density thickness of pixel (*x*, *z*).

complete brightness gradient in the integrating direction or not. Figure 11 represents complete (Figure 11(a)) and incomplete (Figure 11(b)) observation of the schlieren brightness gradient for the laminar flame kernel. Hence, if the schlieren photography in one and more directions is not bounded to an ambient condition which means the schlieren photography is not observed complete target flow including the boundary condition, those directions could be omitted (Figure 11(b) (Top direction)), because this treatment will give better results for noise reduction using multi-path integration. In the present study, paths from negative *z*-direction (7 paths, No. 22 to 28, Figure 10) are not considered in the calculation due to that top area of target flame is not captured by the schlieren photography system, because of the fact that lenses diameters are smaller than target flame height.

4. Results and Discussions

4.1. Quantitative Schlieren Brightness Images

Figure 12 shows quantitative schlieren images of the horizontal brightness deviation from the original experimental data of high-speed turbulent premixed flame. The numbers inserted in each picture express the shooting angles *θ*. The images of Figure 12 are processed by the above-mentioned “inverse process” to obtain quantitative schlieren images of vertical and diagonal brightness deviation. Due to page limitation, as a sample just for one direction for camera No.10, the results are depicted for new horizontal, vertical and diagonal (two directions) deviation brightness of quantitative schlieren images in Figure 13.

4.2. Projections Images of Density Thickness (*Dt*)

For better evaluation of the “multi-path integration” noise reduction technique, in the new set of 20-directional schlieren images obtained from the inverse process, one thousand noises are distributed with various brightnesses and in the

(a) (b)

Figure 11. Sample schlieren images of the laminar flame kernel [4]. (a) Complete observation; (b) Incomplete observation.

Figure 12. Original quantitative schlieren images of target flame (*U* = 10 m/s) [2].

Figure 13. Schlieren brightness gradient for (a1) horizontal, (b1) vertical and (c1)-(d1) diagonals directions, (f) conversion to the new projection (density thickness (*Dt*)) of CT by new “multi-path integration” technique, against (a2)-(d2) noisy density thickness (*Dt*) images by former “single-path integration” technique. (The schlieren images (a1-d1) and density thickness images (a2-d2) are depicted as a sample just for one direction for camera No. 10).

different positions randomly (Figure 13(e)). Then, by using the former technique “single-path integration”, in horizontal (Figure 13(a1)), vertical (Figure 13(b1)) and diagonal directions (Figure 13(c1) and Figure 13(d1)), the projections images of density thickness with the existence of noise effects obtained (Figures 13(a2)-(d2)). Finally, the images of Figure 13(a1) (horizontal schlieren brightness gradient) and Figure 13(b1) (vertical schlieren brightness gradient) are processed by the above-mentioned conversion procedure “multi-path integration” to projection images of density thickness “with decreased noise level” (Figure 13(f)). It is again noteworthy that the new set of projections images of density thickness “with decreased noise level” can be obtained by using schlieren brightness gradient in horizontal and vertical or two perpendicular diagonal directions (Figure 13(f)). As seen in Figure 7(i1), Figure 13(f) the noise effects are reduced to an acceptable level and slightly present around initial noise spots.

4.3. CT-Reconstruction Results

By employing the density thickness images, noisy (Figure 13(a2)) and with decreased noise level (Figure 13(f)) as projections for CT-reconstruction, the 3D instantaneous density distributions of target flame have been successfully obtained for *z* = 8 mm to 40 mm. Figure 14 shows a sample set of one horizontal and one vertical cross-section of 3D-CT distributions of density. The darker parts are related to lower density areas with higher temperature and the brighter parts are related to higher density points with lower temperature. For better comparing new and former techniques, the corresponding density contour diagrams of each cross-section are depicted in Figure 15, as well. The density distributions images from multi-path integration in Figure 14(c2) and Figure 14(c3) are uniform, but the results from single-path integration are streaky Figure 14(c1) and Figure 14(c4). As shown in Figure 15(c4) the density contour diagrams from single-path integration are zigzagged form in contour lines.

Figure 16 provides the quantitative information of CT-reconstructed density distributions on the line of (y, z) = (0 mm, 24 mm). By applying the multi-path integration technique, the density distributions diagrams in the presence of noise have coincided with the absence of noise, but the noise effects clearly

Figure 14. Sample set of horizontal ((a) and (d)) and vertical ((b) and (e)) distributions of CT-reconstructed density distributions of target flame.

Figure 15. Corresponding density distributions contour diagrams of Figure 14.

Figure 16. Superposition of radial density distributions of CT-reconstructed density distributions of target flame for *y* = 0 and *z* = 24 mm.

visible on the single-path integration results. Results demonstrate the importance of the new technique for noise reduction.

5. Conclusions

We applied the Cauchy integration theorem to develop the noise reduction technique. Based on Cauchy’s integration theorem statement, contour integral depends only on the beginning and ending points of the paths and does not depend on the integration path. Therefore, the integration procedure in multi-path directions is employed. In the present study, a new technique “multi-path integration” is proposed for “noise reduction” in projections (density thickness) images of CT (computerized tomography). The schlieren brightness gradient information in two orthogonal directions (e.g. horizontal and vertical or two perpendicular diagonal directions) is required in the “multi-path integration” technique. The 20-directional quantitative schlieren optical system gives only images of schlieren brightness in the horizontal gradient and another 20-directional optical system seems necessary to obtain vertical schlieren brightness gradient, simultaneously. One suitable solution without the necessity of a new optical system is using the “inverse process”. By employing the new “inverse process” technique, the actual experimental data is used for reproducing other directions of the schlieren brightness gradient. This new approach is a simple, efficient and cost-effective solution and can be used for obtaining any optional directions of schlieren brightness gradient.

In this study, a schlieren image with information on the brightness gradient in two orthogonal directions is called a “complex schlieren image”. It has shown an important ability of complex schlieren brightness gradient, which is independence on the path and starting point in the integration process. The results validate that the “multi-path integration” novel technique has the capability of noise reduction in obtaining projections (density thickness) images of CT and the 3D density distributions of CT results.

Acknowledgements

This research work was partly supported by The Nitto Foundation and JSPS Grants-in-Aid for Scientific Research (C)16K06118. The supports are gratefully acknowledged.

Nomenclature

*B*: Brightness of schlieren image [-]

*B _{nf}*: Brightness of schlieren image in no-flame (no disturbance) condition in the test section [-]

*Dt*: Density thickness of deviation density [(kg/m^{3})(m)]; derived value

*Dt** ^{’}*: Density thickness [(kg/m

*Dt*^{*}: Density thickness of deviation density [(kg/m^{3})(m)]; actual value

*f*: Focal length of convergent lens in schlieren system [m]

*K*: Gladstone-Dale constant of air [m^{3}/kg]

*R*: Radius of reconstruction area [m]

*x*, *y*, *z*: Cartesian coordinates system of reconstruction volume [m]

*X*(*θ*)*, Y*(*θ*): Inclined coordinates by angle *θ* [m]

*ΔB*: Deviation brightness on schlieren image [-]* *

*Δs*: Transparent width of the light source [m]

*Δρ*: Deviation density [kg/m^{3}]; derived value

*Δρ*^{*}: Deviation density [kg/m^{3}]; actual value, (
$\Delta {\rho}^{*}={\rho}_{a}^{*}-{\rho}^{*}$ )* *

*θ*: Angle of observation [° (degree)]

*ρ*: Density [kg/m^{3}]; derived value, (
$\rho ={\rho}_{a}^{*}-\Delta \rho $ )

*ρ*^{*}: Density [kg/m^{3}]; actual value

${\rho}_{a}^{*}$ : Ambient density of air [kg/m^{3}]; actual value

*ρ _{ref}*: Reference density [kg/m

*φ*: Diameter [m]

References

[1] Settles, G.S. and Hargather, M.J. (2017) A Review of Recent Developments in Schlieren and Shadowgraph Techniques. Measurement Science and Technology, 28, 1-25.

https://doi.org/10.1088/1361-6501/aa5748

[2] Ishino, Y., Hayashi, N., Bt Abd Razak, I.F. and Saiki, Y. (2015) 3D-CT Measurement of Instantaneous Density Distributions of High-Speed Premixed Turbulent Burner Flames with a Multi-Directional Schlieren Camera (Effects of Flow Velocity on 3D Flame Front Shape). Proceedings of the 1st Thermal and Fluids Engineering Summer Conference, New York, 9-12 August 2015, 121-131.

https://doi.org/10.1615/TFESC1.cbf.012970

[3] Ishino, Y., Hayashi, N., Bt Abd Razak, I.F., Kato, T., Kurimoto, Y. and Saiki, Y. (2015) 3D-CT (Computer Tomography) Measurement of an Instantaneous Density Distribution of Turbulent Flames with a Multi-Directional Quantitative Schlieren Camera (Reconstructions of High-Speed Premixed Burner Flames with Different Flow Velocities). Flow, Turbulence and Combustion, 96, 819-835.

https://doi.org/10.1007/s10494-015-9658-5

[4] Ishino, Y., Hayashi, N., Ishiko, Y., Nazari, A.Z., Nagase, K., Kakimoto, K. and Saiki, Y. (2016) Schlieren 3D-CT Reconstruction of Instantaneous Density Distributions of Spark-Ignited Flame Kernels of Fuel-Rich Propane-Air Premixture. Proceedings of ASME Heat Transfer Summer Conference, Washington DC, 10-14 July 2016, HT2016-7423, 1-10.

https://doi.org/10.1115/HT2016-7423

[5] Ishino, Y., Hayashi, N., Ishiko, Y., Nagase, K., Kakimoto, K., Nazari, A.Z. and Saiki, Y. (2017) 3D Printing of Spark-Ignited Flame Kernels, Experimentally Captured by 3D-Computer Tomography and Multi-Directional Schlieren Photography. Heat Transfer, 139, HT-16-1726.

https://doi.org/10.1115/1.4035583

[6] Nazari, A.Z., Ishino, Y., Motohiro, T., Yamada, R., Ishiko, Y. and Saiki, Y. (2017) Schlieren CT Measurement of 3D Density Distributions of Flame Kernels of Spark-ignited Direct-Injection of Free, Cavity-Guided and Plane-Guided Fuel Jets. Proceedings of the 11th Asia-Pacific Conference on Combustion, Sydney, 10-14 December 2017, 1-4.

http://www.anz-combustioninstitute.org/proceedings.php

[7] Nazari, A.Z., Ishino, Y., Ishiko, Y., Ito, F., Kondo, H., Yamada, R., Motohiro, T.,
Miyazato, Y. and Nakao, S. (2020) Multi-Schlieren CT Measurements of Supersonic
Microjets from Circular and Square Micro Nozzles. Journal of Flow Control, Measurement
& Visualization, 8, 77-101.

https://doi.org/10.4236/jfcmv.2020.83005

[8] Nazari, A.Z., Ishino, Y., Yamada, R., Motohiro, T., Ito, F., Kondo, H., Miyazato, Y. and Nakao, S. (2018) CT (Computer Tomography) Measurement of 3D Density Distributions of High-Speed Premixed Turbulent Flames (Multi-Path Integration Image-Noise Reduction Technique Based on Novel Concept of Complex Brightness Gradient in Quantitative Schlieren Images). Proceedings of the 18th International Symposium on Flow Visualization, Zurich, 26-29 June 2018, 1-13.

https://doi.org/10.3929/ethz-b-000279202

[9] Nazari, A.Z., Ishino, Y., Yamada, R., Motohiro, T., Ito, F., Kondo, H., Miyazato, Y. and Nakao, S. (2019) CT (Computer Tomography) Measurement of 3D Density Distributions of Flame: Obtaining Vertical Gradient Schlieren Brightness from Horizontal Gradient for Image-Noise Reduction. Proceedings of the 4th Thermal and Fluids Engineering Conference, Las Vegas, NV, 14-17 April 2019, 185-194.

https://doi.org/10.1615/TFEC2019.cbf.027746

[10] Huang, C., Gregory, J.W. and Sullivan, J.P. (2007) A Modified Schlieren Technique for Micro Flow Visualization. Measurement Science and Technology, 18, N32-N34.

https://doi.org/10.1088/0957-0233/18/5/N04

[11] Collicott, S.H. and Salyer, T.R. (1994) Noise Reduction Properties of a Multiple-Source Schlieren System. AIAA Journal, 32, 1683-1688.

https://doi.org/10.2514/3.12160

[12] Chaloupka, J.L., Woods, M., Aas, J., Hutchins, J. and Thistle, J.D. (2014) Color Schlieren Imaging with a Two-Path, Double Knife Edge System. Optics Express, 22, 8041-8046.

https://doi.org/10.1364/OE.22.008041

[13] Barry, F.W. and Edelman, G.M. (1948) An Improved Schlieren Apparatus. Journal of the Aeronautical Sciences, 15, 364-365.

https://doi.org/10.2514/8.11588

[14] Hadjadj, A. and Kudryavtsev, A. (2005) Computation and Flow Visualization in High-Speed Aerodynamics. Journal of Turbulence, 6, N16.

https://doi.org/10.1080/14685240500209775

[15] Yokoi, T., Shinohara, H., Hashimoto, T., Yamamoto, T. and Niio, Y. (2000) Implementation and Performance Evaluation of Iterative Reconstruction Algorithms in SPECT: A Simulation Study Using EGS4. Proceedings of the Second International Workshop on EGS, KEK Proceedings, 200-20, 224-234.

http://rcwww.kek.jp/egsconf/proceedings/2iwoegs/yokoi.pdf