Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Stability boundary”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

A weighted Shifted Boundary Method for free surface flow problems

The Shifted Boundary Method (SBM) belongs to the class of unfitted (or immersed, or embedded) finite element methods and was recently introduced for the Poisson, linear advection/diffusion, Stokes, Navier-Stokes, acoustics, and shallow-water equations. By reformulating the original boundary value problem over a surrogate (approximate) computational domain, the SBM avoids integration over cut cells and the associated problematic issues regarding numerical stability and matrix conditioning. Accuracy is maintained by modifying the original boundary conditions using Taylor expansions. Hence the name of the method, that shifts the location and values of the boundary conditions. In this article, we extend the SBM to the simulation of incompressible Navier-Stokes flows with moving free-surfaces, by appropriately weighting its variational form with the elemental volume fraction of active fluid. This approach prevents spurious pressure oscillations in time, which would otherwise be produced if the total active fluid volume were to change abruptly over a time step. In fact, the proposed weighted SBM method induces small mass (i.e., volume) conservation errors, which converge quadratically in the case of piecewise-linear finite element interpolations, as the grid is refined. Finally, we present an extensive set of two- and three-dimensional tests to demonstrate the robustness and accuracy of the method.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Phase transitions in HfO2 probed by first-principles computations

Ever since ferroelectricity was discovered in HfO2, the question of its origin remains controversial. Here, we probe this question using a combination of Landau theory of phase transitions and first-principles computations. In such an approach, the energy landscape associated with the phase transition between cubic and different experimentally demonstrated phases of HfO2 (tetragonal, monoclinic, orthorhombic Pbca, orthorhombic Pnma, and orthorhombic Pca21) is explored using density functional theory calculations. Computations revealed that stabilization of all but orthorhombic Pbca phase is driven by a single unstable zone-boundary antipolar mode X2−. When coupled with zone-center modes (Γ1+ and Γ3+), it stabilizes the tetragonal phase. Coupling with four additional modes (Γ5+, X3−, X5−, X5+) results in the monoclinic phase, which is the ground state of the material. If, however, Γ5+ mode is replaced with Γ4− mode, orthorhombic polar phase Pca21 is stabilized. The application of this framework to examine the effect of electric field on the ferroelectric phase of hafnia reveals that the field of 5 MV/cm is capable of stabilizing ferroelectric phase over the monoclinic one at 0 K.

Physics↗
DOE OSTI2021

Stability Functions in the Stable Surface Layer Derived from the Mellor–Yamada–Nakanishi–Niino (MYNN) Scheme

It is desirable that a surface layer scheme in an atmospheric numerical model is consistent with an atmospheric boundary layer scheme incorporated in the same model. In this study, stability functions based on Monin–Obukhov similarity theory for momentum and heat, φ m and φ h , in the stable surface layer are derived from the Mellor–Yamada–Nakanishi–Niino (MYNN) scheme. The resulting stability functions are approximated by φ m = 1 + 4.8z/L and φ h = 0.74 + 6.0z/L, which can be analytically integrated with respect to height z to obtain momentum and heat fluxes, where L is the Obukhov length. The fluxes thus, obtained are compared with those acquired from stability functions in four previous studies: they are nearly the same as those from two of them, and show better agreement with observational data of the Surface Heat Budget of the Arctic Ocean experiment (SHEBA) over sea ice than those from the other two studies. Detailed comparisons of the results of the MYNN scheme with the SHEBA data suggest that significant variations of the fluxes observed during “winter” when the ice was covered with dry snow may have been caused by those of the surface roughness around the observational site. The stability functions obtained from the MYNN scheme predict that the bulk and flux Richardson numbers approach critical values of 0.26 and 0.21, respectively, in the limit of z/L → ∞. These critical values result from the Kolmogorov hypothesis applied to the turbulent dissipation in the MYNN scheme and are considered to correspond to a transition from Kolmogorov to non-Kolmogorov turbulence.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Anharmonicity and Phase Diagram of Magnesium Oxide in the Megabar Regime

With density functional molecular dynamics simulations, we computed the phase diagram of MgO from 50 to 2000 GPa up to 20000 K. Via thermodynamic integration (TDI), we derive the Gibbs free energies of the B1, B2, and liquid phases and determine their phase boundaries. With TDI and a pseudo-quasi-harmonic approach, we show that anharmonic effects are important and stabilize the B1 phase in particular. As a result, the B1-B2 transition boundary in the pressure-temperature plane exhibits a steep slope. We predict the B1-B2-liquid triple point to occur at approximately T= 10000 K and P = 370 GPa, which is higher in pressure than was inferred with quasi-harmonic methods alone. We predict the principal shock Hugoniot curve to enter the B2 phase stability domain but only over a very small range of parameters. Here, this may render it difficult to observe this phase with shock experiments because of kinetic effects.

3-dimensional systems↗
DOE OSTI2025

Development of an interatomic potential for the Ta–Li system

A new interatomic potential for the Ta–Li system is introduced to facilitate the study of phase stability, mechanical properties and non-equilibrium dynamics after Li implantation in Ta. Here, this potential is based on a generalization of the embedded atom method (GEAM) and includes contributions from embedding energy, explicit two- and three-body interactions, and nonlocal many-body interaction terms. The parameters of the potential are optimized using energies and atomic forces for a wide range of configurations obtained from ab initio density functional theory (DFT) calculations. The potential is rigorously validated across a range of physical properties, including elastic constants, equations of state, phonon dispersion curves, point defect properties, and melting temperatures for different compositions. Although our potential is trained on a small dataset, its accuracy is comparable to that of available machine learning potentials for Li and Ta. Our simulations show that at temperatures below 500 K, Li atoms in Ta–Li alloys form clusters separated by Ta-rich domains, and we find no evidence of ordered phase formation. For Li concentrations below a few percent, Li atoms preferentially segregate to surfaces and grain boundaries. However, in alloys containing more than ~10% Li, the accumulation of Li in symmetric-tilt grain boundaries can lead to one of the following effects: formation of amorphous-like regions, changes in grain boundary structural units, or lateral movement of the grain boundary.

GEAM potential↗
DOE OSTI2026

The Spatial and Temporal Variability of the Clear Convective Boundary Layer at the ARM SGP Supersite

The convective boundary layer (CBL), also known as the mixing layer, constitutes the critical lower segment of the atmosphere that significantly influences daily human activities. The growing demand for precise weather forecasts is driven by the requirements of agriculture, transportation, and routine societal functions. Here, to enhance understanding of the CBL, this study investigates the spatiotemporal variability in the CBL and its controlling factors using four-year Doppler lidar, surface flux, and profiling measurements at five ARM Southern Great Plains sites within a 100 km radius. This investigation utilizes data collected exclusively under clear-sky conditions or scattered low-cloud conditions. Results reveal significant spatial differences in CBL evolutions. Daily mixing layer heights (MLHs) vary up to 1 km (30% of the mean) in late afternoon. There is a clear east–west contrast: western sites (C1, E32, E37) exhibit higher summer MLH (1.9–2.1 km) and vertical velocity variances (1.0–1.2 m 2 s −2 ) than eastern sites (1.6–1.8 km), reversing in winter. Temporally, the MLH peaks at 70% of the sunrise–sunset interval, the lagging heat flux (HF) peaks at 50%; and the seasonal MLH maxima lag the HF by approximately one month, influenced by nighttime PBL (planetary boundary layer) properties. The HF and lower tropospheric stability are the main factors of influence for the CBL, but site-specific dependencies highlight the critical roles of local factors, underscoring the need for including them in CBL modeling.

ARM SGP site↗
DOE OSTI2026

Operational stability of mixed Sn–Pb perovskite solar cells: Mechanisms, mitigation strategies, and perspectives

Mixed Sn–Pb perovskites provide the 1.2–1.3 eV narrow-bandgap absorber needed for high efficiency of all perovskite tandems, but their deployment is limited by operational instability under light. This review synthesizes mechanistic origins and recent mitigation strategies for mixed Sn–Pb perovskite solar cells. The oxidation of Sn 2+ to Sn 4+ , often driven by iodine formation under light and bias, is the primary failure pathway; it creates Sn vacancies, self-doping, and nonradiative loss. Surface/grain-boundary defects, halide migration, reactive oxygen species, and interfacial redox at charge-transport layers, along with hole accumulation from poor band alignment, further accelerate Sn–Pb perovskite degradation. Here we survey recent stability advances across additive chemistry, surface and grain-boundary passivation, buried-interface redesign with modified or alternative hole transport layers, and solvent systems that preserve Sn 2+ and correct Sn–Pb speciation for scalable coating. Together, these recent advances have enabled devices to retain 80%–90% output for hundreds to over a thousand hours. Lastly, we provide our perspectives on further improving the operational stability of Sn–Pb perovskite and solar cells.

