SEARCH · Engineering Papers

Results for “Monte Carlo code”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

74 records · Page 3

All sources NASA NTRS DOE OSTI

NASA NTRSDate not supplied

Softening the Gap between Wöhler and Paris – New Approaches for Fatigue Analysis –

Fatigue analysis tools can vary across industries. For example, automotive engineers often use the Wöhler (S-N) approach to design for safe-life, while aerospace engineers prioritize damage tolerance and inspection intervals, relying instead on crack growth models such as Paris’ law. Although both approaches may deal with the control of cracks in similar materials, their analysis tools and material characterizations are fundamentally distinct. This divide mirrors the classic split between stress-based strength analysis and linear elastic fracture mechanics. However, modern nonlinear models that incorporate material softening, such as cohesive laws, blur this boundary and capture fracture behaviors across scales. This presentation describes the CF23 fatigue model, which uses cohesive softening to link S-N crack initiation with crack propagation rates. CF23 spans the full fatigue spectrum, from initial propagation transients to steady-state growth and threshold conditions, offering a unified framework that bridges Wöhler and Paris-based methodologies. Example applications include fatigue crack propagation transients in adhesive interfaces and skin/stiffener separation.

cohesive elements↗

NASA NTRS1983

Modeling a Li/SOCl2 battery for design purposes

A generalized code applicable to many different electrochemical systems and geometric designs is discussed. The code is to be set up so that physical property data such as thermal conductivity, viscosity, density, and configuration (e.g. physical dimensions) are the input data. Thus, by changing these parameters many different battery configurations can be handled. The outputs, as a function of time and space, are voltage, current, temperature, pressure, velocity, and species concentration.

Ernst, D. W.↗

NASA NTRS1995

Transitional Flow in Thin Tubes for Space Station Freedom Radiator

A two dimensional finite volume method is used to predict the film coefficients in the transitional flow region (laminar or turbulent) for the radiator panel tubes. The code used to perform this analysis is CAST (Computer Aided Simulation of Turbulent Flows). The information gathered from this code is then used to augment a Sinda85 model that predicts overall performance of the radiator. A final comparison is drawn between the results generated with a Sinda85 model using the Sinda85 provided transition region heat transfer correlations and the Sinda85 model using the CAST generated data.

Patrick Loney↗

NASA NTRSDate not supplied

Leveraging High-Level Synthesis to Migrate Motor Control Algorithms From Microcontroller to FPGA

As motor control algorithms become increasingly complex, traditional microcontroller-based implementations are reaching computational limits that prevent the controller from operating at the required speed. This paper presents a novel workflow leveraging High-Level Synthesis (HLS) to migrate motor control algorithms from a microcontroller implementation to a Field-Programmable Gate Array (FPGA) implementation. The proposed approach utilizes the free Vitis HLS software to automatically convert Embedded Coder-generated C code from a Simulink model into Hardware Description Language (HDL) code suitable for FPGA deployment.

FPGA↗

NASA NTRSDate not supplied

Leveraging High-Level Synthesis to Migrate Motor Control Algorithms From Microcontroller to FPGA

FPGA↗

NASA NTRS1996

Minimum Heating Re-Entry Trajectories for Advanced Hypersonic Launch Vehicles

Optimal re-entry trajectories are generated for reusable launch vehicles which minimize: (1) the heat absorbed at the vehicle surface, (2) the lower surface temperature, and (3) the heat absorbed by the internal structure. The approach uses the energy state approximation technique and a finite control volume heat transfer code coupled to a flight path integration code. These trajectories are compared to the optimal re-entry trajectory minimizing the integrated convective heat rate to determine which trajectory produces the minimum internal structural temperatures for a given thermal protection system. Three different thermal protection systems are considered: tile, blanket, and metallic.

Robert Windhorst↗

NASA NTRSDate not supplied

On-Orbit Calibration Assessment of NOAA-21 VIIRS Thermal Emissive Bands and Implications for JPSS-4 VIIRS

