Engineering PapersSearch

SEARCH · Engineering Papers

Results for “Multi-material simulations”

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.

Three-Dimensional, Multi-Phase Unsteady Flow of Water Undergoing Multiple Flow Regimes (Free Surface, Droplet, and Fluid Film)

Upgrades to the Mobile Launcher Rainbird nozzles were implemented for the Artemis II launch to increase water flow rates and improve the mobile launcher deck water coverage while minimizing the water throw distance to avoid impacting the Space Launch System (SLS) rocket nozzles. This analysis describes the use of a three-dimensional, multi-phase, free surface transient computational fluid dynamics (CFD) model to simulate the water flowing out through the Rainbird nozzle and ensure compliance with current Artemis requirements. Flow rates, pressures, velocities, nozzle throw distances and coverage areas were analyzed to ensure that system requirements are met. The completed results were used to make necessary modifications to the Ignition Overpressure Protection and Sound Suppression (IOP/SS) system Rainbird design by reshaping the nozzle used in the system.

CFD

Simulations of Yarn Micro-Mechanics of Woven Heat Shield Materials

Carbon and phenolic fibers are commonly used in ablative thermal protection materials, such as 3-dimensional Mid-Density Carbon Phenolic (3MDCP), a 3D-woven composite comprised of mixed-fiber yarn bundles. Predicting the micro-mechanical response and fracture of twisted yarns composed of brittle and ductile fibers requires a modeling approach that captures per-fiber yielding, fiber fracture, and inter-fiber friction and contact. This work presents an extended bonded particle model (BPM) for discrete element method (DEM) simulation of fiber and yarn mechanics, implemented in LAMMPS. The model builds on the incremental bond formulation of Guo et al. and introduces a piecewise elasto-plastic constitutive law for axial extension, enabling representation of fibers that yield before failure. 3MDCP yarns were constructed using measured fiber radius distributions and helical twist geometry. Tensile simulations of single-ply 3MDCP yarns show good agreement with vender stress–strain results. Fiber breakage models also show details on yarn breakage propagration, centered radially in the yarn. Yarn breakage of multi-ply 3MDCP also matched experimental observations in per-ply breakage; however, predicted yarn breakage strength were found higher than experimental observations.

Discrete Element Method

Simulations of Yarn Micro-Mechanics of Woven Heat Shield Materials (3MDCP)

Carbon and phenolic fibers are commonly used in ablative thermal protection materials, such as 3-dimensional Mid-Density Carbon Phenolic (3MDCP), a 3D-woven composite comprised of mixed-fiber yarn bundles. Predicting the micro-mechanical response and fracture of twisted yarns composed of brittle and ductile fibers requires a modeling approach that captures per-fiber yielding, fiber fracture, and inter-fiber friction and contact. This work presents an extended bonded particle model (BPM) for discrete element method (DEM) simulation of fiber and yarn mechanics, implemented in LAMMPS. The model builds on the incremental bond formulation of Guo et al. and introduces a piecewise elasto-plastic constitutive law for axial extension, enabling representation of fibers that yield before failure. 3MDCP yarns were constructed using measured fiber radius distributions and helical twist geometry. Tensile simulations of single-ply 3MDCP yarns show good agreement with vender stress–strain results. Fiber breakage models also show details on yarn breakage propagration, centered radially in the yarn. Yarn breakage of multi-ply 3MDCP also matched experimental observations in per-ply breakage; however, predicted yarn breakage strength were found higher than experimental observations.

Woven

Pragmatic Stress Prediction on Additively Manufactured Coupons

