Multifactorial Diffusion of Chloride Ions in Concrete considering the Non-linear Voltage Distortion of Other Species

One of the most relevant deterioration mechanisms of reinforced concrete (RC) structures is the corrosion of steel bars by chloride ions. The service life of structures could be extended by monitoring this process. The lack of correlation data between chloride ion diffusion and the electrical field has meant that alternative methods have been developed in this paper. One such method considers diffusivity and the electrical field as non-constant, although they are usually assumed constant, that is, without there being any variability of the involved factors. Therefore, a multi-factorial diffusivity that accounts for the material characteristics and environmental variations is adopted. In addition, an innovative stochastic approach to evaluate the variation of the electrical field while accounting for the non-linear voltage distortion, which considers a multispecies model, has been proposed. More complete 1D/2D electro-mechanical models that are solved by advanced numerical and analytical solutions for electro-mechanical diffusion have been developed. Results show that spatiotemporal numerical models estimate the chloride concentration well, estimating more reliable scenarios to provide information on the safety situation of the structure and make it more sustainable by employing successful preventive treatments.


Introduction
Old, reinforced concrete (RC) constructions can easily be found worldwide (e.g., viaducts, bridges, stadiums, and buildings), which will probably have structural problems in the near future.One of the main mechanisms causing failure in RC structures is the corrosion of steel rebars by chloride ions [1].
Chloride ions reach the steel rebars by migrating throughout the concrete cementitious matrix from the external environment.The corrosion of the steel reduces the area of the rebars, which weakens the resistance and ductility of the structure [2].This can lead to structural instability and pose a risk to life [3].
Structures built in aggressive environments such as industrial sites or near the coast are the most exposed to this kind of corrosion phenomenon [4].Recent studies have shown that the affected structures can be from 1.5 [5] to 4.5 km from the sea [6].This high figure can be correlated with the wind action as analysed in Ou et al. [7], also reaching 6.5 km from the coastline, where a possible correlation between wind and chloride concentration is shown.
The approach to monitoring the chloride diffusion described in this paper combines some of the key parameters of the diffusion, i.e., the diffusion coefficient (or diffusivity) adjusted for an equivalent electric field.Fick's first and second laws describe the traditional diffusion process.Diffusivity has mainly been used as a constant value without considering the effect on its value of factors like water/cement (w/c) ratio, binding chloride, temperature, concrete aging, and concrete deformation [8,9].Only in a few studies is concrete deformation regarded as a factor for diffusivity [10,11].
In Yari et al. [12] and Wang et al. [13], the aggregate volume and size effects are also considered to estimate diffusivity for inhomogeneous concrete.Fu et al. [14] introduces a factor that considers the change in pore structure of concrete due to carbonation.
Recent works [2,15] have focused on some of those aspects of correct diffusivity estimation.In Carrara et al. [15], diffusivity has been separated into four different types (i.e., self-diffusivity, free diffusivity, diffusivity in saturated capillary pores, and solid diffusivity) to calculate which is better to use, whereas in Nguyen et al. [2], diffusivity is considered by using another factor, such as humidity.
All these studies on chloride diffusivity in concrete show how it is necessary to focus on diffusion uncertainties, e.g., material characteristics, models, environmental exposure, and electrical diffusion.Regarding the last parameter, the literature often does not consider electrical diffusion (or electromigration) [1,16], which means the diffusion process is purely mechanical [17].
Multispecies approaches based on coupling Nernst-Planck models have become popular in the last ten years because they account for all the ionic species that participate in transport, not only chloride ions [16,18].Nernst-Planck models include the electric potential gradient as a driven force for ionic transport, which can be used both in modelling chloride tests, such as ASTM C1202 [19], and in chloride self-diffusion in concrete.The electric potential gradient can be described by a linear-ohmic potential and a non-linear voltage distortion [20].
Electro-mechanical chloride ion diffusion can be a tool for structural health monitoring (SHM) systems.Despite the interest in this technique over the last twenty years, there is a gap between the theoretical concept and its application due to low competitiveness and the high cost of sensors [21,22].Recent studies refer to SHM systems as "futuristic ideas" for civil engineering [23] and an "emerging topic" since they make it possible to foresee potential damage and identify the real structural behaviour of structures [24].
In some papers, for instance, the electro-mechanical interactions have been studied by adopting carbon nanotubes to monitor a structure's response under loading up to damage through varying the electrical resistance [25].In Sanchez-Romate et al. [26], some relationships between resistance and deformation were established, showing a quasi-linear trend.
Therefore, this paper sets out a novel study on the electromechanical diffusion of chloride ions using non-constant multi-factorial diffusivity.The internal potential (or selfdiffusion non-linear voltage distortion) is also considered.
The new contributions, which should complete some aspects partially solved by the authors of this paper (see among others [1,[27][28][29]), concern: 1.A multi-factorial diffusivity to combine concrete deformation with the chloride ion diffusion.Given that in the literature [25,30] several studies show a correlation between concrete deformation and the resistivity, it could be possible to find a direct relationship between chloride ion diffusion with resistivity through a diffusivity factor that accounts for the concrete deformation.2. A non-constant equivalent electrical field and mechanical diffusivity are considered when estimating the fluctuation of the chloride concentration in time t and space x.This analysis is developed using advanced numerical solutions given by Hermite polynomials (HPs) [28,31] that account for the variable characteristics of the material and external environmental factors.3. The electrical field correlates with the external ohmic and self-diffusion non-linear voltage distortion.This correlation can be analysed using coupling Nernst-Planck equations.Here, a stochastic analysis was proposed to estimate the variation of the self-diffusion.
From these aspects, it is possible not only to estimate the chloride concentration entirely but also to try to find a correlation between the diffusivity/deformations/electric fields.A correct estimation would avoid useless repairs, treatments, and maintenance strategies with high costs and extended timelines.It is preferable to obtain accurate predictions of the structure response, pre-and post-chloride attacks, during design and construction rather than in-service.

Transport of Ions
The following 2D mechanical model describes the diffusion of chloride ions in concrete structures.The chloride concentration, C(x,t), in direction x (i.e., x 1 , x 2 ), and time t in saturated concrete is governed by partial differential equations (PDEs) called Fick's second law [32]: (1) where ∇C(x,t) is the total chloride concentration gradient, and D cl is the diffusivity.The 1D problem in, let us say, x 1 direction has an analytical solution under the following initial condition (i.c.) a Dirichlet boundary condition (b.c.), respectively: where C 0 and C s are the inner and surface chloride concentration, respectively [33,34].Therefore, the 1D analytical solution is: where erfc(•) is the complementary error function.The total chloride concentration, C, is divided into free chloride dissolved in a pore solution, C f , and chloride binding concentration, C b .As mentioned in [35], the "chloride binding process is very complicated and influenced by various factors, which is generally simplified by the binding adsorption isotherms".Here, the correlations between C f and C b at a given temperature are defined by the non-linear Langmuir isotherms [32,35,36], thus: where ω e is the evaporable water content [37] and α and β are material constants (here α 11.8, β 4.0 [38]).
C b is usually assumed to immediately reach equilibrium since the binding process is faster than the diffusion process, therefore it can be considered constant [39].Also, in [40] there is a linear relation between C and C f under specific conditions that indicates a small contribution of C b .

Inter-Combination of Electric/ Diffusion/Deformation Fields
In this section, three coupled fields are explained: (i) the electric diffusion with the elastic deformation of concrete; (ii) the mechanical diffusion with the elastic deformation of concrete; (iii) the electric diffusion with the mechanical diffusion.
If the first and third fields are known, the second field is often neglected, as already mentioned.The introduction of the three fields, correlated with each other, serves to justify the correlation of the current/diffusion/deformation proposed in this paper.

Coupled Electro-Elastic Fields
The electro-elastic model used to study the composite scheme has been retrieved from the literature [27].The general relationships which linearly couple the electrical and mechanical fields are expressed as: where σ kj and ε kj are the components of the elastic stress and strain tensors, respectively; D j are the components of the electric flux density vector; E i are the components of the electric field vector; c kjil , e kji and kj are the elastic coefficients, piezoelectric coefficients, and the dielectric permittivity coefficients.
The relation between the strain field ε kj and the displacement vectors u k,j , u j,k , and that between E k and the gradient of the electric potential, ∇φ k , are, respectively: The fundamental relationships are further subject to conditions of equilibrium and Gauss's law.Under the assumptions of vanishing body forces and vanishing volume charge densities, these governing equations are: Equations ( 6), (7), subject to assumptions given by Eq. ( 8), determine the steady-state behaviour of a linear electroelastic system.

Coupled Mechanical Diffusion-Elastic Fields
As already mentioned, the diffusivity D cl is the key parameter of the diffusion process.It can be considered constant [10], variable [15], or variable multi-factorial [41].
A complete multi-factorial relation was recently proposed by the authors of this paper in [28,42], and it could be described as: where the factors f 0 (w/c), f 1 (C b ), f 2 (T), f 3 (t), and f 4 (h) are correlated with the water/cement (w/c) ratio, C b concentration, temperature T, age t, humidity h, respectively.These factors are not explained here for the sake of brevity, but they can be found in the literature [10,11,41].Equation ( 9) is an empirical equation that accounts for some external actions (e.g., temperature and humidity) and internal characteristics (e.g., w/c ratio, concrete age) of the used material.In [9], several empirical equations have been listed indicating a possible combination of these factors to estimate a reliable D cl .To the best of the authors' knowledge, Eq. ( 9) could represent a good model for estimating D cl .Diffusivity, D cl 0.34 cm2 /year (Eq.9) a For w/c 0.50 [11,63] b For ω e 0.14 and C f C 0 + C s 1.65 kg/m 3 [10,38,63] c For R 8.314 J/mol/K, activation energy E a 41.80 kJ/mol, reference temperature T 0 296.0K and T 293.0 K (i.e., 20.0 ºC) [10,53] d For reference ager t ref 28.0 days (i.e., 0.08 years), t 30.0 years and speed index m 0.20 [10,64] e For reference humidity h ref 0.75, spread factor n 4.0 and h 0.60 [41,53] f For ε kj 1.50‰ (< 5.0‰ [26]), and B -3.55 × 10 -3 [11,21] Liu et al. ( 2021) [43] presents a possible relationship between electricity and the concrete age.Indeed, it could be possible to correlate the f 3 factor with E k .Similarly, here we use the f 5 factor.Particular attention is needed for this since it should allow the correlation of concrete's elastic deformation ε kj (x) with the chloride diffusion. 1This factor, shown in [10,11,44] is defined as: where B is a coefficient that depends on the material characteristics (Table 1), and experimental laboratory tests have estimated it.In [11], this coefficient has been related to the concrete material regardless of the damage levels.It has been fitted by correlations between strain values and D cl for sound concrete under loadings.It represents the situation where no damage occurs in the elastic state for compressive stresses.
The elastic state implies that the pore structure of the concrete changes with the increasing ε kj (x), and thus D cl values.
It is important to note that Eq. ( 10) refers only to a low linear compressive deformation.This implies that: (i) the effect of the breakage of electrical pathway due to matrix fracture and damage is not prevalent (e.g., ε kj (x) < 5.0‰) [26,45]; (ii) a low ε kj (x) guarantees a linear state [46] during the diffusion, thus also maintaining the validity of Eq. ( 7).

Coupled Electro-Mechanical Diffusion Fields
The electro-mechanical model is based on the Nernst-Planck model, which simultaneously describes several ions' mechanical and electric flux, in which chloride is usually of the most interest [1,16,39,[47][48][49].
By extending Eq. ( 1) and correlating D cl with the absolute value of |E k (x)| as indicated in [39], the Nernst-Planck equation is: where z is the valency of the ion, F is the Faraday constant, R is a universal gas constant, and T is the temperature.The a-parameter, considering Eq. ( 7), is expressed as a [16], which is a function of x-depth due to E k (x).It can be expressed as a constant or variable (e.g., a × sin(x)) parameter.
The analytical solution of Eq. ( 11) in x 1 direction, under the same hypotheses used to obtain Eq. ( 3) plus the aparameter considered constant, is [16,39,50]: Equation ( 11) is usually developed only considering the a-parameter as a constant value, whereas C 0 0 and C s 0 in Eq. (12).In this respect, Eq. ( 12) is shown in full.

1D/2D Constant Diffusion
The exact solution of Eq. ( 1) for D cl 1 by using the Kampé de Fériet polynomials (KDFPs) of non-negative degree n, in direction x 1 and time t (x 1 , t ∈ R), is defined by the series [51,52]: By substituting Eq. ( 13) in Eq. ( 1), the following relation is obtained: i.c.: For 2D constant diffusion, Eq. ( 14) is also valid when t is replaced with x 2 .Therefore the boundary conditions are H n (x 1 , 0) x 1 n and H n (0, x 2 ) 0.

1D/2D Non-constant Diffusion
For D cl 1, Eq. ( 13) is not the exact numerical solution of Eq. ( 1), but it can be used to obtain an approximate solution.
For the 1D model, the ordinary Hermite polynomials (HPs) in x 1 and t are specified by the series [12,52]: To solve Eq. ( 1) and Eq. ( 11) by considering a nonconstant D cl and a-parameter, the following conditions are defined: where x 1,m , and t m are two variables related to position and time, respectively.They must be pre-defined not to generate the arbitrary constants that represent the degrees of freedom associated with i.c. and b.c.; q indicates a multiplier parameter that quantifies the convergence (here, it is assumed q ≈ 20.0 as already validated by the authors of this paper in [28,53]).The 2D model is obtained by adding the following series to Eq. ( 15) [12,52]: To solve Eq. ( 1) and Eq. ( 11) by considering 2D nonconstant D cl and a-parameter, with a fixed t, the conditions are: where x 2,m is a pre-defined variable related to the position like x 1,m .

Materials for Mechanical Diffusion
The diffusion process is divided into mechanical and electrical diffusion.Table 1 shows the adopted parameters for mechanical diffusion for ordinary Portland concrete.References justify some used parameters, others are directly calculated; in this respect, they should assume realistic and reliable values.
The external conditions regard a moderate level of aggressiveness (i.e., structures located between 0.10 and 2.84 km from the coast without direct contact with seawater, providing C s 1.15 kg/m 3 [54]), with temperature T 20 ºC and humidity 60.0% (i.e., h 0.60).The main internal conditions concern a w/c 0.50 and C 0 0.50 kg/m 3 [11,55].
As mentioned in Sect.2.2.2, the f 5 factor needs particular attention since it represents the "bridge" between the elastic field and the electro-mechanical diffusion.In this study, the f 5 factor is estimated as 4.36 using the empirical Eq. ( 10).This factor is difficult to calculate since it correlates the elastic deformation ε kj (x), due to external actions on an RC element, with the diffusion of chloride ions.In the literature, it can take a value of 1.89 in [56], 3.23 in [17] up to 8.95 in [10,11], indicating a significant impact on the final D cl , as also shown in Fig. 1.
Figure 1 shows the weight of each multiplicator factor to obtain D cl in this study.The weights are calculated as a percentage (%) with respect to the reference factor f 0 .
The weights shown in Fig. 1 help quantify the general role of each factor.Here, f 4 has the lowest weight (f 4 (h) 3.04%) for estimating D cl .Obviously, these weights can vary as a function of the adopted parameters, but generally the predominant factors are f 0 and f 5 [10,42].

Materials for Electrical Diffusion
The coupled electro-mechanical diffusion fields in Eqs.(11)(12) need to be solved for the electrical field, for which an equivalent average value was proposed in the present paper.
In concrete, the electrical field E k suffers a further internal variation.The variation of E k in direction x 1 could be formed by an external ohmic potential, which corresponds to the external voltage applied to a sample in a migration test such as the ASTM C1202 [19] and NT492 [57], and a self-diffusion non-linear voltage (or potential) distortion when the pure diffusion occurs.
The method proposed in this paper to calculate an equivalent electric field, E k , to represent the phenomenon corresponding to the non-linear potential distortion of the ionic self-diffusion, was obtained from a multispecies approach, in which several ions are involved (i.e., Na + , K + , and OH − from the concrete pore solution, and Cl − penetrating from the outside of the sample).The voltage distortion in concrete subjected to self-diffusion is formed since when all the species start to move, they cannot move freely.This is because they are charged particles that interact with the other species, generating a non-linear voltage distortion.
The theory of the two potential contributions in the x 1 direction are described by [20]: where c i is the ionic concentration of species i, u i is the mobility of species i in the pore fluid, I is the external total current density.The other parameters have already been explained.The method used to obtain a characteristic value and determine which values this electrical field could vary; many simulations were carried out using the multispecies model reported by [58].In the model, it is assumed that each ion flux is dependent on the other ions; therefore, to maintain the principle of electroneutrality (i.e., ensuring no excess of charge [20]), a non-linear potential distortion is generated even if there is no external E k .Then, in a self-diffusion condition with no external voltage source, although the external  Electric field, E k 0.320 ± 0.174 V/m (Fig. 2) Valency of the ion, z 1.0 [47] a-parameter 0.13 ± 0.07 cm −1 (in Eq. ( 11)) Electric current, I 2.30 × 10 -6 A a a Calculated by I (A s × E k )/ρ ≤ 3.30 × 10 -4 A (risk-free value for humans [60]), with a square section area A s 3.60 × 10 -3 m 2 and a specific resistivity ρ 5.0 × 10 2 ohms × m voltage applied is zero, the internal mobility of the species in concrete generates a distortion of the electric potential.About 1100.0 simulations were carried out, randomly varying the diffusion coefficients of the ionic species [59], the chloride binding capacity, and the concentration of the pore solutions of the species involved.Although the potentials generated by the ions coupling would vary in both time and space, for simplicity, here, only the voltage distortion in the central part of a concrete sample (i.e., at 3.0 cm) for several years was considered.Figure 2 shows E k values by histograms and the fitted probability density function (PDF) by a normal distribution, in which the average value obtained was 0.32 V/m.Table 2 lists the adopted parameters for electrical diffusion.Some parameters are justified by references, equations, or specific studies, as shown in Fig. 2; therefore, they should take realistic and reliable values.It should be noted that although the coupled mechanical diffusion-elastic and electro-mechanical diffusion fields approaches are related only to chloride ions, the multispecies model (OH − , Na + , Fig. 3 The sinusoidal trend of a non-constant spatial a-parameter (a 0.13 cm −1 ); b non-constant temporal D cl (D cl 0.34 cm. 2 /year) K + , and Cl − ) was used here only to obtain a characteristic value for E k .
The key parameter in Table 2 is the electrical field E k since it plays a double role in this study: (i) by conducting a multispecies analysis (Eq. ( 19)) it is possible to achieve a probabilistic estimate of E k thus the a-parameter used to carry out Eqs. (11)(12) to evaluate the more complete spatial/temporal electro-mechanical diffusion; (ii) from E k it is possible to obtain an equivalent electrical current, I, for a determined concrete geometry; thus, a possible correlation could be obtained.
Figure 3 shows the trend of the two key parameters that govern the electro-mechanical diffusion process.In Fig. 3a) it is possible to see how D cl affects the electrical a-parameter by reducing its amplitude.
A sinusoidal trend could simulate the trend of the T and h for the mechanical diffusivity, as shown in [41,42], and for the electrical diffusivity, an alternating current (AC) [60].Both are usually represented by a sinusoidal function, which is why this function was chosen in this study, although other function types could be adopted [53].

