SEARCH · Search NASA

Results for “Finite elements”

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.

SourceResearch connectionFrom yearThrough year

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

At least 235 records · Page 13

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Mesh-based super-resolution of fluid flows with multiscale graph neural networks

A graph neural network (GNN) approach is introduced in this work which enables mesh-based three-dimensional super-resolution of fluid flows. In this framework, the GNN is designed to operate not on the full mesh-based field at once, but on localized meshes of elements (or cells) directly. To facilitate mesh-based GNN representations in a manner similar to spectral (or finite) element discretizations, a baseline GNN layer (termed a message passing layer, which updates local node properties) is modified to account for synchronization of coincident graph nodes, rendering compatibility with commonly used element-based mesh connectivities. Furthermore, the architecture is multiscale in nature, and is comprised of a combination of coarse-scale and fine-scale message passing layer sequences (termed processors) separated by a graph unpooling layer. The coarse-scale processor embeds a query element (alongside a set number of neighboring coarse elements) into a single latent graph representation using coarse-scale synchronized message passing over the element neighborhood, and the fine-scale processor leverages additional message passing operations on this latent graph to correct for interpolation errors. Demonstration studies are performed using hexahedral mesh-based data from Taylor–Green Vortex and backward-facing step flow simulations at Reynolds numbers of 1600 and 3200. Through analysis of both global and local errors, the results ultimately show how the GNN is able to produce accurate super-resolved fields compared to targets in both coarse-scale and multiscale model configurations. Reconstruction errors for fixed architectures were found to increase in proportion to the Reynolds number. Geometry extrapolation studies on a separate cavity flow configuration show promising cross-mesh capabilities of the super-resolution strategy.

Backward-facing step↗

DOE OSTI2024

Processing MCNP Elemental Edit Outputs

The Monte Carlo N-Particle (MCNP) transport code version 6 (also known as MCNP6) has the capability for tracking particles on unstructured mesh (UM) geometry models embedded into constructive solid geometry (CSG) cells. A UM geometry is a collection of elements representing a solid geometry. The first step of MCNP UM modeling is using other software packages to create a finite element mesh representation of a solid 3D geometry. Computer-aided design (CAD) or computer-aided manufacturing (CAM) software is typically used to create a solid geometry model, which is later imported into mesh generation software to create a UM model. The MCNP UM feature was originally designed for models generated by the Abaqus/CAE software. The MCNP code version 6.0 and later can process UM models formatted as Abaqus input files. MCNP can process a UM model consisting of several different element types including linear tetrahedral or hexahedral elements and calculate quantities of interest such as flux and energy deposition at elements. An MCNP UM simulation provides high-fidelity elemental edit (i.e., tally) outputs, which can be further used in multiphysics calculations. The MCNP UM feature was used for multiphysics simulations where quantities of interest calculated by MCNP are used as inputs for heat transfer calculations in Abaqus. MCNP6.3 can produce two types of elemental edit output (EEOUT) file formats: ASCII and HDF5. An EEOUT file type must be requested on an EMBED card while output type (flux or energy deposition) must be requested on an EMBEE card. We wrote Python3 scripts to extract energy deposition values in an ASCII or HDF5 EEOUT file and compute a heat flux profile for an Abaqus heat transfer calculation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

Addendum to SAND2023-09604 Xyce lumped-element transmission line model verification to support Empire-Cable cable SGEMP analyses

This report supplements the Verification of Empire-Cable SAND report by expanding on the use of Xyce to simulate the coupling to a transmission line cable model. While Empire-Cable solves its governing equations on a high-order, finite-element mesh with an an implicit-in-time formulation, Xyce must use a first order graph for the circuit and explicit-in-time approach to be compatible with non-linear electrical device models. Thus, given the different solution methodologies in Xyce as compared to Empire-Cable, the convergence rates are expected to be different but the overall quality of the solution should be the same. The original four canonical problems studied in the Empire-Cable verification report are replicated here running in Xyce using transmission line modeling parameters from the verification report. Overall, agreement between the codes is excellent with Xyce’s convergence rates being limited mostly to first order due to the circuit network approximation of a transmission line being a first order approximation.

