4 March 2015 Laser diode self-mixing interferometry for velocity measurements
Author Affiliations +
A diode laser velocimeter based on laser self-mixing has been developed and characterized as a reliable, precise, comparably cheap, and compact monitor. The resolution of this sensor at different incident angles and for a variety of solid and liquid targets moving at velocities between 0.1 and 50  m/s is presented. This includes a theoretical analysis of the underlying measurement principle, highlighting possibilities to extend the measurement capabilities to even higher velocities by altering the sensor design. Finally, an outlook on future applications of the sensor for detailed studies of supersonic gas jets used in beam diagnostics and atomic physics applications is given.



The field of self-mixing (SM) interferometry has been growing since the mid 1960s, with significant growth in the late 2000s due to better fabrication techniques of semiconductor lasers, which has opened up the possibilities to work with cheaper and higher-quality lasers in a wide range of applications. The SM effect leads to a precise and accurate metrology instrument for different purposes, such as distance measurements,1 vibration detection,2 and velocity characterization for solids,3 liquids,4 and plasmas.5

Gas targets are used for a number of accelerator based experiments, as well as for beam diagnostic purposes. In the case of supersonic gas jets, they combine low internal temperatures with high directionality, and as a result, they are very interesting as experimental targets in different fields of science.6 For the optimization and verification of the properties of such a monitor, it should be characterized with high accuracy, in particular with regard to its velocity and density profile. In these applications, gas jet velocities can be up to 2000m/s and can be inhomogeneously distributed across the jet. However, all currently used methods for such a characterization are either not reliable7,8 or require a powerful laser system.9

At the Cockcroft Institute, a new type of ionization beam profile monitor was devised and developed utilizing a supersonic gas jet. The monitor is based on a neutral gas, shaped into a thin sheet, a fraction of which is ionized when it crosses the beam. Gas ions, created in the jet, are then extracted, allowing two-dimensional beam profile imaging.10,11

The gas jet consists of a high density (109 to 1012particles/cm3) of neutral molecules, which could be molecular nitrogen, argon, or helium. The gas-jet curtain is expected to have a velocity range of 100 to 2000m/s, and the diameter of the jet itself is 1 to 20 mm. The characterization of the gas-jet curtain is important for further calibration of the beam monitor, so the velocimeter should measure the velocity in that range together with the density. Moreover, it is desirable to have a sensor that could be easily integrated into an existing setup with minimum additional optics components due to the complexity of the setup itself. For example, often used methodologies, such as particle image velocimetry,12 particle tracking velocimetry,13 and laser speckle velocimetry,14 which have a similar working principle and are based on analyzing the net displacement of moving particles, and interferometric and absorbing techniques15 are not applicable for use in a gas-jet setup due to the complexity of installing them into an existing system.

The SM scheme can be employed in nearly every situation a traditional interferometer (such as the Michelson or the Mach-Zender interferometers) can. The main advantage of the SM scheme with respect to traditional methods is the possibility of an unambiguous measurement from a single interferometric channel. Sensors based on SM compared to other techniques are compact, have a much lower cost, and are easier to align.

The SM effect implies that the light from the laser is reflected on the target back into the laser cavity. The phenomena can be referred to as an optical feedback, induced-modulated interferometer. In the case of semiconductor lasers, the theory is different from the SM effect in gas (CO2, HeNe) lasers where it was initially observed.16,17 Specifically, the nonlinear nature of the semiconductor active medium coupled to the optical gain and refractive index to the carrier density leads to different equations and solutions resulting in unique effects.

The physics of the SM technique is such that the precise measurement of small changes in the target position is straightforward to perform. As a result, the effect was mainly applied to precise and accurate measurements of vibrations, small displacements, and the study of slow movements. The measured velocities reported in the literature using the SM technique rarely exceed 5m/s (Refs. 3, 18, and 19) with several works mentioning the movement of solid object by 25m/s (Ref. 20) and 55m/s.15

SM interferometry has been applied to flow-velocity measurements often aimed at biological applications, such as blood flow characterization,21 Brownian motion, and biological species dynamic studies, widely presented in the review by Otsuka.4 However, all of the objects considered in these cases of velocity characterization move relatively slowly. Such a high interest in this area has resulted in a number of works with measurements of relatively small velocity from 0.1mm/s to 10cm/s (Refs. 22 and 23) with high accuracy.

Most work has been focused on the precision of such measurements providing a high spatial resolution. Using solid or liquid targets leads to a higher backscattered intensity as compared to gas targets; on the other hand, typical velocities of interest expected for a planar gas jet are at least in the range 0.5<M<4, where M stands for the Mach Number; when translated to m/s, this gives a range of roughly 100 to 2000m/s,24 which has previously not been explored.

The motivation of the present study is to clarify the possibility of measuring high-speed velocities up to 50m/s in the geometry of the gas-jet setup, Secs. 3.1 and 3.2, with solid objects, and the possibility of measuring fluid velocity up to 100cm/s, see Sec. 3.3.



The first analysis of the SM phenomena in the case of a laser diode (LD) was performed by Lang–Kobayashi.25 A useful review of the theory and practical approach to the SM effect in semiconductor lasers is given in the book,18 demonstrating theoretical investigations into SM and developments of various applications.

The working principal of the system is that the light from the laser scatters off the target and the scattered light returns back into the laser cavity. The laser plays the role of a coherent heterodyne receiver and amplifier of the signal, and at the same time, the quadratic detection of the registered signal is inside the laser. The behavior of the system can be described in different ways. The target can be seen as an additional mirror for the laser system, and then the system is a known three-mirror laser, where the laser has an external cavity with a length of the distance between the laser and the target, often called the theory of the equivalent-lens cavities.26 From this point of view, the SM system can be described by a set of field and rate equations obtained by adding an additional external feedback term to the standard equations.25,27 Another approach is to consider the system as a two-mirror laser with a modified complex reflection coefficient or effective reflection coefficient, and use solutions derived for LD systems.3 As a result, the solution is based on replacing the parameters without feedback with parameters that depend on an additional mirror (target).28

