Search NASA⌕ Search

SEARCH · Search NASA

Results for “COMPUTER MODELS”

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.

At least 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Computational Modeling of Graphite Degradation due to Molten Salt Infiltration and Wear

Molten-salt reactors (MSRs) represent a promising next-generation reactor design, with graphite serving as a moderator and/or reflector in several designs. However, due to limited experimental data and operational experience, a technical understanding of the structural integrity of graphite in molten salt environments remains incomplete. This report presents a modeling-based evaluation of graphite degradation in MSR environments, focusing on the effects of salt infiltration in fuel salt-based designs and surface wear in pebble bed reactor designs. The objective of this study is to enhance understanding of the structural integrity challenges posed by these degradation mechanisms and to provide a framework for assessing graphite behavior in MSRs. The first part of the report investigates the phenomenon of molten salt infiltration into graphite. This infiltration occurs when molten salt permeates the interconnected pore structure of the graphite moderator, driven by factors such as pressure differentials and the physical properties of both the salt and graphite. The infiltration process is influenced by characteristics of the pore structure, viscosity of the molten salt, and the interfacial energies between the graphite, salt, and the atmosphere within the graphite pore. Utilizing a coupled multiphysics modeling approach with Grizzly software, the study evaluates the stress induced by internal heat sources due to infiltration, which can lead to structural concerns. This evaluation is crucial for understanding how infiltration affects the mechanical integrity of graphite components in MSRs. The study considers the Molten-Salt Reactor Experiment (MSRE) graphite stringer geometry due to the availability of relevant data. Through detailed finite element analysis, the study examines stress distributions at varying infiltration percentages, revealing that stress levels increase with higher amounts of infiltration. Rare-event simulations, using the parallel subset simulation (PSS) framework, further quantify the failure probabilities under input uncertainties, with a user-specified failure metric. The PSS framework also identifies critical input parameters that significantly affect the stress values, including infiltration amount, thermal conductivity, and power density. Additionally, considering realistic reactor scenarios, the analysis was performed to account for the combined effects of radiation and infiltration, and modeling strategies on how to analyze new reactor designs or new graphite grades are discussed. The second part of the report focuses on wear mechanisms in pebble bed-based MSRs. As graphite fuel pebbles interact with the graphite reflector block, wear can result in material loss and the formation of surface defects, which may act as stress concentrators. A similar multiphysics modeling framework is employed to assess the impact of wear on the structural integrity of graphite components. This study considers a generic fluoride-cooled high-temperature reactor (gFHR) design due to the availability of comprehensive data. Worst-case scenario dimensions of the reflector blocks were analyzed under thermal and radiation conditions. Subsequently, wear in the form of idealized pits and grooves is modeled on the inner surface of the graphite block, with the maximum stress from previous simulations. The simulations show that groove-type defects are more detrimental than pits, leading to higher stress concentrations. Considering worst-case simulation scenarios and experimental wear rates, it was determined that the formation of a surface defect critical enough to affect the stress may not be possible in a gFHR design. Overall, the findings of this research contribute to the development of robust modeling tools for predicting graphite behavior under various operational conditions in MSRs.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1996

Computational Modeling of Multi-Phase/Multi-Species Flows with Applications to Liquid Rocket Engines

Accurate prediction of all physical phenomena in a combustion chamber is essential for better understanding of the system performance. Atomization, evaporation, combustion, chemical kinetics, and turbulence are those processes of great importance that need to be well understood. Processes involving the liquid phase in a combustion chamber will be further complicated under supercritical conditions. More advanced and accurate numerical techniques are required to extend our understanding of the above phenomena. A computer program for multi-species/multi-phase flow was developed for NASA/MSFC in 1992. This code, called Liquid Thrust Chamber Performance (LTCP) program takes an Eulerian- Eulerian approach and is based on the Total Variation Diminishing (TVD) technique with Lax-Friedrichs upwind method. Under the NASA/ASEE SFFP the LTCP code was used to predict the performance characteristics of several engines that were of particular interest to NASA. This code was also successful in a combustion detonation study. Converting the program to the PC platform was accomplished which extends usability and makes it available to a wider range of users. The Eulerian formulation of the liquid phase provides a suitable model that can be extended to include combustion modeling under supercritical conditions. The results have been compared against the ones of other codes and available measured data. The algorithm proved to be robust and efficient for problems with stiff source terms.

Navaz, Homayun K.↗
NASA NTRS2006

Computer Modeling of Thermal Convection in Melts to Explain Glass Formation in Low Gravity and on Earth

