SEARCH · Engineering Papers

Results for “model 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.

At least 1,459 records · Page 81

All sources NASA NTRS DOE OSTI

NASA NTRSDate not supplied

Updates to PATO’s Thermo-Mechanical Model

The solid mechanics module within the Porous material Analysis Toolbox (PATO) material response code was updated to model the stress contribution from both internal pressure and pyrolysis shrinkage as function of material’s decomposition. Incorporating these physical phenomena, alongside the previously included temperature dependent mechanical properties, thermal expansion, and aerodynamic loads, allows for a more accurate thermo-mechanical behavior modeling of Thermal Protection System (TPS) materials, which is critical to assess spallation risks during atmospheric entry. This poster outlines the workflow for obtaining the necessary PATO inputs, derived from dilatometry experiments, which are used to model the stress field in TPS materials due to thermal expansion and pyrolysis shrinkage. Furthermore, to analyze the relevance of the stress contribution from the internal pressure, a study was conducted using TACOT with different permeability values.

TPS↗

DOE OSTI2025

Simulating Muon Capture with MARLEY

Core-collapse supernovae are intense sources of tens-of-MeV neutrinos. Each explosion releases about 1058 neutrinos in a burst lasting tens of seconds. As neutrinos are weakly interacting, they escape more easily compared to photons, allowing them to be observed before the light signal. Understanding low-energy neutrino interactions with argon will be critical for future supernova measurements with the Deep Underground Neutrino Experiment (DUNE). Interpreting future DUNE supernova neutrino data will require a high-quality model of O(10 MeV) neutrino interactions, but the dominant charged-current cross section has never been measured at these low energies for an argon target. We exploit the theoretical similarities between muon capture, which has been measured, and antineutrino scattering to test MARLEY (Model of Argon Reaction Low Energy Yields), the simulation code used by DUNE to model supernova neutrino reactions. I implemented C++ code to add a muon capture mode to MARLEY and compared the predictions of the new simulation to measured isotope branching ratios of the residual nuclei. The preliminary branching ratios predicted by MARLEY running in muon capture mode demonstrate promising baseline performance when compared to experimental data and are expected to improve with further model adjustments.

Wong, Baker Louie↗

NASA NTRS2007

Toward Scientific Numerical Modeling

Ultimately, scientific numerical models need quantified output uncertainties so that modeling can evolve to better match reality. Documenting model input uncertainties and verifying that numerical models are translated into code correctly, however, are necessary first steps toward that goal. Without known input parameter uncertainties, model sensitivities are all one can determine, and without code verification, output uncertainties are simply not reliable. To address these two shortcomings, two proposals are offered: (1) an unobtrusive mechanism to document input parameter uncertainties in situ and (2) an adaptation of the Scientific Method to numerical model development and deployment. Because these two steps require changes in the computational simulation community to bear fruit, they are presented in terms of the Beckhard-Harris-Gleicher change model.

Kleb, Bil↗

NASA NTRS2004

Statistical Analysis of CMC Constituent and Processing Data

Ceramic Matrix Composites (CMCs) are the next "big thing" in high-temperature structural materials. In the case of jet engines, it is widely believed that the metallic superalloys currently being utilized for hot structures (combustors, shrouds, turbine vanes and blades) are nearing their potential limits of improvement. In order to allow for increased turbine temperatures to increase engine efficiency, material scientists have begun looking toward advanced CMCs and SiC/SiC composites in particular. Ceramic composites provide greater strength-to-weight ratios at higher temperatures than metallic alloys, but at the same time require greater challenges in micro-structural optimization that in turn increases the cost of the material as well as increases the risk of variability in the material s thermo-structural behavior. to model various potential CMC engine materials and examines the current variability in these properties due to variability in component processing conditions and constituent materials; then, to see how processing and constituent variations effect key strength, stiffness, and thermal properties of the finished components. Basically, this means trying to model variations in the component s behavior by knowing what went into creating it. inter-phase and manufactured by chemical vapor infiltration (CVI) and melt infiltration (MI) were considered. Examinations of: (1) the percent constituents by volume, (2) the inter-phase thickness, (3) variations in the total porosity, and (4) variations in the chemical composition of the Sic fiber are carried out and modeled using various codes used here at NASA-Glenn (PCGina, NASALife, CEMCAN, etc...). The effects of these variations and the ranking of their respective influences on the various thermo-mechanical material properties are studied and compared to available test data. The properties of the materials as well as minor changes to geometry are then made to the computer model and the detrimental effects observed using statistical analysis software. The ultimate purpose of this study is to determine what variations in material processing can lead to the most critical changes in the materials property. The work I have taken part in this summer explores, in general, the key properties needed In this study SiC/SiC composites of varying architectures, utilizing a boron-nitride (BN)

