Search NASA⌕ Search

SEARCH · Search NASA

Results for “action–angle framework”

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.

394 records · Page 20
All sourcesNASA NTRSDOE OSTI
NASA NTRS2009

General Purpose Data-Driven Monitoring for Space Operations

As modern space propulsion and exploration systems improve in capability and efficiency, their designs are becoming increasingly sophisticated and complex. Determining the health state of these systems, using traditional parameter limit checking, model-based, or rule-based methods, is becoming more difficult as the number of sensors and component interactions grow. Data-driven monitoring techniques have been developed to address these issues by analyzing system operations data to automatically characterize normal system behavior. System health can be monitored by comparing real-time operating data with these nominal characterizations, providing detection of anomalous data signatures indicative of system faults or failures. Data-driven techniques have a number of advantages over other methods for monitoring complex space vehicles. Unlike model-based systems, the developer does not need to understand or encode the internal operation of the system. The knowledge required to monitor the system is automatically derived from archived data from system operation. Unlike rule-based systems, data-driven systems do not require system analysts to define nominal relationships among sensors. Analysts can and often do determine these relationships for a system with few sensors; it is more difficult to analytically determine the nominal relationship among a large number of sensors. Data-driven techniques are not limited to low-dimensional spaces and work as effectively with dozens of parameters as they do with a few. Knowledge bases formed by data-driven techniques are also easy to update. As the operating envelope of the monitored system is expanded, data-driven techniques can be quickly retrained to incorporate the new behavior into the knowledge base. The expertise and time-consuming process of updating a model or rule base to maintain consistency with the new operation is not required. The Inductive Monitoring System (IMS) is a data-driven system health monitoring software tool that has been successfully applied to several aerospace applications. IMS uses a data mining technique called clustering to analyze archived system data and characterize normal interactions between parameters. This characterization, or model, of nominal operation is stored in a knowledge base that can be used for real-time system monitoring or analysis of archived events. System data is compared with the nominal IMS model to produce a measure of how well current system behavior matches the normal behavior defined by the training data. Significant deviations from the nominal system model can provide alerts to system malfunctions or precursors of significant failures. The scope of IMS based data-driven monitoring applications continues to expand with current development activities. Successful IMS deployment in the International Space Station (ISS) flight control room to monitor ISS attitude control systems has led to applications in other ISS flight control disciplines, such as thermal control. It has also generated interest in data-driven monitoring capability for Constellation, NASA's program to replace the Space Shuttle with new launch vehicles and spacecraft capable of returning astronauts to the moon, and then on to Mars. Several projects are currently underway to evaluate and mature the IMS technology and complementary tools for use in the Constellation program. These include an experiment on board the Air Force TacSat-3 satellite, and ground systems monitoring for NASA's Ares I-X and Ares I launch vehicles. The TacSat-3 Vehicle System Management (TVSM) project is a software experiment to integrate fault and anomaly detection algorithms and diagnosis tools with executive and adaptive planning functions contained in the flight software on-board the Air Force Research Laboratory TacSat-3 satellite. The TVSM software package will be uploaded after launch to monitor spacecraft subsystems such as power and guidance, navigation, and control (GN&C). It will analyze data in real-time to demonstrate detection of faults and unusual conditions, diagnose problems, and react to threats to spacecraft health and mission goals. The experiment will demonstrate the feasibility and effectiveness of integrated system health management (ISHM) technologies with both ground and on-board experiments. Initially, the TVSM software will run open loop, providing system health information and recommendations to ground operators, without automatically performing fault-mitigating corrective actions. After the end of the satellite's mission, closed loop tests combining TVSM monitoring and diagnosis with reactive capabilities by the flight software will be performed. In addition to monitoring for long periods of actual operation, the experiment will include fault injection into TacSat-3 data as well as commanded operations to test and evaluate automatic ISHM monitoring and recovery under controlled conditions.

Satellites↗
DOE OSTI2026

Relativistic Magnetohydrodynamic Simulations of Giant Magnetar Bursts

