Fast Camera Research Articles

Exposure of Sn-Wetted W CPS Targets to Simultaneous NBI Beam and High-Power CW Laser Pulses at the High-Heat Flux OLMAT Facility

First experiments are reported of the simultaneous exposure of a number of Sn-wetted W CPSs and a reference W CPS to 100 ms NBI pulses (divertor steady-state loading conditions) and 2 ms long high-energy laser pulses (divertor ELM like loading conditions) at the High-Heat Flux OLMAT facility. The use of a fast-frame imaging camera allows monitoring the onset of particle ejection from the targets during laser pulses and obtaining the corresponding laser heat fluxes as a measure of the resilience of these targets. Fast camera images are used also to determine ejected particle numbers and to estimate their maximum velocities as laser power is increased in order to compare the influence of W CPS structure on these parameters. In addition, the craters resulting from particle ejection are studied for each target with an optical microscope and a scanning electron microscope. Moreover, in-situ W and Sn particle ejection is followed using visible emission spectroscopy and post-exposure W melting after particle ejection is observed using the energy dispersive X-ray method EDX for all the studied targets. This shows that Sn is unable to protect the underlying W substrate from high-energy laser damage, albeit a subsequent refilling of the formed craters with Sn is visible during NBI-only pulses after laser damage. Thus, it is considered that optimization of surface refilling/replenishment with Sn is needed to improve the W substrate protection. From this work, it is also found that the W CPS reference material has a higher laser heat flux threshold for particle ejection than the Sn-wetted targets. Nevertheless, it is important to take into account that in these experiments with laser pulses, the possible beneficial effects of vapor shielding that can take place during particle irradiation at ELMs or disruptions are not present, thus these experiments represent a worst-case scenario.

Development of MHz X-ray phase contrast imaging at the European XFEL.

We report on recent developments that enable megahertz hard X-ray phase contrast imaging (MHz XPCI) experiments at the Single Particles, Clusters, and Biomolecules and Serial Femtosecond Crystallography (SPB/SFX) instrument of the European XFEL facility (EuXFEL). We describe the technical implementation of the key components, including an MHz fast camera and a modular indirect X-ray microscope system based on fast scintillators coupled through a high-resolution optical microscope, which enable full-field X-ray microscopy with phase contrast of fast and irreversible phenomena. The image quality for MHz XPCI data showed significant improvement compared with a pilot demonstration of the technique using parallel beam illumination, which also allows access to up to 24 keV photon energies at the SPB/SFX instrument of the EuXFEL. With these developments, MHz XPCI was implemented as a new method offered for a broad user community (academic and industrial) and is accessible via standard user proposals. Furthermore, intra-train pulse diagnostics with a high few-micrometre spatial resolution and recording up to 128 images of consecutive pulses in a train at up to 1.1 MHz repetition rate is available upstream of the instrument. Together with the diagnostic camera upstream of the instrument and the MHz XPCI setup at the SPB/SFX instrument, simultaneous two-plane measurements for future beam studies and feedback for machine parameter tuning are now possible.

IR Pulsed Laser Ablation of Carbon Materials in High Vacuum

This work aimed to understand how the energy released by short laser pulses can produce different effects in carbon targets with different allotropic states. The IR pulse laser ablation, operating at 1064 nm wavelength, 3 ns pulse duration, and 100 mJ pulse energy, has been used to irradiate different types of carbon targets in a high vacuum. Graphite, highly oriented pyrolytic graphite, glassy carbon, active carbon, and vegetable carbon have exhibited different mass densities and have been laser irradiated. Time-of-flight (TOF) measurements have permitted the evince of the maximum carbon ion acceleration in the generated plasma (of about 200 eV per charge state) and the maximum yield emission (96 μg/pulse in the case of vegetal carbon) along the direction normal to the irradiated surface. The ion energy analyzer measured the carbon charge states (four) and their energy distributions. Further plasma investigations have been performed using a fast CCD camera image and surface profiles of the generated craters to calculate the angular emission and the ablation yield for each type of target. The effects as a function of the target carbon density and binding energy have been highlighted. Possible applications for the generation of thin films and carbon nanoparticles are discussed.

Two-phase flows downstream, upstream and within Plate Heat Exchangers

