Morphological spatial patterns in a reaction diffusion model for metal growth

  • Received: 01 May 2009 Accepted: 29 June 2018 Published: 01 April 2010
  • MSC : 35K57, 62P30, 65M.

  • In this paper a reaction-diffusion system modelling metal growth processes is considered, to investigate - within the electrodeposition context- the formation of morphological patterns in a finite two-dimensional spatial domain. Nonlinear dynamics of the system is studied from both the analytical and numerical points of view. Phase-space analysis is provided and initiation of spatial patterns induced by diffusion is shown to occur in a suitable region of the parameter space. Investigations aimed at establishing the role of some relevant chemical parameters on stability and selection of solutions are also provided. By the numerical approximation of the equations, simulations are presented which turn out to be in good agreement with experiments for the electrodeposition of Au-Cu and Au-Cu-Cd alloys.

    Citation: Benedetto Bozzini, Deborah Lacitignola, Ivonne Sgura. Morphological spatial patterns in a reaction diffusion model for metal growth[J]. Mathematical Biosciences and Engineering, 2010, 7(2): 237-258. doi: 10.3934/mbe.2010.7.237

    Related Papers:

    [1] Guangyu Sui, Meng Fan, Irakli Loladze, Yang Kuang . The dynamics of a stoichiometric plant-herbivore model and its discrete analog. Mathematical Biosciences and Engineering, 2007, 4(1): 29-46. doi: 10.3934/mbe.2007.4.29
    [2] Md Nazmul Hassan, Angela Peace . Mechanistically derived Toxicant-mediated predator-prey model under Stoichiometric constraints. Mathematical Biosciences and Engineering, 2020, 17(1): 349-365. doi: 10.3934/mbe.2020019
    [3] Ming Chen, Menglin Gong, Jimin Zhang, Lale Asik . Comparison of dynamic behavior between continuous- and discrete-time models of intraguild predation. Mathematical Biosciences and Engineering, 2023, 20(7): 12750-12771. doi: 10.3934/mbe.2023569
    [4] Ilse Domínguez-Alemán, Itzel Domínguez-Alemán, Juan Carlos Hernández-Gómez, Francisco J. Ariza-Hernández . A predator-prey fractional model with disease in the prey species. Mathematical Biosciences and Engineering, 2024, 21(3): 3713-3741. doi: 10.3934/mbe.2024164
    [5] Md Nazmul Hassan, Kelsey Thompson, Gregory Mayer, Angela Peace . Effect of Excess Food Nutrient on Producer-Grazer Model under Stoichiometric and Toxicological Constraints. Mathematical Biosciences and Engineering, 2019, 16(1): 150-167. doi: 10.3934/mbe.2019008
    [6] P. Auger, N. H. Du, N. T. Hieu . Evolution of Lotka-Volterra predator-prey systems under telegraph noise. Mathematical Biosciences and Engineering, 2009, 6(4): 683-700. doi: 10.3934/mbe.2009.6.683
    [7] Aniello Buonocore, Luigia Caputo, Enrica Pirozzi, Amelia G. Nobile . A non-autonomous stochastic predator-prey model. Mathematical Biosciences and Engineering, 2014, 11(2): 167-188. doi: 10.3934/mbe.2014.11.167
    [8] Ya Li, Z. Feng . Dynamics of a plant-herbivore model with toxin-induced functional response. Mathematical Biosciences and Engineering, 2010, 7(1): 149-169. doi: 10.3934/mbe.2010.7.149
    [9] Paulo Amorim, Bruno Telch, Luis M. Villada . A reaction-diffusion predator-prey model with pursuit, evasion, and nonlocal sensing. Mathematical Biosciences and Engineering, 2019, 16(5): 5114-5145. doi: 10.3934/mbe.2019257
    [10] Xinru Zhou, Xinmiao Rong, Meng Fan, Josué-Antonio Nescolarde-Selvaa . Stoichiometric modeling of aboveground-belowground interaction of herbaceous plant. Mathematical Biosciences and Engineering, 2019, 16(1): 25-55. doi: 10.3934/mbe.2019002
  • In this paper a reaction-diffusion system modelling metal growth processes is considered, to investigate - within the electrodeposition context- the formation of morphological patterns in a finite two-dimensional spatial domain. Nonlinear dynamics of the system is studied from both the analytical and numerical points of view. Phase-space analysis is provided and initiation of spatial patterns induced by diffusion is shown to occur in a suitable region of the parameter space. Investigations aimed at establishing the role of some relevant chemical parameters on stability and selection of solutions are also provided. By the numerical approximation of the equations, simulations are presented which turn out to be in good agreement with experiments for the electrodeposition of Au-Cu and Au-Cu-Cd alloys.


    1. Introduction

    Al-Ag supersaturated solid solution evolves towards the equilibrium state following the sequence [1,2,3,4]:

    Supersaturated solid solution → Guinier-Preston (GP) zones → metastable γ' phase → equilibrium γ phase

    The Guinier-Preston zones (GP), consisting of silver atom clusters, are coherent with the matrix. The metastable γ' (Ag2Al) phase is semi-coherent with the matrix and the equilibrium γ (Ag2Al) phase, is incoherent with the matrix. Precipitation starts from the formation of GP zones, which are isomorphous with the matrix and, therefore, have a lower interfacial energy than intermediate or equilibrium precipitate phases which possess a distinct crystal structure. In fact, the GP zones are fully coherent with the matrix and therefore have a very low interfacial energy, whereas the γ' (Ag2Al) phase is semi-coherent with the matrix which can only form with high-energy semi-incoherent interfaces. Therefore, despite the fact that the driving force for precipitation of GP zones is less than for the metastable γ' (Ag2Al) phase, the barrier to nucleation is still less, and the zones nucleate most rapidly.

    The effect of microalloying elements on the behaviour of age-hardenable alloys, such as, Al-Ag, is an interesting physical problem addressing the mechanisms of transport and aggregation of the solute. Trace elements have been found to exert a disproportionate influence on the structure and properties of Al alloys compared with the amounts added which may be less than 0.1 molar fraction. Most trace element effects arise because they modify the nucleation and the growth characteristics of the phases which form during precipitation such the GP zones and the metastable γ' phase in the case of Al-Ag alloys [5,6,7,8,9]. It is well known that the formation of the GP zones in aluminium alloys is closely linked with the excess vacancies. A number of model of GP zones precipitation assisted by vacancies has been developed by several authors [10,11,12]. GP zones formation is governed by a transport mechanism of solute atoms by solute atom-vacancy complexes. The high binding energy between tin atoms and vacancies leads to the formation of vacancy-tin atom pairs and silver atom-vacancy-tin atom complexes.

    In a supersaturated Al-Ag alloy, the precipitation of the GP zones and the metastable γ' phase is responsible of a variation of the physical properties of the alloy and particularly, the mechanical properties. The evolution of all these properties is linked with the nature of the precipitate particles and their precipitation kinetics. Our purpose is to study the effect of the addition of traces of tin on the precipitation of the GP zones and the metastable γ' phase in Al-10at%Ag alloy using a method based on hardness measurements and differential scanning calorimetry.


    2. Materials and Method

    Al-10at%Ag (31wt%) and Al-10at%Ag-0.05at%Sn (2wt%) alloys were prepared by melting 99.99%, 99.99% and 99.99% pure aluminum, silver and tin, respectively, under argon protection. After homogenization 15 days at 540 °C and ice water quenching, the alloys are cutinto platelet specimens which are mechanically polished, homogenized 6 hours at 540 °C and quenched into ice water. A P1710 Phillips X-ray ray powder diffractometer, utilizing the monochromatic CuKα1 radiation, is used for the characterization of our alloys. The Vickers microhardness measurements were carried out under a load of 100 g on specimen treated during different times at different aged temperatures (90, 125, 150, 170, and 200 °C) and quenched into ice water. The Vickers hardness measurements were made using a microhardness tester type SHIMADZU provided with a square pyramidal penetrator. The differential scanning calorimetry analysis, using a NETZSCH 200PC DSC microcalorimeter, were performed on disk shape of 4 mm diameter and 3 mm thickness. The thermal cycle applied consists of a heating with a rate of 5 °C·min−1.


    3. Results and Discussion


    3.1. Characterization of the Solid Solutions

    The as quenched solids solutions, Al-Ag and Al-Ag-Sn, are characterized by X-ray diffraction on powder specimen. It is well known that an as quenched substitutional solid solution is in a disordered state. In the case of aluminium based solids solutions, the X-ray ray powder diffraction spectra give peaks diffraction at the analogous positions of those given by the aluminium. The Al, Al-Ag and Al-Ag-Sn powder diffraction spectra show that the successive peaks are at analogous positions (Figure 1) and their integrated intensities vary in the same ratio (Table 1). The lattice parameters of Al, Al-Ag and Al-Ag-Sn, determined using the Nelson-Riley [13] extrapolation function, are 4.0462 Å, 4.0478 Å and 4.0467 Å with an error estimated at 0.05% respectively. In such a case, there is practically no difference between these lattice parameters. This is due to the fact that there is a little difference between the Al and the Ag atomic radius and the Sn, which have a larger atomic radius, is added in traces proportion.

    Figure 1. X-ray diffractograms of Al, Al-Ag and Al-Ag-Sn as quenched solid solutions.
    Table 1. Integrated intensities ratios.
    hkl Al (Ihkl/I111) Al-Ag (Ihkl/I111) Al-Ag-Sn (Ihkl/I111)
    111 100 100 100
    200 50 51 53
    220 21 30 29
    311 20 21 21
    222 6 6 6
     | Show Table
    DownLoad: CSV

    3.2. Hardening Evolution at 90, 125 and 150 °C

    The isotherm curves of hardness, established at 90, 125 and 150 °C show a first step of hardening due to the GP zones precipitation and a second step due to the precipitation of the metastable γ' phase (Figures 2, 3 and 4). The obtained degree of hardening depends on the volume fraction, the structure of the precipitates and the nature of the interface between the metastable phases and the aluminum matrix [10,14,15,16,17,18,19].

    Figure 2. Isotherm hardness curves at 90 °C.
    Figure 3. Isotherm hardness curves at 125 °C.
    Figure 4. Isotherm hardness curves at 150 °C.

    The intermediate bearing corresponds to the metastable equilibrium state of the precipitation of the GP zones during which their volume fraction is maximum. The softening is due to the coarsening of the particles of the γ' phase and to the precipitation of the equilibrium γ phase. These isotherm curves of hardness show that tin atoms retard the establishment of the metastable equilibrium state of the GP zones precipitation by trapping some of quenched in vacancies. Due to the high binding energy between tin atoms and vacancies, tin atom-vacancy pairs and silver atom-vacancy-tin atom complexes are formed, thus reducing the number of free vacancies available to promote the diffusion of Ag atoms for the formation of the GP zones.

    Due to the relatively high binding energy between tin atoms and vacancies, estimated to be about 0.43 eV [20,21] , compared with that of the silver atoms-vacancy binding energy , estimated to be about 0.25 eV [22], tin atom-vacancy pairs [23,24] are formed and, due to the strong interaction between the tin atoms and the vacancies and the interaction between the silver atoms and the vacancies, silver atom-vacancy-tin atom complexes are formed, thus reducing the number of free vacancies available to promote the diffusion of Ag atoms for the formation of the GP zones.

    These curves also show that tin addition stimulates the precipitation of the metastable γ' phase. This is due to the fact that tin atoms are absorbed at the interface γ'-matrix and reduce the interfacial energy required for precipitate nucleation [25,26]. An alternative explanation is that the clustering of tin elements, which have a larger atomic volume than matrix aluminium, would create a compressive volume strain in the matrix, thus attracting a high concentration of vacancies and silver atoms to the interfacial region, providing an excellent condition for the nucleation of the metastable γ' phase [27].


    3.3. Precipitation Kinetics of the GP Zones

    During the precipitation of the GP zones, the transformed fraction, F, which represents the ratio between the volume occupied by the GP zones at a time t and their volume at the metastable equilibrium state, is given by the Merle relatio [28]:

    Hv(t)=FHv(metastableequilibrium)+(1F)Hv(0)

    where Hv(0) is the as quenched hardness, Hv(t) is the hardness of the alloy at the time t during the precipitation of the GP zones, and Hv(metastable equilibrium), the hardness of the alloy at the metastable equilibrium state of the GP zones precipitation. The incubation times, which the determined values by extrapolation varies between 0.01 and 0.02 hours, compared with the necessary times to reach the metastable equilibrium state, are very short and are characteristics of a rapid nucleation in the two both alloys, Al-Ag and Al-Ag-Sn, because of the high supersaturation of the quenched in vacancies (Figures 5, 6 and 7).

    Figure 5. Transformed fraction during GP zones precipitation in Al-Ag and Al-Ag-Sn alloys at 90 °C.
    Figure 6. Transformed fraction during GP zones precipitation in Al-Ag and Al-Ag-Sn alloy at 125 °C.
    Figure 7. Transformed fraction during GP zones precipitation in Al-Ag and Al-Ag-Sn alloy at 150 °C.

    The curves of the variations of ln(ln(1/1 − F)) versus ln(t), show that the growth stage obeys to the JMAK (Johnson-Mehl-Avrami-Kolmogorov) [29,30,31] law, F = 1 - exp[(−kt)n], of the growth controlled by solute atom diffusion, where n and k are the growth parameters (Figures 8, 9 and 10).

    Figure 8. Determination of the growth parameters for Al-Ag and Al-Ag-Sn alloys at 90 °C.
    Figure 9. Determination of the growth parameters for Al-Ag and Al-Ag-Sn alloys at 125 °C.
    Figure 10. Determination of the growth parameters for Al-Ag and Al-Ag-Sn alloys at 150 °C.

    The values of n are characteristics of a heterogeneous precipitation in the two alloys while the values of k, which characterizes the transformation kinetic, show a delay of the precipitation of the GP zones in Al-Ag-Sn alloy (Table 2).

    Table 2. Values of the growth parameters.
    T (°C) Al-Ag Al-Ag-Sn
    n k (s−1) n k (s−1)
    90 0.80 22 × 10−4 0.84 10 × 10−4
    125 0.82 31 × 10−4 1.0 24 × 10−4
    150 1.08 44 × 10−4 0.94 30 × 10−4
     | Show Table
    DownLoad: CSV

    3.3.1. Activation Energy

    The apparent activation energy, Q, is determined using the Arrhenius law k = A·exp(−Q/RT) [32] where A is a constant, T is the temperature and R is the gas constant. The slopes of the variation curves ln(k) = f(1/T) (Figure 11) gives an apparent activation energy in the order of 14 ± 1.4 kJ/mol and 24 ± 2.4 kJ/mol in the Al-Ag and the Al-Ag-Sn alloys respectively corresponding to a difference of 10 kJ/mol (0.1 eV) which confirms that the reaction of GP zones precipitation is slowest in the microalloyed alloy.

    Figure 11. Determination of the difference between activation energy in Al-Ag and Al-Ag-Sn alloys.

    3.4. Precipitation of the Metastable γ' Phase

    At 90, 125 and 150 °C, the isotherm curves of hardness show that tin addition stimulates the precipitation of the metastable γ' phase. The variations of the duration of the GP zones equilibrium metastable state, Δt, and the beginning time of the second step hardening which corresponds to the γ' particles precipitation, tγ', against 1/T obey to an Arrhenius type law A·exp(−Q/RT), where A is a constant, T is the temperature and R is the gas constant (Figures 12 and 13).

    Figure 12. Variation of the duration of the GP zones equilibrium metastable state, Δt.
    Figure 13. Variation of the beginning time of the second step hardening, tγ'.

    The apparent activation energies of the solute atom diffusion, determined from the slopes of these variations curves, show that the precipitation of the γ' phase occurs earlier in the Al-Ag-Sn alloy (Table 3).

    Table 3. Apparent activation energies determined from the duration of the GP zones equilibrium metastable state, Δt variation and from the beginning time of the second step hardening tγ' variation .
    Alloy Al-Ag Al-Ag-Sn
    Q (kJ/mol) (Δt) 50.1 ± 5.0 38.7 ± 3.9
    Q (kJ/mol) (tγ') 48.8 ± 4.9 36.7 ± 3.7
     | Show Table
    DownLoad: CSV

    The difference between these apparent activation energies, ΔQ = QAl-Ag − QAl-Ag-Sn , of the solute atom diffusion in Al-Ag and in Al-Ag-Sn alloys, determined from the Δt variation (Figure 12) and from the tγ' variation (Figure 13) are in the order of 11.4 ± 1.1 kJ/mol (0.11 ± 0.01 eV) and 12.1 ± 1.2 kJ/mol (0.12 ± 0.01 eV) respectively.

    It is explained by the fact that tin atoms are absorbed at the interface γ'-matrix and reduce the interfacial energy required for precipitate nucleation [25,26]. An alternative explanation is that the clustering of tin elements, which have a larger atomic volume than the aluminium matrix, would create a compressive volume strain in the matrix, thus attracting a high concentration of vacancies and silver atoms to the interfacial region, providing an excellent condition for the nucleation of the metastable γ' phase [27].


    3.4.1. Precipitation at 170 °C and 200 °C

    In the added alloy, the isotherm curves of hardness established at 170 °C and 200 °C (Figures 14 and 15) show, in a first stage, a light softening due, probably, to the coarsening and the dissolution of the GP zones as soon as the nuclei of the γ'phase, more thermodynamically stable, appear as observed on the DSC curves (Figure 16). At this stage, the particle size of the γ' phase produces an insufficient hardening [10,14,15,16,17,18,19]. In the second stage, the hardening is due to the growth of the particles of the γ' phase.

    Figure 14. Isotherm curves of hardness at 170 °C.
    Figure 15. Isotherm curves of hardness at 200 °C.

    In the non added alloy, as well as 170 °C and 200 °C, due to the thermal activation, it the GP zones are rapidly formed and their equilibrium metastable state, which is reached rapidly, is followed by their dissolution and the formation of the γ' phase which hardens the alloy as shown in the DSC curves (Figure 16) [10,14,15,16,17,18,19]. The differential scanning calorimetry (DSC) curves where we observe two peaks (Figure 16). The first one, endothermic, between 150 °C and 200 °C for Al-Ag alloy and between 150 °C and 180 °C for Al-Ag-Sn alloy, corresponds to the dissolution of the GP zones. The second one, exothermic, between 240 °C and 325 °C for Al-Ag alloy and between 225 °C and 310 °C for Al-Ag-Sn, corresponds to the precipitation of the metastable γ' phase. The energy associated with the dissolution of the GP zones and that associated with the precipitation of the γ' metastable phase are given by the peak area in the Table 4. The peaks temperature and the peak area are given in the Table 4.

    Figure 16. DSC thermograms of Al-Ag and Al-Ag-Sn alloys at theheating rate of 5 °C/mm.
    Table 4. DSC results.
    Alloy GP zones dissolution γ' phase precipitation
    Peak (°C) Peak area (J/g) Peak (°C) Peak area (J/g)
    Al-Ag 183 0.28 273 1.04
    Al-Ag-Sn 176 0.35 269 2.69
     | Show Table
    DownLoad: CSV

    In results, the presence of tin atoms in the alloy prevents the formation of the GP zones and promotes the precipitation of the γ' phase by trapping vacancies. In fact, at an early stage of aging treatment, due to the high binding energies between the tin atoms and vacancies, tin atom-vacancy and tin atom-vacancy-silver atom clusters are formed [20,21,22,23,24], thus effectively reducing the number of free vacancies available to promote the diffusion of Ag atoms for the formation of the GP zones, and consequently, suppressing the formation of the GP zone and promoting the precipitation of the γ' metastable phase [29,33] .


    4. Conclusion

    In the Al-10at%Ag-0.05at%Sn alloy, tin atoms form tin atom-vacancy pairs and tin atom-vacancy-silver atom complexes which retards the precipitation of the GP zones at low temperatures (90, 125, 150 °C) and suppress their precipitation at higher temperatures (170, 200 °C). At all temperatures, the presence of tin atoms stimulates the precipitation of the metastable γ' phase.


    Conflict of Interest

    The authors declare that there is no conflict of interest regarding the publication of this manuscript.


  • This article has been cited by:

    1. Haidong Fan, Alfonso H.W. Ngan, Kefu Gan, Jaafar A. El-Awady, Origin of double-peak precipitation hardening in metallic alloys, 2018, 111, 07496419, 152, 10.1016/j.ijplas.2018.07.016
    2. Asadullah Madni, Sobia Noreen, Irsah Maqbool, Faizza Rehman, Amna Batool, Prince Muhammad Kashif, Mubashar Rehman, Nayab Tahir, Muhammad Imran Khan, Graphene-based nanocomposites: synthesis and their theranostic applications, 2018, 26, 1061-186X, 858, 10.1080/1061186X.2018.1437920
  • Reader Comments
  • © 2010 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(2997) PDF downloads(535) Cited by(19)

Article outline

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog