Search NASA⌕ Search

SEARCH · Search NASA

Results for “failure modes and effects analysis”

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.

1,680 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Quantifying System Strength From Grid-Forming Resources Using Frequency Scan Approach: Preprint

Current industry practices for quantifying the system strength contribution from grid-forming (GFM) resources to ensure stability of power systems dominated by inverter-based resources (IBRs) are primarily based on iterative electromagnetic transient (EMT) time-domain simulation studies. While feasible, these approaches are resource-intensive, lack scalability and intuition, and might not evaluate the system strength contribution over the entire frequency range of interest. This paper introduces a novel, frequency-domain approach to quantify system strength support provided by a GFM resource using frequency scans. The proposed method uses transfer functions from the grid voltage magnitude (V) and phase (?), respectively, to the reactive (Q) and active power (P) output of a GFM resource for quantifying its contribution to system strength. These transfer functions provide a direct measure of the ability of a GFM resource to behave as a stiff voltage source behind a reactance over a specified frequency range, enabling robust quantification of its system strength contribution. The key innovation of this work is the development of a frequency domain system strength metric called the dynamic short-circuit ratio (dSCR) that is suitable for IBR-dominated power systems and is directly related with the familiar short circuit ratio (SCR) metric. The new metric, dSCR, enables the assessment of system strength contributions from both synchronous machines and converter-based GFM resources using a unified benchmark, which is not possible with the traditional SCR metric. The paper also demonstrates how impedance scans could identify if an unstable condition observed during weak grid conditions is a result of the lack active or reactive power support or both. By leveraging the proposed frequency-domain dSCR metric for quantifying system strength contribution from GFM IBRs, the paper demonstrates targeted mitigation strategies for weak grid instabilities without resorting to repeated, time-consuming time-domain simulations. The result is a scalable and efficient approach to remediate stability challenges in power systems with high shares of IBRs and accelerating the integration of GFM technologies for system strength support in power systems.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2026

EXERGETIC: De-Risking Next-Generation Resilient Geothermal Hybrids via At-Scale Evaluation Using Virtual Emulation Digital Twin Environment for Efficient Operation

The DOE-GTO-funded project, award number 5.1.2.12, entitled "EXERGETIC - De-risking Next Generation Resilient Geothermal Hybrids via at-Scale Evaluation Using a Virtual Emulation Digital Twin Environment for Efficient Operation," advances the solution to these challenges by developing and validating a geothermal co-emulation environment implemented at the National Laboratory of the Rockies (NLR)'s Advanced Research on Integrated Energy Systems (ARIES) platform. This framework enables the de-risking of next-generation geothermal and geothermal hybrid systems through high-fidelity modeling, real-time digital emulation, advanced control strategies, and techno-economic assessment. The project focused on geothermal hybrid configurations that integrate geothermal power plants with concentrated solar power and underground thermal energy storage, enabling enhanced efficiency, flexibility, and grid support capabilities. The main goal of this project was the development of a geothermal digital co-emulation environment to demonstrate the technical and economic value of geothermal hybrid systems and their contribution to grid stability and flexibility. The EXERGETIC framework combined physics-based models, controls, and real assets at ARIES, including digital real-time simulators (DRTS), a 20-MW-scale controllable grid interface (CGI), and a 2-MW conventional generator. Detailed transient models were developed for the key subsystems of a hybrid geothermal plant, including parabolic trough solar collectors, reservoir thermal energy storage (RTES), and a binary Organic Rankine Cycle (ORC) power plant. The ORC model explicitly captured thermal inertia and off-design operation and integrated control strategies to dynamically respond to electric load profiles. The models were validated against published experimental and numerical studies, demonstrating strong agreement and confirming the accuracy and robustness of the modeling approach. The resulting digital twin represents geothermal-solar-storage systems at multiple scales (1 MW to 100 MW) and enables realistic emulation of grid-connected operation. The control architecture allows the geothermal resource to provide stable baseload generation, while solar and stored thermal energy supply flexible, dispatchable support during periods of high demand or variable grid conditions. A key contribution of the EXERGETIC project is the demonstration that geothermal hybrid systems can be designed to be active grid assets rather than passive baseload generators. Using the ARIES platform, the digital twin was evaluated under multiple grid scenarios, including load following, voltage support at the distribution level, and frequency response at the transmission level. Results show that hybrid geothermal systems can respond effectively to dynamic grid conditions, providing inertia-like behavior, primary frequency support, and voltage regulation through coordinated control. In addition to the performance and grid services capability analysis of geothermal and hybrid geothermal systems, the EXERGETIC project also focused on scalability and techno-economic analysis of geothermal hybrid plants. In particular, for the scalability analysis, machine-learning (ML)-based surrogate models were trained using data generated from the geothermal digital twin under different grid-connected scenarios and plant capacities. These ML models demonstrated strong interpolation and extrapolation capabilities across plant sizes, accurately reproducing both steady-state and transient responses with very low errors. Regarding the techno-economic analysis, plant performance results were integrated with cost models for hybrid geothermal systems, and the levelized cost of electricity (LCOE) was used as the main economic metric to evaluate system performance across a range of system capacities, solar shares, solar multiples, and storage durations. Results indicate that economies of scale significantly reduce geothermal LCOE as plant capacity increases, with large-scale systems (25-100 MW) achieving substantially lower costs than small plants. Hybridization with solar thermal energy and storage further improves economic performance by increasing capacity utilization and enabling flexible dispatch. In addition, thermal storage plays a critical role in reducing LCOE by maximizing geothermal, solar, and stored energy resources. In summary, the results from this project demonstrate that geothermal hybrid systems represent a promising alternative for increasing the energy conversion efficiency of geothermal technologies, contributing to the preservation of geothermal resources, and supporting the transition of geothermal plants from traditional baseload resources into flexible, resilient, and cost-competitive energy conversion technologies.