Air-water flows were assessed within Plate Heat Exchangers (PHE) with the aid of fast camera imaging. Tests occurred in transparent setups with three chevron angle arrangements (30o/30o, 30o/60o and 60o/60o), representative of low, in-between and high pressure drop channels. Evaluation upstream the PHE inlet happened with Electrical Capacitance Tomography. Three patterns were tested: bubbly, slug and stratified. The effects of flow direction, superficial fluid velocities, two-phase pattern, and chevron angle arrangement on air-water distributions were assessed. The PHE channel outlet is characterized by intense flow recirculation. Bubble entrapment occurs in the core of the recirculation zones. Energy dissipation processes along the PHE channel flow affect the inlet gaseous content, intensifying the mixing process of air and water phases, particularly at flow distribution areas owing to the occurrence of flow acceleration and deceleration. Bubble distribution is wide since the break-up process is rather heterogeneous. Prediction of the maximum bubble diameter was obtained with a modification to Hinze's model. Coalescence can occur with small liquid superficial velocities. At the exit manifold, the recirculation zones affect the two-phase pipe flow. In addition to swirling decay, two-phase flow features and gravitational forces need to be accounted to determine the necessary pipe length to attain stationary process.

Managing Residual Heat Effects in Femtosecond Laser Material Processing by Pulse-on-Demand Operation

Femtosecond laser processing combines highly accurate structuring with low residual heating of materials, low thermal damage, and nonlinear absorption processes, making it suitable for the machining of transparent brittle materials. However, with high average powers and laser pulse repetition rates, residual heating becomes relevant. Here, we present a study of the femtosecond laser pulse-on-demand operation regime, combined with regular scanners, aiming to improve throughput and quality of processing regardless of the scanner’s capabilities. We developed two methods to define the needed pulse-on-demand trigger sequences that compensate for the initial accelerating scanner movements. The effects of pulse-on-demand operation were studied in detail using direct process monitoring with a fast thermal camera and indirect process monitoring with optical and topographical surface imaging of final structures, both showing clear advantages of pulse-on-demand operation in precision, thermal effects, and structure shape control. The ability to compensate for irregular scanner movement is the basis for simplified, cheaper, and faster femtosecond laser processing of brittle and heat-susceptible materials.

Development of a three-dimensional scintillation detector for pencil beam verification in proton therapy patient-specific quality assurance.

Pencil Beam Scanning proton therapy has many advantages from a therapeutic point of view, but raises technical constraints in terms of treatment verification. The treatment relies on a large number of planned pencil beams (PB) (up to thousands), whose delivery is divided in several low-intensity pulses delivered a high frequency (1kHz in this study). The purpose of this study was to develop a three-dimensional quality assurance system allowing to verify all the PBs' characteristics (position, energy, intensity in terms of delivered monitor unit-MU) of patient treatment plans on a pulse-by-pulse or a PB-by-PB basis. A system named SCICOPRO has been developed. It is based on a 10×10×10cm3 scintillator cube and a fast camera, synchronized with beam delivery, recording two views (direct and using a mirror) of the scintillation distribution generated by the pulses. A specific calibration and analysis process allowed to extract the characteristics of all the pulses delivered during the treatment, and consequently of all the PBs. The system uncertainties, defined here as average value + standard deviation, were characterized with a customized irradiation plan at different PB intensities (0.02, 0.1, and 1 MU) and with two patient's treatment plans of three beams each. The system's ability to detect potential treatment delivery problems, such as positioning errors of the treatment table in this work (1° rotations and a 2mm translation), was assessed by calculating the confidence intervals (CI) for the different characteristics and evaluating the proportion of PBs within these intervals. The performances of SCICOPRO were evaluated on a pulse-by-pulse basis. They showed a very good signal-to-noise ratio for all the pulse intensities (between 2×10-3MU and 150×10-3MU) allowing uncertainties smaller than 580µm for the position, 180keV for the energy and 3% for the intensity on patients treatment plans. The position and energy uncertainties were found to be little dependent from the pulse intensities whereas the intensity uncertainty depends on the pulses number and intensity distribution. Finally, treatment plans evaluations showed that 98% of the PBs were within the CIs with a nominal positioning against 83% or less with the table positioning errors, thus proving the ability of SCICOPRO to detect this kind of errors. The high acquisition rate and the very high sensitivity of the system developed in this work allowed to record pulses of intensities as low as 2×10-3MU. SCICOPRO was thus able to measure all the characteristics of the spots of a treatment (position, energy, intensity) in a single measurement, making it possible to verify their compliance with the treatment plan. SCICOPRO thus proved to be a fast and accurate tool that would be useful for patient-specific quality assurance (PSQA) on a pulse-by-pulse or PB-by-PB verification basis.