Gradual crustal deformation can generate strongly twisted magnetic fields around magnetars, potentially triggering giant flares with total energies exceeding 10 44 erg. In this letter, we present the first relativistic magnetohydrodynamic simulation of a surface shear-driven magnetar eruption, capturing reconnection-driven plasma heating, the ejection of relativistically hot plasma, and the formation of a hot fireball confined within the inner magnetosphere. We find that magnetic reconnection in the equatorial current sheet launches a hot trailing outflow capable of powering the initial spike observed in giant flares, while simultaneously leaving behind a thermally stratified fireball with sufficient thermal energy to produce the pulsating, decaying tail. Together, these features provide a self-consistent physical framework for understanding the observed energetics of magnetar giant flares. The eruption also expels a magnetically dominated giant plasmoid carrying up to ∼9% of the magnetosphere’s total magnetic energy. Furthermore, our simulation demonstrates how the plasmoid drives the formation of a blast wave—an important ingredient in models linking magnetar eruptions to fast radio bursts.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

First multi-institutional systematic comparison of the neutron ambient dose equivalent produced by proton therapy systems

Objective. Isochronous cyclotrons, synchrocyclotrons, and synchrotrons are used to accelerate protons for proton therapy. An accurate measurement of neutron doses generated by these accelerators and associated delivery systems and its clinical relevance requires systematic protocols and proper neutron dosimetry for a meaningful assessment. We present the first comprehensive comparison of neutron ambient dose equivalent (H*(10)) produced by clinically operational proton therapy systems. Approach. Treatment plans with 10 cm modulation-depth and ranges of 10 cm (R10M10) and 25 cm (R25M10) were created to cover a 10 × 10 × 10 cm 3 water target. The pencil beam scanning proton therapy machines studied were: two gantry-mounted synchrocyclotrons (Hyperscan, Mevion, half-gantry), two isochronous cyclotrons (ProBeam, Varian, full-gantry), one isochronous cyclotron (Proteus, IBA, full-gantry), and two synchrotrons (PROBEAT, Hitachi, full- and half-gantry). Proton beams were delivered to 30 × 30 × 40 cm 3 plastic water phantoms. WENDI-II and LUPIN-BF3-NP neutron rem-meters were positioned at three angles (0°, 45°, 90°) relative to the beam direction to measure the neutron H*(10) at distances between 50–300 cm from the isocenter. Main results. H*(10) showed dependence on beam energy, machine type, and measurement location. The highest reading was for the gantry-mounted synchrocyclotron, whereas other systems produced approximately comparable neutron doses. In all cases, the H*(10) reduced with distance from the isocenter. The H*(10) drop at 2 m distance compared to that at 0.5 m was a factor of ∼5 for the gantry-mounted synchrocyclotron whereas in other systems the decrease was a factor of 10. The WENDI-II device suffered from dead-time-associated under-estimation of the dose by a factor of ∼2–3 under the synchrocyclotron beam due to its high dose-per-pulse. However, WENDI-II and LUPIN-BF3-NP results were within reasonable agreement in isochronous cyclotron and synchrotron beams, indicating that both devices are suitable for those systems. Significance. Neutron H*(10) is dependent on various parameters including beam energy, measurement location, as well as machine design. Caution must be exercised in choosing the appropriate neutron-dose-measurement device to be used for low-duty-factor, particularly in high-instantaneous-rate proton delivery systems. By delivering the same volumetric proton dose across different machines, this work provides a benchmark for inter-system comparisons and serves as a foundation for future studies.

LUPIN↗
NASA NTRS2026

Modeling Large Dust Aerosols in the Community Earth System Model Version 2 (CESM2)