15 GEOTHERMAL ENERGY↗
DOE OSTI2026

Mapping the Microhardness of Al Matrix in U-7Mo/Al Dispersion Fuels at Medium and High Burn-up by in situ Nanoindentation

The mechanical degradation of the Al matrix in U-7Mo dispersion fuels was detected after irradiation in previous studies, but more details of the degradation mechanism are not clear yet to clarify the root cause. In this study, microhardness mapping in fresh and irradiated U-7Mo dispersion fuels coated with ZrN, with burn-up of 3.35 × 10 21 and 6.28 × 10 21 f/cm 3 , are measured to reveal the hardening of the Al matrix as a function of distance to fuel particles. Micro cubes were pulled by plasma focused ion beam scanning electron microscopy and tested by in situ nanoindenter with Berkovitch tip. Size effects are measured in medium and high burn-up samples. The higher burn-up specimen exhibits a more pronounced size effect when indent load is lower than 20 mN. Size effect correction models are also calculated by plots fitting, which predicts the true microhardness without size effect. Results shows that the hardness of the Al matrix near the ZrN coating has the highest value, then decreases moving away from the coating and become stable when the distance reaches ∼10 µm. The hardness starts increasing again getting closer to the next fuel particle, which makes the hardness distribution between two fuel particles to be “U” shape.

36 - MATERIALS SCIENCE↗
DOE OSTI2026

Data-Driven Discovery and Experimental Validation of Solvent Polarity Effects on Conjugated Polymer Solution-to-Film Assembly Pathways

Understanding how solvent properties influence the solution-to-film assembly of conjugated polymers remains a critical challenge due to the complex and intertwined nature of polymer–solvent interactions. In this study, we integrate a data-driven framework with experimental validation to identify key parameters influencing the assembly and performance of poly[2,5-(2-octyldodecyl)-3,6-diketopyrrolopyrrole-alt-5,5-(2,5-di(thien-2-yl)thieno[3,2-b]thiophene)] (DPP-DTT) in organic field-effect transistors (OFETs). A machine learning (ML) approach identified the normalized Reichardt polarity parameter (E T N ) as a significant descriptor correlated with DPP-DTT hole mobility (μ). Systematic DPP-DTT devices fabricated using solvents across a wide E T N range revealed that higher E T N solvents yield enhanced μ. To elucidate the structural origins of high μ, we conducted comprehensive analyses using UV–vis–NIR spectroscopy and grazing incidence wide angle X-ray scattering (GIWAXS) measurements. The results revealed that films processed from high E T N solvents exhibit reduced paracrystallinity. By analyzing the solution-state behavior using optical microscopy and solution WAXS, we revealed polymer solubility differences in the various solvents and associated distinct polymer assembly pathways, elucidating why the high E T N solvent produces long-range ordered films. Notably, the high E T N solvent shows a pronounced preference for liquid-crystal (LC)-mediated assembly, providing a mechanistic explanation for the enhanced structural order. Therefore, these results demonstrate that solvent polarity, as evaluated by E T N , serves as an important parameter that plays a significant role in the DPP-DTT assembly pathway and resultant solid-state morphology. This work provides a strategy for integrating data science with experiments to identify critical parameters associated with complex polymer systems and helps guide rational process design for high-performance organic electronics.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Heliostat sizing methodology for concentrating solar thermal industrial process heat projects