Realization of a gas puff imaging system on the Wendelstein 7-X stellarator.

A system for studying the spatiotemporal dynamics of fluctuations in the boundary of the W7-X plasma using the "Gas-Puff Imaging" (GPI) technique has been designed, constructed, installed, and operated. This GPI system addresses a number of challenges specific to long-pulse superconducting devices, such as W7-X, including the long distance between the plasma and the vacuum vessel wall, the long distance between the plasma and diagnostic ports, the range of last closed flux surface (LCFS) locations for different magnetic configurations in W7-X, and management of heat loads on the system's plasma-facing components. The system features a pair of "converging-diverging" nozzles for partially collimating the gas puffed locally ≈135mm radially outboard of the plasma boundary, a pop-up turning mirror for viewing the gas puff emission from the side (which also acts as a shutter for the re-entrant vacuum window), and a high-throughput optical system that collects visible emission resulting from the interaction between the puffed gas and the plasma and directs it along a water-cooled re-entrant tube directly onto the 8 × 16 pixel detector array of the fast camera. The DEGAS 2 neutral code was used to simulate the Hα (656nm) and HeI (587nm) line emission expected from well-characterized gas-puffs of H2 and He and excited within typical edge plasma profiles in W7-X, thereby predicting line brightnesses used to reduce the risks associated with system sensitivity and placement of the field of view. Operation of GPI on W7-X shows excellent signal-to-noise ratios (>100 at 2 Mframes/s) over the field of view for minimally perturbing gas puffs. The GPI system provides detailed measurements of the two-dimensional (radial and poloidal) dynamics of plasma fluctuations in the W7-X edge and scrape-off layer and in and around the magnetic islands outside the LCFS that make up the island divertor configuration employed on W7-X.

On-axis full-field swept-source optical coherence tomography for murine retinal imaging

A full-field swept-source optical coherence tomography (FF-SS-OCT) for in vivo murine retinal imaging is demonstrated. The on-axis FF-SS-OCT system was built in a Mach–Zehnder interferometer configuration employing a tunable laser source with an adjustable sweep rate and sweep range in conjunction with a fast 2D-CMOS camera. A large field retinal (coherent) illumination was accomplished using an imaging interface comprised of a short-focal length imaging lens and a contact lens. The magnification between the camera and retina (spatial sampling) was appropriately chosen to record the microscopic structural features of the retina in the image. A pupil stop was employed in the detection path to reject unwanted backscattering from the mouse eye and other sources and limit aberrations distorting the retinal images. In vivo mouse retinal imaging was performed at a sweep rate of 150 Hz to acquire volumes unaffected by the system vibrations, which predominated at lower frequencies. Operating the FF-SS-OCT at this speed yielded an effective axial scan rate of 20 million A-scans/s and a field of view of 820 × 410 µm (24.12° × 12.06°). High-quality retinal B-scans and enface images of the retina were obtained with the SS-FF-OCT, revealing all major retinal layers and vascular plexuses.

Experimental Study to Visualize a Methane Leak of 0.25 mL/min by Direct Absorption Spectroscopy and Mid-Infrared Imaging

Tunable laser spectroscopy (TLS) with infrared (IR) imaging is a powerful tool for gas leak detection. This study focuses on direct absorption spectroscopy (DAS) that utilizes wavelength modulation to extract gas information. A tunable interband cascade laser (ICL) with an optical power of 5 mW is periodically modulated by a sawtooth injection current at 10 Hz across the methane absorption around 3271 nm. A fast and sensitive thermal imaging camera for the mid-infrared range between 3 and 5.7 µm is operated at a frame rate of 470 Hz. Offline processing of image stacks is performed using different algorithms (DAS-F, DAS-f and DAS-2f) based on the Lambert–Beer law and the HITRAN database. These algorithms analyze various features of gas absorption, such as area (F), peak (f) and second derivative (2f) of the absorbance. The methane concentration in ppm*m is determined on a pixel-by-pixel analysis without calibration. Leak localization for methane leak rates as low as 0.25 mL/min is accurately displayed in a single concentration image with pixelwise sensitivities of approximately 1 ppm*m in a laboratory environment. Concentration image sequences represent the spatiotemporal dynamics of a gas plume with high contrast. The DAS-2f concept demonstrates promising characteristics, including accuracy, precision, 1/f noise rejection, simplicity and computational efficiency, expanding the applications of DAS.

