Analyzing and Modelling the Corrosion Behavior of Ni/Al2O3, Ni/SiC, Ni/ZrO2 and Ni/Graphene Nanocomposite Coatings

A study has been presented on the effects of intrinsic mechanical parameters, such as surface stress, surface elastic modulus, surface porosity, permeability and grain size on the corrosion failure of nanocomposite coatings. A set of mechano-electrochemical equations was developed by combining the popular Butler–Volmer and Duhem expressions to analyze the direct influence of mechanical parameters on the electrochemical reactions in nanocomposite coatings. Nanocomposite coatings of Ni with Al2O3, SiC, ZrO2 and Graphene nanoparticles were studied as examples. The predictions showed that the corrosion rate of the nanocoatings increased with increasing grain size due to increase in surface stress, surface porosity and permeability of nanocoatings. A detailed experimental study was performed in which the nanocomposite coatings were subjected to an accelerated corrosion testing. The experimental results helped to develop and validate the equations by qualitative comparison between the experimental and predicted results showing good agreement between the two.


Introduction
This paper enables researchers to fully understand the complex relationship between the mechanical properties and macroscopic phenomena, including electrochemical and chemical responses of nanocoatings. It is well known that nanomaterial has large surface-to-volume ratio compared to bulk material, however, in the case of nanocoating this role becomes even more significant [1]. A literature survey shows that the classical electrochemical theory does not fully incorporate the influences of mechanical parameters of nanocoatings [2,3].
It has been found in the literature [4,5] that the combined effects of both the mechanical and electrochemical parameters directly influence the electrochemical reactions on the surface of nanocoating. Khan-Nazir [6,7] suggested that electrochemical reactions on metal coating surfaces are stress-dependent. Sometimes, both mechanical and electrochemical effects, coupled together, result in the modified corrosion rates of the nanocoating compared to its stress-free state [8][9][10]. Therefore, it can be established that the size-dependent electrochemical properties of nanocoatings (31 g/L). Four types of nanoparticles with an amount of 20 g/L were ultrasonically dispersed in different solutions. The sizes of nano-particles were: Nano Al 2 O 3 (50-60 nm), Nano SiC (40-50 nm), Nano ZrO 2 (30-40 nm) and Graphene platelets (6-8 nm). The solutions had been magnetically stirred overnight to yield better particle suspension prior to the start of the deposition process. The pH of the solution was adjusted to between 4.0 and 4.5 by sodium hydroxide or diluted sulphuric acid and recorded using a Tecpel pH meter.
The electrodeposition process parameters were kept constant as current density (3 A/dm 2 ) and pulse on-off time (20-80 ms) with a duty cycle of 20%. A high-quality nickel sheet was used as an anode, and a steel circular disc of 80 mm diameter and 8.20 mm thick was used as a cathode.

Accelerated Corrosion Testing
All the coated samples were exposed to salt spray testing for 450 h. Salt spray testing was performed by placing samples in a test chamber which was operated according to the specified conditions in ASTM B117. This test was performed to analyze the relative resistance to corrosion of coated samples when exposed to a salt spray climate at elevated temperature. The chamber was enclosed with a well-maintained temperature of +35 • C, and a continuous indirect spray of 5% aqueous sodium chloride solution was allowed on to the samples at a specified rate.
The appearance of corrosion products was evaluated after the exposure test. Surface characterization of test samples, including crystal orientation and grain size, was performed both pre-exposure and post-exposure using SEM (JEOL USA Inc., Peabody, MA, USA), EDS (JEOL USA Inc., Peabody, MA, USA) and XRD (BRUKER, ,Petaling Jaya, Malaysia). Figure 1 shows the corrosion status of electrodeposited composite coatings after 450 h of salt spray exposure. The images clearly show that the highest surface corrosion was found for the case of Ni/Al 2 O 3 coating which subsequently decreased for the case of Ni/SiC and Ni/ZrO 2 coatings. The lowest surface corrosion was found for the case of Ni/GPL coating. According to the latest studies [29], the first five atomic surface layers of nano coatings have the highest degradation rates because of direct exposure to the environment. Therefore, prominent changes in the surface morphology of post-test samples can be seen from the images.

Surface Corrosion and Porosity
Materials 2017, 10, 1225 3 of 16 overnight to yield better particle suspension prior to the start of the deposition process. The pH of the solution was adjusted to between 4.0 and 4.5 by sodium hydroxide or diluted sulphuric acid and recorded using a Tecpel pH meter. The electrodeposition process parameters were kept constant as current density (3 A/dm 2 ) and pulse on-off time (20-80 ms) with a duty cycle of 20%. A high-quality nickel sheet was used as an anode, and a steel circular disc of 80 mm diameter and 8.20 mm thick was used as a cathode.

