Engineering Papers⌕ Search

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.

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

At least 829 records · Page 46

A numerical study of chemically reacting flow in nozzles

The space station uses small rocket motors, called thrusters, for orientation control. Because of the lack of viable design tools for small rockets, the initial thruster design was basically a very small version of a large rocket motor. Thrust measurements of the initial design were lower than predicted. To improve predictions it was decided to develop a verison of the RPLUS2D reacting flow code for thruster calculations. RPLUS2D employs an implicit finite volume, lower-upper symmetric successive overrelaxation (LU-SSOR) scheme for solving the complete two-dimensional Navier-Stokes equations and species transport equations in a coupled and very efficient manner. The combustion processes are modeled by a 9-species, 18 step finite-rate chemistry model, and the turbulence is simulated by a Baldwin-Lomax algebraic model. The code is extended to handle multiple subsonic inlet conditions where the total mass flow is governed by conditions calculated at the thruster-throat. Results are shown for a thruster design where the overall mixture ratio is hydrogen rich. A calculation of a large area ratio divergent nozzle is also presented.

Vanoverbeke, Thomas J.↗

A numerical study of chemically reacting flow in nozzles

The Space Station uses small rocket motors, called thrusters, for orientation control. Because of the lack of viable design tools for small rockets, the initial thruster design was basically a very small version of a large rocket motor. Thrust measurements of the initial design were lower than predicted. To improve predictions it was decided to develop a version of the RPLUS2D reacting flow code for thruster calculations. RPLUS2D employs an implicit finite volume, lower-upper symmetric successive overrelaxation (LU-SSOR) scheme for solving the complete two-dimensional Navier-Stokes equations and species transport equations in a coupled and very efficient manner. The combustion processes are modeled by a 9-species, 18 step finite-rate chemistry model, and the turbulence is simulated by a Baldwin-Lomax algebraic model. The code is extended to handle multiple subsonic inlet conditions where the total mass flow is governed by conditions calculated at the thruster-throat. Results are shown for a thruster design where the overall mixture ratio is hydrogen rich. A calculation of a large area ratio divergent nozzle is also presented.

Vanoverbeke, Thomas J.↗

Nodal Modeling of Tank Pressurization and Draining using a Multi-Node-Ullage Approach

The purpose of the pressurization system in liquid rocket propulsion is to control the pressure in the gas space of the propellant tank (known as the ullage space) and the propellant mass flowrate to the engine. A mathematical model is required to predict the amount of pressurant necessary to ensure that pressure and temperature levels inside the tank remain within acceptable limits and that the propellant pressure leaving the tank satisfies the net positive suction pressure (NPSP) requirement of the pump feeding the engine. Nodal codes typically model tank pressurization and draining using a single node to represent the ullage and a single node to represent the propellant. As the tank drains, the ullage node grows and the propellant node shrinks. The heat transfer between ullage to wall and ullage to propellant is governed by natural convection. Designers of liquid propulsion systems often use empirical correlations to estimate the “Collapse Factor” which represents the ratio of pressurant required with heat transfer and the amount of pressurant required without heat transfer. A single node ullage model of tank pressurization was developed using GFSSP to compute the collapse factor reasonably well and later was used to model tank pressurization during test firing of the FASTRAC rocket engine. The predicted tank pressure compared well with the test data. In the early 1970’s, pressurization and drain tests with liquid methane were performed at NASA Lewis Research Center in a vacuum chamber. A 5 ft diameter spherical aluminum tank was tested to drain from 95% to 5% full using gaseous helium, hydrogen, nitrogen and methane as pressurant. Tests were conducted with different pressurant inlet temperatures and drain times. Measured data include pressurant requirement, amount of pressurant condensed, and ullage and wall temperatures at various heights in the ullage space at the end of draining. A single node GFSSP model was developed to simulate helium pressurization of the methane tank. Predicted helium consumption was 8-23% less than measured. The average error of the six test cases was 16%. A single ullage node with multiple solid node model was developed using Thermal Desktop. Predicted helium consumption compares with the test data within 2%. This paper describes the development of a GFSSP multi-node ullage model of the test configuration and compares the predicted pressurant consumption for both helium and autogenous pressurization using gaseous methane with experimental as well as TD predictions.

