Search NASA⌕ Search

Engineering topics

Karma, Alain

Publications and source records attributed to Karma, Alain.

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Clarifying the formation of equiaxed grains and microstructural refinement in the additive manufacturing of Ti-Cu

Controlling microstructural evolution in metallic additive manufacturing (AM) is difficult, especially in producing refined as-built grains instead of coarse, directional grains. Traditional solutions involve adding inoculants to AM feedstocks, but titanium (Ti) alloys cannot employ this approach without producing detrimental secondary phases. Ti-Cu (Ti-copper) alloys offer a solution through constitutional supercooling and/or solid state thermal cycling under AM conditions. This work analyzes a compositionally graded directed energy deposition (DED) Ti-Cu build, single-melt laser tracks, and dilatometric heat treatments to evaluate if, when, and by what mechanism(s) microstructural refinement occurs. Refinement by inoculation of unmelted powder particles was also considered. Constitutional supercooling produced no net microstructural refinement as any equiaxed dendrites which form are remelted with new deposition. This finding agreed with solidification modeling of powder bed fusion-laser beam (PBF-LB) and DED builds. Solid state thermal cycling refined microstructures only during ex-situ dilatometric heat treatments, suggesting build parameter optimization is needed to achieve refinement in-situ. Accidental heterogeneous nucleation on unmelted Ti powder, originating from the different thermophysical properties of Ti and Cu, provided the most significant microstructural refinement. This work systematically assesses the microstructural refinement mechanisms of Ti-Cu in AM builds and offers insights into microstructural control in eutectoid alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Microstructure development during rapid alloy solidification

Abstract Solidification processing of structural alloys can take place over an extremely wide range of solid–liquid interface velocities spanning six orders of magnitude, from the low-velocity constitutional supercooling limit of microns/s to the high-velocity absolute stability limit of m/s. In between these two limits, the solid–liquid interface is morphologically unstable and typically forms cellular-dendritic microstructures, but also other microstructures that remain elusive. Rapid developments in additive manufacturing have renewed the interest in modeling the high-velocity range, where approximate analytical theories provide limited predictions. In this article, we discuss recent advances in phase-field modeling of rapid solidification of metallic alloys, including a brief description of state-of-the-art experiments used for model validation. We describe how phase-field models can cope with the dual challenge of carrying out simulations on experimentally relevant length- and time scales and incorporating nonequilibrium effects at the solid–liquid interface that become dominant at rapid rates. We present selected results, illustrating how phase-field simulations have yielded unprecedented insights into high-velocity interface dynamics, shedding new light on both the absolute stability limit and the formation of banded microstructures that are a hallmark of rapid alloy solidification near this limit. We also discuss state-of-the-art experiments used to validate those insights. Graphical abstract

36 MATERIALS SCIENCE↗
DOE OSTI2023

Topology-enhanced mechanical stability of swelling nanoporous electrodes

Abstract Materials like silicon and germanium offer a 10-fold improvement in charge capacity over conventional graphite anodes in lithium-ion batteries but experience a roughly threefold volume increase during lithiation, which challenges ensuring battery integrity. Nanoporous silicon, created by liquid-metal-dealloying, is a potentially attractive anode design to mitigate this challenge, exhibiting both higher capacity and extended cycle lifetimes. However, how nanoporous structures accommodate the large volume change is unknown. Here, we address this question by using phase-field modeling to produce nanoporous particles and to investigate their elastoplastic swelling behavior and fracture. Our simulations show that enhanced mechanical stability results from the network topology consisting of ligaments connected by bulbous, sphere-like nodes. The ligaments forcefully resist elongation while the nodes, behaving like isolated spherical particles, experience large stresses driving fracture. However, being smaller compared to a sphere of the same volume as the entire nanoporous particle, the nodes are more protected against fracture.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Scaling laws for two-dimensional dendritic crystal growth in a narrow channel

Here, we investigate analytically and computationally the dynamics of two-dimensional needle crystal growth from the melt in a narrow channel. Our analytical theory predicts that, in the low supersaturation limit, the growth velocity $\textit{V}$ decreases in time $\textit{t}$ as a power law $V ~t^{–2/3}$, which we validate by phase-field and dendritic-needle-network simulations. Simulations further reveal that, above a critical channel width $Λ ≈ 5l_D$, where $l_D$ is the diffusion length, needle crystals grow with a constant $V < V_s$, where $V_s$ is the free-growth needle crystal velocity, and approaches $V_s$ in the limit $Λ \gg l_D$.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Phase-Field Modeling of Materials Interfaces and Nanostructures