General Methodology
The general methodology of this study is divided into four main steps, as shown in Fig. 4.
Step 0. The possibility evaluated of finding a correlation between the chloride ion diffusion (mechanical and electrical) and elastic deformation of the material.This step examines research studied in previous works [10,28,58,61,62].

Fig. 4 General methodology of this study
Step 1.In this step, the mathematical models are described considering the mechanical and electrical diffusion.Numerical, modified analytical, and multi-factorial diffusivity solutions are developed.The self-diffusion approach is also proposed.
Step 2. The three coupled fields are defined: (i) electroelastic, which correlates the electrical field with the material deformations; (ii) diffusion-elastic, which correlates the mechanical diffusion of chloride ions with the material deformations; (iii) electro-diffusion, which correlates the electrical diffusion with the mechanical diffusion.Thus, the correlation of the current/diffusion/deformation is proposed.
Step 3. Finally, the chloride ions concentrations C in concrete in time t and space x by 1D/2D model are estimated.Several factors (mechanicals, electricals, environmental, constitutive) are included.The results thus provide several scenarios to anticipate the chloride ions attack within RC structures.

Constant Mechanical and Electro-Mechanical Diffusion
This section shows the mechanical and electro-mechanical chloride ion diffusion, implemented by Mathematica software [51].Both key parameters (i.e., a-parameter and D cl , shown in Tables 1-2) for estimating the diffusion process are considered constant.
It is therefore possible to develop analytical solutions and quantify the merit of the results since they represent the exact solutions.Analytical solutions should provide reliable results given that D cl and a-parameter have been estimated using more complete models.
In Fig. 5a, c, for the curve t 0, C(x 1 , t) C 0 ; for this reason, t > 0 is adopted to show a curve.In Fig. 5b, d), the horizontal line corresponds to x 1 0 when C(0, t) C s following Eq.( 2).
Two main aspects can be noted.The first concerns the diffusion velocity of the chloride ions by including an electric field.The diffusivity D cl is amplified by the a-parameter that allows the chloride concentration C to reach the depth x 1 6.0 cm faster (see black lines in Fig. 5).However, as shown in Fig. 3a), this amplification is not very large (i.e., about 0.38-0.300.08 cm/year).
The other aspect is the sensitivity of the a-parameter variation, which is more evident for a higher a-parameter (see Fig. 5c-d)).This indicates that it could be challenging to combine the mechanical and the electro-mechanical diffusion for E k > > 0 since the two curves could be very distant from each other.Thus, the two diffusions could not be compatible with each other for monitoring RC structures.Also, it is noted that the electric diffusion affects the total diffusion from t > 2.5 years.
Results of the electro-mechanical model, shown in Fig. 5, controlled by the a-parameter (a function of the electrical field), include, in a simplified macroscopic way, the complex physicochemical phenomena that could be modelled through multi-species approaches in which the ion-ion effects are included, coupling all the species involved, to the charge neutrality.From Fig. 5 it can also be seen that the mechanical model tends to underestimate the penetration of chlorides with respect to what was obtained using the electro-mechanical models.As was expected, this condition is exacerbated as the electric field increases.
The numerical solutions for the mechanical diffusion (Fig. 6c-e) are very similar to the analytical solution with an error of ~2.0%, whereas for the electro-mechanical solutions (Fig. 6d-f), the approximation is a bit worse (~5.0%).This indicates that Eq. ( 3), as expected, approximates the results quite well, whereas Eq. ( 12) provides worse results, probably due to double erfc(•).However, a difference between 2.0 and 5.0% can still be considered acceptable.Therefore, this also indicates that the q-parameter (Eq.( 16)) adopted for the numerical solutions provides good results.
Also, we can see that numerical solutions have some advantages; for instance, it is possible to plot a curve for x 1 0 cm (non-constant outputs in Fig. 6e-f) for non-constant a-parameter and D cl (shown in the following sections).
Figure 7 shows the 2D chloride electro-mechanical diffusion in x 1 and x 2 direction for a 0.13 ± 0.07 cm −1 and for t {5.0, 25.0} years.The numerical solutions are obtained using Eq. ( 1) and Eq. ( 11) with the conditions shown in Eq. (18).
Although the variation of the a-parameter is not very high, it is possible to see its effect on the chloride concentration in the 2D model.In particular, at t 25.0 years and x 1 x 2 6.0 cm (i.e., where steel bars could be placed), the maximum C with a 0.20 cm −1 is greater than the maximum C for a 0.06 cm −1 of 1.06 (see Fig. 7c-d).
When the chloride ions come from two directions, x 1 and x 2 , the chloride concentration C is expected to be greater than for unidimensional diffusion.In fact, the ratio between the chloride concentration C is always C 2D /C 1D ≥ 1.0.In particular, at 6.0 cm, this ratio is 1.30 for 5.0 years and 1.24 for 25.0 years, in accordance with [53].This ratio decrease is probably due to a saturation phenomenon that reduces the concrete pores.

