Electrical Tomography: A Review of Configurations, and Application to Fibre Flow Suspensions Characterisation

applied sciences Review Electrical Tomography: A Review of Configurations, and Application to Fibre Flow Suspensions Characterisation Pedro Faia 1,* , Rui Silva 2, Maria G. Rasteiro 2and Fernando Garcia 2 1Department of Electric and Computers Engineering, Faculty of Sciences and Technology, CEMMPRE, University of Coimbra, 3030-290 Coimbra, Portugal 2 Department of Chemical Engineering, Faculty of Sciences and Technology, CIEPQPF, University of Coimbra, 3030-790 Coimbra, Portugal; [email protected] (R.S.); [email protected] (M.G.R.); fgar[email protected] (F.G.) *Correspondence: [email protected]; Tel.: +351-239-796-200 Received: 6 March 2020; Accepted: 19 March 2020; Published: 30 March 2020   Featured Application: The depicted tomographic technique, Electrical Tomography, is a non-invasive and remotely operated imaging system that makes it possible to obtain valuable information about multiphase system flow behaviours: The 2D/3D images obtained are of extreme importance not only in the development of prediction models, but also in their real-time monitoring and control. It has high applicability in several industrial environments, such as chemistry, paper, petroleum, water and wastes industries, just to mention some. Abstract: Understanding the behaviour of suspension flows continues to be a subject of great interest considering its industrial relevance, regardless of the long time and effort dedicated to it by the scientific and industrial communities. Information about several flow characteristics, such as flow regimen, relative velocity between phases, and spatial distribution of the phases, are essential for the development of exact models for description of processes involving pulp suspension. Among the diverse non-invasive techniques for flow characterisation that have been reported in the literature for obtaining experimental data about suspension flow in different processes, Electrical Tomography is one of the most interesting, since it presents perhaps the best compromise among cost, portability, and, above all, safety of handling (indeed there is no need to use radiation, which requires special care when using it). In this paper, a brief review and comparison between existing technologies for pulp suspension flow monitoring will be presented, together with their strengths and weaknesses. Emphasis is given to Electrical Tomography, because it offers the above-mentioned compromise and thus was the strategy adopted by the authors to characterise different flow processes (solid–liquid, liquid–liquid, fibres, etc.). The produced portable EIT system is described, and examples of results of its use for pulp suspension flow characterisation are reported and discussed. Keywords: fibre suspensions flow; electrical tomography; complex flows; disperse phase distribution; velocity evaluation 1. Introduction For the obtention of tomographic information, different measurement methods can be followed with respect to the sensing procedure: intrusive, non-intrusive, invasive, and non-invasive. The choice of an approach is the result of a wise choice made by the users. By ‘invasive’ is meant that the sensor surface is parallel to the inner zone of the wall of the channel under study, while the term ‘intrusive’ is used when the sensor itself interacts with the flow being visualised. In opposition non-invasive refers to cases in which the sensor is located on the outside of the wall, without crossing it, while Appl. Sci. 2020,10, 2355; doi:10.3390/app10072355 www.mdpi.com/journal/applsci Appl. Sci. 2020,10, 2355 2 of 30 non-intrusive indicates cases in which the sensor placement forces the pipe wall to be crossed, but without disturbing the flow under inspection [ 1 ]. In recent years, non-invasive procedures have shown a high degree of potential growth in many areas, such as in the chemical industry. When compared with non-intrusive procedures, non-invasive approaches bring benefits to the above-mentioned industry sectors, such as reduction of hazards due to operation of corrosive or radioactive materials, increase in operator safety, etc. Consequently, the non-invasive approach is the imaging strategy most commonly used in process plants. For hard-field tomography, the case must be considered in which there the measured parameter distribution and the sensitivity of the imaged media in the entire 3D domain are independent of one another [ 2 ]; optical, ultrasonic and X-ray are examples of this type of tomography. On the other hand, soft-field tomography refers to cases in which a dependency is found, over the entire volume, between the medium sensitivity and the distribution of the measures of the desired variable [ 3 ]. Electrical tomography is the most relevant member of this class, but it presents an important effort regarding the resolution of the inverse problem on the domain under observation. 2. Pulp Suspensions Imaging Overview Pulp suspensions are characterised as being opaque and intricate mixtures that exhibit a quite heterogeneous distribution of fibres under flow, presenting areas in which flocs are present (the term flocs refers to a high concentration of fibres). In addition, along the entire pulp suspension and above certain concentrations, as a function of the observed shear rates, a continuous fibre network forms. For the purposes of process control nowadays, knowledge regarding flow characteristics is essential, once the suspensions have been subjected to high shear forces. The numbers regarding turnover in the pulp and paper industry confirm the importance of this knowledge, as in some countries they exceed the return from other industrial areas, as is the case in Sweden and Canada. Understanding and controlling the manufacturing process enables the development of newer products with a high innovation content. Additionally, obtaining precise values for the velocity profiles associated with a certain shear rate in a pulp suspension is also of high importance for optimising the flow and, consequently, the energy consumption; for the later, however, there is scarce information reported in the literature. Diverse flow imaging methods, along with particle following velocimetry approaches, have been described in the literature; however, most of these need to be experimentally trained, with data acquisition procedures thus usually lasting quite a long time, and they can typically only be used in non-opaque suspensions. Consequently, it is not possible to use them in pulp suspension flows not belonging to the dilute range. For data collection in pulp suspensions, among the more promising techniques can be found Laser Doppler Anemometry (LDA) (Optomet GmbH, Darmstadt, Germany) [ 4 ], Magnetic Resonance Imaging (MRI) (National Institute of Biomedical Imaging and Bioengineering, Bethesda, MD, USA) [ 5 ], Ultrasonic Velocity Profiling (UVP) (MET-FLOW S.A., Lausanne, Switzerland) [ 6 ], Optical Coherence Tomography (OCT) (PHOENIX TECHNOLOGY GROUP,Pleasanton, CA,USA)[ 7 ], ParticleImageVelocimetry(PIV)(DANTECDYNAMICS,Skovlunde, Denmark) [ 8 ], X-ray tomography (Department of Physics P.O., University of Jyväskylä, Jyväskylän yliopisto, Finland) [ 9 ], Gamma ray densitometry [ 10 ], and Electrical Tomography (ET) (PASI SRL, Torino, Italy) [11]. Magnetic Resonance Imaging is a tomographic technology with a high spatial resolution that can supply information regarding the behaviour of the protons present in a system, which are typically contained in the 1H nuclei of water. MRI can be used, among other applications, to image the motion of water, e.g., expressed in the form of velocity profiles. The major limitations of MRI are the types of particles and the size of the system that can be studied, together with the size, weight and cost of the instrumentation. Only particles containing MR-sensitive nuclei, such as 1H, can be detected [ 12 , 13 ]. The maximum diameter of the system is bounded by the inner diameter of the coil in the magnet. Thus, MRI experiments are typically limited to laboratory-scale fluidised beds or pipes with small Appl. Sci. 2020,10, 2355 3 of 30 diameters [ 14 ]. In summary, the strength of MRI is in studying centimetre-scale systems at spatial resolutions of ~100 µ m. Time-averaged velocity profiles of diluted pulp suspension have been obtained using MRI with fibre concentrations of up to 0.86% (w/w) [ 15 – 20 ]. However, when compared with LDA and UPV, MRI presents some disadvantages, such as longer observation times and equipment cost. Arola et al. [ 20 ] were able, to gather information about the suspension velocity profiles of 0.5% fibres diluted in water (w/w) using MRI in a matter of milliseconds. Also using MRI, under steady flow conditions, and for pulp suspensions with concentrations above 3% (w/w), Seymour et al. [ 21 ] were able to measure time-averaged velocity profiles. Despite these lab trials that have taken place around the world, wide MRI usage in industry is not expected in the coming years. Laser Doppler Anemometry (LDA) is an optical technique used to investigate velocity and turbulence in gas, liquid, and mixed fluids, among other applications [ 22 ]. The basic idea underlying LDA is to measure the velocity of tiny particles transported by the flow. If these particles are small enough, their velocity is assumed to be that of the liquid stream, and LDA thus provides a measure of the local instantaneous velocity, the mean velocity, and the turbulent quantities. Laser anemometers offer unique advantages in comparison with other fluid flow instrumentation [ 22 ]: (1) Non-contact optical measurement—LDA probes the flow with focused laser beams and can determine the velocity without disturbing the flow; (2) Well-defined directional response—the quantity measured by LDA is the projection of the velocity vector in the measuring direction defined by the optical system; (3) High spatial and temporal resolution—the optics of the laser anemometer are able to define a very small measuring volume, thus providing good spatial resolution and allowing for local measurement of velocity; (4) Multi-component and multi-directional measurements—combinations of laser anemometer systems with component separation based on colour, polarisation or frequency shift allow one-, two-, or three-component LDA systems to be put together based on common optical modules. So indeed, LDA has a lot of advantages, but some compromise has to be made when selecting and setting up laser anemometer systems; indeed, among the disadvantages of LDA are the high cost of the equipment, the need for a transparent flow through which the light beams can pass (usually corresponding to a low concentration of the disperse phase), and the fact that they do not give continuous velocity signals. Nevertheless, the special properties of the gas lasers (such as high energy, spatial and temporal coherence, and stability) make this method applicable for solving a large number of problems. LDA has been used to examine low concentrated pulp suspensions containing air bulbs by Ek et al. [ 23 ]; carry out a study of 0.5% fibre suspensions in water (w/w) by Kerekes et al. [ 24 ]; study pulp suspension flow at concentrations between 1.2 and 12 g per litre of water, using transparent model systems, based on the refractive index matching of the liquid and the fibres by Steen [ 25 , 26 ]; and, using the same model system presented by Steen, performing measurements in a stirred tank with fibre concentrations in the range 3–20% (w/w) by both Andersson et al. [27] and by Pettersson and Rasmussen [28,29]. Since Acoustic Doppler devices achieved easy handling and low cost, their use has experienced a high rate of growth in recent years. Doppler effect-based echography equipment allows the fabrication of systems that are able to measure velocity profiles nearly instantaneously [ 30 ]. The Ultrasonic Velocity Profile (UVP) technique was initially limited to the inspection of non-transparent fluids [ 31 ], and to measurements through pipe walls in small-scale systems [ 32 ]. In the UVP technique, a probe, acting firstly as an emitter, produces an ultrasound signal that travels along a certain direction, usually the suspension direction, and later, acting as a receiver, senses the reflexions originated by particles present in the suspension; the time, as well as the Doppler frequency shift between emission and reflections, is measured. Using the value of sound velocity in the fluid, the distance and velocity of the particles originating the echoes (the reflexions), along the suspension direction, are calculated. UVP was originally established for 1D measurements; however, promising 2D velocity profiles have recently been published [ 30 ]. The UVP technique has been used extensively in recent years to obtain instantaneous velocity profiles, but very few studies on pulp suspensions can be found in the literature. Hirsimäki [ 33 ] used an early ultrasonic profiling technique to obtain radial velocity profiles in pulp suspensions with concentrations up to 1% (w/w). Karema and his co-workers [ 34 , 35 ] Appl. Sci. 2020,10, 2355 4 of 30 characterised velocity fluctuations and studied paper formation by fluidisation and reflocculation of wood pulp suspensions with concentrations up to 1% (w/w) using UVP. Wiklund et al. [ 31 , 36 ] obtained instantaneous velocity profiles in a steady laminar flow of 0.5–3% (w/w) and later up to 7.8% (w/w) cellulose pulp suspensions using UVP. Claesson and her co-workers [ 37 ] used UVP to characterise pulp fibre suspensions through a sudden expansion at two different concentrations, 1.8% and 2.8% (w/w), and at diverse velocities, from 1 up to 2.2 m/s: They were able to instantaneously obtain velocity profiles as well as axial real mean square velocity fluctuations in the evaluated semi-concentrated flows just after the expansion. To enable improved flow information for higher concentrations, Kotz é et al. [ 38 ] tested a pulsed velocimetry-based system to measure velocity profiles of cellulose flows for two different concentrations, 6.9% and 8% (w/w), in three different pipe diameters (16, 22.5 and 52.8 mm). Also using a pulsed approach, Hanjiang and colleagues [ 39 ] obtained velocity profiles of fibre flows in a rectangular channel, and investigated the influence of fibre concentration and Reynolds number on the shape of the measured profiles; they used natural wood fibres (with an average length and diameter of 2.3 and 35 µ m, respectively), and the tests were conducted with varying average flow velocities of between 0.16 and 7.0 m/s, and different fibres concentrations between 0.05% and 1.0% (w/w). Optical Coherence Tomography (OCT) is a recent technique that allows the internal structure and the motion behaviour of opaque suspensions to be measured simultaneously and with high spatial resolution. Depending on the technology used in the OCT construction, the axial scan rates can go up to hundreds of kHz [ 7 ]. OCT operating modus is quite similar to ultrasound imaging, with the main difference being its operating frequency; while OCT operates in the light frequency band, ultrasound uses the sound band. The depth of the obtained images depends to a large extent on the optical properties of the domain; consequently, its distance from the wall can vary between micrometres and some millimetres [ 40 ]. This technique was used by Kataja et al. [ 41 ] to characterise the boundary layer of micro-fibrillated cellulose flowing in a straight pipe, for consistencies in the range 0.4–1.6% (w/w). However, in the case of ultrasound systems, which usually operate in the MHz frequency band, the waves are readily transmitted to the domain, making it easy to obtain images of significant depth: nevertheless, at those frequencies they can also suffer from high attenuation, limiting then the used frequencies. Furthermore, it is much more difficult to focus sound waves than light waves [42]. Particle Image Velocimetry (PIV) is a laser measuring technique, as is LDA, which measures the velocity by means of tracking tracers in the suspension. However, PIV measures the velocity by cross-correlation algorithms, allowing the velocity 3D-components or the velocity gradient tensor in the measuring plane to be obtained [ 43 ]. PIV has the advantage of not being a single point measurement technique, and consequently is attracting increasing interest in the evaluation of complex and turbulent flows. The main disadvantage of this technique is the limitation it presents when more opaque domains are the object of visualisation. Sorokin and his co-researchers used PIV to estimate sets of 2D pulp flow velocity vectors [ 8 ]. With recent improvements in digital imaging, PIV also makes it possible to obtain information about the orientation of the fibres, as in Fan et al. [ 44 ]. In the experiments, they used synthetic fibres with diverse aspect ratios (18.8, 25.0 and 37.5) in two different concentrations, 0.1% and 0.2%. They obtained important information about the influence of the impeller speed, fibre aspect ratio and concentration on the field velocity and orientation. Fock and colleagues [ 45 ] used PIV combined with LDA to further understand the mechanisms of plug flow and near-wall behaviour in fibre flow through a pipe at concentrations of up to 4.7% (w/w) with a flow rate of 0.003 m 3 /s; the obtained results showed that the pulp suspension flow in the plug regime can be very inhomogeneous in the region near the pipe wall, particularly at low flow velocities, becoming more structured with an increase in flow velocity. Among the difficult field techniques, X-ray computed tomography is the safest one; it has excellent spatial resolution and image reconstruction is easy, but due to the source and construction of detection components, its operation is slow [ 9 ]. X-ray usage for fibre flow characterisation has rarely been reported in the literature; however, flash X-ray radiography has been used to characterise pulp flows similar to those that can be found in the paper industry. Flash X-ray is an X-ray radiography method Appl. Sci. 2020,10, 2355 5 of 30 that uses an intense burst of radiation that is produced for a short time period, making it possible to capture images of high-speed events that are hidden by, for instance, dust, smoke or light [ 9 ]. Farrington used Flash X-ray [ 46 ] in order to image tungsten fibres with a median diameter of 25 µ m (acting as tracers) in a wood fibre suspension with a concentration of 2% (w/w). He was able to asses fibre mass distribution, confirming that this technique could be used for pulp flow characterisation. Later, Heidel and co-workers [ 47 , 48 ] used Flash X-ray to characterise and quantify gas behaviour in several suspensions containing fibres; they used it to visualise air flow in diverse systems, particularly in one composed of air/water/wood fibres, with fibre consistencies of 0.5%, 1% and 1.5% (w/w), while varying the average injection velocity of the air between 0.1 and 12.5 cm/s and maintaining the remaining mixture injection velocity at a constant rate in the pipe. Gamma ray has also been tried for inspection of pulp suspensions, particularly for acquiring density/void fraction information [ 49 ]. Gamma ray densitometry [ 10 ] uses Gamma radiation, which is an electromagnetic radiation similar to X-rays, but which is produced by the natural decay of radionuclides. It is reliable for the determination of densities, it is inexpensive, and it is portable; however, it does not allow local information to be obtained (the collimated beam of the Gamma ray produces a line averaged value and local information in points of the chord cannot be obtained) [ 50 ]. Also, Gamma rays’ intensities in the domain are quite dependent on the thickness of the pipe wall, and their accuracy depends on the radiation reaching the detector, which is affected by the material and thickness of the pipe wall [ 50 ]. Xie et al. [ 51 ] combined Flash X-ray tomography with Gamma ray densitometry to study flow structures, gas holdup, and the geometry of the gas bubbles in a gas–pulp–liquid system; in their work, low-fibre suspension consistencies were used—up to 1.5% (w/w)—and superficial flow velocities of gas and liquid/pulp suspensions in the ranges 0–26 and 21–51 cm/s, respectively. The study of larger-scale systems can be addressed by Electrical Impedance Tomography (EIT), which makes it possible to obtain good images of the distribution of gas–liquid–solid mixtures [ 52 ]. Low cost, portability and ease of upscaling are some of the advantages of EIT. However, its spatial resolution is lower than that of other tomographic techniques. This is mainly due to the number of matchless conductivity measurements that can be obtained, which is strongly tied to the current injection-voltage acquisition system protocol and to the number of measuring electrodes [ 53 ]. The reconstruction procedure is poorly conditioned and non-linear [ 53 ]. Consequently, the amount of fibres or particles present in the suspension cannot be precisely calculated, with the obtained tomographic images representing the distribution of variation in conductance/impedance, rather than absolute values. To maximise the EIT conductivity distribution images resulting from the reconstruction process, the conductivity difference between the two phases to be distinguished should be large, with larger conductivity differences providing better distinguishability and measurement accuracy. In the pulp and paper industry, EIT has recently been used in applications including the measurement of mixing quality, flow uniformity, and velocity, just to mention a few [11,54–56]. 3. Electrical Tomography Principles and Applications 3.1. Electrical Tomography Due to its low cost, portability and handling safety (because no harmful radiation is used), Electrical Impedance Tomography exhibits strong potential for use in several academic/industrial areas. It has a simple setup, and offers quick and easy operation, making it robust enough to be used in many industrial situations. As has already been mentioned, its main drawback is its low spatial resolution, which is usually confined to between 6% and 20% of the pipe radius [ 53 ]. However, with the development of more advanced algorithms for inverse problem resolution, there have been significant improvements with respect to increasing the achievable spatial resolution. With Electrical Tomography, both quantitative and qualitative data regarding multiphase systems can be obtained for modelling purposes. Appl. Sci. 2020,10, 2355 6 of 30 Non-invasive collection of data by electrical tomography enables the reconstruction of cross-section images, providing material distribution profiles in a suspension flowing through a pipe, or information about possible transient events taking place. The resulting conclusions can then be used in processes redesign or control. In some cases of the real-time imaging of industrial processes, electrical tomography is the most appealing technique, due to its simplicity and low cost. Electrical tomography can be divided in two procedures: Electrical Impedance Tomography, EIT (note that Electrical Resistance Tomography, ERT, can be seen has a particular case of EIT), and Electrical Capacitance Tomography, ECT. They both allow the generation of images that reflect the variation of the conductivity/resistivity or permittivity along the under-observation domain; due to its simplicity of implementation, ERT, which is ideal for pure resistive domains, is the most used [57]. When systems are heterogeneous, and composed of materials possessing diverse electrical characteristics, EIT and ECT can be used in the observation of changing processes; pneumatic and hydraulic transportation, mixing, cyclones, and fluidised beds are only some examples. The number of images generated by time interval is variable: it depends mainly on the electronics and structure of the acquisition system, on the measurement protocol used to gather the data, and on the reconstruction algorithm used for image building. Typically, increasing the data gathering speed would make it possible to obtain data about faster process, but consequently with a higher noise level and reduced image quality of the distribution profiles. Electrical tomography can be used on-line or off-line. For on-line measurements, poor spatial resolution images are usually obtained, consequence of the need to use either fast acquisition procedures or data processing algorithms that only require a short time. Conversely, for off-line measurements, the acquisition times can be customised to adapt to the changes taking place in the flow process, while the processing algorithms can be optimised to the highest spatial resolution possible, because the processing time can be longer. In addition to the reconstruction algorithm used, other features determine the quality of the obtained images, mainly physical ones, such as the number of electrodes, electrode diameter and material, spacing between them, etc. Because in some applications the aim is to obtain information about average quantities, such as void fraction or mean velocity, the developed model’s validation requires images with high spatial resolution. In the following sections, the authors start by addressing the basic principles and operation of Electrical Impedance Tomography, followed by a short review of the use of Electrical Tomography for industrial/process monitoring purposes, with an emphasis on the pulp and paper area. Finally, a succinct description of the authors’ approach to developing high spatial resolution portable EIT equipment is presented; the authors’ EIT system makes it possible not only to inspect media with high conductivity, as are typically found in industry, as well as using a reconstruction procedure that makes it possible to obtain sharper impedance distribution images. The results obtained with this new system using fibre suspensions on a laboratory level pilot rig will be presented. Finally, the incoming challenges for Electrical Tomography are addressed. 3.2. Fundamentals of Electrical Tomography Maxwell [ 58 ] was the first to present a theoretical model connecting the permittivity/conductivity of a two-phase mixture to the volume fraction of one the compounds present in the mixture dispersed in the second. For his calculations, he used small spheres of a certain material possessing equal but small size when compared with the distance between them; additionally, he presumed that in a certain domain, the spheres were distributed in a uniform way along the continuous phase of a second material, and that consequently their presence disturbs a homogeneous electrical field passing through the domain. Extrapolating, by measuring quantities such as conductance/capacitance, between a pair of electrodes placed around the pipe wall, the particle distribution inside that pipe can be obtained by means of the distribution of the electrical properties and through the usage of an adequate mathematic procedure: this mathematic algorithm links the measurements with the particle distribution. The relation between the homogeneity of a multiphase flow and this idea has been Appl. Sci. 2020,10, 2355 7 of 30 used for some time to observe different systems, such as liquid–solid and gas–liquid. For instance, to quantify localised amounts of solids in a gas flow, small capacitive probes have been used [ 59 , 60 ]. In the gas and oil industry, similar probes have been used with analogous objectives [ 61 ]. Typically, in electrical tomography, the excitation sources, either a voltage or a current source, operate at low frequencies (below 5 MHz). Consequently, the systems in which these types of sources are used are described by the governing equations of the electrostatic field. When the flux (or current) lines encounter an interface with different permittivity or conductivity, it gets deflected. In order to allow the mining of the maximum amount of information about the process, the electrodes are usually mounted over the periphery of the pipe at intervals of equal distance. In the case of capacitive systems, the measuring electrodes are installed outside of the process vessel or pipe which is made of a dielectric material (i.e., they are mounted in a non-invasive way). Furthermore, the electrode area must be high, making it possible to detect a large range of capacitance changes. For system measuring Resistivity changes, the sensing electrodes are quite small when compared to those used in capacitive systems. They will typically be in contact with the flow under inspection, mounted level along the inner side of the pipe wall. 3.2.1. Capacitance Tomography The purpose of ECT it to obtain a reconstruction of the dielectric characteristics of a system, using measures of capacitance taken between all existing electrode pairs. In Figure 1, a schematic of a cross-section representation of an ECT system possessing eight electrodes can be seen. If the pipe where the suspension is flowing is made of a conductive material, then the electrodes are internally mounted. If, conversely, the pipe is made of an insulating material, then the electrodes will be located externally. External electrodes are easily conceived, assembled and maintained; additionally, once they have been directly exposed to extreme turbulence, pressures and temperatures, their electrical behaviour remains unchanged for a longer time (contamination by the materials flowing in the tube is not encountered, either). Their greatest disadvantage is that they possess non-linear characteristics; however, the use of adequate correction factors allows them to display almost linear characteristics. The design of internal electrodes is harder and more intricate, because they are exposed as described above to extreme conditions, and may also suffer from corrosion; however, the capacitance changes observed with internal electrodes can be assumed to be directly proportional to the permittivity changes occurring inside the pipe. In ECT, all electrodes are stimulated, one at a time, and all pairs of capacitances between the stimulated electrode and the remaining ones are measured; for an N electrodes system, this produces N(N-1)/2 independent measurements, once capacitance C i,j =C j,i and C i,i , i.e., the self-capacitance, are all ignored. In summary, the measuring protocol operates in a similar way to the source-detector movement in computerised tomography in medical imaging; in ECT, the electrical field rotates around the pipe cross-section in steps, with a step angle of 360 ◦ /N. The number Nof electrodes is chosen based on a balance between the desired spatial resolution and acquisition rate of the images (acquisition rates around 100 frames per second are frequent). Maxwell equations make it possible to establish a relationship between the measured capacitances and the actual spatial distribution of the permittivity. In ECT, only one electrode is stimulated at one time, while the remaining are at a virtual earth potential. Consequently, over tall electrode surfaces, the total electric flux is equal to zero. The capacitance measures can be easily influenced by external parasite capacitances (usually greater than the measured ones), such as earth capacitance, because the inter-electrode capacitance is quite small. To prevent this perturbation, electrode shielding is implemented. With respect to the excitation signal frequency, a value of around 1 MHz is common; for this order of magnitude, the signal wavelength is of some hundreds of meters, greatly exceeding the sensor size by many orders of magnitude. In summary, electrostatic field theory can be described by the electrical potential distribution inside the medium under inspection. Appl. Sci. 2020,10, 2355 8 of 30 Appl. Sci. 2020, 10, x FOR PEER REVIEW 7 of 31 conductivity, it gets deflected. In order to allow the mining of the maximum amount of information about the process, the electrodes are usually mounted over the periphery of the pipe at intervals of equal distance. In the case of capacitive systems, the measuring electrodes are installed outside of the process vessel or pipe which is made of a dielectric material (i.e., they are mounted in a non-invasive way). Furthermore, the electrode area has to be high, making it possible to detect a large range of capacitance changes. For system measuring Resistivity changes, the sensing electrodes are quite small when compared to those used in capacitive systems. They will typically be in contact with the flow under inspection, mounted level along the inner side of the pipe wall. 3.2.1. Capacitance Tomography The purpose of ECT it to obtain a reconstruction of the dielectric characteristics of a system, using measures of capacitance taken between all existing electrode pairs. In Figure 1, a schematic of a crosssection representation of an ECT system possessing eight electrodes can be seen. If the pipe where the suspension is flowing is made of a conductive material, then the electrodes are internally mounted. If, conversely, the pipe is made of an insulating material, then the electrodes will be located externally. External electrodes are easily conceived, assembled and maintained; additionally, once they have been directly exposed to extreme turbulence, pressures and temperatures, their electrical behaviour remains unchanged for a longer time (contamination by the materials flowing in the tube is not encountered, either). Their greatest disadvantage is that they possess non-linear characteristics; however, the use of adequate correction factors allows them to display almost linear characteristics. The design of internal electrodes is harder and more intricate, because they are exposed as described above to extreme conditions, and may also suffer from corrosion; however, the capacitance changes observed with internal electrodes can be assumed to be directly proportional to the permittivity changes occurring inside the pipe. In ECT, all electrodes are stimulated, one at a time, and all pairs of capacitances between the stimulated electrode and the remaining ones are measured; for an N electrodes system, this produces N(N1)/2 independent measurements, once capacitance Ci,j = Cj,i and Ci,i, i.e., the self-capacitance, are all ignored. In summary, the measuring protocol operates in a similar way to the source-detector movement in computerised tomography in medical imaging; in ECT, the electrical field rotates around the pipe crosssection in steps, with a step angle of 360°/N. The number N of electrodes is chosen based on a balance between the desired spatial resolution and acquisition rate of the images (acquisition rates around 100 frames per second are frequent). Maxwell equations make it possible to establish a relationship between the measured capacitances and the actual spatial distribution of the permittivity. In ECT, only one electrode is stimulated at one time, while the remaining are at a virtual earth potential. Consequently, over tall electrode surfaces, the total electric flux is equal to zero. The capacitance measures can be easily influenced by external parasite capacitances (usually greater than the measured ones), such as earth capacitance, because the inter-electrode capacitance is quite small. To prevent this perturbation, electrode shielding is implemented. With respect to the excitation signal frequency, a value of around 1 MHz is common; for this order of magnitude, the signal wavelength is of some hundreds of meters, greatly exceeding the sensor size by many orders of magnitude. In summary, electrostatic field theory can be described by the electrical potential distribution inside the medium under inspection. Figure 1. Schematic representation of the measurement principle of an ECT system (reproduced with permission from Hansen et al. [62]). ERT sensing arrays are less complex than those required by ECT, and the difficulty in their conception arises when the domains present high conductivity [ 63 , 64 ]. Consequently, it can be stated that it is quite adequate for the inspection of non-conductive mixtures containing diverse permittivity’s. ECT produces low-resolution images, but, when dealing with disperse system flows, offers good accuracy with respect to volume fraction estimation [ 65 ]. If in ECT a dual plane system is used, then cross correlating the cross-section average time series enables the evaluation of a particle’s velocity [ 65 ]. 3.2.2. Resistivity/Impedance Tomography The impedance variation distribution of a mixture is the focus of EIT. In EIT, the stimulus used is a sinusoidal electrical signal (current or voltage) with a certain frequency and amplitude. Then, the real and imaginary parts of the impedances are calculated using phase sensitivity measurements [ 66 , 67 ]. The resistive component is obtained by an in-phase measurement with the stimulus signal, while the capacitive component is found by a quadrature of phase measurement with the excitation signal [ 53 ]. The invention credit for EIT is given to John G. Webster, who reported it in a 1978 publication [ 68 ]. Still, the first known experimental application of EIT took place in 1984, for which the authors were Barber and Brown [ 69 ]. For investigating processes where the continuous phase is electrically conductive, either ERT or EIT can be used. Specifically, when the inspected medium is purely resistive, then ERT should be used. Similarly to ECT, in ERT/EIT an electrical current is injected through a pair of electrodes that are placed along the boundary of the domain under study. Consequently, the resulting electrical field distribution in the domain will be conditioned by the material distribution within it [ 70 ]. The potential voltage differences between all pairs of electrodes placed around the domain perimeter, apart from the pair used for exciting the domain, are measured, and used to nourish a non-linear algorithm. This algorithm solves the inverse problem described above and permits the calculation of the unknown conductivity/resistivity distribution along the domain. However, the procedure only ends when all electrodes are used for stimulating the domain, and thus a full cycle has as many projections as the number of electrodes comprising the sensing system (see Figure 2) [ 71 ]. Therefore, in summary, for estimating the material distribution in the domain, an electrical field distribution characterisation is used. In ERT, which is easier to implement [ 67 , 72 ], the aim is to obtain a conductivity/resistivity distribution along the domain. The working principle used in ERT to obtain the data related to a full set of projections, is identical to that described for EIT (requiring the use of all electrodes to introduce the electrical excitation signal in the domain, and to measure all the potential voltage differences between the remaining ones). The single difference is that in ERT, the electrical signal used as stimulus is a continuous current. Therefore, for EIT, the data acquisition cycle only ends when a full rotation of the applied electric field excitation is complete. Consequently, for both EIT and ERT, if the sensing system is composed of Nelectrodes, the number of independent differential voltage measurements is N(N − 3)/2. Typical values for the frequencies used in EIT systems lie in the range 20–150 kHz, so quasi-static conditions can be assumed when a full rotation data is Appl. Sci. 2020,10, 2355 9 of 30 acquired. In recent years, efforts have been directed towards the development of acquisition systems and to the conception of sensors (materials used in their fabrication, dimensions and location of the sensors, etc.). The use of a discretised electrode system is restricted to based flows with constant conductance [ 73 ]. When large bubbles are present in the fluid, those bubbles can cause electrodes to lose contact with the fluid, and consequently, the obtained images of the conductance distribution will be inconsistent. When faced with this type of suspension, one possible solution would be the use of continuous electrode rings instead of a discrete electrode sensing system [ 73 ]. An additional limitation of using discrete electrodes mounted in the pipe wall, the surface of which is parallel to the inner pipe wall (as required by resistive/impedance tomography), is their surface modification over time, particularly when used for evaluation processes that operate under extreme conditions. Temperature is one example of such conditions: temperature changes can dramatically cause fluid conductivity to alter, and consequently fluid temperature compensation mechanisms have to be applied to “correct” current/voltage measurements. Appl. Sci. 2020, 10, x FOR PEER REVIEW 9 of 31 of such conditions: temperature changes can dramatically cause fluid conductivity to alter, and consequently fluid temperature compensation mechanisms have to be applied to “correct” current/voltage measurements. Looking now to the evolution of acquisition systems, recent developments in the electronics field have made it possible to achieve frame acquisition speeds of 1000 frames/s [73,74]. In acquisition systems, a good excitation strategy is also vital, and these are typically current-based ones. However, the use of voltage-based excitation strategies, where a voltage-controlled source (typically with low output impedance) combined with an equal width pulse synthesiser, generates a subset of waveforms used for exciting the domain and to implement demodulation, have attracted increasing attention, since they produce less noisy data [73]. The use of voltage-controlled sources in the domain excitation causes electrical currents inside the domain that increase with the increase in suspension fluid conductivity; this fact is very important for highly conductive suspensions, above 2 S/m, which are quite typical of industrial processes [67,74]. Figure 2. EIT/ERT adjacent injection and measurement protocol for the first (A) and second (B) projections (reproduced with permission from Malmivuo et al. [71]). In addition to the shape and size of the electrodes, and the material from which they are made, when using ERT/EIT, other characteristics must be considered, such as the spacing between electrodes, and for EIT the frequency of the electrical excitation signal. Altogether, these characteristics condition the electric field distribution along the domain and determine the true measuring volume and overall sensitivity [74,75]. However, some of the negative effects of these characteristics can be overcome; for instance, the effect of the electrode dimensions can be reduced by using appropriate calibration procedures [76]. EIT, similar to ECT, can also be used to obtain disperse phase velocity by using a pair of data acquisition planes and applying cross-correlation algorithms to the acquired data [77,78]. However, to properly estimate velocity, a minimum frame rate acquisition of 100 frames per second is needed [78]. Demodulation has to be performed through software to overcome the problem of the low settling time of the low-pass filters. The way electrodes are combined to introduce the electrical signal into the domain is designated by the Injection Protocol. Thorough studies on injection and measurement protocols or patterns can be found in the literature. The most common strategy is referred to as the four-electrode method, where the signal is injected through some electrodes and the voltages are measured from the remaining electrodes. Different patterns for using this strategy can be found in the literature [79,80]: (a) Adjacent Pattern: where the signal is injected through neighbouring or adjacent electrodes and the resulting electrical potential differences are measured in the remaining pairs of electrodes, in a similar adjacent combination; (b) Opposite Pattern: In this pattern, the signal is injected through opposite electrodes and the electrical potential is measured with respect to one reference electrode adjacent to the injecting electrode. In addition to the above-mentioned strategies, which are the most frequently used, others have been employed, such as the Cross or Diagonal Pattern [80], which is a combination of the two described above, and the Trigonometric or Adaptive Pattern [81]. Figure 2. EIT/ERT adjacent injection and measurement protocol for the first ( A ) and second ( B ) projections (reproduced with permission from Malmivuo et al. [71]). Looking now to the evolution of acquisition systems, recent developments in the electronics field have made it possible to achieve frame acquisition speeds of 1000 frames/s [ 73 , 74 ]. In acquisition systems, a good excitation strategy is also vital, and these are typically current-based ones. However, the use of voltage-based excitation strategies, where a voltage-controlled source (typically with low output impedance) combined with an equal width pulse synthesiser, generates a subset of waveforms used for exciting the domain and to implement demodulation, have attracted increasing attention, since they produce less noisy data [ 73 ]. The use of voltage-controlled sources in the domain excitation causes electrical currents inside the domain that increase with the increase in suspension fluid conductivity; this fact is very important for highly conductive suspensions, above 2 S/m, which are quite typical of industrial processes [67,74]. In addition to the shape and size of the electrodes, and the material from which they are made, when using ERT/EIT, other characteristics must be considered, such as the spacing between electrodes, and for EIT the frequency of the electrical excitation signal. Altogether, these characteristics condition the electric field distribution along the domain and determine the true measuring volume and overall sensitivity [ 74 , 75 ]. However, some of the negative effects of these characteristics can be overcome; for instance, the effect of the electrode dimensions can be reduced by using appropriate calibration procedures [ 76 ]. EIT, similar to ECT, can also be used to obtain disperse phase velocity by using a pair of data acquisition planes and applying cross-correlation algorithms to the acquired data [77,78]. However, to properly estimate velocity, a minimum frame rate acquisition of 100 frames per second is needed [ 78 ]. Demodulation has to be performed through software to overcome the problem of the low settling time of the low-pass filters. Appl. Sci. 2020,10, 2355 16 of 30 pulp) and observed that as their size increased to accommodate higher production rates, operating problems related to liquor flow increased (as defined by the authors, liquor flow is the relative flow of liquor through the chip mass). Among the different problems, they found decreased counter current washing efficiency, accelerated corrosion, increased pulp non-uniformity and problems with chip movement. In a column arranged to operate as the wash section of a continuous digester, they used ERT to obtain images of the flow properties. For the tests carried, the packed bed was modelled using HDPE pellets of high-density polyethylene, having diameters of around 4.44 mm and thickness of around 1.91 mm. Regarding the liquid phase, a split between liquor flows was made to obtain ratios of axial (up flow in all cases), Q A , to radial flow, Q R , of 0.5, 1.0 and 2.0. The tests with wood chips showed that liquid flow adjacent to the vessel wall was significantly greater than at the centre of the vessel. This difference was not as significant for the HDPE disks for which the wall-affected zone is much smaller and below the spatial resolution of the ERT system. When a tracer was isokinetically added, ERT reconstructed images below the down comer showed that the entire trace leaved the down comer when Q R was twice the value of Q A ; nevertheless, for a 1:1 flow ratio, some tracer was detected above the screen level. Figure 4presents some images of the measurements conducted. Lee et al. [ 56 ] used Electrical resistance tomography (ERT) to evaluate the uniformity of liquor flow through a porous media in a laboratory model digester under scaled industrial conditions: a 1:15 geometrically scaled vessel, a vessel to particle diameter ratio of 93:1 to minimise wall effects, and close approximation of liquor superficial velocity and pore Reynolds number (pore Reynolds number is based on the capillary model of flow through packed beds and accounts for the porosity, ε , tortuosity, τ , and specific surface area presented to the flow). Their motivation departed from the literature, where the extent of delignification in batch digesters varied as a function of chip location in the vessel. This non-uniformity may be exacerbated by several factors but is commonly attributed to poor and/or non-uniform liquor flow through the reactor. Local interstitial flow velocities were measured using pixel-pixel cross correlation techniques. Tests were conducted to cover the range of typical liquor splits, Q T /Q B (the liquor volumetric flow ratio between the top and bottom of the digester) found in industrial units: 0.5, 1.0, and 2.0. The results showed that creation of uniform zones was readily achieved, although a certain flow rate was needed in the upper section to establish uniformity due to asymmetric liquor addition to the vessel. A mixing region was observed at the screen elevation in the outflow zone from the digester. The location of this interaction zone moved, depending on the relative liquor flux to the top and bottom of the vessel. When local flow velocities were measured using the tomographic data (equal liquor fluxes to the top and bottom circulation zones) no flow crossed the screen level. Thus, a stagnation point existed at the centre of the vessel near the plane of the extraction screens. However, by forcing liquor to pass through the screen plane (by directing the flow from the top to the bottom of the digester and bypassing the circulation screens, or vice versa) during part of the cook, the spatial dependence of kappa number (related to the degree of cooking) should be reduced. Hui and co-workers [ 127 ] studied the effect of non-Newtonian rheology, including the presence of yield stress, characteristic of pulp fibre in mixing operations, where regions of active motion around the impellers can be created with the cavern size affecting the quality of mixing attained. A scaled version of a commercial axial flow impeller was used in a standard side-entering configuration. Measured cavern diameters were compared against model predictions available in the literature. Since the studied pulp suspensions were opaque, cavern size had to be imaged using indirect techniques for which two methods were used, being ERT one of them (the other was ultrasonic Doppler velocimetry, UDV): ERT was the preferred one due to the marked decrease in the time required to acquire data. The shape of the cavern was best approximated as a truncated right-circular cylinder and, as expected, increasing the impeller speed increased the cavern volume. However, the development of the cavern volume with increasing impeller speed was not uniform, which was attributed to the interaction between the cavern and the vessel walls. A model was developed for predicting cavern development in yield stress fluids which included interaction with vessel boundaries. The proposed model predicts well the trend Appl. Sci. 2020,10, 2355 17 of 30 in cavern volume change with increasing impeller speed although the absolute cavern volume is only predicted to within 13%. Appl. Sci. 2020, 10, x FOR PEER REVIEW 16 of 31 from the literature, where the extent of delignification in batch digesters varied as a function of chip location in the vessel. This non-uniformity may be exacerbated by several factors but is commonly attributed to poor and/or non-uniform liquor flow through the reactor. Local interstitial flow velocities were measured using pixel-pixel cross correlation techniques. Tests were conducted to cover the range of typical liquor splits, QT/QB (the liquor volumetric flow ratio between the top and bottom of the digester) found in industrial units: 0.5, 1.0, and 2.0. The results showed that creation of uniform zones was readily achieved, although a certain flow rate was needed in the upper section to establish uniformity due to asymmetric liquor addition to the vessel. A mixing region was observed at the screen elevation in the outflow zone from the digester. The location of this interaction zone moved, depending on the relative liquor flux to the top and bottom of the vessel. When local flow velocities were measured using the tomographic data (equal liquor fluxes to the top and bottom circulation zones) no flow crossed the screen level. Thus, a stagnation point existed at the centre of the vessel near the plane of the extraction screens. However, by forcing liquor to pass through the screen plane (by directing the flow from the top to the bottom of the digester and bypassing the circulation screens, or vice versa) during part of the cook, the spatial dependence of kappa number (related to the degree of cooking) should be reduced. Figure 4. Tomographic images of a liquid with a certain concentration, a liquor, flowing through a packed bed of HDPE pellets at steady flow conditions, obtained using a set of sensing electrode rings located sequentially along a pipe (the planes are numbered from the bottom, starting with number 1). A step reduction in liquid conductivity was made in the flow exiting the down comer. Images display the liquid concentration distribution at steady-state flow ratios of QA/QR of 0.5, 1.0 and 2.0. A total liquid flux of 75 mL/s was used for QA/QR of 0.5 and 2.0, and of 50 mL/s for QA/QR equal to 1.0. In the tomographic diagrams, low conductivity regions are coloured blue and high conductivity regions are coloured red (Reproduced with permission from Ruzinsky et al. [126]). Hui and co-workers [127] studied the effect of non-Newtonian rheology, including the presence of yield stress, characteristic of pulp fibre in mixing operations, where regions of active motion around the impellers can be created with the cavern size affecting the quality of mixing attained. A scaled version of a commercial axial flow impeller was used in a standard side-entering configuration. Measured cavern diameters were compared against model predictions available in the literature. Since the studied pulp suspensions were opaque, cavern size had to be imaged using Figure 4. Tomographic images of a liquid with a certain concentration, a liquor, flowing through a packed bed of HDPE pellets at steady flow conditions, obtained using a set of sensing electrode rings located sequentially along a pipe (the planes are numbered from the bottom, starting with number 1). A step reduction in liquid conductivity was made in the flow exiting the down comer. Images display the liquid concentration distribution at steady-state flow ratios of Q A /Q R of 0.5, 1.0 and 2.0. A total liquid flux of 75 mL/s was used for Q A /Q R of 0.5 and 2.0, and of 50 mL/s for Q A /Q R equal to 1.0. In the tomographic diagrams, low conductivity regions are coloured blue and high conductivity regions are coloured red (Reproduced with permission from Ruzinsky et al. [126]). 4.3. Pulp Suspensions Characterisation Using ET: The Authors’ Approach 4.3.1. Overview of the Developed EIT System and of the Experimental Setup The authors of this paper have also been working on the improvement of a proprietary EIT system used in data acquisition of two-phase flows (solid/liquid). The reconstructed images aim to validate the CFD developed model [ 74 , 112 , 128 ], particularly when applied to the flow characterisation of pulp suspensions. To obtain improved sensitivity, the developed EIT system determines the differences in both the real and imaginary parts of the impedance over the domain. It is a portable device; consequently, it can be used in any location where 16 or 32 electrode planes have been installed on pipes. To avoid the cost and difficulties inherent in the construction of a current source, a voltage source was projected instead, a choice different from what is typical when designing electrical tomography apparatus. Jia et al. [ 67 ], demonstrated that the use of a potential voltage difference instead of current in the domain excitation procedure made it possible to surpass the limitations found when studying more conductive media. The developed EIT system has already been described previously [ 74 ]; nevertheless, together with a block diagram of the system architecture and its main modules, presented in Figure 5, a brief description is given here. A set of input/output boards with USB ports (DAQ boards) are used for interfacing the EIT system with a personal computer (PC). Using the analogue input channels of the DAQ boards, the electrical potentials, resulting from each applied stimulus, are read. The developed Appl. Sci. 2020,10, 2355 18 of 30 EIT system can acquire data at a rate that allows an image reconstruction rate of up to 4000 or 1000 frames per second, depending on the number of electrodes composing the sensing rings, either 16 or 32. To obtain the desired excitation signals, with the desired output amplitude and frequency (up to 25 MHz for this system), a programmable frequency voltage waveform was used. The authors also adopted an excitation strategy based on a voltage signal instead of a current, as has already been mentioned, as the use of this approach resulted in less noisy data [ 73 ]. In parallel with the generation of the stimulus signal, two others were also created, both with unitary magnitude, but where one is in phase with the excitation signal and the other is 90 degrees out of phase. To demodulate the read electrical potentials, each of them is multiplied simultaneously by both in phase and out of phase unitary produced signals; the signals outputted by the multipliers go through low pass filters, enabling only the constant terms to be extracted from them, which are the terms needed to calculate the impedance distribution in the domains. A set of analogue bidirectional multiplexers are used for routing of all signals: the excitation is directed to the desired electrode pair, and the read electric potentials to the multipliers. The digital outputs of the DAQ boards are used to address the multiplexers. This EIT system was used to characterise pulp suspensions in flow tests conducted in a pilot rig described elsewhere [ 129 ] and schematically shown in Figure 6. In short, the test section consisted of a horizontal PVC pipe, 7.5 m in length and 100 mm internal diameter. Appl. Sci. 2020, 10, x FOR PEER REVIEW 18 of 31 Figure 5. Block diagram of the actual tomographic system main modules. Figure 6. Pilot rig schematic. The flow rate could be regulated through the manipulation of two valves located after the pump and was measured by an electromagnetic flow meter. Appropriate lengths were inserted before and after the test section to account for entrance and exit effects. The EIT electrodes rings (constructed in Teflon), see Figure 7, consisted of 32 Titanium-based electrodes (each with a diameter of 5 mm which was previously optimised experimentally), circumferentially and equally spaced, and mounted in a Perspex tube inserted in the rig between the measuring pressure taps. For all tests reported here, an excitation frequency of 10 kHz with 2Vpp amplitude was imposed, and opposite injection and adjacent measuring protocols were used. For image reconstruction, the open-source software EIDORS [130], considering direct differential measured voltages and using a structured Mesh consisting of 2304 linear elements and 1201 nodes, was used. EIDORS implements a non-linear back projection method using a regularised algorithm (Tikhonov regularisation). To solve the forward problem the Complete Electrode model (CEM) was chosen [131]. In the reconstructed images of the conductivity distribution in the pipe cross-section, electrode number 1 corresponds to the top position in the cross-section of the tube; a darker blue colour corresponds to a lower conductivity region, where the concentration of the fibres is higher, while yellow or red colours identify higher conductivities regions, where the fibres are present in lower concentrations. Figure 5. Block diagram of the actual tomographic system main modules. Appl. Sci. 2020, 10, x FOR PEER REVIEW 18 of 31 Figure 5. Block diagram of the actual tomographic system main modules. Figure 6. Pilot rig schematic. The flow rate could be regulated through the manipulation of two valves located after the pump and was measured by an electromagnetic flow meter. Appropriate lengths were inserted before and after the test section to account for entrance and exit effects. The EIT electrodes rings (constructed in Teflon), see Figure 7, consisted of 32 Titanium-based electrodes (each with a diameter of 5 mm which was previously optimised experimentally), circumferentially and equally spaced, and mounted in a Perspex tube inserted in the rig between the measuring pressure taps. For all tests reported here, an excitation frequency of 10 kHz with 2Vpp amplitude was imposed, and opposite injection and adjacent measuring protocols were used. For image reconstruction, the open-source software EIDORS [130], considering direct differential measured voltages and using a structured Mesh consisting of 2304 linear elements and 1201 nodes, was used. EIDORS implements a non-linear back projection method using a regularised algorithm (Tikhonov regularisation). To solve the forward problem the Complete Electrode model (CEM) was chosen [131]. In the reconstructed images of the conductivity distribution in the pipe cross-section, electrode number 1 corresponds to the top position in the cross-section of the tube; a darker blue colour corresponds to a lower conductivity region, where the concentration of the fibres is higher, while yellow or red colours identify higher conductivities regions, where the fibres are present in lower concentrations. Figure 6. Pilot rig schematic. The flow rate could be regulated through the manipulation of two valves located after the pump and was measured by an electromagnetic flow meter. Appropriate lengths were inserted before and Appl. Sci. 2020,10, 2355 19 of 30 after the test section to account for entrance and exit effects. The EIT electrodes rings (constructed in Teflon), see Figure 7, consisted of 32 Titanium-based electrodes (each with a diameter of 5 mm which was previously optimised experimentally), circumferentially and equally spaced, and mounted in a Perspex tube inserted in the rig between the measuring pressure taps. For all tests reported here, an excitation frequency of 10 kHz with 2V pp amplitude was imposed, and opposite injection and adjacent measuring protocols were used. For image reconstruction, the open-source software EIDORS [ 130 ], considering direct differential measured voltages and using a structured Mesh consisting of 2304 linear elements and 1201 nodes, was used. EIDORS implements a non-linear back projection method using a regularised algorithm (Tikhonov regularisation). To solve the forward problem the Complete Electrode model (CEM) was chosen [ 131 ]. In the reconstructed images of the conductivity distribution in the pipe cross-section, electrode number 1 corresponds to the top position in the cross-section of the tube; a darker blue colour corresponds to a lower conductivity region, where the concentration of the fibres is higher, while yellow or red colours identify higher conductivities regions, where the fibres are present in lower concentrations. Appl. Sci. 2020, 10, x FOR PEER REVIEW 19 of 31 Figure 7. Image of one of the 32 Titanium electrodes test rings. 4.3.2. Overview of the Pulp Suspension Flow Imaging For the tests depicted in this section, industrial pulp suspensions of Pine and Eucalyptus fibres, both collected in a local paper mill, were firstly squeezed and washed in water, and the fibres were then resuspended in tap water doped with NaCl with adjusted electrical conductivity. The behaviours of Eucalyptus pulp suspensions (average Fibre length of 0.706 mm) and of Pine pulp suspensions (average Fibre length of 2.18 mm) were compared, as in Figure 8, under the conditions listed in Table 1. Table 1. Table of concentrations and flow velocities tested. Type of Pulp Concentration (wt/wt) Flow Velocities (m/s) Suspension Conductivity (mS·cm−1) Eucalypt 2.35 0.5, 1, 1.5, 2, 2.5, 3 1.460 1.5 0.5, 1, 1.5, 2, 2.5, 3 1.150 1 0.5, 1, 1.5, 2, 2.5, 3 1.128 Pine 2.35 0.5, 1, 1.5, 2, 2.5, 3 1.471 1.5 0.5, 1, 1.5, 2, 2.5, 3 1.461 1 0.5, 1, 1.5, 2, 2.5, 3 1.435 Figure 7. Image of one of the 32 Titanium electrodes test rings. 4.3.2. Overview of the Pulp Suspension Flow Imaging For the tests depicted in this section, industrial pulp suspensions of Pine and Eucalyptus fibres, both collected in a local paper mill, were firstly squeezed and washed in water, and the fibres were then resuspended in tap water doped with NaCl with adjusted electrical conductivity. The behaviours of Eucalyptus pulp suspensions (average Fibre length of 0.706 mm) and of Pine pulp suspensions (average Fibre length of 2.18 mm) were compared, as in Figure 8, under the conditions listed in Table 1. Table 1. Table of concentrations and flow velocities tested. Type of Pulp Concentration (wt/wt) Flow Velocities (m/s) Suspension Conductivity (mS·cm−1) Eucalypt 2.35 0.5, 1, 1.5, 2, 2.5, 3 1.460 1.5 0.5, 1, 1.5, 2, 2.5, 3 1.150 1 0.5, 1, 1.5, 2, 2.5, 3 1.128 Pine 2.35 0.5, 1, 1.5, 2, 2.5, 3 1.471 1.5 0.5, 1, 1.5, 2, 2.5, 3 1.461 1 0.5, 1, 1.5, 2, 2.5, 3 1.435 Before analysing the obtained suspensions flow images, we must reinforce that the tomography’s reflect the voltage differences measured. For fibre concentrations of around 2.35%, the fibre distributions in the cross-pipe section are diverse for the two types of fibres. Concerning the Pine suspension, a plug type flow behaviour is already visible at the lowest velocity, with an increase in the plug (size and Appl. Sci. 2020,10, 2355 20 of 30 intensity of the blue colour) when the velocity increases to 2 m/s (see Figure 8a,c,e). However, a further increase in the velocity up to 3 m/s leads to breakage of the plug, as in Figure 8g), as was expected for this type of long and stifffibre. Appl. Sci. 2020, 10, x FOR PEER REVIEW 20 of 31 Figure 8. EIT images obtained with a concentration of ≅2.35% (w/w) of Pine (left column: (a,c,e,g)) and Eucalyptus (right column: (b,d,f,h)) fibres with increasing flow velocity (0.5 m/s, 1m/s, 2m/s and 3m/s, top to bottom). Before analysing the obtained suspensions flow images, we must reinforce that the tomography’s reflect the voltage differences measured. For fibre concentrations of around 2.35%, the fibre distributions in the cross-pipe section are diverse for the two types of fibres. Concerning the Pine suspension, a plug type flow behaviour is already visible at the lowest velocity, with an increase in the plug (size and intensity of the blue colour) when the velocity increases to 2 m/s (see Figure 8a,c,e). However, a further increase in the velocity up to 3m/s leads to breakage of the plug, as in Figure 8g), as was expected for this type of long and stiff fibre. When analysing the Eucalyptus fibre suspension behaviour, at the lowest velocity the fibres remain more or less dispersed in the cross-pipe section (see Figure 8b), and the plug only starts to be visible at the velocity of 1 m/s, becoming denser at 2 m/s (Figure 8d,f). Similarly to the Pine fibres, a further increase in the Eucalyptus suspension velocity up to 3 m/s leads to plug breakage (Figure 8h). These differences in the reconstructed EIT images for this concentration are related to the physical properties of the fibres, in particular the fibre length and stiffness of the two pulps that originate different interactions between fibres when the suspensions flow in the tube. For a fibre concentration of around 1% (EIT images depicted in Figure 9), it was not possible to collect consistent data for the Eucalyptus fibres, and so Figure 9 depicts EIT images only for Pine fibre suspensions. For these suspensions, and by comparison with Figure 8, it is evident that at lower concentrations and 1 2 3 b) d) f) g) h) e) c) a) Figure 8. EIT images obtained with a concentration of  2.35% (w/w) of Pine (left column: ( a , c , e , g )) and Eucalyptus (right column: ( b , d , f , h )) fibres with increasing flow velocity (0.5 m/s, 1m/s, 2m/s and 3m/s, top to bottom). When analysing the Eucalyptus fibre suspension behaviour, at the lowest velocity the fibres remain more or less dispersed in the cross-pipe section (see Figure 8b), and the plug only starts to be visible at the velocity of 1 m/s, becoming denser at 2 m/s (Figure 8d,f). Similarly to the Pine fibres, a further increase in the Eucalyptus suspension velocity up to 3 m/s leads to plug breakage (Figure 8h). These differences in the reconstructed EIT images for this concentration are related to the physical properties of the fibres, in particular the fibre length and stiffness of the two pulps that originate different interactions between fibres when the suspensions flow in the tube. For a fibre concentration of around 1% (EIT images depicted in Figure 9), it was not possible to collect consistent data for the Eucalyptus fibres, and so Figure 9depicts EIT images only for Pine fibre suspensions. For these suspensions, and by comparison with Figure 8, it is evident that at lower concentrations and lower velocities (in this case 0.5 m/s), Pine fibres migrate and become agglomerated in a central plug, while for a more a more concentrated suspension, the fibres (Figure 8) occupy a larger region of the pipe cross-section. However, at more diluted suspensions, the strength of Appl. Sci. 2020,10, 2355 21 of 30 the interconnection of the fibres is lower and the plug breakage occurs at lower velocities as the velocity increases. Appl. Sci. 2020, 10, x FOR PEER REVIEW 21 of 31 lower velocities (in this case 0.5 m/s), Pine fibres migrate and become agglomerated in a central plug, while for a more a more concentrated suspension, the fibres (Figure 8) occupy a larger region of the pipe crosssection. However, at more diluted suspensions, the strength of the interconnection of the fibres is lower and the plug breakage occurs at lower velocities as the velocity increases. Figure 9. EIT images obtained with a concentration of ≅1% (w/w) of Pine fibres (a–d) for 0.5 m/s, 1 m/s, 1.5 m/s and 2.5 m/s, respectively. It must be emphasised at this point that environmental factors that were not easily kept under control may affect the colour mapping of the reconstructed images. In particular, the temperature changes significantly during continuous testing due to the energy input from the pump. Moreover, since cooling of the system was performed by passing tap water through a tube immersed in the suspension tank, the temperature of this cooling water changes significantly with respect to the seasons, which also affects the temperature in the rig. Changes in the suspension temperature affect the viscosity, but mainly the electric conductivity of the suspension and, consequently, the colour map of the reconstructed images is affected. 4.3.3. Fibre Flow Velocity Estimation By using two EIT planes and cross-correlation, EIT can also be used to determine the velocity profile in a pipe. For example, Mosorov et al. [132] used a best-correlated pixel method to measure both axial and radial velocities simultaneously. Independently of the procedure used, a marker is needed. As a first approach, the introduction of a small fraction of conductive or non-conductive solid spheres in the flow was used. After several attempts without success to distinguish these markers in the EIT images, the authors instead decided to use a trace of a highly concentrated NaCl solution as a marker. For the purpose of introducing this NaCl solution in the flowing fibre suspension, an injection system was developed and added to the existing pilot rig: it was located around 3 m before the first EIT electrode ring. This injector was built in such a way that the flow suffered almost no disturbance. For all the tests reported in this section, two 16 electrode titanium rings A and B were used, 30 cm apart, with ring A closer to the injector inlet, as shown in Figure 10. The procedure here reported for estimating velocity was different from those reported in the literature, but, nevertheless, with some similarities with the one described by Dong et al. [133]. Instead of using the complete set of EIT projections, the authors only made use of the first projection; first, the potential differences of each pair of the first projection with no marker injection in the flow were acquired, afterwards, they were subtracted to the measured voltage differences obtained during the runs with NaCl injection. 1 2 d) b) a) c) Figure 9. EIT images obtained with a concentration of  1% (w/w) of Pine fibres ( a – d ) for 0.5 m/s, 1 m/s, 1.5 m/s and 2.5 m/s, respectively. It must be emphasised at this point that environmental factors that were not easily kept under control may affect the colour mapping of the reconstructed images. In particular, the temperature changes significantly during continuous testing due to the energy input from the pump. Moreover, since cooling of the system was performed by passing tap water through a tube immersed in the suspension tank, the temperature of this cooling water changes significantly with respect to the seasons, which also affects the temperature in the rig. Changes in the suspension temperature affect the viscosity, but mainly the electric conductivity of the suspension and, consequently, the colour map of the reconstructed images is affected. 4.3.3. Fibre Flow Velocity Estimation By using two EIT planes and cross-correlation, EIT can also be used to determine the velocity profile in a pipe. For example, Mosorov et al. [ 132 ] used a best-correlated pixel method to measure both axial and radial velocities simultaneously. Independently of the procedure used, a marker is needed. As a first approach, the introduction of a small fraction of conductive or non-conductive solid spheres in the flow was used. After several attempts without success to distinguish these markers in the EIT images, the authors instead decided to use a trace of a highly concentrated NaCl solution as a marker. For the purpose of introducing this NaCl solution in the flowing fibre suspension, an injection system was developed and added to the existing pilot rig: it was located around 3 m before the first EIT electrode ring. This injector was built in such a way that the flow suffered almost no disturbance. For all the tests reported in this section, two 16 electrode titanium rings A and B were used, 30 cm apart, with ring A closer to the injector inlet, as shown in Figure 10. The procedure here reported for estimating velocity was different from those reported in the literature, but, nevertheless, with some similarities with the one described by Dong et al. [ 133 ]. Instead of using the complete set of EIT projections, the authors only made use of the first projection; first, the potential differences of each pair of the first projection with no marker injection in the flow were acquired, afterwards, they were subtracted to the measured voltage differences obtained during the runs with NaCl injection. Finally, the sum of the module of those differences was computed. For each suspension velocity and fibre concentration tested, several injection trials were performed to assess the reproducibility of this approach. The conditions used for the velocity estimation are summarised in Table 2. Some examples of data obtained with this approach are shown in Figure 11: three sets of NaCl injections for suspensions containing 1% and 2% of Eucalyptus fibres with a velocity of 0.25 m/s (Figure 11a,b). Appl. Sci. 2020,10, 2355 22 of 30 Appl. Sci. 2020, 10, x FOR PEER REVIEW 22 of 31 Figure 10. Schematic of the relative position of the EIT electrodes sections and NaCl injection inlet. Finally, the sum of the module of those differences was computed. For each suspension velocity and fibre concentration tested, several injection trials were performed to assess the reproducibility of this approach. The conditions used for the velocity estimation are summarised in Table 2. Some examples of data obtained with this approach are shown in Figure 11: three sets of NaCl injections for suspensions containing 1% and 2% of Eucalyptus fibres with a velocity of 0.25 m/s (Figure 11a,b). Table 2. Table of fibre concentrations and flow velocities used for the velocity estimation tests. Concentration (wt/wt) Flow Velocities (m/s) 3 0.25, 0.5 2 0.25, 0.5, 0.75, 1 1 0.25, 0.5, 0.75, 1 In that same figure, a graph showing the estimated velocity versus the measured suspension velocity is represented (Figure 11c). A quite reasonable estimation is obtained, and so the authors’ approach was validated. Moreover, as fibre concentration increases, the advance of the NaCl front is slowed due to the presence of more fibres, and also due to the structures that build up among fibres. This can be observed between Figures 11a and b, where for the higher fibre concentration, there is a delay in the detection of the NaCl in the second ring (ring B). Therefore, in the future, it can be used to evaluate those structures as a function of fibre concentration and flow velocity. Figure 10. Schematic of the relative position of the EIT electrodes sections and NaCl injection inlet. Table 2. Table of fibre concentrations and flow velocities used for the velocity estimation tests. Concentration (wt/wt) Flow Velocities (m/s) 3 0.25, 0.5 2 0.25, 0.5, 0.75, 1 1 0.25, 0.5, 0.75, 1 Appl. Sci. 2020, 10, x FOR PEER REVIEW 23 of 31 Figure 11. EIT injection test trials for fibre suspension concentrations of 1% and 2%, at the velocity of 0.25 m/s (a,b), and estimated velocity versus measured velocity for all fibre concentrations (c). 5. Conclusions Theoretical models developed to describe multiphase flows have been developed that have different assumptions requiring experimental verification for validation. For example, reliable data on velocity profiles in concentrated pulp suspensions during pipe flow is scarce in the literature. A multitude of techniques is now available, but a condition of non-disturbance of the flows is required in view of the complexity of multiphase systems. Non-invasiveness is therefore beneficial, and is characteristic of some of these techniques. Laser Doppler Anemometry (LDA), Magnetic Resonance Imaging (MRI), Ultrasonic Velocity Profiling (UVP) and Electrical Tomography were reviewed and compared for their energy requirements, spatial and temporal resolution, duration of observation times, spatial and temporal coherence and stability, cost of equipment, and field and scale of applications. The compromise between the advantages and limitations of these techniques led us to focus on the Electric Tomographic techniques, the principles behind them, and their application in the visualisation of pulp suspension flows. Our efforts were directed towards designing proprietary, low-cost, portable Electric Impedance Tomographic (EIT) hardware. To circumvent the difficulties and costs of designing a current source that suits the goals of the EIT system being designed, the choice was made to move from traditional EIT systems and, instead, to design a voltage source. The EIT system depicted here can be used with 16-electrode or 32-electrode rings. Examples of Eucalyptus and Pine pulp suspension flows are reported here, demonstrating the influence of fibre concentration 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 Injection 1 Injection 2 Injection 1 Injection 2 Ring A detection Ring B detection a) b) c) Figure 11. EIT injection test trials for fibre suspension concentrations of 1% and 2%, at the velocity of 0.25 m/s (a,b), and estimated velocity versus measured velocity for all fibre concentrations (c). In that same figure, a graph showing the estimated velocity versus the measured suspension velocity is represented (Figure 11c). A quite reasonable estimation is obtained, and so the authors’ approach was validated. Moreover, as fibre concentration increases, the advance of the NaCl front is slowed due to the presence of more fibres and to the structures that build up among fibres. This can Appl. Sci. 2020,10, 2355 23 of 30 be observed between Figure 11a,b, where for the higher fibre concentration, there is a delay in the detection of the NaCl in the second ring (ring B). Therefore, in the future, it can be used to evaluate those structures as a function of fibre concentration and flow velocity. 5. Conclusions Theoretical models developed to describe multiphase flows have been developed that have different assumptions requiring experimental verification for validation. For example, reliable data on velocity profiles in concentrated pulp suspensions during pipe flow is scarce in the literature. A multitude of techniques is now available, but a condition of non-disturbance of the flows is required in view of the complexity of multiphase systems. Non-invasiveness is therefore beneficial and is characteristic of some of these techniques. Laser Doppler Anemometry (LDA), Magnetic Resonance Imaging (MRI), Ultrasonic Velocity Profiling (UVP) and Electrical Tomography were reviewed and compared for their energy requirements, spatial and temporal resolution, duration of observation times, spatial and temporal coherence and stability, cost of equipment, and field and scale of applications. The compromise between the advantages and limitations of these techniques led us to focus on the Electric Tomographic techniques, the principles behind them, and their application in the visualisation of pulp suspension flows. Our efforts were directed towards designing proprietary, low-cost, portable Electric Impedance Tomographic (EIT) hardware. To circumvent the difficulties and costs of designing a current source that suits the goals of the EIT system being designed, the choice was made to move from traditional EIT systems and, instead, to design a voltage source. The EIT system depicted here can be used with 16-electrode or 32-electrode rings. Examples of Eucalyptus and Pine pulp suspension flows are reported here, demonstrating the influence of fibre concentration and type, as well as flow velocity, on the distribution of the fibres in the cross-sectional area of a horizontal pipe. The developed Electrical Impedance Tomography system was used to characterise the behaviour of pine and eucalyptus fibre suspension flow through a pipe at different fibre concentrations—1%, 1.5% and 2.35% (w/w)—and at varying flow velocities (between 0.5 and 3 m/s) based on the cross-section images obtained for totally developed flows. Especially for the pine fibre suspensions, it was clearly possible to identify a plug formation, which was probably due to the geometry of the pine fibres, which were longer than the eucalyptus ones. With respect to their usage for the estimation of suspension velocity, eucalyptus fibre suspensions were used in three different concentrations—1%, 2% and 3% (w/w)—and at velocities of between 0.25 and 1 m/s. The real suspension velocity measured by a flowmeter was quite reasonably estimated by the cross correlation of EIT images obtained in two consecutive electrode rings separated from one another by 300 mm. The developed EIT system allows 2D image reconstruction rates of up to 4000 or 1000 frames per second, depending on the number of electrodes used in the acquisition (16 or 32). With respect to the resolution of the obtained cross-section images, it is much higher than previously reported in the literature (which is evident by comparing Figure 3[53] and Figure 8). Thus far, the authors have used the developed equipment not only in the monitoring of fibre suspension flow and the estimation of fibre flow velocity, but also with great success in other diverse two-phase systems (solid–liquid [ 72 ], liquid–liquid [ 111 ], just to mention two). In the next step, the objective will be to upgrade the overall system (signal acquisition and image reconstruction algorithms) to visualise three-phase systems. Considering the results presented here, the usefulness of this technique and equipment is demonstrated and clearly justifies its use at an industrial scale, despite its lower spatial resolution when compared with other tomographic procedures. Author Contributions: All authors contributed equally to this work. All authors have read and agreed to the published version of the manuscript. Funding: This research was funded by COST Action FP1005, Fibre suspension flow modelling—A key for innovation and competitiveness in the pulp & paper industry, whose support is gratefully acknowledged. The authors thank the financial support of Fundaç ã o para a Ci ê ncia e Tecnologia through the project, contract Appl. Sci. 2020,10, 2355 24 of 30 PTDC/EQU-EQU/112388/2009 (COMPETE-01-0124-FEDER-015247), and contracts Pest-C/EME/UI0285/2013 and Pest-C/EQB/UI0102/2013, both financed by FCT/MCTES (PIDDAC) and co-financed by the European Regional Development Fund through the program COMPETE. Conflicts of Interest: The authors declare no conflict of interest. References 1. Steiner, G.; Podd, F. A non-invasive and non-intrusive ultrasonic transducer array for process tomography. In Proceedings of the XVIII IMEKO WORLD Congress on Metrolofy and Sustainable Development, Rio de Janeiro, Brazil, 17–20 September 2006; pp. 17–22. 2. Ijaz, U.Z.; Kim, J.-H.; Khambampati, A.K.; Kim, M.-C.; Kim, S.; Kim, K.-Y. Concentration distribution estimation of fluid through electrical impedance tomography based on interacting multiple model scheme. Flow Meas. Instrum. 2007,18, 47–56. [CrossRef] 3. Murphy, S.C.; York, T.A. Electrical impedance tomography with non-stationary electrodes. Meas. Sci. Technol. 2006,17, 3042–3052. [CrossRef] 4. Durst, F.; Melling, A.; Whitelaw, J.H. Principles and Practice of Laser-Doppler Anemometry; Academic Press: New York, NY, USA, 1981. 5. Chaouki, J.; Larachi, F.; Dudukovic, M.P. Non-invasive tomographic and velocimetric monitoring of multiphase flows. Ind. Eng. Chem. Res. 1997,36, 4476–4503. [CrossRef] 6. Takeda, Y. Velocity Profile Measurement by Ultrasonic Doppler Shift Method. In Proceedings of the National Heat Transfer Conference of the American Society of Mechanical Engineers, Philadelphia, PA, USA, 6–9 August 1989; pp. 155–160. 7. Haavisto, S.; Cardona, M.; Salmela, J.; Powell, R.; McCarthy, M.; Kataja, M.; Koponen, A. Experimental investigation of the flow properties and rheology of complex fluids in pipe flows by hybrid multiscale velocimetry. Exp. Fluids 2017,58, 158. [CrossRef] 8. Sorokin, M.; Strokina, N.; Eerola, T.; Lensu, L.; Karttunen, K.; Kalviainen, H. Image-based characterisation of the pulp flows. Pattern Recognit. Image Anal. 2016,26, 630–637. [CrossRef] 9. Heindel, T. A review of X-ray flow visualisation with applications to multiphase flows. J. Fluids Eng. 2011 , 133. [CrossRef] 10. Schlieper, G. Principles of Gamma ray densitometry. Met. Powder Rep. 2000,55, 20–23. [CrossRef] 11. Yenjaichon, W.; Grace, J.R.; Lim, C.J.; Bennington, C.P.J. Characterisation of gas mixing in water and pulp-suspension flow based on electrical resistance tomography. Chem. Eng. J. 2013 ,214, 285–297. [CrossRef] 12. Nakagawa, M.; Altobelli, S.A.; Caprihan, A.; Fukushima, E.; Jeong, E.K. Non-invasive measurements of granular flows by magnetic resonance imaging. Exp. Fluids 1993,16, 54–60. [CrossRef] 13. Sederman, A.J.; Gladden, L.F.; Mantle, M.D. Application of magnetic resonance imaging techniques to particulate systems. Adv. Powder Technol. 2007,18, 23–38. [CrossRef] 14. Müller, C.R.; Holland, D.J.; Sederman, A.J.; Mantle, M.D.; Gladden, L.F.; Davidson, J.F. Magnetic resonance imaging of fluidized beds. Powder Technol. 2008,183, 53–62. [CrossRef] 15. Li, T.Q.; Ödberg, L. Flow properties of cellulose fiber suspensions flocculated by cationic polyacrylamide. Colloids Surf. A Physiochem. Eng. Asp. 1996,115, 127–135. [CrossRef] 16. Li, T.Q.; Ödberg, L. Studies of Flocculation in Cellulose Fibre Suspensions by NMR Imaging. J. Pulp Pap. Sci. 1997,23, 401–405. 17. Li, T.Q.; Powell, R.L.; Ödberg, L.; McCarthy, M.J.; McCarthy, K.L. Velocity measurements of fiber suspensions by the nuclear magnetic resonance imaging method. Tappi J. 1994,77, 145–149. 18. Li, T.Q.; Seymour, J.D.; Powell, R.L.; McCarthy, M.J.; McCarthy, K.L.; Ödberg, L. Visualization of flow patterns of cellulose fiber suspensions by NMR imaging. Aiche J. 1994,40, 1408–1411. [CrossRef] 19. Li, T.Q.; Weldon, M.; Ödberg, L.; McCarthy, M.J.; Powell, R.L. Pipe flow behaviour of hardwood pulp suspensions studied by NMRI. J. Pulp Pap. Sci. 1995,21, 408–414. 20. Arola, D.F.; Li, T.Q.; Ödberg, L.; Powell, R.L.; McCarthy, M.J. Flow of pulp suspensions through a sudden expansion. Aiche J. 1998,44, 2597–2606. [CrossRef] 21. Seymour, J.D.; Maneval, J.E.; McCarthy, K.L.; McCarthy, M.J.; Powell, R.L. NMR velocity phase encoded measurements of fibrous suspensions. Phys. Fluids A Fluid Dyn. 1993,5, 3010–3012. [CrossRef] Appl. Sci. 2020,10, 2355 25 of 30 22. Slavica, R. Laser Doppler Anemometry and its Application in Wind Tunnel Tests. Sci. Tech. Rev. 2007 ,LVII, 64–76. 23. Ek, R.; Möller, K.; Norman, B. Measurement of velocity and concentration variations in dilute fiber/air suspensions using a laser Doppler anemometer. Tappi J. 1978,61, 49–52. 24. Kerekes, R.J.; Garner, R.G. Measurement of Turbulence in Pulp Suspensions by Laser Anemometry. Trans. Tech. Sect. 1982,CPPA 8, TR53. 25. Steen, M. The Application of Refractive Index Matching for Two-phase Flow Measurements in Turbulent Fibre Suspensions by Laser Doppler Anemometry. Nord. Pulp Pap. Res. J. 1989,4, 236–243. [CrossRef] 26. Steen, M. On Turbulence Structure in Vertical Pipe Flow of Fiber Suspensions. Nord. Pulp Pap. Res. J. 1989 ,4, 244–252. [CrossRef] 27. Andersson, S.; Rasmuson, A. Flow Measurements on a Turbulent Fibre Suspension by Laser Doppler Anemometry. Aiche J. 2000,46, 1106–1119. [CrossRef] 28. Pettersson, J.; Rasmuson, A. LDA Measurements on a Turbulent Gas/Liquid/Fibre Suspension. Can. J. Chem. Eng. 2004,82, 265–274. [CrossRef] 29. Pettersson, J.; Wikström, T.; Rasmuson, A. Near Wall Studies of Pulp Suspension Flow Using LDA. Can. J. Chem. Eng. 2006,84, 422–430. [CrossRef] 30. Takeda, Y. Velocity Profile Measurement by Ultrasonic Doppler Method. Exp. Fluid Sci. 1995 ,10, 444–453. [CrossRef] 31. Wiklund, J.; Pettersson, J.; Stading, M.; Rasmuson, A. A Comparative Study of UVP and LDA Techniques for Pulp Suspensions in Pipe Flow. Aiche J. 2006,52, 484–495. [CrossRef] 32. Ozaki, Y.; Kawaguchi, T.; Takeda, Y.; Hishida, K.; Maeda, M. High time resolution ultrasonic velocity profiler. Exp. Fluid Sci. 2002,26, 253–258. [CrossRef] 33. Hirsimäki, O. Determination of Radial Velocity Profile and Flow Disturbance of Pulp Suspension by Ultrasonic Echo Correlation. Pap. Ja Puu 1978,60, 95–97. 34. Karema, H.; Kataja, M.; Kellomäki, M.; Salmela, J.; Selenius, P. Transient Fluidisation of Fibre Suspension in Straight Channel Flow. In Proceedings of the TAPPI Internation Paper Physics Conference, San Diego, CA, USA, 26–30 September 1999; pp. 369–379. 35. Karema, H.; Salmela, J.; Tukiainen, M.; Lepomäki, H. Prediction of Paper Formation by Fluidisation and Reflocculation Experiments. In Proceedings of the 12th Fundamental Resesearch Symposium, Oxford, UK, 11 October 2001; pp. 559–589. 36. Wiklund, J.; Johansson, M.; Shaik, J.; Fischer, P.; Windhab, E.; Stading, M.; Hermansson, A.-M. In-Line Ultrasound based Rheometry of industrial and model suspensions flowing through pipes. In Proceedings of the Third International Symposium on Ultrasonic Doppler Methods for Fluid Mechanics and Fluid Engineering, Lausanne, Switzerland, 9–11 September 2002; pp. 69–76. 37. Claesson, J.; Rasmuson, A.; Wiklund, J.; Wikström, T. Measurement and analysis of flow of concentrated fiber suspensions through sudden expansion using UVP. Aiche J. 2013,59, 1012–1021. [CrossRef] 38. Kotz é , R.; Wiklund, J.; Haldenwang, R. Optimisation of Pulsed Ultrasonic Velocimetry system and transducer technology for industrial applications. Ultrasonics 2013,53, 459–469. [CrossRef] [PubMed] 39. Xu, H.; Aidun, C. Characteristics of fiber suspension flow in a rectangular channel. Int. J. Multiph. Flow 2005 , 31, 318–336. [CrossRef] 40. Dresxler, W.; Fujimoto, J. Optical Coherence Tomography: Technology and Applications; Springer: Berlin, Germany, 2008. 41. Kataja, M.; Haavisto, S.; Slamela, J.; Letho, R.; Koponen, A. Characterisation of micro-fibrilated cellulose fiber suspension flow using multi scale velocity profile measurements. Nord. Pulp Pap. J. 2017 ,31, 473–482. [CrossRef] 42. Fujimoto, J.; Pistris, C.; Boppart, S.; Brezinski, M. Optical coherence tomography: An emerging techonology for biomedical imaging and optical biopsy. Neoplasia 2000,2, 9–25. [CrossRef] 43. Gao, Q.; Wang, H.; Shen, G. Review on development of volumetric particle image velocimetry. Chin. Sci. Bull. 2013,58, 4541–4556. [CrossRef] 44. Fan, L.; Xu, N. Study on laminar flow containing fibre particles in a stirred tank using DPI. Powder Technol. 2017,320, 420–427. [CrossRef] 45. Fock, H.; Claesson, J.; Rasmusson, A.; Wikström, T. Near wall effect in the plug flow of pulp suspensions. Can. J. Chem. Eng. 2011,89, 1207–1216. [CrossRef]