Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “diffusion behavior”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,081 records · Page 60

Modeling of HTTF test PG-26 using RELAP5-3D and SAM

The High Temperature Test Facility (HTTF) at Oregon State University (OSU) is a scaled integral effects experiment designed to investigate transient behavior in high-temperature gas-cooled nuclear reactors (HTGR) with prismatic fuel and reflector blocks [1]. Several tests have been completed at the HTTF including depressurized conduction cooldown (DCC) and pressurized conduction cooldown (PCC) transients. This summary reports on the analysis of test PG-26 using the INL system code REALP5-3D [2] as well as the ANL system code SAM [3]. Test PG-26 is a progression of the Double Ended Inlet-Outlet Crossover Duct Break transient that is referred to as a DCC [4]. Core initial conditions (i.e., before the DCC started) have been met using low power (<100 kW) and two of ten available electric heater banks. The DCC transient was initiated during the 50th hour of the test. The break valves were opened, and hot helium from the core and cold helium from the reactor cavity simulation tank (RCST) started mixing. The gases flowed in a countercurrent fashion, where the top half of the hot duct contained hot helium that flowed in one direction and cold helium that flowed in the other direction in the bottom half of the duct. After the pressure and density reached equilibrium, the event entered a diffusion mode. The onset of a reverse natural circulation was not observed during the DCC period of the test.

97 MATHEMATICS AND COMPUTING↗

Fluid modeling of low-temperature plasmas

Fluid models are essential for understanding and predicting low-temperature plasma (LTP) behavior in various scientific and industrial settings. This paper provides an introductory tutorial on fluid modeling of LTPs, covering model formulation, implementation, and computational simulations. The tutorial focuses on five main components of the formulation of LTP fluid models: fluid flow, energy, chemistry, electromagnetism, and material properties, as well as in essential aspects of model implementations, including multiscale phenomena, multiphysics coupling, and numerical convergence. Designed for students and early-career researchers, this work offers a practical foundation for developing and using fluid models, from in-house computational codes to commercial software, bridging fundamental theory with real-world applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Tailoring Redox Active Ligands for Probing the Reactivity of Actinides