Spatial Non-Constant Electro-Mechanical Diffusion
The electro-mechanical diffusion for a non-constant aparameter is shown in Fig. 8.These numerical results, Fig. 5 1D mechanical and electro-mechanical chloride diffusion using analytical solutions for D cl 0.34 cm 2 /year and a 0.13 ± 0.07 cm −1 .Curves are plotted for t {0.01, 5.0, 10.0, 15.0, 20.0, 25.0} year and x 1 {0, 1.0, 2.0, 3.0, 4.0, 5.0, 6.0} cm obtained by developing Eq. ( 15) with the condition in Eq. ( 16), account for the variability in the space of the electrical field E k (x 1 ).The curves are plotted for different t.Given that variability is spatial, the C curves plotted in time t could not provide significant alterations in this case.
The curves in Fig. 8a) give higher concentration values than those in Fig. 8b).Different from the results with constant a-parameters (Fig. 5a)), here, with a low a-parameter, the chloride concentration increases a little.This estimated rise, at 6.0 cm for 25.0 years, is ~1.17; therefore, it is possible to affirm that by considering a-parameter as constant (as usually assumed), the results could provide slightly underestimated values.This effect is due to the fluctuating trend shown in Fig. 3a).
These non-expected results could be correlated to the fact that, for a non-constant D cl , the concentration of the chloride ions in a specific section x 1 is not cumulative; in fact, in general, the physical meaning of the chloride diffusion is that the concentration C increases in t with a constant D cl .When D cl is non-constant, this physical meaning changes; therefore, the mathematical influence of the a(x 1 ) with a negative sign in Eq. ( 11) alters the results.Here, the variability in the space of the function a × Sin(x 1 ) provides these changes.

