International Journal of Optics
Volume 2016 (2016), Article ID 3582649, 19 pages
http://dx.doi.org/10.1155/2016/3582649

Research Article

Gas Bubbles Investigation in Contaminated Water Using Optical Tomography Based on Independent Component Analysis Method

Mohd Taufiq Mohd Khairi, Sallehuddin Ibrahim, Mohd Amri Md Yunus, and Mahdi Faramarzi

Department of Control and Mechatronics Engineering, Faculty of Electrical Engineering, Universiti Teknologi Malaysia, 81310 Skudai, Johor, Malaysia

Received 3 September 2015; Revised 18 March 2016; Accepted 13 April 2016

Academic Editor: Augusto Beléndez

Copyright © 2016 Mohd Taufiq Mohd Khairi et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

This paper presents the results of concentration profiles for gas bubble flow in a vertical pipeline containing contaminated water using an optical tomography system. The concentration profiles for the bubble flow quantities are investigated under five different flows conditions, a single bubble, double bubbles, 25% of air opening, 50% of air opening, and 100% of air opening flow rates where a valve is used to control the gas flow in the vertical pipeline. The system is aided by the independent component analysis (ICA) algorithm to reconstruct the concentration profiles of the liquid-gas flow. The behaviour of the gas bubbles was investigated in contaminated water in which the water sample was prepared by adding 25 mL of colour ingredients to 3 liters of pure water. The result shows that the application of ICA has enabled the system to detect the presence of gas bubbles in contaminated water. This information provides vital information on the flow inside the pipe and hence could be very significant in increasing the efficiency of the process industries.

1. Introduction

Information about the flow regime in a process vessel or pipelines is important so as to design an accurate, safe, and low cost conveying system in various applications such as chemical engineering and nuclear engineering [1]. The flow regime identification is also vital in evaluating the performance of a process system since incorrect analysis from the system can lead to the reduction of production rates. Hence, an appropriate instrumentation system is essential in providing vital information on the flow patterns. The tomography technique has been selected extensively in the industrial field as a tool to provide information about the phase and spatial distribution without interrupting the process flow [2]. Before being applied in the industry, the technique has been successfully implemented in the medical field where it is responsible for capturing images of body tissue and detecting tumour [3]. A tomography system consists of several parts such as sensors, signal conditioning circuit, a data acquisition system, and a computer. There are many types of sensors that have been used in process tomography, that is, electrical capacitance, electrical resistance, ultrasonic, and optical. In this paper, an optical sensing technique has been selected to be integrated into the tomography system due to its low cost, being straightforward, and having a better dynamic response [4]. Many researchers have published their result using the optical sensor in order to monitor the movement of gas, liquid, and solid. Idroas et al. [5] have successfully applied an optical tomography system for particle sizing identification using sample of beads and irregular shaped nut. Rzasa and Plaskowski [6] have presented gas bubbles measurement using an optical tomography system in a vertical aeration column. Ibrahim et al. [7] proposed a concentration measurement using halogen bulbs as light projectors to reconstruct images of small and large gas bubbles. This paper describes the analysis of gas bubble detection in a pipeline filled with contaminated water. Investigation on gas bubble behaviour in industry such as food and pharmaceutical is vital since the existence of unwanted bubbles can reduce the quality of product [8]. Previous research on gas bubbles measurement in the tomography field was conducted by Rzasa [9], Jin et al. [10], and Ayob et al. [11] which used pure water as a liquid. Generally, the unwanted gas bubbles are in the form of opaque liquid such as oil and paint. Therefore, this paper proposed a technique for imaging the gas bubbles in which the investigation was carried out using opaque liquid.

2. Mathematical Modelling

Mathematical modelling is an essential part in optical tomography. This section discusses the optical attenuation model and the application of independent component analysis. The first part will discuss the optical attenuation modelling for the presence of gas bubbles in contaminated water. The second part elaborated the application of ICA method to detect the distribution of gas bubbles.

2.1. Optical Attenuation Model

Light experienced attenuation when it traverses from one medium to another medium. This attenuating process is due to scattering and refraction of light when it travels through a different medium such as water and gas. This paper highlighted the absorption and attenuation model and neglected the light scattering effect since it is difficult to perform the modelling due to the random shape [12]. Beer-Lambert’s Law is a popular law regarding the investigation of attenuation and absorption of the light energy. The law can be expressed as

where is the receiving sensor voltage, is the voltage of the transmitter, is the attenuation coefficient of water medium (mm⁻¹), is the attenuation coefficient of air medium (mm⁻¹), is the path length of water (mm), and is the path length of air (mm). The reference attenuation coefficients for water and air are = 0.0287 mm⁻¹ and = 0.0142 mm⁻¹ [12]. The value of pipe diameter is = 100 mm and the length of small bubble is assumed to be = 5 mm. The value of is being set to 5 volts (V) which was taken from the transmitter value programmed using the Peripheral Interface Controller (PIC) Microcontroller. Two cases are investigated: water with no gas bubble presence and water with gas bubble presence.

2.1.1. Water with No Gas Bubble Presence

The arrangement of the transmitter and receiver when there is no gas bubble is illustrated in Figure 1. Light from the transmitter traversed water without the interruption of bubble. The value of the attenuation coefficient for water is entered in (2) along with 100 mm of water path length. The value of the receiver voltage is obtained by

Figure 1: The light path in water when there is no gas bubble.

2.1.2. Water with Gas Bubble Presence

The aim of this modelling is to determine the change of parameters and when a single gas bubble intercepts the transmitting light as shown in Figure 2. In this case, the illustration neglects the reflection and scattering effect when light beam passed through the bubble. The value of the receiving voltage when light passed through bubble and water can be determined as follows:

The values of when there is no gas bubble and when gas bubble is present are compared. The existence of bubble in water resulted in a higher value of as indicated in (5) compared to when no bubble is present as indicated in (3). Light experienced less attenuation when gas bubble is present. When there is a bubble, light needs to pass two dissimilar mediums: water and gas. Hence the exponential value is higher compared to when no bubble exists. The values of and are predicted to increase when more and bigger size bubbles intercepted the light before it reaches the receiver’s side. The illustration and an example of calculation of for double bubble model can be seen in Figure 3. Related equations are given in

Figure 2: The condition of light path in water with gas bubble presence in contaminated water.

Figure 3: The light path when double model is present in contaminated water.

2.2. Independent Component Analysis