Prediction of residual stresses from process parameters for additively manufactured large metal parts is computationally expensive. NASA is currently developing meter-scale parts with direct energy deposition. Practically, the predictive computational methods need to efficiently scale-up to meter-scale parts. Coupled thermal-mechanical multi-physics simulations have been developed with the pragmatic method using ABAQUS, COMSOL Multiphysics, ALE3D software. The residual stresses are a result of the manufacturing process which creates thermal cycling of the build layers. The pragmatic method uses lumped thermal layers for stress predictions to reduce computational costs. The stress predictions as well as deformations of the different codes are compared with each other and with ANSYS Additive using identical material models, boundary and initial conditions. The codes were used to simulate three different geometries: a thin wall, hollow cylinder and twin-cantilever part. The coupon parts were then manufactured with Inconel-625. The residual stresses in these parts were measured using X-ray diffraction as well as neutron beam diffraction at NIST. The stress measurements for the two technologies are compared. The pragmatic stress prediction method enabled predictions of the multi-centimeter scale parts using desktop computer workstations in only a few hours for each coupon. The results of the simulated stress predictions compared favorably with the measured stresses even though thermally lumped layers were employed. Finally, a two-meter scale nozzle was simulated using ANSYS Additive. The simulations were used to examine the build orientation trade-space with respect to resulting geometric deformation. The predicted deformations were compared to measurements of an actual subscale part manufactured with direct energy deposition.

pragmatic method

Hybrid Modeling Study on Grain Evolution in the Metal Welding Process and Its Potential Lunar Application

Metal is most commonly used structural material in a wide range of spacecraft, and welding is the principal method for joining metal components into functional systems. However, conducting welding experiments under extreme environments—such as microgravity or vacuum conditions in space—is prohibitively expensive and experimentally challenging. To overcome these limitations, multi-physics computational welding models provide a cost-effective and versatile alternative. In this work, the authors have developed a coupled thermal (fluid) microstructure simulation framework to model metal welding under varying gravity conditions. The framework integrates a mixed-mode heat transfer formulation (conduction, convection, and radiation) with molten pool fluid dynamics, enabling accurate prediction of temperature fields and weld-pool geometry. A grain growth model is further incorporated to capture the spatial and temporal evolution of microstructure, including grain size distribution and morphological transitions during solidification. This approach provides detailed insight into molten pool evolution and grain-level microstructure development throughout the welding process. By explicitly parameterizing environmental conditions, the model supports extrapolation to off-Earth manufacturing scenarios such as welding on the lunar surface. Tantalum—chosen in this study due to its high melting point, oxidation resistance, and mechanical stability at elevated temperatures—serves as the material system for model demonstration. Beyond Tantalum, the integrated multi-physics framework offers broad applicability for predictive welding simulations of various structural and refractory metals or alloys used in extreme terrestrial or extraterrestrial environments.

kinetic Monte Carlo (SPPARKS)

Elastic-plastic models for multi-site damage

This paper presents recent developments in advanced analysis methods for the computation of stress site damage. The method of solution is based on the p-version of the finite element method. Its implementation was designed to permit extraction of linear stress intensity factors using a superconvergent extraction method (known as the contour integral method) and evaluation of the J-integral following an elastic-plastic analysis. Coarse meshes are adequate for obtaining accurate results supported by p-convergence data. The elastic-plastic analysis is based on the deformation theory of plasticity and the von Mises yield criterion. The model problem consists of an aluminum plate with six equally spaced holes and a crack emanating from each hole. The cracks are of different sizes. The panel is subjected to a remote tensile load. Experimental results are available for the panel. The plasticity analysis provided the same limit load as the experimentally determined load. The results of elastic-plastic analysis were compared with the results of linear elastic analysis in an effort to evaluate how plastic zone sizes influence the crack growth rates. The onset of net-section yielding was determined also. The results show that crack growth rate is accelerated by the presence of adjacent damage, and the critical crack size is shorter when the effects of plasticity are taken into consideration. This work also addresses the effects of alternative stress-strain laws: The elastic-ideally-plastic material model is compared against the Ramberg-Osgood model.

Actis, Ricardo L.

Multi-Agent Swarm State of the Art Report