Accelerated Corrosion Testing
All the coated samples were exposed to salt spray testing for 450 h. Salt spray testing was performed by placing samples in a test chamber which was operated according to the specified conditions in ASTM B117. This test was performed to analyze the relative resistance to corrosion of coated samples when exposed to a salt spray climate at elevated temperature. The chamber was enclosed with a well-maintained temperature of +35 °C, and a continuous indirect spray of 5% aqueous sodium chloride solution was allowed on to the samples at a specified rate.
The appearance of corrosion products was evaluated after the exposure test. Surface characterization of test samples, including crystal orientation and grain size, was performed both pre-exposure and post-exposure using SEM (JEOL USA Inc., Peabody, MA, USA), EDS (JEOL USA Inc., Peabody, MA, USA) and XRD (BRUKER, ,Petaling Jaya, Malaysia). Figure 1 shows the corrosion status of electrodeposited composite coatings after 450 h of salt spray exposure. The images clearly show that the highest surface corrosion was found for the case of Ni/Al2O3 coating which subsequently decreased for the case of Ni/SiC and Ni/ZrO2 coatings. The lowest surface corrosion was found for the case of Ni/GPL coating. According to the latest studies [29], the first five atomic surface layers of nano coatings have the highest degradation rates because of direct exposure to the environment. Therefore, prominent changes in the surface morphology of post-test samples can be seen from the images.  Table 1, the pre-exposure SEM micrographs for all the composite coatings demonstrate compacted fine-deposit grain microstructure while the pre-exposure EDS micrographs for the coatings identify the elements and their concentrations. The SEM results showed that Ni/Al2O3 nano-  Table 1, the pre-exposure SEM micrographs for all the composite coatings demonstrate compacted fine-deposit grain microstructure while the pre-exposure EDS micrographs for the coatings identify the elements and their concentrations. The SEM results showed that Ni/Al 2 O 3 nano-sized particles revealed bigger deposit grain size when compared to the remainder coatings. This also contributed to an increase in the average surface roughness of Ni/Al 2 O 3 coating with a value of 0.21 µm followed by Ni/SiC (0.19 µm), Ni/ZrO 2 (0.06 µm) and Ni/GPL (0.007 µm) coatings respectively. The reason for large surface roughness behavior can be attributed to the existence of bigger particles on the surface. With the decrease in the particle size from Ni/Al 2 O 3 to Ni/GPL, the surface roughness of coating decreased. Also, large particle size in Ni/Al 2 O 3 coatings accounted for large porous structure. The decreasing particle size from Ni/Al 2 O 3 to Ni/GPL coating resulted in the formation of fine coating with minimum porosity. For example, Ni/GPL coating which had the smallest particle size formed a very fine glass-like coating with minimum porosity.