PLATO’s signal and noise budget

ESA’s PLATO mission aims the detection and characterization of terrestrial planets around solar-type stars as well as the study of host star properties. The noise-to-signal ratio (NSR) is the main performance parameter of the PLATO instrument, which consists of 24 Normal Cameras and 2 Fast Cameras. In order to justify, verify and breakdown NSR-relevant requirements the software simulator PINE was developed. PINE models the signal pathway from a target star to the digital output of a camera based on physical models and considers the major noise contributors. In this paper, the simulator’s coarse mode is introduced which allows fast performance analyses on instrument level. The added value of PINE is illustrated by exemplary applications.

SiC plasma characterization produced by 3 ns pulse laser at 1010 W/cm2 intensity

A silicon carbide (SiC) plasma has been produced by a Q-switched laser at 1064 nm, 3 ns pulse and 1010 W/cm2 intensity, in high vacuum to investigate the fundamental mechanisms of laser-semiconductor interaction. Although the high hardness and chemical stability restrain the SiC micromachining, it is a valuable semiconductor for devices working under extreme conditions. The plasma has been characterized by using time-of-flight (TOF) measurements of ions detected by an ion collector ring (ICR) and an ion energy analyzer (IEA). Fast CCD camera, quadrupole mass spectroscopy (QMS) and surface profiler have been employed to provide deep insights into the laser-SiC interaction mechanisms. Ion kinetic energies, charge states, plasma expansion velocity, ablation yield per laser pulse, plasma temperature and density have been measured. Applications of the laser-generated plasma are presented and discussed.

Retraction Note: Fast and robust absolute camera pose estimation with known focal length

Retraction Note: Fast and robust absolute camera pose estimation with known focal length

High-speed 4-dimensional scanning transmission electron microscopy using compressive sensing techniques.

Here we show that compressive sensing allows 4-dimensional (4-D) STEM data to be obtained and accurately reconstructed with both high-speed and reduced electron fluence. The methodology needed to achieve these results compared to conventional 4-D approaches requires only that a random subset of probe locations is acquired from the typical regular scanning grid, which immediately generates both higher speed and the lower fluence experimentally. We also consider downsampling of the detector, showing that oversampling is inherent within convergent beam electron diffraction (CBED) patterns and that detector downsampling does not reduce precision but allows faster experimental data acquisition. Analysis of an experimental atomic resolution yttrium silicide dataset shows that it is possible to recover over 25 dB peak signal-to-noise ratio in the recovered phase using 0.3% of the totaldata. Lay abstract: Four-dimensional scanning transmission electron microscopy (4-D STEM) is a powerful technique for characterizing complex nanoscale structures. In this method, a convergent beam electron diffraction pattern (CBED) is acquired at each probe location during the scan of the sample. This means that a 2-dimensional signal is acquired at each 2-D probe location, equating to a 4-D dataset. Despite the recent development of fast direct electron detectors, some capable of 100kHz frame rates, the limiting factor for 4-D STEM is acquisition times in the majority of cases, where cameras will typically operate on the order of 2kHz. This means that a raster scan containing 256^2 probe locations can take on the order of 30s, approximately 100-1000 times longer than a conventional STEM imaging technique using monolithic radial detectors. As a result, 4-D STEM acquisitions can be subject to adverse effects such as drift, beam damage, and sample contamination. Recent advances in computational imaging techniques for STEM have allowed for faster acquisition speeds by way of acquiring only a random subset of probe locations from the field of view. By doing this, the acquisition time is significantly reduced, in some cases by a factor of 10-100 times. The acquired data is then processed to fill-in or inpaint the missing data, taking advantage of the inherently low-complex signals which can be linearly combined to recover the information. In this work, similar methods are demonstrated for the acquisition of 4-D STEM data, where only a random subset of CBED patterns are acquired over the raster scan. We simulate the compressive sensing acquisition method for 4-D STEM and present our findings for a variety of analysis techniques such as ptychography and differential phase contrast. Our results show that acquisition times can be significantly reduced on the order of 100-300 times, therefore improving existing frame rates, as well as further reducing the electron fluence beyond just using a faster camera.

