Outline
Comptes Rendus

Estimation of the thermal conductivity λ(T,P) of ionic liquids using a neural network optimized with genetic algorithms
Comptes Rendus. Chimie, Volume 19 (2016) no. 3, pp. 333-341.

Abstract

In this study, an artificial neural network was optimized using a genetic algorithm in order to estimate the thermal conductivity of ionic liquids at different temperatures and pressures. Experimental thermal conductivity data of 41 ionic liquids (400 experimental data points) in the range from 0.10 to 0.22 W m−1 K−1 were used to obtain the proposed method for the temperature range of 273–390 K and the pressure range of 100–20,000 kPa. In addition, the molecular mass M and structure of molecules, represented by the number of well-defined groups forming the molecule, were provided as input parameters in order to characterize the different molecules of ionic liquids. A heterogeneous set of ionic liquids includes cations such as imidazolium, ammonium, phosphonium, pyrrolidinium, and pyridinium. It also includes anions such as halides, sulfonates, tosylates, imides, borates, phosphates, acetates, and amino acids. The whole dataset was divided into a training set with 300 experimental data points and a prediction set with 100 experimental data points. Several architectures were studied, and the optimum weights for the network were determined. The results showed that the proposed method to estimate the thermal conductivity of ionic liquids at different temperatures and pressures presented a good accuracy with lower deviations such as AARD less than 0.91% and R2 of 0.9969 for the training set, and AARD less than 0.84% with R2 of 0.9963 for the prediction set.

Supplementary Materials:
Supplementary material for this article is supplied as a separate file:

Metadata
Received:
Accepted:
Published online:
DOI: 10.1016/j.crci.2015.09.010
Keywords: Ionic liquids, Thermal conductivity, Structural groups, Artificial neural network, Genetic algorithm, Transport properties

Juan A. Lazzús 1

