Search NASASearch

SEARCH · Search NASA

Results for “Numerical Algorithms”

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.

226 records · Page 11

Coordinated Natural Gas and Electric Planning: Case Studies of Current Approaches and Practices

This paper examines how natural gas and electric utilities across eight U.S. states and two Canadian provinces are beginning to coordinate historically separate planning processes in response to growing system interdependencies, policy mandates, aging infrastructure, and changing customer energy choices. Electricity planning has long relied on robust integrated resource planning frameworks that weigh numerous objectives, risks, and costs. Natural gas planning, on the other hand, is typically less transparent and more narrowly focused on safety and system integrity. As economic, reliability, and policy drivers place new and shared demands on both systems, jurisdictions and utilities are experimenting with approaches such as coordinated forecasting, non pipeline alternatives, and combined planning pilots. A few have issued regulatory or statutory directives for greater data sharing and methodological alignment. Case studies from British Columbia, California, Colorado, Illinois, Massachusetts, Minnesota, New York, Québec, Rhode Island, and Washington illustrate a wide range of approaches to rethinking siloed planning. The paper identifies several common themes across the jurisdictions examined. It provides observations on the methods, processes, and organizational steps that may be needed in the future to address the challenges being faced by states and utilities. Lastly, it identifies some initial steps that states and utilities can consider if they would like to pursue more integrated, cost-effective, policy-aligned energy system planning.

03 NATURAL GAS

Review of European Microgravity Measurements

AA In a French/Russion cooperation, CNES developed a microgravity detection system for analyzing the Mir space station micro-g-environment for the first time. European efforts to characterize the microgravity (1/9) environment within a space laboratory began in the late seventies with the design of the First Spacelab Mission SL-1. Its Material Science Double Rack was the first payload element to carry its own tri-axial acceleration package. Even though incapable for any frequency analysis, the data provided a wealth of novel information for optimal experiment and hardware design and operations for missions to come. Theoretical investigations under ESA contract demonstrated the significance of the detailed knowledge of micro-g data for a thorough experiment analysis. They especially revealed the high sensitivity of numerous phenomena to low frequency acceleration. Accordingly, the payloads of the Spacelab missions D-1 and D-2 were furnished with state-of-the-art detection systems to ensure frequency analysis between 0.1 and 100 Hz. The Microgravity Measurement Assembly (MMA) of D-2 was a centralized system comprising fixed installed as well as mobile tri-axial packages showing real-time data processing and transmission to ground. ESA's free flyer EURECA carried a system for continuous measurement over the entire mission. All EURECA subsystems and experimental facilities had to meet tough requirements defining the upper acceleration limits. In a French/Russion cooperation, CNES developed a mi crogravity detection system for analyzing the Mir space station micro-g-environment for the first time. An approach to get access to low frequency acceleration between 0 and 0.02 Hz will be realized by QSAM (Quasi-steady Acceleration Measurement) on IML-2, complementary to the NASA system Spacelab Acceleration Measurement System SAMS. A second flight of QSAM is planned for the Russian free flyer FOTON.

H Hamacher

Mechanism of H 2 plasma-enabled reduction of hematite thin films

Hydrogen plasma is gaining significant interest as a promising pathway for direct iron ore reduction and for lowering process temperatures but the reduction mechanism remains poorly understood. In this work we analyzed the plasma and thermal reduction of thin-film hematite (Fe 2 O 3 ) at temperatures below 340 °C using X-ray diffraction and scanning electron microscopy with energy-dispersive X-ray spectroscopy. Plasma reduced the incubation period by an order of magnitude and increased the reduction rate by a factor of 2.6 compared to thermal reduction. Plasma-produced H-atoms facilitate the formation of numerous iron nucleation sites, bypassing the energetically unfavorable dissociative adsorption of H 2 on iron oxide. These iron nuclei can autocatalyze the reduction of the surrounding hematite via a hydrogen spillover mechanism. Our results demonstrate that plasma-derived H-atoms primarily impact the initial nucleation-limited stage. These new insights provide a mechanistic framework that can aid the implementation and optimization of hydrogen plasma-assisted iron oxide reduction at reduced temperatures.

08 HYDROGEN

Recirculating and Sub-Atmospheric Rapid Cycle Amine Swing Bed Testing at Various Metabolic Profiles, Temperatures, and Half-Cycle Times

