Search NASASearch

SEARCH · Search NASA

Results for “Design”

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.

990 records · Page 45

Spinel high-entropy oxides (FeNiCrMnZnX) 3 O 4 (X = Al, mg) as anode materials for high-performance lithium-ion batteries

To address the high cost, cobalt dependency, and resource constraints typical of conventional high-entropy oxide (HEO) anodes, this study reports the successful synthesis of two Co-free, six-component spinel-type HEOs(FeNiCrMnZnAl) 3 O 4 (HEO-Al) and (FeNiCrMnZnMg) 3 O 4 (HEO-Mg), via a sol-gel method. The distinct effects of Al 3+ and Mg 2+ incorporation on the electrochemical performance and lithium storage kinetics were systematically investigated. XRD, Raman, and TEM characterizations confirm that both materials possess a pure spinel phase, uniform particle size, and homogeneous elemental distribution. Notably, electrochemical evaluations reveal that HEO-Al delivers a superior reversible capacity of 480.7 mAh g −1 after 100 cycles at 0.1 A g −1 , and maintains 354.5 mAh g −1 after 1000 long-term cycles at 1 A g −1 , significantly outperforming HEO-Mg. Kinetic analysis indicates that HEO-Al exhibits lower charge transfer resistance, a higher Li + diffusion coefficient, and a pseudocapacitive contribution of up to 82%. Furthermore, these findings demonstrate that Al substitution effectively optimizes the structural stability and lithium storage kinetics of Co-free HEOs, providing a viable strategy for designing low-cost, highly stable HEO anode systems.

Anode materials

Fluorophilic Sigma(σ)‐Lock Self‐Healable Copolymers