Nodal Model↗

Atmospheric Correction of Satellite Imagery Using Modtran 3.5 Code

When performing satellite remote sensing of the earth in the solar spectrum, atmospheric scattering and absorption effects provide the sensors corrupted information about the target's radiance characteristics. We are faced with the problem of reconstructing the signal that was reflected from the target, from the data sensed by the remote sensing instrument. This article presents a method for simulating radiance characteristic curves of satellite images using a MODTRAN 3.5 band model (BM) code to solve the radiative transfer equation (RTE), and proposes a method for the implementation of an adaptive system for automated atmospheric corrections. The simulation procedure is carried out as follows: (1) for each satellite digital image a radiance characteristic curve is obtained by performing a digital number (DN) to radiance conversion, (2) using MODTRAN 3.5 a simulation of the images characteristic curves is generated, (3) the output of the code is processed to generate radiance characteristic curves for the simulated cases. The simulation algorithm was used to simulate Landsat Thematic Mapper (TM) images for two types of locations: the ocean surface, and a forest surface. The simulation procedure was validated by computing the error between the empirical and simulated radiance curves. While results in the visible region of the spectrum where not very accurate, those for the infrared region of the spectrum were encouraging. This information can be used for correction of the atmospheric effects. For the simulation over ocean, the lowest error produced in this region was of the order of 105 and up to 14 times smaller than errors in the visible region. For the same spectral region on the forest case, the lowest error produced was of the order of 10-4, and up to 41 times smaller than errors in the visible region,

Gonzales, Fabian O.↗

Corrections to official ENDF/B Evaluation Releases for SCALE Nuclear Data

The official ENDF/B nuclear data evaluation releases in the past two decades have incrementally incorporated more detailed information, new nuclide evaluations, and very often have improved the accuracy of radiation transport codes when modeling shielding, fission reactors, criticality benchmarks, and fusion systems. However, like any large collaborative data compilation, these releases have all included small errors. This report documents the small corrections to ENDF/B releases ENDF/B-VII.1, ENDF/B-VIII.0, and ENDF/B-VIII.1 that have been applied during nuclear data processing to produce data libraries for the SCALE code system.