Blind source separation (BSS) is a technique which separates a set of individual source signals from their mixtures without any specific knowledge as with very limited information about the source signals [13, 14]. An adaptive system is required to solve this problem. Independent component analysis (ICA) is a method which is extensively used for extracting the mixture signals to individual signal [13]. The method is developed by Jutten and Herault at the end of 1980s for retrieving the original sound from the mixed sound signals in the microphone [15]. The method is different from the Principle Component Analysis (PCA) method where PCA is focused in finding the principle components which consist of linear combination of the observed variables, whereas ICA is concerned with discovering the independent components from the observed variables [16]. ICA has been applied in many areas such as telecommunications [17], image processing [18], gas-liquid phase separation in tomography [14], and biomedical applications [19]. A simple example on the application of ICA is shown in Figure 4. In Figure 4, the voice signals from two speakers are mixed together in each microphone. The crowded voices signals also add up and resulted in noise being produced. By using ICA, the original voice signal for both speakers can be acquired.

Figure 4: A simple example for understanding the ICA concept [13].

The general model of ICA is described as

where is matrix of mixture of source signals, is mixing matrix, and is matrix denoting a source signal. The aim of ICA is to get the value of mixing matrix and/or source signal matrix given only realizations of the observation matrix of . After the ICA process was executed, (10) converts to (11) where is the unmixing matrix and equals :

The illustration in Figure 5 can assist in understanding the concept of ICA technique. After going through the ICA process, the estimation of source signal is produced and it looks similar to the original signal, . Moreover, the ICA has formed the separating matrix which was produced alongside the output signal at the output side.

Figure 5: ICA block diagram.

Among the most widely used ICA algorithms are FastICA and Infomax [20]. FastICA provides recovery of independent sources by employing the higher order statistics and the estimation process is done one by one [13]. The algorithm separates the source signals based on their non-Gaussianity and it finds one independent component at a time instead of solving the mixing matrix [20]. The FastICA algorithm is as follows [21]:

(1)Choose an initial weight vector, .(2)Let

where is a temporary variable used to calculate weight vector, is averaging over all column-vectors of matrix , is a transpose of weight vector, is input data, and is the derivative of nonquadratic used in the contrast function for solving ICA problem.(3)Let

(4)If not converged, go back to (7).

The expression “converged” means the value is either old or new for in the same direction. The advantages of FastICA are that it is parallel and simple for computation process and requires small memory size [22]. For the Infomax algorithm, it is based on maximizing the information transferred in a network of nonlinear units and it is proposed by Bell and Sejnowski [23]. The algorithm is derived from a maximization principle between inputs and nonlinear outputs. The entropy can be described as

is the conditional entropy where the entropy of is conditional on . and are the entropy of and is mutual information. From this equation, maximizing the joint entropy of the outputs amount has led to minimising the mutual information [13]. The disadvantage of Infomax is that the function does not converge properly if step size is wrongly chosen [24].

This research emphasized the separation of the transmitters signals which were mixed in the receivers. In ICA, there is some ambiguity that should be handled properly. Firstly, the separated signals indicated that the amplitude of the signals is not identical as the source signals. Secondly, the separating matrix is formed in the arbitrary row [20]. A similar parameter from input and output that can be investigated in ICA is the frequency value. Therefore, each transmitter is set to have a different duty cycle in order to facilitate the rearrangement of the separating matrix row. An in-depth explanation about the rearrangement procedure is discussed in Section 4. For simplicity, the modelling for the algorithm is described using four transmitters and four receivers as shown in Figure 6 while the actual system consists of eighteen transmitters and eighteen receivers. This is because the explanation of the ICA concept using eighteen pairs of sensors is complex since the matrix equation is too long. Hence, in this section only four pairs of sensor and four absorption coefficients are chosen and are sufficient to show how ICA works. Transmitters 1, 2, 3, and 4 are represented as TX1, TX2, TX3, and TX4, respectively, while receivers 1, 2, 3, and 4 are represented as RX1, RX2, RX3, and RX4, respectively. The attenuation coefficients are represented as , , , and .

Figure 6: Model of four-sensor pair with four pixels.

In this model, each receiver is supposed to detect light from all transmitters. Take, for example, receiver RX1. The four signals received at RX1 have been mixed where each light has its own characteristics such as light length and absorption coefficient. The equation of receiver RX1 voltage can be written as

where is transmitter TX1 voltage, is transmitter TX2 voltage, is transmitter TX3 voltage, and is transmitter TX4 voltage, while is attenuation coefficient for material 1, is attenuation coefficient for material 2, is attenuation coefficient for material 3, and is attenuation coefficient for material 4, whereas is length 1, is length 2, and is length 3. Voltage transmitter () is a source signal when we analyse the light received by receiver. For example, as shown in Figure 6, receiver 1 (RX1) detects four lights from transmitter 1 (TX1), transmitter 2 (TX2), transmitter 3 (TX3), and transmitter 4 (TX4). The distribution of bubbles in a pipeline is estimated based on the parameters of and receiver output voltage as indicated in (16). Equation (12) can be arranged as

The equations for receiver RX2 voltage, receiver RX3 voltage, and receiver RX4 voltage are simplified as

Equations (13) to (14) are transformed into matrix as shown in (15), where is a matrix of mixture of source signals, represented by the voltage receiver (). is denoted as a mixing matrix and voltage transmitter () is signified as a source signal :

ICA is utilized to separate the mixture of source signals where in (16) the matrix positioned of source signals represented as is interchanged with denoted by matrix of mixture of source signals . The separating matrix is formed to replace the mixing matrix where theoretically is obtained by the inversing process of :

In this case, the parameter considered for computing the concentration profile of bubble flow is the matrix . However, the matrix cannot be acquired directly from the ICA method. Hence, the related parameter regarding is supported where the matrix is obtained as an inverse process in order to find the matrix .

3. Hardware Design