Surface Corrosion and Porosity
The post exposure SEM and EDS results revealed highest corrosion products (Fe = 3.0%) appearance in Ni/Al 2 O 3 and smallest corrosion product (Fe = 0.3%) appearance in Ni/GPL coating. This behavior can be directly linked to the porosity level which sets the Oxygen and NaCl electrolyte permeability of coatings. Post exposure micrographs showed that Oxygen and electrolyte permeability of coatings decreased with a decrease in the particle size. The inhibition efficiency was found to increase with decreasing the size of particles showing that Ni/GPL coating had the least corrosion found after exposure. sized particles revealed bigger deposit grain size when compared to the remainder coatings. This also contributed to an increase in the average surface roughness of Ni/Al2O3 coating with a value of 0.21 µm followed by Ni/SiC (0.19 µm), Ni/ZrO2 (0.06 µm) and Ni/GPL (0.007 µm) coatings respectively. The reason for large surface roughness behavior can be attributed to the existence of bigger particles on the surface. With the decrease in the particle size from Ni/Al2O3 to Ni/GPL, the surface roughness of coating decreased. Also, large particle size in Ni/Al2O3 coatings accounted for large porous structure. The decreasing particle size from Ni/Al2O3 to Ni/GPL coating resulted in the formation of fine coating with minimum porosity. For example, Ni/GPL coating which had the smallest particle size formed a very fine glass-like coating with minimum porosity. The post exposure SEM and EDS results revealed highest corrosion products (Fe = 3.0%) appearance in Ni/Al2O3 and smallest corrosion product (Fe = 0.3%) appearance in Ni/GPL coating. This behavior can be directly linked to the porosity level which sets the Oxygen and NaCl electrolyte permeability of coatings. Post exposure micrographs showed that Oxygen and electrolyte permeability of coatings decreased with a decrease in the particle size. The inhibition efficiency was found to increase with decreasing the size of particles showing that Ni/GPL coating had the least corrosion found after exposure.  sized particles revealed bigger deposit grain size when compared to the remainder coatings. This also contributed to an increase in the average surface roughness of Ni/Al2O3 coating with a value of 0.21 µm followed by Ni/SiC (0.19 µm), Ni/ZrO2 (0.06 µm) and Ni/GPL (0.007 µm) coatings respectively. The reason for large surface roughness behavior can be attributed to the existence of bigger particles on the surface. With the decrease in the particle size from Ni/Al2O3 to Ni/GPL, the surface roughness of coating decreased. Also, large particle size in Ni/Al2O3 coatings accounted for large porous structure. The decreasing particle size from Ni/Al2O3 to Ni/GPL coating resulted in the formation of fine coating with minimum porosity. For example, Ni/GPL coating which had the smallest particle size formed a very fine glass-like coating with minimum porosity. The post exposure SEM and EDS results revealed highest corrosion products (Fe = 3.0%) appearance in Ni/Al2O3 and smallest corrosion product (Fe = 0.3%) appearance in Ni/GPL coating. This behavior can be directly linked to the porosity level which sets the Oxygen and NaCl electrolyte permeability of coatings. Post exposure micrographs showed that Oxygen and electrolyte permeability of coatings decreased with a decrease in the particle size. The inhibition efficiency was found to increase with decreasing the size of particles showing that Ni/GPL coating had the least corrosion found after exposure.   sized particles revealed bigger deposit grain size when compared to the remainder coatings. This also contributed to an increase in the average surface roughness of Ni/Al2O3 coating with a value of 0.21 µm followed by Ni/SiC (0.19 µm), Ni/ZrO2 (0.06 µm) and Ni/GPL (0.007 µm) coatings respectively. The reason for large surface roughness behavior can be attributed to the existence of bigger particles on the surface. With the decrease in the particle size from Ni/Al2O3 to Ni/GPL, the surface roughness of coating decreased. Also, large particle size in Ni/Al2O3 coatings accounted for large porous structure. The decreasing particle size from Ni/Al2O3 to Ni/GPL coating resulted in the formation of fine coating with minimum porosity. For example, Ni/GPL coating which had the smallest particle size formed a very fine glass-like coating with minimum porosity. The post exposure SEM and EDS results revealed highest corrosion products (Fe = 3.0%) appearance in Ni/Al2O3 and smallest corrosion product (Fe = 0.3%) appearance in Ni/GPL coating. This behavior can be directly linked to the porosity level which sets the Oxygen and NaCl electrolyte permeability of coatings. Post exposure micrographs showed that Oxygen and electrolyte permeability of coatings decreased with a decrease in the particle size. The inhibition efficiency was found to increase with decreasing the size of particles showing that Ni/GPL coating had the least corrosion found after exposure.  sized particles revealed bigger deposit grain size when compared to the remainder coatings. This also contributed to an increase in the average surface roughness of Ni/Al2O3 coating with a value of 0.21 µm followed by Ni/SiC (0.19 µm), Ni/ZrO2 (0.06 µm) and Ni/GPL (0.007 µm) coatings respectively. The reason for large surface roughness behavior can be attributed to the existence of bigger particles on the surface. With the decrease in the particle size from Ni/Al2O3 to Ni/GPL, the surface roughness of coating decreased. Also, large particle size in Ni/Al2O3 coatings accounted for large porous structure. The decreasing particle size from Ni/Al2O3 to Ni/GPL coating resulted in the formation of fine coating with minimum porosity. For example, Ni/GPL coating which had the smallest particle size formed a very fine glass-like coating with minimum porosity. The post exposure SEM and EDS results revealed highest corrosion products (Fe = 3.0%) appearance in Ni/Al2O3 and smallest corrosion product (Fe = 0.3%) appearance in Ni/GPL coating. This behavior can be directly linked to the porosity level which sets the Oxygen and NaCl electrolyte permeability of coatings. Post exposure micrographs showed that Oxygen and electrolyte permeability of coatings decreased with a decrease in the particle size. The inhibition efficiency was found to increase with decreasing the size of particles showing that Ni/GPL coating had the least corrosion found after exposure.   420) and (422) crystal planes, respectively. There were no significant changes recorded in any of the spectra, however, intensity of the peaks (especially (111) and (200) peaks) after the immersion reduced for all the samples in order to compensate the externally diffused NaCl crystal peaks. This behavior is attributed to the increase in porosity and reduction in the inhibition efficiency of coatings. For the case of Ni/GPL, minimum drop in the intensity peaks in addition to the few peaks for NaCl was observed compared to the other coatings. There were no significant changes recorded in any of the spectra, however, intensity of the peaks (especially (111) and (200) peaks) after the immersion reduced for all the samples in order to compensate the externally diffused NaCl crystal peaks. This behavior is attributed to the increase in porosity and reduction in the inhibition efficiency of coatings. For the case of Ni/GPL, minimum drop in the intensity peaks in addition to the few peaks for NaCl was observed compared to the other coatings. Further investigation showed that the highest % rise in the grain size after exposure was found in the case of Ni/Al2O3 coating, while the lowest % rise was observed in case of Ni/GPL coating, as shown in Figure 3. The large % rise in Ni/Al2O3 occurred due to its open, disordered nature of grain boundaries, meaning that NaCl crystals diffused more rapidly down boundaries leading to more rapid coble creep [30]. Since grain boundaries are regions of high energy, they make excellent sites for the nucleation of salt precipitates. Contrarily, Ni/GPL coating had compact, refined grain structure leading to least diffusion of NaCl crystals.  Further investigation showed that the highest % rise in the grain size after exposure was found in the case of Ni/Al 2 O 3 coating, while the lowest % rise was observed in case of Ni/GPL coating, as shown in Figure 3. The large % rise in Ni/Al 2 O 3 occurred due to its open, disordered nature of grain boundaries, meaning that NaCl crystals diffused more rapidly down boundaries leading to more rapid coble creep [30]. Since grain boundaries are regions of high energy, they make excellent sites for the nucleation of salt precipitates. Contrarily, Ni/GPL coating had compact, refined grain structure leading to least diffusion of NaCl crystals.