Experiments conducted up to this time on glass forming melts in the low gravity environment of space show that glasses prepared in low-g are more chemically homogeneous and more resistant to crystallization than the comparable glasses prepared at 1-g on Earth. This result is somewhat surprising and opposite to the accepted concept on glass formation for a melt. A hypothesis based on "shear thinning" of a melt, a decrease in viscosity with increasing shear stress, is proposed as an explanation for the observed low-gravity results. This paper describes detailed simulation procedures to test the role of thermal convection in introducing shear stress in glass forming melts, using a lithium disilcate melt as a model. The simulation system in its idealized version consists of a cylinder that is heated at one end and cooled at the other with gravity acting in a transverse direction to the thermal gradient. The side wall of the cylinder is assumed to be insulating. The governing equations of motion and energy are solved using variable properties for viscosity (Arrehenius and non-Arrehenius behaviors) and density (constant and temperature dependent). Other parametric variables in the calculations include gravity level and gravity vector orientation. The shear stress in the system are then computed as a function of gravity from the calculated values of maximum melt velocity, and its effect on melt viscosity (shear thinning) is predicted. Also included and discussed are the modeling efforts related to other potential convective processes in glass forming melts and their possible effects on melt viscosity.

Ray, Chandra S.↗
DOE OSTI2018

Computational Modeling of Exciton-Bath Hamiltonians for Light Harvesting 2 and Light Harvesting 3 Complexes of Purple Photosynthetic Bacteria at Room Temperature

Light harvesting 2 (LH2) complex is the primary component of the photosynthetic unit of purple bacteria that is responsible for harvesting and relaying excitons. The electronic absorption line shape of LH2 contains two major bands at 800 and 850 nm wavelength regions. Under low light conditions, some species of purple bacteria replace LH2 with light harvesting 3 (LH3), a variant form with almost the same structure as the former but with distinctively different spectral features. The major difference between the absorption line shapes of LH2 and LH3 is the shift of the 850 nm band of the former to a new 820 nm region. The microscopic origin of this difference has been the subject of some theoretical/computational investigations. However, the genuine molecular level source of such a difference is not clearly understood yet. This work reports a comprehensive computational study of LH2 and LH3 complexes so as to clarify different molecular level features of LH2 and LH3 complexes and to construct simple exciton-bath models with a common form. All-atomistic molecular dynamics simulations of both LH2 and LH3 complexes provide detailed molecular level structural differences of bacteriochlorophylls (BChls) in the two complexes, in particular, in their patterns of hydrogen bonding (HB) and torsional angles of the acetyl group. Time-dependent density functional theory calculation of the excitation energies of BChls for structures sampled from the MD simulations suggests that the observed differences in the HB and torsional angles cannot fully account for the experimentally observed spectral shift of LH3. Potential sources that can explain the actual spectral shift of LH3 are discussed, and their magnitudes are assessed through fitting of experimental line shapes. Furthermore, these results demonstrate the feasibility of developing simple exciton-bath models for both LH2 and LH3, which can be employed for large-scale exciton quantum dynamics in their aggregates.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Computational modeling of coupled mechanical damage and electrochemistry in ternary oxide composite electrodes

Performance degradation of ternary layered oxide cathodes largely originates from their loss of structural integrity in cyclic usage. Mechanical damage, such as intergranular fracture of the active particles, is not only a mechanical cleavage process but also interferes with electrochemical kinetics such as infiltration of liquid electrolyte, surface corrosion of the constituent primary particles, and may eventually isolate the primary grains from the electron conducting network. Here, in this work, we develop a computational framework that integrates electrochemistry of a LiNi x Mn y Co 1−x−y O 2 (NMC) composite cathode with mechanical damage of the active particles. To fully examine the intricate chemomechanical behavior of the electrode, we evaluate the effects of the anisotropic material properties, the influence of mechanical potential on Li transport, and the concurrent intergranular fracture and electrolyte penetration along the grain boundaries upon multiple cycles. Electrolyte infiltration benefits capacity retention but aggravates further mechanical damage by corrosion. Structural failure mostly occurs in the first charging due to the anisotropic mechanical strain between the primary grains, while the resulting damage remains stable in the later few cycles. The results are consistent with experimental observations and the integration of electrochemistry and mechanical failure enables a step further understanding of the complex mechanism of battery degradation.

Battery degradation↗
NASA NTRS2010

A Mass Computation Model for Lightweight Brayton Cycle Regenerator Heat Exchangers