This study presents a method to obtain a heliostat size that minimizes the levelized cost of heat (LCOH) of a heliostat-based concentrating solar thermal system for applications of solar heating for industrial processes at operating temperatures from 565 to 1550°C. The method extends prior work by embedding a routine for system design that obtains near-optimal subsystem sizes, increasing the fidelity of drive cost functions, and adding an optical performance model to supplement the previously developed cost models, which we update to reflect current pricing trends. An illustrative business case is developed for Daggett, California, targeting specified annual thermal energy outputs of 50 to 400 GWh th . Optical performance is modeled using verified estimates from the literature. A surrogate heliostat cost model, derived from commercial heliostat designs and scaled for production volume, installation, and operations and maintenance costs, is used to develop cost functions. Results show that heliostat size strongly affects the LCOH, producing a characteristic U-shaped trend with a robust near-optimal window of 7-20 m 2 ; the heliostat size producing the lowest project cost in our study grows slightly as the project size increases, and is reduced as the operating temperature increases. The findings in this study are consistent with the general trend of smaller heliostats being deployed at existing projects for high-temperature industrial process heat and reflect the significant reduction in power electronics and other per-heliostat costs. The methodology we propose is general and can be tailored to revised cost curves as the technology continues to evolve.

14 SOLAR ENERGY↗
DOE OSTI2026

The impact of aerosol mixing state on immersion freezing: insights from classical nucleation theory and particle-resolved simulations

Immersion freezing, initiated by ice-nucleating particles (INPs) in supercooled aqueous droplets, plays an important role in the formation of ice crystals within clouds. The efficiency of immersion freezing depends strongly on INP composition and, crucially, on the mixing state – how chemical species are distributed across the particle population. Here, we quantify the impact of aerosol mixing state on immersion freezing using a combined theoretical and particle-resolved modeling approach. We derive analytical expressions for the frozen fraction of internally and externally mixed INP populations based on classical nucleation theory, showing that the frozen fraction is sensitive to whether ice-active species are present in all particles or only in a subset of the population. We introduce a multi-species immersion freezing scheme into the particle-resolved model PartMC, using the water activity-based immersion freezing model (ABIFM) to compute freezing probabilities for mixed-composition particles. To improve computational efficiency, we implement a Binned Tau-Leaping algorithm and demonstrate an order-of-magnitude speedup with minimal accuracy loss. Simulations reproduce the analytical trends in limiting cases and extend the analysis to more general aerosol populations, where mixing state continues to exert a substantial control on frozen fraction. Sensitivity analyses across particle size, species type, and cooling condition reveal that the mixing state effect is most pronounced when small amounts of highly efficient INPs are mixed with less efficient materials. These findings underscore the need to represent aerosol mixing state explicitly in models of heterogeneous ice nucleation to reduce uncertainty in cloud-phase partitioning.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Design and Synthesis of Cubic K 3−2 x Ba x SbSe 4 Solid Electrolytes for K–O 2 Batteries

Developing K-ion conducting solid-state electrolytes (SSEs) plays a critical role in the safe implementation of potassium batteries. In this work, a chalcogenide-based potassium ion SSE is reported, K 3 SbSe 4 , which adopts a trigonal structure at room temperature. Single-crystal structural analysis reveals a trigonal-to-cubic phase transition at the low temperature of 50 °C, which is the lowest among similar compounds and thus provides easy access to the cubic phase. The substitution of barium for potassium in K 3 SbSe 4 leads to the creation of potassium vacancies, expansion of lattice parameters, and a transformation from a trigonal phase to a cubic phase. As a result, the maximum conductivity of K 3−2x Ba x SbSe 4 reaches around 0.1 mS cm −1 at 40 °C for K 2.2 Ba 0.4 SbSe 4 , which is over two orders of magnitude higher than that of undoped K 3 SbSe 4 . This novel SSE is successfully employed in a K–O 2 battery operating at room temperature where a polymer-laminated K 2.2 Ba 0.4 SbSe 4 pellet serves as a separator between the oxygen cathode and the potassium metal anode. Effective protection of the K metal anode against corrosion caused by O 2 is demonstrated.

25 ENERGY STORAGE↗
DOE OSTI2026

Computational simulation of asymmetric phase transformation in cracked Li₇La₃Zr₂O₁₂: Variant selection and chemo-mechanical implications

Coupling between microcracks and phase transformation in ion-conducting ceramics can jointly affect mechanical responses and ion transport. In this work we investigate the cubic-to-tetragonal phase transformation in Li₇La₃Zr₂O₁₂ in the presence of a microcrack under hydrostatic loading and quantify its implications for crack-tip stress concentration and Li-ion transport using phase-field and molecular dynamics simulations. The phase transformation exhibits a strong asymmetry between hydrostatic tension and compression. Under tension, the crack edge nucleates one tetragonal variant that amplifies the crack-tip stress intensity and promotes crack opening. Under compression, the crack tip nucleates a different tetragonal variant that enhances the stress-induced crack-closure tendency. Effective Li diffusivity analysis shows faster transport degradation under compression due to accelerated transformation kinetics, exposing a trade-off between mechanical stability and ionic conductivity. These results highlight the intertwined nature of cracking, phase transformation, and ionic transport in ion-conducting oxides and provide mechanistic insights into chemo-mechanical degradation of solid electrolytes.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Field-tunable BKT and quantum phase transitions in spin-$\frac{1}{2}$ triangular lattice antiferromagnet

Quantum magnetism is one of the most active fields for exploring exotic phases and phase transitions. The recently synthesized Na 2 BaCo(PO 4 ) 2 (NBCP) is an ideal material incarnation of the spin-$\frac{1}{2}$ easy-axis triangular lattice antiferromagnet (TLAF). Experimental evidence shows that NBCP hosts the spin supersolid state with a giant magnetocaloric effect. Theory further predicts that magnetic fields can drive NBCP through Berezinskii-Kosterlitz-Thouless (BKT) and other richer quantum phase transitions. However, detecting these transitions is challenging, as they onset at ultralow temperatures near 60 mK and require high magnetization sensitivity. Using a newly developed gradient force magnetometer in a dilution refrigerator, we mapped the magnetic susceptibility phase diagram down to 30 mK. Our results provide a more comprehensive and accurate understanding of BKT melting of spin supersolidity and several field-tunable quantum phase transitions, which establish NBCP as a model platform for frustrated magnetism and highlight potential applications of its giant magnetocaloric effects.

BKT transition↗
DOE OSTI2024

Derivation of low-energy Hamiltonians for heavy-fermion materials

Here, by utilizing a multiorbital periodic Anderson model with parameters obtained from ab initio band structure calculations, combined with degenerate perturbation theory, we derive effective Kondo-Heisenberg and spin Hamiltonians that capture the interaction among the effective magnetic moments. This derivation encompasses fluctuations via both nonmagnetic 4⁢𝑓 0 and magnetic 4⁢𝑓 2 virtual states, and its accuracy is confirmed through comparison with experimental data obtained from CeIn 3 . The significant agreement observed between experimental results and theoretical predictions underscores the potential of deriving minimal models from first-principles calculations for achieving a quantitative description of 4⁢𝑓 materials. Moreover, our microscopic derivation unveils the underlying origin of anisotropy in the exchange interaction between Kramers doublets, shedding light on the conditions under which this anisotropy may be weak compared to the isotropic contribution.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Energy conversion and transport in molecular-scale junctions