Within both ways of describing the system, the modulation of light as a result of the interaction with the target can be mathematically described by



P0 and P are the optical powers without and with feedback, m is a modulation parameter, F is the self-mixing interferometric waveform, τD is the time of flight of the wave, and υ is the frequency of the laser light undergoing feedback. The modulation index m and the shape of F are the function of the feedback parameter Cr3D, so that C depends both on the level of feedback (the portion of light that returns back to the laser cavity), which is defined by the reflectivity coefficient r3 together with coupling optics, and on the target distance D.3,27,28 There are four different types of signals in SM depending on the amount of backscattered light coupled into the cavity, see Fig. 1. In the case of a low amount of back-injected light (weak feedback, C1), the signal is a sinusoidal and the function F is cosine. Gradually increasing the level of feedback, the shape of the signal can be deformed from a sinusoidal (C1) to a nonsymmetrical sinusoidal-like shape (0.1<C<1) and then into an asymmetrical sawtooth-like form (C>1). The orientation of the function F depends on the sign of the phase, so the direction of the target movements can be discriminated.2930.31 Figure 1 demonstrates the example of the SM signals together with the principal of direction determination. In both Figs. 1(a) and 1(b), the two top plots represent the case with the targets approaching the laser, and the two bottom plots represent that of the target moving away from the laser. The signals in the middle with C<1 are close to a sinusoidal shape.

Fig. 1

Experimentally obtained self-mixing signals for two different types of surfaces: (a) mirror and (b) white paper when the translational stage was moving with a velocity of 5.00cm/s toward (top plots) and away (bottom plots) from the laser. The plots are shifted along the y axis, so the comparison of the shapes of self-mixing (SM) signals is possible as a function of the level of the feedback re-entering the cavity. The coupling efficiency is related to the coefficient C. The plots in the middle with C<1 correspond to weak feedback, and the shape of the signal is close to sinusoidal. As the amount of light coupling back into the cavity increases, the shape of the signal distorts and becomes a nonsymmetrical sinusoidal signal with C1. In the case of strong feedback (C>1), which is shown in the top and bottom plots in both (a) and (b), the signal is sawtooth-like. In the case of an asymmetrical shape of the signal, the direction of leaning denotes the direction in which the target was moving. The top plots lean to the left corresponding to when the target is moving toward the laser diode (LD). The bottom plots, where the signal is leaning to the right, are when the target was moving away from the laser. In case of the white paper (b), the signal amplitude is slightly reduced and distorted in shape as a result of speckle effects, though the shape and trends are the same with the mirror (a) as a target.


The velocity can be calculated from the spectrum of the power of the signal. As previously mentioned, the signal is the result of mixing the initial lasing field and the Doppler shifted backscattered light. The power spectrum of the laser changes due to interaction between the lasing field and the small backscattered field, which re-enters the laser cavity with results described by Eq. (1). The Doppler shift of the periodic optical power fluctuation (f) is equal to


where n is the unit vector in the direction in which the harmonic wave of light propagates, λ is the wavelength of the lasing light, and V is the velocity of the target. Equation (2) shows the shift frequency at the moment of scattering. In general, both parameters f and V depend on time and the signal is a frequency modulated signal.32

In the case of mirrors or other surfaces with low roughness, the specular (mirror-like) reflection prevails over the diffuse reflection. As a result, the level of the feedback in the case of alignment is very high and leads to cases with C>1. Using additional filters, it is possible to have a system with weak feedback.

Using a rough surface as a target for SM, it is possible to apply the SM technique without additional alignment, so it can be performed at any angle. Light scatters off a rough surface in all directions; this is in contrast to smooth surfaces where light reflects mainly in one direction. However, due to random interference of light scattered off the rough surface, speckles33 are observed affecting the shape of the SM signal received, causing an error in the measurement.34 An ideal case of an SM signal is a harmonious function with constant amplitude with a narrow spectral line in the spectrum. In the presence of speckles, the amplitude and the phase of the SM signal is modulated, which leads to spectral broadening and additional noises in the signal.35

The speckle phenomena in SM can be used not only for a displacement measurement,36 but also for the surface classification37 as well as for the length of the moving objects.38

When the light travels through a liquid or gas, the interaction between the light and particles randomly changes the phase and amplitude of the scattered light.33 The power spectrum [Eq. (1)] provides information about the scattered particles, such as velocity, along the laser light axis and other properties, such as the concentration in the volume and their size.39


Experimental Results

The SM technique can be applied using any type of laser, and both the target and property of the laser itself should be correctly taken into account. In the case presented here, the experiments were performed using an LD with 650 nm wavelength (L650P007) with a power of <5mW, driven by a constant current. The light was focused onto the target by a lens with numerical aperture 0.26. The signal was received from the built-in photodiode, which is usually used for monitoring the LD. The current was amplified and modified into a voltage signal by a transimpedance amplifier produced by Femto (DHPCA 100). The velocities of the targets were calculated by performing a fast Fourier transformation (FFT) of the voltage signal applying Hamming and Blackman-Harris windows. The spectrum was processed with a smoothing filter, and the peak of the spectra was found by fitting a Gaussian distribution in the frequency domain.