Fornuff, Jonathan↗

DOE OSTI2022

Thermal-Hydraulics Modeling and Simulations of Hot Pool Using the SAS-CFD Coupled Code

The main goal of this activity is to test the dynamic coupling of the SAS4A/SASSYS-1 (SAS) and CFD models, using a recently patched version of the SAS code intended to address an undocumented limitation that hindered the Versatile Test Reactor (VTR) simulation efforts in FY21. As described in previous VTR calculation reports, the undocumented limitation in SAS v5.4 does not allow the user to activate the CFD coupling option during restart calculations. Since the analysts were unaware of this limitation, prior SAS-CFD simulation results for the protected station blackout (PSBO) transient were erroneous. Root-cause analysis was performed to determine the cause of this undocumented limitation in SAS v5.4, the SAS software was updated in a new patch, and the SAS-CFD simulations were repeated with this patched software. The results of the SAS-CFD simulations documented in this report show that the software patch does address the cited issue, and that the patched software indeed supports the activation of the CFD coupling model in restart calculations. The SAS development team will determine the schedule for implementing the patch in an official software release. This report documents updated SAS-CFD simulations of the PSBO transient response in the VTR. The hot pool is modeled with the CFD code STAR-CCM+, which is coupled at the flow boundaries to the SAS model of the primary heat transport system. SAS computes the mass flow rate and temperature at each core subassembly outlet, the thermal insulation cavity bypass, and the IHX inlet windows. CFD in turn computes the absolute pressure and temperature at each of these boundaries. The SAS code will ignore the temperature data at flow boundaries where flow is directed into the hot pool, i.e., at the core subassembly outlets unless flow reversal occurs. Similarly, CFD will ignore temperature data at boundaries where the flow is directed out of the hot pool, i.e., at the IHX inlets except under flow reversal. The focus of this work is to ensure that the SAS software patch addresses the undocumented limitation described in prior VTR calculation reports, rather than the accurate assessment of thermal stratification in the VTR during protected transients. This motivates the development of a new, simplified CFD model with a coarser mesh to accelerate the testing process. The updated model, and simplifying assumptions, are documented in this report. In future work, the thermal stratification assessment should be performed in more detail. The simplified CFD model can be improved by performing grid convergence studies sensitivity studies of turbulence parameters (e.g., Prandtl number, turbulence production and dissipation parameters) on temperature distributions and thermal stratification.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2023

Modeling of fission gas diffusion and release for Gd 2 O 3 doped UO 2

Uranium dioxide (UO 2 ) is the primary nuclear fuel in light water reactors, and its excess neutronic reactivity can be controlled by adding burnable absorbers, such as Gd 2 O 3 . This burnable absorber has a large neutron absorption cross-section, lowering the high reactivity of the reactor's initial fuel load. However, there needs to be more understanding of how added Gd 2 O 3 influences the properties of UO 2 under irradiation. To understand the behavior of defects and fission gas in the UO 2 /Gd 2 O 3 system under irradiation, we use cluster dynamics modeling supported by density functional theory calculations. First, we calculate the formation energies of Gd point and cluster defects, and evaluate the temperature-dependent defect concentrations using the defect formation energies and entropies. We show that Gd is soluble in UO 2 , introducing a negative charge in the system. Using this information, we adapted the cluster dynamics code Centipede to model the influence of Gd on U self-diffusion and Xe diffusion in UO 2 with 10 wt% Gd 2 O 3 . Also, we analyzed the Xe diffusion as a function of Gd 2 O 3 concentration, showing that the Xe diffusivity is decreased, which means that the athermal diffusivity due to electronic stopping persists at higher temperatures. In conclusion, the decrease in Xe diffusion means that more Xe stays in the matrix, decreasing the Xe release, and lowering its influence of fission gas release on the thermomechanical properties of UO 2 .

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2024

Modeling the Progenitor Stars of Observed Type IIP Supernovae

