Simulation Approach for Cathodic Protection Prediction of Aluminum Fin-Tube Heat Exchanger Using Boundary Element Method

Abstract

The multi-galvanic effect of an Al fin-tube heat exchanger was evaluated using polarization tests, numerical simulation, and the seawater acetic acid test (SWAAT). Determination of the polarization state using polarization curves was well correlated with numerical simulations using a high-conductivity electrolyte. However, the polarization results did not match those of the low-conductivity electrolyte due to the lower galvanic effect. Although the polarization state is changed by electrolyte conductivity, the total net current of the tube is decreased in the case of the anodic joint. From SWAAT results, the leakage time of Al fin-tube heat exchanger assembled by anodic joint was longer than the case with cathodic joint.

1. Introduction

Aluminum (Al) is a light metal that has interesting properties for heat exchanger applications (e.g., its low density, high thermal conductivity [1], good corrosion resistance [2], and good mechanical properties [3]) in the HVACR (heating, ventilation, air conditioning, and refrigeration) industry for use in cooling systems and air-ventilated units [4,5,6]. In the mid-1990s, the mechanical assembly of automotive heat exchangers started using brazed Al alloys, and this trend is currently applied in the heat exchangers of air conditioners [7]. In air conditioners, AA 1xxx and AA 3xxx series alloys are used because automotive heat exchangers do not require high mechanical properties and these series have higher thermal transfer efficiency and are economic advantages [8].

To increase the wettability of the melted filler metal at the surface, fluxing is used to remove the natural oxide layer covering the Al surface [9]. Due to the high temperature and cladding with the flux material [10,11], the geometry and microstructure of the joint region, tube, and fin are modified [12,13,14,15]. The microstructural change can influence the corrosion behavior. Thus, the microstructural effects on the corrosion properties after brazing should also be considered. In addition, the Al series used in the tube, fin, and joint (filler metal) are different according to their applications. In this study, AA 1100, AA 3003, and AA 4343 or modified AA 4343 are used in the tube, fin, and filler metal, respectively. For assembly of the fin and tube, each part underwent brazing and the fin-tube was jointed as shown in Figure 1. Although Al is used in all parts of the fin-tube heat exchanger, the different alloying elements in the Al series influence the corrosion properties, resulting in galvanic corrosion [16,17].

An important consideration related to the corrosion protection of Al fin-tube heat exchangers is the galvanic coupling between the fin and tube materials. Generally, the material used in the fin has a more negative potential than the tube materials in order to provide cathodic protection for the tube. However, in the case of the Al fin-tube heat exchangers that are assembled by the brazing process, the brazing joint influences the galvanic coupling; this is called multi-material galvanic corrosion. The basic principles of galvanic corrosion are well established and commonly accepted, but galvanic corrosion in realistic situations, such as with multi-material coupling and complex geometries, is hard to predict. Fortunately, the advance of computational simulations has made it possible to model many complex corrosion situations; thus, computational simulation can be directly used to solve engineering design problems [18]. Computational simulation approaches are one of the most effective methods for corrosion design of a product.

Generally, the outside part of the Al fin-tube heat exchanger is the main corroded region because the acid rain and outside pollutants such as sulphur oxides (SO_x) and nitrogen oxides (NO_x) cause a severe corrosion problem. Thus, a test for outside corrosion is mainly considered for Al fin-tube heat exchanger. One of the corrosion reliability tests of the Al fin-tube heat exchanger is the sea water acetic acid test (SWAAT), which is a cyclic spray test (as opposed to a full immersion test). Therefore, in SWAAT, wet and high humidity environments are produced in turn. This indicates that the corrosion condition of the Al fin-tube heat exchanger is alternatively exposed to the thin-layer electrolyte and full immersion. The conductivities of full immersion electrolytes and thin-layer electrolytes are significantly different, which affects the potential and current distribution of the galvanic coupled metals [19]. Thus, the solution conductivity needs to be taken into account when designing a system that protects against galvanic corrosion. In this study, we evaluate the corrosion of Al fin-tube heat exchangers with different brazing joints in high- and low-conductivity electrolyte conditions based on computational simulations. Also, to compare the result of computational results, the actual SWAAT was conducted for Al fin-tube heat exchanger applied different joint materials. Although it would be not a perfect validation of simulation results, it can be supporting the simulation results for application of actual corrosion reliability test.