The following subsections are divided into three parts, with each subsection devoted to a different investigated target: a mirror with a high reflective index mounted on a translational stage, a diffusive surface (white paper) first mounted on a translational stage and then on a rotational disc, see Figs. 2(a) and 2(b), and, finally, milk diluted in water and also a titanium dioxide (emulsion paint) water solution measured by means of laminar flow out of a hose, see Fig. 2(c). The right part of Fig. 2 demonstrates the projection of the targets’ velocity in different experimental setups. Only that component of the velocity that lies on the laser axis interacts with the laser light.

Fig. 2

Illustration of the layout of the experiments on the left side, and the schematic illustration of the projection of the velocity on the laser axis for each case on the right side. The light from the LD is focused by means of multiple lenses onto a target. Light then scatters off the target and is coupled back into the laser cavity, leading to an SM effect. A signal is produced by an inbuilt photodiode (PD), and amplified and converted into a voltage signal by an amplifier. Three experimental setups were established for testing the SM velocimetry technique: (a) the translation stage; (b) rotating disc; and (c) fluid stream. The translation stage (a) moves the target (mirror or white paper) toward and away from the laser; as a result, the scattered light has a Doppler shift proportional to the velocity vector in the direction of movement. In the case of rotating disc (b) covered with white paper, the direction of the rotation defines the direction of the velocity vector relating to the laser axis. The velocity vector has its projection on the laser axis depending on the incident angle. The fluid (c) disgorges from the nozzle, so the flow is laminar. The laser light shines at an angle to the velocity [αf in (c)], so the measured velocity is the velocity projected onto the laser axis.



Moving Target with Different Reflectivity

First, the experiments with a mirror and with white paper were performed using the same setup, where the target moves back and forth on the laser light axis, see Fig. 2(a). The target was moved by a stage driven by a direct current motor controller, which gives precise velocities in the range of 0.01 to 10.00cm/s. In the case of the mirror, all four feedback regimes, mentioned earlier, were observed and studied depending on how focused/defocused the light was on the target. Figure 1(a) demonstrates different types of SM signals and the distortion of the shapes of the signal depending on the amount of light coupling back into cavity. In the case of moderate (C>1) and weak feedback (0.1<C<1), both the velocity and its direction were measured with high accuracy. The reflectivity coefficient directly influences the feedback parameter and, thus, the performance of the SM technique. The refractive index depends on the real refractive index and light absorption coefficient. Reflectivity is a function of the refractive index and incident angle. The reflectivity of the aluminum mirror at the normal angle for 650 nm is 92%. In order to reduce the amount of light coupling into the laser cavity, density filters were used. Experiments with a mirror, which moves across the transverse axis of the laser beam, were performed to understand the SM technique as a whole, such as the best geometry of the setup, reduction of the noise level, the focusing/defocusing properties of the LD, and the investigation of the system itself. In the case of a well-adjusted setup, the mirror target was moved up to 90 cm from the laser without changing the shape of the signal [Fig. 1(a), regimes with C>1]. The experiments demonstrated that accurate velocity measurements with a precision of >1% can be obtained at a distance up to 90 cm from the target.

With a rough surface target, such as paper, weak or moderate feedback is expected due to a lower reflectivity equal to 60% when normal to the surface for laser light of 650 nm. Typical examples of the signals received with white paper are shown in Fig. 1(b). The signals are similar in shape and the precision of the methods are the same in both cases. However, the amplitude of the signal is less when compared with the same type of signals in the mirror experiments. The size of the laser spot on the white paper required readjustments while the distance between the LD and target was changing. Using the mirror as a target allowed measurements up to 90 cm from the target without changing anything in the setup; however, when using the white paper as a target, readjustments were required periodically up to 50 cm where no signal was detected.

Accurate measurements with a >1.5% precision of the velocity were obtained in the range of 0.01 to 10.00cm/s at a distance up to 50 cm between the LD and the white paper target. The results of the measurements are presented in Fig. 3, which shows both the accuracy and precision of the measurements. The thin straight solid line is a reference velocity, and the measured velocity is shown with associated error bars, which give the range within which the measured velocities were obtained over multiple readings. The dashed lines represent the measured velocity in the SM signal at different distances up to 90 cm from the laser with a mirror as a target, and the solid line shows a measured velocity using the white paper as a target.

Fig. 3

Example of measurements from the experiments using a movable stage, which allows for a precise velocity measurement (straight solid line) in the range of 0.01 to 10.00cm/s with a mirror (dot lines) and white paper (thick solid lines) with an error bar obtained at different distances between the laser and the target. The velocity measurements were carried out up to 90 cm for the mirror and 50 cm for the white paper and showed a better than 1% accuracy and 1.5%, respectively. The error bars indicate the range of the measured velocity obtained in different independent sets of experiments.



Rotating Target with Low Reflectivity (60% at 650 nm)

A rotating disc was then covered with white paper, allowing the limitations of the setup to be investigated at high velocities. The distance between the target and the laser was fixed to 10 cm in the geometry shown in Fig. 2(b). This distance was selected as this is the expected distance when the velocimeter will be installed into the gas-jet setup. The rotational speed ω of the disc can vary from 2π to 300πrad/s which translates into tangential velocities V from 1cm/s to 100m/s.

The white paper has micrometer (40 to 50μm) surface roughness. The power spectrum of the light scattered off rough surfaces is anticipated to be Gaussian with a frequency corresponding to the Doppler shift of the SM system. The reflectivity properties of the white paper are such that 60% of the light, with a wavelength of 650 nm, is expected to scatter. The amount of light scattered in different directions depends on the incident light angle with a maximum at 90 deg.40