1 Departamento de Física, Universidad de La Serena, Casilla 554, La Serena, Chile
@article{CRCHIM_2016__19_3_333_0,
     author = {Juan A. Lazz\'us},
     title = {Estimation of the thermal conductivity {\ensuremath{\lambda}(\protect\emph{T},\protect\emph{P})} of ionic liquids using a neural network optimized with genetic algorithms},
     journal = {Comptes Rendus. Chimie},
     pages = {333--341},
     publisher = {Elsevier},
     volume = {19},
     number = {3},
     year = {2016},
     doi = {10.1016/j.crci.2015.09.010},
     language = {en},
}
TY  - JOUR
AU  - Juan A. Lazzús
TI  - Estimation of the thermal conductivity λ(T,P) of ionic liquids using a neural network optimized with genetic algorithms
JO  - Comptes Rendus. Chimie
PY  - 2016
SP  - 333
EP  - 341
VL  - 19
IS  - 3
PB  - Elsevier
DO  - 10.1016/j.crci.2015.09.010
LA  - en
ID  - CRCHIM_2016__19_3_333_0
ER  - 
%0 Journal Article
%A Juan A. Lazzús
%T Estimation of the thermal conductivity λ(T,P) of ionic liquids using a neural network optimized with genetic algorithms
%J Comptes Rendus. Chimie
%D 2016
%P 333-341
%V 19
%N 3
%I Elsevier
%R 10.1016/j.crci.2015.09.010
%G en
%F CRCHIM_2016__19_3_333_0
Juan A. Lazzús. Estimation of the thermal conductivity λ(T,P) of ionic liquids using a neural network optimized with genetic algorithms. Comptes Rendus. Chimie, Volume 19 (2016) no. 3, pp. 333-341. doi : 10.1016/j.crci.2015.09.010. https://comptes-rendus.academie-sciences.fr/chimie/articles/10.1016/j.crci.2015.09.010/

Version originale du texte intégral

1 Introduction

Ionic liquids (ILs) are salts with a very low melting temperature and typically consist of a large organic cation and an inorganic polyatomic anion [1]. In the last few years, ILs have drawn the attention of the scientific community, and hundreds of studies that involve different aspects of ILs have been published in the scientific literature. From the scientific and industrial points of view, a fundamental understanding of the physico-chemical properties of ILs is needed before their application to several processes. For instance, knowledge of some basic properties can be useful for fluid property estimation, thermodynamic property calculations, and phase equilibrium [2].

Regarding ILs, thermal conductivity is important to be known before obtaining the heat transfer coefficient of fluids that is essential for the design of heat transfer fluid and equipment [3]. There is a large variety of analytical expressions that allow the correlation and prediction of thermal conductivity. In the case of ILs, these expressions are usually based on the use of adjustable parameters of each fluid (correlations) [3–6], or based on the group contribution methods (GCMs). In this context, Gardas and Countinho [7], and Albert and Müller [8] proposed a GCM for predicting the thermal conductivity of imidazolium, pyrrolidinium, and phosphonium ILs with a deviation below the mean. Wu et al. [9] proposed another GCM for predicting the thermal conductivities of ILs, in combination with the Valderrama's group contribution method for critical properties of ILs [10]. However, to the best of the author's knowledge there is no application of GCM for the estimation of the thermal conductivity of ILs as a function of temperature and pressure, as presented here, and certainly there is not any publications on the estimation of thermal conductivity for a heterogeneous set of ILs using a hybrid artificial neural network.

The relationship between the physical and thermodynamic properties is highly nonlinear, and consequently an artificial neural network (ANN) can be a suitable alternative to modelling underlying thermodynamic properties. An ANN is an especially efficient algorithm for approximating any function with a finite number of discontinuities by learning the relationships between the input and output vectors [11]. Therefore, an ANN is an appropriate technique for modelling the nonlinear behaviour of thermophysical properties [2]. In this context, Hezave et al. [12] presented a neural network model for predicting the thermal conductivity of ionic liquids using temperature T, pressure P, molecular weight MW, and melting point temperature Tm of ILs as input parameters.

In this work, thermal conductivities of several ILs at different temperatures and pressures were correlated and estimated using an ANN optimized with a genetic algorithm (GA) [13] in order to update the weights of the network.

2 Neural network and genetic algorithm (ANN+GA)

The most successful and frequently used type of neural network, a multilayer feed-forward neural network with a back-propagation learning algorithm (gradient descent error), was implemented in the study. The ANN consisted of one input layer with N inputs, one hidden layer with q units and one output layer with n outputs. The output of this model can be expressed as [11]:

yn=fn(j=1qWnjfj(i=1NWjixi+Wj0)+Wn0)(1)
where Wnj are the weights between unit j and unit n of input and hidden layers and Wji are the weights between the hidden layer and an output. The activation functions fn(x) and fj(x) are linear or nonlinear. One hidden layer with fj(x) was used as a tangent hyperbolic nonlinear activation function and fn(x) as the linear function in the output layer. For a given set of N inputs, the root mean square error (RMSE) is defined by:
RMSE=s=1N(ysexpyscalc)2N2(2)
where yexp denotes the experimental input and ycalc is the neural network output. This network was trained to minimize RMSE, replacing the gradient descent error with a genetic algorithm (GA).

The GA was first developed by Holland [13] and it was based on the mechanics of natural selection in biological systems. It uses a structure to utilize genetic information for finding new search directions. Major genetic operators that reflect the nature of the evolutionary process are reproduction, crossover and mutation [14].

The GA maintains a population of individuals, whose characteristics are encoded in a fixed-length bit string, modelling the biological genotype [15]. The way these bits represent the phenotype (ontogeny) is at the programmer's direction. As in nature, the genetic material is swapped between the individuals and mutated to produce offspring, with the corresponding changes in phenotypic performance. The crossover operator is an analogue for the recombination of the genetic material as observed in reproduction. Crossover involves splitting the genomic two-parent bit-strings at a given number of locations and then splicing complementary sections of each parent's bit-string in order to form the genotype of the new individual. Crossover occurs with a random probability. The mutation operator simulates natural mutation of DNA. This involves flipping bits in the string in a stochastic manner. Mutation should be fairly infrequent and should be applied following the crossover [14].

The most significant differences among GAs are: i) only the objective function and the corresponding fitness levels influence search directions; ii) GAs use probabilistic transition rules, not deterministic ones; and iii) GAs work on an encoding environment of the parameter set rather than the parameter set itself [15].

3 Database and training

Thermal conductivity at different temperatures and pressures of 41 ionic liquids (400 experimental data points in total) was collected. In this dataset, λ(T,P) properties cover wide ranges: 273–390 K for temperature, 100–20000 kappa for pressure, and 0.10–0.22 W m−1 K−1 for thermal conductivity. This dataset includes cations such as imidazolium, ammonium, phosphonium, pyrrolidinium, and pyridinium. It also includes anions such as halides, sulfonates, tosylates, imides, borates, phosphates, acetates, and amino acids. Fig. 1 shows a general picture of the range of thermal conductivities and ILs considered. These values are of special importance in order to verify that an acceptable range of λ was covered in this study. All data used were chosen from specific databases [16], and corresponded to those claimed by the authors as being experimentally determined [5,6,17–24]. Data available in the literature from theoretical methods, correlations, and extrapolations of any kind were not considered, as well as data whose accuracy was not guaranteed by the authors themselves for any reasons whatsoever (presence of impurities, fluid instability, or problems with the equipment) [25].