In the optical tomography system, eighteen transmitters and eighteen receivers are mounted around a 100-millimeter (mm) diameter pipe by using a sensor jig. The transparent pipe is made from an acrylic material which allows visual observation of the flow process. The transmitters are positioned next to the receivers since the system is modelled based on fan beam light projection. The infrared Light-Emitting Diode (LED) model TSAL6200 is used as the transmitters where it emits an invisible light. The invisible light has a smaller value of wavelength compared to the visible light in which the range is 750 nanometers (nm) and above. The photodiode model BPV10NF is selected due to its ability to detect the infrared light from the transmitter. A frequency of 50 Hertz (Hz) is chosen for the transmitting circuit frequency, and each transmitter is set using different pulse durations for switching time by using a Peripheral Interface Controller (PIC) Microcontroller. Transmitter 1 (TX1) is set to 1 millisecond (ms) and TX2 for 2 ms and followed until 18 ms for TX18. The purpose for doing this is to facilitate the determination of the transmitter’s signal detected by the photodiode where it will be explained in Section 4. The outputs from the receivers are connected to a signal conditioning circuit in order to convert the receiving light into voltage and for voltage amplification. The amplified voltage signals are converted into the digital form using the data acquisition system (DAQ) model U2331 manufactured by Agilent Technologies Inc. It has 3 megasamples per second (MS/s) of sampling rate, 64 single ended/32 differential inputs for analog input channels, and 24-bit programmable TTL input/output channels. The experiment is performed offline. The voltage output from receiver is converted into a digital form using the DAQ and further analysis is performed using the LabVIEW software, which provided the graphical based system for the ICA algorithm. The concentration profiles of the bubble’s condition are displayed on a computer. The total time to obtain the final result of concentration profile is around 1 minute. An illustration of the overall design of the optical tomography system is shown in Figure 7.

(a)

(b)

Figure 7: (a) The overall design of the optical tomography system. (b) The preparation of the apparatus.

4. Experimental Procedure

This section will elaborate the approach to calculate and analyse the result of bubble flow using ICA method. The steps of the experiment are according to these stages:

(1)Each transmitter has a unique value where it is programmed with different pulse duration value. Hence, the switching times are different from each other. It is assumed that there are eighteen signals from transmitters to be detected by each receiving sensor. Then, the ICA method is applied in order to separate the eighteen signals in order to get a single transmitter signal in each output indicator. The pulse duration value for the transmitter can be from 1 ms to 18 ms, which is set in the transmitter circuit.(2)A transmitter’s signal at each output indicator also brings along a row of unmixed matrix . The whole 18 rows for are produced in arbitrary rows and rearrangement of the matrix’s row is necessary according to the order of transmitters pulse duration. For example, a signal with 10 ms is shown in output indicator; hence the value at row 1 will be shifted to row 10. In other words, the rearrangement process for the row depends on the value of pulse duration displayed at each output indicator.(3)The new matrix is formed using the ICA method. The matrix of is emphasized in the steps since it contains an important parameter that has a relationship with the mixing matrix . In the ICA theory, . Hence, the matrix will be inversed to get the value of which consists of matrix 18 × 18 and expansion of . The matrix of is named as the turbidity factor .(4)The 18 columns and 18 rows contained in matrix represented 18 receivers and 18 transmitters, respectively. Each value is represented as where is the transmitter number and is the receiver number. For example, by referring to Table 1, has the value 2.782. It indicates that when light is transmitted from transmitter 2 (TX2) to receiver 1 (RX1), it has a turbidity factor of 2.782.(5)The value at each pixel contained in the cross section of the pipe has to be determined. The value is the affiliation between and light’s path length, where light crossed a pixel. An example of the light path along the pixels can be seen in Figure 8. For light transmitted from transmitter TX2 to receiver RX1, two pixels and are traversed by the light. The light’s lengths in both pixels are measured as 0.4 cm and 0.5 cm, respectively, with the total light’s length as 0.9 cm. The value is computed in order to know the quantity of light that has crossed pixels and . The calculation started from pixel , followed by other pixels contained in the pipe’s cross section by using

where is the pixel at row and column , is the path length of light that has crossed the pixel in centimeter (cm), is the total length from the transmitter to the receiver, and is the turbidity factor. An example for calculating the value for pixel can be shown as

(6)From Figure 8, pixel is intercepted by four lights, which are from transmitters TX2, TX3, TX4, and TX5. By referring to Table 1, the value for light from transmitter TX3 to receiver RX1 is . Then, the value is recalculated using (18) based on its light’s path length. The calculation process for getting the value continued to other light sources, namely, transmitter TX4 and transmitter TX5.(7)The transmitter’s light also could cross pixel when it is transmitted to receiver RX2 until receiver RX18. Hence, it is necessary to recognise whether light from the transmitters has crossed the pixel or not. If it has, the values from all transmitters are totalled up and divide into the number of lights that has passed through pixel . The result from this calculation is named as the average value . The value is the final value to be computed in order to estimate the location of bubbles at certain pixels. Detail of the calculation is shown in Table 2.

Table 1: Inverse matrix of 25% of air opening flow in 25 mL of water.

Row
Column

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18

1
0.369
1.706
0.692
0.077
0.05
0.335
0.139
0.918
0.439
0.041
0.828
0.08
0.248
1.114
1.123
2.595
1.706
0.978

2
2.782
5.009
2.06
0.295
0.039
0.208
0.003
1.384
0.275
0.481
2.333
1.678
1.607
2.722
1.29
7.143
5.219
3.903

3
1.647
3.761
2.121
0.125
0.077
0.097
0.791
1.407
0.381
0.89
0.918
1.247
0.714
2.083
1.345
4.635
4.479
2.775

4
2.812
1.763
0.797
0.114
0.109
0.12
0.183
0.06
0.336
0.007
0.321
0.063
0.641
0.465
0.898
3.436
1.275
1.264

5
0.321
2.998
1.212
0.008
0.002
0.19
0.032
1.099
0.573
0.088
1.181
0.684
0.93
1.829
2.363
4.242
2.499
1.654

6
0.002
1.648
0.746
0.081
0.153
0.443
0.224
1.074
0.654
0.002
1.25
0.852
0.573
1.227
1.394
2.245
2.667
1.212

7
0.07
0.985
0.393
0.284
0.372
0.597
0.123
0.505
0.2
0.653
0.605
0.639
0.355
1.782
1.353
1.436
2.448
0.119

8
2.568
5.099
2.003
0.489
0.395
0.606
0.078
1.302
0.607
1.392
1.466
1.332
0.569
2.464
2.074
7.335
6.613
4.337

9
4.519
7.131
3.013
0.326
0.326
0.671
0.324
1.897
0.004
1.225
2.591
0.988
0.916
3.111
0.617
10.786
8.416
5.695

10
2.934
4.501
1.44
0.354
0.095
0.247
0.06
0.612
0.085
0.736
1.272
0.833
0.869
1.574
0.137
6.034
3.741
2.812

11
2.637
4.568
2.215
0.534
0.266
0.588
0.119
0.526
0.103
1.031
1.181
0.859
0.663
1.693
0.033
6.054
4.261
2.341