Based on a theoretical analysis of convective heat transfer across large internal surface areas, this paper discusses the design implications for generating lightweight gas-gas heat exchanger designs by packaging such areas into compact three-dimensional shapes. Allowances are made for hot and cold inlet and outlet headers for assembly of completed regenerator (or recuperator) heat exchanger units into closed cycle gas turbine flow ducting. Surface area and resulting volume and mass requirements are computed for a range of heat exchanger effectiveness values and internal heat transfer coefficients. Benefit cost curves show the effect of increasing heat exchanger effectiveness on Brayton cycle thermodynamic efficiency on the plus side, while also illustrating the cost in heat exchanger required surface area, volume, and mass requirements as effectiveness is increased. The equations derived for counterflow and crossflow configurations show that as effectiveness values approach unity, or 100 percent, the required surface area, and hence heat exchanger volume and mass tend toward infinity, since the implication is that heat is transferred at a zero temperature difference. To verify the dimensional accuracy of the regenerator mass computational procedure, calculation of a regenerator specific mass, that is, heat exchanger weight per unit working fluid mass flow, is performed in both English and SI units. Identical numerical values for the specific mass parameter, whether expressed in lb/(lb/sec) or kg/(kg/sec), show the dimensional consistency of overall results.

Juhasz, Albert J.↗
NASA NTRS2010

A Mass Computation Model for Lightweight Brayton Cycle Regenerator Heat Exchangers

Based on a theoretical analysis of convective heat transfer across large internal surface areas, this paper discusses the design implications for generating lightweight gas-gas heat exchanger designs by packaging such areas into compact three-dimensional shapes. Allowances are made for hot and cold inlet and outlet headers for assembly of completed regenerator (or recuperator) heat exchanger units into closed cycle gas turbine flow ducting. Surface area and resulting volume and mass requirements are computed for a range of heat exchanger effectiveness values and internal heat transfer coefficients. Benefit cost curves show the effect of increasing heat exchanger effectiveness on Brayton cycle thermodynamic efficiency on the plus side, while also illustrating the cost in heat exchanger required surface area, volume, and mass requirements as effectiveness is increased. The equations derived for counterflow and crossflow configurations show that as effectiveness values approach unity, or 100 percent, the required surface area, and hence heat exchanger volume and mass tend toward infinity, since the implication is that heat is transferred at a zero temperature difference. To verify the dimensional accuracy of the regenerator mass computational procedure, calculation of a regenerator specific mass, that is, heat exchanger weight per unit working fluid mass flow, is performed in both English and SI units. Identical numerical values for the specific mass parameter, whether expressed in lb/(lb/sec) or kg/ (kg/sec), show the dimensional consistency of overall results.

Juhasz, Albert J.↗
NASA NTRS2000

Computer Modeling of Real-Time Dynamic Lighting

Space Station tasks involve procedures that are very complex and highly dependent on the availability of visual information. In many situations, cameras are used as tools to help overcome the visual and physical restrictions associated with space flight. However, these cameras are effected by the dynamic lighting conditions of space. Training for these is conditions is necessary. The current project builds on the findings of an earlier NRA funded project, which revealed improved performance by humans when trained with computer graphics and lighting effects such as shadows and glare.

Maida, James C.↗
DOE OSTI2020

Computational Modeling To Adapt Neutralizing Antibody

Monoclonal antibodies (mAbs) is the leading therapy for viral infections because it provides immediate protection and can be administered at higher levels than in a natural immune response. Finding mAbs that neutralize a broad spectrum of viral targets has proven difficult because many species and strains exist and blanket targeting is a slow and laborious process to experimentally screen 108 variants. A new method is needed to rapidly redesign mAbs for homologous targets. This project speeds up redesign using structure-based computational design to reduce the mAbs search space to a manageable level and screen mutants at a much higher rate than in experiments. Computation will also provide critical knowledge about the fundamental interactions. The project will adapt S230, a human antibody that neutralizes SARS-CoV, to neutralize SARS-COV-2.

60 APPLIED LIFE SCIENCES↗
NASA NTRS1994

Computer modeling of the sensitivity of a laser water vapor sensor to variations in temperature and air speed