Sixty years since the first spacewalk, NASA’s Extravehicular Mobility Unit (xEMU) exploration space suit has undergone numerous design improvements and iterations. As commercial manufacturers develop next-generation systems to support increasingly complex missions, including Moon-to-Mars architecture prototyping and the design of the Martian Portable Life Support System (PLSS), material selection and individual subunit testing under realistic metabolic and pressure conditions are paramount. To determine whether a reusable carbon dioxide (CO 2 ) scrubber design and adsorbent can be utilized in new suit configurations, an all-encompassing test rig is required to compare existing and newly developed technologies. The test rig must be capable of sub-atmospheric testing, calibrated humid CO 2 dosing, and facilitation of the evaluation of sorbents across a wide range of operational requirements. Recently, XploSafe has developed an Extravehicular Mobility Unit (xEMU) testing rig specifically designed to explore CO 2 and humidity control materials for spacesuits. This closed-loop recirculating swing bed (Rapid Cycle Amine Test Rig) was utilized to compare adsorbent utility across a wide range of xEMU operating conditions. Following successful testing with Xplo-SA9T over multiple simulated eight-hour Extravehicular Activity (EVA) tests, improvements to the test stand were considered. Regeneration vacuum supply within the test rig was updated to more closely resemble space vacuum, theoretically increasing half-cycle timing. An improved water vapor delivery system was incorporated to dose a wider calibrated range of humidity over both short and long-duration EVAs. A range of metabolic profiles between 800 and 3000 BTU/h were tested, and the corresponding average CO 2 steady-states were compared. The effect of dosing temperature at the adsorbent bed inlet between 10 and 40 ̊C and the corresponding CO 2 removal efficiency was also evaluated.

John R Tidwell

Atomic faulting drives exceptional toughness in low thermal expansion chromium alloys

Endowing functional properties with mechanical responses in traditional metals has been a frontier topic, akin to transforming base metal into gold. Chromium and its alloys, with their functional deficiencies and limited ductility, serve as typical examples. Herein, we report a Cr96Fe4Ge1.3B1 alloy that unifies low thermal expansion (LTE, αl = 1.79 × 10-6 K-1, 200 − 315 K) with exceptional toughness (240.2 J·cm-3). The enhancement in mechanical responses is primarily attributed to layered Cr2B intermetallic precipitates, which ameliorate interfacial cohesion and simultaneously refine the grain structure. The weakened interlayer interactions within the Cr-B layers facilitate the nucleation and movement of numerous tiny stacking faults in precipitates, efficiently alleviating strain energy and resulting in marked work-hardening ability. Additionally, antiferromagnetic fluctuations in the BCC matrix contribute to the unique LTE behavior. This paves the way for the design of high-performance alloys featuring layered-symmetry precipitates.