2. Boundary Element Method (BEM)

A model of an Al fin-tube heat exchanger is shown in Figure 2; this has electrolyte domains (Ω_bulk and Ω_thin-layer) and is surrounded by the surface of the electrolyte (Γ_n), the surface of the tube part (Γ_a), the surface of the fin part (Γ_b), and the surface of the joint part (Γ_c). The electrolyte conductivity (σ) is uniform in the whole domain and there is no current loss. The potential field in the electrolyte domain (Ω) can be modeled by Laplace’s equation [20]:

$${\nabla }^2\Phi =0$$

(1)

Here, Φ is the electrical potential, which is the potential relative to a reference electrode, such as a saturated calomel electrode (SCE).

Laplace’s equation is calculated using the following boundary conditions:

i = i₀, on Γ_n

(2)

i_a = f_a(Φ_a), on Γ_a

(3)

i_b = f_b(Φ_b), on Γ_b

(4)

i_c = f_c(Φ_c), on Γ_c

(5)

Here, Γ is the entire surface of the electrolyte domain, which includes Γ_n (electrolyte surface), Γ_a (tube surface), Γ_b (fin surface), and Γ_c (joint surface). f_a(Φ_a), f_b(Φ_b), and f_c(Φ_c) are the non-linear functions on the surfaces of the tube, fin, and joint areas, respectively, which represent the experimentally achieved polarization curves. Thus, the boundary element method (BEM) can be used to calculate Laplace’s Equation (1) when the tube, fin, and joint areas are prescribed on the Al surface and their polarization curves are known [21]. Based on this, Φ and i on the whole surface can be determined.

An Al fin-tube heat exchanger model was created using the Rhinoceros 4.0 (Robert McNeel & Associates, Seattle, USA)3D drawing software based on the real shapes and dimensions. The inner parts of the heat exchanger were not modeled to avoid computational errors because outside corrosion is the main consideration in this study. The 3D model of the Al fin-tube heat exchanger, with detailed dimensions, is shown in Figure 3. The joint part is modeled by rectangular shape between the tube and fin to simplify the heat exchanger model. Areas of the tube, fin, and joint are 870, 1115, and 40 mm², respectively. After modeling, the heat exchanger, the 3D model was imported into the program BEASY version 10.0r14, which is BEM-based software(BEASY, Southampton, England). Setting the boundary conditions is an essential step for corrosion simulation. Different electrolyte) conductivity values were applied for comparison of the bulk (0.4 S/m) and thin-layer (0.00004 S/m) electrolytes; this was done because the electrolyte conductivity is generally decreased in thin-layer electrolyte conditions. Although the oxygen concentration and ion transfer parameter can also influence the corrosion rate, we consider the electrolyte conductivity difference because the galvanic effect between the tube, fin, and joint is the main focus in this study. Also, the localized corrosion caused by galvanic corrosion is the important factor for corrosion reliability. However, the corrosion design for a large structure by simulation is the main goal of this study. Thus, not the localized corrosion which is focused on the micro-scale but the uniform corrosion in bulk structure was focused in this study.

3. Experimental Procedures

3.1. Materials and Solution

As mentioned above, AA 1100, AA 3003, and AA 4343 or modified AA 4343 (which was modified to include more Zn to decrease the corrosion potential), are used in the tube, fin, and filler metals, respectively. The chemical composition of each part is listed in

Table 1. Chemical composition of the parts (Al 1100, Al 3003, Al 4343 and modified Al 4343) of Al fin-tube heat exchanger.

Parts	Chemical Composition (wt%)
	Cu	Fe	Si	Zn	Mn
Tube (AA 1100)	0.001	0.200	0.010	-	-
Fin (AA 3003)	0.002	0.210	0.220	0.150	0.640
Cathodic joint (AA 4343)	0.110	0.120	5.380	0.080	0.010
Anodic joint (modified AA 4343)	0.120	0.120	7.730	0.480	0.010

. Cathodic and anodic joints indicate that the higher and lower corrosion potential than the tube materials. Fin material always has lower potential than tubes and joints for the corrosion protection of the tube. The test solution consisted of a seawater acidified solution which is a synthetic sea salt solution made with the addition of glacial acetic acid (pH 3.0), as described by ASTM G85. Test temperature of all the tests is 49 °C.