Temporal Non-constant Electro-Mechanical Diffusion
The influence of time variability is studied in this section.Figure 9 shows the 1D mechanical and electro-mechanical diffusion for a non-constant diffusivity in t.The curves are obtained by developing Eq. ( 15) with the conditions indicated in Eq. (16).
The difference between the two models is almost imperceptible; it is very challenging to quantify the chloride ion diffusion by comparing mechanical and electrical models.At t 25.0 years and x 1 6.0 cm, the electro-mechanical Fig. 6 1D electro-mechanical and mechanical diffusion using numerical solutions (D cl 0.34 cm 2 /year and a 0.13 cm.−1 ) diffusion provides larger values of 1.04 than mechanical diffusion.
Physically, the difference could be small because electrical diffusion is much faster than mechanical diffusion.Moreover, the temporal interval of the diffusivity is different, with the former being studied in seconds and the latter in years.
Evident differences between the chloride concentration in time variability with respect to a constant analysis relate to the oscillations of the curves.In fact, the curves in Fig. 9 are directly affected by the sinusoidal trends shown in Fig. 3b), which could provide values lower than curves for a constant analysis (Fig. 6e-f).In particular, at the point (t 5.0 years, x 1 3.0 cm), the difference is 1.05, i.e., the temporal nonconstant curve provides a smaller value.
Considering that the model is dynamic, it is possible that for the same depth x 1 , the chloride concentration C decreases or increases in t differently for the constant D cl , which increases the concentration C monotonically, as shown in [28].When a sinusoidal diffusivity is composed of more points placed above the constant D cl (horizontal line in Fig. 3b), the chloride concentration by numerical solution is greater than the chloride concentration by exact solution and vice-versa.This can be also noted in Fig. 11 (shown later), where the difference between C by constant and temporal analysis at 10.0 years (i.e., 10.0 ≈ 3π) is greater than at 25.0 years (i.e., 25.0 ≈ 8π).