The Visible Infrared Imaging Radiometer Suite (VIIRS) instrument onboard the NOAA-21 (N21) satellite has been successfully operating for over three years collecting valuable scientific measurements. A large suite of weather forecasting models and scientific research applications are supported with the VIIRS measurements from N21 combined with those acquired from the VIIRS instrument onboard the S-NPP and NOAA-20 spacecraft. Among the 22 VIIRS spectral bands, there are seven thermal emissive bands (TEB) covering the 3.7 to 12.2 µm spectral range at two different spatial resolutions. The VIIRS TEB detectors are calibrated onboard using a blackbody at a controlled temperature and a space view for background signal measurement. In this paper, we discuss the on-orbit performance of the N21 VIIRS TEB using various parameters such as the detector gain, noise, and offset. This on-orbit assessment provides critical insights into instrument behavior and calibration performance that will inform pre-launch testing and post-launch validation strategies for the upcoming JPSS-4 VIIRS mission scheduled for launch in 2027. The methodologies and lessons learned from characterizing the N21 TEB performance will enable more efficient and comprehensive radiometric assessment of the next VIIRS instrument. Specific focus is placed on identifying performance trends and anomaly signatures observed during the N21 commissioning phase that can enhance the JPSS4 instrument checkout procedures and accelerate the transition to routine on-orbit operations. The established performance baselines and uncertainty estimates will serve as an acceptance criterion for the JPSS-4 commissioning activities, ensuring mission readiness and data continuity for operational users.

Amit Angal↗

NASA NTRS2026

Atacama Cosmology Telescope A Measurement of Galaxy Cluster Temperatures Through Relativistic Corrections to the Thermal Sunyaev-Zeldovich Effect

The high electron temperature in galaxy clusters (> 1 keV or > 107 K) leads to corrections at the level of a few percent in their thermal Sunyaev-Zeldovich effect signatures. Both the size and frequency dependence of these corrections, which are known as relativistic temperature corrections, depend upon the temperature of the objects. In this work we exploit this effect to measure the average temperature of a stack of Compton-y selected clusters. Specifically, we apply the “spectroscopic method” and search for the temperature that best fits the clusters’ signal measured at frequencies from 30 to 545 GHz by the Atacama Cosmology Telescope and Planck satellite. We measure the average temperature of clusters detected in the Atacama Cosmology Telescope maps to be 8.5 2.4 keV, with an additional systematic error of comparable amplitude dominated by passband uncertainty. Upcoming surveys, such as the Simons Observatory and CMB-S4, have the potential to dramatically improve upon these measurements and thereby enable precision studies of cluster temperatures with millimeter observations. The key challenge for future observations will be mitigating instrumental systematic effects, which already limit this analysis.

William R Coulton↗

NASA NTRSDate not supplied

Evaluation of Low-Energy Hydrogen Separation Method Using Metal-Organic Frameworks (MOFs) for Closed-Loop ECLSS Air Revitalization (CLEAR)