Different points on the rotating disc move at different tangential velocities V. The motion can be characterized by the rotational speed, which is given by revolutions per second, ν. The tangential velocity V is perpendicular to the diameter of the disc at each point and at a tangent to the circle. The tangential velocity of the point at the radius R of the disc is V=ωR, where ω is 2πν. The laser is turned to an angle α in order to have a vector of velocity, which would not be equal to zero in relation to the laser beam axis. The right part of Fig. 2(b) illustrates this concept, showing the projection of the velocity on the disc depending on whether the rotation is clockwise or anticlockwise when the laser spot is above the point of rotation, i.e., the center of the disc.

The Doppler shift for the SM system can be described by Eq. (2), where the velocity vector of V and the unit vector n can be described differently depending on the position of the laser spot.

The co-ordinate system is defined to be x-y on the surface of the disc with the origin at the center, and z is parallel to the normal vector of the disc surface. In Fig. 4(a), for vectors V=V(0,1,0) and n=(0,cosα,sinα), the Doppler shift is equal to f=2ωRcosα/λ. In the case when the laser is shining on any point on the disc surface [see Fig. 4(b)], the vector V is given by V=V(sinβ,cosβ,0). Thus, the Doppler shift in the most general case is equal to



Fig. 4

The laser light from the LD is focused by means of multiple lenses onto a point on a rotating disc located a distance R from the point of rotation. The origin of x-y-z coordinate system is defined to be at the point of rotation of the disc. n is a unit vector in the direction of propagation of light, and ω is the rotational speed of the disc. (a) In the case of the laser shining onto any point on the x axis, the Doppler shift depends only on the angle α between laser light direction and the disc surface (y axis). (b) At any other point, the Doppler shift depends on the angle β between the velocity vector and y axis.


The variation of the velocity can be achieved by changing the revolution per second (ν), the position of the laser spot on the disc (which will change the radius R), and the angle α.

In order to validate the results of the measured SM system, a photointerrupter circuit based on the OPT101 photodiode amplifier from Texas Instruments and a microcontroller board (ATMega 328) were used for measuring the rotational speed of the motor.

From every measured signal, a corresponding frequency spectrum can be found, and by using Eq. (1), the velocity can be calculated with a reference velocity based on Eq. (2). The uncertainty of the radius ΔR is defined as the size of the laser spot on the paper; all other errors are limited by the instrumentation. The relative error on the reference velocity can be calculated by differentiating Eq. (2) with the angle variation between the laser and the disc giving the most significant uncertainty. The experimental uncertainty of the velocity, based on the uncertainties in the experimentally measured quantities, was calculated according to δVlaser=(δR)2+(δν)2+(δα)2, where δR=ΔR/R, δν=Δν/ν, and δα=tan(α)Δα. The calculated dependency on the velocity uncertainty of the setup from the angle is shown in Fig. 5 with a solid black line with a 5% level of the accuracy by the dashed line. For calculations, it was assumed that the angle could be measured with a 0.5 deg precision. The rotational velocity ω was defined with an accuracy of 1%. The size of the laser spot on the paper depends on the angle α as well being proportional to the cosine of this angle, leading to an increasing relative error of the measurements at small, <10deg, angles due to a broadening of the frequency spectrum. Speckle modulation increases dramatically at angles >80deg, which also results in frequency spectrum broadening.19,35 Considering the trend of the calculated experimental uncertainties, presented in Fig. 5 by a black solid line, the optimal range of angles from an accuracy point of view is between 13 and 77 deg, which gives an error of <5%, and between 25 and 63 deg, which gives an error <3%.

Fig. 5

The theoretical and experimental relative error dependence on the angle between the disc covered with white paper and the laser. The theoretical error was calculated as the uncertainties in the experimentally measured quantities (see Sec. 3.2) and shown by the solid line. The experimental accuracy is shown by the dashed line with round dots, which represents the mean accuracy calculated from multiple readings over a range of velocities. The error bars on the experimental plot represent the range of the measured velocities obtained after set of independent experiments and demonstrate the precision of the results. The experimental results follow the same trend as the theoretical one. In the range of 15 to 85 relative errors were better than 5%, and the most precise measurements were received in the range of 30 to 75.


During the experiment, the following parameters were varied: the position of the laser spot on the disc, the angle α, and the rotational speed of the disc. Figure 5 represents the results of the accuracy, which was achieved during the experiments for different angles, shown by a dashed line. The accuracy of the measurements was defined as the relative error in relation to the reference velocity. For each angle taken in the plot, multiple measurements were taken over a velocity range of 0.5 to 30m/s with at least 10 measurements for each velocity. The points in Fig. 4 represent the average of the uncertainty over the velocity range for each angle. Along with the accuracy of the measurements, the error bars are presented for each point, which represents the precision of the independent measurement within the total accuracy of the results.

The experimental accuracy follows the same trend as the theoretical one in Fig. 5 with sharp improvements of the accuracy starting from an angle of 15 deg with a plateau of 3 to 4% of the accuracy in the range of 25 to 80 deg and with a rise of inaccuracy above 85 deg. The error bar for each point in the plot shows that the precision of the measurements improves within the range of 25 to 75 deg. The calculated velocity is based on the frequency extracted from the signal spectrum. In the case of a movable mirror, the whole surface moves with the same velocity. As a result, the peak of frequency is quite narrow, even if the spot of the laser light is nonpoint-like. Each point on the disc, however, has a different velocity. The laser spot has a finite size, i.e., it is nonpoint-like; this means that for each point on the disc, a different tangential velocity is measured, which results in a distribution of velocities with different weights in the signal measured. This can be seen by a broadening of the signal spectrum.