Brown, Jesse M. [Oak Ridge National Laboratory (OR↗

Thermodynamic Limits of Redox-Based Thermochemical Processes (REDOTHERM)

Solar thermochemical fuel production is a potential pathway for the production of sustain liquid drop-in fuels, which can help decarbonize the aviation and maritime sectors. In an attempt to analyze the commercial viability of this technology, several studies have been conducted, including system and technoeconomic analysis (TEA) modeling. However, most studies to date simply assume a given redox reactor efficiency, which is significantly higher than demonstrated values to date. While it is widely recognized that utilizing a counter-current flow (CF) configuration could increase the redox reactor efficiency, an over-simplification in the thermodynamic modeling may lead to unphysical results which has been included in multiple publications. The fact that the solar redox reactor is the least developed component in the process chain makes it hard to identify technology gaps and evaluate pathways to deployment at scale using this approach. In this work, a thermodynamic model for a moving oxide system has been developed, in a general form that allows to analyze the system for different redox-active materials, under a wide range of operating conditions, for both parallel and countercurrent flows. The model capabilites are demonstrated, and the model's code will be shared as an open-source on GitHub in the next few months.

chemical looping↗

The investigation of tethered satellite system dynamics

The tether control law to retrieve the satellite was modified in order to have a smooth retrieval trajectory of the satellite that minimizes the thruster activation. The satellite thrusters were added to the rotational dynamics computer code and a preliminary control logic was implemented to simulate them during the retrieval maneuver. The high resolution computer code for modelling the three dimensional dynamics of untensioned tether, SLACK3, was made fully operative and a set of computer simulations of possible tether breakages was run. The distribution of the electric field around an electrodynamic tether in vacuo severed at some length from the shuttle was computed with a three dimensional electrodynamic computer code.

Lorenzini, E.↗

Software for Collaborative Engineering of Launch Rockets

The Rocket Evaluation and Cost Integration for Propulsion and Engineering software enables collaborative computing with automated exchange of information in the design and analysis of launch rockets and other complex systems. RECIPE can interact with and incorporate a variety of programs, including legacy codes, that model aspects of a system from the perspectives of different technological disciplines (e.g., aerodynamics, structures, propulsion, trajectory, aeroheating, controls, and operations) and that are used by different engineers on different computers running different operating systems. RECIPE consists mainly of (1) ISCRM a file-transfer subprogram that makes it possible for legacy codes executed in their original operating systems on their original computers to exchange data and (2) CONES an easy-to-use filewrapper subprogram that enables the integration of legacy codes. RECIPE provides a tightly integrated conceptual framework that emphasizes connectivity among the programs used by the collaborators, linking these programs in a manner that provides some configuration control while facilitating collaborative engineering tradeoff studies, including design to cost studies. In comparison with prior collaborative-engineering schemes, one based on the use of RECIPE enables fewer engineers to do more in less time.

Stanley, Thomas Troy↗

A coupled radiative-convective-photochemical model of the stratosphere

A one-dimensional atmospheric model is described. The neutral photochemical processes of the troposphere and stratosphere are presented by a reaction set of 87 kinetic processes and 28 photolytic reactions, with vertical transport being taken into account by an eddy diffusion parameter. The vertical temperature profiles are coupled with the chemical calculations by means of a radiative-convective code. The model has the capability of including effects of multiple scattering and diurnal variations in species concentrations. The numerical solution methodology is described in detail, and current input data are given along with results obtained in a sample calculation.

Boughner, R. E.↗

Modelling of capillary optics as a focussing hard X-ray concentrator

The development of glass microcapillary arrays by Kumakhov and co-workers in the Soviet Union suggests that it is possible to construct novel devices for concentrating and focussing X-ray beams. Hypothetical focussing concentrators made up of capillary arrays have been configured to operate over a broad range of X-ray energy and its characteristics have been studied with a Monte Carlo ray tracing code. These modeling results indicate that the concentrators would have significant effective area up to 100 keV while their dimensions are similar to conventional grazing incidence telescopes that operate in the soft X-ray band.

Gorenstein, P.↗

Numerical Predictions of Damage and Failure in Carbon Fiber Reinforced Laminates Using a Thermodynamically-Based Work Potential Theory

A thermodynamically-based work potential theory for modeling progressive damage and failure in fiber-reinforced laminates is presented. The current, multiple-internal state variable (ISV) formulation, referred to as enhanced Schapery theory (EST), utilizes separate ISVs for modeling the effects of damage and failure. Consistent characteristic lengths are introduced into the formulation to govern the evolution of the failure ISVs. Using the stationarity of the total work potential with respect to each ISV, a set of thermodynamically consistent evolution equations for the ISVs are derived. The theory is implemented into a commercial finite element code. The model is verified against experimental results from two laminated, T800/3900-2 panels containing a central notch and different fiber-orientation stacking sequences. Global load versus displacement, global load versus local strain gage data, and macroscopic failure paths obtained from the models are compared against the experimental results.

composite materials↗

Space weather observations, modeling, and alerts in support of human exploration of Mars

Space weather observations and modeling at Mars have begun but they must be significantly increased to support the future of Human Exploration on the Red Planet. A comprehensive space weather understanding of a planet without a global magnetosphere and a thin atmosphere is very different from our situation at Earth so there is substantial fundamental research remaining. It is expected that the development of suitable models will lead to a comprehensive operational Mars space weather alert (MSWA) system that would provide rapid dissemination of information to Earth controllers, astronauts in transit, and those in the exploration zone (EZ) on the surface by producing alerts that are delivered rapidly and are actionable. To illustrate the importance of such a system, we use a magnetohydrodynamic code to model an extreme Carrington-type coronal mass ejection (CME) event at Mars. The results show a significant induced surface field of nearly 3,000 nT on the dayside that could radically affect unprotected electrical systems that would dramatically impact human survival on Mars. Other associated problems include coronal mass ejection (CME) shock-driven acceleration of solar energetic particles producing large doses of ionizing radiation at the Martian surface. In summary, along with working more closely with international partners, the next Heliophysics Decadal Survey must include a new initiative to meet expected demands for space weather forecasting in support of humans living and working on the surface of Mars. It will require significant effort to coordinate NASA and the international community contributions.

79 ASTRONOMY AND ASTROPHYSICS↗

Optimizing Price-Informed Operation of a Battery Storage System in an Office Building

New prescriptive efficiency requirements amended to US model energy codes historically have been evaluated using an average, blended electricity rate, which obscures demand charges. Post 2019, new measures can be evaluated using a representative time-of-use (TOU) tariff yet more sophisticated analysis methods are needed to consider a variety of TOU tariffs and assess price-informed control. To address these needs, this study couples building prototype simulation model with varying-in-sophistication battery storage operating strategies for different electricity TOU tariffs. The analysis compares the impact of operating a battery storage system following simpler rule-of-thumb methods versus a semi-optimized priced-informed heuristic approach. The investigation demonstrates that the heuristic approach results in greater electricity cost savings and is practical to implement.

Lei, Xuechen↗

An integrated optimum design approach for high speed prop rotors

The objective is to develop an optimization procedure for high-speed and civil tilt-rotors by coupling all of the necessary disciplines within a closed-loop optimization procedure. Both simplified and comprehensive analysis codes are used for the aerodynamic analyses. The structural properties are calculated using in-house developed algorithms for both isotropic and composite box beam sections. There are four major objectives of this study. (1) Aerodynamic optimization: The effects of blade aerodynamic characteristics on cruise and hover performance of prop-rotor aircraft are investigated using the classical blade element momentum approach with corrections for the high lift capability of rotors/propellers. (2) Coupled aerodynamic/structures optimization: A multilevel hybrid optimization technique is developed for the design of prop-rotor aircraft. The design problem is decomposed into a level for improved aerodynamics with continuous design variables and a level with discrete variables to investigate composite tailoring. The aerodynamic analysis is based on that developed in objective 1 and the structural analysis is performed using an in-house code which models a composite box beam. The results are compared to both a reference rotor and the optimum rotor found in the purely aerodynamic formulation. (3) Multipoint optimization: The multilevel optimization procedure of objective 2 is extended to a multipoint design problem. Hover, cruise, and take-off are the three flight conditions simultaneously maximized. (4) Coupled rotor/wing optimization: Using the comprehensive rotary wing code CAMRAD, an optimization procedure is developed for the coupled rotor/wing performance in high speed tilt-rotor aircraft. The developed procedure contains design variables which define the rotor and wing planforms.

Chattopadhyay, Aditi↗

Computational model for predicting the performance of fluid-powered longiudinal piston modules and motors.

A Matlab-based computational model line code has been developed for predicting the performance of fluid-powered longitudinal piston modules and motors. A system of equations has been developed based upon applying the conservation of mass to a control volume defined by a reciprocating piston. The method results in a differential equation that predicts the pressure gradient in the control volume as a function of rotor rotation angle. The model uses an integration scheme to predict the pressure state in the control volume. The orifice equation is used to predict the flowrate in and out of the control volume. From the resulting solution the torque and speed characteristics of the motor configuration can be predicted, optimized, or tailored for a specific requirement.SAND2019-13983 M Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Taylor, Michael↗

Influence of Antarctic and Greenland Continental Shelf Circulation on High‐Latitude Oceans in E3SM

The science objectives of this project are to simulate and understand the impacts of both deep-basin warm-water intrusions and land-ice melt on the continental shelf circulations and sea-ice distributions around the margins of Greenland and Antarctica. As well, the role of subsurface ocean heat from the Atlantic on declining sea-ice cover in the Arctic is explored. Mesoscale processes and fine bathymetry are implicated in cross-shelf property transports around both Greenland and Antarctica. Therefore, we configured and ran an atmospheric reanalysis-forced global ocean/sea-ice simulation on a grid that reduces from 8 km at the Equator to 2 km at the poles (UH8to2) with 60 vertical levels. It was produced using the Energy Exascale Earth System Model ‘‘HiLAT’’ code (E3SMv0-HiLAT) that uses the Parallel Ocean Program (POP) and CICE5 as its ocean and sea-ice components, respectively. Two main UH8to2 simulations were carried out: one for 1975-2009 and the other for July 2016-2020 after it was initialized from a 1/25° data-assimilative ocean/sea-ice prediction system ocean/sea-ice state. The UH8to2 is not coupled to an active land-ice model. Rather, land-ice melt is represented by observationally informed freshwater fluxes (FWFs). Short (multi-year) UH8to2 simulations were conducted to understand sensitivities when Greenland ice sheet (GrIS) melt is released only at the ocean surface or when it is distributed over the upper water column in accordance with fjord melt plume behavior; these cases were compared with a no GrIS melt case. West Greenland continental shelf currents were fastest in the vertical distribution case and an increase in baroclinic conversion at the shelf break associated with increased eddy kinetic energy was found relative to the surface release case. Further, salinity is lower and meltwater volume greater in the eastern Labrador Sea in the vertical distribution case. For the Arctic, the veracity of the UH8to2 was evaluated for 2017-2020 using available observations. Simulated seasonal sea-ice thickness and concentration are realistic, but the ice is unrealistically thin in the central and eastern Arctic in the fall. Comparisons of vertical sections of ocean temperature, salinity, and buoyancy collected from Ice-Tethered Profilers (ITPs) in the eastern Arctic in the fall and winter of 2019/2020 and co-located/concurrent UH8to2 fields show the stratification over the top 100 m of the water column is too low in the model, the simulated mixed layer too deep, and the simulated subsurface Atlantic Water (AW) too warm; these biases may contribute to the sea-ice biases. A model intercomparison study using the UH8to2 and a forced 1/25° regional Arctic ocean/sea-ice (uses the HYbrid Coordinate Ocean Model and CICE5) simulation further investigates the relationship between AW and sea-ice in the eastern Arctic. The models show a mesoscale-rich pulse of Atlantic Water extending into the eastern basin that reaches maximum intensity in late winter of 2018, after which it decreases in strength. Concurrent and co-located sea-ice melt or the inhibition of sea-ice growth is seen and is attributed to halocline mesoscale eddies doming into the mixed layer with convection bringing this heat into the vicinity of the sea-ice.

58 GEOSCIENCES↗

AGR-5/6/7 Final Release-to-Birth Ratio Data Analysis

AGR-5/6/7 is the last of a series of Advanced Gas Reactor (AGR) experiments conducted in the Advanced Test Reactor (ATR) at Idaho National Laboratory (INL) in support of development and qualification of tristructural isotropic (TRISO) low-enriched fuel for use in the high-temperature gas cooled reactor (HTGR). AGR configuration and irradiation conditions are based on prismatic HTGR technology that is distinguished primarily through the use of helium coolant, a low-power-density ceramic core capable of withstanding very high temperatures, and TRISO-coated particle fuel. The AGR tests provide valuable irradiation-performance data to support fuel process development, qualify fuel for normal operation and accident conditions, and support development and validation of fuel performance and fission-product (FP) transport models and codes. Each AGR test consists of multiple independently controlled and monitored capsules containing fuel compacts placed in a graphite cylinder shrouded by a steel shell. Release-to-birth ratios (R/B) for fission-gas isotopes released from each capsule are calculated from release rates, measured by germanium detectors in the Fission Product Monitoring System (FPMS) installed downstream from each capsule, and birth rates calculated using numerical models of FG generation. The R/Bs are a critical measure of the ability of the fuel kernel, the particle coating layers, and the compact matrix to retain fission-gas atoms, preventing their release into the sweep-gas flow, and the impact of initially defective particles and/or particle-coating failures that occur during irradiation. For fission-gas isotopes, particle failure is defined as failure of all coating layers, allowing gaseous fission atoms to escape from a particle. During the first five cycles (162B ? 165A), R/Bs were stable in the 10-8?10-6 range, and no in-pile particle failures were observed, based on the gross gamma counts. The maximum R/B value of around 2 ? 10-6 for Kr-85m resulted from the presence of as-fabricated exposed kernels (based on the high exposed kernel fraction), the dispersed uranium, and high fuel particle temperatures in Capsule 1. Comparison of capsule-measured R/Bs from these early cycles to predictions using the previously developed AGR R/B model demonstrated FG release from the AGR-5/6/7 TRISO fuel was comparable to that of previous experiments. In addition, the Kr-85m R/B per-exposed-kernel values are comparable to R/B values obtained in AGR-3/4 irradiation experiment and four irradiation experiments performed during 1980s: (1) HRB-17/18, (2) COMEDIE-BD1, (3) HFR B1, and (4) HRB-21. In contrast, all measured R/B values are lower than predictions by the commonly used Richards and German models, which are intentionally conservative. A large number of in-pile particle failures occurred in Capsule 1 by the end of Cycle 166A. During the final four cycles (166A ? 168A), apparent damage to the Capsule 1 gas line appeared to cause FG leakage from that capsule into the other four capsules, resulting in an increase in fission gas (FG) detected in the effluent for all capsules. Isolation of the Capsule 1 gas line during the last three cycles also prevented measurement of its FG release. Thus, R/Bs in all capsules after Cycle 166A are highly uncertain because of undefined amount of leakage from Capsule 1, especially for long-lived isotopes. A few hundred in-pile particle failures were estimated for Capsule 1 before the end of Cycle 166A, but the total number of failures is unknown due to the lack of FG release data in the later cycles. Based primarily on evidence from the gross gamma counts during Cycle 168A, approximately 15 particles failed in Capsule 3 and four particles failed in Capsule 2. In-pile failures in Capsule 3 were anticipated because this capsule was designed to operate beyond the HTGR normal operating temperature range. In contrast, no in-pile failures were identified in the top two capsules (4 and 5) based on the absence of the typical spikes in gross gamma counts and low failure estimates using the AGR model, developed in INL/EXT-14-32970, for R/B of the short-lived isotopes (Kr-89 and Xe-137) with minimal leakage from Capsule 1.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

White Paper on Case Study of Safe Installation of Second-Life Energy Storage System

This technical report provides for a case study for the safe installation of a second-life, or repurposed, battery, that has been reconfigured for use as a stationary energy storage systems (ESS). Driven by legislative requirements such as California Senate Bill 615 and projections that retired EV batteries could meet a substantial portion of U.S. grid ESS needs beginning in 2035, the repurposing of EV batteries is anticipated to grow significantly. However, safety concerns arise from the effects of aging, unknown prior usage history, and changes in thermal runaway behavior, which may increase failure risks compared to new batteries. NFPA 855, the predominant U.S. standard for ESS installation, mandates that second-life batteries meet all requirements for new batteries, with repurposers complying with UL 1974 in addition to obtaining UL 9540 and UL 1973 listings. These are certifications that few repurposers have achieved and represent a regulatory barrier to entry for the market as a whole.

47 OTHER INSTRUMENTATION↗