Dust aerosols have a wide size distribution from less than 0.1 to over 100 μm and dominate Earth's atmospheric aerosol mass. However, most Earth system models (ESMs) inadequately represent dust aerosols larger than 10 μm in diameter, limiting the accuracy of the simulated dust cycle and climate impacts. Here, we introduce a new modeling framework that captures the full observed size distribution of dust aerosols, incorporating recent advances into a mineral-resolved version of the Community ESM, while addressing known issues in previous versions. Comprehensive evaluation against diverse observations of bulk dust and component minerals demonstrates that the model reproduces the observed dust cycle across particle sizes. Incorporating the previously unrepresented large-dust fractions substantially alters dust budget estimates, highlighting potential changes in simulated climate impacts and underscoring the importance of comprehensive size-resolved dust modeling. Despite these advancements, uncertainties persist. Our results indicate that a size-dependent reduction in settling velocity is required to reproduce the observed dust size distribution downwind of source regions. Specifically, in the new model, the gravitational settling velocity of dust particles larger than 10 μm in diameter must be reduced by as much as 85% to achieve agreement with observations. This empirical reduction serves as a constraint on physics-based models of dust settling. Future developments should address misrepresented physical processes that hinder accurate modeling of the large dust aerosol transport. Expanding observational data sets covering the full-size distribution is also essential to better constrain the dust cycle and improve the representation of dust optical properties and climate effects.

mineral dust↗
DOE OSTI2026

A Tutorial on Bayesian analysis of linear shock compression data

Gas gun and other shock compression experiments often produce shock wave velocity measurements that are linearly associated with particle velocity. Traditionally, this empirical relationship is quantified with a single Hugoniot curve that is estimated using least squares regression. However, for downstream modeling and simulation tasks, it is often more useful to have multiple Hugoniot curves in the pressure–volume plane that are consistent with the data. We employ Bayesian uncertainty quantification methods as a framework for propagating measurement uncertainty through to model parameters and predictions. Specifically, this Tutorial shows how to sample multiple Hugoniot curves in the pressure–volume plane that are consistent with the shock wave-particle velocity measurements in a two-step Bayesian approach. First, we obtain an analytical expression for the posterior distribution of the linear model parameters using Bayesian linear regression. Second, we propagate samples from the posterior distribution through the Rankine–Hugoniot equations to yield Hugoniot curves in the pressure–volume plane. The procedure is demonstrated with publicly available data on argon, copper, and nickel, and compared against bootstrapping and linear regression. The Bayesian procedure is shown to be interpretable, computationally inexpensive, and less sensitive than an alternative bootstrapping approach to the removal of the point in the copper dataset that has the largest particle velocity. As a Tutorial on Bayesian methodology for the shock compression community, we provide several derivations and explanations that make this paper self-contained, and make all code and data available at github.com/llnl/BALSCD.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Mesoscale Magnetostructural Phase Separation in Fe‐deficient Fe 5 GeTe 2

Two-dimensional van der Waals ferromagnet Fe 5-x GeTe 2 (F5GT) is promising for spintronic applications due to its high Curie temperature, layered structure, and ability to host complex magnetic textures. However, the origin of its sample-dependent magnetic anisotropy remains unclear, hindering control of its magnetic behavior. Here, we use spatially resolved cryogenic scanning transmission electron microscopy (STEM) to correlatively map magnetism, lattice structure, and chemistry across atomic-to-micron scales. We reveal that only mesoscale, not nanoscale, inclusions of a Fe-deficient secondary phase significantly modify magnetic behavior, establishing a previously unrecognized critical length scale. This phase separation, induced by quenching, leads to in-plane magnetic anisotropy, while slow cooling confines separation to a few nanometers and preserves out-of-plane anisotropy. These findings reconcile prior inconsistencies and establish a predictive framework for tuning magnetism in F5GT through thermal processing, with broader implications for controlling anisotropy in other two-dimensional magnetic materials.

2D ferromagnets↗
DOE OSTI2026

Linking international technical specifications for acoustic characterization of marine energy converter sounds with environmental compliance criteria