14 SOLAR ENERGY↗
DOE OSTI2025

Identifying synoptic controls on boundary layer thermodynamic and cloud properties in a regional forecast model

Although most of our understanding of boundary layer cloudiness is based on idealized, subtropical, barotropic marine environments, boundary layer clouds exist across a range of conditions. In this study, we use the Naval Research Laboratory's Coupled Ocean/Atmosphere Mesoscale Prediction System (COAMPS) and an automated cold-front-relative analysis framework to explore the boundary layer structure associated with low clouds across a transect through the cold front of a midlatitude synoptic cyclone. The model credibly captures boundary layer structure in line with conceptual models of midlatitude cyclones and ground-based observations at Graciosa Island in the Azores. The warm sector is conditionally unstable, with clouds that are too shallow and with too little liquid water, compared to cloud property retrievals from satellite and surface-based instruments. The cold-frontal region exhibits convection associated with weak stability and ascent. Northwest of the cold front, the boundary layer is well mixed, deeper, and capped by a strong inversion maintained by large-scale subsidence. Simulated clouds in frontal and post-frontal regions are mostly too thick, with too much liquid water and too little cloud base drizzle. The post-frontal clouds are associated with grid-scale updrafts, which appear to be the model's attempt to represent mesoscale organization of cellular convection typically observed in the cold sector of midlatitude cyclones. The deep, well-mixed post-frontal boundary layers and cloudiness are maintained by strong surface fluxes, as in cold air outbreaks. Our analysis framework serves as a unique approach to model verification, and our results offer insights into differences in boundary layer cloud behavior between subtropical and synoptic cold-sector regimes.