Abstract Type IIP supernovae (SNe IIP) are thought to originate from the explosion of massive stars >10M ⊙ . Their luminosity is primarily powered by the explosion energy and the radioactive decay energy of 56 Co, with the photosphere location regulated by hydrogen recombination. However, the physical connections between SNe IIP and their progenitor stars remain unclear. This paper presents a comprehensive study of SNe IIP and their progenitor stars by using the one-dimensional stellar evolution code,MESA. Our model grids consider the effects of stellar metallicity, mass, and rotation in the evolution of massive stars, as well as the explosion energy and 56 Ni production in modeling supernovae. To elucidate the observed SNe IIP and their origins, we compare their light curves (LCs) with our models. Furthermore, we investigate the impact of stellar parameters on LCs by considering stellar mass, metallicity, rotation, explosion energy, and 56 Ni production. We find that more massive stars exhibit longer plateaus due to increased photon diffusion time caused by massive ejecta. Higher metallicity leads to increased opacity and mass loss of progenitor stars. Rapid rotation affects internal stellar structures, enhancing convective mixing and mass loss, potentially affecting the plateau’s brightness and duration. Higher explosion energy results in brighter but shorter plateaus due to faster-moving ejecta. 56 Ni mass affects late-time luminosity and plateau duration, with larger masses leading to slower declines.

Astronomy & Astrophysics↗

DOE OSTI2025

Integrating Energy-Efficient Computing with Computational Research to Accelerate Energy Technology

NREL's computational sciences center hosts the largest high performance computing (HPC) capabilities dedicated to energy research while functioning as a living laboratory for energy-efficient computing. NREL's HPC capabilities support the research needs of the Department of Energy's Office of Energy Efficiency and Renewable Energy (EERE). In ten years of operation, HPC use in EERE-sponsored research has grown by a factor of 30, including work in electricity generation, energy efficiency, transportation, and energy system modeling. This paper analyzes this research portfolio, providing examples of individual use cases. The paper documents NREL's history of operating one of the world's most energy-efficient data centers while examining pathways to reduce economic and environmental impact beyond reduction of Power Usage Efficiency (PUE). This paper concludes by examining the unique opportunities created for accelerating improvements in data center efficiency created by combining an HPC system dedicated to energy research and a research program in energy-efficient computing.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

Simulation of the PBX 9501 Self-Ignited Heavily Confined Cookoff Test

The PBX 9501 self-ignited thermally damaged explosive (SITX) test was carried out by Holmes et al. A 6-inch diameter sphere made of the PBX 9501 was heated inside a steel confinement vessel until it underwent a thermal explosion. This test provided measured temperature data that was used to define a boundary condition to carry out post-experiment finite element (FE) thermal simulations and comparisons of the FE thermal response with the experimental time to thermal ignition and internal temperature measurements at the center of the PBX 9501 sphere. Finite element models of this experiment were constructed using the heat transfer code Aria. A cookoff model for PBX 9501 was used to simulate the thermal decomposition of PBX 9501 and to predict the time to thermal ignition. An inverse heat conduction scheme was designed and applied to estimate a 1D temperature boundary condition.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

DOE OSTI2023

HydroChrono: An Open-Source Hydrodynamics Package for Project Chrono

In this paper we present the development and verification of HydroChrono, a hydrodynamics package for the Project Chrono physics engine. This package includes the implementation of hydrodynamics equations, the added mass for multibody systems, the development of I/O functions as well as a Python API, and comparison against standard reference cases and other existing tools. HydroChrono provides a flexible, fully open-source solution for simulating wave energy converters (WECs), floating offshore wind turbines (FOWTs) platforms, and other hydrodynamic systems. Here we show, via comparisons with existing tools for benchmark verification cases, that HydroChrono accurately models hydrodynamic forces - making it a useful tool for the design and optimization of these systems. Additionally, the integration of HydroChrono with Project Chrono offers access to finite element modeling capabilities and high-fidelity modelling - with Chrono's existing coupling to CFD and SPH codes. This provides numerical modelers with a multifidelity simulation framework for designing and validating these systems. The development of HydroChrono provides a new, open-source solution for simulating hydrodynamic systems. Its compatibility with other simulation tools enables a more streamlined and efficient design process, advancing the field and providing new opportunities for innovation in this area.

BEM↗

NASA NTRS1977

Prediction of compliant wall drag reduction, part 1

Computer codes developed to test Bushnell's compliant wall drag reduction model are discussed. One code computes the evolution of mean velocity profiles during the period between bursts as forced by an imposed large-scale pressure pulse due to earlier bursts. Another code computes the local stability characteristics of these computed profiles. The programs use Chebyshev polynomials to resolve the normal boundary layer (y) direction and a staggered grid of mesh points to resolve the x direction. Typically, 257 grid points and 33 Chebyshev polynomials are used in the computations.

Orszag, S. A.↗

NASA NTRS1982

Radiative transfer model for heterogeneous 3-D scenes