As new ocean energy technologies emerge and are deployed for testing and operations, sound emissions are a potential concern for environmental effects to marine life. Consistent acoustic measurement and data analysis methods can help promote comparisons of technologies and transferability between project sites. In 2022, acoustic emissions from a prototype scale wave energy converter (WEC) were characterized for a range of environmental conditions and power generation states in the coastal waters off southern California using a set of international technical specifications. Results from the international technical specification analyses were applied to United States regulatory threshold criteria for acoustic impacts to marine mammals and examined in the context of European underwater noise monitoring guidelines. Weighted 24 hour cumulative sound exposure levels SEL24h calculated from the highest power generation state WEC sound pressure levels were often more than 20 dB below threshold criteria for temporary threshold shifts in five relevant marine mammal hearing groups. Following European Union recommendations for analyses and reporting, WEC sound characterization in third octave bands centered at 63 Hz and 125 Hz show clear spatial decay of WEC-generated noise, with more pronounced attenuation at 63 Hz, and a less marked but still detectable gradient at 125 Hz, collectively suggesting a relatively confined acoustic footprint under the observed conditions. The value of the international technical specification approach is highlighted by the isolation of WEC sounds from the surrounding soundscape. This allows for a robust characterization of acoustic emissions through a range of device power generation and sea states. Furthermore, in threshold-based regulatory contexts like the U.S., this facilitates direct evaluation of source contributions, while in broader monitoring frameworks used in the E.U. it provides a reproducible foundation for assessing the contribution of emerging ocean energy technologies to the underwater acoustic environment.

Haxel, Joseph H. (ORCID:0000000273864761)↗
DOE OSTI2026

Reaction Optimization for Enzymatic Deconstruction of Industrially Relevant Nylon Composites

Plastics such as polyamides (PAs) possess unique physicochemical properties that make them indispensable in modern society. However, their energy‐intensive production and challenging end‐of‐life management highlight the urgent need for efficient recycling or remanufacturing solutions. Enzymatic depolymerization offers a promising route toward circular recycling, yet remains constrained by limited enzyme characterization, lack of validation under industrially relevant conditions and substrates, and overall performance. Here, we optimized the reaction conditions for three recently discovered nylon‐degrading enzymes. One of them, Nyl12, achieved product titers with PA6 and PA66 that exceed previously reported values, without enzyme engineering or substrate pretreatment. We further demonstrated the scalability of the process and its application to complex PA‐based materials used in microelectronic components. Analysis of substrate features, including surface area and particle size, revealed key parameters governing enzymatic activity and provided a framework for future pretreatment and process optimization efforts. In combination, these efforts provide a new benchmark for enzymatic nylon recycling.

nylon↗
DOE OSTI2026

Refractory-based thermal energy storage for industrial process heat: one-dimensional modeling, control, and optimization

The variable and weather-dependent output of wind and solar power plants present a substantial challenge for planning and operating electricity-systems, particularly in the absence of cost-effective and dispatchable energy storage technologies. This study investigates a high-temperature, electrically heated, refractory-based thermal energy storage (RTES) system that stores electrical energy as sensible heat in dense ceramic bricks over the 950–1800 °C range. The stored heat can be discharged as a controlled hot-gas stream for industrial heating, fuel substitution in high-temperature processes, or electricity generation. The main novelty is a comprehensive modelling, control, mapping, and optimization framework that integrates one-dimensional transient gas–solid heat transfer, fan-assisted discharge, bypass-flow regulation, reheating logic, fan-power evaluation, insulation-loss assessment, and genetic-algorithm-based design optimization. The model uses feedback from outlet temperature and delivered power to regulate discharge, while a two-stage genetic algorithm optimizes brick-channel geometry, gas-flow operation, and multilayer insulation thicknesses. Storage capacities below 50 MWh and discharge powers of 5–30 MW are analyzed to evaluate hold time, thermal delivery, fan-power penalty, heat loss, state-of-charge evolution, and indicative capital cost. Results demonstrate that optimized and well-insulated refractory-based thermal energy storage units can provide stable, efficient, and repeatable heat delivery over multiple discharge cycles. The generated performance and cost maps support modular refractory thermal energy storage as a practical option for large-scale integration of wind and solar generation and for high-temperature industrial process heat.

25 ENERGY STORAGE↗
DOE OSTI2024

The ReSWARM microgravity flight experiments: Planning, control, and model estimation for on‐orbit close proximity operations