Experimental investigation of electrically enhanced boiling of FC 72 using high-resolution phase-detection diagnostics

This work features the boiling of FC72 on a transparent sapphire substrate. Tests are performed using a thin ITO film (indium tin oxide) coated on the sapphire substrate as the heat source, a fast speed video camera to capture the boiling process using a technique to track the phase (liquid or vapor) in contact with the heating surface, an InfraRed video camera to capture the average temperature on the surface, and a metallic grid to impose an electric field perpendicular to the boiling surface. The maximum average electric field tested in this work is 3.3kV/mm, which led to a 18% increase of the critical heat flux. This study analyses in detail the phase distribution data showing that (1) there is no evidence of microlayer formation, and suggesting that (2) the triple contact line evaporation accounts for approximately 20% of the total heat flux, while (3) the quenching stage accounts for approximately 80%. Finally, phase distribution images are processed to characterize the size of vapor patches, showing that a boiling crisis occurs when the distribution of the vapor patches becomes scale-free, and corroborating the hypothesis that the boiling crisis can be modelled as a bubble interaction instability using a percolation model.

A study of turbulent filaments in the edge plasma of the Wendelstein 7-X stellarator

Filaments are studied by examining fast camera images on the Wendelstein 7-X stellarator. Fast cameras offer a unique perspective, revealing the complex 3D structure of filaments in the entire poloidal cross-section of the plasma. By correlating individual pixels, their location, shape, and movement are analyzed in standard and high-ι configurations. The presence of filaments is not uniform poloidally around. The number of active areas matches the number of magnetic islands in both configurations. Filaments are found to extend to multiple toroidal turns in standard configuration. No time delay is observed between the different toroidal sections. Such behavior is not seen in high-ι configuration. Filaments are observed within and without the edge shear layer, indicated by the direction of their poloidal rotation. Inside the shear layer, their velocity scatters around 1.25 km s−1, accompanied by a lifetime between 80 and 120 µs. Outside, their velocity shows greater absolute values and variance, but still in a few km s−1 range. The similarities and differences between the two configurations are discussed and compared to previous results.

Interpretative 3D MHD modelling of deuterium SPI into a JET H-mode plasma

The pre-thermal quench (pre-TQ) dynamics of a pure deuterium ( D2 ) shattered pellet injection (SPI) into a 3MA / 7MJ JET H-mode plasma is studied via 3D non-linear MHD modelling with the JOREK code. The interpretative modelling captures the overall evolution of the measured density and radiated power. The simulations also identify the importance of the drifts of ablation plasmoids towards the tokamak low field side (LFS) and the impurities in the background plasma in fragment penetration, assimilation, radiative cooling and MHD activity in D2 SPI experiments. It is found that plasmoid drifts lead to an about 70% reduction of the central line-integrated density (compared to a simulation without drifts) in the JET D2 SPI discharge considered. Impurities that pre-exist before the SPI as well as those from possible impurity influxes related to the SPI are shown to dominate the radiation in the considered discharge. With inputs from JOREK simulations, modelling with the Lagrangian particle-based pellet code PELOTON reproduces the deviation of the SPI fragments in the direction of the major radius as observed by the fast camera. This confirms the role of rocket effects and plasmoid drifts in the considered discharge and reinforces the validity of the JOREK modelling. The limited core density rise due to plasmoid drifts and the strong radiative cooling and MHD activity with impurities (depending on their species and concentration) could limit the effectiveness of LFS D2 SPI in runaway electron avoidance and are worth considering in the design of the ITER disruption mitigation system.

Mutual interactions between plasma filaments in a tokamak evidenced by fast imaging and machine learning.

Magnetically confined fusion plasmas are subject to various instabilities that cause turbulent transport of particles and heat across the magnetic field. In the edge plasma region, this transport takes the form of long filaments stretched along the magnetic field lines. Understanding the dynamics of these filaments, referred to as blobs, is crucial for predicting and controlling their impact on reactor performance. To achieve this, highly resolved passive fast camera measurements have been conducted on the COMPASS tokamak. These measurements are analyzed using both conventional tracking methods and a custom-developed machine-learning approach designed to characterize more particularly the mutual interactions between filaments. Our findings demonstrate that up to 18% of blobs exhibit mutual interactions in the investigated area close to the separatrix, at the border between confined and nonconfined plasma. Notably, we present direct observations and radial dependence of blob coalescence and splitting, rapid reversals in the propagation direction of the blob, as well as their dependence on the radial position. The comparison between observations realized with passive and gas puff imaging does not evidence any significant bias due to the use of the latter technique.