Spatial/Temporal Non-constant Electro-Mechanical Diffusion
This section shows a complete model by considering all dependencies, i.e., the variation of the diffusivity in t, x 1, Fig. 7 2D electro-mechanical diffusion in two directions, x 1 and x 2 , using numerical solutions for a .13± 0.07 cm −1 at 5.0 and 25.0 years Fig. 8 1D electro-mechanical diffusion by using numerical solutions for D cl 0.34 cm 2 /year, a 0.06 Sin(x 1 ) cm −1 , and a 0.20 Sin(x 1 ) cm.In the literature, the diffusivity is divided into intrinsic and apparent diffusivity, which are correlated by porosity and chloride binding [15,49].In this study, these aspects are covered in Eq. ( 9) b ratio is calculated as the ratio between C values at x 1 1.0 cm and those at x 1 0 cm and x 2 .Figure 10 shows the electro-mechanical diffusion in the 1D/2D models.
Results in Fig. 10 show that the chloride concentration further increases by considering spatial and temporal variations.In the 2D model, the maximum values are higher than (~1.22) the values shown in Fig. 7c-d).
Figure 11 compares the results of chloride concentration C in x 1 of the four 1D models.It is possible to confirm that the chloride concentration with the constant parameters (dashed curve) is lower than that with spatial-temporal parameters (solid curve).This should indicate that the analytical solutions will not estimate the scenarios correctly.
Figure 12 compares the same four models as Fig. 11 but with some changes.The difference is that the sinusoidal functions shown in Fig. 3 have been translated downwards by 0.34.Thus, they have upper peaks coinciding with D cl .This makes quantifying the fluctuation trends with positive and negative amplitudes possible.This could simulate, for instance, (i) internal convective flux of chloride ions with retrocession of the diffusion, (ii) possible interactions between different species (described in Sect.4.1.2),and (iii) binding effects that decelerate the diffusion.
In this case, the chloride concentration with the constant parameters is higher than that with spatial-temporal parameters.In particular, for a lengthy period (e.g., t 25.0 years), the chloride concentration for constant parameters (i.e., using analytical solutions) is somewhat overestimated.This indicates that analytical solutions should not be used to estimate relatively long-term scenarios.