Abstract On‐orbit close proximity operations involve robotic spacecraft maneuvering and making decisions for a growing number of mission scenarios demanding autonomy, including on‐orbit assembly, repair, and astronaut assistance. Of these scenarios, on‐orbit assembly is an enabling technology that will allow large space structures to be built in situ, using smaller building block modules. However, like many of these scenarios, robotic on‐orbit assembly involves several technical hurdles, such as changing system models. For instance, grappled modules moved by a free‐flying “assembler” robot can cause significant changes in the combined system inertia, which have cascading impacts on motion planning and control portions of the autonomy stack. Further, on‐orbit assembly and other scenarios require collision‐avoiding motion planning, particularly when operating in a “construction site” scenario of multiple assembler robots and structures. Multiple key technologies that address these complicating factors for autonomous microgravity close proximity operations are detailed in this work, in particular: (1) application of global long‐horizon planning, accomplished using offline and online sampling‐based planner options that consider the system dynamics; (2) adaptation of the recently proposed RATTLE information‐aware planning framework for on‐orbit reconfiguration model learning; and (3) connection with robust control tools to provide low‐level control robustness using current system knowledge. These approaches were demonstrated for an autonomous on‐orbit assembly use case by the RElative Satellite sWarming and Robotic Maneuvering (ReSWARM) experiments using NASA's Astrobee robots on the International Space Station. Results of the ReSWARM experiments are provided along with significant operational and implementation detail discussing the practicalities of hardware implementation and unique aspects of working with the Astrobee free‐flyer robots in microgravity. ReSWARM provides a base set of planning and control tools for robotic close proximity operations, demonstrates them in microgravity, and outlines some of the important hardware aspects that future autonomous free‐flyers will need to consider.

Robotics↗
DOE OSTI2026

How Proton Incorporation Reshapes Lattice Dynamics In BaSnO 3 ‐Type Proton Conductors

Proton conduction in acceptor-doped perovskites is fundamentally a vibronic process: mobile H + and D + do not move independently, but dynamically co-vibrate with the surrounding oxygen–metal framework. Direct experimental evidence for this behavior is presented using in situ 119⁢ 𝑆⁢𝑛 nuclear resonance vibrational spectroscopy (NRVS) on hydrated, deuterated, and dry 𝐵⁢𝑎⁢𝑆⁢𝑛 1−𝑥⁢ 𝑌 𝑥 ⁢𝑂 3−𝛿 . Hydration induces systematic redistributions in the Sn-projected phonon density of states (PDOS), including an upshift of the first spectral moment by about 0.4 meV, indicating a stiffening of the extended Sn–O network. H/D isotopic substitution leaves the Sn-projected PDOS largely unchanged, with only subtle isotope-dependent spectral reweighting, demonstrating that protonic degrees of freedom are not localized oscillators but are embedded in collective lattice modes. These results are rationalized using a classical coupled proton–phonon oscillator model that links the observed PDOS variations to changes in effective force constants and vibrational mass terms. The model captures how H + and D + participate in cooperative lattice dynamics rather than forming isolated OH/OD entities. Overall, NRVS probes proton–lattice coupling in ceramic proton conductors and quantitatively describes how protonic defects modulate host lattice dynamics to enable phonon-assisted long-range proton transport.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Kinetic control of phase evolution and defect-mediated recombination in AACVD-grown copper antimony sulphide thin films