The Next-Generation Multi-Agent Swarm (NGS) Study conducted by NASA’s Ames Research Center for NASA’s Space Technology Mission Directorate (STMD) will develop a comprehensive understanding of emerging multi-agent swarm capabilities. The study aims to identify existing swarm capabilities and asses their potential for persistent lunar space situational awareness, surface monitoring, and distributed autonomy demonstrations. A key objective is to inform the design of a next-generation multi-agent swarm that can perform autonomous distributed remote sensing, position, navigation, and timing (PNT) services, automated deployment that leverages autonomy, edge computing, and interoperable networking to enable cooperative operations without the need for immediate human operation. This study will address specific shortfalls identified by STMD, including intelligent multi-agent constellations, autonomy, edge computation, position, navigation, and timing for small spacecraft, small spacecraft propulsion, and space situational awareness (1625, 1438, 1433, 1557, 1431, 1430, 1589). The NASA Ames Mission Design Center (MDC) will provide subject matter expertise to support systems engineering trades, while experts in autonomy and spacecraft swarms in NASA’s Intelligent Systems Division will lead the study and focus on identifying emerging next-generation swarm capabilities. The study objectives include: capturing the current state-of-the-art for multi-agent swarm capabilities, evaluating technologies and creating technology roadmaps, and developing at least one new technology demonstration mission concept. This initial NGS study report surveys the current state of the art in technology areas relevant for the next-generation multi-agent swarm design. Our primary focus is on surveying relevant deployed space systems1, supplemented with selective analysis of relevant proposed missions and technology developments that have yet to fly.

agent

Control of Thermo-Acoustics Instabilities: The Multi-Scale Extended Kalman Approach

"Multi-Scale Extended Kalman" (MSEK) is a novel model-based control approach recently found to be effective for suppressing combustion instabilities in gas turbines. A control law formulated in this approach for fuel modulation demonstrated steady suppression of a high-frequency combustion instability (less than 500Hz) in a liquid-fuel combustion test rig under engine-realistic conditions. To make-up for severe transport-delays on control effect, the MSEK controller combines a wavelet -like Multi-Scale analysis and an Extended Kalman Observer to predict the thermo-acoustic states of combustion pressure perturbations. The commanded fuel modulation is composed of a damper action based on the predicted states, and a tones suppression action based on the Multi-Scale estimation of thermal excitations and other transient disturbances. The controller performs automatic adjustments of the gain and phase of these actions to minimize the Time-Scale Averaged Variances of the pressures inside the combustion zone and upstream of the injector. The successful demonstration of Active Combustion Control with this MSEK controller completed an important NASA milestone for the current research in advanced combustion technologies.

Dzu K Le

Artificial Intelligence for Multi-Mission Planetary Operations

A brief introduction is given to an automated system called the Spacecraft Health Automated Reasoning Prototype (SHARP). SHARP is designed to demonstrate automated health and status analysis for multi-mission spacecraft and ground data systems operations. The SHARP system combines conventional computer science methodologies with artificial intelligence techniques to produce an effective method for detecting and analyzing potential spacecraft and ground systems problems. The system performs real-time analysis of spacecraft and other related telemetry, and is also capable of examining data in historical context. Telecommunications link analysis of the Voyager II spacecraft is the initial focus for evaluation of the prototype in a real-time operations setting during the Voyager spacecraft encounter with Neptune in August, 1989. The preliminary results of the SHARP project and plans for future application of the technology are discussed.

David J Atkinson

Glenn Research Center Propulsion Systems Laboratory 2026 Customer Guide

This guide describes the Propulsion Systems Laboratory (PSL) at the NASA Glenn Research Center. It was written to help customers understand the various components involved in conducting a test program within the PSL. The PSL complex supports two large-engine test cells that simulate altitude flight conditions for a wide range of research and experimental tests. These test cells operate at altitudes up to 90,000 ft and speeds from subsonic to above supersonic. Test points such as pressure, temperature, and Mach number can be set at the engine or test article inlet based on customer requirements. The facility’s support systems include the heated and cooled combustion air systems; altitude exhaust system; hydraulic system; nitrogen, oxygen, and hydrogen systems; thrust measurement system, which includes the facility’s single- and multi-axis thrust stands; inlet system; and electrical systems. In addition to providing a detailed description of PSL systems and capabilities, this guide discusses the facility’s history and past tests and addresses facility safety procedures, pretest requirements, and test operation standards.

Engine Icing Test Facility

A Multi-Tier Autonomous Aerial Architecture for Wildfire Detection, Characterization, and Communication in Infrastructure-Denied Environments