With this project, we aim to further enhance our understanding of fundamental f-element chemistry, including electronic structure and bonding behaviors. Key project goals include the characterization of actinide complexes bearing redox-active ligands, including those of lower-oxidation state uranium, and examining how structural changes to the ligand and coordination sphere affect structure and bonding. Redox-active and redox non-innocent ligands have been used widely in transition metal chemistry, but uranium complexes with redox-active ligands are surprisingly rare. We recently reported the synthesis of the redox-active ligand, “phen-BIAN” (N,N’-bis(iminophenol)acenaphthene), in an investigation of the electronic behavior of uranium and other actinides. This ligand framework borrows features from two classes of Schiff base ligands—the tetradentate O-N-N-O binding pocket from salens, and the redox-active α-diimine unit and backbone from Ar-BIANs (N,N’-bis[(aryl)imino]acenaphthenes, which can accept up to four electrons upon reduction. To continue this work, the electronics of uranyl and thorium complexes will be probed using the reduced forms of these ligands, using these in the preparation of lower-oxidation state and non-oxo uranium complexes. Additional members of the phen-BIAN family will be synthesized, including naphthol and thiol derivative. With these we can characterize the influence of substitutent groups on the redox activity, coordination geometry, and covalent interactions. The work proposed here entails extensive preparation, structural and spectroscopic characterization, of an array of actinide complexes. This will allow us to address probative questions about the nature of actinide bonding, the degree of covalency, the validity of lanthanides as models for the actinides, hard-soft interactions, magnetic interactions between metals in bimetallic complexes and what combination of electronic and steric affects produce an actinide selective ligand. The latter will assist greatly in providing a benchmark against which to compare and evaluate the distinct behaviors exhibited by 5f-element complexes. This will enable us to learn about the differences in binding between the 4f lanthanides and the 5f elements like uranium and plutonium. Such difference in binding can be exploited in separations and the development of new materials.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Final Technical Report on Investigation of Short-range Ordering in Transition Metal Compounds by Diffuse Scattering

Future energy needs and sustainability require new materials with novel properties for such applications as energy production, storage and transport and microelectronics. Quantum materials with several competing interactions at the electronic level offer tremendous opportunity to discover, design and tune properties for such applications. Although it has been recognized that small deviations in atomic positions, driven by the competing electronic interactions, in crystalline materials can have significant impact on properties, it remains a challenge to accurately characterize such deviations (short-range order) due to the lack of advanced instruments and analysis tools to characterize them. This project used the powerful neutron and x-ray scattering instruments recently developed at the DOE user facilities to address this challenge. New methods and efficient analysis tools were employed to uncover the hidden ordering that is behind the unique properties of transition metal compounds. One example of hidden order revealed by this project comes from vanadium (IV) oxide (VO 2 ) and related compounds. The project found that the chemical bonds in VO 2 compete against each other, unlike most crystalline compounds where the chemical bonds cooperate to yield the ordered structure. This helps explain why different measurements yielded different, competing pictures about the nature of VO 2 , which has led to disagreement about where its physical properties come from. Another example from this project of hidden ordering comes from a class of compounds known as condensed Chevrel phases, which are superconducting compounds that are generally believed to be non-magnetic. The neutron scattering methods used in this project revealed evidence of magnetism, which usually does not coexist with superconductivity. In this case, magnetism is believed to be one part of a special type of electronic behavior that can occur in compounds that have one-dimensional bonding character.

36 MATERIALS SCIENCE↗

Corrosion of 316 L stainless steel under the natural circulation of molten NaCl-MgCl 2 salt

The corrosion behavior of 316 L stainless steel (SS) was studied via the natural circulation of molten eutectic NaCl-MgCl 2 salt through a microloop. The post-corrosion tested 316 L SS microloop sections were characterized with microscopy techniques to determine the microstructural and microchemical changes that occurred at the alloy/salt interface. It was found that 316 L SS showed heterogeneous dissolution at the hot-leg, whereas deposition of corrosion products occurred at the cold-leg. For the first time, experimentally obtained molten salt flow-induced corrosion of 316 L SS results were combined with computational thermodynamic-kinetic models to validate the dissolution and deposition in terms of elemental compositional changes at the alloy/salt interface. The thermodynamic-kinetic modeling predicted that the heterogeneous dissolution of Cr and Fe from the hot-leg section of 316 L SS persisted throughout the salt circulation. The model also estimated that, despite Cr deposition starting earlier than Fe, the total redeposition of Fe is expected to be significantly greater than that of Cr over the circulation of salt. Furthermore, the modeling accurately predicted the subsurface enrichment of Mo which is attributed to the reduced Cr activity and the relatively higher diffusion rate of Mo within the alloy matrix. Here, the agreement between modeling and experimental results confirms that Fe chlorides dissolve at the hot-leg and subsequently deposit at the cold-leg due to activity changes driven by the thermal gradient. By contrast, Cr was not detected in the cold-leg deposits, which is attributed to its weaker temperature dependence on activity, limiting its redeposition under these conditions.

36 - MATERIALS SCIENCE↗

Creep behavior of high-entropy alloys: A critical review

High-entropy alloys (HEAs) have been extensively investigated during the last two decades. While substantial progress has been made in understanding their phase stability, microstructure, and deformation mechanisms at room and cryogenic temperatures, the long-term creep behavior (>100 h) of HEAs at high temperatures (>0.6 T m , where T m is the melting temperature) remains relatively underexplored. This knowledge gap is critical, as many engineering applications, including those for power generation and propulsion, require materials with good creep resistance to maintain structural integrity over extended service lifetimes. This review provides a focused and critical assessment of the current understanding of high-temperature deformation and creep behavior of HEAs, with particular attention paid to face-centered cubic HEAs and body-centered cubic refractory HEAs. The underlying deformation mechanisms governing their creep response and the influence of phase stability at elevated temperatures are examined in detail. Recent studies reveal mechanistic differences between HEAs and conventional dilute alloys that do not always lead to improved creep resistance belying their initial promise. Based on these findings, we discuss the challenges in designing HEAs for high-temperature structural applications and outline future research directions that may lead to creep-resistant HEAs.

Creep↗

Tuning Solvation Dynamics of Electrolytes at Their Eutectic Point Through Halide Identity

Deep eutectic solvents (DESs) are regarded as highly promising solvent systems for redox flow batteries. DESs, composed of choline halides (ChX, X = F − , Cl − , Br − , I − ) and ethylene glycol (EG), exhibit distinct physicochemical properties at their eutectic points, including halide-dependent phase behavior, viscosity, polarity, conductivity, and solvation dynamics. In this study, we investigate the effects of the halide identity on the solvation properties of ChX:EG mixtures at varying mol % of ChX salt content. The solvatochromic polarity based on ET(30) measurements indicates higher polarity for larger halides (I − > Br − ) than for smaller halides (Cl − > F − ), which exhibit larger compensating solvation shells. The ionic conductivity follows the trend of the solvent fluidity (the inverse of the viscosity), namely ChCl > ChBr > ChI > ChF, influenced by the ion mobility and solvodynamic radii. Measurements of the liquidus temperatures (T L ) reveal that the system with ChCl exhibits the deepest eutectic point (at ~20 mol % ChCl), while ChBr and ChI have shallower minima at ~10 mol % ChBr and ~3 mol % ChI, respectively. ChF does not display a eutectic transition but instead appears to readily supercool at salt concentrations above 30 mol % ChF. Consistent with the phase transition measurements, femtosecond transient absorption spectroscopy shows that in the ChCl system, the solvation dynamics become faster with an increasing salt concentration up to ~16.67 mol %, after which the dynamics slow down with further increases in the salt content. The ChF-based system exhibits similar behavior, though with slower dynamics. In contrast, the solvation dynamics of the systems containing ChBr and ChI monotonously slow down with an increasing salt concentration, in agreement with the phase transition measurements, which show that the eutectic points occur at low salt concentrations. These measurements suggest that the solvent composition and, in particular, the identity of the halide anion play a significant role in the solvation behavior of these ethylene-glycol-based DESs, offering a foundation for tuning the DES properties for specific applications.

Biochemistry & Molecular Biology↗

Field Testing of Self-Healing Metallic Coatings for Internal Corrosion Protection of Natural Gas Pipelines

Internal corrosion occurs in natural gas pipelines primarily due to the presence of water, carbon dioxide, and hydrogen sulfide. Internal corrosion can eventually result in leakage, cracks, and rupture of the pipeline. The objective of this work is to mitigate internal corrosion in steel pipelines transporting natural gas using cold spray coatings. The corrosion behavior of carbon steel coated with self-healing zinc chromium (ZnCr) and zinc niobium (ZnNb) cold spray coatings was investigated in a natural gas environment. For comparison purposes, hot-dip galvanized steel (HDGS) was tested under the same conditions as cold spray coatings. The field test was performed at the NW Natural gas storage facility in Mist, Oregon. The field test was conducted in a 6-inch diameter pipe transporting wet natural gas out of an underground storage well at 500 psi and 4.4 °C. The coupons were tested for 15 and 32 days under stagnant and flow conditions, respectively. Weight loss method was used to measure the corrosion rate of metallic coatings. Post-corrosion surface characterization was performed on the specimen using a scanning electron microscopy (SEM) equipped with energy dispersive X-ray spectroscopy (EDS). Crystalline phases were determined by X-ray diffraction (XRD). The field-test results confirmed that the metallic coatings provided corrosion protection of carbon steel exposed to wet natural gas under stagnant and flow conditions. The formation of the ZnCO 3 layer on top of ZnNb and ZnCr coatings led to passivation of the coatings which helps to reduce the self-corrosion. These layers form a barrier for diffusion of corrosive species to the surface.

03 NATURAL GAS↗

AGC-2 Specimen Post Irradiation Data Package Report

This report documents results of the post-irradiation examination material property testing of the creep, control, and piggyback specimens from the irradiation creep capsule Advanced Graphite Creep (AGC)-2 are reported. This is the second of a series of six irradiation test trains planned as part of the AGC experiment to fully characterize the neutron irradiation effects and radiation creep behavior of current nuclear graphite grades. The AGC-2 capsule was irradiated in the Idaho National Laboratory Advanced Test Reactor at a nominal temperature of 600°C and to a peak dose of 5 dpa (displacements per atom). One half of the creep specimens were subjected to mechanical stresses (an applied stress of either 13.8, 17.2, or 20.7 MPa) to induce irradiation creep. All post-irradiation testing and measurement results are reported with the exception of the irradiation mechanical strength testing, which is the last destructive testing stage of the irradiation testing program. Material property tests were conducted on specimens from 15 nuclear graphite grades using a similar loading configuration as the first AGC capsule (AGC-1) to provide easy comparison between the two capsules. However, AGC-2 contained an increased number of specimens (i.e., 487 total specimens irradiated) and replaced specimens of the minor grade 2020 with the newer grade 2114. The data reported include specimen dimensions for both stressed and unstressed specimens to establish the irradiation creep rates, mass and volume data necessary to derive density, elastic constants (Young’s modulus, shear modulus, and Poisson’s ratio) from ultrasonic time of flight velocity measurements, Young’s modulus from the fundamental frequency of vibration, electrical resistivity, and thermal diffusivity and thermal expansion data from 100–500°C. No data outliers were determined after all measurements were completed. A brief statistical analysis was performed on the irradiated data and a limited comparison between pre- and post-irradiation properties is presented. A more complete evaluation of trends in the material property changes, as well as irradiation-induced creep due to irradiation, temperature, and applied load on specimens will be discussed in later AGC-2 post-irradiation examination analysis reports.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Reactive particle-tracking solutions to a benchmark problem on heavy metal cycling in lake sediments

Geochemical systems are known to exhibit highly variable spatiotemporal behavior. This may be observed both in non-smooth concentration curves in space for a single sampling time and also in variability between samples taken from the same location at different times. However, most models that are designed to simulate these systems provide only single-solution smooth curves and fail to capture the noise and variability seen in the data. We apply a recently developed reactive particle-tracking method to a system that displays highly complex geochemical behavior. When the method is made to most closely resemble a corresponding Eulerian method, in its unperturbed form, we see near-exact match between solutions of the two models. More importantly, we consider two approaches for perturbing the model and find that the spatially-perturbed condition is able to capture a greater degree of the variability present in the data. This method of perturbation is a task to which particle methods are uniquely suited and Eulerian models are not well-suited. Additionally, because of the nature of the algorithm, noisy spatial gradients can be highly resolved by a large number of mobile particles, and this incurs negligible computational cost, as compared to expensive chemistry calculations.

54 ENVIRONMENTAL SCIENCES↗

Ion irradiation study of lithium silicates for fusion blanket applications

Here, this study reports on the microstructural and compositional evolutions, deuterium release and lithium loss behavior in lithium silicates. Pellets of biphasic Li 4 SiO 4 and Li 2 SiO 3 were fabricated using hot pressing of powders. Sequential irradiation with Si + , He + and D + ions was performed up to 773 K to emulate one-year 6 Li burnup in 6 Li 4 SiO 4 inside the SlimCS DEMO. Crystalline Li 4 SiO 4 phase was not observed in the irradiated depth region of all the samples in this study. In spite of full amorphization in the near-surface region, Li 2 SiO 3 irradiated to a high dose at 773 K remained crystalline in the damage peak region, suggesting that Li 2 SiO 3 is more irradiation resistant to amorphization than Li 4 SiO 4 . Radiolysis may be primarily responsible for silicate decomposition and amorphization. Data from this study also show that D release from the pellet is very efficient during ion irradiation at 773 K, which is accompanied with a significant Li loss. Thermal annealing for 10 min at 773 K for the room-temperature irradiated pellet leads to a nearly complete D release without an observable Li loss. A concept of Ni coating on tritium breeder materials is discussed with supporting data to minimize Li loss without a significant impact on tritium diffusion and release.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Efficient quadrature rules for finite element discretizations of nonlocal equations

Here, in this paper, we design efficient quadrature rules for finite element (FE) discretizations of nonlocal diffusion problems with compactly supported kernel functions. Two of the main challenges in nonlocal modeling and simulations are the prohibitive computational cost and the nontrivial implementation of discretization schemes, especially in three-dimensional settings. In this work, we circumvent both challenges by introducing a parametrized mollifying function that improves the regularity of the integrand, utilizing an adaptive integration technique, and exploiting parallelization. We first show that the “mollified” solution converges to the exact one as the mollifying parameter vanishes, then we illustrate the consistency and accuracy of the proposed method on several two- and three-dimensional test cases. Furthermore, we demonstrate the good scaling properties of the parallel implementation of the adaptive algorithm and we compare the proposed method with recently developed techniques for efficient FE assembly.

97 MATHEMATICS AND COMPUTING↗

Efficient quadrature rules for finite element discretizations of nonlocal equations

In this paper we design efficient quadrature rules for finite element discretizations of nonlocal diffusion problems with compactly supported kernel functions. Two of the main challenges in nonlocal modeling and simulations are the prohibitive computational cost and the nontrivial implementation of discretization schemes, especially in three-dimensional settings. In this work we circumvent both challenges by introducing a parametrized mollifying function that improves the regularity of the integrand, utilizing an adaptive integration technique, and exploiting parallelization. We first showthat the “mollified” solution converges to the exact one as the mollifying parameter vanishes, then we illustrate the consistency and accuracy of the proposed method on several two- and three-dimensional test cases. Furthermore, we demonstrate the good scaling properties of the parallel implementation of the adaptive algorithm and we compare the proposed method with recently developed techniques for efficient finite element assembly.

97 MATHEMATICS AND COMPUTING↗

Effect of microphase separation on the limiting current density in hybrid organic-inorganic copolymer electrolytes

Hybrid organic-inorganic block copolymer electrolytes are of interest to enable batteries containing lithium metal anodes. The conductive block is a standard polymer electrolyte of poly(ethylene oxide) and the mechanically rigid block is an inorganic poly(acryloisobutyl polyhedral oligomeric silsesquioxane) polymer. Here, in this paper, we compare a poly(acryloisobutyl polyhedral oligomeric silsesquioxane)-b-poly(ethylene oxide)-b-poly(acryloisobutyl polyhedral oligomeric silsesquioxane) (POSS-PEO-POSS) triblock copolymer and a poly(ethylene oxide)-b- poly(acryloisobutyl polyhedral oligomeric silsesquioxane) (PEO-POSS) diblock copolymer mixed with lithium bis(trifluoromethanesulfonyl)imide salt. We have experimentally measured the limiting current density in lithium symmetric cells containing hybrid organic-inorganic electrolytes at 90 °C. The cells were polarized at a large range of applied current density. The diblock copolymer electrolyte exhibited a clear plateau in cell potential at all current densities below the limiting current density. At low applied current density, the triblock copolymer electrolyte also exhibited a clear plateau in cell potential. At currents approaching the limiting current density, the triblock copolymer electrolyte exhibited an underdamped potential profile. The cell potential did not reach a plateau at current densities above the limiting current in both systems. The diblock and triblock copolymer electrolytes were fully characterized using electrochemical methods to determine the ionic conductivity, cation current fraction, salt diffusion coefficient, and open circuit voltage as a function of salt concentration. Cell potential and salt concentration as functions of position in the cell at various current densities were calculated using Newman's concentrated solution theory. The theoretical limiting current density was calculated to be the current density at which salt is depleted at the cathode. We see quantitative agreement between experimental measurements and theoretical predictions for the limiting current density in the diblock copolymer electrolyte which has an ordered structure at all salt concentrations, while the experimental limiting current density is lower than the theoretical prediction for the triblock copolymer electrolyte, which exhibits a disordered morphology at high salt concentrations.

25 ENERGY STORAGE↗

Machine learning molecular dynamics simulations toward exploration of high-temperature properties of nuclear fuel materials: case study of thorium dioxide

Predicting materials properties of nuclear fuel compounds is a challenging task in materials science. Their thermodynamical behaviors around and above the operational temperature are essential for the design of nuclear reactors. However, they are not easy to measure, because the target temperature range is too high to perform various standard experiments safely and accurately. Moreover, theoretical methods such as first-principles calculations also suffer from the computational limitations in calculating thermodynamical properties due to their high calculation-costs and complicated electronic structures stemming from f-orbital occupations of valence electrons in actinide elements. Here, we demonstrate, for the first time, machine-learning molecular-dynamics to theoretically explore high-temperature thermodynamical properties of a nuclear fuel material, thorium dioxide. The target compound satisfies first-principles calculation accuracy because f-electron occupation coincidentally diminishes and the scheme meets sampling sufficiency because it works at the computational cost of classical molecular-dynamics levels. We prepare a set of training data using first-principles molecular dynamics with small number of atoms, which cannot directly evaluate thermodynamical properties but captures essential atomistic dynamics at the high temperature range. Then, we construct a machine-learning molecular-dynamics potential and carry out large-scale molecular-dynamics calculations. Consequently, we successfully access two kinds of thermodynamic phase transitions, namely the melting and the anomalous λ transition induced by large diffusions of oxygen atoms. Furthermore, we quantitatively reproduce various experimental data in the best agreement manner by selecting a density functional scheme known as SCAN. Our results suggest that the present scale-up simulation-scheme using machine-learning techniques opens up a new pathway on theoretical studies of not only nuclear fuel compounds, but also a variety of similar materials that contain both heavy and light elements, like thorium dioxide.

74 ATOMIC AND MOLECULAR PHYSICS↗

Nanoconfinement of Carbon Dioxide within Interfacial Aqueous/Ionic Liquid Systems

Nanoporous, gas-selective membranes have shown encouraging results for the removal of CO 2 from flue gas, yet the optimal design for such membranes is often unknown. Therefore, we used molecular dynamics simulations to elucidate the behavior of CO 2 within aqueous and ionic liquid (IL) systems ([EMIM][TFSI] and [OMIM][TFSI]), both confined individually and as an interfacial aqueous/IL system. We found that within aqueous systems the mobility of CO 2 is reduced due to interactions between the CO 2 oxygens and hydroxyl groups on the pore surface. Within the IL systems, we found that confinement has a greater effect on the [EMIM][TFSI] system as opposed to the [OMIM][TFSI] system. Paradoxically, the larger and more asymmetrical [OMIM]+ molecule undergoes less efficient packing, resulting in fewer confinement effects. Finally, free energy surfaces of the nanoconfined aqueous/IL interface demonstrate that CO 2 will transfer spontaneously from the aqueous to the IL phase.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Origin of Rapid Delithiation In Secondary Particles Of LiNi 0.8 Co 0.15 Al 0.05 O 2 and LiNi y Mn z Co 1− y − z O 2 Cathodes

Abstract Most research on the electrochemical dynamics in materials for high‐energy Li‐ion batteries has focused on the global behavior of the electrode. This approach is susceptible to misleading analyses resulting from idiosyncratic kinetic conditions, such as surface impurities inducing an apparent two‐phase transformation within LiNi 0.8 Co 0.15 Al 0.05 O 2 . Here, nano‐focused X‐ray probes are used to measure delithiation operando at the scale of secondary particle agglomerates in layered cathode materials during charge. After an initial latent phase, individual secondary particles undergo rapid, stochastic, and largely uniform delithiation, which is in contrast with the gradual increase in cell potential. This behavior reproduces across several layered oxides. Operando X‐ray microdiffraction (‐XRD) leverages the relationship between Li content and lattice parameter to further reveal that rate acceleration occurs between Li‐site fraction ( x Li ) ≈0.9 and ≈0.5 for LiNi 0.8 Co 0.15 Al 0.05 O 2 . Physics‐based modeling shows that, to reproduce the experimental results, the exchange current density ( i 0 ) must depend on x Li , and that i 0 should increase rapidly over three orders of magnitude at the transition point. The specifics and implications of this jump in i 0 are crucial to understanding the charge‐storage reaction of Li‐ion battery cathodes.

25 ENERGY STORAGE↗

Characterization of electromigration-induced short-range stress development in Al(0.25 at. % Cu) conductor line

Scanning x-ray microbeam topography and fluorescence experiments were conducted in situ to study the electromigration behavior of a 0.5 μm thick, 10 μm wide, and 200 μm long Al(0.25 at. % Cu) conductor line with 1.5 μm-thick SiO2 passivation on a single crystal Si substrate. The strain sensitivity of x-ray topography measurement allowed detailed examination of the electromigration-induced stress distribution and evolution in the conductor line in response to the depletion of Cu solute early in the electromigration process. Upon electromigration at 0.4 MA/cm2 and 303 °C, a short-range stress gradient was quickly induced by Al migration in the Cu-depleted cathode region to counteract further Al flow. The stress gradient was fully developed during the 5.3 h incubation time, extending over the critical Blech length of about 66 μm from the cathode end. Plastic deformation then occurred at the downstream end of the Cu-depleted region. The preferential electromigration of Cu did not cause detectable stress change outside the Cu-depleted region, except for the significant stress development from the Al2Cu precipitation at the anode end which appeared to initiate the fracture in the passivation. Preliminary finite difference modeling was undertaken to simulate the experimental observations, from which important parameters dictating electromigration in Al(Cu) line were extracted: an apparent effective valence of −5.6 and −1.9 for Cu and Al in Al(Cu), respectively, and a critical Cu concentration of 0.16 at. % above which Al grain boundary diffusion is effectively blocked.

Physics↗