Copper antimony sulphide (CAS) is a multinary chalcogenide semiconductor in which small deviations from stoichiometry can drive phase competition and strong defect-mediated modulation of optoelectronic properties. However, the systematic roles of copper precursor fraction and deposition time in governing phase evolution, off-stoichiometry, and recombination dynamics in aerosol-assisted chemical vapour deposition (AACVD)-grown undoped CAS thin films remain insufficiently understood. In this work, undoped CAS thin films were deposited by AACVD using Cu(dedtc)2 and Sb(dedtc)3 single-source precursors at 550 °C and a carrier gas flow rate of 150 sccm, while the Cu(dedtc)2 mole fraction (x = 0.15 – 0.55) and deposition time (1 – 2 hours) were systematically varied to probe how growth kinetics influence phase composition, microstructure, and defect-mediated optical properties without post-deposition annealing or extrinsic doping. Increasing copper precursor content drives phase evolution toward tetrahedrite-dominant CAS films at intermediate Cu(dedtc)2 mole fractions, with the film deposited at x = 0.35 exhibiting the strongest tetrahedrite character within the parameter space examined. The films are also copper-rich, antimony-poor, and sulphur-deficient, consistent with off-stoichiometric growth and intrinsic defect formation, plausibly including copper interstitials, Cu-on-Sb antisites, and sulphur vacancies. These growth-dependent compositional deviations are accompanied by tunable indirect optical bandgaps of approximately 1.60 – 2.18 eV and weak visible photoluminescence governed by defect-mediated recombination. Time-resolved photoluminescence reveals bi-exponential decay behaviour with lifetimes of approximately 0.1 – 3.6 ns, with emission dominated by slower donor–acceptor pair recombination and a smaller contribution from faster trap-assisted pathways. Collectively, these results establish an explicit kinetic process–structure–defect–property relationship for AACVD-grown CAS thin films and provide a growth–structure–property framework relevant to future optimization of CAS-based optoelectronic and energy materials.

36 MATERIALS SCIENCE↗
DOE OSTI2026

SAM Code Enhancements for Fission Product Tracking of Noble Gases and Metals in MSRs

This report documents fiscal year 2026 enhancements to the System Analysis Module (SAM) for modeling fission product transport in liquid-fueled molten salt reactors (MSRs). The work advances three principal areas: noble gas transport, noble metal deposition, and user interface improvements. The noble gas transport capability integrates drift-flux gas transport, Henry’s law two-film interphase mass transfer with pressure-based nucleation suppression, Knudsen-regime pore diffusion into porous graphite with a conjugate salt-graphite interface constraint, built-in material properties, five Sherwood-number mass transfer correlations including three derived from high-fidelity NekRS simulations, and xenon-135 reactivity feedback through SAM’s point-kinetics model. This work also presents a comprehensive verification test suite, including new analytically verified cases for pressure-dependent onset of interphase gas transfer in a stagnant vertical pipe, a postulated FLiBe-graphite Xe extraction permeator, a gravity riser with a fission-product source, and a descending pipe with gas redissolution driven by hydrostatic pressure. A machine learning framework for bubble rise velocity prediction in molten salt systems is developed and benchmarked on molten-salt and diverse aqueous bubble datasets. The best-performing fine-tuned transfer-learning networks achieve an 82% reduction in RMSE relative to the Clift correlation, and is implemented directly in SAM. The noble metal transport capability is developed, including a liquid-wall deposition model and a gas-surface flotation mechanism that transfers insoluble particles entrained by sparging gas to wetted structures. Verification tests and demonstration cases cover the surface deposition, flotation efflux, and flotation shedding. Finally, a new [SpeciesTransport] input structure replaces positional global vectors with selfcontained, order-independent, named species blocks, simplifies the specification of multiphase species and decay chains, and remains fully compatible with existing SAM input files. Together, these developments improve the physical fidelity, verification basis, and usability of SAM for system-level analyses of fissionproduct behavior in MSRs.

Mui, Travis (ORCID:0000000303736470)↗
DOE OSTI2026

Experimental observation and integrated modelling of proton-beryllium fusion in He and D plasmas at JET