Molecular-scale junctions (MSJs) have been considered the ideal testbed for probing physical and chemical processes at the molecular scale. Due to nanometric confinement, charge and energy transport in MSJs are governed by quantum mechanically dictated energy profiles, which can be tuned chemically or physically with atomic precision, offering rich possibilities beyond conventional semiconductor devices. While charge transport in MSJs has been extensively studied over the past two decades, understanding energy conversion and transport in MSJs has only become experimentally attainable in recent years. As demonstrated recently, by tuning the quantum interplay between the electrodes, the molecular core, and the contact interfaces, energy processes can be manipulated to achieve desired functionalities, opening new avenues for molecular electronics, energy harvesting, and sensing applications. This Review provides a comprehensive overview and critical analysis of various forms of energy conversion and transport processes in MSJs and their associated applications. We elaborate on energy-related processes mediated by the interaction between the core molecular structure in MSJs and different external stimuli, such as light, heat, electric field, magnetic field, force, and other environmental cues. Key topics covered include photovoltaics, electroluminescence, thermoelectricity, heat conduction, catalysis, spin-mediated phenomena, and vibrational effects. Furthermore, the review concludes with a discussion of existing challenges and future opportunities, aiming to facilitate in-depth future investigation of promising experimental platforms, molecular design principles, control strategies, and new application scenarios.

Charge transport↗
DOE OSTI2026

BLEECAM™ (Benchmarking Life Cycle Environmental, Economic, and Social Metrics for Critical and Advanced Minerals and Materials) [SWR-25-125]

The National Laboratory of the Rockies' (NLR) Benchmarking Life Cycle Environmental, Economic, and Social Metrics for Critical and Advanced Minerals and Materials (BLEECAM™) is an open-source, integrated decision-support tool for evaluating the impacts, risks, and trade-offs across U.S. and global materials supply chains. Funded by the U.S. Department of Energy, BLEECAM supports supply chain and market analysis. The tool integrates multi-objective supply chain optimization, system dynamics, network design, lifecycle assessment, techno-economic modeling, and social impact assessment methods to evaluate how supply chains evolve over time, geography, and deployment scenarios. BLEECAM also supports analysis related to energy infrastructure, data centers and digital infrastructure, advanced manufacturing, and other sectors that depend on critical materials.