Nanostructured materials offer unique properties for a wide range of energy applications. Among various known processing methods, liquid metal dealloying (LMD)—the selective dissolution of a base alloy element into a metallic melt—has emerged as a powerful technique to produce a new class of nano-/meso-scale open porous and bicontinuous composite structures with ultra-high interfacial area. LMD has recently been complemented by the advent of vapor phase dealloying (VPD), a novel technique that exploits the selective evaporation of one element from a parent alloy containing elements with very different vapor pressures, thereby enabling to fabricate open nanoporous structures of various elements from less-noble metals to inorganic elements regardless of their chemical activity without requiring high LMD temperatures and chemical etching. Together LMD and VPD have greatly expanded the scope of dealloying techniques, limited by traditional electrochemical means to noble metals, and boosted the design of new functional and structural materials that combine a wide variety of elements. Topologically-connected open porous structures with ultra-high surface area allow mass transport within the structure while preserving structural integrity, enabling them to serve as catalysts, fuel cells, supercapacitors, or high-capacity battery materials. Bicontinuous composite structures in turn can display high strength and high ductility or superior radiation-damage resistance due to the ultra-high interface area between interpenetrating solid phases. This research program makes use of state-of-the-art computational methods to understand at a basic level the self-organizing dealloying process with main focus on dealloying kinetics and interfacial pattern formation at the dealloying front controlling initial structure size, topology, and phase compositions. Phase-field simulation studies of LMD focus on solid solutions, line compounds, and intermetallic systems that can form ternary composites by nucleation and growth of a new phase. Studies of VPD employ phase-field modeling and a hybrid method combining a kinetic Monte Carlo (KMC) model of evaporation and surface diffusion with molecular dynamics for vapor-phase transport inside nanopores. Simulations explore mechanisms of interface- and diffusion-controlled dealloying kinetics, both observed in VPD but not fundamentally understood. In addition, our newly developed multi-physics phase-field approach of large-volume-change phase transformations is being used to model novel 3D anode geometries based on dealloyed nanoporous structures including novel sandwiched graphene/Si/silica for highrate Li ion battery. Those studies are aimed at elucidating geometric design principles that improve mechanical stability. We expect this research to enhance the capability to tailor nano-/mesoscale structures for a wide range of energy-related materials applications and to yield further advances in computational methodologies that benefit a broad materials research community.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2022

Fragmentation of faceted crystalline wires

The reduction of surface energy is a well-known driving force for the fragmentation of nanowires. Fragmentation can limit the service life of those small low-dimensional building blocks or can be controlled to induce beneficial shape changes. For isotropic surface energy, fragmentation is analogous to the classic Rayleigh-Plateau instability of liquid jets. However, commonly synthesized crystalline nanowires have strongly anisotropic surface energies and exhibit facets coinciding with cusps in the Wulff plot. Depending on growth orientation, different fragmentation behaviors have been seen in nanowires. Using phase-field simulations, we show that fragmentation of faceted nanowires with cubic crystal symmetry may occur by a finite-amplitude nonlinear instability, as opposed to a Rayleigh-Plateau-like linear instability, depending on nanowire growth orientation. Here, we carry out a weakly nonlinear analysis based on sharp interface theory to characterize the faceted nanowire shape corresponding to a surface-energy saddle point. The analysis predicts that the minimum amplitude of a periodic shape perturbation to trigger fragmentation increases with cusp strength but decreases inversely proportionally to the perturbation wavelength for long wavelengths, in good quantitative agreement with phase-field predictions. The results provide the theoretical foundation to predict nanowire stability as a function of length and surface energy in diverse applications.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Microstructural pattern formation during liquid metal dealloying: Phase-field simulations and theoretical analyses

In recent years, liquid metal dealloying has emerged as a promising material processing method to generate micro- and nanoscale bicontinuous or porous structures. Most previous studies focused on the experimental characterization of the dealloying process and on the properties of the dealloyed materials, leaving the theoretical study incomplete to fully understand the fundamental mechanisms of the liquid metal dealloying process. In this paper, we use theoretical models and phase-field simulations to clarify the kinetics and pattern formation during liquid metal dealloying. Our investigation starts from a theoretical analysis of the 1D dissolution of a binary precursor alloy, which reveals that the 1D dissolution process involves two regimes. In the first regime, due to the low solubility of one of the elements in the melt, it accumulates at the solid-liquid interface, which reduces the dissolution kinetics. In the second regime, the interface kinetics reaches a stationary regime where both elements of the precursor alloy dissolve into the melt. Previous works revealed that in the early dealloying stage, the dealloying front is destabilized by an interfacial spinodal decomposition, which triggers the formation of interconnected ligaments. Here we extend this line of work by proposing a linear stability analysis able to predict the initial length scale of the ligaments formed in the initial stage of the dealloying. Combining this analysis with the 1D dissolution model proposed here enables us to better understand the initial conditions (composition of the precursor alloy and the melt) leading to a planar dissolution without interface destabilization. Finally, we report a strong influence of solid-state diffusion on dealloying. Although the solid-state diffusivity is four to five orders of magnitude smaller than in the liquid phase, it is found to affect both dissolution kinetics and ligament morphologies.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Vulnerable window of yield strength for swelling-driven fracture of phase-transforming battery materials