Fig. 1

Thermal conductivity as a function of the molecular mass of all ionic liquids used in this study. (a) Total mass distribution, (b) cation distribution, and (c) anion distribution.

Additionally, molecular mass M and molecular structures, represented by the number of well-defined groups forming the molecule, were provided as input parameters in order to characterize the different molecules of ILs. A leave-25%-out cross-validation method was used to estimate the predictive capabilities of the proposed method for entrances such as temperature T, pressure P, molecular mass M, and 28 structural groups. Thus, the whole dataset was divided into a training set with 300 experimental data points, and a prediction set with 100 experimental data points. Table S1 shows the values of experimental thermal conductivity considered in this study (see Supplementary data). Training and prediction sets were selected considering that molecules were decomposed into fragments and all fragments were present with adequate frequency in the training dataset for structural groups. The value associated with the structural group was defined as 0 when the group does not appear in the substance, and as k, when the group appeared k-times in the substance [26]. Table 1 shows all 31 input parameters considered in this study, and contains minimum and maximum input values, and the number of input occurrences in the datasets.

Table 1

Input considered in the proposed method.

No.InputXminXmaxNo. occurrence
Training setPrediction setTotal set
1T273.15390.00300100400
2P10020000300100400
3M66.43515.13300100400
4
Imidazolium(+)
0119942241
5
Pyridinium(+)
0118927
6
Pyrrolidinium(+)
017916
7
Ammonium(+)
0130737
8
Phosphonium(+)
01463379
9–H(+)0121741258
10CH3(+)14300100400
11–CH2(+)128300100400
12CH3(−)0474579
13–CH2(−)010561470
14>CH–(−) [–CH–](−)02511263
15>C<(−) [>C–](−)0123528
16–COO–(−)01909
17–COOH(−)01421254
18–OH(−)0130737
19–O–(−) [–O](−)02651479
20O(−)01909
21–CN(−)039918
22>N–(−)016449113
23NH2(−)02422668
24–SO2(−)0211154165
25–SH(−)01707
26CF3(−)039440134
27–CF2(−)0316016
28–F(−)0610425129
29–Cl(−)01707
30–P(−)01591675
31–B(−)0163972

Subsequently, these input parameters were normalized using the following equation [11]:

xi=(XiXimin)2XimaxXimin1(3)
where Xi is the input variable i, and Ximin and Ximax are the smallest and largest values of the data, respectively. Thus, the net inputs (N) are calculated for the hidden neurons coming from the input neurons as shown in Eq. (1), and then a GA is used to obtain the optimum weights for the feed-forward ANN.

The step-by-step for optimum weights calculation using GA is described as follows:

  • 1) Initial weights in the ANN are randomly generated (initial population). M-chromosomes are generated randomly in order to serve as the initial population. Each chromosome represents all initial weights and biases in the ANN, which are optimized by the GA.
  • 2) Chromosome fitness is evaluated based on the ANN performance. In this case, the fitness function is defined as the root mean square error (RMSE).
  • 3) Fitness function value of each individual in the population is evaluated. The lower the RMSE, the higher the probability for chromosomes to be passed down to the next generation. The best individual chromosomes are selected for mating. The selection is repeated until the number of individuals in the mating pool is the same as the number of individuals in the population [26].
  • 4) Two individuals are selected randomly from the mating pool as the parent, and random gene selection is used as a crossover point [26]. The transfer of the parent genes to the right of the crossover point generates two new individuals, ending of the single-point crossover [27].
  • 5) A mutation operator is applied in order to maintain diversity within the population. Since the initial weights of the ANN could take any values between zero and one, the mutation is conducted by switching random genes. The approximate optimum solutions can be found quickly in order to set up the mutation rate as a parameter to control mutation probability. The mutation rate is a real number between 0 and 1 [28].
  • 6) Finally, the best chromosome with the minimum RMSE is chosen and set as the initial weights for the next period in the ANN. The process is repeated until the convergence toward a population of individuals encoding the set of weights that keep the error at a desired level occurs.

