
This work aims to model the influence of biofuels on localized "pitting" corrosion that occurs at the bottom of atmospheric storage tanks. To achieve this purpose, an electro-chemical phase-field model is proposed to include the extra chemical reaction due to the presence of organic acids in an electrolyte solution. The resulting set of nonlinear coupled partial differential equations is numerically integrated by means of finite element methods with a twofold aim: tracking the evolution of the metal/electrolyte interface and predicting the corrosion rates observed when either single or multiple interacting pits are formed in the bottom of a carbon steel tank. The results obtained in the case of single pit, which exhibited a good quantitative agreement with recent experimental data, can be summarized as follows: the presence of organic acids led to higher corrosion rates in comparison with conventional fuels; the corrosion rate is a two-stage process; the dependence of the pit depth as a function of time; and the solid potential, which can be successfully described via a double power law. For multiple interacting pits, the larger corrosivity associated to biofuels was further amplified and the long-time behavior of pit growth gave rise to a "band" behavior, with the major role being played by the number of pits rather than the initial spacings among them. Thus, the proposed model can be employed as a sophisticated tool to predict and quantify the real hazards associated with the release of pollutants in the environment, as well as to optimize the maintenance strategies based on an improved risk-based inspection planning.
Citation: Hossein Moradi, Gabriele Grifò, Maria Francesca Milazzo, Edoardo Proverbio, Giancarlo Consolo. Modeling localized corrosion in biofuel storage tanks[J]. Mathematical Biosciences and Engineering, 2025, 22(3): 677-699. doi: 10.3934/mbe.2025025
[1] | Christoph Sadée, Eugene Kashdan . A model of thermotherapy treatment for bladder cancer. Mathematical Biosciences and Engineering, 2016, 13(6): 1169-1183. doi: 10.3934/mbe.2016037 |
[2] | Heiko Enderling, Alexander R.A. Anderson, Mark A.J. Chaplain, Glenn W.A. Rowe . Visualisation of the numerical solution of partial differential equation systems in three space dimensions and its importance for mathematical models in biology. Mathematical Biosciences and Engineering, 2006, 3(4): 571-582. doi: 10.3934/mbe.2006.3.571 |
[3] | Maria Vittoria Barbarossa, Christina Kuttler, Jonathan Zinsl . Delay equations modeling the effects of phase-specific drugs and immunotherapy on proliferating tumor cells. Mathematical Biosciences and Engineering, 2012, 9(2): 241-257. doi: 10.3934/mbe.2012.9.241 |
[4] | Avner Friedman, Harsh Vardhan Jain . A partial differential equation model of metastasized prostatic cancer. Mathematical Biosciences and Engineering, 2013, 10(3): 591-608. doi: 10.3934/mbe.2013.10.591 |
[5] | Jie Hu, Yangyang Wang, Yuwei Ma, Jiangxiong Wei, Qijun Yu . Identification on acidification damage of external anode system induced by impressed current cathodic protection for reinforced concrete. Mathematical Biosciences and Engineering, 2019, 16(6): 7510-7525. doi: 10.3934/mbe.2019377 |
[6] | Krzysztof Fujarewicz, Krzysztof Łakomiec . Adjoint sensitivity analysis of a tumor growth model and its application to spatiotemporal radiotherapy optimization. Mathematical Biosciences and Engineering, 2016, 13(6): 1131-1142. doi: 10.3934/mbe.2016034 |
[7] | Ting-Hao Hsu, Tyler Meadows, LinWang, Gail S. K. Wolkowicz . Growth on two limiting essential resources in a self-cycling fermentor. Mathematical Biosciences and Engineering, 2019, 16(1): 78-100. doi: 10.3934/mbe.2019004 |
[8] | Yue Xia, Lijuan Chen, Vaibhava Srivastava, Rana D. Parshad . Stability and bifurcation analysis of a two-patch model with the Allee effect and dispersal. Mathematical Biosciences and Engineering, 2023, 20(11): 19781-19807. doi: 10.3934/mbe.2023876 |
[9] | Alexandra Shyntar, Thomas Hillen . Mathematical modeling of microtube-driven regrowth of gliomas after local resection. Mathematical Biosciences and Engineering, 2025, 22(1): 52-72. doi: 10.3934/mbe.2025003 |
[10] | Alexander S. Bratus, Svetlana Yu. Kovalenko, Elena Fimmel . On viable therapy strategy for a mathematical spatial cancer model describing the dynamics of malignant and healthy cells. Mathematical Biosciences and Engineering, 2015, 12(1): 163-183. doi: 10.3934/mbe.2015.12.163 |
This work aims to model the influence of biofuels on localized "pitting" corrosion that occurs at the bottom of atmospheric storage tanks. To achieve this purpose, an electro-chemical phase-field model is proposed to include the extra chemical reaction due to the presence of organic acids in an electrolyte solution. The resulting set of nonlinear coupled partial differential equations is numerically integrated by means of finite element methods with a twofold aim: tracking the evolution of the metal/electrolyte interface and predicting the corrosion rates observed when either single or multiple interacting pits are formed in the bottom of a carbon steel tank. The results obtained in the case of single pit, which exhibited a good quantitative agreement with recent experimental data, can be summarized as follows: the presence of organic acids led to higher corrosion rates in comparison with conventional fuels; the corrosion rate is a two-stage process; the dependence of the pit depth as a function of time; and the solid potential, which can be successfully described via a double power law. For multiple interacting pits, the larger corrosivity associated to biofuels was further amplified and the long-time behavior of pit growth gave rise to a "band" behavior, with the major role being played by the number of pits rather than the initial spacings among them. Thus, the proposed model can be employed as a sophisticated tool to predict and quantify the real hazards associated with the release of pollutants in the environment, as well as to optimize the maintenance strategies based on an improved risk-based inspection planning.
The global transition from traditional fossil fuels to biofuels marks a significant step toward reducing greenhouse gas emissions and promoting sustainable energy practices. Biofuels, which are derived from renewable biological resources, offer an environmentally friendly alternative to conventional petroleum-based fuels [1,2]. However, this shift introduces new challenges, particularly regarding the storage and handling of biofuels. Atmospheric Storage Tanks (ASTs), typically made by carbon steel [3], are essential infrastructure for the biofuel industry. However, they are susceptible to corrosion, which is a major concern for their long-term integrity, as it can lead to the release of hazardous materials causing, in turn, serious damages to the environment and living beings.
Corrosion is the deterioration of a material due to a chemical or electrochemical interaction with its environment [4]. Localized corrosion is a specific type of corrosion that occurs at discrete sites on a metal surface rather than uniformly. This can result in the formation of pits, crevices, or other defects [5]. Pitting corrosion, which is a form of localized corrosion resulting in the formation of small holes [6], is especially prevalent in tank bottoms. One of the challenges that green industries face is the increased susceptibility of ASTs to localized corrosion when storing biofuels [7]. The bottom of an AST is particularly vulnerable to these corrosion processes due to several factors: water accumulation from product contamination, condensation, or rainwater infiltration; the presence of contaminants within the stored product; and stagnant conditions in certain areas of the tank bottom that speed up the process of localized corrosion [8]. In particular, water can dissolve various substances present in the fuel, such as low molecular weight carboxylic acids, carbon dioxide, and chlorides, which are derived from biomass decomposition and other processes [9]. Unlike traditional fossil fuels, biofuels introduce unique challenges due to the presence of corrosion-causing byproducts such as organic acids, aqueous carbonic acids, and chlorides (depending on the feedstock). The term "organic acid" generally refers to a wide variety of organic compounds that contain the carboxyl (COOH) functional group, which plays a crucial role in the corrosion mechanism [9]. Biofuels, such as Fatty Acid Methyl Esters (FAME) and Palm Oil Mill Effluent (POME), are primarily composed of organic fats [10,11,12]. Hydrolysis, oxidation, and microbial activity are the primary mechanisms contributing to the deterioration of FAME and POME. These degradation processes result in the formation of free fatty acids, glycerol, and other corrosive compounds, thus leading to an increase in the total acid number of the biofuel [7,12]. Among the major corrodents identified in biofuel streams, there are low molecular weight organic acids, such as formic acid, acetic acid, and propionic acid, as well as carbon dioxide and chlorides. These stem from biomass decomposition, mineral acids used for hydrolysis, and chlorides found in algae grown in salt water. While some of these chemicals can be removed before further processing, others are either generated at high temperatures or remain uneconomical to eliminate. This issue necessitates the use of corrosion-resistant materials in biofuel processing units [9]. The presence of these corrosive compounds significantly contributes to the development of localized corrosion within ASTs. However, it is crucial to clarify that carbon steel does not exhibit passivity in neutral or low-pH environments [13]. Pitting corrosion in carbon steel ASTs is primarily attributed to material heterogeneities, weld-induced defects, and non-uniform water distribution rather than a classical pitting corrosion mechanism observed in passivated materials [14]. The presence of organic acids in biofuels can contribute to the overall corrosion rates by reducing the pH of accumulated water; however, the primary mechanism of localized attack remains associated with material imperfections and environmental factors.
To ensure the integrity of ASTs and prevent potential leaks, it is essential to regularly monitor the thickness of their bottom sections. This typically involves comprehensive inspections that are conducted every 10 years or more [15,16]. However, discrete thickness measurements cannot exactly determine the maximum corrosion depth of the bottom of storage tanks, where materials may exhibit localized corrosion in the form of pits. To assure the reliability of a given material under a certain environmental condition, quantitative predictions of the local corrosion rate as a function of exposure time and under various environmental influences should be made, which must include the most significant chemical reactions therein occurring.
This work constitutes an attempt to fill this gap, as it aims to develop a predictive mathematical tool for the estimation of corrosion rates in the bottom of ASTs. To achieve this purpose, the model needs to track the time evolution of the interface between the metallic tank and the electrolyte solution, which is generated from the deposition of water in the bottom of the tank. In building such a framework, it becomes crucial to identify those underlying mechanisms which, once the pit initiation stage has been concluded, rule the spatio-temporal propagation of such corrosion pits.
A consolidated framework, named the so-called phase-field (PF), has been used to account for continuously changing complex interfaces and to track the pit growth and the mutual interaction that arises when more localized corrosion sites form and grow in close proximity [17]. This method, commonly known as the "diffuse interface" model, is a mathematical framework that is employed to study interfacial phenomena [18]. The PF description has been employed to describe the evolution of an interface by means of an auxiliary order parameter that undergoes a smooth variation between two phases (metal and electrolyte, in our context) and whose dynamics is ruled by an Allen-Cahn balance equation. This approach avoids the complexities associated with the need of explicitly tracking moving boundaries.
Additionally, modeling the corrosion in the bottom of fuel tanks needs to be complemented by a set of Nernst-Planck equations to describe the transport of ionic species balance of metal ions [17,19,20] and by a Poisson equation which rules the spatial distribution of the electrostatic potential.
In this paper, a relatively simple PF model is employed to qualitatively capture the most relevant electro-chemical reactions that occur at the metal/electrolyte interface in the bottom of a biofuel tank. Indeed, to keep the complexity of the proposed model at a reasonable level, some additional effects, such as those associated with stress corrosion cracking, mechanically- and microbiologically-assisted corrosion, are disregarded here.
The novelties here are introduced are manifold, as they encompass the possibility to numerically investigate how the electrolyte solution is affected by the presence of a biofuel, which additional chemical reactions can be responsible for the larger corrosion rates typically observed in experiments on biofuels in comparison with traditional fuels, how the pit growth rate can be modeled, which dynamics among multiple interacting pits takes place, which characteristic stages the pit undergoes during its growth and which model parameters govern such dynamics.
The paper is organized as follows: in Section 2, the generalized electro-chemical PF model of pitting corrosion, which includes the additional chemical reaction due to the presence of organic acid in the biofuel, is presented; in Section 3, results of numerical investigations are addressed to characterize the time evolution of the pit depth and the corrosion rate observed in configurations with single and multiple interacting pits; and in Section 4, concluding remarks and some future research directions are provided.
To investigate the evolution of the metal/electrolyte interface and predict the corrosion rate caused by single and multiple interacting localized corrosion areas on carbon steel, a PF-based numerical framework is presented, which will be named as pits for simplicity.
Unlike sharp interface models, which require tracking a free boundary, the PF method introduces an auxiliary field variable, φ, defined across the entire domain Ω. This variable undergoes smoothly transitions between the two physical states at the interface, mimicking the two different phases (solid and liquid). In the model presented here, φ = 1 corresponds to the solid phase (metal) while φ = 0 corresponds to the liquid phase (NaCl-based electrolyte solution). A smooth variation of φ occurs within the interfacial region, where the pitting corrosion front evolves. Consequently, the free boundary problem is formulated as a partial differential equation of the Allen-Cahn type, which describes the spatiotemporal evolution of the PF variable. This method allows us to effectively manage complex geometries and to handle topological changes [18,19].
To describe the localized corrosion process in detail, let us consider a 2D domain composed of a metal phase, such as carbon steel, which is exposed to a NaCl solution simulating water solution that can be found on the AST bottom [21]. A schematic of the electrolyte/metal interface, together with some pictures of real biofuel tank bottoms that reveals the presence of circular pits, is shown in Figure 1. Once the dissolution process starts, it leads to the release of metal cations into the electrolyte. The general form of the anodic reaction is represented by the following:
Fe→Fe2++2e−, | (1) |
where Fe represents the corroding metal, and e- refers to electrons released during the process. To reduce the model complexity, it is assumed that the corrosion process solely takes place within the pit area. To this aim, outside the pit area, an inactive layer is introduced to separate the solid phase from the liquid one.
To model the evolution of corrosion waves at the interface between the electrolyte and metal, it is necessary to specifically consider the presence of the following anionic and cationic species: Fe2+, FeOH+, Na+, Cl-, H+, and OH-. Moreover, additional chemicals need to be included because of the extra reactions associated with the specific nature of biofuels that are mainly composed of methyl esters produced through the esterification of oils. When esters interact with water, they break down into organic acids and alcohol [22], thus leading to more organic acid formation. In turn, these acids dissociate in the aqueous environment, thus affecting the acidity of the mixture and lowering the pH.
The above additional chemical reactions will be specified in Section 2.2. In detail, the following subsections tackle different topics. Section 2.1 focuses on the theoretical framework, based upon the PF model, which is used to describe the localized corrosion process. Section 2.2 describes the transport of dilute ionic species, thus providing a detailed analysis of the key chemical reactions involved in the investigated phenomenon of biofuel-induced corrosion. Section 2.3 addresses the spatial distribution of electrostatic potential.
To describe the localized corrosion that occurs in the bottom of biofuel tanks, it is necessary to generalize the framework used with success in the literature [17,19,20]. The time evolution of the PF variable φ is described by an Allen-Cahn equation [23], which reads as follows:
∂φ∂t=−L(∂Ψ∂φ−τ∇2φ), | (2) |
where L denotes the interface kinetics coefficient, Ψ represents the total electrochemical energy, and τ is the gradient energy coefficient. Moreover, the coefficient L can be correlated with the corrosion current density i [19,20], or the overpotential η [17] in the case of activation-controlled corrosion, through the following:
ξ=Li=L0i0, | (3) |
where i0 represents the exchange current density, and L0 represents the interface kinetics coefficient when the overpotential η vanishes. This latter parameter is defined as η = φ0 -φl -Eeq, where φ0 represents the applied (solid) potential, φl represents, the electrostatic potential, and Eeq represents the equilibrium potential (assumed to be zero, for simplicity) [17].
To estimate the corrosion current density i, the Butler-Volmer equation is used under the assumption that the only significant reaction at the dissolution interface is the anodic one:
i=i0exp(anMFηRgT), | (4) |
in such a way that the corrosion current density i takes the form of a Tafel-type [21,24]. In (4), the parameter a represents the anodic charge transfer coefficient, nM represents the charge number of the metal, F is the Faraday constant, Rg is a gas constant, and T is the temperature. According to the above assumptions, L becomes the following:
L=L0exp(anMF(ϕ0−ϕl)RgT). | (5) |
The total electrochemical energy Ψ is an integral of the corresponding energy density ψ over a 2D domain Ω, enclosing the electrolyte solution and the metal layers. The above energy gathers the following contributions:
Ψ=∫Ω(ψ(α)+ψ(β)+ψ(γ))dΩ, | (6) |
where ψ(α), ψ(β), and ψ(γ) represent the chemical, interfacial, and electrostatic free energy densities, respectively. The chemical free energy density can be considered as ψ(α) = ψ(α, d)+ψ(α, s), where ψ(α, d) is the energy for metal dissolution and ψ(α, s) is the energy stored in the dilute solution. In line with the Kim-Kim-Suzuki (KKS) model proposed in [25], each material point is considered as a mixture of solid and liquid phases. Given that the two phases coexist with different concentrations, while maintaining equal chemical potentials, the energy for metal dissolution ψ(α, d) can be expressed as follows:
ψ(α,d)=h(φ)ψ(α,d)S+[1−h(φ)]ψ(α,d)L, | (7) |
where the degradation function is defined as h(φ) = -2φ3 + 3φ2. This function is such that h(φ = 0) = 0 and h(φ = 1) = 1, which identify the fully corroded and the undamaged metal regions, respectively.
The other functions that appear, namely ψ(α,d)S=A(cS−cSe)2 and ψ(α,d)L=A(cL−cLe)2, represent the chemical free energy density terms associated with the concentrations of the solid phase, cS, and the liquid phase, cL, respectively. The quantities involved are as follows: cSe=csolid/csolid=1, which represents the normalized equilibrium concentration for the solid phase; cLe=csat/csolid, which represents the normalized equilibrium concentration for the liquid phase; A, which represents the free energy density parameter (assumed to be identical for both solid and liquid phases); csolid, which represents the concentration of atoms in the metal; and csat, which represents the average saturation concentration. Moreover, additional constraints are given in line with the KKS model [25]. In detail, ∂ψ(α,d)S∂cS=∂ψ(α,d)L∂cL and, by introducing the normalized concentration of metal as cM=c0M/c0Mcsolidcsolid, being c0M the ion metal concentration, and requiring that cM=1 inside the metal and cM=0 in the electrolyte solution (far from the evolving interface), it is possible to express the metal concentration as follows:
cM=h(φ)cS+[1−h(φ)]cL. | (8) |
Finally, by combining all the above assumptions, the energy for metal dissolution can be described by the following equation:
ψ(α,d)=A[cM−h(φ)(cSe−cLe)−cLe]2. | (9) |
On the other hand, the energy stored in the dilute solution is defined as follows:
ψ(α,s)=n∑i=1ciRgT(lnci−1)+ciμ0i, | (10) |
where ci represents the concentration of the i-th ionic species, and μ0i represents the reference chemical potential.
The interfacial energy density, ψ(β), is given by the sum of double-well potential energy and the energy corresponding to the PF gradient:
ψ(β)=g(φ)ω+τ2|∇φ|2, | (11) |
where τ is the gradient energy coefficient, and ω is the height of the double-well potential g(φ) defined as g(φ) = φ2(1 − φ)2. In the present PF corrosion model, these quantities are correlated to the interface energy per area σ and its thickness l via the following [17]:
σ=√τω18,l=√8τω. | (12) |
The electrostatic energy density, ψ(γ) , is expressed as a function of the charge density as follows:
ψ(γ)=Fϕl(nMcMcsolid+n∑i=1nici), | (13) |
where ni is the charge number of the i-th species.
The transport of each aqueous chemical species in the dilute electrolytic solution, with a concentration ci, is ruled by the mass conservation law and consists of three parts [21,24,26]: the migration of ions under concentration gradients, the migration under electrostatic potential gradients, and the rate of depletion or production by chemical reaction. It reads as follows:
∂ci∂t+∇⋅Ni=Πi | (14) |
where Ni and Πi are the total flux in an electrolyte and the chemical reaction term of the i-th species, respectively. The total flux includes phenomena of diffusion (according to Fick's law) and migration (dictated by the electrolyte potential), whereas convection is neglected here, as usual in the literature [17,19,20].
In order to define the Nernst-Planck equations that govern the transport of ionic species, the saline solution that mimics the water collected on the bottom of AST is modeled here as a NaCl-based solution that contains the following six ionic species: c1 ≡ cM = [Fe2+], c2 = [FeOH+], c3 = [Na+], c4 = [Cl-], c5 = [H+], and c6 = [OH-]. This approach holds true when a traditional fuel is stored into the AST. In this case, the most relevant chemical reactions therein occurring are as follows [17,21,26]:
Fe2+(aq)+H2O(aq)k1f⇌k1bFeOH+(aq)+H+(aq) | (15) |
H2O(aq)k2f⇌k2bH+(aq)+OH−(aq) | (16) |
where kb and kf represent the rate constants for the forward and backward reactions, respectively.
Alternatively, when AST is used to store biofuels, it is necessary to account for additional factors related to the occurrence of extra chemical reactions. In detail, biofuels primarily composed of methyl esters R - COOR' produced through the esterification of oils, such as palm oil with methanol, can undergo hydrolysis in the presence of water (Ester Hydrolysis) [27]. Here, "R" denotes a straight or branched chain. Ester hydrolysis is a critical chemical reaction to consider when evaluating the stability and corrosive effects of biodiesel, especially in the context of its interaction with metallic materials [12]. This reaction breaks down the ester bonds, thus leading to the formation of organic acids R - COOH and alcohols R' - OH:
R - COOR'(aq)+H2O(aq)⇌R - COOH(aq)+R' - OH(aq), | (17) |
where R' represents an organyl group. This equation represents the reaction of an ester with water to produce an aliphatic (fatty) acid and an alcohol. In turn, these organic acids dissociate in the electrolyte, thus generating additional ions:
R - COOH(aq)k3f⇌k3bH + (aq)+R - COO−(aq). | (18) |
The dissociation of organic acids contributes to the overall acidity of the electrolyte by introducing additional H+ ions. This increased acidity can accelerate localized corrosion by promoting the dissolution of the metal at the anodic sites [22].
The above chemical reaction constitutes the most significant generalization introduced to describe the different corrosive behavior exhibited by biofuels in AST. This leads to include the presence of two additional species in our framework: c7 = [R-COOH] and c8 = [R-COO-].
By including all the above considerations, and accounting for the concentrations of the mentioned eight ions, the overall Nernst-Planck equations (14) specialize into the following:
∂cM∂t−∇⋅DM∇[cM−h(φ)(cSe−cLe)−cLe]−∇⋅{[1−h(φ)]DMcMRgTnMF∇ϕl}=ΠMcsolid | (19) |
∂ci∂t−∇⋅{[1−h(φ)]Di∇ci}−∇⋅{[1−h(φ)]DiciRgTniF∇ϕl}=Πi(i=2,..,8), | (20) |
where DM and Di are the diffusion coefficients of the metal cation and of the i-th species, respectively. According to (15), (16), and (18), the reaction terms for each involved species, take the following form:
ΠM=ΠFe2+=−k1fc1+k1bc2c5Π2=ΠFeOH+=k1fc1−k1bc2c5Π3=ΠNa+=0Π4=ΠCl−=0Π5=ΠH+=k1fc1−k1bc2c5+k2f−k2bc5c6+k3fc7−k3bc5c8Π6=ΠOH−=k2f−k2bc5c6Π7=ΠR - COOH=−k3fc7+k3bc5c8Π8=ΠR - COO - =k3fc7−k3bc5c8. | (21) |
Notice that sodium Na+ and chloride Cl- ions are assumed not to directly participate in chemical reactions (Π3 = Π4 = 0). However, they play a crucial role in facilitating the transfer of current within the pit by contributing to the overall ionic conductivity of the electrolyte [24].
The distribution of electrostatic potential φl can be estimated through a Poisson-type equation. This equation relates the local variations in the electrostatic potential to the distribution of electric charge within a specific region. It reads as follows [17]:
∇⋅(θ∇ϕl)=nMFcsolid∂φ∂t, | (22) |
where θ is the electric conductivity of the electrolyte.
System (2), (19), (20), and (22) rules that the evolution of localized corrosion in the bottom of biofuel tanks is in the form of a Nernst-Planck-Poisson system [28], which is rewritten for the sake of clarity as follows:
∂φ∂t=−L(∂Ψ∂φ−τ∇2φ), | (23.1) |
∂cM∂t−∇⋅DM∇[cM−h(φ)(cSe−cLe)−cLe]−∇⋅{[1−h(φ)]DMcMRgTnMF∇ϕl}=ΠMcsolid, | (23.2) |
∂ci∂t−∇⋅{[1−h(φ)]Di∇ci}−∇⋅{[1−h(φ)]DiciRgTniF∇ϕl}=Πi(i=2,..,8), | (23.3) |
∇⋅(θ∇ϕl)=nMFcsolid∂φ∂t, | (23.4) |
where cM = [Fe2+], c2 = [FeOH+], c3 = [Na+], c4 = [Cl-], c5 = [H+], c6 = [OH-] and, due to the presence of the biofuel, c7 = [R-COOH] and c8 = [R-COO-]. The reactions terms ΠM and Πi which appear in (23.2) and (23.3) are given in (21).
The complexity of the above system prevents the possibility of addressing further analytical investigations. Therefore, to capture some aspects of qualitative behavior, it is necessary to rely on numerical investigations.
To highlight the main role played by biofuels, a distinction is made between two different frameworks: classical fuel and biofuel. In the former case, organic acids are absent and thus the concentrations of species c7 = [R-COOH] and c8 = [R-COO-] are set to zero. In this case, eight partial differential equations (six balance laws for the species c1, ..., c6, together with those for the PF variable φ and the electrostatic potential φl) describe the key physical phenomena involved. In the latter one, organic acids are present in the electrolyte solution; therefore, it is necessary to account for the chemical reaction (18) and to include two further balance laws for the species R-COOH and R-COO-, coupled with the previously-mentioned ones.
The set of PDEs is numerically integrated by means of finite-elements framework of COMSOL Multiphysics® [29]. To keep the computational cost at a reasonable level, a narrow region of interest was selected at the interface between a layer mimicking the NaCl electrolyte solution having a thickness Hl = 0.6 μm and another layer representing the top surface of the carbon steel fuel tank (the solid phase), having a thickness Hs = 0.4 μm (Figure 1).
Moreover, in order to allow the corrosion process to solely take place within the pit area, an "inactive" layer, having a thickness Hin = 1 µm, is introduced as a meshless layer separating the solid (metal) phase from the liquid (electrolyte) one (see Figure 1(c)). By means of this mathematical artifice, it is intended to describe a layer which does not affect any chemical reaction, does not undergo degradation and, thus, does not participate in the localized corrosion process.
It is also assumed that the pit initiation stage has already occurred [30]; therefore, the processes here solely numerically investigated the pit growth. During this stage, as known, the development of some of the initiated pits occurs, the rate of growth starts to increase and, possibly, the interaction among different pits takes place. In particular, pits with a radius rpit = 10 µm constitute the initial condition for all the numerical simulations. This choice allows us to improve the convergence of the numerical algorithm and is in line with the literature [17,19].
To accurately assess the corrosion rates, the simulations are carried out over a time window of 30,000 sec and considers a computational domain having a width W = 0.6 μm.
Table 1 outlines the initial conditions for various field variables within the system. These conditions represent the state of the system at the beginning of the simulation. The variables include the PF variable (φ), the metal ion concentration cM, various ionic species ci, and the electric potential φl. The liquid and solid phases are differentiated, with specific values assigned to each phase for each variable. Table 2 details the boundary conditions for the PF variable, the metal ion concentration, and other ionic species, imposed at the top surface of the electrolyte layer (x = 1000 μm) and at the bottom of the metal layer (x = 0 μm). Moreover, zero flux boundary conditions are applied to the lateral boundaries. The additionally parameter values used in all simulations are reported in Table 3.
Field variable | Liquid phase | Solid phase |
φ | 0 | 1 |
cM/csolid = [Fe2+]/csolid | 0 | 1 |
c2 = [FeOH+] | 0 mol/m3 | 0 mol/m3 |
c3 = [Na+] | 493.65 mol/m3 | 0 mol/m3 |
c4 = [Cl-] | 1 mol/m³ | 0 mol/m3 |
c5 = [H+] | 1.5849E-2 mol/m³ | 0 mol/m3 |
c6 = [OH-] | 6.3809E-7 mol/m³ | 0 mol/m3 |
c7 = [R-COOH] | 331.01 mol/m³ | 0 mol/m3 |
c8 = [R-COO-] | 416.71 mol/m³ | 0 mol/m3 |
φl | 0 V | 0 V |
Field variable | Top surface of electrolyte layer (x = 1000 μm) | Bottom of the metal layer (x = 0 μm) |
φ | 0 | 1 |
cM/csolid = [Fe2+]/csolid | 0 | 1 |
c2 = [FeOH+] | 0 mol/m3 | 0 mol/m3 |
c3 = [Na+] | 493.65 mol/m3 | 0 mol/m3 |
c4 = [Cl-] | 1 mol/m³ | 0 mol/m3 |
c5 = [H+] | 1.5849E-2 mol/m³ | 0 mol/m3 |
c6 = [OH-] | 6.3809E-7 mol/m³ | 0 mol/m3 |
c7 = [R-COOH] | 331.01 mol/m³ | 0 mol/m3 |
c8 = [R-COO-] | 416.71 mol/m³ | 0 mol/m3 |
ϕl | 0V | ∂ϕ/∂ϕ∂x=0∂x=0V/m |
Parameter | Description | Value | Unit of measure |
A | Free energy density curvature | 5.35E7 | N/m² |
α | Anodic charge transfer coefficient | 0.5 | - |
csolid | Average concentration of metal | 1.43E5 | mol/m³ |
csat | Average saturation concentration | 5100 mol/m³ | mol/m³ |
DM, Di | Diffusion coefficient | 8.5E-10 | m²/s |
K1 | Chemical equilibrium constant | 1.625E-4 | mol/m³ |
K2 | Chemical equilibrium constant | 1.0113E-8 | mol²/m⁶ |
K3 | Chemical equilibrium constant | 0.028184 | mol/m³ |
k1f | Forward reaction rate | 0.08125 | 1/s |
k2f | Forward reaction rate | 5.0565E-5 | mol/(m³·s) |
k3f | Forward reaction rate | 1409.2 | 1/s |
k1b | Backward reaction rate | 500 | m³/(s·mol) |
k2b | Backward reaction rate | 5000 | m³/(s·mol) |
k3b | Backward reaction rate | 50,000 | m³/(s·mol) |
nM | Average charge number: [Fe2+] | 2.1 | - |
nc2 | Average charge number: [FeOH+] | 1.1 | - |
nc3 | Average charge number: [Na+] | 1 | - |
nc4 | Average charge number: [Cl-] | -1 | - |
nc5 | Average charge number: [H+] | 1 | - |
nc6 | Average charge number: [OH-] | -1 | - |
nc7 | Average charge number: [R-COOH] | 0 | - |
nc8 | Average charge number: [R-COO-] | -1 | - |
T | Temperature | 300 | K |
ω | Height of the double-well potential | 2.4E7 | N/m² |
Rg | Gas constant | 8.314 | J/(mol·K) |
τ | Gradient energy coefficient | 7.5E-5 | N |
θ | Electric conductivity | 1E7 | S/m |
The numerical values reported in these tables are extracted from the literature [17,21,24]. The only exception is constituted by the values of the initial concentrations of R-COOH and R-COO- reported in Table 1, which are not generally available in the literature. In the following, the procedure to deduce their values is outlined.
Let us introduce the concept of the partition coefficient κ, which quantifies the distribution of a solute between two immiscible phases [31]. In our case, κ will quantify the distribution of the organic acid between the aqueous electrolyte (where it dissociates) and the palmitic oil phase:
κ=[organic acids]oil[organic acids]water. | (24) |
Moreover, as outlined in [17], the equilibrium constants K1 associated with reaction (15) quantifies as follows:
K1=k1fk1b=c2c5cMcsolid. | (25) |
Similarly, K2 characterizes the equilibrium of reaction (16) and is determined by the following:
K2=k2fk2b=c5c6. | (26) |
In the specific case of acid dissociation, where a weak acid forms its conjugate base, the equilibrium constant is given by Ka = [R-COO-][H+]/[R-COOH], and is referred as the acid dissociation constant. On the other hand, pKa is related to the acid dissociation constant through pKa = -log(Ka) [32,33,34].
According to the reaction (18), the dissociation equilibrium constant K3 ≡ Ka of the organic acid in the electrolyte solution is given by the following:
K3=k3fk3b=10−pKa. | (27) |
Therefore, taking the partition coefficient for acetic acid between oil and water given in [35] κ = 0.075 into account, and extracting the mass concentration of acetic acid in oil so through the knowledge of typical molar concentration and molar mass [36],
so = 0.246molL⋅65.05gmol=16.0023gL, | (28) |
it is possible to deduce the following:
sw=soκ=16.0023gL0.075=213.364gL. | (29) |
Then, let us consider a mixture containing 0.15% water and 99.85% oil. On the total amount, the weight percentage of acetic acid in water Ww and the weight percentage of acetic acid in oil Wo are given by Ww = 1.96% and Wo = 98.04%, respectively.
Therefore, the solubility of acetic acid in the mixture saa is determined by considering the solubility in both water and oil:
saa=Ww⋅sw+Wo⋅so=19.88gL. | (30) |
Once the solubility has been determined, considering the molar mass of acetic acid PM(acetic_acid) = 60.052 g/mol, the concentration of the undissociated form [R-COOH] is calculated by the following:
[R-COOH]=saaPM(acetic\_acid)=331.01molm3. | (31) |
Finally, considering that the pKa value for acetic acid is 4.7, as reported in [37], and that the initial pH of the solution is 4.8 (see Table 1), the concentration of the dissociated form [R-COO-] is given by the following:
[R-COO−]=10−pKa[R-COOH]10−pH=416.71molm3. | (32) |
Numerical investigations have been performed by following this criterion. First, the corrosion process is modeled under the simplistic assumption that only a single pit forms and evolves. In such an analysis, an attempt is made to extract a suitable predictive law and highlight the most relevant differences between classical fuels and biofuels. Then, the analysis solely focuses on the biofuel, and the most relevant features that occur when two or three pits interact and evolve on the metal surface in close proximity are described.
It is expected that the most relevant physical quantities that contribute to the evolution of a single pit are as follows: time t and solid potential φ0. To investigate this dependence in detail, several numerical simulations have been carried out by considering a computational time window t = [0, 3e4] sec, and by varying the solid potential in the range φ0 = [0.3, 0.42] V. To accelerate observations on pit growths within a computationally feasible timeframe, the chosen solid potential values are intentionally larger than those that might be encountered in real-world applications.
Then, a nonlinear regression analysis on these numerical data is performed, for both classical and biofuels, in order to find the best fit for the function depth(t, φ0), which describes the time evolution of the bottom of the pit (i.e., the worst situation). This analysis has led to the following expression:
depth(t,ϕ0)=B0(tt∗)B1(ϕ0ϕ0∗)B2 | (33) |
where B0, B1, and B2 are fitting parameters, whereas t* and φ0* are normalization coefficients. In detail, B0 is a characteristic depth which also determines the pit growth rate, and B1 and B2 describe the power laws for time and solid potential, respectively. The normalization coefficients are arbitrarily set in such a way that t* = 31,536,000 sec, which corresponds to the period of 1 year, whereas φ0* = 0.42 Volt is the largest value considered for solid potential.
For classical fuel, the best-fit values of these parameters are as follows:
B0=22803.19μm,B1=0.64178,B2=9.0612, |
whereas, for biofuel,
B0=25505.81μm,B1=0.65476,B2=9.23315. |
Among the possible functional dependencies that can be considered to fit a real function of two-variables, it should be mentioned that this law appears to be attracting since it encloses just three fitting parameters and, above all, the power laws have been already used with success in the literature, especially to describe the time evolution of pit growth [30,38].
The numerical fits reported in Eq (33) are depicted in Figure 2 through surfaces for classical fuel (panel (a)) and biofuel (panel (b)). In the same figures, the results of numerical simulations are represented by solid lines. A regression analysis for both classical fuel and biofuel yielded an adjusted R-square value of 0.997, thus indicating a very strong fit of the model to the data, as it also appears from the close agreement reported in the figures. From the close agreement between the numerical and analytical data, it can be argued that the model in Eq (33) appears to be quite suitable to predict both short-time and long-time behavior (i.e., from the conclusion of the pit initiation stage to the achievement of quasi-steady-state).
Additionally, the results depicted in this figure clearly indicate that the pit depth observed in biofuels is always larger than that reported for classical fuels, which is in qualitative agreement with literature [22].
Moreover, from the inspection of the above fitting data, it is worth noticing that the parameter providing the most relevant contribution to the higher corrosion rate exhibited by biofuel is the prefactor B0, which exhibits an increase of about 12% with respect to that of a classical fuel. Interestingly, all the other parameters exhibit negligible differences (smaller than 2%), so that the power laws describing time and potential dependencies may be ruled by the same function. Noteworthy, this enables the possibility to characterize the typology of a fuel solely by means of a single parameter, that is, the role played by additional chemical reactions due to the presence of organic acids in the electrolyte solution may be enclosed in the sole coefficient B0.
In our view, the results discussed here represent one of the most relevant achievements of this work as the possibility to manage a best fit function allows us to address some quantitative predictions on the time required to achieve a given pit depth once the solid potential is known. This constitutes a more than helpful tool for green industries to optimize the protocol adopted to test the integrity of bottom tanks.
For instance, it is possible to address a quantitative comparison with experimental data reported by several groups [12,22,27,39,40,41], where a quite wide range of corrosion rates due to biofuel in contact with carbon steel was reported, from 1 μm/year to about 2 mm/year. For instance, by considering the value of pit depth equal to 17 μm reported after 8 hours of laboratory observations (that roughly corresponds to our computational time window of 30.000 sec) (see Figure 28 in [40]), it is possible to get an indirect measure of the potential φ0 = 0.3125 V by means of Eq (33). Then, by considering this value of applied potential, the same Eq (33) allows us to estimate the final pit depth after 1 year, which is 1.66 mm, which is again in line with the above literature.
In Figure 3(a), the time evolution of a single pit for the two fuels is compared for different values of solid potentials φ0. From this figure, it can be seen that the larger corrosivity of biofuel increases with the rise in solid potential. Apart from the time evolution of the pit depth, it is interesting to evaluate its time derivative, which is representative of the corrosion rate and provides significant information on the different stages occurring during the corrosion process. Figure 3(b) shows the behavior obtained for a specific applied potential (φ0 = 0.4 V), where the behavior associated to the classical fuel is given as the solid line while the one exhibited by biofuel is given via the dashed line.
Two key pieces of information might be extracted from this figure. First, the higher corrosivity of the biofuel manifests itself from the very beginning. Second, the time evolution of the rate of corrosion, in the case of a single pit, is a two-steps process characterized by an initial sudden increase up to a given maximum, followed by a slow decrease over time. This pattern, which appears to be independent of the fuel type, suggests that it might be a general and common feature of the time evolution of localized corrosion. It is also interesting to notice that this two-steps process was also observed in the literature of corrosion of metal in biodiesel [41] and even in other contexts (e.g., marine corrosion) [30]. The initial rapid growth is typically associated with favorable conditions (i.e., large supply of chemical agents) that support a high rate of corrosion. However, this rate reduces to a quasi-steady state situation that is supposed to continue indefinitely and whose precise underlying mechanisms are still unclear.
From now on, the investigations will solely focus on the behavior of the biofuel. Additionally, for simplicity but without loss of generality, the analysis is limited to a solid potential of φ0 = 0.4V. The behavior associated with other values of the potential is qualitatively unchanged.
In Figure 4, the time evolutions of the pit depth (panel (a)) and corrosion rate (panel (b)) are depicted for frameworks where a single, double, and triple pits are initially formed. In simulations with multiple pits, the initial condition consists of pits which are equal-in-size (radius rpit = 10 μm) and located at a distance d = 20 µm from each other.
In the single pit scenario (dash-dotted line) already discussed in the previous figure, the maximum depth achieved at the final simulation time is approximately 180 µm. Notably, the presence of multiple pits significantly impacts the corrosion process. Indeed, in the double pit configuration (dashed line), the final depth is around 200 µm and reaches about 220 µm in the case of triple pits (solid line) (i.e., an increase of about 20% (according to this parameter set)). These observations strongly suggest that the interaction between adjacent pits exacerbates the corrosion process, thus leading to a more rapid increase in the pit depth. The analysis of the corrosion rate (panel (b)) revealed that the maximum corrosion rate obtained in the multiple pit configuration (double and triple), about 13.5 nm/s, is about 40% larger than that reported in the single pit case, 8.5 nm/s. Two key differences may be pointed out between the two multiple-pits configurations: in the triple one, the peak is slightly delayed and, after it is achieved, the decreasing trend of corrosion rate appears to be slower than that observed in the double pit.
A superficial analysis of Figure 4 might erroneously lead to the conclusion that the corrosion process in the presence of multiple pits is a two-stage process similar to the one observed in Figure 3(b). However, deeper investigations reveal that it is not. To prove it, let us consider a double-pit configuration with an initial spacing d = 20 μm and track the time evolution of the corrosion rate as well as the spatial profiles of the PF variable φ at each time instant, with a particular emphasis at those at which the corrosion rate exhibits a change in slope. The results of this analysis are reported in Figure 5.
In detail, the simulation starts with the initial (t = 0 s) configuration with two equal pits, see point P1 in panel (b). As the simulation time progresses, at P2 (t = 750 s), the corrosion rate reaches the value of about 10 nm/s and the curve exhibits a kink. At this point, the rate of corrosion has indeed risen. The reason for that is that a change of morphology is taking place: the boundaries of the two pit areas are now in contact, as shown in Figure 5(c). During the subsequent time instants, from P2 to P3, the initially distinct pits progressively merge together and form a unique pit with a larger volume, as can be appreciated in Figure 5(d), in line with both experiments [42] and other numerical simulations [19,43]. Here, the maximum corrosion rate of about 13.5 nm/s is obtained. Starting from P3, the evolution is the one typical of a single pit.
To gain further insight into the process of corrosion due to multiple pits, let us now vary the initial spacing between them. In the case of double pits, d is considered in the range (20,160) μm, while for triple pits, d is applied in the range (10, 70) μm. These values are chosen to be compatible with the considered width of each layer and the simulation time window. Figure 6 reports the overall results (panel (a) for double, panel (b) for triple). In the same figure, small squares denote the time instant at which the boundary of the adjacent pit areas become in touch with each other. The results reveal that the closest pits get in touch earlier, as expected, and thus the corresponding pit depth exhibits a faster increase. However, these two combined observations don't imply that the pit depth achieves the largest value at the final simulation time. Indeed, since the corrosion rate starts monotonically decreasing once the maximum is reached, it happens that, after a substantial long time (in our case, the upper bound of the simulated time window), the pit depth observed for d = 20 μm appears to be lower than those observed for larger spacings. This behavior holds for double and triple pits.
The above considerations may be further supported and enriched by addressing an overall comparison among all the time evolutions of the pit depth for single, double, and triple pit configurations, for different values of initial spacings between the pits, which is in line with Figure 6. The results of this analysis are reported in Figure 7.
Here, the dashed line depicts the pit depth evolution for a single pit and serves as a baseline for the comparative analysis, as it represents the lower limit for the corrosion process. Dashed and solid lines denote the behavior for double and triple pits, respectively. The evidence that the increase of interacting sources enhances the pit depth is not the only conclusion suggested by this figure. In fact, these results reveal that the long-time behavior associated to a given number of interacting pits, independently of the considered initial spacings between pits, forms a "band", which appears disjointed with respect to the band formed by a different number of interacting pits. In other words, the long-time dynamics results to be mostly determined by the number of interacting pits rather than the initial spacings among pits. In our view, this is another intriguing feature to be considered in the design of the maintenance system of fuel tanks.
AST bottoms are critical components susceptible to corrosion due to the accumulation of sediments and exposure to biofuels. This study employed a PF model to investigate the influence of biofuels on pitting corrosion at the bottom of ASTs made by carbon steel. To emphasize the role played by biofuel, the model presented here also accounted for the additional chemical reaction which described the dissociation of organic acids which are present in the electrolyte solution. Thus, the overall system included an equation for each of the eight ionic species, one for the PF variable and one for the electrostatic potential. The obtained nonlinear and coupled system was numerically integrated to describe the spatio-temporal evolution observed during the stage of pit growth. In particular, the dynamics of single and multiple interacting pits were inspected, and results agree with some experimental observations.
The most relevant achievements contained in this manuscript can be summarized as follows:
● A best-fit law describing the dependence of the pit depth as a function of time and solid potential was deduced. This function, apart from being characterized by just three free parameters (B0, B1, B2), yielded the most significant advantage of highlighting that the presence of organic acids in the electrolyte solution may be enclosed in the sole coefficient B0. Thus, this fitting parameter can be used to identify and characterize the typology of biofuel.
● It is confirmed that biofuels exhibit a larger corrosivity in comparison with traditional fuels.
● As far as single, isolated, pits are concerned, the pit growth rate is a two-steps process.
● When multiple pits are formed in close proximity, their interaction yields a threefold contribution:
○ it enhances the pit growth;
○ the pit growth consists of three distinct stages; and
○ the time evolution of pit growth may be classified in terms of separated "bands", within which the curves representative of different initial spacings among pits fall.
These latter observations served to stress the importance of considering pit interactions during inspections and risk assessments for ASTs. At the same time, they suggested that while the proximity of multiple pits becomes crucial for the corrosion of the bottom of the tank, as it may potentially led to enhanced localized breaches and leaks, the number of such interacting pits plays a relevant role rather than their initial spacings.
It is believed that the findings described here may contribute to acquiring a better understanding on the influence of biofuels and pit interactions on corrosion behavior at the tank bottom. Additionally, they may help to develop more robust design strategies for AST bottoms, to optimize the inspection procedures to target areas susceptible to multiple interacting pits, and to select appropriate materials to ensure the long-term integrity and reliability of these critical storage systems.
Future research directions could include the following: modeling additional contributions arising from mechanically- and microbiologically-assisted corrosion and the inclusion of stochastic effects describing the initial formation of several, randomly dislocated, pits.
The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.
Authors gratefully acknowledge fruitful discussions with Prof. C. Cui and E. Martinez-Paneda. This research was funded by MUR (Italian Ministry of University and Research) PNRR - Missione 4, Componente 2, Investimento 1.1 – Bando Prin 2022, Decreto Direttoriale no. 104 del 02-02-2022, through PRIN2022 Project PRIN_202248TY47 ''Modelling complex biOlogical systeMs for biofuEl productioN and sTorAge: mathematics meets green industry (MOMENTA)" CUP J53D23003580006; by INdAM-GNFM and by INAIL within the BRIC/2021 ID = 3 project DRIVERS.
The authors declare there is no conflict of interest.
[1] |
L. Cherwoo, I. Gupta, G. Flora, R. Verma, M. Kapil, S. K. Arya, et al., Biofuels an alternative to traditional fossil fuels: A comprehensive review, Sustain. Energ. Technol. Assess. , 60 (2023), 103503. https://doi.org/10.1016/j.seta.2023.103503 doi: 10.1016/j.seta.2023.103503
![]() |
[2] |
M. V. Rodionova, R. S. Poudyal, I. Tiwari, R. A. Voloshin, S. K. Zharmukhamedov, H. G. Nam, et al., Biofuel production: Challenges and opportunities, Int. J. Hydrogen Energ. , 42 (2017), 8450–8461. https://doi.org/10.1016/j.ijhydene.2016.11.125 doi: 10.1016/j.ijhydene.2016.11.125
![]() |
[3] | API Standard 650, 13th Ed., Available from: https://www.api.org/products-and-services/standards/important-standards-announcements/standard650. |
[4] | N. Perez, Electrochemical corrosion, in Electrochemistry and Corrosion Science, Springer International Publishing, Cham, (2016), 1–23. https://doi.org/10.1007/978-3-319-24847-9_1 |
[5] | D. A. Shifler, Localized corrosion, in Corrosion in Marine Environments, John Wiley & Sons, Inc., (2022), 63–121. https://doi.org/10.1002/9781119788867.ch4 |
[6] | J. R. Galvele, Pitting corrosion, in Treatise on Materials Science and Technology (ed. J. C. Scully), Elsevier, (1983), 1–57. https://doi.org/10.1016/B978-0-12-633670-2.50006-1 |
[7] |
L. N. Komariah, S. Arita, B. E. Prianda, T. K. Dewi, Technical assessment of biodiesel storage tank: a corrosion case study, J. King Saud. Univ. Eng. Sci., 35 (2023), 232–237. https://doi.org/10.1016/j.jksues.2021.03.016 doi: 10.1016/j.jksues.2021.03.016
![]() |
[8] |
M. F. Milazzo, G. Ancione, P. Bragatto, E. Proverbio, A probabilistic approach for the estimation of the residual useful lifetime of atmospheric storage tanks in oil industry, J. Loss Prevent. Proc. Ind. , 77 (2022), 104781. https://doi.org/10.1016/j.jlp.2022.104781 doi: 10.1016/j.jlp.2022.104781
![]() |
[9] | C. A. Shargay, K. Moore, M. West, NACE-2014-3729, Corrosion 2014, San Antonio, Texas, USA, 2014. |
[10] |
K. A. Zahan, M. Kano, Biodiesel production from palm oil, its by-products, and mill effluent: A review, Energies, 11 (2018), 2132. https://doi.org/10.3390/en11082132 doi: 10.3390/en11082132
![]() |
[11] |
M. Rostami, S. Raeissi, M. Mahmoodi, M. Nowroozi, Liquid–Liquid equilibria in biodiesel production, J. Am. Oil Chem. Soc. , 90 (2013), 147–154. https://doi.org/10.1007/s11746-012-2144-5 doi: 10.1007/s11746-012-2144-5
![]() |
[12] | A. Groysman, Corrosion in Systems for Storage and Transportation of Petroleum Products and Biofuels: Identification, Monitoring and Solutions, Springer Dordrecht, 2014. https://doi.org/10.1007/978-94-007-7884-9 |
[13] |
H. Parangusan, J. Bhadra, N. Al-Thani, A review of passivity breakdown on metal surfaces: influence of chloride- and sulfide-ion concentrations, temperature, and pH, Emerg. Mater. , 4 (2021), 1187–1203. https://doi.org/10.1007/s42247-021-00194-6 doi: 10.1007/s42247-021-00194-6
![]() |
[14] |
Y. Tan, Understanding the effects of electrode inhomogeneity and electrochemical heterogeneity on pitting corrosion initiation on bare electrode surfaces, Corros. Sci. , 53 (2011), 1845–1864. https://doi.org/10.1016/j.corsci.2011.02.002 doi: 10.1016/j.corsci.2011.02.002
![]() |
[15] | EEMUA Publication 159 Digital, Available from: https://www.eemua.org/Products/Publications/Digital/EEMUA-Publication-159.aspx. |
[16] | API 580 - Risk Based Inspection, 2023. Available from: https://www.api.org/products-and-services/individual-certification-programs/certifications/api580. |
[17] |
C. Cui, R. Ma, E. Martínez-Pañeda, Electro-chemo-mechanical phase field modeling of localized corrosion: theory and COMSOL implementation, Eng. Comput. , 39 (2023), 3877–3894. https://doi.org/10.1007/s00366-023-01833-8 doi: 10.1007/s00366-023-01833-8
![]() |
[18] |
W. J. Boettinger, J. A. Warren, C. Beckermann, A. Karma, Phase-field simulation of solidification, Annu. Rev. Mater. Res. , 32 (2002), 163–194. https://doi.org/10.1146/annurev.matsci.32.101901.155803 doi: 10.1146/annurev.matsci.32.101901.155803
![]() |
[19] |
W. Mai, S. Soghrati, R. G. Buchheit, A phase field model for simulating the pitting corrosion, Corros. Sci. , 110 (2016), 157–166. https://doi.org/10.1016/j.corsci.2016.04.001 doi: 10.1016/j.corsci.2016.04.001
![]() |
[20] |
C. Cui, R. Ma, E. Martínez-Pañeda, A phase field formulation for dissolution-driven stress corrosion cracking, J. Mech. Phys. Solids, 147 (2021), 104254. https://doi.org/10.1016/j.jmps.2020.104254 doi: 10.1016/j.jmps.2020.104254
![]() |
[21] |
S. M. Sharland, P. W. Tasker, A mathematical model of crevice and pitting corrosion—Ⅰ. The physical model, Corros. Sci. , 28 (1988), 603–620. https://doi.org/10.1016/0010-938X(88)90027-3 doi: 10.1016/0010-938X(88)90027-3
![]() |
[22] |
A. T. Hoang, M. Tabatabaei, M. Aghbashlo, A review of the effect of biodiesel on the corrosion behavior of metals/alloys in diesel engines, Energy Sources Part A, 42 (2020), 2923–2943. https://doi.org/10.1080/15567036.2019.1623346 doi: 10.1080/15567036.2019.1623346
![]() |
[23] |
J. W. Cahn, J. E. Hilliard, Free energy of a nonuniform system. Ⅰ. interfacial free energy, J. Chem. Phys. , 28 (1958), 258–267. https://doi.org/10.1063/1.1744102 doi: 10.1063/1.1744102
![]() |
[24] |
S. M. Sharland, A mathematical model of crevice and pitting corrosion—Ⅱ. The mathematical solution, Corros. Sci. , 28 (1988), 621–630. https://doi.org/10.1016/0010-938X(88)90028-5 doi: 10.1016/0010-938X(88)90028-5
![]() |
[25] |
S. G. Kim, W. T. Kim, T. Suzuki, Phase-field model for binary alloys, Phys. Rev. E, 60 (1999), 7186–7197. https://doi.org/10.1103/PhysRevE.60.7186 doi: 10.1103/PhysRevE.60.7186
![]() |
[26] |
S. M. Sharland, C. P. Jackson, A. J. Diver, A finite-element model of the propagation of corrosion crevices and pits, Corros. Sci. , 29 (1989), 1149–1166. https://doi.org/10.1016/0010-938X(89)90051-6 doi: 10.1016/0010-938X(89)90051-6
![]() |
[27] |
L. M. Baena, J. A. Calderón, Effects of palm biodiesel and blends of biodiesel with organic acids on metals, Heliyon, 6 (2020), e03735. https://doi.org/10.1016/j.heliyon.2020.e03735 doi: 10.1016/j.heliyon.2020.e03735
![]() |
[28] |
B. E. McNealy, J. L. Hertz, Extended Poisson–Nernst–Planck modeling of membrane blockage via insoluble reaction products, J. Math. Chem. , 52 (2014), 430–440. https://doi.org/10.1007/s10910-013-0270-4 doi: 10.1007/s10910-013-0270-4
![]() |
[29] | COMSOL Multiphysics® v. 6.2 (2023), COMSOL AB, Stockholm, Sweden. Available from: http://www.comsol.com. |
[30] |
J. Bhandari, F. Khan, R. Abbassi, V. Garaniya, R. Ojeda, Modelling of pitting corrosion in marine and offshore steel structures – A technical review, J. Loss Prevent. Proc. Ind. , 37 (2015), 39–62. https://doi.org/10.1016/j.jlp.2015.06.008 doi: 10.1016/j.jlp.2015.06.008
![]() |
[31] |
M. A. Reinsel, J. J. Borkowski, J. T. Sears, Partition coefficients for acetic, propionic, and butyric acids in a crude oil/water system, J. Chem. Eng. Data, 39 (1994), 513–516. https://doi.org/10.1021/je00015a026 doi: 10.1021/je00015a026
![]() |
[32] | J. E. McMurry, R. C. Fay, Chemistry, 6th ed., Pearson Education, 2011. |
[33] |
J. Drzymala, Chemical equilibria in the oleic acid-water-NaCl system, J. Colloid Interface Sci. , 108 (1985), 257–263. https://doi.org/10.1016/0021-9797(85)90259-0 doi: 10.1016/0021-9797(85)90259-0
![]() |
[34] | P. Atkins, J. de Paula, R. Friedman, Physical Chemistry: Quanta, Matter, and Change, OUP Oxford, 2014. https://doi.org/10.1093/hesc/9780199609819.001.0001 |
[35] |
B. Rangarajan, A. Havey, E. A. Grulke, P. D. Culnan, Kinetic parameters of a two-phase model for in situ epoxidation of soybean oil, J. Am. Oil Chem. Soc. , 72 (1995), 1161–1169. https://doi.org/10.1007/BF02540983 doi: 10.1007/BF02540983
![]() |
[36] |
D. Gamaralalage, O. Sawai, T. Nunoura, Degradation behavior of palm oil mill effluent in Fenton oxidation, J. Hazard. Mater., 364 (2019), 791–799. https://doi.org/10.1016/j.jhazmat.2018.07.023 doi: 10.1016/j.jhazmat.2018.07.023
![]() |
[37] | National Center for Biotechnology Information, PubChem Compound Summary for CID 985, Palmitic Acid, 2025. Available from: https://pubchem.ncbi.nlm.nih.gov/compound/Palmitic-Acid. |
[38] |
Y. Katano, K. Miyata, H. Shimizu, T. Isogai, Predictive model for pit growth on underground pipes, Corrosion, 59 (2003), 155–161. https://doi.org/10.5006/1.3277545 doi: 10.5006/1.3277545
![]() |
[39] |
F. Pessu, R. Barker, A. Neville, Understanding pitting corrosion behavior of X65 carbon steel in CO2-saturated environments: The temperature effect, Corrosion, 72 (2015), 78–94. https://doi.org/10.5006/1338 doi: 10.5006/1338
![]() |
[40] |
F. Pessu, R. Barker, A. Neville, CO2 corrosion of carbon steel: The synergy of chloride ion concentration and temperature on metal penetration, Corrosion, 76 (2020), 957–974. https://doi.org/10.5006/3583 doi: 10.5006/3583
![]() |
[41] |
M. A. Fazal, A. S. M. A. Haseeb, H. H. Masjuki, Corrosion mechanism of copper in palm biodiesel, Corros. Sci. , 67 (2013), 50–59. https://doi.org/10.1016/j.corsci.2012.10.006 doi: 10.1016/j.corsci.2012.10.006
![]() |
[42] |
M. Ghahari, D. Krouse, N. Laycock, T. Rayment, C. Padovani, M. Stampanoni, et al., Synchrotron X-ray radiography studies of pitting corrosion of stainless steel: Extraction of pit propagation parameters, Corros. Sci. , 100 (2015), 23–35. https://doi.org/10.1016/j.corsci.2015.06.023 doi: 10.1016/j.corsci.2015.06.023
![]() |
[43] |
T. Q. Ansari, Z. Xiao, S. Hu, Y. Li, J. L. Luo, S. Q. Shi, Phase-field model of pitting corrosion kinetics in metallic materials, NPJ Comput. Mater. , 4 (2018), 38. https://doi.org/10.1038/s41524-018-0089-4 doi: 10.1038/s41524-018-0089-4
![]() |
Field variable | Liquid phase | Solid phase |
φ | 0 | 1 |
cM/csolid = [Fe2+]/csolid | 0 | 1 |
c2 = [FeOH+] | 0 mol/m3 | 0 mol/m3 |
c3 = [Na+] | 493.65 mol/m3 | 0 mol/m3 |
c4 = [Cl-] | 1 mol/m³ | 0 mol/m3 |
c5 = [H+] | 1.5849E-2 mol/m³ | 0 mol/m3 |
c6 = [OH-] | 6.3809E-7 mol/m³ | 0 mol/m3 |
c7 = [R-COOH] | 331.01 mol/m³ | 0 mol/m3 |
c8 = [R-COO-] | 416.71 mol/m³ | 0 mol/m3 |
φl | 0 V | 0 V |
Field variable | Top surface of electrolyte layer (x = 1000 μm) | Bottom of the metal layer (x = 0 μm) |
φ | 0 | 1 |
cM/csolid = [Fe2+]/csolid | 0 | 1 |
c2 = [FeOH+] | 0 mol/m3 | 0 mol/m3 |
c3 = [Na+] | 493.65 mol/m3 | 0 mol/m3 |
c4 = [Cl-] | 1 mol/m³ | 0 mol/m3 |
c5 = [H+] | 1.5849E-2 mol/m³ | 0 mol/m3 |
c6 = [OH-] | 6.3809E-7 mol/m³ | 0 mol/m3 |
c7 = [R-COOH] | 331.01 mol/m³ | 0 mol/m3 |
c8 = [R-COO-] | 416.71 mol/m³ | 0 mol/m3 |
ϕl | 0V | ∂ϕ/∂ϕ∂x=0∂x=0V/m |
Parameter | Description | Value | Unit of measure |
A | Free energy density curvature | 5.35E7 | N/m² |
α | Anodic charge transfer coefficient | 0.5 | - |
csolid | Average concentration of metal | 1.43E5 | mol/m³ |
csat | Average saturation concentration | 5100 mol/m³ | mol/m³ |
DM, Di | Diffusion coefficient | 8.5E-10 | m²/s |
K1 | Chemical equilibrium constant | 1.625E-4 | mol/m³ |
K2 | Chemical equilibrium constant | 1.0113E-8 | mol²/m⁶ |
K3 | Chemical equilibrium constant | 0.028184 | mol/m³ |
k1f | Forward reaction rate | 0.08125 | 1/s |
k2f | Forward reaction rate | 5.0565E-5 | mol/(m³·s) |
k3f | Forward reaction rate | 1409.2 | 1/s |
k1b | Backward reaction rate | 500 | m³/(s·mol) |
k2b | Backward reaction rate | 5000 | m³/(s·mol) |
k3b | Backward reaction rate | 50,000 | m³/(s·mol) |
nM | Average charge number: [Fe2+] | 2.1 | - |
nc2 | Average charge number: [FeOH+] | 1.1 | - |
nc3 | Average charge number: [Na+] | 1 | - |
nc4 | Average charge number: [Cl-] | -1 | - |
nc5 | Average charge number: [H+] | 1 | - |
nc6 | Average charge number: [OH-] | -1 | - |
nc7 | Average charge number: [R-COOH] | 0 | - |
nc8 | Average charge number: [R-COO-] | -1 | - |
T | Temperature | 300 | K |
ω | Height of the double-well potential | 2.4E7 | N/m² |
Rg | Gas constant | 8.314 | J/(mol·K) |
τ | Gradient energy coefficient | 7.5E-5 | N |
θ | Electric conductivity | 1E7 | S/m |
Field variable | Liquid phase | Solid phase |
φ | 0 | 1 |
cM/csolid = [Fe2+]/csolid | 0 | 1 |
c2 = [FeOH+] | 0 mol/m3 | 0 mol/m3 |
c3 = [Na+] | 493.65 mol/m3 | 0 mol/m3 |
c4 = [Cl-] | 1 mol/m³ | 0 mol/m3 |
c5 = [H+] | 1.5849E-2 mol/m³ | 0 mol/m3 |
c6 = [OH-] | 6.3809E-7 mol/m³ | 0 mol/m3 |
c7 = [R-COOH] | 331.01 mol/m³ | 0 mol/m3 |
c8 = [R-COO-] | 416.71 mol/m³ | 0 mol/m3 |
φl | 0 V | 0 V |
Field variable | Top surface of electrolyte layer (x = 1000 μm) | Bottom of the metal layer (x = 0 μm) |
φ | 0 | 1 |
cM/csolid = [Fe2+]/csolid | 0 | 1 |
c2 = [FeOH+] | 0 mol/m3 | 0 mol/m3 |
c3 = [Na+] | 493.65 mol/m3 | 0 mol/m3 |
c4 = [Cl-] | 1 mol/m³ | 0 mol/m3 |
c5 = [H+] | 1.5849E-2 mol/m³ | 0 mol/m3 |
c6 = [OH-] | 6.3809E-7 mol/m³ | 0 mol/m3 |
c7 = [R-COOH] | 331.01 mol/m³ | 0 mol/m3 |
c8 = [R-COO-] | 416.71 mol/m³ | 0 mol/m3 |
ϕl | 0V | ∂ϕ/∂ϕ∂x=0∂x=0V/m |
Parameter | Description | Value | Unit of measure |
A | Free energy density curvature | 5.35E7 | N/m² |
α | Anodic charge transfer coefficient | 0.5 | - |
csolid | Average concentration of metal | 1.43E5 | mol/m³ |
csat | Average saturation concentration | 5100 mol/m³ | mol/m³ |
DM, Di | Diffusion coefficient | 8.5E-10 | m²/s |
K1 | Chemical equilibrium constant | 1.625E-4 | mol/m³ |
K2 | Chemical equilibrium constant | 1.0113E-8 | mol²/m⁶ |
K3 | Chemical equilibrium constant | 0.028184 | mol/m³ |
k1f | Forward reaction rate | 0.08125 | 1/s |
k2f | Forward reaction rate | 5.0565E-5 | mol/(m³·s) |
k3f | Forward reaction rate | 1409.2 | 1/s |
k1b | Backward reaction rate | 500 | m³/(s·mol) |
k2b | Backward reaction rate | 5000 | m³/(s·mol) |
k3b | Backward reaction rate | 50,000 | m³/(s·mol) |
nM | Average charge number: [Fe2+] | 2.1 | - |
nc2 | Average charge number: [FeOH+] | 1.1 | - |
nc3 | Average charge number: [Na+] | 1 | - |
nc4 | Average charge number: [Cl-] | -1 | - |
nc5 | Average charge number: [H+] | 1 | - |
nc6 | Average charge number: [OH-] | -1 | - |
nc7 | Average charge number: [R-COOH] | 0 | - |
nc8 | Average charge number: [R-COO-] | -1 | - |
T | Temperature | 300 | K |
ω | Height of the double-well potential | 2.4E7 | N/m² |
Rg | Gas constant | 8.314 | J/(mol·K) |
τ | Gradient energy coefficient | 7.5E-5 | N |
θ | Electric conductivity | 1E7 | S/m |