Surface Stress and Corrosion Rate
The residual stresses in all the coatings were measured by XRD and Reuss model using the computer program, Stress [16,31,32], while the corrosion rate measurement was performed by EIS measurements using potentiostat setup. From Figure 4 it was established that the nano particle size affects both the residual stresses and the corrosion rate of nickel-based composite coatings. The larger particles of Ni/Al2O3 impeded the lateral growth of the nickel layer, producing a weak (111) plane compared to the rest of the coatings, which resulted in high tensile residual stresses and large porosity. Large tensile stress with high porosity led to the large openings of micro pores in Ni/Al2O3, providing easy pathways for NaCl ions to diffuse through the pores. This behavior, commonly known as 'diffusion creep' resulted in a high intergranular corrosion rate of Ni/Al2O3. Similar tensile behavior was also found in the case of Ni/SiC and Ni/ZrO2, but rather in a decreasing order. On the other hand, the finer particles of Ni/GPL resulted in compressive stresses and minimum porosity, ultimately leading to the minimum corrosion rate.

Surface Stress and Corrosion Rate
The residual stresses in all the coatings were measured by XRD and Reuss model using the computer program, Stress [16,31,32], while the corrosion rate measurement was performed by EIS measurements using potentiostat setup. From Figure 4 it was established that the nano particle size affects both the residual stresses and the corrosion rate of nickel-based composite coatings. The larger particles of Ni/Al 2 O 3 impeded the lateral growth of the nickel layer, producing a weak (111) plane compared to the rest of the coatings, which resulted in high tensile residual stresses and large porosity. Large tensile stress with high porosity led to the large openings of micro pores in Ni/Al 2 O 3 , providing easy pathways for NaCl ions to diffuse through the pores. This behavior, commonly known as 'diffusion creep' resulted in a high intergranular corrosion rate of Ni/Al 2 O 3 . Similar tensile behavior was also found in the case of Ni/SiC and Ni/ZrO 2 , but rather in a decreasing order. On the other hand, the finer particles of Ni/GPL resulted in compressive stresses and minimum porosity, ultimately leading to the minimum corrosion rate.

Surface Stress and Corrosion Rate
The residual stresses in all the coatings were measured by XRD and Reuss model using the computer program, Stress [16,31,32], while the corrosion rate measurement was performed by EIS measurements using potentiostat setup. From Figure 4 it was established that the nano particle size affects both the residual stresses and the corrosion rate of nickel-based composite coatings. The larger particles of Ni/Al2O3 impeded the lateral growth of the nickel layer, producing a weak (111) plane compared to the rest of the coatings, which resulted in high tensile residual stresses and large porosity. Large tensile stress with high porosity led to the large openings of micro pores in Ni/Al2O3, providing easy pathways for NaCl ions to diffuse through the pores. This behavior, commonly known as 'diffusion creep' resulted in a high intergranular corrosion rate of Ni/Al2O3. Similar tensile behavior was also found in the case of Ni/SiC and Ni/ZrO2, but rather in a decreasing order. On the other hand, the finer particles of Ni/GPL resulted in compressive stresses and minimum porosity, ultimately leading to the minimum corrosion rate.

Mathematical Model
A set of mechano-electrochemical equations was developed for the nanocomposite coating. The nanocomposite coating model consists of the 3D surface layers adhered to a 3D core ( Figure 5). It is assumed that the surface layers and the core are linearly elastic and mechanically isotropic in addition to homogenous stresses. The <100>, <010> and <001> lattice directions represent the x, y and z coordinates respectively. This model develops, in a step-wise manner, the thermodynamic equations for: • Mechanical parameters including surface stress, porosity, permeability, particle size and grain size of the composite coating. • Corrosion rate of the composite coating.

Mathematical Model
A set of mechano-electrochemical equations was developed for the nanocomposite coating. The nanocomposite coating model consists of the 3D surface layers adhered to a 3D core ( Figure 5). It is assumed that the surface layers and the core are linearly elastic and mechanically isotropic in addition to homogenous stresses. The <100>, <010> and <001> lattice directions represent the x, y and z coordinates respectively. This model develops, in a step-wise manner, the thermodynamic equations for: • Mechanical parameters including surface stress, porosity, permeability, particle size and grain size of the composite coating.

•
Corrosion rate of the composite coating. The focus of study is on the elastic properties of a composite coating with typical threedimensional (3D) structure to simplify the theoretical analysis. In the study of surface thermodynamic properties of composite materials, the 3D surface layer is assumed to be 1 nm thick. This assumption follows the fact that during the relaxation process of nanocomposite coating, only the first five atomic surface layers (shown as yellow in Figure 5) encounter normal relaxation, while the rest of the atomic layers remain undeformed [29,33]. After this, the surface thickness is assumed to be 1 nm to model the size-dependent properties of the coating.

Mechanical Parameters Modelling
Usually, 3D surfaces are modelled using the diffusive and interphase surface approaches [34]. The diffusive interface represents the concentration gradient of electrolyte [35]. The interface is treated as a thermodynamic phase by the interphase approach with no significant variation of properties along the interface. As the atoms within the 1 nm thin surface layer face a very distinct local environment (e.g., corrosive) compared to the core, the mechanical response and physical properties of the surface will be different from those of the core.
Consider a composite coating which has debonded from the stress-free substrate; the coating will reach its minimum energy after relaxing and attaining equilibrium. There are two types of relaxation processes, including normal and parallel relaxation [9,36]. The normal relaxation leads to the creation of an eigenstress σ in the surface layer, while the normal relaxation leads to creation of a stress (or initial strain ε ) in the core to balance the surface stress. As the surface is bonded to the core, both the surface and core should undergo the same magnitude of deformation. The focus of study is on the elastic properties of a composite coating with typical three-dimensional (3D) structure to simplify the theoretical analysis. In the study of surface thermodynamic properties of composite materials, the 3D surface layer is assumed to be 1 nm thick. This assumption follows the fact that during the relaxation process of nanocomposite coating, only the first five atomic surface layers (shown as yellow in Figure 5) encounter normal relaxation, while the rest of the atomic layers remain undeformed [29,33]. After this, the surface thickness is assumed to be 1 nm to model the size-dependent properties of the coating.

Mechanical Parameters Modelling
Usually, 3D surfaces are modelled using the diffusive and interphase surface approaches [34]. The diffusive interface represents the concentration gradient of electrolyte [35]. The interface is treated as a thermodynamic phase by the interphase approach with no significant variation of properties along the interface. As the atoms within the 1 nm thin surface layer face a very distinct local environment (e.g., corrosive) compared to the core, the mechanical response and physical properties of the surface will be different from those of the core.
Consider a composite coating which has debonded from the stress-free substrate; the coating will reach its minimum energy after relaxing and attaining equilibrium. There are two types of relaxation processes, including normal and parallel relaxation [9,36]. The normal relaxation leads to the creation of an eigenstress σ o s in the surface layer, while the normal relaxation leads to creation of a stress (or initial strain ε ini ) in the core to balance the surface stress. As the surface is bonded to the core, both the surface and core should undergo the same magnitude of deformation. If the coating state after normal relaxation is taken as the reference configuration, the total potential energy of a coating is given by [37]: where h s is the surface layer thickness and h is the overall thicknesses, ε is the strain (bi-axial), L o is the length and width of the coating, u s is the energy density of the surface layer and u c is the energy density of the core. Both u s and u c can be written as [37]: where u o s is the strain-free energy density of the surface layer, u o c is strain-free energy density of the core, Y c is the bulk biaxial Young's modulus; Y s is the 3D surface biaxial Young's modulus and σ o s is the 3D surface eigenstress.
Substituting Equations (2) and (3) into Equation (1) give the potential energy of a coating as: At the equilibrium, the energy minimization requires ∂ U(ε) ∂ε ε=ε ini = 0, which yields the generalized Young-Laplace Equation [38] to describe the mechanical force balance between the surface layer and the core: where F ini c = (h − 2h s )Y c ε ini is the force on core and F ini s = h s (σ o s + Y s ε ini ) is the force on surface. The solution of Equation (5) gives the initial strain which is induced by the surface eigenstress as, Equation (6) can be used to calculate the biaxial surface stress σ s of the coating by using Hooke's law as, It should be noted that the biaxial surface stress of the coating depend on the thickness h of the coating. Surface stress and Young's modulus can be used to calculate the surface porosity P s of coating as [39].
where γ s is the stress sensitivity coefficient [40]. Surface porosity of coating can be related to the surface permeability k s of coating as [41].
where β and m are material constants for coating. Permeability and Porosity can be used to calculate the nano particle size (diameter) D p using modified form of Darcy's law as [42] The relation between average maximum grain size D max with the particle radius D p 2 and its volume fraction (f ) in the nano composite coatings is given as [43]; Average maximum grain size is the average size of grains after deformation has occurred, while the deposit grain size D o is the size of grain at the time of deposition. It is noteworthy that Equation (11) was derived by assuming that particles would not be separated from grain boundaries even when the moving grain boundaries (due to external pressure) drag particles during the relaxation stage.