Fig. 2 shows a block diagram of the ANN+GA developed in this study. The full methodology was programmed in C++. In a GA, the number of individuals, the crossover operator, the crossover probability, the mutation operator, and the mutation probability summarize the parameters for synchronizing their application in a given problem. An exhaustive trial-and-error procedure was applied for tuning GA parameters. Table 2 shows the selected parameters for the ANN+GA algorithm.

Fig. 2

Flow diagram for training of the ANN using GA.

Table 2

Parameters used in the hybrid ANN+GA technique.

SectionParameterValue
ANNNN-typefeed-forward
No. hidden layers1
Transfer function (hidden)tansig
Transfer function (output)linear
No. iterations1500
Normalization range[−1, 1]
Weight range[−35, 35]
Minimum error1e−4
GANo. individuals100
Crossover operatortwo point
Crossover probability0.8
Mutation operatorbinary
Mutation probability0.02
Objective functionRMSE

Using the above methodology, several network architectures were tested to select the most accurate topology. The most basic architecture usually used for this type of application involves a neural network consisting of three layers [29]. The number of hidden neurons should be sufficient for ensuring an adequate representation of the information in the data used for training the network [2]. There is no specific approach for determining the number of neurons of the hidden layer (NHL), but many alternative combinations are possible. Therefore, the optimum number of neurons was determined by systematically adding neurons and evaluating the RMSE of the sets during the learning process [29]. Fig. 3 shows the RMSE found when correlating λ(T,P) as a function of the number of neurons in the hidden layer. As observed in this figure, the optimum number of neurons in the hidden layer is between 4 and 6. The network giving the lowest deviation during training was the one with 31 parameters in the input layer, 5 neurons in the hidden layer, and one neuron in the output layer. The RMSEs of this architecture were 0.00199 during training, and 0.00197 during prediction, respectively.

Fig. 3

Root mean square error (RMSE) found in correlating the λ(T,P) of ILs as a function of the number of neurons in the hidden layer (NHL) using ANN+GA. Training step (●) and prediction step (◊).

4 Results and discussion

Once the best architecture was determined, the optimum weights required to carry out the estimate λ(T,P) of ILs were obtained. Table 3 shows the optimum weights for the ANN+GA 31-5-1.

Table 3

Optimum weights obtained by ANN+GA for model 31-5-1.

InputWij12345
T1−0.783160.00703−0.179600.261680.04082
P20.674040.01776−0.28596−0.081770.02026
M30.49567−0.42633−0.336460.51141−1.92883
Imidazolium(+)4−0.16750−1.164202.856271.423020.53367
Pyridinium(+)52.239310.576466.45933−0.866210.11607
Pyrrolidinium(+)6−0.755671.19110−0.72874−1.65053−0.62319
Ammonium(+)7−0.431860.060060.10105−1.37138−0.54707
Phosphonium(+)8−0.81134−0.47999−9.264921.617492.96438
–H(+)9−0.93286−0.955189.891400.154190.19737
CH3(+)10−0.76023−0.64633−8.345960.688901.17553
–CH2(+)110.72921−2.0156834.540030.97577−2.75101
CH3(−)12−0.357010.05285−0.080030.553030.21510
–CH2(−)131.66628−0.28268−0.15251−0.10097−0.15836
>CH–(−) [–CH–](−)14−0.32900−0.328690.25780−0.75351−0.34744
>C<(−) [>C–](−)15−0.36681−0.15418−0.07819−0.26981−0.01303
–COO–(−)160.804871.64832−0.74100−1.342530.33933
–COOH(−)170.34297−0.55203−0.215940.932110.86835
–OH(−)18−0.046080.007960.798310.09337−0.41559
–O–(−) [–O](−)19−0.916700.790960.748440.230562.02816
O(−)20−0.27203−0.119961.04748−1.99322−0.72868
–CN(−)210.683361.58632−0.89114−1.318130.43514
>N–(−)22−0.18281−1.47621−1.277672.788031.16030
NH2(−)23−1.181640.24742−0.035500.13463−0.52892
–SO2(−)24−0.31328−2.36592−0.259110.78516−1.13372
–SH(−)25−0.012350.03879−0.276300.87636−0.11776
CF3(−)26−1.00636−1.26141−1.997340.771690.57938
–CF2(−)27−0.368831.13458−2.08071−0.36091−1.08332
–F(−)281.22737−0.579600.294771.417481.66512
–Cl(−)291.06402−0.389470.04274−0.52829−0.41591
–P(−)30−0.604280.083092.65385−0.69970−0.41548
–B(−)31−0.38003−0.20115−2.50992−0.15801−0.06608
W0j1.52286−0.415960.31931−0.68584−2.44430
Wnj0.29242−13.790250.20755−13.640707.61938
Wn07.26828

