Next Article in Journal
Influence of 0.5% Ag Addition on Low-Cycle Fatigue Behavior of Hot-Extruded Al-5Cu-0.8Mg-0.15Zr-0.2Sc Alloy Subjected to Peak-Aging Treatment
Next Article in Special Issue
Effect of Nitrogen on the Corrosion Resistance of 6Mo Super Austenitic Stainless Steel
Previous Article in Journal
Effect of Magnet Alternate Stirring on the Internal Quality of Sn-Pb Alloy
Previous Article in Special Issue
Probing Localised Corrosion Inhibition of AA2024-T3 by Integrating Electrode Array, SVET, SECM, and SEM-EDS Techniques
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Experimental Investigation and Modeling: Considerations of Simultaneous Surface Steel Droplets’ Evaporation and Corrosion

by
Marius Ciprian Ilie
1,2,
Timur Vasile Chiş
3,
Ioana Maior
4,
Cristian Eugen Răducanu
1,*,
Iuliana Mihaela Deleanu
1,
Tănase Dobre
1,5,* and
Oana Cristina Pârvulescu
1
1
Chemical and Biochemical Engineering Department, National University of Science and Technology Politehnica Bucharest, 1-7 Gheorghe Polizu St., 011061 Bucharest, Romania
2
TotalEnergies Marketing Romania SA, 4 Vasile Alecsandri St., 010639 Bucharest, Romania
3
Doctoral School Associate, Petroleum-Gas University of Ploiești, 39 Bucharest Blvd., 100680 Ploiesti, Romania
4
Inorganic Chemistry, Physical Chemistry and Electrochemistry Department, National University of Science and Technology Politehnica Bucharest, 1-7 Gheorghe Polizu St., 011061 Bucharest, Romania
5
Technical Sciences Academy of Romania, 26 Dacia Boulevard, 030167 Bucharest, Romania
*
Authors to whom correspondence should be addressed.
Metals 2023, 13(10), 1733; https://doi.org/10.3390/met13101733
Submission received: 30 August 2023 / Revised: 5 October 2023 / Accepted: 11 October 2023 / Published: 12 October 2023

Abstract

:
The present work focuses on the problem of steel surface corrosion as a kinetic expression when water droplets are repeatedly deposited and evaporated on/from its surface. This process, together with the rainwater film corrosion process, belongs to the theoretical foundations of the problem of atmospheric corrosion. It was considered that the formation of water droplets on surfaces is a random but repetitive process, as well as the fact that experimental and theoretical observations show that the droplet corrosion front of a metal surface is located in its zone circumference. We thus aimed to establish how the corrosion process evolves on a steel plate when many drops are deposited and removed repeatedly. An experimental setup and working procedure were used to obtain data characterizing the simultaneous process of steel surface corrosion and water droplet evaporation. For natural convection conditions with a variable relative humidity and temperature environment, an extensive data set consisting of the dynamics of individual droplet evaporation coupled simultaneously with the corrosion of the steel surface under the droplet was obtained. The mathematical models for evaporation and corrosion under the droplet have the same dynamic transfer surface for water evaporation and oxygen supply in the droplet. An approach for determining this surface depending on the momentary droplet mass was considered. Several simultaneous measurements of evaporation–corrosion dynamics were used to calibrate the coupled models, which were then used to show their compatibility with experimental data.

1. Introduction

Metal losses, especially steel losses, caused by the corrosion of metal structures located in the atmospheric environment are impressive [1,2]. The latest estimation of World Corrosion Organization (WCO) considers that the direct cost of worldwide corrosion is about of EUR 1.3 and 1.4 trillion, which is equivalent to 3.8% of the global Gross Domestic Product (GDP) [3]. Public data on metal corrosion costs show that it is an average of 3.4% of the global GDP, distributed by economic sectors, of which the industrial sector is at a level of 2% (USD 1600 billion out of a GDP of USD 80,000 billion). It is known that the surface chemical process due to oxygen brought from the air into the water that comes periodically as rain or as droplets on the metal surface structures cannot be interrupted. However, it can be controlled in such a way that the operation of metal structures in an atmospheric environment is as safe as possible. The atmospheric corrosion of steels is a complex process with evolution at the molecular scale which is influenced by an extremely large number of factors, starting with the steel structure, the geographical positioning of the steel structure, the design and installation of the metal structure, the frequency and flow rate of the rains, the frequency and intensity of the deposition of drops on the surface, the presence of ionizable pollutants in the water, the non/existence of an anti-corrosion protection solution, etc. In the approach to atmospheric corrosion of weathering steels, there is a quasi-recognized consideration regarding the expression of the thickness of the corroded steel [4,5]. It starts from the acceptance of the simplest model of corrosion, according to which water forms a diffusion film to move towards the surface of the steel where it reaches and causes the maintenance of anodic and cathodic processes as a surface reaction [6]. The result of this model can be seen in Equation (1) which shows the time dependence, expressed in years, of the thickness of the corroded steel (δst), which, as a convention, is given in μm. The experimentally determinable constants A and n take into account the influence of the factors mentioned above on the dynamics A as well as the superficial formation of the rust layer, which reduces the intensity of corrosion n.
δ s t = A   τ n
Improvements have been made to Equation (1) to take into account excess moisture in the environment, the presence of sulfur dioxide (acid rain) or the presence of Cl ions in the corrosion environment [7]. We show that most countries in the world have participated, since 1980, in a program to monitor the atmospheric corrosion of steels, so that data on the parameters of this model are available [8]. Now, the problem of atmospheric corrosion of steels is treated much more complexly, starting mainly from the fact that water reaches the steel structures through rain and through condensation in drops (dew), which determines two main cases of corrosion, namely film corrosion and corrosion drops. Both cases of atmospheric corrosion can be characterized by approaches which consider them to be cases of simultaneous transfer of momentum and mass [9,10], respectively, and of simultaneous transfer of heat and mass [11,12] which are associated with the electrochemical process that takes place at the corrosion surface. In the case of droplet corrosion, their appearance on the condensation surface is a relatively fast process, under 10–15 min [13], if the condensation conditions are met [14], as is shown by Equation (2), where ps is the water saturation pressure, ts and tg represent the surface and air temperature, respectively, Ur gives the air relative humidity and tdew is the dew point temperature.
  p s t s p s t g U r   or   t s t d e w U r
Being a condensation process, its speed is controlled by the heat transfer coefficient, which for condensation has extremely high values (5000–7000 w/m2·grd). A coalescence of the small droplets on the steel surface, which is also fast, leads to a droplet size distribution estimated to follow a normal distribution, between 0.2 mm and 6 mm. Evaporation of a droplet from a surface is a simultaneous mass and heat transfer process with its duration depending on droplet size, on temperature difference ts-tg and on the fact that, here, the heat transfer coefficient is low (below 100 w/m2·grd), corresponding to weak natural convection or laminar gas flow along the evaporation surface. In other words, considering its duration in relation to the surface corrosion process, droplet evaporation is much more important than deposition.
Considering these remarks, we focused in this study on the analysis of the kinetics within the corrosion process of the droplets, upon their evaporation. Many of the current approaches to modeling droplet corrosion refer to the Evans droplet model [12,15,16], which, in its current form [15], characterizes the onset of hemispherical droplet corrosion using unsteady-state species concentration field equations with reactions occurring on the surface at the cathodic and anodic (droplet-based) sites. For small hemispherical drops, the absence of movement in the droplet can be accepted [17], while for large drops this fact is no longer acceptable because here the evaporation of the droplet is accompanied by the Marangoni effect [18,19,20]. For small droplet corrosion, the COMSOL solution of the Evans model [17], where it uses only an anodic site on the surface and a cathodic site near the droplet hemisphere, shows both that the current density curves inside the droplet are arcs, starting from the anode and closing at the cathodic site, and that here, in the cathodic site (the area where the drop closes to the surface), the concentration of dissolved species reaches maximum values. With the same limit conditions, the mentioned model shows about the same result, and, in the case of large drops, the reports are limited to a duration of 600 s. Additionally, in this case, it is shown that the rust is deposited in the marginal area of the drop, while the photographs of the situation after the evaporation of the drop show the imprint of the rust all over its original location [16]. It should be noted that the model in question accepts the hemispherical shape of the droplet [16,17] which implies a wetting angle of the water surface of 90 degrees.
Depending on the metal surface, this angle is between 73 and 84 degrees [21,22], with about 80 degrees for water–steel. For large droplets with a hemispheric diameter in the range of 0.5–5 mm, formed by dew deposition and coalescence on the surface of steel or other metal, which corrode with simultaneous evaporation, there are many problems that are still the focus of research. Thus, it is not clear whether the Marangoni flow inside the droplet causes it to flatten into a cylinder or whether during evaporation the base of the droplet remains unchanged, grows or shrinks. In addition, there are few data on the average corrosion rate during droplet evaporation and how it is influenced by the presence of corrosion accelerators (Cl for example). Given the repeatability of surface dew deposition, it must be shown whether repeated formation of the droplet on a rust trace accelerates or decelerates corrosion. The current paper focuses on these issues or, more correctly, on some aspects of these issues. Specifically, a new approach was used to experimentally investigate the simultaneous process of droplet evaporation and sub-droplet corrosion on the steel surface, the purpose of which was to validate the models developed to describe this phenomenon. The droplet evaporation model and the sub-droplet corrosion model are coupled by having the same transfer surface, which decreases with time. A new solution was thus proposed to express the dynamics of the transfer surface for evaporation, and for oxygen supply in the droplet as a function of the momentary droplet mass, respectively, while several new elements were considered for the two coupled models, the most significant being the decreasing circulation due to the Marangoni phenomenon.