Wildfire response depends on how fast an ignition can be confirmed and located, especially in remote regions where ground-based communication and monitoring may be limited. Geostationary sensors provide frequent observations but at kilometer-scale resolution, which is too coarse to resolve small fires in remote terrain. Ground camera networks require sightlines and infrastructure that back-country areas lack. To address these limitations, this work proposes a Multi-Tier Autonomous Wildfire Intelligence System that combines wide-area monitoring with targeted, high-resolution sensing. A solar-powered high-altitude long endurance (HALE) platform operating at approximately 60,000 ft provides persistent wide-area thermal and optical surveillance, running onboard edge inference to screen candidate ignitions and reduce false positives and downlink bandwidth. When a candidate ignition is detected, low-altitude uncrewed aircraft systems (UAS) can be deployed to conduct localized observations, including high-resolution imaging and atmospheric measurements such as wind and plume observation. By combining persistent detection with local sensing, the proposed architecture is designed to provide first responders with timely, high-resolution information about fire location and behavior to aid in emergency decision making.

Wildfire management, UAS, drones

An Overview of CMIP5 and CMIP6 Simulated Cloud Ice, Radiation Fields, Surface Wind Stress, Sea Surface Temperatures and Precipitation over Tropical and Subtropical Oceans

The potential links between ice water path (IWP), radiation, circulation, sea surface temperature (SST) and precipitation over the Pacific and Atlantic Oceans resulting from the falling ice radiative effects (FIREs) are examined from present day model outputs of CMIP5 and CMIP6. The latter is divided into two subsets with (SON6) and without FIREs (NOS6) as more models with FIREs are included in CMIP6 than in CMIP5. Improvement in floating cloud ice (~20 g m-2) is noticeable over convective regions in CMIP6 relative to CMIP5. The inclusion of FIREs in SON6 subset may contribute to reduce biases of overestimated outgoing longwave radiation and downward surface shortwave and overestimated reflected shortwave at the top of the atmosphere (TOA) by magnitudes of 4?8 W m-2 over convective regions against CERES, compared to NOS6 subset. The reduced biases in radiative fluxes in convective regions stabilize the atmosphere and lead to circulation, SST, cloud and precipitation changes over the trade-wind regions, as seen from improved radiative fluxes (4?15 W m-2), surface wind stress biases, SST (0.2?0.8 K) and precipitation (1 mm day-1) biases. The significant improvement from NOS6 to SON6 leads to improved multi-model means for CMIP6 relative to CMIP5 for radiation fields over the trade wind regions but the degradation over convective zones is attributed to NOS6 subset. The results suggest that other sources of uncertainty and deficiencies in climate models may play significant roles for reducing discrepancies although FIREs, via radiation-circulation coupling, may be one of the factors that help to reduce regional biases.

Jui-Lin F Li

Space Nuclear Propulsion for Space Applications: Multi-Mission Nuclear Thermal Propulsion Vehicle Design

This NASA Technical Memorandum (TM) is intended to capture design work that the Advanced Concepts Office (ACO) at NASA’s Marshall Space Flight Center (MSFC) performed on behalf of the Space Nuclear Propulsion (SNP) Project during an SNP-driven design and mission analysis initiative in FY2025. Specifically, this TM focuses on several spacecraft design concepts intended to explore the feasibility of nuclear thermal propulsion (NTP) for use in a selection of notional, unmanned space missions. While references to the missions will be made in this TM, the primary focus will be on the spacecraft design concepts themselves, the ground rules and assumptions used for designing the spacecraft, a description of the various subsystems and their performance parameters, and the methodologies used to design them. The TM will also highlight areas of future work and additional considerations that may need to be taken into account when evaluating the feasibility of NTP space systems for other applications. The report will discuss the results from ACO’s analysis as they pertain to the different vehicle concepts that were evaluated. Each discipline section in this report will discuss the discipline-specific ground rules and assumptions, the methodology and design approaches employed, and present the results along with any relevant discussion, and highlight areas of future work where applicable.

Mitchell A Rodriguez

Advances in Design Capabilities for Planetary Missions from the NASA Entry Systems Modeling and Instrumentation Portfolio