Khalifa, SherifA. [National Laboratory of the Rock↗
DOE OSTI2025

Theoretical perspective on phase stability and polarization switching in ferroelectric hafnia

Fluorite-based ferroelectric materials are revolutionizing the application space of polar semiconductors. These materials leverage robust polarization of extremely thin films, compatibility with the silicon chip processing, and decades of manufacturing experience to enable a new generation of ferroelectric memory devices. As a new paradigm for ferroelectrics, understanding of phase transitions and the switching mechanism in fluorites is essential both for advancing applications and for fundamental science. In this article, we outline the recent progress that has been made to understand the relative phase stability, phase transition and order parameter coupling, ferroelectric switching through unique nucleation and growth processes, and how defects affect these phases and processes. The main challenges, opportunities, and next steps for leveraging these materials for next-generation devices are reviewed.

Condensed Matter Physics↗
DOE OSTI2023

The UCLA Cosmochemistry Database

Abstract The UCLA Cosmochemistry Database was initiated as part of a data-rescue and -storage project aimed at archiving a variety of cosmochemical data acquired at University of California, Los Angeles (UCLA). The data collection includes elemental compositions of extraterrestrial materials analyzed by UCLA cosmochemists over the last five decades. The analytical techniques include atomic absorption spectrometry (AAS) and neutron activation analysis (NAA) at UCLA. The data collection is stored on the Astromaterials Data System (Astromat). We provide both interactive tables and downloadable datasheets for users to access all data. The UCLA Cosmochemistry Database archives cosmochemical data that are essential tools for increasing our understanding of the nature and origin of extraterrestrial materials. Future studies can reference the data collection in the examination, analysis, and classification of newly acquired extraterrestrial samples.

Science & Technology - Other Topics↗
DOE OSTI2023

Four Years of Type Ia Supernovae Observed by TESS: Early-time Light-curve Shapes and Constraints on Companion Interaction Models

Abstract We present 307 type Ia supernova (SN) light curves from the first 4 yr of the Transiting Exoplanet Survey Satellite mission. We use this sample to characterize the shapes of the early-time light curves, measure the rise times from first light to peak, and search for companion star interactions. Using simulations, we show that light curves must have noise <10% of the peak flux to avoid biases in the early-time light-curve shape, restricting our quantitative analysis to 74 light curves. We find that the mean power-law index t β 1 of the early-time light curves isβ 1 = 1.93 ± 0.57, and the mean rise time to peak is 15.7 ± 3.5 days. The underlying population distribution forβ 1 may instead consist of a Gaussian component with mean 2.29, width 0.34, and a long tail extending to values less than 1.0. We find that the data can rarely distinguish between models with and without companion interaction models. Nevertheless, we find three high-quality light curves that tentatively prefer the addition of a companion interaction model, but the statistical evidence for the companion interactions is not robust. We also find two SNe that disfavor the addition of a companion interaction model to a curved power-law model. Taking the 74 SNe together, we calculate 3σupper limits on the presence of companion signatures to control for orientation effects that can hide companions in individual light curves. Our results rule out common progenitor systems with companions having Roche lobe radii >31R ⊙ (separations >5.7 × 10 12 cm, 99.9% confidence level) and disfavor companions having Roche lobe radii >10R ⊙ (separations >1.9 × 10 12 cm, 95% confidence level). Lastly, we discuss the implications of our results for the intrinsic fraction of single degenerate progenitor systems.

Astronomy & Astrophysics↗
DOE OSTI2026

Sierra/SD - Users Manual - 5.32

Sierra/SD provides a massively parallel implementation of structural dynamics finite element analysis, required for high-fidelity, validated models used in modal, vibration, static and shock analysis. This document provides a user's guide to the input for Sierra/SD. Details of input specifications for the different solution types, output options, element types and parameters are included.

Beale, Dagny Marie [Sandia National Laboratories (↗
DOE OSTI2026

Anion-dependent phase behavior of methylimidazolium-based ionic liquids mixed with water: Correlation between local molecular structure and mesoscale behaviors

The mesoscopic phase behavior of decylmethylimidazolium (C 10 mim) ionic liquids (ILs) bearing three monovalent anions—thiocyanate (SCN − ), nitrate (NO 3 − ), and chloride (Cl − )—mixed with water at relatively high IL contents (50–95 wt%) was investigated. Small-angle and wide-angle X-ray scattering (SAXS/WAXS) were employed to follow the evolution of both local and mesoscale structures across this composition range. In the absence of water, C 10 mimSCN and C 10 mimCl behaved as disordered liquids, whereas C 10 mimNO 3 spontaneously formed a hexagonally ordered cylindrical mesophase and displayed a sticky-solid macroscopic appearance. Upon addition of water, C 10 mimSCN remained a viscous liquid and only weakly ordered lamellar domains were observed. This limited ordering is attributed to the lack of hydrogen bonding and weak interaction energy between SCN − anions and C 10 mim + cations. In contrast, the trigonal-planar NO 3 − and point-like Cl − anions promoted the formation of well-defined hexagonal mesophases up to 35–45 wt% water. The formation of hydrogen bonding of the two anions with imidazolium ring protons likely enabled the creation of compact ion clusters that effectively exclude water molecules from the immediate vicinity of the IL aggregates. These findings demonstrate that the interaction energy between ion pairs dictate IL-water interactions and therefore control the transition from disordered liquids to ordered mesophases in IL/water mixtures with high IL contents. In conclusion, the combined SAXS/WAXS analysis reveals a correlation between local intermolecular structure and the emergence of mesoscopic order, providing a systematic framework for tailoring mesoscale structures in alkylimidazolium-based IL/water systems.

36 MATERIALS SCIENCE↗
DOE OSTI2025

A rigorous framework for an improved Messinger/Myers model of ice accretion under conditions of variable property and unsteady aircraft icing

We analyse the Messinger/Myers model by critically evaluating simplifying assumptions through a rigorous formulation of the rime ice accretion process. We explore the effects of both constant and variable ice density and thermal conductivity, along with the effects of sublimation from the ice surface. The effects of key factors such as droplet impact rate, ambient temperature relative to the freezing temperature and the temperature difference between the ambient air and the airfoil surface are examined. Under these varying conditions, the present rigorous formulation is used to assess the significance of unsteady effects, variable ice properties and sublimation. We observe that the Myers model performs remarkably well in certain icing situations and analyse the reasons for this strong performance. We also show that partially relaxing the model’s assumptions can lead to poorer performance. The Myers model can lead to overprediction of ice surface temperature and correspondingly underprediction of transition time under conditions of relatively weak sublimation and surface cooling. A modified Myers model is presented, which can be used to recover near-perfect results under widely varying icing conditions of relevance. This article is part of the theme issue ‘Heat and mass transfer in frost and ice’.

Science & Technology - Other Topics↗