The State-of-the-Art (SOA) air revitalization architecture onboard the International Space Station (ISS) recovered approximately 50% of the oxygen (O 2 ) from metabolic carbon dioxide (CO 2 ) via the Sabatier process from 2011 to 2017. O 2 recovery is currently constrained by the limited availability of reactant hydrogen (H 2 ) preventing complete conversion of CO 2 to H 2 O. Increasing O 2 recovery within Closed-Loop ECLSS is essential to reducing resupply mass for long-duration manned missions; specifically focusing on water (H 2 O) which supplies H 2 for Sabatier via water electrolysis. Past ground test endeavors at Marshall Space Flight Center (MSFC) have attempted to recover H 2 from Sabatier-produced CH 4 using technologies such as carbon vapor deposition (CVD) and plasma pyrolysis. The byproducts of these technologies can act as a catalyst poison or reactor deadload to the Sabatier reactor. Hydrogen separation techniques must be utilized to maintain the Sabatier catalyst during gas recycling and must be scalable, non-energy intensive, and safe to operate in a habitation setting. Research indicated that metal-organic frameworks (MOFs) could meet these criteria and were tested for their capability to capture the various carbon-based gaseous products of CVD and plasma pyrolysis such as acetylene (C 2 H 2 ), ethylene (C 2 H 4 ), ethane (C 2 H 6 ), and carbon monoxide (CO) which would purify the hydrogen gas stream passing through the MOF. A sub-scale adsorption column was developed by Marshall Space Flight Center to test three MOF candidates against a synthetic gas mixture comprised of process-relevant carbonous gases and hydrogen to evaluate the separation capability of the MOFs. The results of the hydrogen separation capability, isothermal desorption capability, and demonstrated cyclic reuse of the MOF are presented in this paper.

Kagen Crawford↗

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↗

NASA NTRS2026

Global Model Estimates of Atmospheric Al, Ca, Fe, Si, and Ti from Dust and Non-Dust Aerosols Informed by EMIT Surface Mineralogy and Evaluated Against Observations

Atmospheric deposition of micro-nutrients like Fe has been shown to be important for ocean biogeochemistry. The largest source of atmospheric Fe and other elements (e.g., Ca, Al, Si, and Ti) is desert dust, although there are significant non-dust sources in some regions. However, past estimates of these elements have been substantially uncertain due to limited information about the composition of the desert source regions. Here we use elemental distributions estimated from new Earth Surface Mineral Dust Source Investigation (EMIT) observations, which provide mineralogical composition at the surface of the Earth based on imaging spectroscopy measurements from the International Space Station. We add in other sources of these elements (anthropogenic and natural) and compare to a compilation of available surface concentration data from stations over land and from shipborne observations. Our results suggest that the modeled distribution is similar to available observations, but discrepancies still exist in both natural desert dust regions as well as regions dominated by anthropogenic sources. Global budgets for the elements Ca, Al, Fe, Si, and Ti suggest that desert dust remains the dominant source for these elements but anthropogenic or volcanic sources are also important for these elements. Changes in elemental distributions since preindustrial times were also estimated.

aerosols↗

NASA NTRS2026

The Contaminant Footprint of Landed Spacecraft: Toward an Inventory and Modelling Framework

All spacecraft generate and carry contaminants, i.e., unwanted and potentially harmful material. When a spacecraft lands and operates in near-vacuum, as onto Earth’s Moon, it introduces contaminants into its environment that may compromise mission science objectives and engineering performance. Contamination of solar system bodies may irrevocably degrade targets of unique value to planetary scientists, for instance, as lunar landed spacecraft introduce propellant effluents into the otherwise pristine ice of the Moon’s permanently shadowed regions. NASA’s planetary protection discipline seeks to ensure that solar system bodies are not contaminated, for scientific purposes, by terrestrial material (i.e., forward contamination). This interest aligns with planetary science interest in mitigating the transport of terrestrial contaminants onto solar system bodies and in controlling types of contamination that could compromise the scientific value of samples or measurements. NASA, and other entities that practice planetary science, have compelling and multidisciplinary interests in the preservation of special regions and sampling sites of high scientific value – including lunar permanently shadowed regions (PSRs) – from inadvertent contamination by any spacecraft, and in understanding the contamination of such regions by all spacecraft. Organic molecular contamination here represents a primary threat. Organic molecules will be introduced to solar system bodies by nominal landed spacecraft and crew processes – including by the action of descent and ascent engines; natural materials outgassing; and crew environmental and life support system sources. Molecular contaminants can also travel in the free-molecular sense at global scale across near-vacuum bodies, including into regions where they may be permanently trapped. This presentation will address a high-level study to identify sources of contaminants – in particular, organic material – generated by landed spacecraft along with the transport vectors by which these contaminants can reach sites of scientific interest on bodies like the Moon. A vision for an integrated modeling framework for the organic contamination footprint of spacecraft missions, individually and collectively, will also be described and presented along with initial conclusions related to organic molecular transport.