12
1.282
3.864
1.875
0.413
0.537
0.764
0.101
0.453
0.186
0.717
0.911
0.657
0.534
1.25
0.518
4.831
3.379
2.337

13
1.018
2.277
1.045
0.259
0.206
0.586
0.087
0.013
0.362
0.622
0.334
0.524
0.693
0.846
0.385
2.59
1.54
1.07

14
2.523
3.722
1.544
0.228
0.079
0.347
0.141
0.611
0.217
1.141
1.227
0.951
0.76
2.282
0.543
5.418
3.854
2.447

15
0.728
2.226
0.936
0.335
0.011
0.093
0.009
1.954
0.252
0.266
1.098
0.733
0.439
1.708
1.698
2.726
2.803
1.828

16
1.02
1.797
0.933
0.088
0.711
0.184
0.143
0.411
0.121
0.548
0.036
0.226
1.185
0.876
0.347
2.357
3.662
1.622

17
1.164
0.334
0.131
0.12
0.159
0.976
0.235
0.464
0.373
1.086
0.175
0.105
0.725
0.198
1.029
0.802
1.385
0.677

18
0.243
4.09
1.562
0.768
0.805
0.812
1.465
2.311
0.716
1.571
0.867
1.576
0.25
2.802
2.89
4.063
6.35
3.268

Table 2: The calculation of average value for pixels and in 25% of air opening flow in 25 mL of colour ingredient.

Pixels
Turbidity factor () of pixel
Pixel’s value ()
Total
Number of pixels
Average value ()

+ + + + + + + + + + + + + + + + + + + + + + + +
0.1107 + 1.236 + 0.244 + 0.2295 + 0.027 + 0.00016 + 0.0049 + 0.138 + 0.231 + 0.1467 + 0.133 + 0.0534 + 0.034 + 0.0615 + 0.0205 + 0.035 + 0.052 + 0.025 + 0.3656 + 0.0706 + 0.0013 + 0.018 + 0.4843 + 0.58 + 0.697
5.0000
25
0.2000

+ + + + + + + + + + + + + + + + + + + + +
1.546 + 0.366 + 0.2296 + 0.011 + 0.5009 + 0.5961 + 0.076 + 0.0067 + 0.00066 + 0.0055 + 0.00011 + 0.044 + 0.0083 + 0.0143 + 0.07 + 0.0519 + 0.0706 + 0.164 + 0.0896 + 0.6053 + 0.5798 + 0.7000
5.5022
22
0.2501

Figure 8: The analysis of transmitting light to receiver RX1.

5. Result and Discussion

This section presents the concentration profiles of gas bubbles flow in the vertical pipe filled with contaminated water. Five types of bubble flowing conditions are selected which are a single bubble flow, double bubble flow, and 25%, 50%, and 100% of air opening flow. Table 3 shows the gas bubble condition and concentration profile for a single bubble and double bubbles flow in contaminated water. The sample of contaminated water is prepared by adding 25 mL of green food colouring to the pipe which contains 3 liters of pure water filled in the pipe. The limitation of this research is that the experiment cannot be performed in volume higher than 25 mL such as 30 mL or 35 mL of green contaminated water. This is due to the fact that the ICA cannot separate the source signals properly since the water sample is too dark.

Table 3: The gas bubble condition and concentration profile for a single bubble and double bubble in contaminated water.

Gas bubble flow
Gas bubble condition
Concentration profile

A single bubble

Double bubbles

5.1. 25% of Air Opening

The experiments were conducted by opening up a quarter of the valve to allow the 25% volume of gas entering upwards in the vertical pipe as shown in Figure 9. The transmitter circuit is switched on, before the LabVIEW software is run to apply the ICA algorithm. The value is computed using the Matlab software by entering the inverse matrix acquired by the ICA process.

Figure 9: The condition of 25% of air opening.

Table 4 shows the result of value for 25% of air opening flow. The location of bubbles is estimated by analysing the value of . By referring to the mathematical model that has been discussed in considering the presence of gas bubbles, the value of receiver voltage is getting high compared to in homogeneous phase flow. The high value of is related to which varied when gas bubbles exist. The properties of are the outcome for matrix and have influenced the calculation of parameter. The range of values of is considered from 0.1 and above. Therefore, the bubbles are supposed to be located in the pixels which contain at least 0.1. The result in Table 4 is transformed to concentration profile as shown in Figure 10 in order to visualize bubble location. The pink colour in the concentration profile indicates that some bubbles are inclined to rise towards to the border of the pipe although the bubble injector is installed at the centre of the pipe base. The illustration for the bubble’s movements can be seen in Figure 11. Other than that, the high value of pixels located near the pipe’s periphery is probably due to the existence of receivers around the pipe’s periphery. Hence, most of the light is likely to cross the pixels around the pipe’s periphery. Although the average value is calculated, the value at the pipe’s periphery is still higher compared to other locations.

Table 4: value for 25% of air opening flow.

0.000
0.000
0.000
0.000
0.000
0.166
0.190
0.181
0.125
0.159
0.200
0.250
0.589
0.000
0.000
0.000
0.000
0.000

0.000
0.000
0.000
0.194
0.480
0.129
0.147
0.156
0.124
0.131
0.131
0.107
0.132
0.270
0.446
0.000
0.000
0.000

0.000
0.000
0.093
0.211
0.155
0.111
0.098
0.079
0.067
0.088
0.092
0.067
0.116
0.130
0.162
0.303
0.000
0.000

0.000
0.211
0.211
0.140
0.106
0.106
0.100
0.098
0.068
0.082
0.071
0.067
0.073
0.075
0.052
0.110
0.110
0.000

0.000
0.229
0.159
0.098
0.096
0.105
0.099
0.069
0.076
0.073
0.072
0.089
0.059
0.088
0.059
0.068
0.069
0.000

0.154
0.140
0.128
0.095
0.077
0.102
0.082
0.089
0.089
0.095
0.076
0.059
0.075
0.052
0.041
0.045
0.049
0.015

0.124
0.121
0.095
0.110
0.111
0.077
0.091
0.090
0.089
0.074
0.071
0.063
0.083
0.061
0.054
0.044
0.043
0.062

0.140
0.079
0.097
0.099
0.099
0.084
0.079
0.062
0.130
0.059
0.093
0.049
0.070
0.058
0.047
0.039
0.048
0.057