Free boundary equilibrium determination of HL-3 discharges

In 2022–2023, HL-3 (previously known as HL-2M) realized the divertor discharge after its first plasma campaign. To analyze equilibrium configurations of HL-3, a free boundary equilibrium code FEEQS, based on the finite element method, has been used in recent experiments. Free boundary equilibrium determinations are carried out for discharges performed in both limiter and divertor configurations. Compared with the experimental diagnostic data from fast camera photos and saturated ion current signals.etc, various time evolution of plasma parameters are well reproduced, such as plasma shape, currents in the poloidal field coils, strike-point displacement. The use of FEEQS successfully supported the operations in HL-3, and physics study as well, which laid a foundation for the application of free boundary equilibrium integration simulation in the next step.

A new automatic sugarcane seed cutting machine based on internet of things technology and RGB color sensor.

Egypt is among the world's largest producers of sugarcane. This crop is of great economic importance in the country, as it serves as a primary source of sugar, a vital strategic material. The pre-cutting planting mode is the most used technique for cultivating sugarcane in Egypt. However, this method is plagued by several issues that adversely affect the quality of the crop. A proposed solution to these problems is the implementation of a sugarcane-seed-cutting device, which incorporates automatic identification technology for optimal efficiency. The aim is to enhance the cutting quality and efficiency of the pre-cutting planting mode of sugarcane. The developed machine consists of a feeding system, a node scanning and detection system, a node cutting system, a sugarcane seed counting and monitoring system, and a control system. The current research aims to study the pulse widths (PW) of three-color channels (R, G, and B) of the RGB color sensors under laboratory conditions. The output PW of red, green, and blue channel values were recorded at three color types for hand-colored nodes [black, red, and blue], three speeds of the feeding system [7.5 m/min, 5 m/min, and 4.3 m/min], three installing heights of the RGB color sensors [2.0 cm, 3.0 cm, and 4.0 cm], and three widths of the colored line [10.0 mm, 7.0 mm, and 3.0 mm]. The laboratory test results s to identify hand-colored sugarcane nodes showed that the recognition rate ranged from 95% to 100% and the average scanning time ranged from 1.0 s to 1.75 s. The capacity of the developed machine ranged up to 1200 seeds per hour. The highest performance of the developed machine was 100% when using hand-colored sugarcane stalks with a 10 mm blue color line and installing the RGB color sensor at 2.0 cm in height, as well as increasing the speed of the feeding system to 7.5 m/min. The use of IoT and RGB color sensors has made it possible to get analytical indicators like those achieved with other automatic systems for cutting sugar cane seeds without requiring the use of computers or expensive, fast industrial cameras for image processing.

Performance characteristics of the airlift pump under vertical solid–water–gas flow conditions for conveying centimetric-sized coal particles

In this study, the installation of an airlift pump with inner diameter of 102 mm and length of 5.64 m was utilized to consider the conveying process of non-spherical coal particles with density of 1340 kg/m3 and graining 25–44.5 mm. The test results revealed that the magnitude of increase in the solid transport rate due to the changes in the three tested parameters between compressed air velocity, submergence ratio, and feeding coal possibility was not the same, which are stand in range of 20%, 75%, and 40%, respectively. Hence, creating the optimal airlift pump performance is highly dependent on submergence ratio. More importantly, we measured the solid volume fraction using the method of one-way valves in order to minimize the disadvantages of conventional devices, such as fast speed camera and conductivity ring sensor. The results confirmed that the volume fraction of the solid phase in the transfer process was always less than 12%. To validate present experimental data, the existing empirical correlations together with the theoretical equations related to the multiphase flow was used. The overall agreement between the theory and experimental solid delivery results was particularly good instead of the first stage of conveying process. This drawback can be corrected by omitting the role of friction and shear stress at low air income velocity. It was also found that the model developed by Kalenik failed to predict the performance of our airlift operation in terms of the mass flow rate of the coal particles.