Currently, there is disagreement among existing methods of determining atmospheric water vapor concentration at dew-points below -40 C. A major source of error is wall effects which result from the necessity of bringing samples into the instruments. All of these instruments also have response times on the order of seconds. NASA Langley is developing a water vapor sensor which utilizes the absorption of the infrared radiation produced by a diode laser to estimate water vapor concentration. The laser beam is directed through an aircraft window to a retroreflector located on an engine. The reflected beam is detected by an infrared detector located near the laser. To maximize signal to noise, derivative signals are analyzed. By measuring the 2f/DC signal and correcting for ambient temperature, atmospheric pressure and air speed (which results in a Doppler shifting of the laser beam), the water vapor concentration can be retrieved. Since this is an in situ measurement there are no wall effects and measurements can be made at a rate of more than 20 per second. This allows small spatial variations of water vapor to be studied. In order to study the sensitivity of the instrument to variations in temperature and air speed, a computer program which generated the 2f, 3f, 4f, DC and 2f/DC signals of the instrument as a function of temperature, pressure and air speed was written. This model was used to determine the effect of errors in measurement of the temperature and air speed on the measured water vapor concentration. Future studies will quantify the effect of pressure measurement errors, which are expected to be very small. As a result of these studied, a retrieval algorithm has been formulated, and will be applied to data taken during the PEM-West atmospheric science field mission. Spectroscopic studies of the water vapor line used by the instrument will be used to refine this algorithm. To prepare for these studies, several lasers have been studied to determine their output frequency range and power.

Tucker, George F.↗
DOE OSTI2020

Computational modeling and neutron imaging to understand interface shape and solute segregation during the vertical gradient freeze growth of BaBrCl:Eu

In this work we apply continuum models to analyze phase change, heat transfer, fluid flow, solute transport, and segregation in order to understand prior neutron imaging observations of the vertical gradient freeze growth of Eu-doped BaBrCl. The models provide a rigorous framework in which to understand the mechanisms that are responsible for the complicated evolution of interface shape and dopant distribution in the growth experiment. We explain how a transition in the solid/liquid interface shape from concave to convex is driven by changes in radial heat transfer caused by furnace design. We also provide a mechanistic explanation of how dynamic growth conditions and changes of the flow structure in the melt result in complicated segregation patterns in this system. A growth pause caused by controller lock-up is shown to result in a band of solute depletion in accordance with classical theory. However, changing flow patterns during growth result in a non-monotonic axial distribution of solute that cannot be explained by simple application of classical segregation models. We assert that the approach presented here, namely the use of rigorous models in conjunction advanced diagnostics, such as neutron imaging, provides an exciting path forward for process optimization and control, accelerating the incremental advances that have, in the past, typically relied on empiricism, experience, and intuition.

36 MATERIALS SCIENCE↗
NASA NTRS2007

Computational Modeling Develops Ultra-Hard Steel

Glenn Research Center's Mechanical Components Branch developed a spiral bevel or face gear test rig for testing thermal behavior, surface fatigue, strain, vibration, and noise; a full-scale, 500-horsepower helicopter main-rotor transmission testing stand; a gear rig that allows fundamental studies of the dynamic behavior of gear systems and gear noise; and a high-speed helical gear test for analyzing thermal behavior for rotorcraft. The test rig provides accelerated fatigue life testing for standard spur gears at speeds of up to 10,000 rotations per minute. The test rig enables engineers to investigate the effects of materials, heat treat, shot peen, lubricants, and other factors on the gear's performance. QuesTek Innovations LLC, based in Evanston, Illinois, recently developed a carburized, martensitic gear steel with an ultra-hard case using its computational design methodology, but needed to verify surface fatigue, lifecycle performance, and overall reliability. The Battelle Memorial Institute introduced the company to researchers at Glenn's Mechanical Components Branch and facilitated a partnership allowing researchers at the NASA Center to conduct spur gear fatigue testing for the company. Testing revealed that QuesTek's gear steel outperforms the current state-of-the-art alloys used for aviation gears in contact fatigue by almost 300 percent. With the confidence and credibility provided by the NASA testing, QuesTek is commercializing two new steel alloys. Uses for this new class of steel are limitless in areas that demand exceptional strength for high throughput applications.

Source record↗
DOE OSTI2023

FY23 Progress on Computational Modeling of the Water-Based NSTF