0.121
0.070
0.068
0.091
0.074
0.092
0.094
0.069
0.077
0.061
0.074
0.056
0.064
0.052
0.048
0.044
0.029
0.024

0.194
0.077
0.069
0.071
0.077
0.098
0.033
0.127
0.017
0.093
0.056
0.061
0.070
0.054
0.034
0.041
0.035
0.035

0.157
0.070
0.051
0.073
0.089
0.040
0.083
0.058
0.097
0.066
0.074
0.045
0.072
0.045
0.038
0.032
0.032
0.052

0.099
0.065
0.062
0.074
0.049
0.078
0.075
0.066
0.076
0.081
0.060
0.055
0.068
0.048
0.047
0.028
0.026
0.055

0.064
0.067
0.054
0.053
0.062
0.069
0.072
0.083
0.072
0.066
0.087
0.065
0.071
0.047
0.050
0.034
0.031
0.030

0.000
0.076
0.068
0.045
0.068
0.064
0.077
0.053
0.105
0.079
0.058
0.081
0.079
0.043
0.067
0.041
0.050
0.000

0.000
0.069
0.094
0.072
0.061
0.041
0.065
0.062
0.045
0.104
0.049
0.074
0.088
0.061
0.052
0.039
0.027
0.000

0.000
0.000
0.074
0.081
0.068
0.067
0.054
0.048
0.054
0.056
0.085
0.064
0.055
0.087
0.073
0.013
0.000
0.000

0.000
0.000
0.000
0.083
0.147
0.098
0.052
0.070
0.055
0.036
0.088
0.092
0.058
0.132
0.056
0.000
0.000
0.000

0.000
0.000
0.000
0.000
0.000
0.069
0.075
0.140
0.082
0.049
0.075
0.148
0.088
0.000
0.000
0.000
0.000
0.000

Figure 10: The concentration profile for 25% of air opening.

Figure 11: Bubbles movement upwards from the bottom.

5.2. 50% of Air Opening

Half of the valve is opened in this experiment as in Figure 12(a) to ensure that more bubbles are inserted into the pipe than the first experiment. The concentration profile for the medium bubble flow is shown in Figure 12(b). The result of value is higher compared to the 25% of air opening where the range for bubble is from 0.12 to 0.5. It is due to the light’s energy experiencing more attenuation when the bubbles interrupted the light path as modelled in (4) to (6). Moreover, the bubbles also became close to each other when more bubbles went upward since the space around the pipe is limited and compact. Thus, the collision between bubbles occurred and generated the additional energy which resulted in the value of increasing.

(a)

(b)

Figure 12: (a) The condition of 50% of air opening. (b) The concentration profile for 50% of air opening.

5.3. 100% of Air Opening

In this experiment, the valve is fully open and the bubble flow condition is shown in Figure 13(a). The range of bubble is from 0.2 to 0.5, and this represents an increase compared to medium bubble flow. The location of bubble can be seen in Figure 13(b) which shows that most of the area in the cross section of the pipe is occupied by gas bubbles. Thus, all sensors in the system can detect the presence of bubble in high bubble flow condition. This condition also creates more collision between the bubbles than medium bubble flow since the free space in the pipe has been reduced by the addition of gas bubbles. More collisions of bubbles have influenced the reading in terms of the increasing value of .

(a)

(b)

Figure 13: (a) The condition of 100% of air opening. (b) The concentration profile for 100% of air opening.

6. Conclusion

In this paper, an optical tomography system using the ICA method has been presented to estimate the distribution of bubbles in a vertical pipeline. The system proved its ability to obtain information on the gas bubbles concentration profile in the pipe without disturbing the process flow. The concentration profiles for the bubble flow quantities are investigated under five different flows conditions: a single bubble, double bubbles, 25% of air opening, 50% of air opening, and 100% of air opening flow rates. The ICA method has a great potential in improving the process efficiency by visualizing the concentration profiles inside the pipelines.

Appendix

See Tables 5–9.

Table 5: value of 25 mL.

0.000
0.000
0.000
0.000
0.000
0.264
0.316
0.203
0.147
0.220
0.211
0.173
0.204
0.000
0.000
0.000
0.000
0.000

0.000
0.000
0.000
0.083
0.291
0.161
0.154
0.106
0.094
0.091
0.089
0.114
0.130
0.088
0.070
0.000
0.000
0.000

0.000
0.000
0.121
0.201
0.091
0.088
0.092
0.079
0.065
0.068
0.059
0.048
0.073
0.068
0.072
0.107
0.000
0.000

0.000
0.055
0.125
0.105
0.083
0.077
0.067
0.065
0.058
0.087
0.052
0.047
0.058
0.061
0.054
0.084
0.059
0.000

0.000
0.137
0.078
0.049
0.079
0.055
0.060
0.072
0.042
0.053
0.055
0.062
0.042
0.079
0.070
0.049
0.108
0.000

0.358
0.093
0.054
0.048
0.061
0.059
0.050
0.070
0.062
0.090
0.067
0.035
0.076
0.034
0.058
0.092
0.094
0.115

0.091
0.062
0.042
0.051
0.056
0.044
0.040
0.067
0.068
0.061
0.047
0.062
0.054
0.085
0.075
0.054
0.099
0.156

0.067
0.065
0.046
0.068
0.048
0.036
0.042
0.038
0.054
0.072
0.056
0.050
0.053
0.062
0.072
0.082
0.070
0.134

0.060
0.066
0.051
0.055
0.056
0.034
0.038
0.046
0.049
0.050
0.060
0.071
0.039
0.061
0.051
0.077
0.099
0.200

0.057
0.055
0.058
0.053
0.060
0.041
0.025
0.053
0.020
0.066
0.050
0.073
0.057
0.053
0.084
0.047
0.068
0.192

0.066
0.061
0.070
0.057
0.048
0.034
0.048
0.053
0.060
0.049
0.058
0.037
0.048
0.060
0.057
0.054
0.063
0.137

0.094
0.050
0.074
0.053
0.042
0.037
0.042
0.056
0.064
0.041
0.040
0.051
0.049
0.045
0.057
0.079
0.075
0.232

0.046
0.045
0.053
0.051
0.067
0.036
0.042
0.046
0.056
0.047
0.052
0.028
0.053
0.062
0.059
0.065
0.077
0.099

0.000
0.046
0.059
0.062
0.041
0.051
0.026
0.040
0.057
0.041
0.058
0.060
0.053
0.054
0.064
0.049
0.099
0.000