The Entry Systems Modeling project (ESM) is supported by both the NASA Space Technology and the Science Mission Directorates and focuses on developing simulation tools and validated models for characterizing the performance of entry systems tailored to planetary destinations across the Solar System. ESM is organized into six technical capability areas that together address all relevant factors related to spacecraft entry, as well as some aspects of descent: Thermal Protection System (TPS) Materials; Aerothermodynamics; Entry & Descent Vehicle Dynamics; Guidance, Navigation, and Control; Vehicle Systems Analysis; and Advanced Tools and Numerical Methods. Development within the capability areas is undertaken explicitly with a focus on transition and infusion to science missions, human exploration missions, and commercial space activities. The present talk details developments that specifically impact science missions, including simulation tool capabilities that aid in mission design and model development to understand entry system performance at a given destination. Examples of the successful infusion and transition of such project outcomes to science missions also are provided. Several simulation tool development efforts within ESM have resulted in new design capabilities for missions. One such outcome is improved toolsets for mission trajectory and concept of operations design. Specifically, an initiative to couple a leading tool for entry, ascent/descent, and orbital trajectory optimization (Program to Optimize Simulated Trajectories II or POST2) to those used within the Agency for interplanetary trajectory optimization (Copernicus and Monte) has made substantial progress, with the outcomes to date promising to allow efficient trajectory optimization across mission phases. Additionally, toolchains for the evaluation of vehicle performance during entry and descent have been developed that allow assessment of multi-dimensional aeroheating on detailed vehicle geometries, characterization of deployment and inflation of parachutes, and assessment of vehicle dynamic stability during descent. These capabilities are achieved by coupling diverse sets of physics together – material response, computational fluid dynamics, radiation, and vehicle dynamics – to suitably describe complex entry and descent phenomena. Several model development and validation efforts for specific destinations and entry regimes also are underway within the ESM project. For instance, new experimental capabilities to validate radiation models at low densities/high altitudes recently have been established with project support, specifically the Low-Density Shock Tube (LDST) at the NASA Ames Research Center Electric Arc Shock Tube (EAST) facility. The LDST is being leveraged to develop improved models of shock layer kinetics and radiation in Titan atmospheres, while future studies will be conducted in the LDST and the existing high velocity shock tube to provide validation data for radiation models of Venus, Ice Giants, and Mars atmospheres. Models describing the aerothermal and thermo-structural performance of Thermal Protection System (TPS) materials has been another focus, with multiscale modeling activities on-going for the two leading TPS materials applicable to a range of entry conditions and science missions: the Phenolic-Impregnated Carbon Ablator (PICA) and woven materials like 3D Mid-Density Carbon Phenolic (3MDCP). A continual effort is made to infuse and transition outcomes from ESM simulation tool and model development activities into relevant science missions. Significant progress has been made on this front, with missions such as Dragonfly, DAVINCI, and Mars Missions benefitting from project outcomes. The groundwork also is being laid to provide insights into forward looking missions to Gas/Ice Giants as well as for potential sample returns.

Justin Haskins

Summary of Research Report

Ten papers, published in various publications, on buckling, and the effects of imperfections on various structures are presented. These papers are: (1) Buckling mode localization in elastic plates due to misplacement in the stiffner location; (2) On vibrational imperfection sensitivity on Augusti's model structure in the vicinity of a non-linear static state; (3) Imperfection sensitivity due to elastic moduli in the Roorda Koiter frame; (4) Buckling mode localization in a multi-span periodic structure with a disorder in a single span; (5) Prediction of natural frequency and buckling load variability due to uncertainty in material properties by convex modeling; (6) Derivation of multi-dimensional ellipsoidal convex model for experimental data; (7) Passive control of buckling deformation via Anderson localization phenomenon; (8)Effect of the thickness and initial im perfection on buckling on composite cylindrical shells: asymptotic analysis and numerical results by BOSOR4 and PANDA2; (9) Worst case estimation of homology design by convex analysis; (10) Buckling of structures with uncertain imperfections - Personal perspective.

Isaac Elishakoff