3.2. Potentiodynamic Tests

Polarization data for the different parts of the Al fin-tube heat exchanger (i.e., AA 1100, AA 3003, and AA 4343 or modified AA 4343) were needed to conduct simulations. Potentiodynamic polarization tests were conducted by a conventional three-electrode cell. A purified carbon rod was used as the counter electrode and a saturated calomel electrode (SCE) was used as the reference electrode. The tested specimen for carrying out the polarization test of each part was extracted from the Al fin-tube heat exchanger which was machined by a micro-cutting machine (The CUTLAM^®micro 1.1, LAMPLAN, Gaillard, France). The working electrode was abraded by a series of abrasive papers (from 220 to 600 grit), rinsed ultrasonically with ethanol, and dried with nitrogen (N₂) gas. The specimen was covered with silicone rubber, leaving an area of 25 mm² unmasked. The prepared specimen was then exposed to the test solution for 1 h. Potentiodynamic polarization tests were carried out using a Bio-Logic VSP-300 potentiostat. The potential range was from −0.3 V_SCE vs. the open-circuit potential (OCP) to −0.4 V_SCE at a scan rate of 0.166 mV/s. All polarization tests were repeated a minimum of three times to ensure accuracy.

3.3. Sea Water Acetic Acid Test (SWAAT)

As a corrosion reliability test of fin-tube heat exchanger, SWAAT has been used in the field of Al fin-tube heat exchanger. Thus, to compare the corrosion reliability of Al fin-tube heat exchanger applied anodic and cathodic joint and validate the corrosion resistance of suggested joint materials, SWAAT was conducted in this study. SWAAT was produced based on ASTM G85 which utilizes the following cycle: 30-minute spray followed by a 90-minute soak at above 98% relative humidity and 49 °C. Whole Al fin-tube heat exchangers were tested to evaluate the tube leakage time caused by corrosion degradation. To record the leakage time caused by corrosion degradation, the tube was filled with air at the inlet part and the pressure gage was installed at the outlet part. Until the pressure gage was set to 5 MPa, the air was filled. The leakage of the Al fin-tube exchanger was determined from the decrease of a pressure gage during the SWAAT.

4. Results and Discussion

4.1. Polarization Curves

Figure 4 shows the polarization curves of the tube (AA 1100), fin (AA 3003), cathodic joint (AA 4343), and anodic joint (modified AA 4343) in the seawater acetic acid solution. Parameters such as the anodic and cathodic Tafel slopes (β_a and β_c), corrosion current density (i_corr), and corrosion potential (E_corr) were determined from the polarization curves and are listed in

Table 2. Potentiodynamic polarization test results of different parts of the Al fin-tube heat exchanger (Al 1100, Al 3003, Al 4343, and modified Al 4343).

Parts	βa (V/decade)	βc (V/decade)	icorr (μA/mm2)	Ecorr (VSCE)
Tube (AA 1100)	0.021 $\pm$ 0.003	0.886 $\pm$ 0.01	0.109 $\pm$ 0.05	−0.716 $\pm$ 0.02
Fin (AA 3003)	0.010 $\pm$ 0.005	0.223 $\pm$ 0.01	2.351 $\pm$ 0.25	−0.792 $\pm$ 0.02
Cathodic joint (AA 4343)	0.109 $\pm$ 0.03	0.748 $\pm$ 0.09	0.129 $\pm$ 0.06	−0.657 $\pm$ 0.01
Anodic joint (modified AA 4343)	0.018 $\pm$ 0.002	0.859 $\pm$ 0.05	0.349 $\pm$ 0.09	−0.746 $\pm$ 0.03

. Although i_corr is different for each part of the Al fin-tube heat exchanger due to the alloying and microstructure effects [16,17,22,23,24], the actual current density in the Al-fin tube heat exchanger should vary due to the galvanic corrosion between the three metals. E_corr exhibits the following decreasing order: cathodic joint > tube > anodic joint > fin. This means that the fin acts as a sacrificial anode for both of the assembled cases; cathodic joint > tube > fin and tube > anodic joint > fin. Thus, the effect of different tri-metal galvanic couplings on the corrosion reliability should be considered for designing Al fin-tube heat exchangers.