Although F-Containing molecules and macromolecules are often used in molecular biology to increase the binding with Lewis acidic groups by introducing favorable C−F dipoles, there is virtually no experimental evidence and limited understanding of the nature of these interactions, especially their role in synthetic polymeric materials. These studies elucidate the molecular origin of inter- and intra-Chain interactions responsible for self-healing of F-Containing copolymers composed of pentafluorostyrene and n-butyl acrylate units (p(PFS/nBA). Guided by dynamic surface oscillating force (SOF) and spectroscopic measurements supported by molecular dynamics (MD) simulations, these studies show that the reformation of σ-σ orbitals in −C−F of PFS and CH 3 CH 2 − of nBA units enables the recovery of entropic energy via fluorophilic-σ-lock van der Waals forces when PFS/nBA molar ratios are ~50/50. The strength of these interactions determined experimentally for self-healable PFS/nBA compositions is in the order ~0.3 kcal/mol which primarily comes from fluorophilic-σ-lock (~70 %) contributions. These interactions are significantly diminished for non-self-healable counterparts. Strongly polarized −C−F σ orbitals create lateral dipolar forces enhancing the affinity towards −C−H orbitals, facilitating energetically favorable interactions. Entropic recovery driven by non-Covalent bonding offers a valuable tool in designing materials with unique functionalities, particularly self-healable batteries and energy storage devices.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

Micrometer: Micromechanics transformer for predicting full field mechanical responses of heterogeneous materials

Predicting mechanical responses of heterogeneous materials across scales remains a significant challenge. Traditional computational methods often struggle with complex and multiscale nature of these materials, limiting their effectiveness in real-world applications. Here, in this paper, we introduce Micrometer, a vision transformer based deep learning model designed to predict full field mechanical responses of heterogeneous materials, bridging the gap between computer vision and solid mechanics problems. We show that Micrometer, trained on a large-scale high-resolution dataset of 2D fiber-reinforced composites, can achieve state-of-the-art performance in predicting microscale strain fields across a wide range of material properties and loading conditions. Our model demonstrates accuracy and computational efficiency in applications such as computational homogenization and multiscale modeling, reducing computational time by up to two orders of magnitude compared to conventional numerical solvers while maintaining less than 1 % errors in predicting macroscale stress fields. Furthermore, we showcase Micrometer’s adaptability through transfer learning experiments on new materials with limited data, highlighting its potential to tackle diverse scenarios in computational solid mechanics. These results represent a significant step towards AI-driven innovation in materials science, addressing the limitations of traditional numerical methods and paving the way for more efficient simulations of heterogeneous materials across various industrial applications.

Composite materials

Observation of Extraordinary Vibration Scatterings Induced by Strong Anharmonicity in Lead‐Free Halide Double Perovskites

Abstract Lead‐free halide double perovskites provide a promising solution for the long‐standing issues of lead‐containing halide perovskites, i.e., the toxicity of Pb and the low stability under ambient conditions and high‐intensity illumination. Their light‐to‐electricity or thermal‐to‐electricity conversion is strongly determined by the dynamics of the corresponding lattice vibrations. Here, the measurement of lattice dynamics is presented in a prototypical lead‐free halide double perovskite(Cs 2 NaInCl 6 ). The quantitative measurements and first‐principles calculations show that the scatterings among lattice vibrations at room temperature are at the timescale of ≈1 ps, which stems from the extraordinarily strong anharmonicity in Cs 2 NaInCl 6 . Further the degree of anharmonicity of each type of atom is quantitatively characterized in the Cs 2 NaInCl 6 single crystal, which stems from the interatomic forces, and demonstrate that this strong anharmonicity is synergistically contributed by the bond hierarchy, the tilting of the NaCl 6 and InCl 6 octahedral units, and the rattling of Cs + ions. Consequently, the crystalline Cs 2 NaInCl 6 possesses an ultralow thermal conductivity of ≈0.43 W mK −1 at room temperature, and a weak temperature dependence ofT −0.41 . These findings uncovered the underlying mechanisms behind the dynamics of lattice vibrations in double perovskites, which can largely benefit the design of optoelectronics and thermoelectrics based on halide double perovskites.

Chemistry

Performance of Heterostructural TaC/AlGaN Schottky Diodes Based on First Principles Electronic Structure Properties

Advances in ultra-wide bandgap materials, such as high Al-content AlxGa1-xN (AlGaN), are essential for next generation power electronics, but the requirement for lattice matched substrates is currently a significant obstacle. Recently, conductive TaC has emerged as a promising virtual substrate for AlGaN heteroepitaxy, with wurtzite (0001) AlxGa1-xN lattice-matched to rocksalt (111) TaC at x ~ 0.5. Thus, understanding and controlling the electronic properties of the TaC/AlGaN interface is key for developing technological applications based on TaC/AlGaN devices. Using density functional theory and electronic structure calculations, we here investigate TaC/Al0.5Ga0.5N interfaces, where we include explicit alloy models in the slab calculations. We predict the Schottky barrier height and the electric field discontinuity resulting from interface charges. Considering all possible combinations of (Ta, C) substrate termination, (Al/Ga, N) nucleation, and (Al/Ga, N) polarity, we construct a chemical potential phase diagram to identify the stable interfaces that can be accessed through variation of the synthesis conditions. The predicted interface electronic properties are implemented in device performance simulations to demonstrate a practical design for a strain-free, high-efficiency TaC/AlGaN Schottky diode with a low barrier height and without interface charges, underscoring the potential of TaC as a substrate for ultra-wide bandgap devices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC

Advanced Model Development for Large Eddy Simulation of Oxy-Combustion and Supercritical Carbon Dioxide Power Cycles

A joint experimental and numerical study is performed to observe the characteristics of a supercritical carbon dioxide turbulent mixing layer in the presence of strong nonlinearities in the thermodynamic and transport properties. A bespoke experimental setup is designed and employed for this purpose and provides insight into macroscopic mixing behavior. The mixing is experimentally observed using two techniques: shadowgraphy and spontaneous Raman scattering. Qualitative and quantitative intensity fields obtained via these techniques yield instantaneous and mean density data. Spanwise temperature data is also collected using analogue resistance temperature detectors. These measurements are used to quantify the level of mixed material within the field. The experimental data are supplemented by a companion high-fidelity numerical study. The numerical results are obtained through fully resolved, three-dimensional direct numerical simulation. The numerical dataset permits observation of the near-field mixing characteristics, which are difficult to measure experimentally due to the rapid dynamics and sharp thermophysical gradients in this area. Qualitative field visualizations are presented, followed by quantitative mixed material results and observations regarding thermodynamic property trends at select locations within the field. One-dimensional spectra of the turbulent kinetic energy and solenoidal dissipation are provided to observe the spectral characteristics of the flow. Reynolds stress anisotropy is analyzed graphically through anisotropy invariance maps (Lumley triangles). The mixing quantification, spectral data and anisotropy analysis of a flow at these thermodynamic conditions represent the main outcomes of the work.

20 FOSSIL-FUELED POWER PLANTS

Optimizing Porous Transport Layer Porosity for Proton Exchange Membrane Water Electrolysis

An empirical model is presented that describes anode-side losses related to porous transport layer (PTL) morphology in proton exchange membrane water electrolysis (PEMWE). The model is based on an advanced voltage breakdown analysis that links various overpotentials to PTL morphology. Custom Ti PTLs, spanning uncommonly low porosities (22 - 31%), were fabricated and analyzed with X-ray CT to obtain pore and particle size distributions. Particle size distributions were consistent across samples with an average particle diameter of 12.0?..mu..m, whereas average pore diameters ranged from 6.0 to 7.0?..mu..m. The PTLs were tested in standard PEMWE cell assemblies with anode catalyst loadings of 0.1 mgIr cm-2 to obtain polarization curves, electrochemical impedance spectra, and augmented Tafel analysis. The PTL-dependent anode side losses were deconvoluted and assigned to excess utilization, concentration, ion transport resistance, and electrical contact resistance overpotentials. The data and model reveal an optimal 20 - 28% PTL porosity region where utilization and contact resistance overpotentials are minimized without triggering concentration and ion transport losses related to water deprivation. The optimal PTL porosity depends on the operating current density and is demonstrated at realistic PEMWE water flow rates to establish PTL design guidance for operation at scale.

08 HYDROGEN

Rewet Temperature Correlations for Liquid-Nitrogen Boiling Pipe Flows Across Varying Flow Conditions and Orientations

In many convective liquid-vapor phase change heat transfer engineering applications, cryogenic fluids are widely used in industrial processes, spacecraft and cryosurgery systems, and so on. For example, cryogens are usually used as liquid fuels such as liquid hydrogen and oxygen in the rocket industry, liquid nitrogen (LN2) and helium are frequently used to cool superconducting magnetic device for medical applications. In these systems, proper transport, handling, and storage of cryogenic fluids are of extreme importance. Among all the cryogenic transport processes performed in room temperatures, quenching, also termed chilldown, is a unavoidable initial, transient phase-change heat transfer process that brings the system down to the cryogenic condition. The Leidenfrost temperature or rewet temperature that signals the end of film boiling is practically considered the completion point of a quenching process. Therefore, rewet temperature has been considered the most important parameter for the engineering design of cryogenic thermal management systems. As most of the previous correlations for predicting the Leidenfrost temperature and the rewet temperature have been basically developed for water, they are shown to disagree with recent liquid nitrogen pipe chilldown experiments in upward and downward flow directions over a wide range of flow rates, pressures, and degrees of inlet subcooling. In addition to a complete review of the literature, two new correlations are presented in this work, one based on bubble growth and another based on the theoretical maximum limit of superheat. Each correlation performs well over the entire data set.

liquid nitrogen

Smart building HVAC control challenge: experience and solutions from the ADRENALIN project

A smart building HVAC control competition crowdsourced and compared algorithms on fair and equal ground using the standardized BOPTEST framework. The competition attracted 138 participants, but only 9% submitted valid solutions for the final stage, highlighting the complexity of advanced HVAC control design. The winning solutions showed significant potential to reduce energy use and cost by shifting demand, without compromising occupant comfort. Across scenarios, thermal energy cost reductions of 36–76% relative to a baseline, were achieved. In peak heat periods, the cost reduction leveraged limited energy use reduction (0–15%), but more significant energy price reduction (34–62%). This shows smart controls' ability to avoid as much as possible consumption during the morning peak hours, when spot prices are tendentially the highest. Hosting the competition has highlighted challenges in creating competitions that both are fair and promotes solutions that are transferable to real life implementation.

BOPTEST

Guidelines for Schlieren Systems at Langley Research Center

The original Langley Working Paper (LWP 448) published on July 27th, 1967, provided guidance to NASA Langley Research Center personnel on how to set up conventional path-integrated schlieren flow visualization systems and what pitfalls could be expected with such setups. The guidance and information contained in the document continues to be used for schlieren setups at NASA Langley Research Center to this day. The most commonly used method of flow visualization in supersonic wind tunnels is the schlieren method, or method of striae. At the Langley Research Center, this method is so commonly used that it has come to be regarded somewhat as merely another piece of instrumentation which is to be maintained by service people, and which is too mysteriously complicated to be touched by tunnel personnel. This should not be the case. High-quality results from a schlieren system can only be obtained by having personnel on the tunnel staff who are well-informed on the subject and who have the time and the interest to do a good job of adjustment. This paper is a discussion in practical terms, of some of the points to be considered and the pitfalls that may be encountered in designing, setting up and using a working schlieren system. It is directed specifically toward helping specify and install a Z-type mirror system. It is hoped that the material will be of assistance to facilities throughout the Center and NASA and that users of the schlieren method will be inspired to the point of training people in the adjustment of the systems in use at their facilities.

Langley Working Papers

Peptoid-Based Nanosheets Exhibiting Broad Antiviral Activity Against Enveloped RNA Viruses

Enveloped RNA viruses, such as Influenza A (H1N1) and Sindbis virus, pose persistent global health threats due to their high mutation rates, efficient transmission, and frequent drug resistance. By mimicking host cell membrane receptors, multivalent virus inhibitors can block viral attachment, making them promising broad-spectrum antiviral agents. However, most of existing antivirals are often limited by strain specificity, short-lived efficacy, and toxicity. Here, we introduce a broad-spectrum antiviral platform based on highly tunable and biocompatible two-dimensional nanomembranes (2DNMs) self-assembled from amphiphilic peptoids, operating via a non-genomic, mutation-insensitive mechanism. By varying peptoid sequence, we design and synthesize over twenty different 2DNMs with various surface charge and high density of viral-attachment ligands (VALs). The self-assembled architecture of these stable 2DNMs provides cooperative noncovalent multivalent binding to virus particles that result in effective inhibition of viral infection. Screening of variants identified three leads that potently suppressed Influenza A (H1N1) and Sindbis virus infection across median tissue culture infectious dose (TCID50), plaque, RT–qPCR, and immunofluorescence assays, while maintaining >90% cell viability. These nanosheets significantly reduced infectious titers, viral RNA replication, and intracellular viral protein expression, indicating inhibition at early stages of viral entry and propagation. The sequence programmability, chemical robustness, and mutation-insensitive antiviral activity distinguish 2DNMs from traditional antivirals and positions them as a versatile materials platform for antiviral coatings, protective barriers, and prophylactic biomedical applications.

Influenza A virus

Preliminary Assessment for the Electric Load Shifting Potential of Integrating Thermal Energy Storage with Heat Pumps in Residential Buildings in Texas of United States

The widespread adoption of electric-driven heat pumps for heating and cooling is expected to significantly increase electric demand. Cooling electric demand will result in an increase in peak hours electric loads, placing additional strain on the grid during peak hours. This challenge, combined with the current rapidly increasing demand from data centers, exacerbates the electric demand duck curve problem. Integrating thermal energy storage (TES) with heat pump can shift electric use for heating and cooling from peak to off-peak hours of the electric grid, which can help flatten the daily electric demand profile of a building. This paper presents a novel design that integrates heat pump with TES (HP-TES), which uses phase change materials (PCM). Heat pump charges TES by melting or freezing PCM during off-peak hours when there is no thermal demand from the building. TES is then discharged (i.e., by freezing or melting PCM) during peak hours to provide a more favorable heat source or heat sink for the heat pump to meet the thermal demand of the building with lower electricity use than conventional air-source heat pumps. Computer simulations were developed to predict the performance of HP-TES applied to a typical single-family house in the US. The building-level simulation results were scaled up to preliminarily assess the aggregated impacts of deploying HP-TES across all single-family houses in Texas of the United States, including reduction of peak demand of the electric grid.

Anees, Fady [ORNL]

Microstructure, mechanical, and thermal properties of compositionally complex (Hf,Zr,Nb,Ti)B 2 ‒LaB 6 ceramics

Novel compositionally complex borides, (Hf,Zr,Nb,Ti)B 2 and (Hf,Zr,Nb,Ti)B 2 ‒LaB 6 , were fabricated using spark plasma sintering process. (Hf,Zr,Nb,Ti)B 2 ‒LaB 6 exhibits a dual-phase microstructure, in which (Hf,Zr,Nb,Ti)B 2 is a primary phase with the hexagonal structure and LaB 6 is a secondary phase with a cubic structure. The mechanical properties of both (Hf,Zr,Nb,Ti)B 2 and (Hf,Zr,Nb,Ti)B 2 ‒LaB 6 are comparable, with a combination of high hardness and moderate fracture toughness. Thermal diffusivity and conductivity of (Hf,Zr,Nb,Ti)B 2 are much lower than the individual transition metal borides but are significantly increased by the addition of LaB 6 . Herein, it is implied that the thermal properties of boride ceramics can be controlled through the appropriate design of principal metal element compositions.

borides

Gear Test Assembly - Experimental Testing and Analysis of Gears and Bearings - FY2026

The Gear Test Assembly (GTA) is an experimental test apparatus built and installed in the Mechanisms and Engineering Test Loop (METL) at Argonne National Laboratory (ANL). While designed to accommodate a variety of components intended for use in advanced compact in-vessel transfer machine, the focus of GTA been the testing of large radial spur gears. The performance of the components used in GTA motivates the choices made in the design, components, and material choices for the forthcoming METL test articles; the Gripper Test Assembly (GrTA), and the Bearing Test Assembly (BTA). To date, GTA has completed nine experimental testing campaigns and achieved over 21 million revolutions under prototypic sodium fast reactor loads and operating conditions. This report provides an overview and results for the two most recent experiments Campaign #8 and Campaign #9. `Campaign #8 was the first campaign to make use of a set of nickel alloy radial spur gears, made from Hastelloy C-22HS, and combined with heat-treated tapered roller bearings. The bearings and torque profile used allowed Campaign #8 to be directly compared with Campaign #1 where Inconel 718 gears were used instead. Campaign #8 attained 1,058,880 shaft revolutions within 23% of the 1,314,855 revolutions achieved in Campaign #1. These results demonstrate that the heat-treated tapered roller bearings show consistent increased lifetimes, approximately three times longer than for non heat-treated taper rollers, and that the particular material of the gear does not play a substantial role in performance. Similar to Campaign #7 where bearing failure resulted in fragmentation and transport of bearing material into the gear teeth led to tooth damage in the Inconel gears, a similar process led to the damage of the Hastelloy gears. Subsequent NDE analysis showed that the damage was only limited to the tooth surface and does not extend into the gear interiors. Based on the NDE reports, the Hastelloy gears can be reconditioned and put back into service for future campaigns as was done for the Inconel gears. Reconditioned Inconel 718 gears damaged in Campaign #7 were returned to service in Campaign #9 and paired with non heat-treated tapered roller bearings. In Campaign #9 a novel set of operating conditions were implemented whereby the gears and bearings were cycled between hot, submerged refueling conditions where fuel handling operations occurred and cold, dry, and inert conditions where the sodium was drained and the vessel was allowed to cool simulating long, ex-vessel storage of the fuel handling machine. A total of eight of these operating cycles were achieved before a bearing failure occurred after the accumulation of 475,692 shaft revolutions were achieved. During Campaign #9, one of the drive motors was damaged and over the course of a 240 day standby period, the entire GTA drive system was upgraded from a 480VAC system to a 240VAC system that provided the opportunity to consolidate wall-space in the METL facility and provide more space for additional test articles being installed in the coming year. The over 475,000 shaft revolutions achieved in Campaign #9 eclipses the previous record for lifetime of non heat-treated tapered roller bearings held by Campaign #3 of 392,000 shaft revolutions. Further analysis and future campaigns hope to shed light on the cause of this approximately 17% increase in lifetime.

42 ENGINEERING

Phase Stability and Electrochemical Performance of La-Site-Doped Li6La3Zr0.5Nb0.5Ta0.5Hf0.5O12 High-Entropy Garnets

We investigate La-site substitution in the high-entropy garnet Li6La3Zr0.5Nb0.5Ta0.5Hf0.5O12 (LLZNTH) using Ba2+, Sr2+, and Sm3+ to elucidate how dopant governs phase stability, Li-site distribution, and electrochemical behavior. X-ray diffraction shows that Sr2+ is incorporated homogeneously into the garnet lattice, whereas the larger Ba2+ and smaller Sm3+ ions partially exceed the structural tolerance, generating secondary phases. Nevertheless, the Sm-doped composition (x = 0.05) exhibits the highest room-temperature ionic conductivity (2.7 × 10–4 S cm–1). Neutron powder diffraction reveals that Sm substitution drives a redistribution of Li+ from the tetrahedral 24 d sites into the higher-mobility 96 h positions, enhancing the connectivity of the three-dimensional Li-ion migration network. A Sm-doping series (x = 0.01–0.05) further shows that only sufficiently high Sm levels induce this redistribution, whereas lower concentrations retain Li arrangements similar to the undoped garnet. Critical current density measurements demonstrate that La-site dopants also influence interfacial stability against Li metal, underscoring a trade-off between bulk transport enhancement and mechanical robustness. Collectively, these findings reveal that in high-entropy garnets improved ionic conductivity can originate not only from phase-pure structures but also from targeted modification of the Li sublattice, even when accompanied by secondary phases, offering a compositional design principle for garnet electrolytes.

Li, Chang [Mechanical Engineering, School of Scien

Carbon source–driven metabolic and regulatory remodeling defines phenomic states in Lipomyces starkeyi

Lipomyces is a genus of oleaginous yeasts with potential for contributing to reliable biomanufacturing supply chains. However, progress in advanced strain designs and engineering efforts are still constrained by a lack of understanding of the underlying molecular drivers of Lipomyces phenotypes. To address this gap, we collected a suite of multi-omic data to dissect how carbon source availability reshapes the metabolic network, lipid allocation, and regulatory architecture of Lipomyces starkeyi. We observed that glucose promotes biosynthetic and proliferative processes supported by abundant energy and carbon intermediates, xylose enhances redox-balancing mechanisms centered on the pentose phosphate pathway, and glycerol activates respiratory metabolism, ß-oxidation, and the glyoxylate cycle. Lipid species distributions remained consistent in both nitrogen replete and depleted conditions across the carbon sources, indicating robust production mechanisms. Regulatory protein identification and network analysis revealed glycerol-driven respiratory growth favors regulatory programs integrating stress tolerance, redox balance, and lipid-associated metabolism, whereas xylose growth activates compensatory transcriptional responses aimed at maintaining mitochondrial function. Nitrogen limitation modulates the strength of these responses but does not fundamentally alter their direction, reinforcing carbon source as the dominant driver of regulatory architecture. Taken together, this data enhances the understanding of Lipomyces molecular rearrangements and provides a foundation for further development of predictive phenotypic tools in this genus.

Biotechnology

A Full-Induction Magnetohydrodynamics Solver for Liquid Metal Fusion Blankets in Vertex-CFD

Multiphysics modeling of liquid metal fusion blankets, which produce tritium and convert energy of neutrons created via fusion reactions into heat, is crucial for predicting performance, ensuring structural integrity, and optimizing energy production. While traditional blanket modeling of liquid metal flows during normal steady operating conditions commonly employs the inductionless approximation of the magnetohydrodynamics (MHD) equations, transient scenarios, when the plasma-confining magnetic field varies on millisecond time scales, require a full-induction MHD approach that dynamically evolves the magnetic field via the time-dependent induction equation. This paper presents the formulation, implementation, and initial verification of a full-induction MHD solver integrated within the open-source Vertex-CFD framework, which aims to achieve tight multiphysics coupling, a flexible software design enabling easy extension and addition of physics models, and performance portability across computing platforms. The solver utilizes finite element spatial discretization, implicit Runge–Kutta time integration, and an inexact Newton method to solve the resulting discrete nonlinear system, leveraging Trilinos packages for efficient computation. Verification against selected benchmark problems demonstrates accuracy and robustness of the solver. Furthermore, when the solver is applied to an idealized blanket model in 2.5D and full 3D, results obtained with Vertex-CFD are in good agreement with recently published quasi-2D simulations. These findings establish a computational foundation for future simulations of transient MHD phenomena in liquid metal blankets with Vertex-CFD, and open avenues for future extensions and performance optimizations.

Endeve, Eirik [ORNL] (ORCID:0000000312519507)

Time-History Statistics of Soot Formation in A Model Gas Turbine Combustor

Soot formation is a complex dynamic and intermittent process determined by properties of the fuel, combustor design, and combustor operation. Although the major steps in soot formation (i.e., formation of precursors, inception, growth and evolution) are similar for a variety of carbonaceous fuels, applications, and operating conditions, it remains unclear when the temporal transition between these steps occurs. An engineering prediction tool coupled with computational fluid physics (CFD), therefore needs to accurately model all these complex steps. To develop such a model, we propose the time-history concept for understanding the time dependency of soot formation as a function of local properties (i.e., temperature, velocity, local fuel air ratio, etc.). We continue our previous work with modeling the DLR aero-combustor [1] with our updated in-house CFD code, Open National Combustion Code (OpenNCC), that now includes a Multiple Time-Scale Flamelet Progress Variable approach and a the semi-empirical two-equation soot model. We injected massless tracer particles upstream of the injector region of the combustor to collect time-history statistics of the solution variables. The correlations between the collected statistics with respect to the experimental soot volume fraction data showed that time-history effect of certain flow variables, including turbulent kinetic energy (TKE), and multiple species is indeed important for soot formation. We then conducted a time-history based correlation analysis to determine the key species and the concentration ranges critical for soot formation (C6H5-based nucleation, acetylene-based surface growth, and oxidation with OH and O2). Based on the time-history correlation coefficient (THCC) analysis, we propose possible modifications to improve the current two-equation model.

LES