Figure 6 demonstrates typical waveforms and their FFTs in the case of an accuracy <5% [Fig. 6(a), 1.4% accuracy] and >5% [Fig. 6(b), 6.7%]. The signals were chosen such that the peaks of the FFTs are relatively close to each other for comparison. The Doppler shifts are 3.41 MHz for Fig. 6(a) and 3.23 MHz for Fig. 6(b), which correspond to tangential velocities V of 17.6 and 10.6m/s, respectively, while the angle between the laser axis and the surface of the disc [see Fig. 2(b)] was 51 and 9 deg, respectively. The disc, when covered with white paper, leads to an intensity modulation due to speckle effects, which can be seen in the signal plots in Fig. 6. The signals in Fig. 6 have a shape typical for weak feedback. The setup allows calculation and verification of the SM method with velocities of up to 50m/s. The only limitation of the experiment found so far is from the electrical circuit parameters (bandwidth of the amplifier). As a result, the high velocities could be measured only with a small angle and, as a consequence, with low accuracy. The experimental results on the velocity measurements with white paper as a target are shown in Fig. 7. Using a disc as a target allows the SM velocimeter to be validated. The accuracy is dominated by the angle between the target and the laser, see Fig. 5. The accuracy of the reference velocity was obtained based on the uncertainties of variables used for calculations. The measured velocity shown in Fig. 7 is an average of repeated measurements with a combined error. Starting from 20m/s, only large values of α were used, close to the normal of the disc surface, due to the low bandwidth of the transimpedance amplifier in the readout electronics. The precision of the results are, on average, better than 6%.

Fig. 6

Examples of typical SM signals and their fast Fourier transformations received from the setup with a rotating disc for the case of a relative error (a) <5% and equal to 1.4%, and (b) >5% and equal to 6.7%. The reference and measured velocities were (a) 17.6m/s and (b) 10.6m/s. The angles between the velocity vector and laser axis were (a) 51 deg and (b) 9 deg, leading to Doppler frequencies of 3.4 and 3.23 MHz, respectively. The waveforms demonstrate the intensity modulation effect due to nonpoint-like laser spot and speckle effects leading to a broadening of the frequency spectrum. In (b), the spectrum is expected to be broader due to the changing of the size of the laser spot on the disc.


Fig. 7

Experimental results from measurements of a rotating disc covered with white paper. The difference in the resulting and estimated velocity is within the limits of uncertainty. The limitation of the velocity was due to the bandwidth of the amplifier used for the SM velocimeter. The increasing uncertainty at high speeds, equal to and higher than 20m/s, is caused by the angle, 90deg, used for the measurement.



Fluid as a Target

The next step in the experiment was to check the possibility of measuring the velocity in the case of smaller scale targets in the case of low-level feedback.

To measure the velocity of liquids, an experimental setup was established which allows measurements to be compared with analytical and numerical solutions for verification and benchmarking of the SM velocimeter.

Laminar flow was achieved by discharging liquid from a nozzle of a tube pointed horizontally into the atmosphere. The liquid flow outlet velocity from the nozzle depends on the pump power. The laser light was focused close to the nozzle, 0.4 cm from the outlet; and the direction of the velocity, αf in Fig. 2(c), was at an angle of 28 deg with respect to the laser light axis as seen in Fig. 2(c). The right part of the schematic configuration represents the position of the velocity vector with respect to the laser axis.

A mixture of water with milk was used for fluid velocity measurements at first. Milk was chosen as it consists of an emulsion of fat in water, where 87% of the milk itself is water, and only 4% is fat globules. The refractive index for the milk at 650 nm wavelength is 1.45 and the light absorption coefficient is <105,41 which resulted in reflectivity of 3.5%.

In the first set of experiments, the velocity of the diluted milk was measured in the range between 1 and 20cm/s. The distance between the LD and the target was fixed to 4 cm. The measurements were performed for 5% concentration of milk. The fat globules’ shape is close to spherical, with the average size of the molecules being 3 to 6μm.41

In the second set of experiments, the velocity of water seeded with titanium dioxide was measured with the same wavelength of 650 nm. A refractive index of 1.49 was calculated for the mixture based on a measurement of the critical angle from a measurement of the angle of total internal reflection. The setup was used for studying smaller-scale seeds, 0.5 to 1μm, allowing lower densities of the concentration, <1% of the volume, to be probed.

Since both the milk and titanium dioxide solutions appeared white, the visible light is scattered equally. Light scatters off moving particles in the fluid (the initial particles or seeds) in all directions with some angular distribution of scattering efficiencies corresponding to the properties of the scattered particles.42

The SM technique allows measurements of the velocities of the solution of up to 1m/s. An example of the SM signal obtained during the experiment can be seen in Fig. 8. The signal can be characterized as a signal with weak feedback. Depending on the optics used and the concentration of the seeds, the direction of the moving target was able to be defined only when the type of the SM signal corresponded to the weak feedback with C>0.1.

Fig. 8

Experimental setup for flow velocity measurements. The liquid flow discharges from the nozzle into the atmosphere using a controlled pump system, which allows a reference velocity to be defined. The laser light was focused on the center of the liquid flow at a distance of 0.4 cm from the nozzle outlet at 28 deg.


The velocities of fluids up to 1m/s were successfully measured with high accuracy, as shown in Fig. 9(b). The velocity was found from the averaged spectrum of 20 power spectra for each velocity on the plot. The dashed line is the reference velocity of the fluids. The differences in the resulting and estimated velocity of the fluid were within the limits of uncertainty of measurement and calculations.

Fig. 9