Abstract Despite numerous experimental and theoretical investigations of the mechanical behavior of high-capacity Si and Ge Li-ion battery anodes, our basic understanding of swelling-driven fracture in these materials remains limited. Existing theoretical studies have provided insights into elasto-plastic deformations caused by large volume change phase transformations, but have not modeled fracture explicitly beyond Griffith’s criterion. Here, we use a multi-physics phase-field approach to model self-consistently anisotropic phase transformation, elasto-plastic deformation, and crack initiation and propagation during lithiation of Si nanopillars. Our computational results reveal that fracture occurs within a “vulnerable window” inside the two-dimensional parameter space of yield strength and fracture energy and highlight the importance of taking into account the surface localization of plastic deformation to accurately predict the magnitude of tensile stresses at the onset of fracture. They further demonstrate how the increased robustness of hollow nanopillars can be understood as a direct effect of anode geometry on the size of this vulnerable window. Those insights provide an improved theoretical basis for designing next-generation mechanically stable phase-transforming battery materials undergoing large volume changes.

25 ENERGY STORAGE↗
DOE OSTI2020

Configurational stability of a crack propagating in a material with mode-dependent fracture energy – Part II: Drift of fracture facets in mixed-mode I+II+III

In earlier papers (Leblond et al., 2011; 2019), we presented linear stability analyses of the coplanar propagation of a crack loaded in mixed-mode I+III, based on a “double” propagation criterion combining (Griffith, 1920)’s energetic condition and (Goldstein and Salganik, 1974)’s principle of local symmetry. The difference between the two papers was that in the more recent one, the local value of the critical energy-release-rate was no longer considered as a constant, but heuristically allowed to depend upon the ratio of the local mode III to mode I stress intensity factors. This led to a much improved, qualitatively acceptable agreement of theory and experiments, for the “threshold” value of the ratio of the unperturbed mode III to mode I stress intensity factors, above which coplanar propagation becomes unstable. In this paper, the analysis is extended to the case where a small additional mode II loading component is present in the initially planar configuration of the crack, generating a small, general kink of this crack from the moment it is applied. The main new effect resulting from presence of such a loading component is that the instability modes present above the threshold must drift along the crack front during its propagation. This prediction may be useful for future theoretical interpretations of a number of experiments where such a drifting motion was indeed observed.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Universality and Stability Phase Diagram of Two-Dimensional Brittle Fracture

The two-dimensional oscillatory crack instability, experimentally observed in a class of brittle materials under strongly dynamic conditions, has been recently reproduced by a nonlinear phase-field fracture theory. Here, we highlight the universal character of this instability by showing that it is present in materials exhibiting widely different near crack tip elastic nonlinearity, and by demonstrating that the oscillations wavelength follows a universal master curve in terms of dissipation-related and nonlinear elastic intrinsic length scales. Moreover, we show that upon increasing the driving force for fracture, a high-velocity tip-splitting instability emerges, as experimentally demonstrated. The analysis culminates in a comprehensive stability phase diagram of two-dimensional brittle fracture, whose salient properties and topology are independent of the form of near tip nonlinearity.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Pattern formation during electrochemical and liquid metal dealloying

Dealloying has evolved from a problematic corrosion process to a versatile tool for scalable fabrication of nanostructured metals. While the original, and majority of, work in the area has focused on electrochemical dealloying, a powerful variation of dealloying—liquid metal dealloying—has recently gained popularity. This process relies on a melt to carry out selective dissolution, replacing the traditional electrolyte solution. While electrolytes and molten metals are both suitable dealloying media, they can lead to very different morphologies. In this article, we compare and contrast what is known about the microscale physics and chemistry controlling microstructural evolution in electrochemical and liquid metal dealloying. Here, we conclude that the core phenomenology of porosity evolution—a competition between dissolution and interface diffusion—is similar in both dealloying processes, but that the relative magnitudes of these two processes control interfacial pattern formation.

McCue, Ian↗
DOE OSTI2017