2. Materials and Methods

2.1. Experimental Setup and Procedures

The corrosion processes from one drop, complemented with total Fe+2 to Fe+3 oxidation and with metal oxide precipitation, is expressed [8,9,10] by Equation (3):
2 F e s + 3 2 O 2   l + z + 3 H 2 O 2 F e O H 3 z H 2 O F e 2 O 3 z H 2 O + 3 H 2 O .
If we consider a sufficiently large number of drops deposited on a steel surface, then, at the end of their evaporation, the mass of the plate will increase. This increase in mass allows the determination of the corrosion rate (Equation (4)) and, by relating to the droplet surface, the specific corrosion flow rate (Equation (5)). Herem is the increase in plate mass, τ corresponds to the droplet evaporation time, Ap represents the average droplet corrosion surface and using Mi (I = Fe, Fe2O3 ·zH2O) gives the molecular mass of species i.
G m = Δ m τ 2 M F e M F e 2 O 3 z H 20
N F e = G m n A p
Tracking the dynamics of droplet evaporation brings information about the kinetics of this process related to the external environmental conditions, mainly its temperature and relative humidity.
Repeating this process several times can show the extent to which the rust deposited on the surface is influencing the process. Replacing the water droplets with NaCl-containing droplets will show how the corrosion process is accelerated. The experimental setup is shown in Figure 1, and we mention the precision balance (Kern, Germany) coupling via USB, Kern Software SCD—4.2.4 Pro (Kern, Germany) Temperature and Humidity Data Logger (USB connetion). The weathering steel plate was 2 mm thick, 200 mm wide and 300 mm long. The composition of the steel is given in Table 1. The water used as the corrosion medium had an electrical conductivity of 100 mS/cm) and a pH of ~6.7, which makes it close to clean rainwater. In the accelerated corrosion tests, the NaCl content in the water was 1 g/L.
The atmospheric experiment parameters had values located in different ranges since the cumulative duration of the experiments was more than 9700 min: 18–29 °C for temperature, 27–60% for relative humidity and 6–17 °C for dew point.
The working sequences for a corrosion test include:
  • Determining the mass of the steel plate;
  • Loading the plate with droplets in the preset positions;
  • Starting the recording the momentary mass of the plate with droplets and the state of the moist air parameters near the surface (at ~8 cm from it);
  • Taking digital photographs of the plate with droplets at test starting, during the test and at the end of each test experiment;
  • Ascertaining the drying of the droplets and saving the recordings in order to process them.
Table 2 contains further data for a complete description of the experimental corrosion tests.
If we refer to a droplet corrosion test, we have, experimentally, the dynamics of water evaporation from a droplet, the dynamics of the moisture content near the droplet, photo images that show the shape of the droplet during evaporation as well as the pressure, temperature and moisture content from the air from the evaporation medium. At the end of the test, the mass of rust associated with the mass plate, and therefore associated with each droplet, is available. The entire plan of the experimental investigation considered that its development takes place under natural convection conditions. Thus, each of the 30 tests reported in the paper were performed under conditions in the presence of moisture content and temperature different to the external air environment, recorded for each experiment. Then, each experiment started from the existence of a previous layer of rust on the surface. The duration between two experiments was, as shown in the results tables, random, but always allowed the newly deposited rust layer to reach the structural state of the previous one. The effective time of an experiment, strongly dependent on air humidity and temperature, was slightly modified, as is shown in the results from the data tables in the results section, by the mass of the droplet placed in the fixed location.

2.2. Mathematical Modeling