Experimental results from measuring the velocity of fluids. (a) Example of the signal received from the setup with a titanium oxide water solution as the target. The reference and measured velocity was 0.55m/s. (b) The results of the measurements of the velocity up to 1m/s were measured with an accuracy better than 4%. Velocities in the range of 1 to 20cm/s were carried out for milk diluted in water to a 5% concentration, shown by round dots. Higher velocities, up to 1m/s, shown in the graph are obtained from measurements of the titanium oxide water solution with a <1% concentration in volume and were also measured with accuracy >4%.



Conclusion and Outlook

A diode laser velocimeter based on the laser SM method has been developed as an easy-to-build and compact solution when compared to alternative measurement techniques. First tests with solid targets and fluids were carried out and the technique seems a promising way for a complete characterization of moving objects. The laser velocimeter could provide in-detail information about velocity and density of liquids and, eventually, a gas jet with micrometer resolution. Such a sensor would allow for unambiguous measurements from a single interferometric channel and could be installed even in radiation-exposed environments.

The laser velocimeter is a self-aligning device based on the SM method, where the laser is both a transmitter and receiver of the signal. It should be noted that laser SM is usually used for measurements of low velocities and vibrations. The theoretical analysis presented here shows the possibility to extend these measurement capabilities to higher velocities by altering the design.

Table 1 presents the results of the velocity measurements for different targets. In the experiments, the bandwidth of the transimpedance amplifier was 15 MHz, which was the main limitation in measuring higher velocities. Measurements of high, up to 50m/s, velocities with white paper on a disc were achieved with an uncertainty better than 5%. The fluid velocity measurement showed the possibility to measure a signal with a low-level feedback with promising results of measuring high velocities, up to 1m/s, which were not covered in previous studies.

Table 1

Results of the experiment with different targets with different reflectivity.

Target, n reflectivity for 650 nmVelocity, m/sError in the measurement
Mirror, n=92%Up to 0.1 (limitation from the moving setup)<1%
White paper, n=60%Up to 0.1 (limitation from the moving setup)<1.5%
White paper on a disc, n=60%Up to 50 (limitation from the electronic part of setup)<5%
Fluid, milk diluted in water at a concentration of 5% in volume, n=3.5%Up to 0.2 (limitation is under investigation)<5%
Fluid, titanium dioxide water solution at a concentration of <1% in volume, n=4%Up to 1 m/s (limitation is under investigation)<4%

An initial calibration and verification of the accuracy of the method was performed on a setup that gave access to velocities of up to 50m/s. In order to verify the possibility of using the velocimeter as a flow sensor, measurements of local velocities were performed on fluids with externally imposed flow profiles and different densities. The first experiments allowed velocity measurements of fluids up to 1m/s.

First experimental results show that the velocity of a white paper target and fluids can be measured with an accuracy of >5% and micrometer spatial resolution over a velocity range from 0.5mm/s to 50m/s. The sensor’s capability to measure even higher target velocities will be subject to further studies.

For the next step, it is planned to apply the above method to characterize liquids with different densities and velocity distributions, as well as ultrasonic gases. The extension of the methodology to these particular cases will have to overcome two major challenges. The first is that the intensity of the backscattered light may not be large enough, and second, the sensor may not be able to cope with velocities on the order of the speed of sound. Taking into account optical properties and density of the gas jet together with collimation optics for the LD, calculations show that, potentially, seeding techniques will need to be applied to increase the reflectivity of the supersonic jet, and this will be part of future investigations.

Successful use of SM strongly depends on the quantity of feedback coupled back into the laser cavity. It is stated15 that for the weak feedback regime to be attained, a level of feedback of 90dB is usually required. The required level of feedback also depends on the particular type of laser used and on the particular optics. The quality of the detectable signal also depends on the electrostatic insulation of the system and the pickup noise of the electronics. It is, therefore, necessary to set up an experiment that allows the attenuation in the external cavity leading to the different regimes to be measured, using the particular elements and external conditions of the specific laboratory. One particular area of interest in this study is to find the minimum level of concentration (the number of particles) required to be detected together with the maximum speed that can be measured. This will provide insight into the applicability of the method for targets with lower density.


This project has received funding from the European Union’s Seventh Framework Programme for research, technological development and demonstration under Grant Agreement No. 289191. Work supported by HGF and GSI under contract VH-NG-328 and STFC under the Cockcroft Institute Core Grant No. ST/G008248/1. The authors would like to thank Dr. Massimiliano Putignano for his help with the initial setup, as well as Professor Gaetano Scamarcio and his group from the University of Bari for their introduction to the self-mixing technique.


1. S. Ottonelli et al., “Laser-self-mixing interferometry for mechatronics applications,” Sensors 9(5), 3527–3548 (2009).SNSRES0746-9462 http://dx.doi.org/10.3390/s90503527 Google Scholar

2. L. Scalise et al., “Self-mixing laser diode velocimetry: application to vibration and velocity measurement,” IEEE Trans. Instrum. Meas. 53(1), 223–232 (2004).IEIMAO0018-9456 http://dx.doi.org/10.1109/TIM.2003.822194 Google Scholar

3. T. Bosch, N. Servagent and S. Donati, “Optical feedback interferometry for sensing application,” Opt. Eng. 40(1), 20–27 (2001).OPEGAR0091-3286 http://dx.doi.org/10.1117/1.1330701 Google Scholar

4. K. Otsuka, “Self-mixing thin-slice solid-state laser metrology,” Sensors 11(2), 2195–2245 (2011).SNSRES0746-9462 http://dx.doi.org/10.3390/s110202195 Google Scholar

5. I. J. Rasiah, “Improved Ashbya-Jephcott interferometer for temporal electron density measurements in plasmas,” Rev. Sci. Instrum. 65(5), 1603–1605 (1994).RSINAK0034-6748 http://dx.doi.org/10.1063/1.1144899 Google Scholar

6. J. Ullrich et al., “Recoil-ion and electron momentum spectroscopy: reaction-microscopes,” Rep. Prog. Phys. 66(9), 1463 (2003).RPPHAG0034-4885 http://dx.doi.org/10.1088/0034-4885/66/9/203 Google Scholar

7. S. C. Morris and J. F. Foss, “Transient thermal response of a hot-wire anemometer,” Meas. Sci. Technol. 14(3), 251 (2003).MSTCEP0957-0233 http://dx.doi.org/10.1088/0957-0233/14/3/302 Google Scholar

8. M. Nishida, T. Nakagawa and H. Kobayashi, “Density contours of a supersonic freejet,” Exp. Fluids 3(3), 181–183 (1985).EXFLDU0723-4864 http://dx.doi.org/10.1007/BF00280458 Google Scholar

9. M. Raffel et al., Particle Image Velocimetry: A Practical Guide, pp.259–388, Springer, New York (2007). Google Scholar

10. V. Tzoganis and C. P. Welsch, “A non-invasive beam profile monitor for charged particle beams,” Appl. Phys. Lett. 104(20), 204104 (2014).APPLAB0003-6951 http://dx.doi.org/10.1063/1.4879285 Google Scholar

11. V. Tzoganis, A. Jeff and C. P. Welsch, “Gas dynamics considerations in a non-invasive profile monitor for charged particle beams,” Vacuum 109, 417–424 (2014).VACUAV0042-207X http://dx.doi.org/10.1016/j.vacuum.2014.07.009 Google Scholar

12. R. J. Adrian and J. Westerweel, Particle Image Velocimetry, pp.165–240, Cambridge University Press, New York (2011). Google Scholar

13. N. A. Buchmann et al., “Ultra-high-speed 3D astigmatic particle tracking velocimetry: application to particle-laden supersonic impinging jets,” Exp. Fluids 55(11), 1842 (2014). http://dx.doi.org/10.1007/s00348-014-1842-1 Google Scholar

14. N. A. Fomin, Speckle Photography for Fluid Mechanics Measurements, pp.129–148, Springer, New York (1998). Google Scholar

15. G. Giuliani et al., “Laser diode self-mixing technique for sensing applications,” J. Opt. A Pure Appl. Opt. 4(6), S283–S294 (2002).JOAOF81464-4258 http://dx.doi.org/10.1088/1464-4258/4/6/371 Google Scholar

16. M. J. Rudd, “A laser Doppler velocimeter employing the laser as a mixer-oscillator,” J. Phys. E Sci. Instrum. 1(7), 723 (1968). 0022-3735 http://dx.doi.org/10.1088/0022-3735/1/7/305 Google Scholar

17. J. H. Churnside, “Laser Doppler velocimetry by modulating a CO2 laser with backscattered light,” Appl. Opt. 23(1), 61–66 (1984).APOPAI0003-6935 http://dx.doi.org/10.1364/AO.23.000061 Google Scholar

18. D. M. Kane and K. A. Shore, Unlocking Dynamical Diversity: Optical Feedback Effects on Semiconductor Lasers, John Wiley & Sons, Chichester (2005). Google Scholar

19. W. Huang et al., “Effect of angle of incidence on self-mixing laser Doppler velocimeter and optimization of the system,” Opt. Commun. 281(6), 1662–1667 (2008).OPCOB80030-4018 http://dx.doi.org/10.1016/j.optcom.2007.11.029 Google Scholar

20. X. Raoul et al., “A double-laser diode onboard sensor for velocity measurements,” IEEE Trans. Instrum. Meas. 53(1), 95–101 (2004).IEIMAO0018-9456 http://dx.doi.org/10.1109/TIM.2003.821483 Google Scholar

21. M. Norgia, A. Pesatori and L. Rovati, “Self-mixing laser Doppler spectra of extracorporeal blood flow: a theoretical and experimental study,” IEEE Sens. J. 12(3), 552–557 (2012).ISJEAZ1530-437X http://dx.doi.org/10.1109/JSEN.2011.2131646 Google Scholar

22. M. Nikolic et al., “Self-mixing laser Doppler flow sensor: an optofluidic implementation,” Appl. Opt. 52(33), 8128–8133 (2013).APOPAI0003-6935 http://dx.doi.org/10.1364/AO.52.008128 Google Scholar

23. L. Campagnolo et al., “Flow profile measurement in microchannel using the optical feedback interferometry sensing technique,” Microfluid. Nanofluid. 14(1–2), 113–119 (2013).1613-4982 http://dx.doi.org/10.1007/s10404-012-1029-0 Google Scholar

24. M. Putignano and C. P. Welsch, “Numerical study on the generation of a planar supersonic gas-jet,” Nucl. Instrum. Methods Phys. Res. A 667, 44–52 (2012).NIMRD90167-5087 http://dx.doi.org/10.1016/j.nima.2011.11.054 Google Scholar

25. R. Lang and K. Kobayashi, “External optical feedback effects on semiconductor injection laser properties,” IEEE J. Quantum Electron. 16(3), 347–355 (1980).IEJQA70018-9197 http://dx.doi.org/10.1109/JQE.1980.1070479 Google Scholar

26. W. T. Silfvast, Laser Fundamentals, pp.402–433, Cambridge University Press, New York (2004). Google Scholar

27. G. A. Acket et al., “The influence of feedback intensity on longitudinal mode properties and optical noise in index-guided semiconductor lasers,” IEEE J. Quantum Electron. 20(10), 1163–1169 (1984).IEJQA70018-9197 http://dx.doi.org/10.1109/JQE.1984.1072281 Google Scholar