Yu, Chengyi [University of Science and Technology

Modeling for Battery Prognostics

For any battery-powered vehicles (be it unmanned aerial vehicles, small passenger aircraft, or assets in exoplanetary operations) to operate at maximum efficiency and reliability, it is critical to monitor battery health as well performance and to predict end of discharge (EOD) and end of useful life (EOL). To fulfil these needs, it is important to capture the battery's inherent characteristics as well as operational knowledge in the form of models that can be used by monitoring, diagnostic, and prognostic algorithms. Several battery modeling methodologies have been developed in last few years as the understanding of underlying electrochemical mechanics has been advancing. The models can generally be classified as empirical models, electrochemical engineering models, multi-physics models, and molecular/atomist. Empirical models are based on fitting certain functions to past experimental data, without making use of any physicochemical principles. Electrical circuit equivalent models are an example of such empirical models. Electrochemical engineering models are typically continuum models that include electrochemical kinetics and transport phenomena. Each model has its advantages and disadvantages. The former type of model has the advantage of being computationally efficient, but has limited accuracy and robustness, due to the approximations used in developed model, and as a result of such approximations, cannot represent aging well. The latter type of model has the advantage of being very accurate, but is often computationally inefficient, having to solve complex sets of partial differential equations, and thus not suited well for online prognostic applications. In addition both multi-physics and atomist models are computationally expensive hence are even less suited to online application An electrochemistry-based model of Li-ion batteries has been developed, that captures crucial electrochemical processes, captures effects of aging, is computationally efficient, and is of suitable accuracy for reliable EOD prediction in a variety of operational profiles. The model can be considered an electrochemical engineering model, but unlike most such models found in the literature, certain approximations are done that allow to retain computational efficiency for online implementation of the model. Although the focus here is on Li-ion batteries, the model is quite general and can be applied to different chemistries through a change of model parameter values. Progress on model development, providing model validation results and EOD prediction results is being presented.

Prognostics

CFD Simulation of the Dosing Behavior within the Atomic Layer Deposition Feeding System

The effective operation of atomic layer deposition (ALD) feeding system is the premise of realizing specific ALD processes. In the present work, a detailed computational fluid dynamics (CFD) model of the feeding system has been developed and validated, which accounts for the roles of ALD valves and manifolds. A numerical simulation of the compressible fluid flow and heat/mass transfer within the feeding system was conducted. The dosing amounts and the spatiotemporal distributions of the precursors can be accurately predicted using the CFD model, as validated by experimental results. Different precursors, operating conditions, and structures of the feeding system were simulated and analyzed to examine the operating flexibility of the feeding system. The simulation results can be adopted as the upstream boundary conditions for simulations of the ALD process in the reaction chamber. The substrate-scale simulation indicates that the effect of the feeding system on the film deposition is highly related to the surface kinetics of ALD. The present work can serve as a guide for the development and optimization of different ALD-based processes via proper operation and even the design of the feeding system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

Multilevel Conditional Disturbance Rejection Control for Satellite Attitude Tracking

Recently, the conditional disturbance rejection controller (CDRC) was proposed to improve control performance by leveraging disturbances that have beneficial effects on the system. However, it only considers disturbances acting on the state variable directly influenced by them. Although fast convergence of this state can be achieved with the CDRC, it may unintentionally affect the convergence of the output (i.e., the primary state). Here, in this article, a multilevel CDRC is proposed to enhance satellite attitude control performance by accounting for the effect of disturbances on both attitude (output) and angular velocity. The extended state observer is employed to estimate the lumped disturbance, including the modeling errors and external disturbances. Then, a backstepping-based controller with the multilevel disturbance rejection law (ML-DRL) is designed for attitude tracking. The ML-DRL is developed to improve the control performance by using a disturbance with a damping effect on both attitude and velocity. Faster convergence of attitude and velocity can be achieved by conditionally compensating for the disturbance. The stability of the proposed control method is analyzed by demonstrating that the errors are bounded as time tends to infinity. The attitude control performance of the proposed method is evaluated through numerical examples conducted using the MATLAB/Simulink Multibody tool.

Active disturbance rejection control (ADRC)

High-Order Methods in NASA’s Next Generation of Computational Fluid Dynamics Tools

The missions of the National Aeronautics and Space Administration (NASA) routinely produce unique requirements and challenges for development and application of Computational Fluid Dynamics (CFD) methods. NASA presently embodies four distinct Mission Directorates: Aeronautics Research, Exploration Systems, Science, and Space Operations. These missions generate requirements for systems that operate in a wide variety of environments. They range from the high-speed flight of aerodynamically optimized vehicles operating in the earth’s atmosphere to spacecraft designed for missions that don’t favor aerodynamic optimization, some operating in the atmosphere of planets and planetary moons such as Mars and Venus or Saturn’s moon Titan. Systems supporting these vehicles, such as rocket and jet propulsion, reaction control systems, fluid and thermal transfer systems, etc. can also generate their own unique set of flow phenomena that challenge today’s CFD methodology. Through the NASA Engineering and Safety Center (NESC), NASA annually conducts state-of-the-discipline assessments in fifteen distinct engineering disciplines. These assessments are performed by the NASA Technical Fellows that lead Technical Discipline Teams (TDT) of recognized experts in these fifteen areas. In the Aerosciences discipline, three topics have been identified as the top challenges for the discipline: aero-plume interaction prediction, unsteady separated flows, and aerothermodynamic prediction. These challenge areas are defined by the Agency’s high-risk projects and problems on which the NESC is requested to perform independent testing, analysis, and assessments. When viewed as a whole, these tests, analyses, and assessments provide a clear view of the recurring technical challenges facing Agency engineers and researchers and can be used to guide future research and technology development. The present state-of-the-art in the application of CFD at NASA is the use of Reynolds-Averaged Navier- Stokes (RANS) solvers, primarily executed in a steady-state mode of operation. In isolated cases, Unsteady RANS (URANS) solvers have been employed when steady RANS solutions produce poorly converging or oscillating results or in cases, such as aeroelastic analysis, which require unsteady aerodynamic simulation. For most traditional external and internal aerodynamic flows, structured overset grids or unstructured grids are employed to minimize geometric modeling and grid generation times. Grid adaptation, primarily as a series of coarse-grain intermediate processing steps is also seeing use on particularly complex flow problems and configurations. In the case of aerothermodynamic flows, engineers have been forced to continue to employ structured grid techniques as the present unstructured grid methodology has proven inadequate in the prediction of surface heating. In the area of aero-plume interaction modeling, two-gas, frozen chemistry simulation is generally the state-of-the- art, with some production solvers capable of predicting flows with only a single gas component. Prediction of flows falling into the afore-mentioned top Aerosciences technical challenges have severely stressed the present state-of-the-art in CFD prediction and for some problems, such as unsteady separated flows and aero-plume interaction cases, engineers have begun employing Large Eddy Simulation (LES) and Hybrid RANS/LES techniques. In some isolated aero-propulsion interaction cases, chemically reacting flow simulations have been applied. These methods are highly evolutionary and engineers have little experience in their application, so they cannot be heavily relied upon in today’s application environment. Therefore, this leads one to muse over which numerical technologies will be included in the CFD tools that will be employed 30 years in the future. This presentation will describe specific technical problems that have stressed NASA’s traditional CFD methods to their breaking point and will link these issues to the Agency’s top Aerosciences technical challenges. The discussion will then shift to the characteristics of future CFD solvers that will be required to attack these challenges and how these characteristics differ from the present state-of-the- art. High-order methods certainly appear to have a place in the development of future CFD tools and some of the physical characteristics of our most challenging problems suggest that high-order methods are the only way to effectively solve them. But there are some relatively severe implementation issues that face these methods, particularly in the area of general applicability and robust operation as an engineering tool. Desired characteristics of next-generation CFD solvers will be discussed and the author’s view of which emerging numerical technologies might be employed to address these attributes will also be presented

David M Schuster

Advances in Entry Systems Modeling: Modeling Summer Visit 2025 Technical Reports

This NASA Technical Memorandum presents the collective work carried out during the Modeling Summer Visit (MSV) 2025, a collaborative research program between NASA Ames Research Center and three leading European research institutions: the University of Bordeaux, the von Karman Institute (VKI), and Ecole Polytechnique. Held from July 7th to August 1st, 2025, the NASA ARC MSV program brought together 42 researchers from Europe to work alongside NASA staff on 22 projects in entry systems modeling. The program was organized by Bruno Dias, Sergio Fraile-Izquierdo, and Jeremie Meurisse from Analytical Mechanics Associates, Inc., in collaboration with Marc Massot (Ecole Polytechnique), Jean Lachaud (University of Bordeaux), and Thierry Magin (VKI). The MSV program traces its origins to an initiative started by Nagi N. Mansour in 2018, with this 2025 edition representing a consolidation and expansion of that original vision. The projects were organized into six thematic groups: Porous Media, Numerical Methods, Plasma/Rarefied Methods, Experimental Reconstruction, Flow and Material Coupling Methods, and Transition and High-Speed Flows. Throughout the four-week program, participants presented their work in three public presentations (July 7th, July 21st, and August 1st), allowing for continuous feedback and discussion of their evolving research.

Plasmas

RG-CAT: Detection pipeline and catalogue of radio galaxies in the EMU pilot survey

Abstract We present source detection and catalogue construction pipelines to build the first catalogue of radio galaxies from the 270$\rm deg^2$pilot survey of the Evolutionary Map of the Universe (EMU-PS) conducted with the Australian Square Kilometre Array Pathfinder (ASKAP) telescope. The detection pipeline uses Gal-DINO computer vision networks (Gupta et al. 2024, PASA, 41, e001) to predict the categories of radio morphology and bounding boxes for radio sources, as well as their potential infrared host positions. The Gal-DINO network is trained and evaluated on approximately 5 000 visually inspected radio galaxies and their infrared hosts, encompassing both compact and extended radio morphologies. We find that the Intersection over Union (IoU) for the predicted and ground-truth bounding boxes is larger than 0.5 for 99% of the radio sources, and 98% of predicted host positions are within$3^{\prime \prime}$of the ground-truth infrared host in the evaluation set. The catalogue construction pipeline uses the predictions of the trained network on the radio and infrared image cutouts based on the catalogue of radio components identified using theSelavysource finder algorithm. Confidence scores of the predictions are then used to prioritiseSelavycomponents with higher scores and incorporate them first into the catalogue. This results in identifications for a total of 211 625 radio sources, with 201 211 classified as compact and unresolved. The remaining 10 414 are categorised as extended radio morphologies, including 582 FR-I, 5 602 FR-II, 1 494 FR-x (uncertain whether FR-I or FR-II), 2 375 R (single-peak resolved) radio galaxies, and 361 with peculiar and other rare morphologies. Each source in the catalogue includes a confidence score. We cross-match the radio sources in the catalogue with the infrared and optical catalogues, finding infrared cross-matches for 73% and photometric redshifts for 36% of the radio galaxies. The EMU-PS catalogue and the detection pipelines presented here will be used towards constructing catalogues for the main EMU survey covering the full southern sky.

Astronomy & Astrophysics

Dynamics and Structure of Weakly-Strained Flames In Normal- and Micro-Gravity

Strained laminar flames have been systematically studied, as the understanding of their structure and dynamic behavior is of relevance to turbulent combustion. Most of these studies have been conducted in opposed-jet, stagnation-type flow configurations. Flame studies in stagnation flows also allow for the determination of fundamental flame properties such as laminar flame speeds and extinction states that can be also conveniently modeled in detail. Studies at high strain rates are important in quantifying and understanding the response of vigorously-burning flames under conditions in which the transport time scales become comparable to the chemical time scales. Studies of weakly-strained flames can be of particular interest for all stoichiometries. For example, the laminar flame speeds for any equivalence ratio, phi, can be accurately determined by using the counterflow technique only if measurements are obtained at very low strain rates. Furthermore, near-limit flames can be only stabilized by weak strain rates. Previous studies have shown that weakly-burning flames are particularly sensitive to chain mechanisms, thermal radiation, and unsteadiness. The stabilization and study of weakly-strained flames is complicated by the presence of buoyancy that can render the flames unstable to the point of extinction. Such instabilities are caused either by the induced natural convection or the fact that higher density fluid can find itself on top of a lower density fluid. Thus, the use of microgravity (mu-g) becomes essential in order to provide meaningful insight into this important combustion regime. In view of the foregoing considerations, the main objectives of the program are to: (1) Experimentally determine the laminar flame speed at near-zero strain rates; (2) Experimentally determine the extinction limits of near-limit flames; (3) Experimentally determine the response of near-limit flames to unsteadiness and heat loss; (4) Introduce Digital Particle Image Velocimetry (DPIV) in microgravity; (5) Conduct detailed numerical simulations of the experiments; and (6) Provide physical insight into the underlying physico-chemical mechanisms.

Hai Zhang

Timing the Flames: Geostationary Satellite Detection of Diurnally Shifting Stubble Burning in Northwestern India

Post-monsoon open-field stubble burning in northwestern (NW) India—a key agricultural region known as the “breadbasket”—is a longstanding practice used to clear fields. Satellite observations spanning over two decades have revealed significant upward trends in crop production, vegetative greenness, and the frequency of post-harvest fires, with this last contributing to hazardous air quality during the peak burning season (mid-October to mid-November). Since 2022, thermal anomaly data from Aqua-MODIS and SNPP-VIIRS sensors have shown a sharp decline in reported fire events—an observation that contrasts starkly with the concurrent rise in regional aerosol loading detected from space. This apparent discrepancy became particularly pronounced in 2024–2025, prompting a closer examination using high-temporal-resolution imagery from the Advanced Meteorological Imager (AMI) on the geostationary satellite GEO-KOMPSAT-2A. These observations revealed a clear spike in fire-related signals occurring around and after 4:00 p.m. local time, i.e., outside the typical noon to 2:00 p.m. detection window of the MODIS and VIIRS. A fire detection algorithm exploiting the fire-sensitive shortwave-infrared 3.8 μm signal and its contrast to 11.2 μm infrared observations is designed to adopt AMI observations and applied to its multi-year observations (2019–2025). The resulting fire dataset unambiguously shows a gradual shift in stubble burning activity toward the late afternoon hours beginning in 2022 which is underreported by polar-orbiting satellites. The orbital drift of NASA’s MODIS sensor on the Aqua platform allows detection of some of the gradually shifting fires during afternoon hours, but the MODIS still misses a large number of fires occurring around and after 4 p.m. The AMI’s relatively coarse spatial resolution (~4 km), a consequence of its slant viewing geometry over NW India, imposes inherent limitations on quantifying the full extent of fire occurrences. The operational air quality forecasting models currently assimilate satellite fire detections predominantly captured during early afternoon overpasses of the MODIS and VIIRS. The temporal shift in fire activity complicates such forecast, leading to a substantial underestimation of emissions. Intense stubble burning and the resulting air pollution highlight the need for effective crop residue management practices for mitigating the frequency of open biomass burning and thereby reducing episodic degradation of air quality and its associated public health and economic impacts.

post-monsoon stubble burning; northwestern India;

Laser-heated diamond anvil cell synthesis and recovery of metastable MnSb 2 and YbZn 2 for post-synthesis transport studies

The creation and exploration of new materials under extreme pressure–temperature conditions has become increasingly reliant on laser-heated diamond anvil cell (LHDAC) techniques, which provide direct access to previously unexplored regions of multinary phase diagrams. Whereas numerous high-pressure phases have been identified in situ, systematic recovery and post-synthesis physical property characterization of these materials remain significant challenges. In this work, we describe the setup and implementation of an LHDAC-based synthesis and recovery workflow and demonstrate its application to metastable MnSb 2 and YbZn 2 phases. Synchrotron x-ray diffraction and spatial mapping confirm dominant formation of the targeted phases, whereas laboratory-based refinement quantifies phase fractions despite intrinsic microstrain and minor secondary phases. High-pressure transport measurements on recovered samples reveal pressure-tunable electronic instabilities in both systems. In MnSb 2 , pressure suppresses two high-temperature magnetic ordering anomalies, observed in transport, by ∼5 GPa and, for higher pressures, induces a new low-temperature feature that increases with further pressure increase. In hexagonal high-pressure YbZn 2 , an electronic reconstruction emerges at ∼11 GPa, characterized by semiconducting-like behavior from ∼30 to 300 K and a broad low-temperature coherence crossover near 30 K. Our results establish LHDAC synthesis not only as a structural discovery tool but also as an experimental platform for investigating correlated quantum states stabilized far from equilibrium thermodynamic conditions.

Huyan, S. [Iowa State University, Ames, IA (United

Design and fabrication of ion traps for low RF power dissipation

Large surface-electrode ion traps with multiple trapping regions and junctions are a natural approach to scaling trapped ion quantum computers, supporting the connectivity and ion counts necessary for complex quantum algorithms. However, a major hurdle in this scaling is on-chip power dissipation from the applied RF voltage, which increases at a rate between linear and cubic relative to trap size, depending on whether the losses are dielectric or Ohmic. Here, we present two versions of a trap with features designed to reduce both types of RF power dissipation. The first variant contains a raised RF electrode that increases the electrode–ground distance to reduce capacitance. Different DC voltage sources are demonstrated on this trap to show that technical noise before the filter remains the dominant source of voltage noise and therefore motional heating. The second variant additionally includes a method for removing dielectric from beneath the RF electrode to further reduce dielectric losses. These traps were demonstrated at room temperature with 40 Ca + ions. In conclusion, the similar heating rates and heating rate axial frequency dependencies between 2.4 and 3.0 MHz illustrate that this dielectric modification is not detrimental to trap performance.

Sterk, J. D. [Sandia National Laboratories (SNL-NM

Analysis of Impact Induced Damage and its Effect on Structural Integrity of Space Flight Composite Overwrapped Pressure Vessels

The objective of this research work has been to provide analytical background and support to the ongoing experimental program at NASA, White Sands Test Facility, involving testing composite overwrapped pressure vessels (COPV) for impact damage and cyclic pressurization. Preliminary theoretical basis, including the governing equations for a shallow shell subjected to internal pressure, has been established. Effects of the Griffith type cracks on the structural integrity of the cylindrical vessel were evaluated by methods of Fracture Mechanics. The results indicate that the effective mass of the pressure vessel is an important factor influencing the response to impact events. We also have found that the material properties of the target, contained in the constitutive equations of the composite attached to the Aluminum liner, dominate the impact event in the low velocity range, the material properties become less important, while the target mass distribution and the impactor mass become more significant as the velocity of the impactor increases. Therefore, at high-velocity impact it is not only the kinetic energy of the impactor but also its mass which has a significant effect on the dynamics of the event, and consequently on the induced damage. This work also suggests a methodology for an assessment of the rate of loading effects on the degradation of the material toughness associated with a high-velocity impact where the rate effects become significant. To model the rate dependence of the material response a viscoelastic-plastic constitutive equations were assumed, and on this basis predictions are made regarding the rate dependent material resistance curve. Other dynamic phenomena associated with the impact event have been treated in the framework of the Computational Mechanics using the courtesy of Prof. P. Guebelle and his graduate student at University of Illinois at Urbana-Champaign who have an access to a super-fast computer located on their campus. Finally, the guidelines for a follow-up research program are provided in the body of this report. They address three major areas: theoretical research, numerical studies, and further experimental work.

Michael P Wnuk

Tree tensor network hierarchical equations of motion based on time-dependent variational principle for efficient open quantum dynamics in structured thermal environments

In this work, we introduce an efficient method, TTN-HEOM, for exactly calculating the open quantum dynamics for driven quantum systems interacting with highly structured bosonic baths by combining the tree tensor network (TTN) decomposition scheme with the bexcitonic generalization of the numerically exact hierarchical equations of motion (HEOM). The method yields a series of quantum master equations for all core tensors in the TTN that efficiently and accurately capture the open quantum dynamics for non-Markovian environments to all orders in the system–bath interaction. These master equations are constructed based on the time-dependent Dirac–Frenkel variational principle, which isolates the optimal dynamics for the core tensors given the TTN ansatz. The dynamics converges to the HEOM when increasing the rank of the core tensors, a limit in which the TTN ansatz becomes exact. We introduce TENSO, tensor equations for non-Markovian structured open systems, as a general-purpose Python code to propagate the TTN-HEOM dynamics. We implement three general propagators for the coupled master equations: two fixed-rank methods that require a constant memory footprint during the dynamics and one adaptive-rank method with a variable memory footprint controlled by the target level of computational error. We exemplify the utility of these methods by simulating a two-level system coupled to a structured bath containing one Drude–Lorentz component and eight Brownian oscillators, which is beyond what can presently be computed using the standard HEOM. Our results show that the TTN-HEOM is capable of simulating both dephasing and relaxation dynamics of driven quantum systems interacting with structured baths, even those of chemical complexity, with an affordable computational cost.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

A GPU-based Approach for Turbomachinery Application

Accurate modeling of coolant airflows, which form protective films over turbine blades, is essential for designing fuel-efficient and environmentally sustainable gas turbine engines. Excessive coolant reduces thermal efficiency, while insufficient coolant leads to blade overheating and causes damage. Therefore, precise prediction of flow field interactions with cooling air is critical for optimizing turbine performance. This study numerically investigates the cooling effectiveness of purge and film cooling flows within a high-pressure turbine (HPT) rotor using Large Eddy Simulation (LES). The study utilizes NASA Glenn Research Center’s Glenn-HT solver. The simulation models ethe conditions of the Penn State University START rotating rig. A high-fidelity structured mesh comprising up to 800 million cells is employed to resolve high-Reynolds number flow (Re ≈ 350,000) and to capture intricate secondary flow structures, including tip leakage and purge-induced vortices. Film cooling effectiveness computations are highly sensitive to boundary conditions at the cooling holes and to grid resolution. Even with well-resolved grids and included plena, strong mixing challenges traditional eddy viscosity models. To address this, a simplified configuration is simulated: a truncated row of shaped holes on the suction side near the leading edge and a row on the pressure side, both fed from internal plena while the purge slot and tip clearance are also modeled. Two isothermal LES cases are conducted at two distinct wall temperatures, which yield the adiabatic wall temperature and the heat transfer coefficient. The definition and means of computation of the effectiveness is discussed in this paper. The simulations reveal detailed three-dimensional unsteady flow features, including coherent vortical structures and secondary flows originating from the purge cavity. Film cooling effectiveness and Nusselt number distributions are presented for both the blade surface and tip, highlighting regions of elevated heat transfer and complex thermal behavior. These findings underscore the importance of high-resolution LES and realistic boundary conditions in capturing the dynamics of purge and film cooling, offering valuable insights for improving turbine blade design and thermal management strategies.

Gas Turbine