Validated integrated modelling of JET ITER-like wall experiments in which fusion performance is driven by reactions between fast ions and intrinsically present metal wall impurities is presented. A steady-state L-mode plasma with dominant proton-beryllium fusion and neutron yields of up to ≈ 6·10 13 s -1 is developed in He and D, via radiofrequency heating of a H minority. The fusion drive is unambiguously confirmed by the neutral particle analyser, fast ion loss detector, and γ-ray diagnostics. Experiments are analysed via an integrated modelling framework, developed to model the two-stage proton beryllium-fusion chain and produce high-fidelity fusion product source terms. The modelling chain comprises TRANSP and JETTO for plasma core modelling, LOCUST for full orbit product tracking and collisional slowing-down, DRESS to resolve two- and three-body fusion kinematics, and MCNP for neutron transport calculations. Modelling shows that the primary 9 Be(p,n) 9 B reaction is the dominant neutron emitter at naturally present concentrations of beryllium in these experiments. The yield contribution of secondary reactions between fusion products and beryllium, 9 Be(d,n) 10 B and 9 Be(α,n) 12 C, is found to be negligible. The proton-deuteron knock-on effect in D plasmas is modelled, which is calculated to contribute ≈ 25% to the total neutron yield. For both He and D discharges the total computed neutron rates match fission chamber (FC) measurements within the combined experimental and computational uncertainty, with an average discrepancy of ≈ ± 20%. Realistic proton-beryllium neutron sources are propagated through JET’s MCNP neutron transport model which shows that 235 U FCs’ response is sensitive to p–Be source changes, with up to ≈ 10% variation compared to a D–D neutron source. We show that the high-energy tail of the fast proton minority can be studied with multi-foil neutron activation. The framework is also applied to the study of interactions between fast protons and boron impurities, of relevance to ITER. We calculate that in JET conditions a significant alpha source with DT-like energies could be generated through 11 B(p, α)2α fusion, and detected via γ-emission in secondary interactions between fast alphas and boron. The work represents an important step towards validating predictive integrated modelling capabilities for non-standard fusion reactions.

JET↗
NASA NTRS1997

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↗
DOE OSTI2026

Metal–bicarbonate ion pairing in alkaline aqueous solutions from multilevel embedded correlated wavefunction theory and molecular dynamics

In this work, we examine ion-pairing mechanisms of bicarbonates in alkaline aqueous solutions with the divalent metal ions most abundantly present in seawater, namely, Ca 2+ and Mg 2+ . We employ a rare-event enhanced sampling approach within first-principles molecular dynamics to explore regions of phase space spanning solvent-shared to contact ion pairs. Second-order Møller–Plesset perturbation theory (MP2) corrections are subsequently applied in an embedding framework (EMB) to refine the electronic structure of stationary states and associated reaction barriers along the free-energy profiles while retaining the extended solvent effects at the density functional theory (DFT) level. Ca 2+ –HCO 3 − was previously hypothesized to exist in a solvent-shared ion pair (SSHIP) by DFT studies with an endergonic contact ion pair (CIP) formation; however, our EMB-MP2 refinement of the DFT ion-pairing pathways reveals that Ca 2+ and HCO 3 − form a virtually barrier-free CIP in alkaline solutions, with even more energetic ease than the widely studied Ca–CO 3 ion pair. We find qualitative agreement between DFT and EMB-MP2 for Mg 2+ —unlike Ca 2+ , Mg 2+ refuses to shed its strong hydration shell, thereby preferring a SSHIP state with a significant activation barrier to crossover to the CIP forms—a trait reminiscent of ion pairing in Mg–CO 3 and closely related to the kinetic limitations underlying the famous subject of the dolomite problem. Our study highlights the importance of improved electronic structure descriptions of liquids, modeled as a condensed phase of matter lacking in long-range crystalline order. It also strongly suggests that Ca 2+ –HCO 3 − CIPs are likely precursors involved in prenucleation of CaCO 3 mineral formation in seawater.