28. W. M. Wang et al., “Self-mixing interference inside a single-mode diode laser for optical sensing applications,” J. Lightwave Technol. 12(9), 1577–1587 (1994).JLTEDG0733-8724 http://dx.doi.org/10.1109/50.320940 Google Scholar

29. E. T. Shimizu, “Directional discrimination in the self-mixing type laser Doppler velocimeter,” Appl. Opt. 26(21), 4541–4544 (1987).APOPAI0003-6935 http://dx.doi.org/10.1364/AO.26.004541 Google Scholar

30. S. K. Ozdemir, S. Shinohara and H. Yoshida, “Effect of linewidth enhancement factor on Doppler beat waveform obtained from a self-mixing laser diode,” Opt. Rev. 7(6), 550–554 (2000).1340-6000 http://dx.doi.org/10.1007/s10043-000-0550-7 Google Scholar

31. Y. Fan et al., “Improving the measurement performance for a self-mixing interferometry-based displacement sensing system,” Appl. Opt. 50(26), 5064–5072 (2011).APOPAI0003-6935 http://dx.doi.org/10.1364/AO.50.005064 Google Scholar

32. G. Plantier et al., “Real-time tracking of time-varying velocity using a self-mixing laser diode,” IEEE Trans. Instrum. Meas. 53(1), 109–115 (2004).IEIMAO0018-9456 http://dx.doi.org/10.1109/TIM.2003.821488 Google Scholar

33. J. W. Goodman, Speckle Phenomena in Optics: Theory and Applications, Roberts & Company, Englewood, Colorado (2007). Google Scholar

34. A. Magnani and M. Norgia, “Spectral analysis for velocity measurement through self-mixing interferometry,” IEEE J. Quantum Electron. 49(9), 765–769 (2013).IEJQA70018-9197 http://dx.doi.org/10.1109/JQE.2013.2273562 Google Scholar

35. R. Kliese and A. D. Rakić, “Spectral broadening caused by dynamic speckle in self-mixing velocimetry sensors,” Opt. Express 20(17), 18757–18771 (2012).OPEXFF1094-4087 http://dx.doi.org/10.1364/OE.20.018757 Google Scholar

36. M. Norgia, S. Donati and D. D’Alessandro, “Interferometric measurements of displacement on a diffusing target by a speckle tracking technique,” IEEE J. Quantum Electron. 37(6), 800–806 (2001).IEJQA70018-9197 http://dx.doi.org/10.1109/3.922778 Google Scholar

37. S. K. Ozdemir et al., “Compact optical instrument for surface classification using self-mixing interference in a laser diode,” Opt. Eng. 40(1), 38–43 (2001).OPEGAR0091-3286 http://dx.doi.org/10.1117/1.1331268 Google Scholar

38. S. K. Ozdemir et al., “Simultaneous measurement of velocity and length of moving surfaces by a speckle velocimeter with two self-mixing laser diodes,” Appl. Opt. 38(10), 1968–1974 (1999).APOPAI0003-6935 http://dx.doi.org/10.1364/AO.38.001968 Google Scholar

39. C. Zakian, M. Dickinson and T. King, “Particle sizing and flow measurement using self-mixing interferometry with a laser diode,” J. Opt. A Pure Appl. Opt. 7(6), S445 (2005).JOAOF81464-4258 http://dx.doi.org/10.1088/1464-4258/7/6/029 Google Scholar

40. R. B. Miles, W. R. Lempert and J. N. Forkey, “Laser Rayleigh scattering,” Meas. Sci. Technol. 12(5), R33 (2001).MSTCEP0957-0233 http://dx.doi.org/10.1088/0957-0233/12/5/201 Google Scholar

41. M.-C. Michalski, V. Briard and F. Michel, “Optical parameters of milk fat globules for laser light scattering measurements,” Lait 81(6), 787–796 (2001).LAITAG0023-7302 http://dx.doi.org/10.1051/lait:2001105 Google Scholar

42. M. Kerker, The Scattering of Light, and Other Electromagnetic Radiation, pp.27–96, Academic Press, New York (1969). Google Scholar


Alexandra S. Alexandrova is a PhD student and Marie Curie fellow at the University of Liverpool, Physics Department, since October 2012. She received her specialist physics and engineering degree from the National Research Nuclear University MEPhI, Russia, in 2012 in concentrated flux of radiation: matter interaction and laser physics. Her research interests are laser physics, interferometry, and laser physics application to accelerator physics.

Vasilis Tzoganis is a PhD student of the University of Liverpool, Physics Department, and currently at RIKEN, Japan. He received his degree in physics and MSc in electronics from the Physics Department of the University of Ioannina, Greece, in 2009 and 2011, respectively. His research focuses on noninvasive diagnostics for particle accelerators.

Carsten P. Welsch is a full professor and head of the accelerator physics group at the University of Liverpool. He studied physics and economics at UC Berkeley and the University of Frankfurt, where he completed his PhD. His career then took him to RIKEN in Japan, the Max Planck Institute for Nuclear Physics in Germany, and CERN in Switzerland. His research interests include beam diagnostics, accelerator design and optimization, novel acceleration concepts, medical applications, and antimatter studies.

© The Authors. Published by SPIE under a Creative Commons Attribution 3.0 Unported License. Distribution or reproduction of this work in whole or in part requires full attribution of the original publication, including its DOI.
Alexandra S. Alexandrova, Alexandra S. Alexandrova, Vasilis Tzoganis, Vasilis Tzoganis, Carsten P. Welsch, Carsten P. Welsch, } "Laser diode self-mixing interferometry for velocity measurements," Optical Engineering 54(3), 034104 (4 March 2015). https://doi.org/10.1117/1.OE.54.3.034104 . Submission:

Back to Top