Generally, in tri-metal galvanic couples that include three metals (A₁, A₂, and A₃) with different corrosion potentials in the order of E_A1 > E_A2 > E_A3, the polarization of the middle potential metal (E_A2) can be defined by the mixed potential theory [25,26]. If the tri-metal galvanic coupled metals (A₁, A₂, and A₃) comply with the Tafel system (a linear system), the polarization state of the middle potential metal can be determined by the Tafel superposition relation. Thus, the tri-metal galvanic couple can be investigated by considering the A₁-A₂ and A₂-A₃ couples. Then, the net current (I_net) on the A₂ metal can be determined as I_A1-A2 $-$ I_A2-A3, where I_A1-A2 and I_A2-A3 are currents flowing from A₁ to A₂ and from A₂ to A₃, respectively. This means that A₂ acts as an anode in the case of a positive I_net and as a cathode in the case of a negative I_net.

Based on the above determination made via the polarization state method, the polarization states of the tube in the two assembled cases (case 1: fin-tube-cathodic joint and case 2: fin-tube-anodic joint) were determined. Schematic polarization curves of the tube, fin, and joint are shown in Figure 5, which reflect the area of the 3D model used to interpret the polarization state; this is done because the ratio of the cathode-to-anode surface areas is an important factor in galvanic corrosion. In case 1, the intersection between the Tafel slopes of the cathodic joint and the tube cannot occur. Thus, the galvanic currents on the anode parts, which are the anodic currents on the tube and fin (I_T-CJ and I_F-T), are not increased after galvanic coupling according to the mixed-potential theory. Alternatively, the galvanic currents on the cathodic parts, which are the cathodic currents on the cathodic joint and tube (I_CJ-T and I_T-F), are increased. Consequently, I_net of the tube in case 1 is negative, acting as a cathode. In case 2, the Tafel slopes of the anodic joint and tube make an intersection. The galvanic currents on the anode parts, which are the anodic currents on the anodic joint and fin (I_AJ-T and I_F-T), are increased only in the anodic joint and tube couple. The galvanic currents on the cathode parts, which are the cathodic currents on the tube (I_T-AJ and I_T-F), are increased in case 2. The value of I_net for the tube in case 2 is also negative, acting as a cathode. In short, the polarization state of the tube in both cases was cathodic; however, this determination did not consider the distance between the anode and cathode or the electrolyte conductivity. Thus, a more detailed determination is needed to design Al fin-tube heat exchangers adequately.

4.2. Corrosion Simulation

Figure 6 and Figure 7 show the 3D models of the Al-fin tube heat exchanger with the cathodic and anodic joints, respectively, in high-conductivity electrolyte (0.4 S/m). The potential and current density ranges of the Al-fin tube heat exchanger assembled with a cathodic joint (Figure 6) were –775.2 mV_SCE to −774.9 mV_SCE and −2.21 $\times$ 10⁻⁷ A/mm² to 3.52 $\times$ 10⁻⁷ A/mm², respectively. While those of the Al-fin tube heat exchanger assembled with an anodic joint (Fig. 7) were −786.2 mV_SCE to –786.1 mV_SCE and −4.09 $\times$ 10⁻⁷ A/mm² to 1.37 $\times$ 10^-7 A/mm², respectively. In both cases, the anodic current was observed only on the fin part and the current density increased near the tube and joint parts. This is due to the increased galvanic effect between the metals which had different corrosion potentials. The potential distribution was not broad (almost single potentials of −775 mV_SCE and −786 mV_SCE) in the anodic potential range of the fin. This means that the fin part acted as a sacrificial anode to the tube and joint (anodic and cathodic joints) in the high-conductivity electrolyte, regardless of the polarization state of the joint. This result is similar to the determination of the polarization state using the polarization curves.

The total net currents of each part of cathodic or anodic joints in high-conductivity electrolyte were −8.87 $\times$ 10⁻⁶ A (cathodic joint), −5.93 $\times$ 10⁻⁵ A (tube), 6.81 $\times$ 10⁻⁵ A (fin) and −1.09 $\times$ 10⁻⁶ A (anodic joint), −5.94 $\times$ 10⁻⁵ A (tube), 6.05 $\times$ 10⁻⁵ A (fin). The total net current is related to the amount of current flowing from the anode or cathode. A decrease in the total net current simply indicates a decrease of the corrosion rate (according to Faraday’s law). Thus, the decrease of the total net current of the fin in the anodic joint implies increased corrosion life, which means that an anodic joint extends the corrosion life of the Al-fin tube heat exchanger compared to the cathodic joint in the high-conductivity electrolyte.