This report summarizes the system-level modeling effort by Argonne National Laboratory (Argonne) of the Natural convection Shutdown heat removal Test Facility (NSTF) in FY23. As a continuation of the modeling effort from FY22, this year’s work focuses on improving the RELAP5-3D model developed previously for two-phase flow simulations. The RELAP5-3D model is updated to more accurately capture the heat loss experienced by the facility. The updated model is compared against experimental data for benchmarking purposes of the RELAP5-3D input model. By correctly accounting for heat loss, the updated RELAP5-3D model can now predict the two-phase baseline case more accurately. The onset and the duration of instability are captured well by the model. Furthermore, analyses are performed to better understand the instability mechanism experienced by the flow where the expansion of the boiling boundary in the chimney is studied in details and the fundamental frequencies of the oscillations are obtained. The updated RELAP5-3D model is further compared against four fault conditions, namely the reduction of riser header inlet flow area, depletion of system inventory, blocked riser channels, and static boiling scenario. For each fault condition, minor modifications and tuning are performed to improve the predictions of the model. The purpose of the analyses is to investigate the capability of RELAP5-3D in predicting complex two-phase flows in possible accident scenarios in actual Reactor Cavity Cooling System (RCCS). Overall, the model is able to capture the behaviors and trends of these fault conditions relatively well. Some discrepancies remain between the experimental data and the predictions, many of which are likely due to the differences in the predicted and experimental vapor generation rate. Future work will focus on the continued development of the current RELAP5-3D input model of the NSTF to both improve the accuracy of the model’s predictive capability and continue supporting the experimental program needs. The mutually beneficial relationship between analysis and experimental efforts has become integral to the parent NSTF program, and the greater objective to fully understand and accurately predict the heat removal performance of a full scale RCCS concept.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

Computational Modeling of CO2 Capture By Novel PIM-RU in a Fluidized Bed Riser

The interest in carbon capture, storage, and utilization (CCUS) has increased significantly in the past few decades as it can help mitigate the threat of global warming caused by substantial increase in CO2 emissions due to anthropogenic activities. Although several CO2 capture technologies have been developed, porous solid sorbents, which adsorb CO2 by physisorption, are considered promising candidates for post-combustion CO2 capture because of the easier recovery of adsorbed CO2 and high material stability. Some prominent types of porous solid sorbents are metal-organic frameworks (MOFs), Zeolite, mesoporous silica, and polymer-based sorbents (e.g., polymers with intrinsic microporosity or PIM). Oak Ridge National Laboratory (ORNL) recently developed a novel PIM-based sorbent, referred to here as PIM-RU. The main objective of this work is to investigate the CO2 capture performance of this sorbent using computational fluid dynamics (CFD). While process configuration and reactor design are open questions, a fluidized bed riser was selected as the contactor type for this study.

Aziz, Hossain↗
DOE OSTI2024

Procedure for the Computational Models created for the DOE/NRC Collaboration for Criticality Safety Support for Commercial-Scale HALEU Fuel Cycles Project (DNCSH)

This procedure outlines the application models developed as part of the US Department of Energy (DOE) and Nuclear Regulatory Commission (NRC) collaboration for criticality safety support for commercial scale HALEU fuel cycle and transportation (DNCSH) project, an effort authorized by the US Congress. The project is a joint effort among the DOE, NRC, numerous national laboratories, and private enterprises with project management from Oak Ridge National Laboratory.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS2006

Computer Modeling of Protocellular Functions: Peptide Insertion in Membranes

Lipid vesicles became the precursors to protocells by acquiring the capabilities needed to survive and reproduce. These include transport of ions, nutrients and waste products across cell walls and capture of energy and its conversion into a chemically usable form. In modem organisms these functions are carried out by membrane-bound proteins (about 30% of the genome codes for this kind of proteins). A number of properties of alpha-helical peptides suggest that their associations are excellent candidates for protobiological precursors of proteins. In particular, some simple a-helical peptides can aggregate spontaneously and form functional channels. This process can be described conceptually by a three-step thermodynamic cycle: 1 - folding of helices at the water-membrane interface, 2 - helix insertion into the lipid bilayer and 3 - specific interactions of these helices that result in functional tertiary structures. Although a crucial step, helix insertion has not been adequately studied because of the insolubility and aggregation of hydrophobic peptides. In this work, we use computer simulation methods (Molecular Dynamics) to characterize the energetics of helix insertion and we discuss its importance in an evolutionary context. Specifically, helices could self-assemble only if their interactions were sufficiently strong to compensate the unfavorable Free Energy of insertion of individual helices into membranes, providing a selection mechanism for protobiological evolution.

Rodriquez-Gomez, D.↗
DOE OSTI2019

High Fidelity Computational Model for Fluidized Bed Experiments (Final Report)

MFiX was integrated with preconditioners and linear solver packages in Trilinos via MFIX, Fortran, C, and C++ wrappers. The MFIX wrapper interpreted the structure of the matrix and vector from MFiX while the Fortran wrapper transfers this information to a low-level language, C. C-wrapper transforms the memory semantic between Fortran and C++ language. C++ wrapper passes the matrix and vector to a Trilinos package and gets the solution. This solution was transferred to MFiX via C++, C, Fortran, MFIX wrapper. The framework is used to integrate first, second and third generation linear solvers in Trilinos with MFiX.

97 MATHEMATICS AND COMPUTING↗