The accuracy of the chosen network was verified between the calculated values of λ and the experimental data from the literature, by using the average relative absolute deviation for each data point (|%Δλ|) and for the total set (AARD). These deviations were calculated as follows:

|%Δλ|=|λcalcλexpλexp|100(4)
AARD=100Ndatai=1N|λcalcλexpλexp|i(5)

Table S1 shows the results obtained for the correlation set and prediction set, respectively (see, Supplementary data). Table 4 shows the results obtained from the proposed GCM for all different types of ILs; it shows the following ranges of AARD in the prediction of λ(T,P): ammonium < pyridinium < imidazolium < phosphonium < 1%, and pyrrolidinium < 3%. All IL-types presented a maximum error lower than 4%, and all the correlation coefficients R2 are greater than 0.9. Table 5 shows a summary of the deviations of all ILs using the proposed ANN+GA method. The results showed that the ANN+GA can provide an accurate estimation of the thermal conductivity of several ILs: AARD lower than 0.91% for the 300 data points used in the training set, and an AARD lower than 0.84% for the other 100 data points used in the prediction set. The AARD is a little higher than 0.8% for all datasets (400 data points of several ILs), and the deviation is lower than 5% for all data points of the database.

Table 4

Deviations obtained with the proposed ANN+GA in the estimation of thermal conductivity of ILs.

IL classNdataΔλ (W m−1 K−1)ΔT (K)ΔP (kPa)AARD|%Δλ|min|%Δλ|maxR2
Training set
Imidazolium1990.12–0.22273−390100−200001.270.004.960.9961
Ammonium300.12–0.17273−3531010.210.000.980.9997
Pyridinium180.16–0.17294−334100−200000.240.000.550.9934
Phosphonium460.12–0.16282−3551010.190.000.640.9997
Pyrrolidinium70.11–0.13293−3531010.220.030.430.9937
Prediction set
Imidazolium420.13–0.21273−390100−200000.680.032.220.9987
Ammonium70.16–0.17298−3531010.090.000.160.9981
Pyridinium90.16–0.17294−334100−200000.380.001.080.9896
Phosphonium330.13–0.16286−3531010.850.003.420.9878
Pyrrolidinium90.12–0.13293−3331012.671.793.530.9288
Table 5

Summary of deviations in the estimation of thermal conductivity of ILs using the ANN+GA method.

StatisticsCorrelation setPrediction setTotal set
Ndata300100400
AARD0.910.840.89
|%Δλ|min0.000.000.00
|%Δλ|max4.963.534.96
|%Δλ|<5%300100400
|%Δλ|>5%000
R20.99690.99630.9967

Fig. 4 shows a comparison between experimental (solid line) and calculated (dots) values by the proposed ANN+GA method to estimate λ(T,P). This figure shows a comparison between correlated and experimental values of thermal conductivity. For the training set, the correlation coefficient R2 was 0.9969, and the curve slope (m) was 0.9947 (expected to be 1.0). Also, the figure shows a comparison for the prediction set between calculated and experimental values of λ. In this case, the correlation coefficient R2 is 0.9963 and m (also expected to be 1.0) is 0.9926. For the total set, the correlation coefficient R2 is 0.9967 and m is 0.9945. This figure ratifies the discussion presented above.

Fig. 4

Accuracy in the estimation of thermal conductivity of ILs with the proposed ANN+GA method. (○) training set, and (×) prediction set.

Note that for the estimation of λ(T,P = 101.325 kappa), the prediction set presents a very low deviation (AARD = 0.88%) with an |%Δλ|max of 3.53%. In particular, the dataset of 1-butyl-1-methylpyrrolidinium bis[(trifluoromethyl)sulfonyl]imide [20] showed a highest AARD of 2.92%. In addition, the prediction set also showed a very low AARD of 0.66% with a |%Δλ|max of 2.22% for the estimation of the thermal conductivity at different temperatures and pressures λ(T,P). In this case, the prediction set presents an AARD of 0.40% for pressures higher than the atmospheric pressure (P > 101.325 kappa).

Recently a number of methods have been proposed for the estimation of thermal conductivity of ILs [30] (see, the Introduction section). It is worth pointing out that all these methods obtained their results from different databases with different correlation and prediction sets (or training and prediction sets), and based on different methodologies and the results cannot be compared directly with one another. However, a comparison can be made for the selected datasets of common ILs for all methods. Table 6 shows a comparison between the proposed neural network method for the λ estimation of ILs and other methods reported in the literature such as group contribution methods (GCM), generalized correlations, and quantitative structure–property relationships (QSPR). This table contains 20 datasets of ILs at atmospheric pressure and several temperatures. Only four methods can be completely compared based on these results. Albert−Müller's GCM [8] showed an AARD of 2.19, Wu's method [9] presented an AARD of 1.89, Lazzús's QSPR model [30] resulted in an AARD of 2.09, and the proposed neural network method showed an AARD of 1.35. On the one hand, Gardas–Coutinho's GCM [7] does not perform for any ILs due to the group division principle, and it cannot be extended to others ILs [30]. On the other hand, Wu's method [9] and Shojaee's method [3] depend on the knowledge of other properties of ILs. Therefore, the proposed neural network method presented better accuracy and more advantages than other methods to predict λ of ILs at atmospheric pressure and different temperatures.

Table 6

Comparison between the proposed neural network method and other methods reported in the literature for the λ(T) estimation of ILs.

Ionic liquidΔT (K)NdataGardas−Coutinho
GCM [7]
AARD
Albert−Müller
GCM [8]
AARD
Wu et al.
Correlation [9]
AARD
Shojaee et al.
Correlation [3]
AARD
Lazzús
QSPR [30]
AARD
This work
ANN+GA
Ref.
[C2mim][NTf2]273.15–353.1596.414.025.846.543.693.83[6]
[C2mim][EtSO4]293.00–353.0070.152.001.309.343.500.87[20]
273.15–353.1592.531.1011.160.890.74[6]
[C2mim][CH3COO]273.17–353.1590.850.863.320.25[6]
[C2mim][DCA]273.17–353.1591.760.5814.720.800.31[6]
[C2mim][C(CN)3]273.17–353.1590.330.741.520.27[6]
[C4mim][OTf]293.00–353.0070.151.710.371.920.551.68[20]
293.00–353.0073.513.093.823.001.74[21]
[C4mim][NTf2]293.00–353.0071.262.651.003.293.212.38[20]
[C4mim][PF6]293.00–353.0070.980.881.552.400.960.26[21]
294.90–335.1030.691.401.010.34[5]
[C6mim][NTf2]293.00–353.0071.401.691.973.662.653.03[20]
273.15–353.1595.261.885.341.000.78[6]
[C6mim][BF4]293.00–353.00723.094.551.200.760.722.03[21]
[C6mim][PF6]293.00–353.0074.322.071.512.274.001.21[21]
294.10–335.2030.951.781.350.95[5]
[C8mim][PF6]295.10–335.2031.202.073.023.150.49[5]
[C8mim][NTf2]293.00–353.0071.093.413.011.010.46[20]
[C4mpyrr][NTf2]293.00–323.0040.192.650.776.052.982.92[20]
293.00–333.0055.844.7610.912.522.47[21]
AARD3.472.191.895.872.091.35

In other comparison, using Shojaee's method [3] at different pressures (100−20000 kappa) resulted in an AARD of 3.47% for three ILs, while the proposed neural network method showed an AARD of 0.64%. Table 7 shows the ILs and ranges of temperatures and pressures used in this comparison. Table 7 clearly shows that the proposed neural network method offered better accuracy and more advantages in the λ estimation of ILs at several temperatures and pressures only by knowing the molecular structure of IL.

Table 7

Comparison between the proposed neural network method and other method reported in the literature for the λ(T,P) estimation of ILs.

Ionic liquidΔT (K)ΔP (kPa)NdataShojaee et al.
Correlation [3]
AARD
This work
ANN+GA
AARD
Ref.
[C4mim][PF6]294.90–335.10100–2000093.160.47[5]
[C6mim][PF6]294.10–335.20100–2000092.970.95[5]
[C8mim][PF6]294.10–335.20100–2000094.290.49[5]
AARD3.470.64