Gas Dynamics↗

NASA NTRS2024

Expansion of Check-Cases for 6DOF Simulation

This is the Appendix containing a description of the solution for Case 1 in the assessment, “Expansion of Check-Cases for 6DOF Simulation”. For cases of spherical gravity, it is possible to provide a two-body solution without recourse to numerical integration and thus it is accurate to machine precision. Python code for a Keplerian Propagator (propagate.py) which produced a reference trajectory for Case 1 is provided in this appendix. There is also code for generating test cases which was used as an independent verification of the propagator. This is a high-level description of the algorithm employed. The documentation of each function includes implementation details, including equations for each task.

Modeling↗

NASA NTRS2020

Evaluation of Fatigue Damage Accumulation Functions for Delamination Initiation and Propagation

The present report follows on the cohesive fatigue damage model methodology proposed in NASA-TP-2018-219838. In that report, an empirical function describing the incremental damage due to cyclic loading was used to calculate fatigue damage within a cohesive formulation. The form of the function was developed such that, when integrated at a constant stress amplitude from no damage to failure, it produces a life versus load response that is consistent with an S-N diagram. Therefore, the parameters of the damage model could be obtained by fitting the model predictions to an S-N diagram. The finite element analyses performed demonstrate that the cohesive fatigue accumulation function provides a link between the S-N diagram that describes crack initiation, and the Paris law that characterizes the rate of crack propagation. However, when the model was proposed, it was not known whether the form of the damage accumulation function associated with a desired S-N diagram is unique and, if not, if the link between S-N and the Paris law is unique and independent of the fatigue function selected. In the effort described herein, several alternative forms of the damage function that reproduce the desired features of S-N diagrams were found and evaluated. The effects of each of these functions on the predicted parameters of the Paris law and the propagation threshold are discussed. The results indicate that the predicted exponent m of the Paris law is indeed independent of the damage accumulation function. However, different functions predict different values for the pre-factor C of the Paris law. Therefore, the proper damage accumulation function must be selected by comparison with experiments. One of the new damage accumulation functions proposed herein was found to be particularly useful for analysis because of the ease with which the model parameters can be determined with a minimal amount of experimental information. The effectiveness of the proposed methodology and damage function was demonstrated by conducting analyses of a double cantilever beam test, a mixed-mode bending test, and a three-point bending test of a skin/doubler specimen. The results indicate that the same set of model parameters can provide accurate predictions of the rate of fatigue crack propagation for a variety of material interfaces, mode mixities, load levels, and stress ratios.

DKIN/Stiffener Debonding↗

NASA NTRSDate not supplied

Turbo-Design: Open-Source Radial Equilibrium Turbomachinery Solver: Part I - Turbines

Advances in 3D Geometrical Designs and Cooling have played a significant role in improving the efficiency of turbomachinery. However, these advancements must be effectively translated back to the modeler. Machine learning can facilitate this transition. Specifically, machine learning–based loss models can bridge the gap between 3D and 1D designs, enabling modelers not only to predict velocity triangles but also to extract additional geometric features. Currently, the design tools used at NASA have not been updated to support such integration—until now. TurboDesign is an open-source, Python-based framework that replaces TD2 (LEW-11029-1) and AXOD2 (LEW-16323-1), both of which are radial equilibrium solvers for axial turbines. The goal of this update is to enable the integration of machine learning loss models into radial equilibrium equations. Additionally, TurboDesign is designed to support radial machines. This paper presents the governing equations, the assumptions underlying the code, the integration of legacy loss models, an example of machine learning model integration, and a validation comparison with CFD. All code, tutorials, and documentation are available at: https://www.github.com/nasa/turbo-design

Radial Equilibrium↗