Figure 8 and Figure 9 show the 3D models of the Al-fin tube heat exchanger using cathodic and anodic joints in low-conductivity electrolyte (0.00004 S/m), respectively. Contrary to the results of the high-conductivity electrolyte, the potential distribution was clearly separated in each part of the Al fin-tube heat exchanger in both cases (anodic and cathodic joints). This means that the galvanic effect among the tube, fin, and joint was significantly decreased due to the low conductivity [19]. Also, these results indicate that the cathodic protection effect between exchanger components was very low, so that sufficient protection could not be obtained at the regions that are far away from the junction. The cathodic current was observed in the case of the cathodic joint, while anodic and cathodic currents were found in the case of the anodic joint. In the tube and fin, only the anodic current was shown, which implies that the current flow is too low to change the polarization state from anodic to cathodic. However, the anodic current was slightly decreased near the anodic joint (Figure 9c) because of a small current effect from the anodic joint. Consequently, the total net current was slightly lower in the case of the anodic joint (6.205 $\times$ 10⁻³ A) than in the case of the cathodic joint (6.329 $\times$ 10⁻³ A). However, according to the simulation results, the anodic joint was more favorable for extending the corrosion life of the Al fin-tube heat exchanger (especially in the tube part) in both conductivity environments. This can be verified by practical corrosion testing.

4.3. SWAAT Test

To compare the corrosion lifetimes of Al fin-tube heat exchangers assembled with cathodic or anodic joints, SWAAT was conducted until leakage of the tube occurred. The polarization states of the Al fin-tube heat exchangers assembled by cathodic or anodic joints and the corresponding leakage time are listed in

Table 3. Sea water acetic acid test (SWAAT) results and polarization states of Al fin-tube heat exchangers assembled with cathodic or anodic joints.

	Polarization State FROM SIMULATION		Leakage Time (day)
	High conductivity	Low conductivity
Cathodic joint	Cathode	Cathode	56 days
Tube	Cathode	Anode
Fin	Anode	Anode
Anodic joint	Cathode	Anode/Cathode	80 days
Tube	Cathode	Anode
Fin	Anode	Anode

. Although the polarization state was changed according to the electrolyte conductivity, the anodic joint was efficient in extending the corrosion life because the anodic joint either decreased the net current on the fin or it protected the tube by serving as a sacrificial anode. In short, based on the corrosion simulation results, it was revealed that the anodic joint can increase the corrosion life. In SWAAT, an increase in the corrosion life was observed on the Al fin-tube assembled with an anodic joint. The leakage time in the case of the anodic joint was increased by about 42% relative to the lifetime with the cathodic joint. This suggests that the numerical simulation results correlated well with the corrosion test of the product. Thus, corrosion simulations can be helpful for the design of Al fin-tube assembly in low and high-conductivity electrolytes.

5. Conclusions

In this study, the multi-galvanic effect in Al fin-tube heat exchangers assembled with cathodic or anodic joints was investigated using a polarization determination method and numerical simulation in the presence of electrolytes with high and low conductivities. SWAAT was also conducted to verify the simulation results. When determining the polarization state from the polarization curves, the tube and joint were found to be cathodes, while the fin acted as an anode, regardless of the corrosion potential of the joint. These results are in good agreement with the numerical simulation results in the high-conductivity environment. However, discrepancies were observed in the low-conductivity environment because the galvanic effect between the fin, tube, and joint was small. Although the polarization state changed according to the electrolyte conductivity, the total net current of the tube, which is related to the corrosion rate, was lower in the case of the anodic joint. Thus, the Al fin-tube heat exchanger assembled with the anodic joint was superior to the exchanger assembled with the cathodic joint. This result was verified by SWATT, and the leakage time of the Al fin-tube heat exchanger assembled with the anodic joint was 42% longer than that of the exchanger assembled with the cathodic joint. Due to the good correlation with the practical corrosion test results, numerical simulation of the multi-galvanic situation can be applied to improve the corrosion design of products.