Modelling of mass and heat transfer phenomena is an effective tool in the design, control and optimization of chemical processes as well as in understanding their mechanisms [23,24]. As a descriptive model, droplet evaporation from a steel surface is a simultaneous mass and heat transfer problem. Here, the water vapor leaves the surface of the droplet in the adjacent environment (mass transfer) and the heat required for this process is brought from the environment (heat transfer). Under the experimental conditions of our test, the mass and heat transport mechanisms correspond to natural convection. The peculiarity of this problem comes from the fact that the drop changes its shape and volume. This fact additionally acts on the corrosion process by changing the specific flow of oxygen transferred to the corrosion surface, and by changing the surface under corrosion, respectively, and, consequently, by reducing rust deposits. Figure 2 shows three possible models for the evolution of the droplet shape in the evaporation process.
Regardless of the shape, if it is identifiable, as shown in Figure 3, then its analytical expression (y = f(x), x = φ(y)) becomes possible from the knowledge of the instantaneous volume [22], as shown by Equation (6). The surface produced by rotating y = f(x) and x = φ(y) with respect to the x and y axes, respectively, [25] is the surface of the air drop interphase (Equation (7)).
V τ 0 x = π a τ a τ f x 2 d x V τ o y = π 0 h τ φ y 2 d y
S τ o x = 2 π a τ a τ f x 1 + f I x 2           S τ 0 y = 2 π 0 h τ φ y 1 + φ I y 2   d y
Our observations, justified with the presentation of the results, support that, for the evaporation of droplets on a steel surface, the initial shape and its evolution correspond to a paraboloid (Equation (8)). Here R is the initial radius of the droplet deposited on the surface and hτ shows the distance from the center of the droplet to its boundary.
x 2 R 2 + y 2 h 2 1 = 0
The assumption that the droplet evaporates keeping its original trail allowed the analysis by modeling the two processes in the droplet. Thus, Figure 4 shows the simultaneous phenomenology of mass and heat transfer during droplet evaporation. Figure 5 shows the transfer of O2 inside the droplet, and its consumption through the reaction on the surface, resulting in the generation of rust. We mention that inside the droplet there is an internal flow determined by the association of the Marangoni phenomenon with its evaporation [25,26,27]. Figure 5 schematically represents this flow. When falling outside, over the boundary layer, the experimental conditions characterize an environment, at most, in natural convection. Under these conditions, the vapors leaving the surface of the drop pass by diffusion over the inert (air) through the boundary layer of thickness δ. On any direction, z, normal to the surface of the drop, the specific flow rate is given by the differential expression in Equation (9), where c is the total molar concentration and Dv represents the air vapor diffusion coefficient.
N v = c D v 1 y v d y v d z
Figure 4 shows the vapor mole fraction at the droplet surface and outside the boundary layer, so that the integration of Equation (9) gives the integral expression of the vapor flow (Equation (10)).
N v = c D v δ y B m y v t s y v δ
Relation (11) is written to highlight the driving force and the mass transfer coefficient in the form (12). Since c, Dv, δ and even yBm in (10) depend on tg and ts, or more correctly on those differences, it was considered, for kg, as dependent on the heat transfer driving force tgts. So,
k g = D v δ y B m = k g t g t s
N v = c k g t g t s y v t s y v δ .
A balance for the duration from τ to τ + , when denoting the interphase transfer surface (evaporation surface) S(τ), yields Equation (12). Here S(τ) depends on the instantaneous mass of the droplet (mp) as follows from the coupling of Equations (6)–(8), where we put Vτ0y = m/ρa. With this consideration, Equation (13) can be written as Equation (14).
d m p d τ = N v M w S τ
d m p d τ = α f m p N v M w
If we had expressions or values for the mass transfer coefficient for the interface temperature as well as for the function f(mp), then the momentary mass of the droplet becomes analytically expressible. Table 3 shows the calculation of the transfer surface, for the droplet masses of interest in the present work, so that α and f(m) from Equation (14) can be identified. In more detail from Equation (8) it is expressed as x = φ(y) and since R is fixed and mp is chosen, in the range of interest of our droplet mass, then solving the equation m p = ρ w V τ o y = π ρ v 0 h τ φ y 2 d y (Equation (6)) gives the h τ , value, which then immediately leads to S or S(τ), as is shown in Equation (7), where it is correlated by regression polynomial after the momentary particle.
To be functional, the model, in terms of the mass dynamics of the evaporating water droplet at the surface, needs a solution to calculate the surface temperature of the water droplet and also to account for the mass transfer coefficient from the droplet to the surrounding air. In this sense, if the heating of the droplet is neglected, then the temperature ts results from the fact that the specific heat flow brought by convection to the droplet (Equation (15)) is equal to the specific heat flow due to vaporization (Equation (16)). It is identified in Equation (15) that the value of the heat transfer coefficient is a function of the air temperature difference and the surface drop. So αg = αg (tgts). In Equation (16) the latent heat of vaporization of water was denoted by rv.
q g = α g   t g t s = α g t g t s   t g t s
q v = N v M w r v
The mentioned equality of qg and qv, coupled with Nv relationship, occurs for ts in Equation (17).
α g t g t s   t g t s = c k g t g t s y v t s y v δ M w r v
The transfer of mass and heat, respectively, from and to the droplet occurs through the boundary layer of the droplet, as shown in Figure 4. Thus, the analogy of transfer phenomena [28,29] can be used to express kg (tgts) using αg (tgts). This is given by Equation (18), where ρg is the air density and cpg represents its specific heat coefficient. Coupling Equation (18) with Equation (19) leads to the expression of the droplet surface temperature by the conditional Equation (19). For yv (ts) Equation (20) is used, where p is the air pressure and A, B, C are the Antoine constants for expressing the water saturation vapor pressure.
k g t g t s = α g t g t s   ρ g c p g
t s = t g c ρ g c p g y v t s y v δ M w r v   ,       y v t s > y v δ
y v t s = 10 A B t s + C p
Regarding Equation (19), it should be specified that y depends on the relative air humidity (φ) and the saturation pressure of water vapor at its temperature tg. It should also be said that, after ts, this relationship is a transcendent equation so to analytically or graphically raise the dependence ts = ts (φ,tg), a calculation program is required.
To complete the model of droplet evaporation by natural convection, i.e., the model showing how the droplet mass evolves over time as a function of initial radius, shape, air relative humidity and air temperature, a relation expressing values of the heat transfer coefficient for this heat transfer mechanism is needed. A relation from the literature [30,31] was thus chosen, here specified in the form of Equation (21), where the coefficients a, m and lc have the specifications shown in Table 4.
α g = α g t g t s = a t g t s l c m
For the modeling of the corrosion process at the droplet–steel interface, it is considered to progress more slowly than the evaporation process. Thus, it may be close to steady state evolution. Therefore, it can be considered to be quasi-stationary. According to Figure 2, the transfer of oxygen through the droplet, from the concentration of cO2d to the concentration of cO2s, corresponding to the surface reaction, is considered the most important. Relation (22) expresses the specific flow rate of oxygen. This becomes (23) by expressing cO2s from the equality of oxygen transferred with oxygen consumed by the reaction surface.
N O 2 l = N O 2 = k l c O 2 d c O 2 s
N O 2 l = k l k r s a k l + k r s a c O 2 d
In Equation (23), cO2d can be considered as the equilibrium of concentration of oxygen and the droplet surface temperature. The mass transfer coefficient kl, given by adapting the literature [32,33], takes into account the fact that the driving force for the Marangoni flow inside the droplet is represented by the surface tension difference between the water at the droplet surface, σ(ts), and at the solid surface, (σ(tp) ≈ σ(tg), as shown by Equation (24). For the apparent reaction constant, krsa, Equation (25) [10] was considered. Here the value of 2.3 × 10−5 m/s [10] was assumed for krs (oxygen surface reaction constant when the steel surface is not rusted (i.e., new, this is a novelty regarding the model from [10] previously)). In regard to the surface reaction rate constant, we show that it has a rather complex meaning. Thus, seeing the corrosion process as a process with a chemical kinetic in which the solid reacts with the limiting reactant (oxygen in water in the case of pure corrosion) and then it is influenced by local surface structure issues (local crystallinity, intra-granular inclusions). As these structural irregularities are distributed, the surface is seen to have average properties. As a result, the reaction rate constant of the surface corrosion process refers to a surface seen with average unevenness. The rust thickness, δr, is linearly dependent on the mass of rust deposited at the solid interface of the droplet (Equation (26)). For the oxygen diffusion coefficient through the rust layer, the range of values is 0.2 × 10−9–2 × 10−9 m2/s [10].
k l = a σ t g σ t s η = 3.75 × 10 5 σ t g σ t s η
1 k r s a = 1 k r s + δ r D O 2 e f
δ r = m r u π R 0 2 ρ r u
The momentary flow rate of consumed oxygen is given by Equation (27) where S(τ) can be expressed as shown in Table 1. Based on Equation (3), the mass of iron dissolved under the droplet (Equation (28),) and the mass of deposited rust under the droplet are obtained, respectively, as is the contour of its trail (Equation (28)). As shown above, the rust mass was measured for each experimental test. Its dynamics, resulting from the measurements performed, can be used to validate the developed models, as well as to promote a global kinetic relationship of droplet corrosion, as shown by Equation (1).
d m O 2 d τ = S τ N O 2 l M O 2 = α f m p N O 2 l M O 2 ,           τ = 0     m O 2 = 0
m F e = 4 3 m O 2 M F e M O 2
m r u = 2 3 m O 2 M r u M O 2
It is worth noting that the droplet evaporation and surface corrosion models are coupled in that the values of variables and parameters from one model are found in the other. We thus have a unitary model of evaporation and corrosion in the droplet. Extending this model to other cases of external droplet flow involves the use of new characteristic relations for the heat and mass transfer coefficients for droplet evaporation.

3. Results

For droplets deposited on weathering steel plates we show that the experimental investigation followed droplet evaporation and steel corrosion by measuring:
  • The water mass dynamics of the drops from the plate;
  • The state of the evaporating drops’ shape;
  • The air parameters near to the plate with the evaporating drops;
  • The increase in plate mass due to rust deposition from corrosion process.
Table 2 shows that in each test 100 drops were deposited on the steel plate. The aggregate droplets deposited on the plate at the beginning of each test were characterized, as shown in Figure 6, by mean diameter and standard deviation. It is also not difficult to appreciate from Figure 6 the paraboloid shape of these drops. The same procedures were used for droplet traces, which remain on the plate after evaporation. It was found, as can be seen in the data below, that the diameter of the trace is between 5 and 10% larger than that of the drop.
Figure 7 and Figure 8 show frames from a filmed experiment step of the corrosion evaporation tests. Figure 7 presents the steps from blank steel plate to the end of droplet evaporation, allowing us to observe the evolution from a water droplet to a shape on the steel plate in a considered first such test experiment. Frames shown in Figure 8 present the same evolution, but for a third considered experiment, which starts with a rusted shape on the steel plate already. Both Figure 7 and Figure 8 come to support that, during evaporation, the droplet shrinks while keeping its diameter very close, practically the same, to the original diameter of the droplet. It is stated in Table 3, that in the building of functions f (mp) this fact, now supported experimentally, was used: namely that the drop shrinks while keeping its initial diameter.
Considering the conclusion after analyzing Figure 7 and Figure 8 we can refer to the differences between the initial diameter of the droplet, established in Figure 6, and the diameter of the droplet trail after the droplet has evaporated. So, this difference in mean diameter and standard deviation shows that there is an extension of preferential rust deposition at the edge of the drop contour, or that the deposition of a new drop into an old trace was not centered enough. Referring to the Evans corrosion model [12], we appreciate that there was a tendency for rust at the edge of the drop contour.
Table 5, Table 6, Table 7 and Table 8 and Figure 9, Figure 10, Figure 11 and Figure 12 present primary data for four evaporation–corrosion tests when the plate surface is made of steel, with its composition shown in Table 1. It is noted that the presence of NaCl in the corrosion medium (column F from Table 7 and Table 8 compared to the same column from Table 4 and Table 5) significantly changes the mass of rust obtained at the end of the test. If we return to the corrosion pattern, it shows, first of all, a significant change in the value of the surface reaction rate constant, which expresses Fe dissolution. Doubling or tripling of the corrosion rate in the presence of Cl is reported in many works [34,35], so our data are in agreement with them.
From Table 5, Table 6, Table 7 and Table 8, and the corresponding figures, supported by all 30 evaporation-corrosion tests in the work, it can be seen that the size of the initial drop, the relative humidity and the temperature of air near the steel surface and the rust loading of location for droplets deposition are important factors that determine the evaporation time and the corrosion rate when the droplets deposited on the surface are similar to rainwater. The chlorine and sodium ions’ presence in drops, in this case at a concentration level of 1 g/L NaCl, does not change the dynamics of evaporation instead, as previously mentioned, instead it causes the increasing of corrosion specific flow rate from 0.01–0.012 mgFe·cm−2h−1 (Table 5 and Table 6, column F) to 0.05–0.06 mgFe·cm−2h−1 (Table 7 and Table 8, column F).
Overall, these tables and figures (each of them) contain all the necessary data, (evaporated water, and the deposited rust mass as dynamics relatable to a single drop, respectively) so that they can be used in testing our developed model.

4. Discussion

In a comprehensive sense, model testing has a calibration part and a validation part. The calibration part involved several sequences, namely:
  • Completing the model with additional data (temperature dependence of dissolved oxygen concentration in water droplet, temperature dependence of water surface tension, densities, molecular masses, etc.) and conditions’ initials;
  • The numerical transposition of the model with the micro sequences: (a) the choice of the parameters of the model that require calibration, namely the coefficient α and the power m in relation (21), the coefficient a in Equation (24), krs and the relation DO2efin (Equation (25)), respectively; (b) selection, from the beginning of the investigation, of the tests with all their data, which are used in the calibration (Table 5 and Table 6 for corrosion in water droplets, respectively, and Table 8 for corrosion with droplets containing NaCl); (c) setting an option regarding air parameters’ (relative humidity and temperature according to Table 5, Table 6, Table 7 and Table 8) use in the model, i.e., as functions of time or as average values; (d) the effective expression of the numerical model from the import of the test data file to the Runge–Kutta integration of the differential equations that provides the dynamics of the mass of the droplet and the dynamics of the mass of rust deposited under the droplet, respectively;
  • The effective use of the numerical model in order to establish values for certain parameters, and strategies for expressing others, respectively.
Table 9 summarizes the model calibration results. Referring to what is presented here, including the reporting of these data by those in the specialized literature, there are some comments to make. The use of mean values for relative humidity (φ) and air temperature (tg) near the plate is supported by the recorded data (Table 5, Table 6, Table 7 and Table 8, column G and H), which show extremely small time variations. The values identified for the surface reaction rate constant (krs) and for the effective diffusion coefficient of oxygen through the rust layer (DO2ef) are consistent with those from film corrosion [9] when evaporation–corrosion occurs in water droplets. The higher values of krsin corrosion by NaCl water droplets show the intensification of the anodic reaction on the steel surface, as reported in other works [34,35,36]. The high value of DO2efin in this case indicates a fairly permeable structure of the rust formed in the corrosion process, as well as that here the transport of oxygen is facilitated by the action of the Cl ions. For the constant α in Equation (24) the model worked well with a higher value compared to those found in the literature (3.75 × 10−5).
If to the presentation of Equation (21) it is added that, in addition to heat transport by natural convection, we have an additional heat transport via vapors leaving the surface, then the value chosen for the constant m is supported as well as the fact that a can be a function of φ and tg.
With the calibrated model, all the corrosion tests performed were simulated with the constant a, from Equation (21), being modified to obtain the best model reproduction of the experimental dynamics showing the mass of the evaporating droplet. For the dynamics of rust deposited during droplet evaporation, we monitored whether, at the end of each test simulation, the rust mass was close to the recorded one (see Table 5, Table 6, Table 7 and Table 8 column F scaled to 100 (number of drops)).
The result of these simulations is concentrated in Figure 13, Figure 14, Figure 15, Figure 16 and Figure 17 and Table 10, which characterize, in comparison with those specified above, each individual test.
The relative deviation between the experimental values and those according to the model, presented in Table 10 for the dynamics of the droplet mass and for the dynamics of the deposited rust mass, respectively, was calculated according to the available data, as it appears in Equations (30) and (31), respectively. In Equation (30) ni is the number of time steps since model integration for the i evaporation-corrosion test.
ε m p M i = 1 n i j = 1 n i m p e x τ j j = 1 n i m p t τ j j = 1 n i m p e x τ j       100             i = 1 ,   2 30
ε m r u M i = m r u e x i m r u t i m r u e x i     100           i = 1 , 2 30
Before making a more detailed presentation of what Figure 13, Figure 14, Figure 15, Figure 16 and Figure 17 and Table 10 show, we can state that we have here a good coverage of the experimental data from those produced by the developed and calibrated model. We want to highlight from the model calibration, and even motivated by it, that for a coefficient in Equation (21) a dependence on the relative humidity of the air and its temperature is expected. We find this dependence in Table 10, which shows that at high relative humidity and temperature, above 27 °C is high (tests 1, 2, 26, 30), and that, at low relative air humidity and temperature, around 24 °C a is the lowest values (tests 10, 13, 14, 21, 27). Calculation of correlation coefficients for a vs. φ and a vs. tg/20 finds the values 0.7551 and −0.2900, respectively, showing that a can, in the limit, be linearly related to φ and can be considered independent of tg. Figure 18 supports these results. The line a vs. φ in this figure is given by Equation (32). Adding Equation (32) to Table 9 means that all parameters of the evaporation–corrosion model are known.
a = 2.423 + 8.839 φ
The dynamics of the evaporation of a single droplet, having variable mass under variable environmental conditions (Table 10), are described well, even very well, by the developed model, as shown by the values of the relative deviations experiment model (Table 10, column εmpM) together with graphical representations from Figure 13, Figure 14 and Figure 15.
Considering the model curve and the experimental curve for each evaporation test (Figure 13, Figure 14 and Figure 15), we note their good correspondence, which supports the quality of the model. From 30 sets of curves, only in the case of test 27, and maybe test 25, does it appear that the model–experiment agreement is borderline. Figure 13, Figure 14 and Figure 15, as well as Figure 16 and Figure 17, in conjunction with the data in Table 10, show that the evaporation time of the drop depends on its size and the relative humidity and temperature of the air. The plots in Figure 17 and Figure 18b, showing the dynamics of the rust mass from each of the 30 experiments, were drawn using the data provided by the model. The slope of these lines, which represent the corrosion rate, is almost the same for tests 1–27 where the water droplet contained no added ions. The nearly threefold increase in this slope when we have 1 g/L NaCl in the droplet (tests 28, 29 and 30) shows the intensification of the corrosion rate by suppressing the cathodic reactions in the corrosion pile. The representation in Figure 18, on the right side, shows the strong agreement between the dynamics of the rust mass deposited under the drop, obtained experimentally and by modeling. Thus, the approach of modeling the corrosion process through a droplet oxygen transfer model with chemical reaction on the steel surface is very well supported. Here too we can observe the above-mentioned change in the evolution of the slope when the composition of the droplet is changed by the addition of ions.
Looking at some possibilities for extending the evaporation–corrosion model, we mention several situations. Given that dewdrops deposited on a surface form in a distribution of sizes and dimensions, the model can simulate cases that are then averaged by the average probability of each size and dimension in the considered distribution. If there are no natural convection conditions during evaporation, then new expressions for the heat transfer coefficient must be added to the model (the case of Equation (21) in the current work). For low temperatures, obviously above the freezing point of water, the model can work, provided that the temperature dependence of the surface reaction rate constant is identified, which may be of the Arrhenius type.

5. Conclusions

Droplet corrosion on the steel surface was analyzed experimentally and by modeling.
A pilot laboratory setup was developed along with the working procedure so that the surface corrosion dynamics of steel could be characterized when droplets of a corrosive medium evaporate from the surface.
The corrosion medium was water with an electrical conductivity and pH close to rainwater, and, for a limited number of water tests, NaCl was added.
In order to analyze the experimental evaporation–corrosion data, a complex mathematical model based on simultaneous heat and mass transfer and simultaneous oxygen transfer with surface reaction was considered.
The expression, in the evaporation–corrosion model, of the evaporation surface for the droplet size was solved by integral analysis with respect to the droplet shape.
It has been shown that the evaporation of the drop occurs with the preservation of the original trace and that the rust resulting from corrosion is deposited on this trace.
In the calibration numerical model, values for the basic parameters of the model were identified, i.e., surface reaction rate constant at 4.5 × 10−5 m/s, oxygen diffusion coefficient through the rust layer at 9.1 × 10−10 m2/s, m constant from the Equation (21) to 0.33, respectively, and α from Equation (24) to 2.63 × 10−4 m2·s (Table 9). When NaCl was present in the corrosion medium, at the concentration level of 1 g/L, the surface reaction rate constant increased 21 times. The same increase was identified for the oxygen diffusion coefficient through the crust.
For all 30 evaporation–corrosion tests, reported in Table 10, a good, even very good, coverage of the experimental results of droplet evaporation and corrosion dynamics was sustained qualitatively by graphic representations from Figure 14, Figure 15, Figure 16 and Figure 17 and quantitatively by identification of an acceptable range for the mean relative deviation experiment model (εmpM r and εmruM, respectively).
Excellent model–experiment agreement was obtained regarding the dynamics of the rust mass associated with a droplet.
In the case of evaporation-droplet corrosion with NaCl content, the involvement of Cl- ions in the anodic corrosion process, led to a strong change in the reaction rate constant and the oxygen diffusion coefficient through the rust layer values, so that we could obtain, for the dynamics of the mass deposited by the rust, a good agreement between the model and the experiment.
It was shown that the developed model can be adapted so that they can serve to simulate surface corrosion when on the droplets have an accepted size distribution or when natural convection is replaced by forced convection.

Author Contributions

Conceptualization, M.C.I., T.D. and C.E.R.; methodology, M.C.I. and T.D.; validation, T.D., O.C.P., I.M., T.V.C. and I.M.D.; formal analysis, T.D., O.C.P., I.M., T.V.C. and I.M.D.; investigation, M.C.I. and T.D.; data curation, M.C.I., T.D. and C.E.R.; writing—original draft preparation, T.D. and M.C.I.; writing—review and editing, O.C.P. and C.E.R.; supervision, T.D. All authors have read and agreed to the published version of the manuscript.

Funding

This work has been funded by the European Social Fund from the Sectoral Operational Programme Human Capital 2014–2020, through the Financial Agreement with the title “Training of PhD students and postdoctoral researchers in order to acquire applied research skills—SMART”, Contract no. 13530/16.06.2022—SMIS code: 153734.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Gerhardus, K.H.; Michel, B.H.P.; Neil, T.G.; Virmani, P.Y.; Payer, J.H. Corrosion Costs and Preventives Strategies in the United States. Available online: https://impact.nace.org/documents/ccsupp.pdf (accessed on 3 August 2023).
  2. Ke, W.; Li, Z. Survey of Corrosion Cost in China and Preventive Strategies. Corros. Sci. Technol. 2008, 7, 259–264. [Google Scholar]
  3. Chico, B.; Fuente, D.; Díaz, J.; Simancas, J.; Morcillo, M. Annual Atmospheric Corrosion of Carbon Steel Worldwide. An Integration of ISOCORRAG, ICP/UNECE and MICAT Databases. Materials 2017, 10, 601. [Google Scholar] [CrossRef] [PubMed]
  4. Hou, W.; Liang, C. Atmospheric corrosion prediction of steels. Corrosion 2004, 60, 313–322. [Google Scholar] [CrossRef]
  5. Díaz, I.; Cano, H.; Chico, B.; De la Fuente, D.; Morcillo, M. Some clarifications regarding literature on atmospheric corrosion of weathering steels. Int. J. Corros. 2012, 2012, 812192. [Google Scholar] [CrossRef]
  6. Forolius, A.Z. The Kinetics of Anodic Dissolution of Iron in High Purity Water. Corros. Eng. 1979, 28, 10–17. [Google Scholar] [CrossRef] [PubMed]
  7. Klinesmith, D.E.; McCuen, R.; Albrecht, H.P. Effect of environmental conditions on corrosion rates. J. Mater. Civil. Eng. 2007, 19, 121–129. [Google Scholar] [CrossRef]
  8. Morcillo, M.; Chico, B.; Díaz, I.; Cano, H.; De la Fuente, D. Atmospheric corrosion data of weathering steels. A review. Corros. Sci. 2013, 77, 6–24. [Google Scholar] [CrossRef]
  9. To, D. Evaluation of Atmospheric Corrosion in Steels for Corrosion Mapping in Asia, Chapter 2. Ph.D. Thesis, Graduate School of Engineering, Yokohama National University, Yokohama, Japan, 2017. [Google Scholar]
  10. Ilie, M.C.; Maior, I.; Raducanu, C.E.; Deleanu, I.M.; Dobre, T.; Parvulescu, O.C. Experimental Investigation and Modeling of Film Flow Corrosion. Metals 2023, 13, 1425. [Google Scholar] [CrossRef]
  11. Tsuru, T.; Tamiya, K.I.; Nishikata, A. Formation and growth of microdroplets during the initial stage of atmospheric corrosion. Electrochim. Acta 2004, 49, 2709–2715. [Google Scholar] [CrossRef]
  12. Sainz-Rosales, A.; Ocampo-Lazcarro, X.; Hernández-Pérez, A.; González-Gutiérrez, A.G.; Larios-Durán, R.E.; Ponce de León, C.; Walsh, C.F.; Bárcena-Soto, M.; Casillas, N. Classic Evans’s Drop Corrosion Experiment Investigated in Terms of a Tertiary Current and Potential Distribution. Corros. Mater. Degrad. 2022, 3, 270–280. [Google Scholar] [CrossRef]
  13. Narhe, R.D.; González-Vinas, W.; Beysens, D.A. Water condensation on zinc surfaces treated by chemical bath deposition. Appl. Surf. Sci. 2010, 256, 4930–4933. [Google Scholar] [CrossRef]
  14. Lawrence, G.M. The relationship between relative humidity and the dewpoint temperature in moist air: A simple conversion and applications. Bull. Am. Meteorol. Soc. 2005, 86, 225–233. [Google Scholar] [CrossRef]
  15. Evans, U.R. Mechanism of rusting. Corros. Sci. 1969, 9, 813–821. [Google Scholar] [CrossRef]
  16. Chen, C.; Mansfeld, F. Potential distribution in the Evans drop experiment. Corros. Sci. 1997, 39, 409–413. [Google Scholar] [CrossRef]
  17. Venkatraman, M.S.; Cole, I.S.; Gunasegaram, D.R.; Emmanuel, B. Modelling Corrosion of a Metal under an Aerosol Droplet. Mater. Sci. Forum 2010, 654–656, 1650–1653. [Google Scholar] [CrossRef]
  18. Tejaswi, J.; Wang, Z.; Askounis, A.; Orejon, D.; Harish, S.; Takata, Y.; Mahapatra, P.S.; Pattamatta, A. Evaporation kinetics of pure water drops: Thermal patterns, Marangoni flow, and interfacial temperature difference. Phys. Rev. E 2018, 98, 052804. [Google Scholar]
  19. Hua, H.; Larson, R.G. Marangoni Effect Reverses Coffee-Ring Depositions. J. Phys. Chem. B 2006, 110, 7090–7096. [Google Scholar] [CrossRef]
  20. Yuri, O.P. Evaporative deposition patterns: Spatial dimensions of the deposit. Phys. Rev. E 2005, 71, 036313. [Google Scholar]
  21. Somlyai-Sipos, L.; Baumli, P. Wettability of Metals by Water. Metals 2022, 12, 1274. [Google Scholar] [CrossRef]
  22. Xu, J.; Li, B.; Lian, J.; Ni, J.; Xiao, J. Wetting behaviors of water droplet on rough metal substrates. Adv. Eng. Res. 2016, 9, 341–345. [Google Scholar]
  23. Dobre, T.; Pârvulescu, O.C.; Stoica-Guzun, A.; Stroescu, M.; Jipa, I.; Al Janabi, A.A. Heat and mass transfer in fixed bed drying of non-deformable porous particles. Int. J. Heat Mass Transf. 2016, 103, 478–485. [Google Scholar] [CrossRef]
  24. Dobre, T.; Zicman, L.R.; Pârvulescu, O.C.; Neacşu, E.; Ciobanu, C.; Drăgolici, F.N. Species removal from aqueous radioactive waste by deep-bed filtration. Environ. Pollut. 2018, 241, 303–310. [Google Scholar] [CrossRef] [PubMed]
  25. Duda, I.; Gradinaru, S. Calculation of areas in parametric coordinates. In Integral Calculus with Applications; Romania of Tomorrow Publishing House: Bucharest, Romania, 2007; pp. 394–452. ISBN 978-973-725-824-3. [Google Scholar]
  26. Xiongjiang, Y.; Jinliang, X. Does sunlight always accelerate water droplet evaporation? Appl. Phys. Lett. 2020, 116, 253903. [Google Scholar]
  27. Van Gaalen, R.T.; Diddens, C.; Wijshoff, H.M.A.; Kuerten, J.G.M. Marangoni circulation in evaporating droplets in the presence of soluble surfactants. J. Colloid Interface Sci. 2021, 584, 622–633. [Google Scholar] [CrossRef] [PubMed]
  28. Girard, F.; Antoni, M.; Sefiane, K. On the effect of Marangoni flow on evaporation rates of heated water drops. Langmuir 24. Am. Chem. Soc. 2008, 17, 9207–9210. [Google Scholar]
  29. Perry, H.D.; Green, W.D.; Malloney, O.J. Perry’s Chemical Engineering Handbook, 7th ed.; McGraw—Hill: San Francisco, CA, USA, 1997; pp. 78–80. [Google Scholar]
  30. Steeman, H.J.; Janssens, A.; De Paepe, M. On the applicability of the heat and mass transfer analogy in indoor air flows. Int. J. Heat Mass Transf. 2009, 52, 1431–1442. [Google Scholar] [CrossRef]
  31. Bird, B.R.; Stewart, E.W.; Lightftoot, N.E. Transport Phenomena, 2nd ed.; John Wiley & Sons: Hoboken, NJ, USA, 2002; pp. 442–443. [Google Scholar]
  32. Leca, A.; Pop, M.G.; Prisecaru, I.; Neaga, C.; Zidaru, G.H.; Musatescu, V.; Isbasoiu, E.C. Calculation Guide. Tables, Nomograms and Thermotechnical Formulas; Technical Publishing House: Bucharest, Romania, 1987; Volume 2, p. 71. [Google Scholar]
  33. Davies, J.T.; Rideal, E.K. Chapter: Mass Transfer across Interfaces. In Interfacial Phenomena; Academic Press: New York, NY, USA, 1963. [Google Scholar]
  34. Dobre, T.; Stoica, A. Mass Transfer Intensification; Electra Publishing House: Florence, Italy, 2017; Chapter 7. [Google Scholar]
  35. Ahmed, S.; Hou, Y.; Lepkova, K.; Pojtanabuntoeng, T. Investigation of the effect chloride ions on carbon steel in closed environments at different temperatures. Corros. Mater. Degrad. 2023, 4, 364–381. [Google Scholar] [CrossRef]
  36. Morcillo, M.; Chico, B.; Alcantara, J.; Diaz, I.; Simancas, J.; De la Fuente, D. Atmospheric corrosion of mild steel in chloride-rich environments. Questions to be answered. Mater. Corros. 2015, 66, 882–892. [Google Scholar] [CrossRef]
Figure 1. Experimental laboratory setup for corrosion in water droplets: (1) black steel plate for corrosion tests; (2) precision balance; (3) water droplets to corrode; (4) balance protection enclosure; (5) temperature humidity dew point Data Logger; (6) laboratory meteorological micro station; (7) filming device; (8) data registration and processing system; (9) stand; (10) clamp stand.
Figure 1. Experimental laboratory setup for corrosion in water droplets: (1) black steel plate for corrosion tests; (2) precision balance; (3) water droplets to corrode; (4) balance protection enclosure; (5) temperature humidity dew point Data Logger; (6) laboratory meteorological micro station; (7) filming device; (8) data registration and processing system; (9) stand; (10) clamp stand.
Metals 13 01733 g001
Figure 2. Possible models for the evolution of droplets shape during their evaporation, hemisphere shape: (A) descending hemisphere; (B) spherical cap with trace preservation; (C) spherical cap without trace preservation. Rust is brown color.
Figure 2. Possible models for the evolution of droplets shape during their evaporation, hemisphere shape: (A) descending hemisphere; (B) spherical cap with trace preservation; (C) spherical cap without trace preservation. Rust is brown color.
Metals 13 01733 g002
Figure 3. Droplet shape in y-x graphic representation.
Figure 3. Droplet shape in y-x graphic representation.
Metals 13 01733 g003
Figure 4. Representation of mass and heat transfer when evaporating the droplet from the steel surface (boundary layer limit as dotted line).
Figure 4. Representation of mass and heat transfer when evaporating the droplet from the steel surface (boundary layer limit as dotted line).
Metals 13 01733 g004
Figure 5. The process of oxygen transfer and corrosion inside the droplet with Marangoni flow (boundary layer limit as dotted line).
Figure 5. The process of oxygen transfer and corrosion inside the droplet with Marangoni flow (boundary layer limit as dotted line).
Metals 13 01733 g005
Figure 6. Determination of the average diameter of the water droplet and its associated standard deviation (Test 1, d0 = 7.71 ± 0.8 mm).
Figure 6. Determination of the average diameter of the water droplet and its associated standard deviation (Test 1, d0 = 7.71 ± 0.8 mm).
Metals 13 01733 g006
Figure 7. Frames from droplet evaporation–corrosion on a steel surface in Test 1, showing that droplet shrinkage occurs with droplet diameter maintained: (φ = 36.5%, d0 = 7.2 mm): (a) blank black steel plate, starting droplet’s deposition; (b) droplet already in position; (c) evaporation evolution after 21 min; (d) end of evaporation and shape formation; (d) bottom-third deposition of the drop, φ = 43.5%, d0 = 7.2 mm.
Figure 7. Frames from droplet evaporation–corrosion on a steel surface in Test 1, showing that droplet shrinkage occurs with droplet diameter maintained: (φ = 36.5%, d0 = 7.2 mm): (a) blank black steel plate, starting droplet’s deposition; (b) droplet already in position; (c) evaporation evolution after 21 min; (d) end of evaporation and shape formation; (d) bottom-third deposition of the drop, φ = 43.5%, d0 = 7.2 mm.
Metals 13 01733 g007
Figure 8. Frames from droplet evaporation–corrosion on a steel surface in Test 3, showing that droplet shrinkage occurs with droplet diameter maintained: (φ = 43.5%, d0 = 7.2 mm): (a) starting droplet’s deposition onto rusted old trails; (b) droplet already in position; (c) evaporation evolution after 21 min; (d) end of evaporation with the new layer of rust.
Figure 8. Frames from droplet evaporation–corrosion on a steel surface in Test 3, showing that droplet shrinkage occurs with droplet diameter maintained: (φ = 43.5%, d0 = 7.2 mm): (a) starting droplet’s deposition onto rusted old trails; (b) droplet already in position; (c) evaporation evolution after 21 min; (d) end of evaporation with the new layer of rust.
Metals 13 01733 g008
Figure 9. Droplet shape and diameter during their evaporation from the steel surface in evaporation-corrosion Test 1: (a) initial; (b) after 120 min; (c) after 210 min (end of test). Mean diameter and standard deviation are shown in Table 5.
Figure 9. Droplet shape and diameter during their evaporation from the steel surface in evaporation-corrosion Test 1: (a) initial; (b) after 120 min; (c) after 210 min (end of test). Mean diameter and standard deviation are shown in Table 5.
Metals 13 01733 g009
Figure 10. Droplet shape and diameter during their evaporation from the steel surface in evaporation-corrosion Test 4: (a) initial; (b) after 120 min; and (c) after 210 min (end of test). Mean diameter and standard deviation are shown in Table 6.
Figure 10. Droplet shape and diameter during their evaporation from the steel surface in evaporation-corrosion Test 4: (a) initial; (b) after 120 min; and (c) after 210 min (end of test). Mean diameter and standard deviation are shown in Table 6.
Metals 13 01733 g010
Figure 11. Droplet shape and diameter during their evaporation from the steel surface in the evaporation-corrosion test 28: (a) initial; (b) after 120 min; (c) after 330 min (end of test). Mean diameter and standard deviation are shown in Table 7).
Figure 11. Droplet shape and diameter during their evaporation from the steel surface in the evaporation-corrosion test 28: (a) initial; (b) after 120 min; (c) after 330 min (end of test). Mean diameter and standard deviation are shown in Table 7).
Metals 13 01733 g011
Figure 12. Droplet shape and diameter during their evaporation from the steel surface in the fourth evaporation-corrosion test 30: (a) initial; (b) after 120 min; (c) after 210 min (end of test). Mean diameter and standard deviation are shown in Table 8).
Figure 12. Droplet shape and diameter during their evaporation from the steel surface in the fourth evaporation-corrosion test 30: (a) initial; (b) after 120 min; (c) after 210 min (end of test). Mean diameter and standard deviation are shown in Table 8).
Metals 13 01733 g012
Figure 13. Dynamics of droplets mass during their evaporation for tests 1 to 10 (T1…T10): (a) T1–red; T3—blue; T5—green; T7—black; T9—brown; (b) T2–red; T4–blue; T6–green; T8–black; T10–brown; continuous curves–model, dashed curves–experimental.
Figure 13. Dynamics of droplets mass during their evaporation for tests 1 to 10 (T1…T10): (a) T1–red; T3—blue; T5—green; T7—black; T9—brown; (b) T2–red; T4–blue; T6–green; T8–black; T10–brown; continuous curves–model, dashed curves–experimental.
Metals 13 01733 g013
Figure 14. Dynamics of droplets mass during their evaporation for tests 11 to 20 (T11…T20): (a) T11—red; T13—blue; T15—green; T17—black; T19—brown; (b) T12—red; T14—blue; T16—green; T18—black; T20—brown; continuous curves—model, dashed curves—experimental.
Figure 14. Dynamics of droplets mass during their evaporation for tests 11 to 20 (T11…T20): (a) T11—red; T13—blue; T15—green; T17—black; T19—brown; (b) T12—red; T14—blue; T16—green; T18—black; T20—brown; continuous curves—model, dashed curves—experimental.
Metals 13 01733 g014
Figure 15. Dynamics of droplets mass during their evaporation for tests 21 to 30 (T21…T30): (a) T21—red; T23—blue; T25—green; T27—black; T29—brown; (b) T22—red; T24—blue; T26—green; T28—black; T30—brown; Continuous curves—model, dashed curves—experimental.
Figure 15. Dynamics of droplets mass during their evaporation for tests 21 to 30 (T21…T30): (a) T21—red; T23—blue; T25—green; T27—black; T29—brown; (b) T22—red; T24—blue; T26—green; T28—black; T30—brown; Continuous curves—model, dashed curves—experimental.
Metals 13 01733 g015
Figure 16. Dynamics of deposited rust mass during droplets evaporation-corrosion, from one droplet, experimental for tests T1…T20: (a) T1 to T10; (b) T11 to T20; continuous curves—model, dashed curves—experimental.
Figure 16. Dynamics of deposited rust mass during droplets evaporation-corrosion, from one droplet, experimental for tests T1…T20: (a) T1 to T10; (b) T11 to T20; continuous curves—model, dashed curves—experimental.
Metals 13 01733 g016
Figure 17. Dynamics of deposited rust mass during droplets evaporation—corrosion: (a) in experimental tests T121…T30; (b) cumulated mass rust vs. effective corrosion time; line—model, dashed line—experimental).
Figure 17. Dynamics of deposited rust mass during droplets evaporation—corrosion: (a) in experimental tests T121…T30; (b) cumulated mass rust vs. effective corrosion time; line—model, dashed line—experimental).
Metals 13 01733 g017
Figure 18. The evolution of a coefficient of Equation (21) in respect to the relative humidity of the air (a) and its temperature (b).
Figure 18. The evolution of a coefficient of Equation (21) in respect to the relative humidity of the air (a) and its temperature (b).
Metals 13 01733 g018
Table 1. Black weathering steel sheet composition used in water droplet corrosion research (according to manufacturer) [10].
Table 1. Black weathering steel sheet composition used in water droplet corrosion research (according to manufacturer) [10].
ElementComposition (wt%)
Manganese (Mn)0.166
Phosphorus (P)0.028
Sulfur (S)0.028
Carbon (C)0.206
Chromium (Cr)0.078
Molybdenum (Mo)0.114
Vanadium (V)0.003
Silicon (Si)0.004
Copper (Cu)0.082
Nickel (Ni)0.088
Titanium (Ti)0.004
Iron (Fe)0.199
Table 2. Experimental data characterizing laboratory droplet corrosion tests.
Table 2. Experimental data characterizing laboratory droplet corrosion tests.
C.N.Experimental FactorsValuesObservations
1Number of drops on the plate100Selected
2The initial mass of the drops on the plate (g)4.5–15.5Selected
3Average droplet volume (μL)45–55Computed
4The initial shape of the dropSpherical capObserved
5Electrical conductivity and water pH (μS/cm)100Selected
6Number of successive tests for corrosion in water26Selected
7NaCl concentration (g/L) in water (accelerated corrosion)1Selected
8Number of accelerated corrosion tests4Selected
Table 3. Calculation of droplet evaporation surface value as a function of droplet mass for the case of paraboloid cap with constant cap radius.
Table 3. Calculation of droplet evaporation surface value as a function of droplet mass for the case of paraboloid cap with constant cap radius.
ParametersSurface Values as Function of Droplet Mass
R (cm)0.30
mp (g)0.0600.0500.0400.0300.0200.0100.005
h or hτ (cm)0.3120.2760.2350.1830.1330.0690.035
S or S(τ) (cm2)0.5880.5280.4530.3780.2860.2830.283
α f ( m p ) α = π R 2     f m p = 0.908 + 7.08 m p + 221.2   m p 2
R (cm)0.35
mp (g)0.0600.0500.0400.0300.0200.0100.005
h or hτ (cm)0.2520.2280.1890.1470.1010.0520.022
S or S(τ) (cm2)0.5950.5300.4580.3850.3850.3850.385
α f ( m p ) α = π R 2     f m p = 1.037 7.79 m p + 277.5   m p 2
R (cm)0.40
mp (g)0.0600.0500.0400.0300.0200.0100.005
h or hτ (cm)0.2170.1850.1520.1160.0780.0400.020
S or S(τ) (cm2)0.6010.5290.5050.5030.5030.5030.503
α f ( m p ) α = π R 2     f m p = 1.043 5.98 m p + 135.7   m p 2
R (cm)0.45
mp (g)0.0600.0500.0400.0300.0200.0100.005
h or hτ (cm)0.1790.1570.1230.0930.0620.0310.016
S or S(τ) (cm2)0.6660.6560.6420.6360.6360.6360.636
α f ( m p ) α = π R 2     f m p = 1.005 0.933 m p + 27.5   m p 2
Table 4. Values for the parameters of Equation (21).
Table 4. Values for the parameters of Equation (21).
C.N.ParameterVertical SurfaceHorizontal Upper SurfaceHorizontal Lower Surface
1 a 1.4201.3201.520
2 m 0.2500.2500.330
3 l c 0.1120.0550.055
Table 5. Experimental data characterizing the evaporation of droplets from the steel surface and its corrosion (Test 1, d0 = 7.78 ± 0.31 mm, du = 8.05 ± 0.45 mm).
Table 5. Experimental data characterizing the evaporation of droplets from the steel surface and its corrosion (Test 1, d0 = 7.78 ± 0.31 mm, du = 8.05 ± 0.45 mm).
ABCDEFGHIJKLM
10482.15487.250.000.0527.753.917.526.055.0757.60.0
230486.350.9027.852.917.30.5
360485.301.0527.953.317.51.0
490484.001.3027.952.817.31.5
5120483.500.5028.052.117.22.0
6150482.800.7027.852.817.22.5
7180482.200.6027.852.917.33.0
8210482.200.0027.853.217.33.5
(A) current number; (B) time (min); (C) initial plate mass (g); (D) current plate mass (g); (E) evaporated water (g); (F) resulted rust mass (g); (G) air temperature (°C); (H) air humidity (%); (I) dew point temperature (°C); (J) external air temperature (°C); (K) external air humidity (%); (L) air environment pressure (mmHg); (M) cumulated time (h).
Table 6. Experimental data characterizing the evaporation of droplets from the steel surface, and its corrosion (Test 4, d0 = 7.79 ± 0.54 mm, du = 8.38 ± 0.55 mm).
Table 6. Experimental data characterizing the evaporation of droplets from the steel surface, and its corrosion (Test 4, d0 = 7.79 ± 0.54 mm, du = 8.38 ± 0.55 mm).
ABCDEFGHIJKLM
10482.30486.850.000.0527.243.213.62.044.075681.0
230484.352.5027.442.813.681.5
360484.100.2527.542.513.682.0
490482.801.3027.741.913.582.5
5120482.350.0027.940.913.383.0
6150482.350.0027.840.913.383.5
(A) current number; (B) time (min); (C) initial plate mass (g); (D) current plate mass (g); (E) evaporated water (g); (F) resulted rust mass (g); (G) air temperature (°C); (H) air humidity (%); (I) dew point temperature (°C); (J) external air temperature (°C); (K) external air humidity (%); (L) air environment pressure (mmHg); (M) cumulated time (h).
Table 7. Experimental data characterizing the evaporation of droplets from the steel surface and its corrosion (Test 28, d0 = 8.70 ± 0.44 mm, du = 8.93 ± 0.50 mm).
Table 7. Experimental data characterizing the evaporation of droplets from the steel surface and its corrosion (Test 28, d0 = 8.70 ± 0.44 mm, du = 8.93 ± 0.50 mm).
ABCDEFGHIJKLM
10484.45501.600.000.2528.532.210.332457599711.5
230498.003.6028.331.79.89712.0
360497.300.7028.233.410.59712.5
490495.401.9028.332.610.39713.0
5120492.353.0528.132.29.99713.5
6150490.202.1528.132.19.99714.0
7180488.401.8028.031.79.69714.5
8210487.650.7528.031.89.69715.0
9240486.301.3527.932.49.89715.5
10270485.400.9027.833.110.19716.0
11300484.850.5527.733.610.29716.5
12330484.850.0027.634.110.39717.0
(A) current number; (B) time (min); (C) initial plate mass (g); (D) current plate mass (g); (E) evaporated water (g); (F) resulted rust mass (g); (G) air temperature (°C); (H) air humidity (%); (I) dew point temperature (°C); (J) external air temperature (°C); (K) external air humidity (%); (L) air environment pressure (mmHg); (M) cumulated time (h).
Table 8. Experimental data characterizing the evaporation of droplets from the steel surface and its corrosion (Test 30, d0 = 9.05 ± 0.47 mm, du = 9.32 ± 0.49 mm).
Table 8. Experimental data characterizing the evaporation of droplets from the steel surface and its corrosion (Test 30, d0 = 9.05 ± 0.47 mm, du = 9.32 ± 0.49 mm).
ABCDEFGHIJKLM
10485.20496.150.000.3028.937.81330357579769.0
230493.103.0529.137.613.19769.5
360491.201.9029.237.713.39770.0
490489.351.8529.036.912.79770.5
5120488.151.2028.837.112.79771.0
6150486.701.4528.737.112.69771.5
7180485.501.2028.637.112.59772.0
8210485.500.0028.536.712.29772.5
(A) current number; (B) time (min); (C) initial plate mass (g); (D) current plate mass (g); (E) evaporated water (g); (F) resulted rust mass (g); (G) air temperature (°C); (H) air humidity (%); (I) dew point temperature (°C); (J) external air temperature (°C); (K) external air humidity (%); (L) air environment pressure (mmHg); (M) cumulated time (h).
Table 9. Values of some model parameters and strategy from Equation (21).
Table 9. Values of some model parameters and strategy from Equation (21).
C.N.CaseDataφtg k r s  
(Equation (25))
D O 2 e f
(Equation (23))
a (Equation (21))m
(Equation (21))
α
(Equation (24))
1Water dropletTable 5
Table 6
meanmean4.5 × 10−5 (m/s)9.1 × 10−10 (m2/s)f(φ,tg)0.332.63 × 10−4 (m2·s)
2Water droplet with NaClTable 8meanmean9.5 × 10−4 (m/s)5.1 × 10−9 (m2/s)f(φ,tg)0.332.63 × 10−4 (m2·s)
Table 10. Corrosion tests and comparison of experimental results with those according to the model by relative deviations (εmpM, εmruM).
Table 10. Corrosion tests and comparison of experimental results with those according to the model by relative deviations (εmpM, εmruM).
Testmp0
(g)
φ
(/)
tg (°C)a
(20)
εmpM
(%)
εmruM
(%)
τt
(h)
τc
(h)
τp
(h)
10.05000.53127.87.66−8.14−9.313.53.53.5
20.04800.45428.27.69−7.65−2.853.06.530.5
30.03600.35428.35.1811.312.942.59.057.5
40.04550.42727.65.9518.621.942.511.583.5
50.04900.32827.75.2315.184.393.014.5134
60.04800.30627.45.41−11.5015.552.517.0161
70.05350.31228.25.3314.148.652.519.5187
80.05350.35028.95.190.71711.023.022.5214
90.06300.38325.65.4115.295.624.026.51658
100.07000.33926.34.391.81−1614.531.01687
110.09150.47921.46.808.452.355.036.02892
120.10400.44121,57.16−12.743.834.040.02920
130.10200.30124.94.3918.004.795.045.05517
140.11500.26125.43.79−0.962.455.050.05645
150.12300.95025.84.65−13.424.685.055.05575
160.13100.29924.94.88−1.392.925.060.06252
170.13200.34923.96.26−3.214.865.064.06282
180.12400.36123.26.25−3.497.195.569.56584
190.12600.30621.36.97−3.317.965.074.56603
200.13500.32622.95.21−4.418.726.080.56663
210.14600.31223.15.010.489.665.586.06903
220.16700.49820.95.311.759.927.593.57014
230.13900.49820.26.993.049.717.0100.57654
240.15000.48920,06.30−6.2011.247.5108.07684
250.16700.48720.96.0516.7512.898.5116.57719
260.16500.57720,37.265.1611.088.5125.07789
270.20800.49520.85.1917.878.959.5134.57820
280.12500.32628.15.9515.1812.275.5140.07860
290.15500.33728.46.2115.2813.274.0144.07884
300.10900.37228.86.86−19.9817.313.5147.57908
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ilie, M.C.; Chiş, T.V.; Maior, I.; Răducanu, C.E.; Deleanu, I.M.; Dobre, T.; Pârvulescu, O.C. Experimental Investigation and Modeling: Considerations of Simultaneous Surface Steel Droplets’ Evaporation and Corrosion. Metals 2023, 13, 1733. https://doi.org/10.3390/met13101733

AMA Style

Ilie MC, Chiş TV, Maior I, Răducanu CE, Deleanu IM, Dobre T, Pârvulescu OC. Experimental Investigation and Modeling: Considerations of Simultaneous Surface Steel Droplets’ Evaporation and Corrosion. Metals. 2023; 13(10):1733. https://doi.org/10.3390/met13101733

Chicago/Turabian Style

Ilie, Marius Ciprian, Timur Vasile Chiş, Ioana Maior, Cristian Eugen Răducanu, Iuliana Mihaela Deleanu, Tănase Dobre, and Oana Cristina Pârvulescu. 2023. "Experimental Investigation and Modeling: Considerations of Simultaneous Surface Steel Droplets’ Evaporation and Corrosion" Metals 13, no. 10: 1733. https://doi.org/10.3390/met13101733

APA Style

Ilie, M. C., Chiş, T. V., Maior, I., Răducanu, C. E., Deleanu, I. M., Dobre, T., & Pârvulescu, O. C. (2023). Experimental Investigation and Modeling: Considerations of Simultaneous Surface Steel Droplets’ Evaporation and Corrosion. Metals, 13(10), 1733. https://doi.org/10.3390/met13101733

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop