Possible Interactions and Interferences of Copper , Chromium , and Arsenic during the Gasification of Contaminated Waste Wood †

A considerable proportion (about 64%) of biomass energy is produced from woody biomass (wood and its wastes). However, waste wood (WW) is very often contaminated with metal(loid) elements at concentrations leading to toxicity emissions and damages to facilities during thermal conversion. Therefore, procedures for preventing and/or alleviating the negative impacts of these elements require further development, particularly by providing informative and supportive information regarding the phase transformations of the metal(loid)s during thermal conversion processes. Although it is well known that phase transformation depends on different factors such as elements’ vaporization characteristics, operational conditions, and process configuration; however, the influences of reaction atmosphere composition in terms of interactions and interferences are rarely addressed. In response, since Cu, Cr, and As (CCA-elements) are the most regulated elements in woody biomass, this paper aims to explore the possible interactions and interferences among CCA-elements themselves and with Ca, Na, S, Cl, Fe, and Ni from reaction atmosphere composition perspectives during the gasification of contaminated WW. To do so, thermodynamic equilibrium calculations were performed for Boudouard reaction (BR) and partial combustion reaction (PCR) with temperature ranges of 0–1300 ◦C and 0–1800 ◦C, respectively, and both reactions were simulated under pressure conditions of 1, 20, and 40 atm. Refinement of the occurred interactions and interferences reveals that Ni-As interactions generate dominant species As2Ni5 and As8Ni11, which increase the solid–gaseous transformation temperature of As. Moreover, the interactions between Ca and Cr predominantly form C3Cr7; whereas the absence of Ca leads to Cr2Na2O4 causing instability in the Cr phase transformation.


Introduction
Waste wood (WW) is the term used for wood material that has previously been used for various purposes and is now being treated as a waste material; in the context of this study the waste wood is being utilized for energy production [1].It comprises a wide variety of wood materials primarily from industrial and commercial activities and from construction and demolition operations [2,3].Since it is collected from diverse sources, WW may contain both physical and chemical contaminants [1,4].Physical contaminates can simply be separated by sorting or using mechanical techniques.However, the majority of chemical contaminants are usually linked to the metal(loid) elements contained in wood Energies 2018, 11, 1966 2 of 17 that has been treated with preservatives, paintings, coatings, and other related activities, and, of course, mechanical separation is not applicable.
Essentially, gasification process involves the conversion of solid carbonaceous fuels into combustible gases composed of a mixture of N 2 , H 2 O, CH 4 , H 2 , CO 2 , CO, and very small quantities of low molecular weight hydrocarbons and contaminants including ash, tar, and carbon particles [5,6].In general, biomass gasification consists of four processing steps: drying, pyrolysis of dried biomass particles (de-volatilization), partial oxidation of pyrolysis gases and/or char, and char gasification (reduction) [7], as depicted in Figure 1 [8].Technically, gasification is an exothermic partial oxidation of biomass operating at high temperatures between 800 and 1300 • C with the presence of gasifying mediums such as carbon dioxide, steam, and air (or a mixture of air and steam) [7,[9][10][11].Further, the compositions and properties of the gasification products vary according to the operational conditions, gasifier type, and biomass materials.
Energies 2018, 11, x FOR PEER REVIEW 2 of 17 However, the majority of chemical contaminants are usually linked to the metal(loid) elements contained in wood that has been treated with preservatives, paintings, coatings, and other related activities, and, of course, mechanical separation is not applicable.Essentially, gasification process involves the conversion of solid carbonaceous fuels into combustible gases composed of a mixture of N2, H2O, CH4, H2, CO2, CO, and very small quantities of low molecular weight hydrocarbons and contaminants including ash, tar, and carbon particles [5,6].In general, biomass gasification consists of four processing steps: drying, pyrolysis of dried biomass particles (de-volatilization), partial oxidation of pyrolysis gases and/or char, and char gasification (reduction) [7], as depicted in Figure 1 [8].Technically, gasification is an exothermic partial oxidation of biomass operating at high temperatures between 800 °C and 1300 °C with the presence of gasifying mediums such as carbon dioxide, steam, and air (or a mixture of air and steam) [7,[9][10][11].Further, the compositions and properties of the gasification products vary according to the operational conditions, gasifier type, and biomass materials.Currently, because of its efficiency and flexibility [12], gasification is vastly employed for generating renewable energy and waste disposal.However, when conducting gasification, metal(loid)s in biomass can cause severe human health and environmental problems as well as significant technical damages when they present at high concentrations.Hence, to avoid and/or limit the risks of metal(loid) elements, the transformation behaviours of these elements across the different phases require further investigating.
A range of studies [13][14][15][16][17] have focused on the distribution of the metal(loid)s between gaseous and solid phases during gasification.These studies report that the partitioning of metal(loid) elements can be directly influenced by a combination of several factors such as feedstock characteristics, operational conditions, and process agents.In fact, phase transformation temperatures of metal(loid)s cannot be accurately determined using traditional analytical tools/instruments because the generated species cannot be easily identified, and hence measured.As such, the mass distribution of metal(loid) elements during gasification is a quite complex problem.Fortunately, thermodynamic equilibrium modelling based on Gibbs energy minimization has effectively been used for yielding the best predictions of the behaviours of metal(loid) elements, and in turn allowing for precise evaluation of the speciation during the thermochemical processes [18][19][20][21][22][23][24][25].
Jiang et al. [18] studied the solid-gaseous transition of Zn, Pb, Ni, Mn, Mg, Fe, Cu, Cr, Co, Cd, As, and Al in five biomass samples during typical gasification scenarios.This study reports that the most volatilized elements are As, Cd, Pb, and Zn, whereas Cu, Co, Ni, and Mn are moderately volatilized, and the phase transformation temperatures of Mg, Fe, Cr, and Al are > 1200 °C.Froment et al. [19] assessed the volatilization and condensation of certain elements in woody biomass during gasification with temperatures ranging from 500 °C to 1500 °C, and pressure conditions from 1 to 10 bar.The calculation results showed that Zn, Pb, Se, Sb, S, N, Hg, F, Cd, and Cl completely volatilized Currently, because of its efficiency and flexibility [12], gasification is vastly employed for generating renewable energy and waste disposal.However, when conducting gasification, metal(loid)s in biomass can cause severe human health and environmental problems as well as significant technical damages when they present at high concentrations.Hence, to avoid and/or limit the risks of metal(loid) elements, the transformation behaviours of these elements across the different phases require further investigating.
A range of studies [13][14][15][16][17] have focused on the distribution of the metal(loid)s between gaseous and solid phases during gasification.These studies report that the partitioning of metal(loid) elements can be directly influenced by a combination of several factors such as feedstock characteristics, operational conditions, and process agents.In fact, phase transformation temperatures of metal(loid)s cannot be accurately determined using traditional analytical tools/instruments because the generated species cannot be easily identified, and hence measured.As such, the mass distribution of metal(loid) elements during gasification is a quite complex problem.Fortunately, thermodynamic equilibrium modelling based on Gibbs energy minimization has effectively been used for yielding the best predictions of the behaviours of metal(loid) elements, and in turn allowing for precise evaluation of the speciation during the thermochemical processes [18][19][20][21][22][23][24][25].
Energies 2018, 11, 1966 3 of 17 Jiang et al. [18] studied the solid-gaseous transition of Zn, Pb, Ni, Mn, Mg, Fe, Cu, Cr, Co, Cd, As, and Al in five biomass samples during typical gasification scenarios.This study reports that the most volatilized elements are As, Cd, Pb, and Zn, whereas Cu, Co, Ni, and Mn are moderately volatilized, and the phase transformation temperatures of Mg, Fe, Cr, and Al are > 1200 • C. Froment et al. [19] assessed the volatilization and condensation of certain elements in woody biomass during gasification with temperatures ranging from 500 to 1500 • C, and pressure conditions from 1 to 10 bar.The calculation results showed that Zn, Pb, Se, Sb, S, N, Hg, F, Cd, and Cl completely volatilized as either elemental or species, whereas, Ca, Si, Al, P, Cu, Mn, Ti, Fe, Cr, Mg, and Ba remained as oxides in the condensed phase, except copper which condensed as Cu 3 As.Kramb et al. [20] modelled the behavior of As during the gasification of CCA-treated wood, and reported that about 99.6% of As can be captured.Thompson and Argent [21] categorized the trace elements in a PRENLFO gasifier into "virtually immobile (Cr, Ba, Ni, V), somewhat mobile (Mo, Be, Cu), mobile (Sb, B) and almost totally mobilized (As, Cd, Pb, Sn, Zn)".Kilgallon et al. [22] investigated the impacts of both pressure conditions and chlorine/sulphur on element phase transitions.Reed et al. [23] reported the speciation prediction of Cd, As, B, Co, Be, Pb, Cu, Sb, V, Sn, and Zn in fuel gas from a gasifier.They explicitly partitioned these elements into three groups: condensed phase (Cu, Be, V, Zn and Co; at T 50 900 K), gas phase (As without Ni, B without Ca, Sb, and Sn), and condensed phase (As with Ni, B with Ca, Cd, and Pb; at T 50 < 900 K).Liu et al. [24] also reported the speciation of Ni, As, Pb, Cd, Se, and Sb during underground coal gasification with specific agent injection and pressure.A study by Bradshaw et al. [25] showed that the elements Hg, As, Pb, Sb, Cd, Sn, Se, and Te, and alkali metals remained in fuel gas, and the species of Pb, Sb, and Cd can pass through the path of fuel gas forming surface deposits.
Although the modelling of phase transformation yields beneficial information (e.g., elemental vaporization and condensation, optimized operational conditions) about the distribution of metal(loid) elements, factors including the interactions and interferences among these elements cannot be ignored.In particular, metal(loid)s can interact and/or interfere, leading to the generation of new species, which may induce remarkable influence on the volatilization of metal(loid)s.At present, to the best of the authors' knowledge, there has been no single study towards the assessment of occurred interferences between metal(loid)s during thermochemical processes.In fact, only a study by K. Svoboda [26] investigated the instantaneous interferences of both HF and HCl on H 2 S sorption.Also, two studies evaluated the effect of interactions on elemental behaviour during coal combustion.Contreras et al. [27] evaluated the interactions between Cd, As, Sb, and Hg, and pointed out that some species can be dominant under specific conditions.For example, As-Cd interaction forms Cd 3 (AsO 4 ) 2 , which reduces the concentration of gaseous oxides.Yan et al. [28] showed that interactions among Se, Hg, and As lead to the generation of HgSe (g) and AsSe (g) as dominant species under reduced conditions.
In this paper, thermodynamic equilibrium calculations were carried out to explore the possible interactions and interferences between Cu, Cr, and As (CCA-elements) themselves and with Ca, Na, S, Cl, Fe, and Ni in contaminated WW during two reactions of gasification: Boudouard reaction (BR) and partial combustion reaction (PCR).The calculations were performed using MTDATA software (Version 5.10, NPL, London, UK) to predict the phase transformation temperatures and speciation of the considered metal(loid)s.Moreover, contaminated WW samples with complete composition data were selected from the ECN PHYLLIS2 database [29] for the modelling.

Thermodynamic Equilibrium Modelling
Thermodynamic equilibrium modellings for Boudouard reaction (BR) and partial combustion reaction (PCR) were carried out using the multiphase module of the MTDATA software coupled with the Scientific Group Thermodata Europe (SGTE) database.With each modelling, the calculations were conducted to study the effects of the occurred interactions and interferences on the fate of the CCA-elements as consequence of different contents of CCA-elements and of Ca, Na, S, Cl, Fe, and Ni in WW samples.
For model simulations, a primary search of the SGTE database was carried out to reveal the relevant species of CCA-elements.A complete list of the adopted species are given in Table 1.The best fit probability distributions of the composition data of the selected WW samples (#2712, #679, #1779, #1364, #871, #1448, #3498, #2901, #2900) [29] are estimated through the @Risk analysis simulation package (Version 6.1, Palisade Corp., Ithaca, NY, USA), and the estimation results are exhibited in Table 2.As an illustrative example, the statistical analysis of the moisture content is demonstrated in Figure 2. The elements to be introduced to the MTDATA were grouped into major elements (Ar, C, N, O, and H) and minor elements (As, Cr, Cu, Co, Cl, Ca, Ni, S, Na, and Fe).Note that the mass ranges of the minor elements, Table 3, were assumed inconsistent with the information provided by the ECN PHYLLIS2 database [29] concerning the contaminated woody biomass.Moreover, the operational conditions and process agents for each reaction are presented in Table 4, in which the required oxygen for PCR and the required carbon dioxide for BR are determined using Equation (1) [30] and Equation (2) [7], respectively, as follows: Note that the air composition (under standard dry at atmospheric pressure) consists of 1.28% Ar, 23.2% O 2 , 0.046% CO 2 , and 75.47%N 2 .Biomass feed rates for both reactions were calculated by means of the higher heating value (HHV d ) from Equation (3) [30].

Equilibrium Composition of CCA-Elements
Initially, the modelling of PCR was based on composition data in Table 2 (the mean values) and process parameters of PCR in Table 4 in order to predict the equilibrium partitioning of the CCA-elements.The simulation results show that, at temperatures ≤1350 • C, about 100% of Cu remained in solid phase as a dominant species of CaCu resulting from Ca-Cu interaction (Figure 3).When the temperature increases from 1350 to 1800 • C, new gaseous species CuH (0.47-4.49%),Cu 2 (0.02-0.43%), and CuCl (5.5 × 10 −3 -5.544 × 10 −1 %) are formed due to Cu-Cl and Cu-H interactions.At temperatures between 1450 and 1600 • C, Cu is proportionally distributed as elemental between gaseous phase (4.29-94.52%)and solid phase (88.84-33.57%).
elements.The simulation results show that, at temperatures ≤1350 C, about 100% of Cu remained in solid phase as a dominant species of CaCu resulting from Ca-Cu interaction (Figure 3).When the temperature increases from 1350 C to 1800 C, new gaseous species CuH (0.47-4.49%),Cu2 (0.02-0.43%), and CuCl (5.5 × 10 −3 -5.544 × 10 −1 %) are formed due to Cu-Cl and Cu-H interactions.At temperatures between 1450 C and 1600 C, Cu is proportionally distributed as elemental between gaseous phase (4.29-94.52%)and solid phase (88.84-33.57%).The predicated distribution of Cr is shown in Figure 4. Obviously, at temperatures ≤ 1550 C, Cr remained in solid phase as C3Cr7 (100%), and this was mainly due to C-Cr interaction.Additionally, at temperatures between 1550 C and 1800 C, the mass percentage of C3Cr7 decreased from 100% to 67.34%, while at the same time it was released as elemental Cr (3.16-32.48%)and as a new set of gaseous species generated by Cr-Cl, Cr-H, Cr-N, and Cr-O interactions; the total mass percentage of these gaseous species was less than 0.2%.Figure 5 shows the equilibrium distribution of As.It is clear that at temperatures ≤ 1400 C, As is predicated as stable solid species AsNa3 (100%) and As2Ni5 (99.99-99.15%)as consequence of As- The predicated distribution of Cr is shown in Figure 4. Obviously, at temperatures ≤1550 • C, Cr remained in solid phase as C 3 Cr 7 (100%), and this was mainly due to C-Cr interaction.Additionally, at temperatures between 1550 and 1800 • C, the mass percentage of C 3 Cr 7 decreased from 100% to 67.34%, while at the same time it was released as elemental Cr (3.16-32.48%)and as a new set of gaseous species generated by Cr-Cl, Cr-H, Cr-N, and Cr-O interactions; the total mass percentage of these gaseous species was less than 0.2%.
elements.The simulation results show that, at temperatures ≤1350 C, about 100% of Cu remained in solid phase as a dominant species of CaCu resulting from Ca-Cu interaction (Figure 3).When the temperature increases from 1350 C to 1800 C, new gaseous species CuH (0.47-4.49%),Cu2 (0.02-0.43%), and CuCl (5.5 × 10 −3 -5.544 × 10 −1 %) are formed due to Cu-Cl and Cu-H interactions.At temperatures between 1450 C and 1600 C, Cu is proportionally distributed as elemental between gaseous phase (4.29-94.52%)and solid phase (88.84-33.57%).The predicated distribution of Cr is shown in Figure 4. Obviously, at temperatures ≤ 1550 C, Cr remained in solid phase as C3Cr7 (100%), and this was mainly due to C-Cr interaction.Additionally, at temperatures between 1550 C and 1800 C, the mass percentage of C3Cr7 decreased from 100% to 67.34%, while at the same time it was released as elemental Cr (3.16-32.48%)and as a new set of gaseous species generated by Cr-Cl, Cr-H, Cr-N, and Cr-O interactions; the total mass percentage of these gaseous species was less than 0.2%.Figure 5 shows the equilibrium distribution of As.It is clear that at temperatures ≤ 1400 C, As is predicated as stable solid species AsNa3 (100%) and As2Ni5 (99.99-99.15%)as consequence of As- Figure 5 shows the equilibrium distribution of As.It is clear that at temperatures ≤1400 • C, As is predicated as stable solid species AsNa 3 (100%) and As 2 Ni 5 (99.99-99.15%)as consequence of As-Na interaction and As-Ni interaction, respectively.As the temperature increases from 1400 to 1800 • C, the mass percentage of As decreases rapidly in solid phase, and it reforms as new gaseous species predominately by interactions.Among them, the dominant ones are As 2 (1.90-46.00%),AsH (0.24-7.13%), and As (5-46.55%).Indeed, AsH is generated by As-H interaction.
Na interaction and As-Ni interaction, respectively.As the temperature increases from 1400 C to 1800 C, the mass percentage of As decreases rapidly in solid phase, and it reforms as new gaseous species predominately by interactions.Among them, the dominant ones are As2 (1.90-46.00%),AsH (0.24-7.13%), and As (5-46.55%).Indeed, AsH is generated by As-H interaction.To highlight the influence of increased pressure conditions on phase transition temperatures of CCA-elements, the reaction pressure was set to 20 atm and to 40 atm.Not surprisingly, Cr remained completely in condensed phase under both pressures, as shown in Figure 6a.Whereas, under 20 atm, the phase transformation temperatures of Cu and As increased by 300 C and 250 C, respectively, compared to those under 1 atm.Also, under 40 atm, the phase transformation temperatures of Cu and As, respectively, increased by 350 C and 300 C compared to those predicated under 1 atm, as depicted in Figures 6b,c.

Influence of Elemental Contents
Since the elemental content is one of the most important factors that directly affects species formation, different element content in WW samples were simulated to identify the possible interactions and interferences in terms of speciation during PCR, such that the element concentration were introduced to the model as presented in Table 3.Note here that the interactions and interferences that lead to appreciable effects on elemental phase transformation temperature and species formation were ranked as significant.

Influences on CCA-elements
In this case, the concentrations of each CCA-elements were individually changed.The simulation results show that the CCA-elements contents exhibit almost no influence on their phase transformation temperatures, but only concentration shifts are observed as shown in Figure 7a-c.Nevertheless, an exception is that when As concentration is relatively high (5 mole/h), its solidgaseous transformation temperature decreased by 400 C because the Ni contained becomes less than the limit to interact with As, and hence captures 100% of the As.Additionally, As2Ni5 (solid) and AsN (gaseous) are the dominant species of As.To highlight the influence of increased pressure conditions on phase transition temperatures of CCA-elements, the reaction pressure was set to 20 atm and to 40 atm.Not surprisingly, Cr remained completely in condensed phase under both pressures, as shown in Figure 6a.Whereas, under 20 atm, the phase transformation temperatures of Cu and As increased by 300 and 250 • C, respectively, compared to those under 1 atm.Also, under 40 atm, the phase transformation temperatures of Cu and As, respectively, increased by 350 and 300 • C compared to those predicated under 1 atm, as depicted in Figure 6b,c.

Influence of Elemental Contents
Since the elemental content is one of the most important factors that directly affects species formation, different element content in WW samples were simulated to identify the possible interactions and interferences in terms of speciation during PCR, such that the element concentration were introduced to the model as presented in Table 3.Note here that the interactions and interferences that lead to appreciable effects on elemental phase transformation temperature and species formation were ranked as significant.

Influences on CCA-elements
In this case, the concentrations of each CCA-elements were individually changed.The simulation results show that the CCA-elements contents exhibit almost no influence on their phase transformation temperatures, but only concentration shifts are observed as shown in Figure 7a-c.Nevertheless, an exception is that when As concentration is relatively high (5 mole/h), its solid-gaseous transformation temperature decreased by 400 • C because the Ni contained becomes less than the limit to interact with As, and hence captures 100% of the As.Additionally, As 2 Ni 5 (solid) and AsN (gaseous) are the dominant species of As.   8.Note that, in all the simulations, the CCA-elements concentrations were kept fixed (as listed in Table 2) while the concentrations of other minor elements were considered as in Table 3.  8.Note that, in all the simulations, the CCA-elements concentrations were kept fixed (as listed in Table 2) while the concentrations of other minor elements were considered as in Table 3.

Copper
As presented in Figure 8, Cu interacts with Ca and Cl; however, all Cl-Cu interactions release only a single gaseous species-ClCu.The solid and gaseous species formed by Ca-Cu interactions are listed in Table 5.The analysis shows that the dominant species in solid phase are CaCu, CaCu and CaCu5, and CaCu when Ca is present, respectively at concentrations of 250, 500, and 1000 mole/h.On the other hand, the presence of Ca at low concentrations would not be sufficient to interact with Cu; alternatively, it interacts with O and Cr, forming CaCr2O4 and CaO, respectively.Further, it is noted that the volatilization of Cu is not affected by the various concentrations of Ca, which is similar to Cl-Cu interactions.

Chromium
The analysis reveals that Ca plays an important role in Cr distribution across the solid and gaseous phases (Figure 9 and Table 6).That is, Na and Cr interact only when the concentration of Ca is ≤ 100 mole/h, leading to the formation of instable solid specie Cr2Na2O4 between 1100 and 1550 C.CaCr2O4 (solid) is the most likely to be generated by Ca-Cr interactions when Ca is present at 100 mole/h, thereby about 100% of Cr remains in solid phase at temperatures ≤ 1550 C.As the Ca content increases (> 100 mole/h), Ca-Cr interferences and C-Cr interactions take place simultaneously,

Copper
As presented in Figure 8, Cu interacts with Ca and Cl; however, all Cl-Cu interactions release only a single gaseous species-ClCu.The solid and gaseous species formed by Ca-Cu interactions are listed in Table 5.The analysis shows that the dominant species in solid phase are CaCu, CaCu and CaCu 5 , and CaCu when Ca is present, respectively at concentrations of 250, 500, and 1000 mole/h.On the other hand, the presence of Ca at low concentrations would not be sufficient to interact with Cu; alternatively, it interacts with O and Cr, forming CaCr 2 O 4 and CaO, respectively.Further, it is noted that the volatilization of Cu is not affected by the various concentrations of Ca, which is similar to Cl-Cu interactions.

Chromium
The analysis reveals that Ca plays an important role in Cr distribution across the solid and gaseous phases (Figure 9 and Table 6).That is, Na and Cr interact only when the concentration of Ca is ≤100 mole/h, leading to the formation of instable solid specie Cr 2 Na 2 O 4 between 1100 and 1550 • C. CaCr 2 O 4 (solid) is the most likely to be generated by Ca-Cr interactions when Ca is present at 100 mole/h, thereby about 100% of Cr remains in solid phase at temperatures ≤1550 • C. As the Ca content increases (> 100 mole/h), Ca-Cr interferences and C-Cr interactions take place simultaneously, forming C 2 Cr 3 and C 3 Cr 7 in solid phase at temperatures up to 1550 • C. A most apparent observation in this case is the concentration shift in gaseous phase of Cr.

Arsenic
When Ni content is present at various concentrations, all Ni-As interactions generate solid specie As 2 Ni 5 , by which the phase transformation temperature of As (Figure 10) increases to about 450 • C. In contrast, the absence of Ni leads to the formation of As 2 (gaseous) and AsNa 3 (solid) as dominant species, as shown in Figure 11.Similarly, Na-As interactions yield AsNa 3 .While the absence of Na shows no influence on the As phase formation; it enhances Ni to capture As, generating As 2 Ni 5 (100-99.15%) up to 1400 • C, as illustrated in Figure 12.Also, S tends to interact with As, yielding, at very low mass percentage, gaseous specie AsS between 1750 and 1800 • C.
When Ni content is present at various concentrations, all Ni-As interactions generate solid specie As2Ni5, by which the phase transformation temperature of As (Figure 10) increases to about 450 C.In contrast, the absence of Ni leads to the formation of As2 (gaseous) and AsNa3 (solid) as dominant species, as shown in Figure 11.Similarly, Na-As interactions yield AsNa3.While the absence of Na shows no influence on the As phase formation; it enhances Ni to capture As, generating As2Ni5 (100-99.15%) up to 1400 C, as illustrated in Figure 12.Also, S tends to interact with As, yielding, at very low mass percentage, gaseous specie AsS between 1750 C and 1800 C.When Ni content is present at various concentrations, all Ni-As interactions generate solid specie As2Ni5, by which the phase transformation temperature of As (Figure 10) increases to about 450 C.In contrast, the absence of Ni leads to the formation of As2 (gaseous) and AsNa3 (solid) as dominant species, as shown in Figure 11.Similarly, Na-As interactions yield AsNa3.While the absence of Na shows no influence on the As phase formation; it enhances Ni to capture As, generating As2Ni5 (100-99.15%) up to 1400 C, as illustrated in Figure 12.Also, S tends to interact with As, yielding, at very low mass percentage, gaseous specie AsS between 1750 C and 1800 C.When Ni content is present at various concentrations, all Ni-As interactions generate solid specie As2Ni5, by which the phase transformation temperature of As (Figure 10) increases to about 450 C.In contrast, the absence of Ni leads to the formation of As2 (gaseous) and AsNa3 (solid) as dominant species, as shown in Figure 11.Similarly, Na-As interactions yield AsNa3.While the absence of Na shows no influence on the As phase formation; it enhances Ni to capture As, generating As2Ni5 (100-99.15%) up to 1400 C, as illustrated in Figure 12.Also, S tends to interact with As, yielding, at very low mass percentage, gaseous specie AsS between 1750 C and 1800 C.Assessing CCA-elements Reactivity All the minor elements were introduced to the model at the same concentration (0.25 mole/h) to assess the reactivity of the CCA-elements.In this case, the resulted speciation of CCA-elements are demonstrated in Figure 13a-c.Assessing CCA-elements Reactivity All the minor elements were introduced to the model at the same concentration (0.25 mole/h) to assess the reactivity of the CCA-elements.In this case, the resulted speciation of CCA-elements are demonstrated in Figure 13a-c As can be seen in Figure 13a, the behavior of Cu between 0 and 700 C indicates that Cu interacts favorably with As and S, forming the solid species AsCu3 (81.89-78.24%)and Cu2S (100%), respectively.At temperatures above 700 C, Cu is redistributed as elemental and a set of gaseous As can be seen in Figure 13a, the behavior of Cu between 0 and 700 • C indicates that Cu interacts favorably with As and S, forming the solid species AsCu 3 (81.89-78.24%)and Cu 2 S (100%), respectively.At temperatures above 700 • C, Cu is redistributed as elemental and a set of gaseous species, and its phase transformation temperature becomes 1200 • C, which is less than about 11% of that obtained with equilibrium composition (Figure 3).The predictive calculations, as shown in Figure 13b, originate the Cr solid species CoCr 2 O 4 (100% between 0-300 • C), Cr 2 FeO 4 (100% between 350-500 • C), and CaCr 2 O 4 (100% between 550-1450 • C) to the Co-O-Cr, Fe-O-Cr, and Ca-O-Cr interactions, respectively.Figure 13c shows As species formation.Intuitively, AsCu 3 (27.27%between 0-400 • C) and As 8 Ni 11 (72.73% between 0-800 • C) are the most probable solid species from Cu-As and Ni-As interactions, respectively.
In this view, the distribution profiles of CCA-elements during gasification offer important trend towards their emissions control.More specifically, the presence of Ni and S enhances As and Cu capture in ash, respectively.Also, Cu and As interactions cause both of them to remain in solid phase, and the elements Ca, Co, and Fe are favored to Cr retention.

Under Boudouard Reaction (BR)
The modelling of BR was carried out using the composition data and BR process parameters, respectively, as reported in Tables 2 and 4. The simulation scenarios were considered similar to those for the case of PCR.On the whole, the obtained results were similar to those observed with the PCR model, except that at equilibrium partitioning, about 99.99% of the CCA-elements remained in solid phase at temperatures less than 1300 • C; also no gaseous species were formed under pressure conditions of 20 and 40 atm.Table 7 demonstrates the phase and species formation of CCA-elements at equilibrium composition under 1 atm.

Conclusion
This paper introduces two predictive models for PCR and BR to understand the behavior of CCA-elements in terms of interactions and interferences during the gasification of contaminated WW.The refinement of the predicated results from PCR modelling indicates that most of the occurred Ni-As interactions led to about 100% of As remaining in solid phase at temperatures less than 1400

Figure 1 .
Figure 1.Schematic of the processes in a typical gasifier.

Figure 1 .
Figure 1.Schematic of the processes in a typical gasifier.

Figure 2 .
Figure 2. (a) Probability distribution of moisture content; (b) p-value; and (c) interval confidence of probability distribution.

Figure 7 .
Figure 7. Resulting phase transformations of CCA-elements by changing the concentrations of each element individually: (a) copper; (b) chromium; (c) arsenic.

Figure 8 .
Figure 8. Interactions and interferences between the CCA-elements and each of Ca, Na, S, Cl, Fe, and Ni.

Figure 8 .
Figure 8. Interactions and interferences between the CCA-elements and each of Ca, Na, S, Cl, Fe, and Ni.
Other interactions, specifically, S-Cr and Cl-Cr interactions, exhibit no influence on Cr volatilization and solely yield gaseous species CrS, ClCr and Cl 2 Cr, respectively.Energies 2018, 11, x FOR PEER REVIEW 12 of 17 forming C2Cr3 and C3Cr7 in solid phase at temperatures up to 1550 C.A most apparent observation in this case is the concentration shift in gaseous phase of Cr.Other interactions, specifically, S-Cr and Cl-Cr interactions, exhibit no influence on Cr volatilization and solely yield gaseous species CrS, ClCr and Cl2Cr, respectively.

Figure 9 .
Figure 9. Effects of Ca on Cr phase formation.

Figure 9 .
Figure 9. Effects of Ca on Cr phase formation.

Figure 10 .
Figure 10.Arsenic distribution according to the presence/absence of Ni.

Figure 11 .
Figure 11.Arsenic speciation in the absence of Ni.

Figure 12 .
Figure 12.Arsenic speciation in the absence of Na.

Figure 10 .
Figure 10.Arsenic distribution according to the presence/absence of Ni.

Figure 10 .
Figure 10.Arsenic distribution according to the presence/absence of Ni.

Figure 11 .
Figure 11.Arsenic speciation in the absence of Ni.

Figure 12 .
Figure 12.Arsenic speciation in the absence of Na.

Figure 11 .
Figure 11.Arsenic speciation in the absence of Ni.

Figure 10 .
Figure 10.Arsenic distribution according to the presence/absence of Ni.

Figure 11 .
Figure 11.Arsenic speciation in the absence of Ni.

Figure 12 .
Figure 12.Arsenic speciation in the absence of Na.Figure 12. Arsenic speciation in the absence of Na.

Figure 12 .
Figure 12.Arsenic speciation in the absence of Na.Figure 12. Arsenic speciation in the absence of Na.

Figure 13 .
Figure 13.Speciation of CCA-elements when all minor elements are present at the same concentrations: (a) copper; (b) chromium; (c) arsenic.

Figure 13 .
Figure 13.Speciation of CCA-elements when all minor elements are present at the same concentrations: (a) copper; (b) chromium; (c) arsenic.
• C. It also indicates that Ca-Cr interferences form solid specie C 3 Cr 7 , but the absence of Ca causes instability in Cr distribution at temperatures from 1100 to 1550 • C. Whereas Co-O-Cr, Fe-O-Cr, and Ca-O-Cr interactions reduced the Cr volatilization as consequence of the generated solid species CoCr 2 O 4 , Cr 2 FeO 4 , and CaCr 2 O 4 , respectively.Most importantly, in both PCR and BR, the presence of Na and or Ni played a vital role for As remaining in solid phase at temperatures ≤800 • C. Eventually, these modellings are valuable for mitigating the negative impacts of toxicity emissions associated with CCA-elements; particularly, by systematically blending different types of contaminated WW prior to gasification.Author Contributions: All authors contributed to this work.S.B.A.: Writing-Original Draft Preparation; S.T.W. and Y.J.: Supervision and Guidance.

Table 1 .
The selected Cu, Cr, and As (CCA-elements) species.

Table 2 .
Statistical analysis results of the composition data.

Table 7 .
Phase transformations and speciation of the CCA-elements.