Engineering topics
Publications and source records attributed to Heo, Tae Wook.
In this quarterly report, we Computed the driving force for the water splitting reaction as a function of reaction progress and continued Development of the numerical method for water splitting reactions of H 2 production on the particle surface within 3D porous microstructures.
Metal hydrides are serious contenders for materials-based hydrogen storage to overcome constraints associated with compressed or liquefied H 2 . Their ultimate performance is usually evaluated using intrinsic material properties without considering a systems design perspective. An illustrative case with startling implications is (LiNH 2 +2LiH). Using models that simulate the storage system and associated fuel cell of a light-duty vehicle (LDV), the performance of the bulk hydrides is compared with a nanoscaled version in porous carbon (PC), (LiNH 2 +2LiH)@(6-nm PC). Using experimental material properties, the simulations show that (LiNH 2 +2LiH)@(6-nm PC) counterintuitively has higher usable gravimetric and volumetric capacities than the bulk counterpart on a system basis despite having lower capacities on a materials-only basis. Nanoscaling increases the thermal conductivity and lowers the desorption enthalpy, which consequently increases heat management efficiency. In a simulated drive cycle for fuel cell-powered LDV, the fuel cell is inoperable using bulk (LiNH 2 +2LiH) as the storage material but completes the drive cycle using the nanoscale material. Further, these results challenge the notion that nanoscaling incurs mass and volume penalties. Instead, the synergistic nanoporous host-hydride interaction can favorably modulate chemical and heat transfer properties. Moreover, a co-design approach considering application-specific tradeoffs is essential to accurately assess a material's potential for real-world hydrogen storage.
The effective diffusivity of ionic species in multiphase materials is critical for the design and function of composite materials for electrochemical energy storage. In practice, effective diffusivity depends sensitively not only on the intrinsic diffusivities of constituting materials but also on their topological arrangement; nevertheless, these coupled contributions are oversimplified in most analytical models. Here, we combine atomistically informed mesoscale modeling and machine learning (ML) analysis to unravel how such features affect effective diffusivity in two-phase composites. Using the Li 7 La 3 Zr 2 O 12 -LiCoO 2 composite solid-state battery cathode as a model system, we compute effective diffusivity for 600 distinct dense polycrystalline microstructures with different topological configurations of grains, grain boundaries, and heterointerfaces. We verify that in addition to atomic-scale variabilities, microstructural feature diversity can significantly impact effective transport properties. Across the ensemble of test microstructures, this often results in bimodal distributions of effective diffusivity that encompass two qualitatively distinct operating mechanisms, which we identify via flux analysis. An ML approach reveals that the most critical determining factors for effective diffusivity are the connectivity of bulk phases and their heterointerfaces. The role of ionic mobility at the heterointerfaces is also discussed. These insights highlight the combined importance of microstructure and interface engineering in tuning the transport properties of ionic species in composite materials. In conclusion, our framework can also be extended for understanding generic microstructure-property relationships in other complex multiphase materials.
Fast charging and stable lithium-ion batteries (LIBs) require prompt Li ion transport between electrodes. 3D porous separators can potentially improve ion transport at both the microscale and macroscale levels. Additive manufacturing (AM) is an ideal tool to produce and study 3D porous separators of LIBs, although little work has been reported in this area. Here, we demonstrate porous structural control in layer-by-layer printed hexanediol diacrylate (HDDA) separators using a projection micro stereolithography (PμSL) AM technology and present their corresponding battery performance in lithium nickel manganese cobalt oxide (NMC)/graphite full cells. The composition of photocurable resin determines the porous morphology (e.g., bicontinuous and droplet-like), while the settings of printing layer thickness and number of exposures can be used to adjust the porosity (e.g., from 27% to 63%). Optimal resin and print conditions are determined, and the resulting separator exhibits superior Li ion conductivity and therefore higher rate performance and longer cycling life. Compared with a Celgard 2325 separator, the optimized porous HDDA separator increases the 10C-rate capacity by 34% and decreases the capacity degradation rate at the 1C-rate by 4.9%. In conclusion, this work paves the way for future AM of 3D separators for high-energy, high-power, and long-cycling-life LIBs.
Grain boundaries (GBs) profoundly influence the properties and performance of materials, emphasizing the importance of understanding the GB structure and phase behavior. As recent computational studies have demonstrated the existence of multiple GB phases associated with varying the atomic density at the interface, we introduce a validated, open-source GRand canonical Interface Predictor (GRIP) tool that automates high-throughput, grand canonical optimization of GB structures. While previous studies of GB phases have almost exclusively focused on cubic systems, we demonstrate the utility of GRIP in an application to hexagonal close-packed titanium. We perform a systematic high-throughput exploration of tilt GBs in titanium and discover previously unreported structures and phase transitions. In low-angle boundaries, we demonstrate a coupling between point defect absorption and the change in the GB dislocation network topology due to GB phase transformations, which has important implications for the accommodation of radiation-induced defects.
The architectural design of electrodes offers new opportunities for next-generation electrochemical energy storage devices (EESDs) by increasing surface area, thickness, and active materials mass loading while maintaining good ion diffusion through optimized electrode tortuosity. However, conventional thick electrodes increase ion diffusion length and cause larger ion concentration gradients, limiting reaction kinetics. We demonstrate a strategy for building interpenetrated structures that shortens ion diffusion length and reduces ion concentration inhomogeneity. This free-standing device structure also avoids short-circuiting without needing a separator. The feature size and number of interpenetrated units can be adjusted during printing to balance surface area and ion diffusion. Starting with a 3D-printed interpenetrated polymer substrate, we metallize it to make it conductive. This substrate has two individually addressable electrodes, allowing selective electrodeposition of energy storage materials. Using a Zn//MnO 2 battery as a model system, the interpenetrated device outperforms conventional separate electrode configurations, improving volumetric energy density by 221% and exhibiting a higher capacity retention rate of 49% compared to 35% at temperatures from 20 to 0 °C. Our study introduces a new EESD architecture applicable to Li-ion, Na-ion batteries, supercapacitors, etc.
Naturally occurring coatings on aluminum metal, such as its oxide or hydroxide, serve to protect the material from corrosion. Understanding the conditions under which these coatings mechanically fail is therefore expected to be an important aspect of predictive models for aluminum component lifetimes. To this end, we develop and apply a molecular dynamics (MD) modeling framework for conducting tension tests that is capable of isolating factors governing the mechanical strength as a function of coating chemistry, defect morphology, and variables associated with the loading path. We consider two representative materials, including γ-Al 2 O 3 and γ-Al(OH) 3 (i.e., oxide and hydroxide), both of which form readily as aluminum surface coatings. Our results indicate that defects have a significant bearing on the strength of aluminum oxide, with grain boundaries serving to reduce the strain at failure from ε zz = 0.300 to 0.219, relative to perfect single crystal. Our simulations also predict that porosity lowers the elastic stiffness and yield strength of the oxide. Relative to perfect crystal, we find porosity factors of 5%, 10% and 20% decrease the yield stress by 26%, 36% and 53%, respectively. MD predicts that perfect hydroxide and oxide single crystal have respective strains at failure of 0.08 and 0.31 under tensile uniaxial strain loading, and that the corresponding yield stresses are respectively 1.6 and 11.1 GPa. These data indicate that the hydroxide is substantially more susceptible to mechanical failure than the oxide. Our results, coupled with literature findings that indicate hot and humid conditions favor formation of hydroxide and defective oxide coatings, indicate the potential for a complicated dependence of aluminum corrosion susceptibility and stress corrosion cracking on aging history.
Explore the source record for details and available documents.
Here we established a comprehensive modeling approach for investigating the microstructure-aware effective thermal conductivity ($κ_{eff}$) for porous microstructures containing solid particles and gaseous pores. Our approach combines the mesoscale computational modeling framework and the semi-analytical method, allowing for efficient prediction of $κ_{eff}$ for realistic porous microstructures, while considering complicated microstructural thermal conduction pathways effectively in the prediction. We used the diffuse-interface mesoscale computational model to generate extensive simulated $κ_{eff}$ data for realistic digital representations of microstructures with wide ranges of porosity ($f_p$), thermal conductivity of the gas phase ($κ_g$), and thermal conductivity of the solid phase ($κ_s$). From the simulated data, we identified two property variation regimes for $κ_{eff}$: (1) a slow $κ_{eff}$ increase for $κ_s ~ κ_g$; and (2) a faster $κ_{eff}$ increase for $κ_s \gg κ_g$. To capture the key features of the relationship between the microstructure and $κ_{eff}$, we derived a semi-analytical model by introducing structure and intensification factors. The two new factors incorporate the calibrated effective contribution of the solid volume with $κ_s$ and additional interfacial effects into the prediction of $κ_{eff}$, respectively, allowing for consideration of parallel, serial, and interfacial conduction mechanisms effectively. Using the selected simulation data, we quantified key model parameters within the semi-analytical model and verified that the parameterized model exhibits excellent agreement with simulated $κ_{eff}$ for the entire range of the parameter space.
Bare aluminum metal surfaces are highly reactive, which leads to the spontaneous formation of a protective oxide surface layer. Because many subsequent corrosive processes are mediated by water, the structure and dynamics of water at the oxide interface are anticipated to influence corrosion kinetics. Here, using molecular dynamics simulations with a reactive force field, we model the behavior of aqueous aluminum metal ions in water adsorbed onto aluminum oxide surfaces across a range of ion concentrations and water film thicknesses corresponding to increasing relative humidity. We find that the structure and diffusivity of both the water and the metal ions depend strongly on the humidity of the environment and the relative height within the adsorbed water film. Aqueous aluminum ion diffusion rates in water films corresponding to a typical indoor relative humidity of 30% are found to be more than 2 orders of magnitude slower than self-diffusion of water in the bulk limit. Connections between metal ion diffusivity and corrosion reaction kinetics are assessed parametrically with a reductionist model based on a 1D continuum reaction–diffusion equation. Our results highlight the importance of incorporating the properties specific to interfacial water in predictive models of aluminum corrosion.
Amorphous sulfides are among the most promising candidates for solid electrolytes (SEs) owing to their excellent deformability and acceptable Li-ion conductivity (σion) at room temperature. However, the complex atomic structure of these amorphous materials without long-range ordering results in a lack of structural understanding and difficulties in tuning material properties. Here, in this study, we performed structural mapping of glassy sulfide SEs composed of various anion clusters using a combinatorial atomic level analysis of synchrotron X-ray-based pair distribution function (PDF) and reverse Monte-Carlo (RMC) methods, demonstrating the potential to tune σ ion in glass SEs. First, we prepared a new glassy sulfide with a mixed anion framework of two halogens (Br and I) as atomic anions, with PS 4 3- molecular anions in a Li 2 S–P 2 S 5 -based glass SE. At a specific Br content, [(Li 2 S) 0.658 (LiI 0.9 LiBr 0.1 ) 0.342 ] 0.825 [P 2 S 5 ] 0.175 recorded a σ ion of 2.27 mS cm -1 , the highest value for any glass sulfides reported to date; however, its elastic modulus was still suppressed to 14.48 GPa. PDF and RMC calculations successfully provided structural mapping of anion clusters, including two halogens. Molecular dynamics simulations of each composition confirmed that flexible coordination caused by the rattling of small polarizable Br ions in the mixed halogens of glassy SEs contributed to the superior σ ion . Our results may provide new insights into the design of superior glassy SEs that play key roles in all-solid-state batteries requiring fast Li-ion conduction and high deformability.
Hydrogen incorporation in native surface oxides of metal alloys often controls the onset of metal hydriding, with implications for materials corrosion and hydrogen storage. A key representative example is titania, which forms as a passivating layer on a variety of titanium alloys for structural and functional applications. These oxides tend to be structurally diverse, featuring polymorphic phases, grain boundaries, and amorphous regions that generate a disparate set of unique local environments for hydrogen. Here, we introduce a workflow that can efficiently and accurately navigate this complexity. First, a machine learning force field, trained on ab initio molecular dynamics simulations, was used to generate amorphous configurations. Density functional theory calculations were then performed on these structures to identify local oxygen environments, which were compared against experimental observations. Second, to classify subtle differences across the disordered configuration space, we employ a graph-based sampling procedure. Finally, local hydrogen binding energies and hopping kinetics are computed using exhaustive density functional theory calculations on representative configurations. Here, we leverage this methodology to show that hydrogen binding energetics are described by local oxygen coordination, which in turn is affected by stoichiometry, and form the basis of hopping kinetics and diffusion. Together these results imply that hydrogen incorporation and transport in TiO x can be tailored through compositional engineering, with implications for improving the performance and durability of titanium-derived alloys in hydrogen environments.
Abstract The garnet‐type phase Li 7 La 3 Zr 2 O 12 (LLZO) attracts significant attention as an oxide solid electrolyte to enable safe and robust solid‐state batteries (SSBs) with potentially high energy density. However, while significant progress has been made in demonstrating compatibility with Li metal, integrating LLZO into composite cathodes remains a challenge. The current perspective focuses on the critical issues that need to be addressed to achieve the ultimate goal of an all‐solid‐state LLZO‐based battery that delivers safety, durability, and pack‐level performance characteristics that are unobtainable with state‐of‐the‐art Li‐ion batteries. This perspective complements existing reviews of solid/solid interfaces with more emphasis on understanding numerous homo‐ and heteroionic interfaces in a pure oxide‐based SSB and the various phenomena that accompany the evolution of the chemical, electrochemical, structural, morphological, and mechanical properties of those interfaces during processing and operation. Finally, the insights gained from a comprehensive literature survey of LLZO–cathode interfaces are used to guide efforts for the development of LLZO‐based SSBs.
Despite the promise of TiFe-based alloys as low-cost solid-state hydrogen storage materials with mild operating conditions and reasonable hydrogen capacity, their initial hydrogenation process is difficult, hindering broad utilization. The effect of alloying element on the initial hydrogenation kinetics of TiFe alloys, TiFe 0.9 M 0.1 (M = V, Cr, Fe, Co and Ni), was evaluated by analyzing changes to the passivating surface oxide layer that inhibits hydrogen permeation, as well as the ease of initial-stage hydrogen absorption into the underlying matrix. X-ray photoelectron spectroscopy and atom probe tomography revealed key variations in surface oxide compositions and thinning of the passivating oxide layer compared to pure TiFe, which suggests suppressed oxide growth by alloying-induced elemental redistribution. At the same time, density functional theory calculations predicted exothermic formation of hydride nuclei when alloying with V or Cr, as well as a reduced nucleation barrier when alloying with Co or Ni. Overall, these results are consistent with the observed experimental trend of the activation kinetics. We propose that improvements in activation kinetics of TiFe with alloying arises from the combined effect of reduced passivating oxide thickness and easier hydride nucleation, offering a starting point for alloy design strategies towards further improvement.
Industrial and public interest in hydrogen technologies has risen strongly recently, as hydrogen is the ideal means for medium to long term energy storage, transport and usage in combination with renewable and green energy supply. In a future energy system, the production, storage and usage of green hydrogen is a key technology. Hydrogen is and will in future be even more used for industrial production processes as a reduction agent or for the production of synthetic hydrocarbons, especially in the chemical industry and in refineries. Under certain conditions material based systems for hydrogen storage and compression offer advantages over the classical systems based on gaseous or liquid hydrogen. This includes in particular lower maintenance costs, higher reliability and safety. Hydrogen storage is possible at pressures and temperatures much closer to ambient conditions. Hydrogen compression is possible without any moving parts and only by using waste heat. In this paper, we summarize the newest developments of hydrogen carriers for storage and compression and in addition, give an overview of the different research activities in this field.
We report FeTi-based hydrides have recently re-attracted attention as stationary hydrogen storage materials due to favorable reversibility, good sorption kinetics and relatively low costs compared to alternative intermetallic hydrides. Employing the OpenCalphad software, the thermodynamics of the (FeTi) 1-x H x (0 ≤ x ≤ 1) system were assessed as a key basis for modeling hydrogenation of FeTi-based alloys. New thermodynamic data were acquired from our experimental pressure-composition-isotherm (PCI) curves, as well as first-principles calculations utilizing density functional theory (DFT). The thermodynamic phase models were carefully selected based on critical analysis of literature information and ab-initio investigations. Key thermodynamic properties such as dissociation pressure, formation enthalpies and phase diagrams were calculated in good agreement to our performed experiments and literature-reported data. This work provides an initial perspective, which can be extended to account for higher-order thermodynamic assessments and subsequently enables the design of novel FeTi-based hydrides. In addition, the assessed thermodynamic data can serve as key inputs for kinetic models and hydride microstructure simulations.
Engineering thermophysical properties of metal hydrides nanocomposites is crucial for effective thermal management during hydrogen storage reactions; however, the effect of microstructure on thermal transport mechanisms is still unclear. Here, we employed an integrated experiment-modeling approach to investigate microstructural factors that determine the effective thermal conductivity of individual reduced graphene oxide-magnesium (rGO/Mg) nanocomposites and their packed bed. Here the effective thermal conductivity of the rGO/Mg nanocomposite packed bed was measured by using guarded hot-plate method under various atmospheric conditions (i.e., vacuum, Ar and He). A microstructure-aware mesoscopic modeling revealed that anisotropy of the effective thermal conductivity of individual rGO/Mg nanocomposites plays an important role in determining the effective thermal conductivity of their packed bed. The validated mesoscopic model also disclosed a nontrivial interplay between the intrinsic rGO properties and the extrinsic composite structural features. Finally, quantitative sensitivity analysis based on the modeling framework is used to provide practical engineering guidance for controlling the thermal transport within nanocomposite packed beds.