Current/Diffusion/Deformation Correlations
Table 3 compares the results of this analysis with some studies published in the literature.The comparison focuses on the a-parameter, D cl , and chloride concentrations (defined as a ratio between the concentration at x 1 1.0 cm and x 1 0 cm, i.e., ratio < 1.0).This ratio represents the available range to compare all studies since it could be difficult to standardize all results if a broad range is considered.
Although the results shown in Table 3 are provided by analyses that differ slightly from each other, it is still possible to discuss some aspects.This is because the key parameters of the electro-mechanical diffusion are E k , thus a-parameter, and D cl .Some results underestimate the chloride concentration (i.e., ratio < 0.5) concerning the analytical solutions that would provide, at x 1 1.0 cm, a value of ratio ≈ 0.80.This is probably due to a poor calibration of the a-parameter in the published studies.For this, a stochastic approach proposed in this paper should provide better values for the a-parameter to reliably monitor chloride diffusion in RC structures.
Finally, in the literature [26,45], there are some relations between ε kj and the variation of I or resistance where a quasilinear relation is shown.Similarly, some correlations between the a-parameter vs. f 5 and I vs. x 1 are estimated in Fig. 13.
The correlation in Fig. 13a) indicates a possible equilibrium between the concrete deformation ε kj , the factor f 5 , and the electric field quantified by the a-parameter.These curves could indicate how much it is necessary to vary ε kj thus f 5 with respect to the a-parameter to individuate diffusivity D cl , therefore a same chloride concentration C in concrete.Consequently, these bi-linear curves show a pair of values that rapidly estimate D cl using Eq. ( 9) and concentration C using Eq.(12). Figure 13b) shows a possible correlation between I and x 1 .Under some conditions of the concrete these curves should indicate A s and ρ (see Table 2) adopted in this paper, what is their current value at a certain depth x 1 .The I values can be used to estimate the a-parameter (i.e., I → E k → a).The upper curve is defined as I + σ and the lower curve as I-σ, with a standard deviation σ 0.174.Figure 13 tries to correlate the current/diffusion/deformation, which is one goal of this paper.In this way, some values should be obtained to correlate electrical parameters (I and a) with mechanical/geometrical parameters (f 5 and x 1 ).In [16,61], the conductivity, which is directly proportional to the electrical current I [24], increases as C increases.This could indicate that the conductivity and I increase as the electrical diffusivity increases.

