Microsoft Word - 202-1279 vetted IIUM Engineering Journal, Vol. 12, No. 6, 2011: Special Issue in Science and Ethics Basuki et al. 117 MATHEMATICAL MODELING AND NUMERICAL SOLUTION OF IRON CORROSION PROBLEM BASED ON CONDENSATION CHEMICAL PROPERTIES BASUKI W. 1 , ARIF F. 2 AND TRI R. 2 1 Department of Mathematics, Institut Teknologi Sepuluh Nopember,Surabaya, Indonesia. 2 Department of Mathematics & Industrial Technology, Jamber University-Wijaya Kusuma Surabaya University ,Jamber-Surabaya, Indonesia. b_widodo@matematika.its.ac.id ABSTRACT: Corrosion process is a natural phenomenin that occurs with the various metals, where the corrosion process in electrochemical environment can be explained by using a galvanic cell. The iron corrosion process is based on the acidity degree (pH) of a condensation, iron concentration and condensation temperature of electrolyte. These are applied at electrochemistry cell. The iron corrosion process at this electrochemical cell also is able to generate electrical potential and electric current during the process. This paper considers how to build a mathematical model of iron corrosion based on electrical potential and electric current. The mathematical model further is solved using the finite element method. The iron corrosion model is built based on the iron concentration, condensation temperature, and iteration time applied. In the electric current density model, the current is generated between the cathode and anode poles and the iteration time applied. Whereas the potential electric model, is based on the electric potential and the iteration time applied. The numerical results show that the part of iron sample, attacked by corrosion, is a function of time, depth difference, iron concentration, condensation temperature and the sum of reduced mass during corrosion process. ABSTRAK: Proses kakisan adalah kes semulajadi yang berlaku di pelbagai logam, di mana proses kakisan dalam elektrokimia boleh diterangkan dengan menggunakan sel galvani. Proses pengaratan besi adalah berdasarkan kepada tahap keasidan (pH) daripada pemeluapan, kepekatan besi dan suhu pemeluapan elektrolit. Sifat-sifat ini digunakan pada sel elektrokimia. Proses pengaratan besi pada sel elektrokimia ini juga mampu menjana arus dan voltan elektrik. Kertas ini menerangkan cara-cara membina model matematik untuk kakisan besi, arus dan voltan elektrik. Model matematik ini diselesaikan dengan menggunakan kaedah unsur terhingga. Model kakisan besi dibina berdasarkan kadar kepekatan besi, suhu pemeluapan, dan lelaran. Ketumpatan arus elektrik berdasarkan jumlah arus elektrik yang ada tiang-tiang katod dan anod tiang pada setiap lelaran. Manakala model potensi elektrik, berdasarkan pada kadar awal voltan elektrik pada masa lelaran. Bacaan experimen menunjukkan bahawa bahagian anod tetap keras walaupun dikakis dan dipengaruhi oleh kepekatan besi dan suhu pemeluwapan dan jumlah jisim. KEYWORDS: iron corrosion; concentration of iron; electrochemical cell; finite element method. IIUM Engineering Journal, Vol. 12, No. 6, 201 1. INTRODUCTION Iron is one of the most important m Corrosion of iron happened naturally. Rustiness can be described as a material ruination by action which embossing area, where resistance of a material on the rustiness depends on the numerous variables such [1]. The important aspect in a material rustiness is mass transport characteristic, that due to evacuation of temperature, migration and diffusion. The distribution character from the effects which are generated in an object depends on force sys itself object. A real rustiness case is a galvanic process with metal that works as an anode, a cathode and as one conductor for the electrons. The quantitative simulation of iron rustiness relates to possibility that deterioration forecasting of conglomerate stones are contain the iron item, thus here pointed to build a rustiness model as possible as free of context in physically and can be interfaced to a model of chemical reaction. This paper focuses on a mathematical modeling of an electro chemistries reaction that ta during corrosion process and the model then solved using finite element method, in which it is shown geometrically as an electrochemistry cell process. reduction-oxidation reaction in electrochemistry cell by using a pa based on the law of physical chemistry. 2. GOVERNING EQUATION The important aspect in a material rustiness is mass transport characteristics, which are occurred because of temperature, migration and diffusion evacuation. These are appropriate with condensation mass i.e.: Jk = –Dk Ck–C where : Dk = diffusion constant Ck = concentration in Molar (moles/liter) zk = charge number F = Faraday’s constant = the electric potential in the electrolyte R = gas constant T = absolute temperature This condensation mass further is transported by electrical field on molecular diffusion process [2]. For each type, the equation of transporta of accumulator acid, i.e.: + .[(– g Journal, Vol. 12, No. 6, 2011: Special Issue in Science and Ethics 118 INTRODUCTION Iron is one of the most important metals that has many benefits to the human iron happened naturally. Rustiness can be described as a material ruination by action which embossing area, where resistance of a material on the rustiness depends such as material properties and environmental characteris The important aspect in a material rustiness is mass transport characteristic, that evacuation of temperature, migration and diffusion. The distribution character from the effects which are generated in an object depends on force system characteristic and A real rustiness case is a galvanic process with metal that works as an anode, a cathode and as one conductor for the electrons. The quantitative simulation of iron rustiness relates to possibility that deterioration forecasting of conglomerate stones are contain the iron item, thus here pointed to build a rustiness model as possible as free of context in physically and can be interfaced to a model of chemical reaction. This paper focuses on a mathematical modeling of an electro chemistries reaction that ta during corrosion process and the model then solved using finite element method, in which it is shown geometrically as an electrochemistry cell process. This paper illustrates the oxidation reaction in electrochemistry cell by using a partial differential model based on the law of physical chemistry. GOVERNING EQUATION The important aspect in a material rustiness is mass transport characteristics, which are occurred because of temperature, migration and diffusion evacuation. These are ropriate with condensation mass Jk of species k, as mentioned in Planck- Ck (1) = concentration in Molar (moles/liter) = the electric potential in the electrolyte This condensation mass further is transported by electrical field on molecular diffusion process [2]. For each type, the equation of transportation can be written in a form ).Ck]–Dk 2 Ck=Sk (2) : Special Issue in Science and Ethics Basuki et al. human society. iron happened naturally. Rustiness can be described as a material ruination by action which embossing area, where resistance of a material on the rustiness depends as material properties and environmental characteristics The important aspect in a material rustiness is mass transport characteristic, that occur evacuation of temperature, migration and diffusion. The distribution character from tem characteristic and A real rustiness case is a galvanic process with metal that works as an anode, a cathode and as one conductor for the electrons. The quantitative simulation of iron rustiness relates to possibility that deterioration forecasting of conglomerate stones are contain the iron item, thus here pointed to build a rustiness model as possible as free of context in physically and can be interfaced to a model of chemical reaction. This paper focuses on a mathematical modeling of an electro chemistries reaction that taken place during corrosion process and the model then solved using finite element method, in which illustrates the rtial differential model The important aspect in a material rustiness is mass transport characteristics, which are occurred because of temperature, migration and diffusion evacuation. These are -Nernst’s law, This condensation mass further is transported by electrical field on molecular tion can be written in a form IIUM Engineering Journal, Vol. 12, No. 6, 201 where Sk is ions production or ruination of type Fe +2 ion because of rust forming. In this corrosion process case, it is taken place in the interface between electrodes and condensation of electrolyte. Electrical tension can be obtained by considering the total electrical current and it disregards the neutrality electric condensation of electrolyte. Newman (1996) explains that the power of the electrical current condensation mass constant change from iron. It can be formulated using the Faraday’s law: A statement of electro-neutrality is that the electric current density is non is: . i = 0 Inserting the Expression (3) in Equation (4), we obtain: . ( ) = 0 By using the Planck–Nernst’s law (1), we then end up with an equation for the electric potential, i.e.: Notice that when the concentration gradients may be neglected, Eq. (6) then reduces to the Laplace’s equation, i.e.: 2 = 0. 3. BOUNDARY CONDITIONS By analyzing the electrochemistry, iron rustiness proc oxidation by the oxygen that comes from the air. Figure electrolyses with (Fe) iron electrode. In this condensation there are some species, for example H + and OH - ions obtained from the result of H ion that come from electrode ionization. Fe OH - ions peripatetic towards positive pole [4]. (a) g Journal, Vol. 12, No. 6, 2011: Special Issue in Science and Ethics 119 is ions production or ruination of type k, for example in elimination case from ion because of rust forming. In this corrosion process case, it is taken place in the interface between electrodes and condensation of electrolyte. Electrical tension can be obtained by considering the total electrical current and it disregards the neutrality electric condensation of electrolyte. Newman (1996) explains that the power of the electrical current I ,in condensation of the electrolyte, is condensation mass constant change from iron. It can be formulated using the Faraday’s neutrality is that the electric current density is non-divergent, that (4) pression (3) in Equation (4), we obtain: ) = 0 (5) Nernst’s law (1), we then end up with an equation for the electric = – ∑ Dk 2 Ck (6) Notice that when the concentration gradients may be neglected, Eq. (6) then reduces to the (7) BOUNDARY CONDITIONS By analyzing the electrochemistry, iron rustiness process is event of iron metal at comes from the air. Figure 1 shows that H2O condensation electrolyses with (Fe) iron electrode. In this condensation there are some species, for ions obtained from the result of H2O ionization as solvent and Fe ion that come from electrode ionization. Fe 2+ ions peripatetic towards negative pole and ions peripatetic towards positive pole [4]. (b) : Special Issue in Science and Ethics Basuki et al. , for example in elimination case from ion because of rust forming. In this corrosion process case, it is taken place in the Electrical tension can be obtained by considering the total electrical current and it disregards the neutrality electric condensation of electrolyte. Newman (1996) explains that ,in condensation of the electrolyte, is flown in the condensation mass constant change from iron. It can be formulated using the Faraday’s divergent, that Nernst’s law (1), we then end up with an equation for the electric Notice that when the concentration gradients may be neglected, Eq. (6) then reduces to the ess is event of iron metal O condensation electrolyses with (Fe) iron electrode. In this condensation there are some species, for ionization as solvent and Fe 2+ ions peripatetic towards negative pole and IIUM Engineering Journal, Vol. 12, No. 6, 201 Fig. 1: (a) Location of the anode and cath During the corrosion process takes place, it will entangle the reaction of reduction oxidation (redox) in electrolytic cell by using a partial di physical chemistry. Newman (1996) explains that the only reaction takes place in anode is disconnection of metal. In this case, reaction of (Fe) iron can be written as: (A)Fe → Fe 2+ where E 0 A is the standard potential of (A) reaction. The current density is produced by (A) reaction and it can be written by the Butler iA = i0A . where : i0A = exchange current density of (A) reaction, th which the net current at the electrode is null (this value is obtained experimentally), άA = transfer coefficient for (A) reaction (this value is obtained experimentally, although a good estimate is 1/2), zA = number of electrons participating in (A) reaction, ήA = overvoltage = äö − EA electrode and the electrolyte (soon after the electronic double layer) and E electrode potential at zero current. As only the species Fe 2+ mass flux will be non-zero only. Hence, once the density current is evaluated from Eq. (2), we can use the Faraday’s law to ev The fluxes are then used as boundary condition at the anode for the transport equation. Further, the cathode pole, between the easiest matters experiences reduction is molecule so that happened reduction to H g Journal, Vol. 12, No. 6, 2011: Special Issue in Science and Ethics 120 Location of the anode and cathode in a corrosion process (b) and a one dimensional galvanic cell. During the corrosion process takes place, it will entangle the reaction of reduction oxidation (redox) in electrolytic cell by using a partial differential model, based on law of physical chemistry. Newman (1996) explains that the only reaction takes place in anode is disconnection of metal. In this case, reaction of (Fe) iron can be written as: + 2e-, E 0 A = - 0,44 (oxidation) (8) is the standard potential of (A) reaction. The current density is produced by (A) reaction and it can be written by the Butler-Volmer’s equation as: (9) exchange current density of (A) reaction, that is the current density evaluated , in which the net current at the electrode is null (this value is obtained experimentally), transfer coefficient for (A) reaction (this value is obtained experimentally, although a good estimate is 1/2), number of electrons participating in (A) reaction, A, where äö is the variation in potential between the electrode and the electrolyte (soon after the electronic double layer) and E electrode potential at zero current. 2+ are involved in the anode reaction, we can assume that its zero only. Hence, once the density current is evaluated from Eq. (2), we can use the Faraday’s law to evaluate the mass flux: Jfe2 + = JH+ = 0 ; on x = 0. (10) The fluxes are then used as boundary condition at the anode for the transport equation. Further, the cathode pole, between the easiest matters experiences reduction is molecule so that happened reduction to H2O: : Special Issue in Science and Ethics Basuki et al. ) and a one- During the corrosion process takes place, it will entangle the reaction of reduction- fferential model, based on law of physical chemistry. Newman (1996) explains that the only reaction takes place in anode is is the standard potential of (A) reaction. The current density is produced by (A) at is the current density evaluated , in which the net current at the electrode is null (this value is obtained experimentally), transfer coefficient for (A) reaction (this value is obtained experimentally, , where äö is the variation in potential between the electrode and the electrolyte (soon after the electronic double layer) and EA is the are involved in the anode reaction, we can assume that its zero only. Hence, once the density current is evaluated from Eq. (2), The fluxes are then used as boundary condition at the anode for the transport equation. Further, the cathode pole, between the easiest matters experiences reduction is H2O IIUM Engineering Journal, Vol. 12, No. 6, 201 (B) 4H + + O 2 +4e The cathodes current density can once again be evaluated according to the Butler Volmer’s equation, by considering the (B) reaction only: iB = - i0B where : i0B = exchange current density of (B) reaction, άB = transfer coefficient for (B) reaction, zB = number of electrons participating in (B) reaction, ήB = over voltage, with the negative sign it means that the flow of current proceeds from the electrode to the electrolyte, on the reverse of direction with respect to the axis. The boundary conditions for the Transport Equations (2) are obtained as before, by considering that only the H+ species are involved in the cathode reaction, and therefore only its mass flux will be non- 4. INITIAL CONDITIONS Every condensation has degree of acidity (hydrogen ion exponent) influences to the corrosion process of a metal. Svante Arhenius comments that acid is a compound, which if it is dissolved in water then it will yield (H alkaline is a compound dissolved in water will yield (OH small fraction from (H2O) water molecules are dissociated in H amount are determined by the constant equilibrium from dissociation reaction , that is: H2O H Concentration of (H + ) Hydronium ion and (OH condensation generally very low but hardly determines characters of condensation, especially condensation in water. Sorensen (1868 concept to avoid usage of a real sm the function of negative logarithm from the concentration of H condensation, and formulated as follows: pH = – log [CH with CH + is concentration from ion H equilibrium of pure water, applies: g Journal, Vol. 12, No. 6, 2011: Special Issue in Science and Ethics 121 +4e - → 2H2O; E 0 B = 1,23 V. (11) The cathodes current density can once again be evaluated according to the Butler Volmer’s equation, by considering the (B) reaction only: , (12) exchange current density of (B) reaction, transfer coefficient for (B) reaction, number of electrons participating in (B) reaction, over voltage, with the negative sign it means that the flow of current proceeds from the electrode to the electrolyte, on the reverse of direction with respect to the axis. The boundary conditions for the Transport Equations (2) are obtained as before, by hat only the H+ species are involved in the cathode reaction, and therefore -zero: JFe2+ = 0 JH + = ; on x = L (13) INITIAL CONDITIONS Every condensation has degree of acidity (hydrogen ion exponent) influences to the corrosion process of a metal. Svante Arhenius comments that acid is a compound, which if it is dissolved in water then it will yield (H + ) Hidronium ion. While alkaline is a compound dissolved in water will yield (OH - ) Hydroxide ion [4]. As known, a O) water molecules are dissociated in H + and OH amount are determined by the constant equilibrium from dissociation reaction , that is: O H + + OH (14) ) Hydronium ion and (OH - ) Hydroxide ion in a liquid condensation generally very low but hardly determines characters of condensation, especially condensation in water. Sorensen (1868-1939) further proposes pH and pOH concept to avoid usage of a real small number [4]. According to Sorensen pH and pOH are the function of negative logarithm from the concentration of H + and OH condensation, and formulated as follows: H + ] dan pOH = – log [COH - ] (15) is concentration from ion H + and COH - is concentration of OH - equilibrium of pure water, applies: : Special Issue in Science and Ethics Basuki et al. The cathodes current density can once again be evaluated according to the Butler– with the negative sign it means that the flow of current proceeds from the electrode to the The boundary conditions for the Transport Equations (2) are obtained as before, by hat only the H+ species are involved in the cathode reaction, and therefore different and influences to the corrosion process of a metal. Svante Arhenius comments that acid is a ) Hidronium ion. While xide ion [4]. As known, a and OH - ions in the amount are determined by the constant equilibrium from dissociation reaction , that is: ) Hydroxide ion in a liquid condensation generally very low but hardly determines characters of condensation, 1939) further proposes pH and pOH all number [4]. According to Sorensen pH and pOH are and OH - ions in a ion where at IIUM Engineering Journal, Vol. 12, No. 6, 201 pH + pOH = 14 Those concentration of the ions expressed in a set of Molar (M), where the molarity expresses that the various solu has pH = 7, so that concentration from both of good ions of H Molar. Smaller pH value of condensation hence level of the solution acidity excelsior and so do on the contrary. 5. RESULTS AND DISCUSSION Finite element method is one of approach method of numeric that bases on the problems at every part of element that is named by finite element. Every problem that exists will be finalized with the quadratic approach, where element method has form of matrix equation: where: [K],[Kt] = assemble matrix {r} = vector from node magnitude which [R] = assemble style parameter vector Global matrix equation hereinafter will be finalized to time, where in this solution it will be used approach of finite difference with the Crack global matrix compiler can be w If it is taken è = 1/2 hence will be obtained Crack pattern: Global matrix equation therefore will have a form of approach to time as follows: g Journal, Vol. 12, No. 6, 2011: Special Issue in Science and Ethics 122 + pOH = 14 (16) Those concentration of the ions expressed in a set of Molar (M), where the molarity expresses that the various solute moles in every 1 condensation liter ( mol/liter). Pure water has pH = 7, so that concentration from both of good ions of H + and OH - ions equal to 10 Molar. Smaller pH value of condensation hence level of the solution acidity excelsior and RESULTS AND DISCUSSION Finite element method is one of approach method of numeric that bases on the problems at every part of element that is named by finite element. Every problem that exists will be finalized with the quadratic approach, where the problem solving form finite element method has form of matrix equation: (17) where: = vector from node magnitude which unknown = assemble style parameter vector Global matrix equation hereinafter will be finalized to time, where in this solution it will be used approach of finite difference with the Crack-Nicholson’s pattern. Equation of global matrix compiler can be written down in the form of finite difference as: If it is taken è = 1/2 hence will be obtained Crack-Nicholson’s pattern, with approach Global matrix equation therefore will have a form of approach to time as follows: : Special Issue in Science and Ethics Basuki et al. Those concentration of the ions expressed in a set of Molar (M), where the molarity te moles in every 1 condensation liter ( mol/liter). Pure water ions equal to 10 -7 Molar. Smaller pH value of condensation hence level of the solution acidity excelsior and Finite element method is one of approach method of numeric that bases on the problems at every part of element that is named by finite element. Every problem that the problem solving form finite Global matrix equation hereinafter will be finalized to time, where in this solution it Nicholson’s pattern. Equation of ritten down in the form of finite difference as: Nicholson’s pattern, with approach Global matrix equation therefore will have a form of approach to time as follows: IIUM Engineering Journal, Vol. 12, No. 6, 2011: Special Issue in Science and Ethics Basuki et al. 123 The numerical solution of those equations which will be visualized by using MATLAB 7.0 program. The result is shown in Fig. 2: (a) (b) (c) IIUM Engineering Journal, Vol. 12, No. 6, 2011: Special Issue in Science and Ethics Basuki et al. 124 Fig. 2: Iron Corrosion Process based on difference Time (a), Concentration of Iron (b), and Temperature condensation (c). Based on the obtained figure above, it can be determined that the relation of the longer time required. Therefore, the chances that the happening of corrosion is substantially increase, and so do on the contrary. The concentration influence of iron towards the iron corrosion process is greater than the applied iron concentration. Therefore, the possibility that the happening of iron corrosion increasingly increases and conversely. Third factor that has an effect on is condensation temperature used, ever greater of temperature used hence possibility that the happening of corrosion exactly increasingly minimizes. The statement chimes in with Planck-Nernst law about mass transport and First Faraday’s law where mass that’s yielded or reduced compares straight to time and concentration of metal that’s applied but proportional inversely with condensation temperature. Numerical result of the electrical potential equation is also be visualized by using MATLAB 7.0 program and it is depicted in Fig. 3: (a) (b) IIUM Engineering Journal, Vol. 12, No. 6, 2011: Special Issue in Science and Ethics Basuki et al. 125 Fig. 3: Electric potential based to time difference (a) and electric potential beginning (b). Based on the above picture, i.e. Fig. 3, which obtained for determinable electrical tension the relation of that longer time required hence level of electrical tension value of mean emerging at every node would increasingly increases, and so do on the contrary. The beginning electric potential influence that’s used smaller that’s applied hence electric potential average value also increasingly minimizes, and so does on the contrary. The statement chimes in with Planck-Nernst’s law about mass transport, where potential electric is compared straight to time. The electrical current also arises during the corrosion process at electrochemistry cell, in which the numerical results of those equations will be visualized in the form of figures by using the MATLAB 7.0 program. The results are shown in Fig. 4: (a) (b) IIUM Engineering Journal, Vol. 12, No. 6, 2011: Special Issue in Science and Ethics Basuki et al. 126 Fig. 4: Electric current based to time difference (a) and condensation temperature (b). Based on the Fig. 4 that are obtained is known that longer time required hence electric current that emerges at every node increasingly declines, and so do on the contrary. Temperature factor also influences to electrical current value that emerges during corrosion process to take place that is ever greater of condensation temperature used hence electric current happened at every node also increasingly increased where both of those relationships chimed in with Faraday’s first law and mass transport Planck-Nernst’s law that is inversely proportional electric current to the time and it compares straight with the condensation temperature. 6. CONCLUSION Based on the numerical simulations which have made so it can be concluded that the inferential part of gristle metal part of functioning metal as anode support, in this case there is in the beginning of metal. If the time is longer applied hence part of metal experiencing impairment concentration by the end of iteration would increasingly long in meaning of gristle area corrosion to be longer so that metal mass that reduce by the corrosion more and more and so do on the contrary. The concentration factor of initial iron also influences to corrosion process, that is value excelsior concentration of iron used hence impairment concentration of iron by the end of iteration also increasingly long in meaning of area which gristle corrosion to be broader, so that number of masses that’s yielded for the concentration of larger ones has number of bigger reduced masses and so do on the contrary. Temperature factor also influences to iron corrosion that is more high condensation temperature excelsior used in condensation of electrolyte hence degradation of concentration of increasingly short, in meaning of gristle area corrosion to be more minimizes, so that mass reduced that’s yielded to increasingly minimizes and so do on the contrary. At the electrical potential increasingly of time hence the potential electrical average value also increases, while ever greater of beginning electric potential that’s applied hence electric potential value which emerges at every node and electric potential average value also increases, and so do on the contrary. At electric current increasingly increasing of time hence electric current value which emerges that’s increasingly minimizes, while IIUM Engineering Journal, Vol. 12, No. 6, 2011: Special Issue in Science and Ethics Basuki et al. 127 growing of condensation temperature of electrolyte used hence electric current value which emerges to increasingly, and so do on the contrary. ACKNOWLEDGEMENT We would like to thank the Institut Teknologi Sepuluh Nopember-Surabaya Indonesia, which give us a grant and chance to disseminate an intellectual output of our research. REFERENCES [1] Shifler D. A, “Corrosion Science”, Physical Science, Vol. 47, 2005, pp. 2335–2352. [2] Atkins P.W., Physical Chemistry, Oxford University Press, Oxford, 1990 [3] Newman, J. S., Electrochemical Systems, Prentice-Hall, Englewood Cliffs, New Jersey, 1996. [4] Sudarmo U., Chemistry, Erlangga, Jakarta, 2006.