Corrosion Rate Modelling
The electrochemical reaction of the surface of coating is written as: where M Z+ and M are the metal oxidised and metal reduced forms and Z represents the number of transferred electrons. From the anodic sites on the coating surface the ions are dissolved in the electrolyte and the current flows between the sites and the electrolyte. For coating in a homogenous electrolyte solution, the current simultaneously faces the influence of two external factors: mechanical and electrical. Considering this, it is possible to write the electrochemical potential of the corresponding particles by using the modified form of Gibbs-Duhem expression as [4,26]: where µ ∆P=0 is the chemical potential for zero pressure, V m is the coating molar volume, T is the temperature, Φ is the electrical potential, a is the activity coefficient, µ o represents the standard potential (when a = 1) and F is the Faraday constant [44]. In equilibrium, for Equation (12), the Gibbs free energy ∆G = 0, therefore using Equation (13) ∆G can be written as, where Φ C is the electrical potential of the coating, Φ sol is the electrolyte potential. From Equation (14), the equilibrium potential of the surface reaction corresponding to pressure ∆P is given as [4,45,46]: When any mechanical deformation occurs on the coating, the coating surface will undergo a change in chemical/mechano-electrochemical potential, as per Equation (13), denoted as ∆PV m .
The current density due to the surface electrochemical reactions accompanying surface stress can be written by using well known Butler-Volmer expression as [4,47]: where, where j o is the equilibrium exchange current density corresponding to equilibrium potential ↔ Φ e (η = 0) and α is the charge transfer coefficient. The chemical potential of coating can be changed by applying both tensile and compression stresses, also referred to as hydrostatic stress tensor ∆P. Surface stress is the main factor for the electrochemical reaction rate on the coating surface. By utilizing the biaxial surface stress (σ s ) of the coating in Equation (7), it is possible to obtain the pressure in terms of stress tensor as ∆P = 2σ s /3.
By substituting ∆P into Equations (16)- (18), the current density of the coating can be written as, The corrosion rate of metal electrode depends on the current density of the electrochemical reaction of the coating. According to Faraday's law, there is a linear relationship between the corrosion rate, and the corrosion current density ↔ J .
where corrosion rate is given in mm/year, K 1 is the constant given as 3.27 × 10 −3 (mm·g)/(µA·cm·year), ρ is the density in g/cm 3 and EW is the equivalent weight in grams.