0.000
0.048
0.053
0.065
0.050
0.060
0.029
0.045
0.044
0.037
0.038
0.065
0.037
0.057
0.057
0.084
0.089
0.000

0.000
0.000
0.047
0.089
0.049
0.051
0.054
0.046
0.046
0.035
0.049
0.052
0.056
0.040
0.054
0.039
0.000
0.000

0.000
0.000
0.000
0.082
0.090
0.064
0.042
0.049
0.062
0.044
0.048
0.053
0.061
0.061
0.013
0.000
0.000
0.000

0.000
0.000
0.000
0.000
0.000
0.025
0.070
0.132
0.129
0.073
0.077
0.049
0.024
0.000
0.000
0.000
0.000
0.000

Table 6: value of a single bubble flow.

0.000
0.000
0.000
0.000
0.000
0.095
0.193
0.123
0.115
0.092
0.087
0.081
0.206
0.000
0.000
0.000
0.000
0.000

0.000
0.000
0.000
0.090
2.730
0.152
0.116
0.096
0.069
0.089
0.082
0.056
0.064
0.250
0.630
0.000
0.000
0.000

0.000
0.000
0.080
0.116
0.149
0.135
0.095
0.067
0.077
0.049
0.069
0.073
0.098
0.141
0.211
0.140
0.000
0.000

0.000
0.122
0.178
0.081
0.129
0.111
0.099
0.088
0.085
0.083
0.034
0.057
0.105
0.077
0.109
0.123
0.090
0.000

0.000
0.125
0.114
0.111
0.149
0.115
0.111
0.077
0.048
0.058
0.071
0.130
0.077
0.059
0.064
0.118
0.099
0.000

0.172
0.099
0.104
0.082
0.114
0.109
0.095
0.100
0.087
0.330
0.116
0.074
0.134
0.083
0.090
0.089
0.131
0.048

0.092
0.087
0.108
0.078
0.090
0.102
0.080
0.060
0.103
0.077
0.047
0.115
0.104
0.089
0.114
0.124
0.101
0.221

0.136
0.059
0.094
0.142
0.052
0.086
0.078
0.102
0.054
0.095
0.067
0.099
0.067
0.076
0.056
0.078
0.167
0.253

0.217
0.087
0.079
0.133
0.077
0.121
0.066
0.073
0.088
0.064
0.056
0.109
0.125
0.120
0.131
0.144
0.122
0.181

0.249
0.118
0.099
0.061
0.141
0.077
0.038
0.077
0.043
0.114
0.138
0.132
0.128
0.065
0.092
0.104
0.129
0.223

0.229
0.101
0.108
0.114
0.120
0.075
0.126
0.121
0.075
0.055
0.116
0.088
0.106
0.085
0.045
0.129
0.109
0.219

0.560
0.175
0.164
0.139
0.102
0.138
0.096
0.078
0.126
0.122
0.110
0.062
0.068
0.087
0.110
0.082
0.113
0.548

0.411
0.218
0.122
0.104
0.139
0.146
0.161
0.112
0.070
0.090
0.053
0.056
0.057
0.084
0.103
0.100
0.115
0.268

0.000
0.268
0.223
0.092
0.106
0.120
0.127
0.049
0.113
0.145
0.080
0.079
0.099
0.146
0.117
0.153
0.274
0.000

0.000
0.163
0.220
0.252
0.111
0.102
0.114
0.101
0.065
0.110
0.064
0.094
0.118
0.134
0.124
0.168
0.152
0.000

0.000
0.000
0.139
0.176
0.175
0.134
0.101
0.125
0.085
0.085
0.096
0.123
0.121
0.120
0.159
0.049
0.000
0.000

0.000
0.000
0.000
0.140
0.198
0.179
0.157
0.173
0.158
0.107
0.161
0.150
0.154
0.183
0.050
0.000
0.000
0.000

0.000
0.000
0.000
0.000
0.000
0.115
0.207
0.430
0.721
0.526
0.304
0.170
0.058
0.000
0.000
0.000
0.000
0.000

Table 7: value of double bubble flow.

0.000
0.000
0.000
0.000
0.000
0.196
0.508
0.109
0.070
0.119
0.089
0.555
0.158
0.000
0.000
0.000
0.000
0.000

0.000
0.000
0.000
0.091
0.195
0.122
0.126
0.087
0.061
0.067
0.069
0.059
0.075
0.090
0.086
0.000
0.000
0.000

0.000
0.000
0.134
0.221
0.096
0.106
0.078
0.062
0.072
0.045
0.063
0.067
0.055
0.048
0.055
0.030
0.000
0.000

0.000
0.139
0.191
0.122
0.075
0.081
0.061
0.058
0.050
0.079
0.044
0.056
0.046
0.063
0.050
0.056
0.017
0.000

0.000
0.218
0.135
0.085
0.068
0.089
0.060
0.058
0.042
0.048
0.044
0.053
0.025
0.068
0.038
0.048
0.058
0.000

0.185
0.161
0.124
0.071
0.071
0.089
0.044
0.067
0.065
0.041
0.078
0.044
0.052
0.035
0.049
0.059
0.065
0.044

0.121
0.153
0.081
0.065
0.095
0.077
0.044
0.038
0.065
0.028
0.041
0.065
0.039
0.046
0.052
0.039
0.068
0.091

0.154
0.108
0.090
0.095
0.068
0.059
0.062
0.050
0.155
0.037
0.040
0.039
0.051
0.050
0.050
0.045
0.046
0.076

0.140
0.111
0.102
0.082
0.075
0.041
0.150
0.043
0.053
0.039
0.043
0.040
0.069
0.058
0.046
0.040
0.052
0.085

0.204
0.115
0.089
0.060
0.084
0.049
0.037
0.073
0.021
0.055
0.052
0.073
0.067
0.078
0.036
0.043
0.057
0.105

0.213
0.135
0.073
0.066
0.085
0.045
0.042
0.050
0.058
0.074
0.061
0.046
0.070
0.081
0.071
0.044
0.041
0.072

0.284
0.117
0.114
0.108
0.063
0.066
0.062
0.080
0.065
0.094
0.040
0.064
0.050
0.055
0.080
0.068
0.058
0.118

0.222
0.149
0.075
0.083
0.094
0.046
0.063
0.056
0.041
0.057
0.093
0.063
0.079
0.084
0.073
0.122
0.104
0.091