NASA NTRS1991

Preparation of the Ice Certification of the Dornier 328 Regional Airliner By Numerical Simulation and By Ground Test

The Dornier 328, a new regional airliner, is to be qualified to FAR25/JAR25 requirements for operation into known icing conditions. All lifting surfaces are pneumatically deiced. Impingement limits were determined using a validated computational fluid dynamic (CFD) code. Wind tunnel tests with a model of the whole aircraft and the isolated empennage equipped with artificial ice shapes showed a degradation of handling characteristics due to ice accretion. Extensive two and three dimensional CFD calculations evaluated ice accretion at the deflection elevator horn. Icing tunnel tests were performed with a model of the horizontal tail with icing similitude fulfilled. The Dornier Do328 powerplant air induction system has successfully completed a program of icing tests in support of aircraft certification for operation into known icing conditions. Potential flight icing conditions were extensively analyzed using a CFD potential flow code and 3-D particle tracking routine to identify critical flight icing conditions. The tunnel test program verified the analytically predicted critical ice accretion surfaces of the air induction system, the adequacy of the ice protection provisions and demonstrated compliance with the applicable (JAR25) regulatory material.

D Welte↗

NASA NTRSDate not supplied

Lessons Learned From A Cryogenic Motor Configuration Study for Aircraft Propulsion

Cryogenic electric machines are a promising technology for multi megawatt electric aircraft drivetrains. Superconducting stators that make use of low AC loss superconducting wires are considered to be a key component for these cryogenic machines. NASA has undertaken 6-year technical challenge to address this need by designing a 5 MW superconducting motor and cryogenic drive and developing multi-MW prototypes. This poster shares the lessons learned from a completed design study that evaluated several configurations of a 5 MW cryogenic motor for this application. The lessons are drawn from a qualitative study of 16 configurations and a quantitative study of 8 down selected configurations. The quantitative study involved using a motor design code based on analytical calculations (electromagnetic, thermofluid, mechanical, and motor sizing) to determine each configuration’s Pareto front of total efficiency and total specific power. This poster builds upon preliminary results that were previously presented by (1) reporting final results based on a revised design code for each configuration, including a considerable increase in fidelity for a novel radial flux motor and (2) making a down selection of the configuration for future detailed design and demonstration and documenting the rationale for the selection.

Superconductors↗

NASA NTRS2026

Updated global and regional trends of stratospheric ozone profiles

We present updated evaluation of stratospheric ozone profile trends in the 60° S–60°N latitude range using long-term ground-based and satellite climate data records, as well as simulations by chemistry-climate models. The trends are evaluated using the LOTUS (Long-term Ozone Trends and Uncertainties in the Stratosphere) regression model. Analyses of satellite data confirm the statistically significant positive ozone trends in the period 2000–2024 in the upper stratosphere of ~1–3% per decade, with larger trends at mid-latitudes compared to the tropics. The trends are slightly positive or close to zero in the middle stratosphere, and mostly negative, -1 to -2% per decade, in the lower stratosphere, but they are not statistically significant. The morphology and magnitude of ozone trends are similar to previous analyses (2000–2020 trends). Ozone trends in 2000–2024 predicted by chemistry-climate model simulations are in good agreement with combined satellite trends. In the upper stratosphere, models predict a slightly stronger ozone recovery than observations. In the lower stratosphere, both models and satellite observations report negative trends in the tropics, while modelled ozone trends are slightly positive at mid-latitudes. Ozone profile trends over several stations estimated from ground-based records capture the same overall vertical pattern of ozone trends as merged gridded satellite datasets. Analyses of regional ozone profile trends in 2003–2024 using merged satellite datasets confirmed the previous observations of a longitudinal structure in ozone trends in the NH mid-latitude stratosphere, with positive trends over Scandinavia and negative trends over Siberia. However, the magnitude of this dipole-like structure is reduced compared to previous analyses.

trends↗