Modelling Results and Discussion
In order to evaluate the intrinsic surface elastic properties of nanocomposite coatings, molecular dynamics (MD) simulations were performed [48]. MD simulations of Ni/Al 2 O 3 , Ni/SiC, Ni/ZrO 2 and Ni/GPL composite coatings face-centered-cubic crystals utilized the LAMMPS code [49] along with the embedded-atom technique [50]. Furthermore, molecular statics (MS) configuration using the conjugate gradient procedure was used during simulations [51].
A stress-free model of nanocomposite coating with atoms positioned within a lattice was created. The simulations of parameters applying boundary conditions were performed only along x and y directions with free surfaces in the z direction. Unreformed volume, length and width L 0 and coating thickness h were evaluated using MD simulations [48]. The unreformed volume of coating multiplied by crystal density was used to evaluate the weight of atoms. The range of coating thickness varied between 100 nm to 5 µm along with 8 × 8 unit cells both along x and y coordinates. The coating relaxation was performed in dual step, i.e., at first normal relaxation along z direction and then parallel relaxation along all three directions, in order to accomplish the minimal energy requirement [48,52]. After parallel relaxation and achieving the equilibrium, the energy reaches a minimum level with length of coating as L ini . Subsequently, the coating's initial strain ε ini was calculated by using ε ini = ln(L ini /L 0 ) . MD simulations were performed for attaining the 3D intrinsic surface elastic parameters of 3D surface eigenstress σ o s , 3D surface biaxial Young's modulus Y s and 3D core biaxial Young's modulus Y c respectively. These parameters are shown in Table 2, subject to assumption that the surface layer thickness h s = 1 nm. Substituting the values in Table 2 into the developed mechano-electrochemical equations, the corrosion rate for all the coatings can be simulated, as shown in Figure 6. The predicted and measured corrosion rates from experiments were compared to validate the model reliability for all the coatings, as shown in Figure 6. The predicted results of corrosion in NaCl environments are in agreement with the experimental results. In the simulation, the model slightly over-predicts the corrosion rates for some data points, which makes some points in this graph deviate from the measured data.   The predicted corrosion rates of Ni/Al2O3, Ni/SiC, Ni/ZrO2 and Ni/GPL with respect to exposure time t are shown in Figure 7. It is clear that the identical corrosion rate curvatures can be observed for all the coatings corresponding to the very short exposure time (i.e., t <25 h). After then, the corrosion rate curvatures reach a maximum and then reduce slowly with passing exposure time, especially in the case of Ni/ZrO2 and Ni/GPL. The reason for this reduction in corrosion rate after 25 h is linked to the increasing grain size with passing time. Lowest reduction in corrosion rate was found in the case of Ni/Al2O3 because the % rise in the grain size for Ni/Al2O3 was largest (also evident from Figure 3). However, the largest reduction in corrosion rate was found in the case of Ni/GPL because the % rise in the grain size for Ni/GPL was smallest. The larger grain size in Ni/Al2O3 resulted in open and disordered nature of grain boundaries, which caused much rapid diffusion of NaCl crystals down boundaries, leading to higher corrosion rate. While Ni/GPL had compact and refined grain structure causing least diffusion of NaCl crystals leading to lowest corrosion rate. The predicted corrosion rates of Ni/Al 2 O 3 , Ni/SiC, Ni/ZrO 2 and Ni/GPL with respect to exposure time t are shown in Figure 7. It is clear that the identical corrosion rate curvatures can be observed for all the coatings corresponding to the very short exposure time (i.e., t <25 h). After then, the corrosion rate curvatures reach a maximum and then reduce slowly with passing exposure time, especially in the case of Ni/ZrO 2 and Ni/GPL. The reason for this reduction in corrosion rate after 25 h is linked to the increasing grain size with passing time. Lowest reduction in corrosion rate was found in the case of Ni/Al 2 O 3 because the % rise in the grain size for Ni/Al 2 O 3 was largest (also evident from Figure 3). However, the largest reduction in corrosion rate was found in the case of Ni/GPL because the % rise in the grain size for Ni/GPL was smallest. The larger grain size in Ni/Al 2 O 3 resulted in open and disordered nature of grain boundaries, which caused much rapid diffusion of NaCl crystals down boundaries, leading to higher corrosion rate. While Ni/GPL had compact and refined grain structure causing least diffusion of NaCl crystals leading to lowest corrosion rate.
The predicted corrosion rates and percentage surface porosity as a function of coating thickness for Ni/Al 2 O 3 , Ni/SiC, Ni/ZrO 2 and Ni/GPL are shown in Figure 8a. It can be seen from Figure 7 that the porosity of extremely thin coatings (less than about 100 nm) is very high as it is directly affected by substrate texture [53,54]. While at greater thicknesses the slope of the porosity-thickness curve is controlled by parameters relevant to the deposit itself. Between these two regimes is a sharp, well-marked transition region in which the porosity of the deposit falls extremely rapidly. The porosity at any given thickness is affected by the deposit grain size. As Ni/Al 2 O 3 has the largest deposit grain size, its porosity-thickness curve accounts for highest porosity and a corresponding highest corrosion rate. While Ni/GPL has the smallest deposit grain size its porosity-thickness curve accounts for lowest porosity and a corresponding lowest corrosion rate. The predicted corrosion rates of Ni/Al2O3, Ni/SiC, Ni/ZrO2 and Ni/GPL with respect to exposure time t are shown in Figure 7. It is clear that the identical corrosion rate curvatures can be observed for all the coatings corresponding to the very short exposure time (i.e., t <25 h). After then, the corrosion rate curvatures reach a maximum and then reduce slowly with passing exposure time, especially in the case of Ni/ZrO2 and Ni/GPL. The reason for this reduction in corrosion rate after 25 h is linked to the increasing grain size with passing time. Lowest reduction in corrosion rate was found in the case of Ni/Al2O3 because the % rise in the grain size for Ni/Al2O3 was largest (also evident from Figure 3). However, the largest reduction in corrosion rate was found in the case of Ni/GPL because the % rise in the grain size for Ni/GPL was smallest. The larger grain size in Ni/Al2O3 resulted in open and disordered nature of grain boundaries, which caused much rapid diffusion of NaCl crystals down boundaries, leading to higher corrosion rate. While Ni/GPL had compact and refined grain structure causing least diffusion of NaCl crystals leading to lowest corrosion rate.  The predicted corrosion rates and percentage surface porosity as a function of coating thickness for Ni/Al2O3, Ni/SiC, Ni/ZrO2 and Ni/GPL are shown in Figure 8a. It can be seen from Figure 7 that the porosity of extremely thin coatings (less than about 100 nm) is very high as it is directly affected by substrate texture [53,54]. While at greater thicknesses the slope of the porosity-thickness curve is controlled by parameters relevant to the deposit itself. Between these two regimes is a sharp, wellmarked transition region in which the porosity of the deposit falls extremely rapidly. The porosity at any given thickness is affected by the deposit grain size. As Ni/Al2O3 has the largest deposit grain size, its porosity-thickness curve accounts for highest porosity and a corresponding highest corrosion rate. While Ni/GPL has the smallest deposit grain size its porosity-thickness curve accounts for lowest porosity and a corresponding lowest corrosion rate. Figure 8b shows the increasing corrosion rate corresponding to decreasing permeability as a function of coating thickness for Ni/Al2O3, Ni/SiC, Ni/ZrO2 and Ni/GPL. This behavior can be directly linked to the porosity level which sets the Oxygen and NaCl electrolyte permeability of coating. As the porosity of coating decreases with the increasing coating thickness, the permeability of Oxygen and electrolyte also decreases with the increasing thickness. Like porosity, the permeability also directly depends on the deposit grain size. Therefore, Ni/Al2O3, which has largest deposit grain size accounts for the highest permeability and corrosion rate, while Ni/GPL, which has the smallest deposit grain size accounts for the lowest permeability and corrosion rate.   Figure 8b shows the increasing corrosion rate corresponding to decreasing permeability as a function of coating thickness for Ni/Al 2 O 3 , Ni/SiC, Ni/ZrO 2 and Ni/GPL. This behavior can be directly linked to the porosity level which sets the Oxygen and NaCl electrolyte permeability of coating. As the porosity of coating decreases with the increasing coating thickness, the permeability of Oxygen and electrolyte also decreases with the increasing thickness. Like porosity, the permeability also directly depends on the deposit grain size. Therefore, Ni/Al 2 O 3 , which has largest deposit grain size accounts for the highest permeability and corrosion rate, while Ni/GPL, which has the smallest deposit grain size accounts for the lowest permeability and corrosion rate.
The predictions in Figure 9 show that, for all the coatings, the corrosion rate increases with increasing average maximum grain size where the average maximum grain size is calculated based on the assumption of spherical grains. The increasing average maximum grain size results in the formation of open, disordered nature of grain boundaries. The grain boundaries make excellent sites for the nucleation of salt precipitates leading to intergranular corrosion. As intergranular corrosion is caused by the precipitation of salt in grain boundaries, the corrosion rate is affected by the volume fraction of precipitated salt per unit grain boundary area. It can be seen that for a given average maximum grain size, the highest corrosion rate was observed in the case of Ni/Al 2 O 3 , while the lowest was observed in the case of Ni/GPL. The reason for the associated corrosion rates is that Ni/Al 2 O 3 accounts for the largest deposit grain size while Ni/GPL accounts for the smallest deposit grain size at the time of deposition. The predictions in Figure 9 show that, for all the coatings, the corrosion rate increases with increasing average maximum grain size where the average maximum grain size is calculated based on the assumption of spherical grains. The increasing average maximum grain size results in the formation of open, disordered nature of grain boundaries. The grain boundaries make excellent sites for the nucleation of salt precipitates leading to intergranular corrosion. As intergranular corrosion is caused by the precipitation of salt in grain boundaries, the corrosion rate is affected by the volume fraction of precipitated salt per unit grain boundary area. It can be seen that for a given average maximum grain size, the highest corrosion rate was observed in the case of Ni/Al2O3, while the lowest was observed in the case of Ni/GPL. The reason for the associated corrosion rates is that Ni/Al2O3 accounts for the largest deposit grain size while Ni/GPL accounts for the smallest deposit grain size at the time of deposition.  Figure 10 shows the relation between predicted corrosion rates as a function of applied biaxial strain. During the deformation of coating under the applied compressive strain, the corrosion rate showed a decreasing trend. While contrarily under the applied tensile strain the corrosion rate showed an increasing trend. The corrosion rate decreases under compressive deformation as pores contract, therefore hindering salt diffusion towards intergranular boundaries [44,55]. On the other hand, the corrosion rate increases under tensile deformation as more pores open up, therefore permitting salts to easily diffuse [56]. Ni/Al2O3 had the highest corrosion rate during compressive and tensile deformations. It is noteworthy that at zero biaxial strain the corrosion rates become minimal, because with no-deformation in the coating system the pores return to normal under stressfree conditions.  Figure 10 shows the relation between predicted corrosion rates as a function of applied biaxial strain. During the deformation of coating under the applied compressive strain, the corrosion rate showed a decreasing trend. While contrarily under the applied tensile strain the corrosion rate showed an increasing trend. The corrosion rate decreases under compressive deformation as pores contract, therefore hindering salt diffusion towards intergranular boundaries [44,55]. On the other hand, the corrosion rate increases under tensile deformation as more pores open up, therefore permitting salts to easily diffuse [56]. Ni/Al 2 O 3 had the highest corrosion rate during compressive and tensile deformations. It is noteworthy that at zero biaxial strain the corrosion rates become minimal, because with no-deformation in the coating system the pores return to normal under stress-free conditions. contract, therefore hindering salt diffusion towards intergranular boundaries [44,55]. On the other hand, the corrosion rate increases under tensile deformation as more pores open up, therefore permitting salts to easily diffuse [56]. Ni/Al2O3 had the highest corrosion rate during compressive and tensile deformations. It is noteworthy that at zero biaxial strain the corrosion rates become minimal, because with no-deformation in the coating system the pores return to normal under stressfree conditions.

Conclusions
A series of combined mechano-electrochemical equations for predicting the corrosion rate of nanocomposite coatings have been developed, which combines the mechanical properties and electrochemical responses of the nanocomposite coatings. The equations take into account the most important parameters in the corrosion of nanocomposite coatings: surface stresses, grain size, porosity, permeability, equilibrium potential and the anodic current density. An extended form of the popular Butler-Volmer expression, combined with the Gibbs-Duhem expression, was used for analyzing the dependence of the corrosion rate on intrinsic mechanical parameters, such as surface stress, surface elastic modulus, surface porosity, permeability and grain size. The predictions showed that the corrosion rate of the Ni/Al 2 O 3 , Ni/SiC, Ni/ZrO 2 and Ni/GPL nanocomposite coatings increased with increasing deposit grain size due to increase in surface stress, porosity and permeability. A model experiment, using accelerated corrosion testing, was performed in conjunction with the analysis. The experiment showed good qualitative agreement with the trends predicted by the theory.