0.000
0.206
0.131
0.098
0.070
0.078
0.056
0.064
0.057
0.055
0.075
0.057
0.055
0.078
0.084
0.078
0.115
0.000

0.000
0.138
0.118
0.145
0.089
0.100
0.050
0.058
0.047
0.040
0.057
0.062
0.054
0.057
0.070
0.098
0.050
0.000

0.000
0.000
0.084
0.131
0.104
0.082
0.061
0.063
0.071
0.068
0.062
0.068
0.054
0.076
0.121
0.049
0.000
0.000

0.000
0.000
0.000
0.091
0.094
0.123
0.071
0.056
0.044
0.062
0.071
0.076
0.108
0.104
0.039
0.000
0.000
0.000

0.000
0.000
0.000
0.000
0.000
0.082
0.131
0.094
0.113
0.096
0.149
0.129
0.028
0.000
0.000
0.000
0.000
0.000

Table 8: value for 50% of air opening flow.

0.000
0.000
0.000
0.000
0.000
0.062
0.224
0.497
0.573
0.770
0.875
1.261
3.137
0.000
0.000
0.000
0.000
0.000

0.000
0.000
0.000
0.233
0.403
0.266
0.365
0.387
0.249
0.345
0.436
0.296
0.410
0.752
0.892
0.000
0.000
0.000

0.000
0.000
0.340
0.366
0.210
0.274
0.284
0.297
0.241
0.266
0.201
0.176
0.200
0.255
0.419
0.931
0.000
0.000

0.000
0.139
0.360
0.200
0.153
0.276
0.257
0.220
0.190
0.241
0.184
0.169
0.188
0.187
0.145
0.340
0.648
0.000

0.000
0.190
0.295
0.193
0.286
0.153
0.312
0.166
0.145
0.187
0.091
0.134
0.146
0.143
0.152
0.153
0.301
0.000

0.330
0.355
0.245
0.169
0.227
0.231
0.181
0.218
0.185
0.164
0.112
0.165
0.118
0.090
0.134
0.131
0.064
0.032

0.406
0.240
0.168
0.192
0.228
0.157
0.186
0.118
0.166
0.116
0.118
0.106
0.079
0.118
0.098
0.058
0.088
0.059

0.349
0.277
0.258
0.187
0.159
0.185
0.155
0.160
0.153
0.088
0.115
0.069
0.090
0.106
0.043
0.065
0.066
0.052

0.455
0.301
0.154
0.212
0.174
0.127
0.139
0.100
0.087
0.081
0.072
0.094
0.085
0.072
0.045
0.091
0.045
0.076

0.704
0.236
0.200
0.130
0.108
0.153
0.072
0.143
0.035
0.110
0.073
0.092
0.089
0.075
0.048
0.079
0.051
0.059

0.670
0.171
0.206
0.160
0.151
0.107
0.111
0.111
0.114
0.066
0.107
0.044
0.069
0.040
0.079
0.069
0.047
0.028

0.595
0.232
0.165
0.100
0.107
0.089
0.138
0.087
0.101
0.075
0.076
0.071
0.052
0.061
0.092
0.048
0.049
0.048

0.341
0.186
0.157
0.163
0.105
0.108
0.134
0.118
0.094
0.082
0.077
0.071
0.082
0.070
0.105
0.051
0.067
0.033

0.000
0.240
0.219
0.158
0.102
0.122
0.104
0.099
0.096
0.096
0.116
0.052
0.073
0.085
0.121
0.098
0.115
0.000

0.000
0.252
0.223
0.196
0.136
0.116
0.088
0.094
0.076
0.100
0.064
0.104
0.101
0.091
0.070
0.043
0.027
0.000

0.000
0.000
0.077
0.232
0.217
0.134
0.103
0.075
0.090
0.075
0.075
0.059
0.108
0.099
0.110
0.031
0.000
0.000

0.000
0.000
0.000
0.151
0.564
0.188
0.105
0.092
0.075
0.047
0.039
0.035
0.068
0.129
0.069
0.000
0.000
0.000

0.000
0.000
0.000
0.000
0.000
0.130
0.097
0.226
0.104
0.053
0.045
0.043
0.048
0.000
0.000
0.000
0.000
0.000

Table 9: value for 100% of air opening flow.

0.000
0.000
0.000
0.000
0.000
0.055
0.141
0.188
0.258
0.506
0.450
0.490
0.682
0.000
0.000
0.000
0.000
0.000

0.000
0.000
0.000
0.080
0.156
0.136
0.170
0.238
0.180
0.207
0.204
0.253
0.453
0.253
0.315
0.000
0.000
0.000

0.000
0.000
0.059
0.144
0.161
0.212
0.177
0.127
0.201
0.195
0.186
0.193
0.304
0.290
0.406
1.011
0.000
0.000

0.000
0.344
0.207
0.204
0.151
0.161
0.160
0.130
0.127
0.094
0.133
0.180
0.175
0.267
0.238
0.413
0.246
0.000

0.000
0.508
0.171
0.140
0.112
0.097
0.105
0.107
0.099
0.066
0.120
0.180
0.078
0.274
0.217
0.299
0.287
0.000

0.501
0.318
0.245
0.213
0.198
0.238
0.198
0.165
0.146
0.175
0.183
0.158
0.211
0.166
0.237
0.354
0.383
0.362

0.207
0.158
0.133
0.154
0.193
0.135
0.170
0.179
0.193
0.115
0.104
0.126
0.139
0.185
0.230
0.300
0.432
0.685

0.246
0.136
0.184
0.179
0.146
0.149
0.158
0.121
0.065
0.121
0.212
0.227
0.190
0.252
0.188
0.183
0.323
0.478

0.240
0.133
0.141
0.168
0.126
0.144
0.148
0.134
0.182
0.156
0.162
0.142
0.201
0.176
0.251
0.274
0.342
0.756

0.282
0.165
0.129
0.172
0.162
0.154
0.102
0.132
0.120
0.197
0.244
0.215
0.177
0.130
0.212
0.323
0.345
1.076

0.336
0.200
0.186
0.154
0.133
0.101
0.138
0.182
0.159
0.126
0.187
0.153
0.241
0.238
0.176
0.226
0.336
0.704

0.377
0.169
0.170
0.146
0.114
0.190
0.143
0.175
0.170
0.217
0.172
0.176
0.144
0.155
0.247
0.186
0.385
1.158

0.309
0.190
0.157
0.152
0.162
0.119
0.168
0.139
0.099
0.144
0.202
0.191
0.172
0.181
0.238
0.385
0.339
0.513