In addition to these results, a comparison was made with a neural network with standard back-propagation (BPNN), and similar architecture (31-5-1) and database. The BPNN showed an AARD of 2.5%, and a |%Δλ|max higher than 15%. Fig. 5 shows a |%Δλ| found during the prediction of thermal conductivity of ILs using the BPNN versus a |%Δλ| obtained using the proposed ANN+GA method. In other applications relative to BPNN, Hezave et al. [12] presented a neural network model to predict thermal conductivity of ionic liquids. In Hezave's neural network model a total of 209 data points from 21 different ionic liquids were used to train and test the BPNN, and the optimum architecture was determined to be 4-13-1, with input parameters such as temperature T, pressure P, molecular weight MW, and melting point temperature Tm. A considerable disadvantage of Hezave's neural network model is the use of melting point temperature of ILs. Note that, for organic compounds, Pérez Ponce et al. [31] showed an AARD higher than 7.5% for several methods reported during years 1987–2006; additionally, they presented an AARD for three commercial softwares for predicting Tm with an AARD higher than 12% and maximum deviations higher than 45%. For ILs, Huo et al. [32] presented a group contribution method (GCM) for estimating the melting point of imidazolium-ILs with an AARD of 5.9% and maximum deviations higher than 32%. Coutinho et al. [33] presented a critical review on predictive methods for the estimation of thermophysical properties of ILs, with R2 from 0.6 to 0.9 for the prediction of Tm. Thus, the reliability of the input data is never established, making interpretation of deviations between the Hezave's neural network model and experimental data of thermal conductivity of ILs impossible to be made. In contrast, low deviations from the proposed ANN+GA method (an AARD a little higher than 0.8% and a |%Δλ|max below than 5%) indicate that it can estimate λ(T,P) of ILs with low deviations and can be relied on an accuracy of 99%. All these results represent a big increase in accuracy for predicting this important property, and show that the application of the proposed ANN+GA method was crucial. An important observation that is worth pointing out is the influential effects of the structural groups in the correlation and prediction of the thermal conductivity. Considering that the big differences in the chemical structure and physical properties of the ILs used in the study place additional difficulties that the proposed ANN+GA method was able to handle.

Fig. 5

Deviations found in the prediction of λ(T,P) of ILs using ANN+GA (○) and BPNN (Ж).

The incorporation of a GA for the optimization of the neural network weights has positive effects on architecture reducing the number of neurons in the hidden layer while controlling the selection of network connection weightings in a wide numeric range (see, Table 2). Note that, traditional optimization techniques such as the back-propagation learning algorithm can also determine the number of network parameters, such as network connection weightings; however, traditional optimization techniques are not able to control the network parameter optimizations. In contrast, the GA was able to solve this problem.

5 Conclusions

In this study, the thermal conductivity of several ILs at different temperatures and pressures was correlated and estimated using an artificial neural network optimized with a genetic algorithm in order to update the weights of the network.

Based on the results and discussion presented in this study, the following conclusions are obtained: (i) the great differences in the chemical structure, and the physical properties of the ionic liquids considered in the study impose additional difficulties to the problems that the proposed ANN+GA method has been able to handle; (ii) the results show that the ANN+GA method can estimate the thermal conductivity of ionic liquids at different temperatures and pressures with low deviations. The consistency of the method has been checked using experimental values of thermal conductivity and comparing them with the calculated values by the proposed method; (iii) the values calculated using the proposed model are believed to be accurate enough for engineering calculations, for generalized correlations and for equations of state methods, among other uses; and (iv) the incorporation of a GA for the optimization of the neural network weights has positive effects on architecture reducing the number of neurons in the hidden layer while controlling the selection of network connection weightings in a wide numeric range.

Acknowledgements

The author thanks the Direction of Research of the University of La Serena (DIULS), through the research project PT13144. Special thanks go to the Department of Physics of the University of La Serena (DFULS) for the special support that made possible the preparation of this paper.


References

[1] H. Ohno Electrochemical Aspects of Ionic Liquids, John Wiley & Sons, New Jersey, 2005

[2] J.A. Lazzús Int. J. Thermophys., 30 (2009), pp. 883-909

[3] S.A. Shojaee; S. Farzam; A.Z. Hezave; M. Lashkarbolooki; S. Ayatollahi Fluid Phase Equilib., 354 (2013), pp. 199-206

[4] C. Frez; G.J. Diebold; C.D. Tran; S. Yu J. Chem. Eng. Data, 51 (2006), pp. 1250-1255

[5] D. Tomida; S. Kenmochi; T. Tsukada; K. Qiao; C. Yokoyama Int. J. Thermophys., 28 (2007), pp. 1147-1160