Instability in dynamic fracture and the failure of the classical theory of cracks

Cracks, the major vehicle for material failure, undergo a micro-branching instability at ~40% of their sonic limiting velocity in three dimensions. Recent experiments showed that in thin systems cracks accelerate to nearly their limiting velocity without micro-branching, until undergoing an oscillatory instability. Despite their fundamental importance, these dynamic instabilities are not explained by the classical theory of cracks, which is based on linear elasticity and an extraneous local symmetry criterion to predict crack paths. We develop a two-dimensional theory for predicting arbitrary paths of ultrahigh-speed cracks, which incorporates elastic nonlinearity without extraneous criteria. We show that cracks undergo an oscillatory instability controlled by small-scale, near crack-tip, elastic nonlinearity. This instability occurs above an ultrahigh critical velocity and features an intrinsic wavelength proportional to the ratio of the fracture energy to the elastic modulus, in quantitative agreement with experiments. As a result, this ratio emerges as a fundamental scaling length assumed to play no role in the classical theory of cracks, but shown here to strongly influence crack dynamics.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Elastically mediated interactions between grain boundaries and precipitates in two-phase coherent solids

In this work, we investigate analytically and numerically the interaction between grain boundaries and second phase precipitates in two-phase coherent solids in the presence of misfit strain. Our numerical study uses amplitude equations that describe the interaction of composition and stress [R. Spatschek and A. Karma, Phys. Rev. B 81, 214201 (2010)] and free-energies corresponding to two-dimensional hexagonal and three-dimensional BCC crystal structures that exhibit isotropic and anisotropic elastic properties, respectively. We consider two experimentally motivated geometries where (i) a lamellar precipitate nucleates along a planar grain boundary that is centered inside the precipitate, and (ii) a circular precipitate nucleates inside a grain at a finite distance to an initially planar grain boundary. For the first geometry, we find that the grain boundary becomes morphologically unstable due to the combination of long-range elastic interaction between the grain boundary and compositional domain boundaries, and shear-coupled grain boundary motion. We characterize this instability analytically by extending the linear stability analysis carried out recently [P.-A. Geslin, Y.-C. Xu, and A. Karma, Phys. Rev. Lett. 114, 105501 (2015)] to the more general case of elastic anisotropy. The analysis predicts that elastic anisotropy hinders but does not suppress the instability. Simulations also reveal that, in a well-developed nonlinear regime, this instability can lead to the break-up of low-angle grain boundaries when the misfit strain exceeds a threshold that depends on the grain boundary misorientation. For the second geometry, simulations show that the elastic interaction between an initially planar grain boundary and an adjacent circular precipitate causes the precipitate to migrate to and anchor at the grain boundary.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Grain growth competition during thin-sample directional solidification of dendritic microstructures: A phase-field study

Here we present the results of a comprehensive phase-field study of columnar grain growth competition in bi-crystalline samples in two dimensions (2D) and in three dimensions (3D) for small sample thicknesses allowing a single row of dendrites to form. We focus on the selection of grain boundary (GB) orientation during directional solidification in the steady-state dendritic regime, and study its dependence upon the orientation of two competing grains. In 2D, we map the entire orientation range for both grains, performing several simulations for each configuration to account for the stochasticity of GB orientation selection and to assess the average GB behavior. We find that GB orientation selection depends strongly on whether the primary dendrite growth directions have lateral components (i.e. components perpendicular to the axis of the temperature gradient) that point in the same or opposite directions in the two grains. Here, we identify a range of grain orientations in which grain selection follows the classical description of Walton and Chalmers. We also identify conditions that favor unusual overgrowth of favorably-oriented dendrites at a converging GB. We propose a simple analytical description that reproduces the average GB orientation selection from 2D simulations within statistical fluctuations of a few degrees. In 3D, we find a similar GB orientation selection as in 2D when secondary branches grow in planes parallel and perpendicular to the sample walls. Remarkably, quasi-2D behavior is also observed even when those perpendicular sidebranching planes are rotated by a finite azimuthal angle about the primary dendrite growth axis as long as the absolute values of those azimuthal angles are equal in both grains. In contrast, when the absolute values of those azimuthal angles differ markedly, we find that unusual overgrowth events at a converging GB are promoted by a high azimuthal angle in the least-favorably-oriented grain. We also find that diverging GBs can be strongly affected by those azimuthal angles, while converging GBs exhibit a weak dependence on those angles. Finally, for diverging GBs, GB orientation is also strongly affected by the relative signs of the lateral components of the primary dendrite growth directions in both grains.

36 MATERIALS SCIENCE↗