0.000
0.327
0.194
0.146
0.140
0.170
0.177
0.165
0.202
0.125
0.154
0.155
0.152
0.280
0.337
0.193
0.443
0.000

0.000
0.060
0.173
0.160
0.158
0.201
0.118
0.184
0.138
0.156
0.151
0.124
0.178
0.301
0.222
0.331
0.207
0.000

0.000
0.000
0.131
0.266
0.154
0.134
0.126
0.148
0.127
0.139
0.165
0.116
0.270
0.246
0.343
0.174
0.000
0.000

0.000
0.000
0.000
0.146
0.184
0.168
0.174
0.158
0.183
0.151
0.155
0.234
0.392
0.287
0.088
0.000
0.000
0.000

0.000
0.000
0.000
0.000
0.000
0.112
0.190
0.201
0.256
0.342
0.636
0.570
0.190
0.000
0.000
0.000
0.000
0.000

Competing Interests

The authors declare that there are no competing interests regarding the publication of this paper.

Acknowledgments

The authors wish to acknowledge the assistance of the Ministry of Higher Education, Malaysia, and Universiti Teknologi Malaysia under the GUP Research Vote 05H67 for providing the funds and resources in carrying out this research.

References

C. Yan, J. Zhong, Y. Liao, S. Lai, M. Zhang, and D. Gao, “Design of an applied optical fiber process tomography system,” Sensors and Actuators B: Chemical, vol. 104, no. 2, pp. 324–331, 2005.

D. L. George, J. R. Torczynski, K. A. Shollenberger, and T. J. O'Hern, “Validation of electrical-impedance tomography for measurements of material distribution in two-phase flows,” International Journal of Multiphase Flow, vol. 26, no. 4, pp. 549–581, 2000.

T. Dyakowski, L. F. C. Jeanmeure, and A. J. Jaworski, “Applications of electrical tomography for gas-solids and liquid-solids flows—a review,” Powder Technology, vol. 112, no. 3, pp. 174–192, 2000.

R. Abdul Rahim, C. K. Thiam, and M. H. F. Rahiman, “An optical tomography system using a digital signal processor,” Sensors, vol. 8, no. 4, pp. 2082–2103, 2008.

M. Idroas, R. Abdul Rahim, M. H. Fazalul Rahiman, M. N. Ibrahim, and R. G. Green, “Design and development of a CCD based optical tomography measuring system for particle sizing identification,” Measurement, vol. 44, no. 6, pp. 1096–1107, 2011.

M. R. Rzasa and A. Plaskowski, “Application of optical tomography for measurements of aeration parameters in large water tanks,” Measurement Science and Technology, vol. 14, no. 2, pp. 199–204, 2003.

S. Ibrahim, M. A. M. Yunus, R. G. Green, and K. Dutton, “Concentration measurements of bubbles in a water column using an optical tomography system,” ISA Transactions, vol. 51, no. 6, pp. 821–826, 2012.

R. M. Detsch and R. N. Sharma, “The critical angle for gas bubble entrainment by plunging liquid jets,” The Chemical Engineering Journal, vol. 44, no. 3, pp. 157–166, 1990.

M. R. Rzasa, “The measuring method for tests of horizontal two-phase gas-liquid flows, using optical and capacitance tomography,” Nuclear Engineering and Design, vol. 239, no. 4, pp. 699–707, 2009.

H. Jin, Y. Lian, S. Yang, G. He, and Z. Guo, “The parameters measurement of air–water two phase flow using the electrical resistance tomography (ERT) technique in a bubble column,” Flow Measurement and Instrumentation, vol. 31, pp. 55–60, 2013.

N. M. N. Ayob, M. H. F. Rahiman, Z. Zakaria, S. Yaacob, and R. A. Rahim, “Detection of small gas bubble using ultrasonic transmission-mode tomography system,” in Proceedings of the IEEE Symposium on Industrial Electronics & Applications (ISIEA '11), pp. 165–170, IEEE, Penang, Malaysia, October 2010.

A. R. Daniels, Dual modality tomography for the monitoring of constituent volumes in multi-component flows [Ph.D. thesis], Sheffield Hallam University, 1996.

G. R. Naik and D. K. Kumar, “An overview of independent component analysis and its applications,” Informatica, vol. 35, pp. 63–81, 2011.

Y. Xu, H. Wang, Z. Cui, F. Dong, and Y. Yan, “Separation of gas-liquid two-phase flow through independent component analysis,” IEEE Transactions on Instrumentation and Measurement, vol. 59, no. 5, pp. 1294–1302, 2010.

C. Jutten and J. Herault, “Blind separation of sources, part I: an adaptive algorithm based on neuromimetic architecture,” Signal Processing, vol. 24, no. 1, pp. 1–10, 1991.

J. Chen and X. Z. Wang, “A new approach to near-infrared spectral data analysis using independent component analysis,” Journal of Chemical Information and Computer Sciences, vol. 41, no. 4, pp. 992–1001, 2001.

J. P. Huang and J. Mar, “Combined ICA and FCA schemes for a hierarchical network,” Wireless Personal Communications, vol. 28, no. 1, pp. 35–58, 2004.

S. Fiori, “Overview of independent component analysis technique with an application to synthetic aperture radar (SAR) imagery processing,” Neural Networks, vol. 16, no. 3-4, pp. 453–467, 2003.

C. J. James and C. W. Hesse, “Independent component analysis for biomedical signals (topical review),” Physiological Measurement, vol. 26, no. 1, pp. R15–R39, 2005.

G. Wang, Q. Ding, and Z. Hou, “Independent component analysis and its applications in signal processing for analytical chemistry,” TrAC—Trends in Analytical Chemistry, vol. 27, no. 4, pp. 368–376, 2008.

A. Hyvärinen, “Fast and robust fixed-point algorithms for independent component analysis,” IEEE Transactions on Neural Networks, vol. 10, no. 3, pp. 626–634, 1999.

A. Hyvärinen and E. Oja, “Independent component analysis: algorithms and applications,” Neural Networks, vol. 13, no. 4-5, pp. 411–430, 2000.

A. J. Bell and T. J. Sejnowski, “An information-maximization approach to blind separation and blind deconvolution,” Neural Computation, vol. 7, no. 6, pp. 1129–1159, 1995.

J. H. Garvey, Independent component analysis by entropy maximization (Infomax) [M.S. thesis], Naval Postgraduate School, Monterey, Calif, USA, 2007.