A general mathematical framework for simulating processes in heterogeneous 3-D scenes is presented. Specifically, a model was designed and coded for application to radiative transfers in vegetative scenes. The model is unique in that it predicts (1) the directional spectral reflectance factors as a function of the sensor's azimuth and zenith angles and the sensor's position above the canopy, (2) the spectral absorption as a function of location within the scene, and (3) the directional spectral radiance as a function of the sensor's location within the scene. The model was shown to follow known physical principles of radiative transfer. Initial verification of the model as applied to a soybean row crop showed that the simulated directional reflectance data corresponded relatively well in gross trends to the measured data. However, the model can be greatly improved by incorporating more sophisticated and realistic anisotropic scattering algorithms

Kimes, D. S.↗

NASA NTRS1986

Photovoltaic conversion of laser power to electrical power

Photovoltaic laser to electric converters are attractive for use with a space-based laser power station. This paper presents the results of modeling studies for a silicon vertical junction converter used with a Nd laser. A computer code was developed for the model and this code was used to conduct a parametric study for a Si vertical junction converter consisting of one p-n junction irradiated with a Nd laser. These calculations predict an efficiency over 50 percent for an optimized converter.

Walker, G. H.↗

NASA NTRS1984

Stress analysis of Pathfinder-2 models

Stress analysis of the Pathfinder 2 fighter model was performed to determine a method for leading- and trailing-edge attachment that gives acceptable stress levels. Structural modeling of the wing was done using the finite element code SPAR. For the models studied, one ordinary lap joint was found to be satisfactory for the leading-edge flap, however, the alternating surface segmented lap joint method fo attachment was necessary for the trailing-edge flap to obtain acceptable stress levels.

Mehrotra, S. C.↗

NASA NTRS1988

Numerical investigations in three-dimensional internal flows

In order to develop an understanding of flow fields that exist in sidewall compression inlet models, an explicit, time-accurate full Navier-Stokes code was used to model them. Both laminar and turbulent boundary layers were numerically simulated. These simulations were examined as they were influenced by Mach number, changes in sidewall compression angle, hypersonic viscous interaction effects, and thick entering boundary layers. Cause and effect relationships were established. The numerical simulations were compared with the limited data available to aid in understanding both the simulations and the experiments.

Rose, William C.↗

NASA NTRS1992

Issues and approach to develop validated analysis tools for hypersonic flows: One perspective

Critical issues concerning the modeling of low-density hypervelocity flows where thermochemical nonequilibrium effects are pronounced are discussed. Emphasis is on the development of validated analysis tools. A description of the activity in the Ames Research Center's Aerothermodynamics Branch is also given. Inherent in the process is a strong synergism between ground test and real-gas computational fluid dynamics (CFD). Approaches to develop and/or enhance phenomenological models and incorporate them into computational flow-field simulation codes are discussed. These models have been partially validated with experimental data for flows where the gas temperature is raised (compressive flows). Expanding flows, where temperatures drop, however, exhibit somewhat different behavior. Experimental data for these expanding flow conditions are sparse; reliance must be made on intuition and guidance from computational chemistry to model transport processes under these conditions. Ground-based experimental studies used to provide necessary data for model development and validation are described. Included are the performance characteristics of high-enthalpy flow facilities, such as shock tubes and ballistic ranges.

Deiwert, George S.↗

NASA NTRS1995

Combustion system CFD modeling at GE Aircraft Engines

This viewgraph presentation discusses key features of current combustion system CFD modeling capabilities at GE Aircraft Engines provided by the CONCERT code; CONCERT development history; modeling applied for designing engine combustion systems; modeling applied to improve fundamental understanding; CONCERT3D results for current production combustors; CONCERT3D model of NASA/GE E3 combustor; HYBRID CONCERT CFD/Monte-Carlo modeling approach; and future modeling directions.

Burrus, D.↗

NASA NTRS2009

Design and Analysis Tools for Supersonic Inlets

Computational tools are being developed for the design and analysis of supersonic inlets. The objective is to update existing tools and provide design and low-order aerodynamic analysis capability for advanced inlet concepts. The Inlet Tools effort includes aspects of creating an electronic database of inlet design information, a document describing inlet design and analysis methods, a geometry model for describing the shape of inlets, and computer tools that implement the geometry model and methods. The geometry model has a set of basic inlet shapes that include pitot, two-dimensional, axisymmetric, and stream-traced inlet shapes. The inlet model divides the inlet flow field into parts that facilitate the design and analysis methods. The inlet geometry model constructs the inlet surfaces through the generation and transformation of planar entities based on key inlet design factors. Future efforts will focus on developing the inlet geometry model, the inlet design and analysis methods, a Fortran 95 code to implement the model and methods. Other computational platforms, such as Java, will also be explored.

Slater, John W.↗