Eissner, Jordan M. [Univ. of Kansas, Lawrence, KS ↗
DOE OSTI2023

Mechanistic origins of corrosion protection of aluminum alloys by graphene/polyetherimide nanocomposite coatings

Light-weighting vehicular components through adoption of light-metal structural alloys holds promise for reducing the fuel consumption of internal combustion engine vehicles and increasing the range of battery electric vehicles. However, the alloyed microstructure and surface precipitates of aluminum alloys render these materials susceptible to corrosion under modest excursions from neutral pH. Traditional chromium-based anodic passivation layers are subject to increasingly stringent environmental regulations, whereas options for sacrificial cathodic films are sparse for electropositive metals. While hybrid nanocomposite coatings have shown initial promise, mechanistic underpinnings remain poorly understood. Here, a fully imidized polyetherimide (PEI) resin is utilized as the continuous phase with inclusion of unfunctionalized exfoliated graphite (UFG). A comprehensive investigation of the mechanisms of corrosion protection reveals key fundamental design principles underpinning corrosion inhibition. First, strong interfacial adhesion, which for PEI is facilitated by binding of imide carbonyl moieties to Lewis acidic sites on Al surfaces. Second, the miscibility of ion-impervious nanoscopic UFG fillers and stabilization of a substantial interphase region at UFG/PEI boundaries that result in minimizing the free volume at the filler/polymer interface. Finally, extended tortuosity of ion diffusion pathways imbued by the below-percolation-threshold 2D fillers. These three design principles help govern and modulate ion transport from electrolyte/coating interfaces to the coating/metal interface and are crucial for the extended preservation of barrier properties. The results suggest an approach to systematically activate multiple modes of corrosion inhibition through rational design of hybrid nanocomposite coatings across hard-to-abate sectors where light metal alloys are likely to play an increasingly prominent role.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Evaluation of idealized large-eddy simulations performed with the Weather Research and Forecasting model using turbulence measurements from a 250 m meteorological mast

Abstract. We investigate the ability of the Weather Research and Forecasting model to perform large-eddy simulation of canonical flows. This is achieved through comparison of the simulation outputs with measurements from sonic anemometers on a 250 m meteorological mast located at Østerild, in northern Denmark. Østerild is on a flat and rough area, and for the predominant wind directions, the atmospheric flow can be considered to be close to homogeneous. The idealized simulated flows aim at representing atmospheric boundary layer turbulence under unstable, neutral, and stable stability conditions at the surface, which are statistically significant conditions observed at Østerild. We found that the resolved fields from the simulations appear to have the characteristics of the three stability regimes. Vertical profiles of observed mean wind speeds and direction are well reproduced by the simulations, with the largest differences under near-neutral conditions, where the effect of the subgrid-scale model is evident on the vertical wind shear close to the surface. Vertical profiles of observed eddy fluxes are also well reproduced by the simulations, with the largest differences for the three velocity component variances under stable stability conditions, although nearly always within the observed variability. With regards to turbulent kinetic energy, we find good agreement between observations and simulations at all vertical levels. Simulated and observed velocity spectra match very well and show very similar behavior with height and with atmospheric stability within the low-frequency interval; at the effective resolution, the simulated spectra show the typical drop-off of finite differences. Our findings demonstrate that these idealized simulations reproduce the characteristics of atmospheric stability regimes often observed at a high turbulent and flat site within a direction sector, where the air flows over nearly homogeneous land.

17 WIND ENERGY↗
DOE OSTI2023

Thermal and Radiation Stability in Nanocrystalline Cu

Nanocrystalline metals have presented intriguing possibilities for use in radiation environments due to their high grain boundary volume, serving as enhanced irradiation-induced defect sinks. Their promise has been lessened due to the propensity for nanocrystalline metals to suffer deleterious grain growth from combinations of irradiation and/or elevated homologous temperature. While approaches for stabilizing such materials against grain growth are the subject of current research, there is still a lack of central knowledge on the irradiation–grain boundary interactions in pure metals despite many studies on the same. Due to the breadth of available reports, we have critically reviewed studies on irradiation and thermal stability in pure, nanocrystalline copper (Cu) as a model FCC material, and on a few dilute Cu-based alloys. Our study has shown that, viewed collectively, there are large differences in interpretation of irradiation–grain boundary interactions, primarily due to a wide range of irradiation environments and variability in materials processing. We discuss the sources of these differences and analyses herein. Then, with the goal of gaining a more overarching mechanistic understanding of grain size stability in pure materials under irradiation, we provide several key recommendations for making meaningful evaluations across materials with different processing and under variable irradiation conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Boundary Layer Climatology at ARM Southern Great Plains

Operational since 1992, the Atmospheric Radiation Measurement (ARM) southern great plains (SGP) site at Oklahoma, USA has become a reference research site for meteorological studies. Due to an open data policy the ARM data are used by researchers all over the world. In this report, we review the long-term climatology of the atmospheric boundary layer, SGP instrumentations, the site and some site-specific atmospheric conditions which potentially effect wind turbines within the region. As the atmospheric boundary layer is bounded and influenced by the surface, observations of surface radiation components and heat fluxes are crucial in understanding land-atmosphere interactions. The entrainment of air, updrafts, downdrafts and boundary layer height characteristics is needed for understanding the structure and growth of the atmospheric boundary layer. Therefore, measurements from both surface in-situ and remote sensing observations at SGP provide an overall climatology and their interactions from surface up to the boundary layer. Measurements from a 60 m meteorological tower, surface flux stations, disdrometers, soil temperature and moisture flux plates, coherent Doppler lidar, Raman lidar, radiosondes, and satellite data at SGP central facility were analyzed. All the measurements were generally within a few square kilometers of each other at the central facility. This report focuses on data from January 2010 to June 2020 at SGP central facility. The various sections describe the ARM SGP site and surrounding wind turbines; in-situ and remote sensing instrumentation used in the report; provides mathematical equations to analyze fluxes, turbulence and other boundary layer parameters; a climatological analysis of surface winds, fluxes and thermodynamic parameters for several years; an analysis of observed winds in the framework of Monin-Obukhov Similarity Theory; an analysis of the boundary layer winds and direction from a Doppler lidar; multi-year turbulence estimates through the boundary layer from a Doppler lidar; atmospheric boundary layer water vapor and relative humidity profiles from Raman lidar; cloud base height and boundary layer height from multiple sensors and satellite data; and finally site specific atmospheric conditions, such as nocturnal low-level jets. Diurnal, seasonal and yearly variations of surface, sub-surface and boundary layer quantities, such as wind speed, direction, temperature, atmospheric stability, soil temperature, and various atmospheric fluxes at SGP showed distinct trends useful for focused modeling studies. The applicability of surface similarity theory on ARM SGP data is also evaluated, which showed northerly flows are aligned with MO theory estimates compared to southerly flows. Boundary layer winds and direction profiles for several years from a Doppler lidar shows a consistent presence of a nocturnal low-level jet and predominant southerly wind directions through the boundary layer at SGP. The inter-annual variability at SGP is low (<3.5%), with a mean annual wind speed of approximately 7 m s -1 at 100 m above ground level. Boundary layer turbulence and moisture transport from Doppler and Raman lidars are evaluated, which provides evidence of increased water vapor mass flux into the great plains during nocturnal low-level jets. The moisture flux from nocturnal low-level jets is observed to be maximum during summer periods. A novel machine learning algorithm is implemented to accurately estimate the planetary boundary layer height, providing further insights into growth and destruction of the convective boundary layer height during various seasons and land-atmosphere conditions. The high frequency of low-level clouds during winter, spring and fall seasons is validated using the multi-sensor array and satellite estimates of cloud top height. Satellite vegetative fraction data provides insight into seasonal surface roughness and vegetation variability around SGP site.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Plasma Wall Interaction with 3-D Plasma Boundaries

The interaction of the edge plasma and the material surfaces is one of the most critical challenges on the path to harness fusion power as new, fundamental energy source. This challenge typically combines the thrust to reach high density, low temperature (detached) plasmas in front of the divertor target plates as well as understanding the plasma material interaction (PMI) in particular in this regime. The combination of both research thrusts represents an extraordinarily challenging subject encompassing spatial scales spanning nanometers to meters in all states of matter and across a broad energy range. Modeling capabilities, which help to interpret data from nowadays experiments and enable extrapolation to future devices are urgently required. This is in particular true for toroidal magnetic confinement devices with three-dimensional (3D) plasma boundaries. Such plasma boundary geometries occur in tokamaks, when small amplitude magnetic perturbations are used to stabilize the unruly edge plasma or in stellarators, that are inherently 3D plasma confinement devices. In this project, the impact of 3D plasma boundaries on the plasma material interaction (PMI) was assessed. This work focused on plasma boundary conditions, in which high-density conditions at the material surfaces yield mitigation of the otherwise immense heat and particle loads that these materials would see. These so-called high recycling and eventually detached plasma regimes are of great interest for future reactor operation. In the project, key features that are unique to 3D boundaries were explored in comparison to canonically assumed axisymmetric plasma edge situations in tokamaks. In particular, the relevance of the 3D boundary situation in the extrapolation to the plasma boundary solution at ITER, the next step fusion energy experiment under construction as a multi-national, world-wide large-science experiment in southern France, has been explored. The EMC3-EIRENE plasma edge fluid and kinetic neutral transport code has been advanced to cope with the challenging and unprecedented conditions in the ITER boundary plasma including 3D fields that are planned to be used to suppress harmful edge instabilities, the so-called edge localized modes. This is a vital integration challenge for ITER and the results from this grant have provide a leading capability for this assessment. It was shown that the detachment process in a 3D edge solution for ITER follows the recycling regimes that are known from axisymmetric solutions, but that multiple plasma exhaust channels connected to the material surfaces are established which feature individual recycling characteristics. Because these channels touch the material surfaces in the divertor in a 3D geometry, the compatibility with the plasma material interaction (PMI), including erosion and impurity generation has been found to be an important part of the integration challenge. To address this, the fully 3D plasma material interaction code ERO2 has been adapted to these ITER specific geometries and a homogeneous mixing model was implemented, that allows to consider the mixing of Be and Was used at ITER in the PMI modeling. This model enhancement has been used to study non-local migration of Be in the JET ITER like wall configuration and it has been shown that with this model such complex migration processes in ITER relevant plasma shapes and with ITER relevant plasma boundary conditions can be addressed. The combined modeling approach using EMC3-EIRENE as a plasma boundary transport code and the ERO2 specialized PMI model will be an asset for the continued preparations of ITER operation as well as for Fusion Pilot Plant efforts that have emerged in the U.S. during the evolution of this grant. The predictive capability of this numerical tool has been validated at the DIII-D US national fusion facility. Here, dedicated plasma edge diagnostics were implemented to measure the impurity household around a 3D edge plasma during ELM suppression by 3D fields. Dedicated experiments with local material probes using these diagnostics and the state-of-the-art suite of boundary measurements at DIII-D have shown that the 3D perturbation of the plasma edge that is excreted by such 3D control fields yield a perturbation of the plasma boundary flux structure and hence also of the resulting PMI. The 3D boundary plasma is composed out of helical magnetic flux channels that intersect the divertor targets at an angle relative to the main guiding field, i.e., the toroidal magnetic field component of the tokamak. A similar effect has been measured as well on limiter surfaces during the startup campaign at the new stellarator experiment Wendelstein 7-X. These experiments ad initial analysis with the ERO plasma material interaction model, suggested that the place of erosion for a given particle from the surface and its re-deposition can be different in such 3D field geometries yielding potentially a significant level of net-erosion. This is not the case for axisymmetric solutions, where it was shown in the past that the eroded particles are effectively re-deposited into gaps produced by erosion at the same position and hence the net-erosion levels are small. For ITER, the quest to suppress the ELMs and at the same time maintain the integrity of the divertor is an issue, which these fundamental findings will help to resolve. The coupling of this work to the extrapolation in the ITER program has been addressed by both the PI and the lead numerical scientist being ITER Science Fellows in the duration of the contract and forward. A second focus in the exploration of 3D boundary effects on tokamaks and stellarators has been set on the measurement of helium exhaust features with such 3D fields. This is important because He represents the ash of the fusion process and needs to be exhausted. It was shown that 3D field application compatible with suppression of ELMs yields an increase of the helium exhaust performance. The ratio of the effective helium confinement time over the energy confinement time was reduced by almost 50% which demonstrated that the impact of helium accumulation in the plasma core with respect to the confinement of energy to sustain the fusion reaction is significantly improved with such 3D control fields. It was shown that this is the case for tokamaks as well as stellarators. At the Large helical Device in Japan, a similar enhancement of the helium exhaust features when small amplitude additional 3D fields were applied was measured. This is an important additional function of 3D field application and its impact on ITER is presently being studied in combination with investigations of helium exhaust in 3D field geometries of stellarator devices.

3D plasma edge transport↗
DOE OSTI2025

Grain boundary zirconia-modified garnet solid-state electrolyte

Here, we report a method for promoting electrochemical stability in garnet Li 6.4 La 3 Zr 1.4 Ta 0.6 O 12 solid-state electrolyte based on a composite two-phase oxide–oxide microstructure. Grain boundary precipitation of the controlled distribution of amorphous zirconium oxide microparticles is achieved through the addition of reactive tantalum carbide. During ambient-atmosphere sintering, the carbide decomposes through an in situ reaction, the ‘extra’ Ta substituting for Zr within the Li 6.4 La 3 Zr 1.4 Ta 0.6 O 12 lattice. Density functional theory (DFT) calculations identify a thermodynamically favourable reaction path and show how substituting Ta 5+ at Zr 4+ sites affects the crystal structure as well as bulk ionic and electronic conductivities. Quantitative stereology highlights that zirconia also acts as a sintering aid, reducing compact porosity. Cryogenic focused-ion-beam scanning electron microscopy and fractography analysis of cycled solid-state electrolytes illustrates that near-universally observed intergranular Li-metal dendrite propagation is suppressed by the two-phase microstructure, favouring transgranular dendrites instead. Importantly, DFT demonstrates that compared with the Li 6.4 La 3 Zr 1.4 Ta 0.6 O 12 surface, the zirconium oxide surface per se is less electronically conductive and does not trap excess electrons to reduce Li ions. This is a key reason for the substantial improvement in the electrochemical properties over the single-phase baseline.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Disentangling Electronic Transport and Hysteresis at Individual Grain Boundaries in Hybrid Perovskites via Automated Scanning Probe Microscopy

Underlying the rapidly increasing photovoltaic efficiency and stability of metal halide perovskites (MHPs) is the advancement in the understanding of the microstructure of polycrystalline MHP thin film. Over the past decade, intense efforts have been aimed at understanding the effect of microstructures on MHP properties, including chemical heterogeneity, strain disorder, phase impurity, etc. It has been found that grain and grain boundary (GB) are tightly related to lots of microscale and nanoscale behavior in MHP thin films. Atomic force microscopy (AFM) is widely used to observe grain and boundary structures in topography and subsequently to study the correlative surface potential and conductivity of these structures. For now, most AFM measurements have been performed in imaging mode to study the static behavior; in contrast, AFM spectroscopy mode allows us to investigate the dynamic behavior of materials, e.g., conductivity under sweeping voltage. However, a major limitation of AFM spectroscopy measurements is that they require manual operation by human operators, and as such only limited data can be obtained, hindering systematic investigations of these microstructures. Here, we designed a workflow combining the conductive AFM measurement with a machine learning (ML) algorithm to systematically investigate grain boundaries in MHPs. The trained ML model can extract GBs locations from the topography image, and the workflow drives the AFM probe to each GB location to perform a current–voltage (IV) curve automatically. Then, we are able to have IV curves at all GB locations, allowing us to systematically understand the property of GBs. Using this method, we discovered that the GB junction points are less conductive, potentially more photoactive, and can play critical roles in MHP stability, while most previous works only focused on the difference between GB and grains.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Impact of Vibrational Nonequilibrium on the Supersonic Mode in Hypersonic Boundary Layers

The supersonic mode in hypersonic boundary layers has recently been shown to achieve a larger peak growth rate than the traditional second mode in particular circumstances relevant to free flight. In this environment, high levels of aerodynamic heating are encountered such that incorporation of chemical and vibrational nonequilibrium effects can become critical. The impact of thermal nonequilibrium (TNE) on the supersonic mode in particular has not been thoroughly and systematically investigated. This work uses thermochemical nonequilibrium (TCNE) direct numerical simulation (DNS) and TCNE linear stability theory (LST) to obtain a more complete investigation of the supersonic mode using both nonequilibrium and frozen gas models. The simulation is Mach 5 flow over a 1-mm-nose-radius axisymmetric cone 1 m in length, representative of experimental conditions in ground test facilities. Both LST and DNS results indicate that TNE effects are destabilizing compared with frozen flow and can lead to an N-factor difference of approximately +1 at the end of the cone. LST results also indicate that TNE effects are stabilizing to the second and supersonic modes in comparison with thermal equilibrium. The supersonic mode appeared to be slightly more affected by TNE than the second mode. In conclusion, collectively, the findings in this work suggest that even in ground test facilities where TCNE effects are relatively small compared with free flight, incorporating nonequilibrium models in both the mean flow and stability solvers can have a practical effect on hypersonic boundary-layer transition prediction analyses.

42 ENGINEERING↗
DOE OSTI2023

High-temperature active oxidation of nanocrystalline silicon-carbide: A reactive force-field molecular dynamics study

Flexible woven SiC ceramics are prone to accelerated fiber embrittlement under high temperature oxidation in dynamic oxygen environments. The nanocrystalline structure of the constituent fibers impacts the reaction kinetics and phase transformations during active oxidation. However, fundamental understanding and quantification of grain boundary effects on oxidation behavior in nanocrystalline SiC remain elusive when temperatures exceed 1500 K. This study deploys large-scale molecular dynamics simulations with a reactive force-field to elucidate the complex roles of atomic oxygen reservoir conditions and grain size on oxidation kinetics and the nature of oxides produced in both monocrystalline and nanocrystalline 3C-SiC between 1100 K and 2000 K. The simulations with dynamically replenished oxygen provide good agreement with oxidation kinetics and activation energies for the monocrystalline Si(100) and C(100) orientations published in the available literature. This study reveals that, by contrast, nanocrystalline SiC samples exhibit two distinct oxidation kinetics with a transition point at 1500 K due to surface melting, which is supported by experimental evidence. The introduction of a grain-boundary network produces a two-fold decrease in oxidation activation energies compared to monocrystalline SiC below 1500 K. Above 1500 K, however, the activation energies rise substantially due to the formation of a liquid Si phase at the SiC/Si oxide interface. Further, it is shown that the stability of the interfacial liquid phase is promoted by incoherent grain boundaries in the crystalline SiC. These findings are important for the deployment of nanocrystalline SiC fibers in advanced thermal protection systems for high-temperature applications.

36 MATERIALS SCIENCE↗