[6] A.P. Fröba; M.H. Rausch; K. Krzeminski; D. Assenbaum; P. Wasserscheid; A. Leipertz Int. J. Thermophys., 31 (2010), pp. 2059-2077

[7] R.L. Gardas; J.A.P. Coutinho AIChE J., 55 (2009), pp. 1274-1290

[8] J. Albert; K. Müller Ind. Eng. Chem. Res., 53 (2014), pp. 17522-17526

[9] K.-J. Wu; C.-X. Zhao; C.-H. He Fluid Phase Equilib., 339 (2013), pp. 10-14

[10] J.O. Valderrama; W.W. Sanga; J.A. Lazzús Ind. Eng. Chem. Res., 47 (2008), pp. 1318-1330

[11] M.T. Hagan; H.B. Demuth; M. Beal Neural Network Design, PWS Publishing Company, Boston, USA, 1996

[12] A.Z. Hezave; S. Raeissi; M. Lashkarbolooki Ind. Eng. Chem. Res., 51 (2012), pp. 9886-9893

[13] J. Holland Adaptation in Natural and Artificial Systems, University of Michigan Press, Michigan, USA, 1975

[14] L. Davis Handbook of Genetic Algorithms, Van Nostrand Reinhold, New York, USA, 1991

[15] J.A. Lazzús; M. Rivera Rev. Mex. Fís., 59 (2013), pp. 577-583

[16] Q. Dong; C.D. Muzny; A. Kazakov; V. Diky; J.W. Magee; J.A. Widegren; R.D. Chirico; K.N. Marsh; M. Frenkel J. Chem. Eng. Data, 52 (2007), pp. 1151-1159

[17] Q.-L. Chen; K.-J. Wu; C.-H. He J. Chem. Eng. Data, 58 (2013), pp. 2058-2064

[18] A.G.M. Ferreira; P.N. Simoes; A.F. Ferreira; M.A. Fonseca; M.S.A. Oliveira; A.S.M. Trino J. Chem. Thermodyn., 64 (2013), pp. 80-92

[19] R.L. Gardas; R. Ge; P. Goodrich; C. Hardacre; A. Hussain; D.W. Rooney J. Chem. Eng. Data, 55 (2010), pp. 1505-1515

[20] R. Ge; C. Hardacre; P. Nancarrow; D.W. Rooney J. Chem. Eng. Data, 52 (2007), pp. 1819-1823

[21] C.A. Nieto de Castro; M.J.V. Lourenco; A.P.C. Ribeiro; E. Langa; S.I.C. Vieira J. Chem. Eng. Data, 55 (2010), pp. 653-661

[22] D. Tomida; S. Kenmochi; K. Qiao; T. Tsukada; C. Yokoyama Fluid Phase Equilib., 340 (2013), pp. 31-36

[23] D. Tomida; S. Kenmochi; T. Tsukada; K. Qiao; Q. Bao; C. Yokoyama Int. J. Thermophys., 33 (2012), pp. 959-969

[24] M.E. Van Valkenburg; R.L. Vaughn; M. Williams; J.S. Wilkes Thermochim. Acta, 425 (2005), pp. 181-188

[25] J.A. Lazzús Fluid Phase Equilib., 313 (2012), pp. 1-6

[26] Z.-G. Che; T.-A. Chiang; Z.-H. Che Int. J. Innovative Comput., 7 (2011), pp. 5839-5850

[27] A. Meng; L. Ye; D. Roy; P. Padilla Comput. Educ., 49 (2007), pp. 1205-1223

[28] E. Eiben; R. Hinterding; Z. Michalewicz IEEE Trans. Evolut. Comput., 3 (1999), pp. 124-141

[29] J.A. Lazzús; I. Salfate; S. Montecinos Neural Netw. World, 24 (2014), pp. 601-617

[30] J.A. Lazzús; G. Pulgar-Villarroel J. Mol. Liq., 211 (2015), pp. 981-985

[31] A.A. Pérez Ponce; I. Salfate; G. Pulgar-Villarroel; L. Palma-Chilla; J.A. Lazzús J. Eng. Thermophys., 22 (2013), pp. 226-235

[32] Y. Huo; S. Xia; Y. Zhang; P. Ma Ind. Eng. Chem. Res., 48 (2009), pp. 2212-2217

[33] J.A.P. Coutinho; P.J. Carvalho; N.M.C. Oliveira RSC Adv., 2 (2012), pp. 7322-7346


Comments - Policy