Sharma, Vidushi [Princeton Plasma Physics Laborato↗
NASA NTRS2014

High Efficiency Centrifugal Compressor for Rotorcraft Applications

The report "High Efficiency Centrifugal Compressor for Rotorcraft Applications" documents the work conducted at UTRC under the NRA Contract NNC08CB03C, with cost share 2/3 NASA, and 1/3 UTRC, that has been extended to 4.5 years. The purpose of this effort was to identify key technical barriers to advancing the state-of-the-art of small centrifugal compressor stages; to delineate the measurements required to provide insight into the flow physics of the technical barriers; to design, fabricate, install, and test a state-of-the-art research compressor that is representative of the rear stage of an axial-centrifugal aero-engine; and to acquire detailed aerodynamic performance and research quality data to clarify flow physics and to establish detailed data sets for future application. The design activity centered on meeting the goal set outlined in the NASA solicitation-the design target was to increase efficiency at higher work factor, while also reducing the maximum diameter of the stage. To fit within the existing Small Engine Components Test Facility at NASA Glenn Research Center (GRC) and to facilitate component re-use, certain key design parameters were fixed by UTRC, including impeller tip diameter, impeller rotational speed, and impeller inlet hub and shroud radii. This report describes the design effort of the High Efficiency Centrifugal Compressor stage (HECC) and delineation of measurements, fabrication of the compressor, and the initial tests that were performed. A new High-Efficiency Centrifugal Compressor stage with a very challenging reduction in radius ratio was successfully designed, fabricated and installed at GRC. The testing was successful, with no mechanical problems and the running clearances were achieved without impeller rubs. Overall, measured pressure ratio of 4.68, work factor of 0.81, and at design exit corrected flow rate of 3 lbm/s met the target requirements. Polytropic efficiency of 85.5 percent and stall margin of 7.5 percent were measured at design flow rate and speed. The measured efficiency and stall margin were lower than pre-test CFD predictions by 2.4 percentage points (pt) and 4.5 pt, respectively. Initial impressions from the experimental data indicated that the loss in the efficiency and stall margin can be attributed to a design shortfall in the impeller. However, detailed investigation of experimental data and post-test CFD simulations of higher fidelity than pre-test CFD, and in particular the unsteady CFD simulations and the assessment with a wider range of turbulence models, have indicated that the loss in efficiency is most likely due to the impact of unfavorable unsteady impeller/diffuser interactions induced by diffuser vanes, an impeller/diffuser corrected flow-rate mismatch (and associated incidence levels), and, potentially, flow separation in the radial-to-axial bend. An experimental program with a vaneless diffuser is recommended to evaluate this observation. A subsequent redesign of the diffuser (and the radial-to-axial bend) is also recommended. The diffuser needs to be redesigned to eliminate the mismatching of the impeller and the diffuser, targeting a slightly higher flow capacity. Furthermore, diffuser vanes need to be adjusted to align the incidence angles, to optimize the splitter vane location (both radially and circumferentially), and to minimize the unsteady interactions with the impeller. The radial-to-axial bend needs to be redesigned to eliminate, or at least minimize, the flow separation at the inner wall, and its impact on the flow in the diffuser upstream. Lessons were also learned in terms of CFD methodology and the importance of unsteady CFD simulations for centrifugal compressors was highlighted. Inconsistencies in the implementation of a widely used two-equation turbulence model were identified and corrections are recommended. It was also observed that unsteady simulations for centrifugal compressors require significantly longer integration times than what is current practice in industry.

Gorazd Medic↗
DOE OSTI2026

Integrated Simulation of Weld Residual Stress Evolution and Crack Propagation Using XFEM

Nuclear power plant components operate in environments that promote multiple degradation mecha- nisms, several of which involve crack initiation and growth. An ongoing effort in the U.S. Department of Energy’s Nuclear Energy Advanced Modeling and Simulation (NEAMS) program is developing a general capability within the Multiphysics Object Oriented Simulation Environment (MOOSE) framework for simulating three-dimensional crack growth under a range of driving conditions, including fatigue, stress corrosion cracking (SCC), brittle fracture, and stress-relaxation cracking. This report demonstrates an end-to-end workflow that uses this capability to model weld-residual-stress-driven SCC in the J-groove weld of a pressurized-water reactor control rod drive mechanism penetration in the vessel head. The workflow consists of a thermomechanical welding simulation with temperature-dependent plasticity, followed by cooldown to ambient conditions, and a restart of the simulation using the MOOSE extended finite element method (XFEM) module to propagate a three-dimensional crack through the residual stress field. New welding capabilities were developed to properly initialize newly activated elements in the weld region, and robustness improvements were made to the mesh-based algorithm for defining cutting planes in the 3D XFEM algorithm, allowing it to handle complex crack fronts and stress fields. Together these advances allowed the simulated SCC crack to grow from an initial elliptical flaw in the weld, across the weld, through the tube wall, and almost to the triple point (where the weld, tube, and reactor pressure vessel head intersect) over roughly 36 years of simulated service. These results demonstrate a workflow that can be extended to fully three-dimensional welding simulations and more complex crack interaction problems.

42 - ENGINEERING↗