Conclusions
The diffusion of electric-mechanical chloride ions in RC elements by non-constant multi-factorial diffusivity has been developed using advanced numerical solutions.The external ohmic and self-diffusion non-linear voltage distortion of multispecies have also been considered.2)) The main conclusions are: 1.A multi-factorial diffusivity D cl (Eq.( 9)) has been defined to account for several factors that include the variable characteristics and external environmental context of the material (i.e., w/c ratio, binding chloride, temperature, concrete age, humidity, elastic deformation).Some of these factors are usually neglected when estimating the chloride concentration in RC elements; however, it is shown that they have an important role with a weight ranging from 3.0 to 99.0%.The elastic deformation factor, f 5 , represents the "bridge" factor to justify the proposed current/diffusion/deformation model (Sect.2.2). 2. The electrical field variation was estimated as the sum of the external ohmic and self-diffusion non-linear voltage distortion (Eq.( 19)).The former corresponds to the external voltage applied to the concrete.In contrast, the latter is formed when the charged species start to move.The multispecies approach has thus been carried out under ~1100.0simulations in an alternative probabilistic way.The a-parameter of the Nernst-Planck equation is estimated as 0.320 ± 0.174 V/m, which plays an important part in electrical-mechanical diffusion.3. The analytical equation for electro-mechanical chloride ion diffusion (Eq.( 12)) has been adapted to quantify the contribution of C 0 present in the concrete.Numerical solutions provided by HPs have been implemented to develop innovative results in 1D/2D directions.It has been shown that by considering a non-constant D cl and a-parameter (i.e., spatial-temporal analysis), the chloride concentration in an RC element increases by up to 1.06 kg/m 3 (i.e., close to a moderate level of aggressiveness [54]).In general, it is difficult to predict the real contribution of the electro-mechanical mechanisms.Some correlations were also estimated, e.g., f 5 vs. aparameter and I vs. x 1 .

Fig. 1
Fig. 1 Weight of each factor with respect to f0 in the mechanical diffusivity

Fig. 2
Fig.2The electrical field, E k , developed for self-diffusion conditions (V/m): stochastic data (histograms) and PDF normal distribution (red line)

Fig. 11 Fig. 12
Fig. 11 Comparison between four models: a for t 10.0 years and b for t 25.0 years

Fig. 13
Fig.13 Equilibrium between a) f 5 vs. a-parameter and b) I vs. x 1 (the upper curve is defined as I + σ and the lower curve as I-σ, with σ 0.174 (Table2))

Table 1
Adopted parameters for the mechanical diffusion

Table 2
Adopted parameters for electrical diffusion

Table 3
Comparison between some published studies and this analysis *Estimated value a