42 ENGINEERING↗

DOE OSTI2025

Blockage effects in the chemotaxis of diffusiophoretic particles

Transport mechanisms at the micro- and nano-scale play an essential role in regulating intracellular organization. Recent work indicates that directed motion of constituents inside cells can emerge through diffusiophoretic transport, in which colloidal particles move under the influence of chemical gradients. Here, we examine how blockers—passive or actively consuming—reshape those gradients and thereby influence the motion of diffusiophoretic particles. By combining analytical solutions with finite element simulations, we first show that a single blocker can distort a background gradient enough to create or eliminate stagnation points, significantly modifying particle transport. We then introduce a second, explicitly sized blocker at one of these stagnation points and measure how its finite radius alters the diffusiophoretic velocity field for a test particle. Even moderate changes in the second blockers size can cause noticeable shifts in the substrate distribution, highlighting the importance of accounting for explicit particle radii under crowded or consumption-driven conditions. Our findings underscore that subtle geometric variations—such as the radii and positions of two or more blockers—can profoundly affect diffusiophoretic motion, providing a more complete picture of how blocking and crowding phenomena shape intracellular transport.

Song, Zehao [Northwestern Univ., Evanston, IL (Uni↗

DOE OSTI2025

Half-closed discontinuous Galerkin discretisations

Here we introduce the concept of half-closed nodes for nodal discontinuous Galerkin (DG) discretisations. Unlike more commonly used closed nodes in DG, where on every element nodes are placed on all of its boundaries, half-closed nodes only require nodes to be placed on a subset of the element's boundaries. The effect of using different nodes on DG operator sparsity is studied and we find in particular for there to be no difference in the sparsity pattern of the Laplace operator whether closed or half-closed nodes are used. On quadrilateral/hexahedral elements we use the Gauss-Radau points as the half-closed nodes of choice, which we demonstrate is able to speed up DG operator assembly in addition to leverage previously known superconvergence results. We also discuss in this work some linear solver techniques commonly used for Finite Element or discontinuous Galerkin methods such as static condensation and block-based methods, and how they can be applied to half-closed DG discretisations.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Alloying effects on deformation induced microstructure evolution in copper

In this work, we investigated the effects of alloying elements on plastic deformation and microstructure evolution in polycrystalline copper (Cu) and Cu alloyed with 1 wt. % lead (Cu-1 %Pb). These materials were selected due to the size mismatch between Cu and Pb, with the latter forming precipitates at grain boundaries. Multi-modal characterization techniques, including neutron diffraction, electron backscatter diffraction (EBSD), and transmission electron microscopy (TEM), along with finite element simulations were employed to study the deformation behavior across multiple length scales. While both Cu and Cu-1%Pb exhibited similar macroscale response and final deformation textures, both dislocation line profile analysis and TEM revealed increased dislocation density in deformed Cu-1%Pb specimens. The presence of lead precipitates also significantly affected local plastic deformation during compression, with their influence diminishing with increasing strain. These results demonstrate the complex relationships between alloying elements, plastic deformation, microstructural evolution, and material behavior under load. The insights gained from this multi-scale and multi-technique approach contribute to the fundamental understanding of microstructural evolution in immiscible alloys and are valuable for tailoring the properties of structural materials for specific engineering applications.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Griffin: A MOOSE-based reactor physics application for multiphysics simulation of advanced nuclear reactors

Griffin is a Multiphysics Object-Oriented Simulation Environment (MOOSE) based reactor physics application for multiphysics simulations of advanced reactor designs jointly developed by Idaho National Laboratory and Argonne National Laboratory. This paper summarizes the motivation, significance, architecture, design, and features of Griffin. Griffin offers flexible and extensible features to address the challenges associated with advanced reactor designs. These features range from fundamental particle transport to specific reactor physics tasks. The features cover a wide range including on-the-fly and traditional two-step cross-section generation methods, steady-state and transient transport solvers suitable for both heterogeneous and homogeneous models, high-fidelity depletion where thousands of isotopes can be tracked and low-fidelity depletion characterized by burnup, etc. The most fundamental aspect that sets Griffin apart from other reactor analysis codes is that it is developed based on the MOOSE framework. A modular development approach is strongly enforced, with multiphysics being an essential element considered since the beginning of Griffin’s development. Griffin links various MOOSE physics modules and couples to other MOOSE-based applications and non-MOOSE-based applications for multiphyiscs simulations. Griffin includes three modules: ISOXML for preparing and managing multigroup cross sections, radiation transport for solving the neutron transport equation, and reactor analysis for user-oriented reactor physics analysis functionalities. Griffin uses various finite element methods for spatial discretization, multigroup approximation for energy discretization and discrete ordinates method, spherical harmonics expansion method, and diffusion approximation for streaming direction discretization to solve the neutron transport equation. Griffin’s flexibility is evidenced through Griffin’s various applications to fast reactor, high-temperature reactor, pebble bed reactor, molten salt reactor, and microreactor designs. Griffin development follows the software quality assurance procedure for MOOSE-based applications and with software requirements consistent with the ASME NQA-1 standard. Griffin has been adopted into the reactor analysis system for the U.S. NRC and is in use at U.S. companies, universities and national laboratories.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2026

C‐Cracking of Brittle‐Material Spheres From Eccentric Hertzian Contact

Contact of brittle-material spheres can occur while they are manufactured, handled, or used in operation as roller elements or in deformably processed composite material. Statistically, the directional contact between two spheres will nearly always be oblique or eccentric in an unconfined and mobile sphere population. Although analyses of oblique contact exist, the portrayal of its maximum (tensile) first principal stress ( S 1 ) field is lacking. This information is needed for improved judgment of the prospect of crack initiation. Given these, the S 1 due to eccentric contact and subsequent crack initiation was analyzed using finite element analysis (FEA) and corroborative demonstration of c-crack creation. The FEA shows that an asymmetric S 1 field is created about the contact patch and is caused by superimposed shear intrinsic to eccentric contact. A crack will initiate if the maximum S 1 exceeds the material's tensile strength, and its field will cause the crack to propagate and arrest to form a C-shaped crack, or c-crack. C-crack initiation trends with lower applied forces, with greater amount of contact eccentricity, higher coefficient of friction, and smaller spherical radius. Finally, these metrics deserve recognition, especially when c-cracked brittle-material spheres are used in a system whose functionality or reliability is potentially compromised by the c-cracks.

Hertzian contact↗

DOE OSTI2024

Observation and Modeling of Dynamic Fracture Behaviors of Battery Cell Under Impact Loading Using Enhanced Representative Volume Element Concept

The burgeoning electric automobile industry has increased interest in battery safety. Battery cells experience significant mechanical stress during operation, including the impact of accidents and vibrations from driving. The potential for thermal runaway reactions in battery cells raises safety concerns. Although numerous researchers have defined the dynamic behavior of battery cells and proposed numerical models to describe it, few studies have focused on the high-strain rate mechanical impact phase correlated with the onset of fracture. In this study, we describe the dynamic behavior of pouch battery cells and propose a modeling method to study their mechanical failure under impact situations. Impact tests are conducted at various velocities and heights. To overcome numerical issues commonly encountered under rapid deformation scenarios, a new finite element model is developed based on the representative volume element model. The proposed approach efficiently simulates continuous crack propagation and brittleness behavior during impact by permitting the individual behavior of the cell components. Therefore, engineers can reliably design safer electric vehicle battery cells by measuring the properties of the cell components.

ENERGY STORAGE↗

DOE OSTI2025

Characterization of W production during ICRF operations: experiments and modeling

For successfully heating plasma with waves in the ion cyclotron range of frequencies (ICRFs), mitigating impurity production is just as crucial as maximizing power coupling, especially in high-Z environments. ICRF can effectively deposit energy on ions, modify turbulence-driven transport, and enhance fusion reaction efficiency, but only when its power coupling has minimal impact on impurity production. To do so, one must rely on a toroidal array of at least three active elements excited with appropriate phasing and power ratio to reduce the currents induced on the antenna frame below levels critical for physical sputtering. In contrast to classic two-strap antennas, which are optimized for dipole phasing with equal power on both straps, three-strap antennas in ASDEX Upgrade (AUG)—but also four-strap antennas in JET, Alcator C-Mod, SPARC and ITER—offer the possibility to act also on the power ratio between the central and outer straps. With optimal settings, impurity production can be reduced substantially, making the ICRF compatible with the high-Z wall. This paper explores the characteristics of the AUG three-strap antennas in terms of impurity production, as well as the key role of plasma composition in this process. Numerical simulations were performed using SSWICH and Petra-M (finite element codes) to quantify impurity production and compare with experimental results. Energies of ions falling on antenna limiters (measured with probes) are well predicted by both codes. These tools are then used to further describe the source of the impurity, namely the gross erosion of tungsten from an ICRF antenna, for different plasma mixtures. Results are also compared to spectroscopy data. Ultimately, we show that deleterious effects of the ICRF on plasma surface interactions will be weaker in plasmas containing larger fractions of highly ionized heavier low-Z impurity, which is typically relevant for experiments relying on impurity seeding.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2024

Enhancing the cooling performance of thermocouples: a power-constrained topology optimization procedure

Abstract Heat pumping through thermoelectric devices has many advantages over traditional cooling. However, their current efficiency is a limiting factor in their implementation. In this paper, we approach the non-convex topology optimization of thermoelectrical elements for cooling applications through the method of moving asymptotes (MMA) to improve their cooling capabilities per watt usage. The optimization problem is defined for a given power budget, aiming for the minimum temperature with a known heat pumping need. The introduction of power as a constraint justifies the introduction of the voltage gradient across the thermocouple as a design variable to maintain the thermoelectrical device in its optimum power-to-heat extraction ratio. To better understand the convergence of this non-convex problem, we present a two-variable analytical thermoelectric optimization model. This example provides information on how to select the penalty parameters used to scale the three material coefficients involved in the problem to obtain lower objective values and better convergence using MMA. The analytical model shows the non-convexity of the problem and provides the recommendation to use penalization coefficients of the form $$p_k=p_{\sigma }>p_{\alpha }=1$$ p k = p σ > p α = 1 for the thermal conductivity, electrical conductivity, and Seebeck coefficients. We tested these penalization coefficients through optimizations of a model based on the 1MC10-031 commercial thermoelectric-cooler (TEC) using the finite element method (FEM). These penalization coefficients provided local minima without the need for volume constraints. With this procedure, we found designs that provided temperatures close to 10 degrees lower using 60% less semiconductor material volume compared to the initial design.

Gutiérrez, G. Reales↗

DOE OSTI2025

Tailoring composition and deformation modes at the microstructural level for next generation low-cost high-strength austenitic stainless steels

The objective of this project is to enable deliberate development of cost-effective, hydrogen resistant alloys by establishing detailed relationships specific to the effects of alloy composition, short-range order (SRO), and microsegregation in the presence of hydrogen on the transition between homogeneous deformation and localized plasticity in shear bands. In collaboration with the International Institute for Carbon-Neutral Energy Research, I2CNER, at Kyushu University in Japan, we conceptualized, designed, and manufactured four austenitic alloys that maintain corrosion resistance and ensure lower cost relative to baseline commercial alloys. The mechanical properties and deformation modes of the novel alloys (KU alloys) were assessed in the presence of hydrogen (H). Correlations between composition and performance revealed that two of the KU alloys are suitable replacements for 316 steel, while another is a viable replacement for 304 steel at room temperature. We found that, in the presence of other austenite stabilizing elements namely Mn and N, replacing Ni with Cu does not lead to martensite formation as has been previously reported.1–3 Furthermore, we found that the addition of Cu leads to an earlier onset of multiple slip resulting in an relative earlier onset of a higher work hardening rate (WHR). Greater understanding of the relationships between alloy composition and SRO required the development of a novel advanced electron diffraction methodology to characterize SRO in complex FCC alloys. This innovative approach, which combines fluctuation and correlation analyses of diffuse-scattering signals, successfully differentiated between SRO and long-range ordering (LRO). Further investigations into annealed austenitic stainless steels could provide insights into manipulating SRO and its effects on material properties. Atomistic simulations provided understanding of SRO behavior that was difficult to capture experimentally. This project created the first spin cluster expansion model that is able to capture and describe SRO effects in Fe-Ni-Cr FCC alloys, accounting for the non-negligible effects of magnetism. An automated computational workflow was established to provide reliable predictions of SRO in Fe-Ni-Cr austenitic alloys, both with and without the presence of H atoms. Analysis of the propensity for SRO in Fe-Ni-Cr alloys revealed that H tends to cluster with specific, well-defined SRO domains. The computational framework is general purpose and can be extended to realistic stainless steels across diverse composition ranges. With confidence that SRO is possible in austenitic stainless steels, we developed a discrete dislocation finite element code to understand the interaction of dislocations with SRO in the presence of H. By incorporating H effects on the dislocation emission and SRO stress field we show that the critical stress for the dislocation pileup to breakthrough the SRO domain decreases in the presence of H, which directly contributes localized deformation at the macroscale. Through the simulation of a uniaxial tension test, we demonstrated that H-induced weakening of SRO stress field and H-enhanced dislocation emission can lead to the onset of shear localization at lower macroscopic strains. As a whole, this project identified three novel alloys that show improvements in performance and cost efficiency for H-facing applications by studying correlations between alloy chemistry and deformation behavior. We also made significant advancements to experimental and computational methodologies necessary to study the chemistry and distribution of SRO across a range of alloys, which in turn allowed us to demonstrate how deformation mechanisms change due to the contributions of SRO in austenitic alloys in the presence of H. The combined advancements in fundamental understanding with novel alloy development in this project has increased the viability of next generation H-technologies for the broader public through accessible low-cost alloys and accelerated development towards future H-infrastructure.

08 HYDROGEN↗

DOE OSTI2026

Mechanical Behavior of Neutron Irradiated Refractory Multi-Principal Element Alloys Processed via Spark Plasma Sintering

The search for advanced materials capable of withstanding the extreme conditions of Generation IV reactors is a critical area in materials science research. These reactors operate under severe environments, including high temperatures, corrosion, stress, and irradiation damage. Consequently, there is a need for innovative alloy systems to ensure the reliability and longevity of proposed Generation IV reactor components. Refractory multi-principal-element alloys (RMPEA) have emerged as a promising candidate due to their exceptional properties. These alloys, characterized by their composition of multiple principal elements in near-equiatomic ratios, exhibit superior resistance to irradiation damage, reduced void swelling, enhanced microstructural stability, and minimal irradiation-induced hardening. While initial studies on RMPEAs have shown promising results, most research has been limited to thin films, nanocrystalline microstructures, and ion irradiation, which do not accurately represent the behavior of bulk materials. To address this gap, our research focused on the neutron irradiation of bulk RMPEAs. We aim to conduct comprehensive post-irradiation examinations (PIE) of RMPEAs irradiated at the Advanced Test Reactor at Idaho National Laboratory. The RMPEAs were synthesized using spark plasma sintering (SPS) with mechanically alloyed metallurgical powder. The RMPEA specimens are a MoNbTi alloy system with additions of -Zr, and -ZrV. Furthermore, PIE consisted of mechanical testing and advanced materials characterization. The mechanical testing consisted of sub-sized tensile testing, micro- and nano- indentation. Microstructural characterization included scanning electron microscopy and transmission electron microscopy. Mechanical testing coupled with advanced microscopy techniques provides insight into phase morphology and its effects on the mechanical properties of the RMPEA specimens. The results indicate that both pristine and irradiated RMPEA specimens exhibited brittle behavior during tensile testing, which can be attributed to their heterogeneous microstructure. The SPS manufacturing process did not include any post treatment, which resulted in a heterogeneous microstructure. Energy-dispersive X-ray spectroscopy revealed the presence of intermetallic such as laves phases within the microstructure. Specifically, Ti-rich precipitates were observed in the MoNbTi specimen, while Mo-rich precipitates were found in the MoNbTiZrV specimen. Nano-hardness testing of pristine samples showed that the laves phases exhibited higher hardness values compared to the matrix phase, suggesting that precipitate hardening is likely the dominant hardening mechanism in these specimens. The results from this work will be used to build a finite element model to predict mechanical behavior of future MPEA compositions. Thus, enabling for a streamlined approach to developing novel MPEAs for the nuclear industry.

36 - MATERIALS SCIENCE↗

DOE OSTI2024

Multilaminate Energy Storage Films from Entropy‐Driven Self‐Assembled Supramolecular Nanocomposites

Abstract Composite materials comprising polymers and inorganic nanoparticles (NPs) are promising for energy storage applications, though challenges in controlling NP dispersion often result in performance bottlenecks. Realizing nanocomposites with controlled NP locations and distributions within polymer microdomains is highly desirable for improving energy storage capabilities but is a persistent challenge, impeding the in‐depth understanding of the structure–performance relationship. In this study, a facile entropy‐driven self‐assembly approach is employed to fabricate block copolymer‐based supramolecular nanocomposite films with highly ordered lamellar structures, which are then used in electrostatic film capacitors. The oriented interfacial barriers and well‐distributed inorganic NPs within the self‐assembled multilaminate nanocomposites effectively suppress leakage current and mitigate the risk of breakdown, showing superior dielectric strength compared to their disordered counterparts. Consequently, the lamellar nanocomposite films with optimized composition exhibit high energy efficiency (>90% at 650 MV m −1 ), along with remarkable energy density and power density. Moreover, finite element simulations and statistical modeling have provided theoretical insights into the impact of the lamellar structure on electrical conduction, electric field distribution, and electrical tree propagation. This work marks a significant advancement in the design of organic–inorganic hybrids for energy storage, establishing a well‐defined correlation between microstructure and performance.

Li, He↗

DOE OSTI2025

Architected Liquid Crystal Elastomer Lattices with Programmable Energy Absorption

Architected LCE lattices are fabricated with flow-induced alignment via direct ink writing and systematically characterized their shape morphing, stiffness, and energy absorption behavior across strain rates spanning six orders of magnitude from 10 −3 to 10 3 s −1 . It is shown that architected liquid crystal elastomer (LCE) lattices exhibit superior energy absorption compared to their non-mesogenic (silicone) counterparts. Importantly, the LCE-to-silicone energy absorption ratios are up to 18-fold higher at the highest strain rate tested. A finite element model that captures their shape-morphing response is developed, which exhibits excellent agreement with the experimental observations. In conclusion, the work opens new avenues for designing and fabricating LCE lattices with programmable alignment, shape morphing, and mechanics.

active lattices↗

DOE OSTI2024

Size Effect of Negative Capacitance State and Subthreshold Swing in Van der Waals Ferrielectric Field–Effect Transistors

Analytical calculations corroborated by the finite element modeling show that thin films of Van der Waals ferrielectrics covered by a 2D-semiconductor are promising candidates for the controllable reduction of the dielectric layer capacitance due to the negative capacitance (NC) effect emerging in the thin films. The NC state is conditioned by energy-degenerated poly-domain states of the ferrielectric polarization induced in the films under incomplete screening conditions in the presence of a dielectric layer. Calculations performed for the FET-type heterostructure “ferrielectric CuInP 2 S 6 film—2D-MoS 2 single-layer—SiO 2 dielectric layer” reveal the pronounced size effect of the multilayer capacitance. Derived analytical expressions for the electric polarization and multilayer capacitance allow to predict the thickness range of the dielectric layer and ferrielectric film for which the NC effect is the most pronounced in various Van der Waals ferrielectrics, and the corresponding subthreshold swing becomes much less than the Boltzmann's limit. Obtained results can be useful for the size and temperature control of the NC effect in the steep-slope ferrielectric FETs.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Proximity Ferroelectricity in Compositionally Graded Structures

Proximity ferroelectricity is a novel paradigm for inducing ferroelectricity in a non-ferroelectric polar material, such as AlN or ZnO that are typically unswitchable with an external field below their dielectric breakdown field. When placed in direct contact with a thin switchable ferroelectric layer (such as Al 1-x Sc x N or Zn 1-x Mg x O), they become a practically switchable ferroelectric. Using the thermodynamic Landau-Ginzburg-Devonshire theory, in this work, we perform the finite element modeling of the polarization switching in the compositionally graded AlN-Al 1-x Sc x N, ZnO-Zn 1-x Mg x O, and MgO-Zn 1-x Mg x O structures sandwiched in both a parallel-plate capacitor geometry as well as in a sharp probe-planar electrode geometry. We reveal that the compositionally graded structure allows the simultaneous switching of spontaneous polarization in the whole system by a coercive field significantly lower than the electric breakdown field of unswitchable polar materials. The physical mechanism is the depolarization electric field determined by the gradient of chemical composition “x”. The field lowers the steepness of the switching barrier in the otherwise unswitchable parts of the compositionally graded AlN-Al 1-x Sc x N and ZnO-Zn 1-x Mg x O structures. In the MgO-like regions of the compositionally graded MgO-Zn 1-x Mg x O structure, a shallow double-well free energy potential emerges. Proximity ferroelectric switching of the compositionally graded structures placed in the probe-electrode geometry occurs due to nanodomain formation under the tip. We predict that a gradient of chemical composition “x” significantly lowers effective coercive fields of the compositionally graded AlN-Al 1-x Sc x N and ZnO-Zn 1-x Mg x O structures compared to the coercive fields of the corresponding multilayers with a uniform chemical composition in each layer. A tip-induced switching further lowers the coercive field, enabling control of ferroelectric domains in otherwise unswitchable compositionally graded structures, which can provide nanoscale domain control for memory, actuation, sensing, and optical applications.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Transient Chemo‐Mechanical Model of Lithium Plating Impacted by External Pressure

Abstract Lithium anodes show great promise in commercial applications, but are hindered by lithium plating and dendrite growth which cause safety concerns during long‐term cell operation. Stack pressure is experimentally observed to improve cell lifetime; however, the relationship between stress and lithium deposition has remained difficult to elucidate. In this work, a transient, 3D, finite‐element model of the evolution of a lithium anode due to stripping and plating is developed. The evolution of a microscale protrusion on the anode surface is tracked over one charge‐discharge cycle with respect to stack pressure and lithium yield strength. Lithium plastic deformation, nonconformal anode‐separator contact, and separator porosity effects are accounted for. Over the course of several hours of stripping/plating, the anode surface evolves to a similar morphology under pressure regardless of the initial conditions due to lithium plastic deformation and hardening. The rate of this evolution highly depends on the applied pressure and assumed lithium yield strength.

25 ENERGY STORAGE↗