SEARCH · Engineering Papers
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.
This note describes the basic levels of Hamiltonian approximation (for straight-axis elements) that are used in the symplectic beam dynamics modeling code ImpactX. These include: purely linear models, models based on a paraxial (chromatic) expansion of the Hamiltonian, and models based on the exact nonlinear Hamiltonian. The focus here is on paraxial (chromatic) models.
At a former uranium mill site where tailings have been removed, prior work has determined several potential ongoing secondary uranium sources. These include locations with uranium sorbed to organic carbon, uranium in the unsaturated zone, and uranium associated with the presence of gypsum. To better understand uranium mobility controls at the site, four single-well push–pull tests (with a drift phase) were completed with the goal of deriving aquifer flow and contaminant transport parameters for inclusion in a future sitewide reactive transport model. This goes beyond the traditional use of a constant sorption distribution coefficient (K d ) and allows for the evaluation of alternative remedial injection fluids, which can produce variable K d values. Dispersion was first removed from the resulting data to determine possible reactions before conducting reactive transport simulations. These initial analyses indicated the potential need to include cation exchange, uranium sorption, and gypsum dissolution. A reactive transport model using multiple layers to account for partially penetrating wells was completed using the PHT-USG reactive transport modeling code and calibrated using PEST. The model results quantify the hydraulic conductivity and dispersion parameters using the injected tracer concentrations. Uranium sorption, cation exchange, and gypsum dissolution parameters were quantified by comparing the simulated versus observed geochemistry. All simulations required some cation exchange and calcite equilibrium, and one simulation required gypsum dissolution to improve the model fit for calcium and sulfate. Uranium sorption parameters were not strongly influenced by the other parameter values but were highly influenced by uranium concentrations during the drift phase, with possible kinetic rate limitations. Thus, a future recommendation for such push–pull tests is to collect more geochemical data during the drift phase. The final uranium sorption parameters were within the range of values determined from prior column testing. The flow and transport parameters derived from these single-well push–pull tests will provide initial parameters for any future sitewide reactive transport model.
Fuel performance modeling status update and potential model improvements overview of TRISO fuel performance modeling codes PARFUME/BISON, AGR experiment support, potential modeling improvements, BISON smeared cracking model, BISON fission product source term, and AGR-3/4 reirradiation heating test.
Numerical modeling of the Pulsed Inductive Thruster exercising the magnetohydrodynamics code, MACH2 aims to provide bilateral validation of the thruster's measured performance and the code's capability of capturing the pertinent physical processes. Computed impulse values for helium and argon propellants demonstrate excellent correlation to the experimental data for a range of energy levels and propellant-mass values. The effects of the vacuum tank wall and massinjection scheme were investigated to show trivial changes in the overall performance. An idealized model for these energy levels and propellants deduces that the energy expended to the internal energy modes and plasma dissipation processes is independent of the propellant type, mass, and energy level.
A study has been performed on a computer code modeling an aircraft wake vortex spacing system during final approach. This code represents an initial engineering model of a system to calculate reduced approach separation criteria needed to increase airport productivity. This report evaluates model sensitivity toward various weather conditions (crosswind, crosswind variance, turbulent kinetic energy, and thermal gradient), code configurations (approach corridor option, and wake demise definition), and post-processing techniques (rounding of provided spacing values, and controller time variance).
The Honeywell Uncertified Research Engine (HURE), a research version of a turbofan engine that never entered production, was tested in the NASA Propulsion System Laboratory (PSL), an altitude test facility at the NASA Glenn Research Center. The PSL is a facility that is equipped with water spray bars capable of producing an ice cloud consisting of ice particles, having a controlled particle diameter and concentration in the air flow. To develop the test matrix of the HURE, numerical analysis of flow and ice particle thermodynamics was performed on the compression system of the turbofan engine to predict operating conditions that could potentially result in a risk of ice accretion due to ice crystal ingestion. The goal of the test matrix was to have ice accrete in two regions of the compression system: region one, which consists of the fan-stator through the inlet guide vane (IGV), and region two which is the first stator within the high pressure compressor. The predictive analyses were performed with the mean line compressor flow modeling code (COMDES-MELT) which includes an ice particle model. The HURE engine was tested in PSL with the ice cloud over the range of operating conditions of altitude, ambient temperature, simulated flight Mach number, and fan speed with guidance from the analytical predictions. The engine was fitted with video cameras at strategic locations within the engine compression system flow path where ice was predicted to accrete, in order to visually confirm ice accretion when it occurred. In addition, traditional compressor instrumentation such as total pressure and temperature probes, static pressure taps, and metal temperature thermocouples were installed in targeted areas where the risk of ice accretion was expected. The current research focuses on the analysis of the data that was obtained after testing the HURE engine in PSL with ice crystal ingestion. The computational method (COMDES-MELT) was enhanced by computing key parameters through the fan-stator at multiple span wise locations, in order to increase the fidelity with the current mean-line method. The Icing Wedge static wet bulb temperature thresholds were applicable for determining the risk of ice accretion in the fan-stator, which is thought to be an adiabatic region. At some operating conditions near the splitter-lip region, other sources of heat (non-adiabatic walls) were suspected to be the cause of accretion, and the Icing Wedge was not applicable to predict accretion at that location. A simple order-of-magnitude heat transfer model was implemented into the COMDES-MELT code to estimate the wall temperature minimum and maximum thresholds that support ice accretion, as observed by video confirmation. The results from this model spanned the range of wall temperatures measured on a previous engine that experienced ice accretion at certain operating conditions.
The Honeywell Uncertified Research Engine (HURE), a research version of a turbofan engine that never entered production, was tested in the NASA Propulsion System Laboratory (PSL), an altitude test facility at the NASA Glenn Research Center. The PSL is a facility that is equipped with water spray bars capable of producing an ice cloud consisting of ice particles, having a controlled particle diameter and concentration in the air flow. To develop the test matrix of the HURE, numerical analysis of flow and ice particle thermodynamics was performed on the compression system of the turbofan engine to predict operating conditions that could potentially result in a risk of ice accretion due to ice crystal ingestion. The goal of the test matrix was to have ice accrete in two regions of the compression system: region one, which consists of the fan-stator through the inlet guide vane (IGV), and region two which is the first stator within the high pressure compressor. The predictive analyses were performed with the mean line compressor flow modeling code (COMDES-MELT) which includes an ice particle model. The HURE engine was tested in PSL with the ice cloud over the range of operating conditions of altitude, ambient temperature, simulated flight Mach number, and fan speed with guidance from the analytical predictions. The engine was fitted with video cameras at strategic locations within the engine compression system flow path where ice was predicted to accrete, in order to visually confirm ice accretion when it occurred. In addition, traditional compressor instrumentation such as total pressure and temperature probes, static pressure taps, and metal temperature thermocouples were installed in targeted areas where the risk of ice accretion was expected. The current research focuses on the analysis of the data that was obtained after testing the HURE engine in PSL with ice crystal ingestion. The computational method (COMDES-MELT) was enhanced by computing key parameters through the fan- stator at multiple span wise locations, in order to increase the fidelity with the current mean-line method. The Icing Wedge static wet bulb temperature thresholds were applicable for determining the risk of ice accretion in the fan-stator, which is thought to be an adiabatic region. At some operating conditions near the splitter-lip region, other sources of heat (non-adiabatic walls) were suspected to be the cause of accretion, and the Icing Wedge was not applicable to predict accretion at that location. A simple order-of-magnitude heat transfer model was implemented into the COMDES-MELT code to estimate the wall temperature minimum and maximum thresholds that support ice accretion, as observed by video confirmation. The results from this model spanned the range of wall temperatures measured on a previous engine that experienced ice accretion at certain operating conditions.
This Modeling Archive is in support the publication “Intensified Soil Moisture Extremes Decrease Soil Organic Carbon Decomposition: A Mechanistic Modeling Analysis” (Liang et al., 2021). Here we provide model code, inputs, outputs and evaluation datasets for the Microbial ENzyme Decomposition (MEND) model for the Missouri Ozarks AmeriFlux eddy covariance measurement site (MOFLUX) near Ashland, Missouri USA. The MEND model was developed with explicit representation of microbial and enzyme pools to mechanistically simulate the role of microbial organisms and extracellular enzymes in soil organic carbon (SOC) decomposition. Long-term SOC dynamics under intensified moisture extremes are studied using the MEND model that is parameterized with 11 years of measurements from the MOFLUX forest. The model explicitly represents microbial dormancy and resuscitation, different types of SOC-degrading enzymes, and how they vary with changes in soil moisture (Wang et al. 2015, 2019). A combination of two levels of frequency and severity of soil moisture, as well as a control with normal interannual variability, are used to simulate a range of moisture scenarios over 100 years. The code of Microbial-ENzyme Decomposition (MEND) as well as the input and output data are included in the archive. A user’s manual (MEND_Readme.pdf) is included with instructions for compiling and running the model to simulate soil organic carbon decomposition under various moisture scenarios. This dataset contains the modelling archive contained within a compressed (*.zip) file, a file-level metadata file in comma separate (*.csv) format, and two instructional files in PDF (*.pdf) format.
Mathematical models and computer codes based on these models, which allow prediction of the product distribution in chemical reactors for converting gaseous silicon compounds to condensed-phase silicon were developed. The following tasks were accomplished: (1) formulation of a model for silicon vapor separation/collection from the developing turbulent flow stream within reactors of the Westinghouse (2) modification of an available general parabolic code to achieve solutions to the governing partial differential equations (boundary layer type) which describe migration of the vapor to the reactor walls, (3) a parametric study using the boundary layer code to optimize the performance characteristics of the Westinghouse reactor, (4) calculations relating to the collection efficiency of the new AeroChem reactor, and (5) final testing of the modified LAPP code for use as a method of predicting Si(1) droplet sizes in these reactors.
The purpose of this report is to summarize and document the work done to enable a NASA CFD code to model laminar-turbulent transition process on an isolated turbine blade. The ultimate purpose of the present work is to down-select a transition model that would allow the flow simulation of a variable speed power turbine to be accurately performed. The flow modeling in its final form will account for the blade row interactions and their effects on transition which would lead to accurate accounting for losses. The present work only concerns itself with steady flows of variable inlet turbulence. The low Reynolds number k- model of Wilcox and a modified version of the same model will be used for modeling of transition on experimentally measured blade pressure and heat transfer. It will be shown that the k- model and its modified variant fail to simulate the transition with any degree of accuracy. A case is thus made for the adoption of more accurate transition models. Three-equation models based on the work of Mayle on Laminar Kinetic Energy were explored. The three-equation model of Walters and Leylek was thought to be in a relatively mature state of development and was implemented in the Glenn-HT code. Two-dimensional heat transfer predictions of flat plate flow and two-dimensional and three-dimensional heat transfer predictions on a turbine blade were performed and reported herein. Surface heat transfer rate serves as sensitive indicator of transition. With the newly implemented model, it was shown that the simulation of transition process is much improved over the baseline k- model for the single Reynolds number and pressure ratio attempted; while agreement with heat transfer data became more satisfactory. Armed with the new transition model, total-pressure losses of computed three-dimensional flow of E3 tip section cascade were compared to the experimental data for a range of incidence angles. The results obtained, form a partial loss bucket for the chosen blade. In time the loss bucket will be populated with losses at additional incidences. Results obtained thus far will be discussed herein.
Evolutionary solar models are computed using a new stellar evolution code, MOSEC (Modular Stellar Evolution Code). This code has been designed with carefully controlled truncation errors in order to achieve a precision which reflects the increasingly accurate determination of solar interior structure by helioseismology. A series of models is constructed to investigate the effects of the choice of equation of state (OPAL or MHD-E, the latter being a version of the MHD equation of state recalculated by the author), the inclusion of helium and heavy-element settling and diffusion, and the inclusion of a simple model of mixing associated with the solar tachocline. The neutrino flux predictions are discussed, while the sound speed of the computed models is compared to that of the sun via the latest inversion of SOI-NMI p-mode frequency data. The comparison between models calculated with the OPAL and MHD-E equations of state is particularly interesting because the MHD-E equation of state includes relativistic effects for the electrons, whereas neither MHD nor OPAL do. This has a significant effect on the sound speed of the computed model, worsening the agreement with the solar sound speed. Using the OPAL equation of state and including the settling and diffusion of helium and heavy elements produces agreement in sound speed with the helioseismic results to within about +.-0.2%; the inclusion of mixing slightly improves the agreement.
The ionosphere of Mars is rather poorly understood. The only direct measurements were performed by the Viking 1 and 2 landers in 1976, both of which carried a Retarding Potential Analyzer. The RPA was designed to measure ion properties during the descent, although electron fluxes were estimated from changes in the ion currents. Using these derived low-energy electron fluxes, Mantas and Hanson studied the photoelectron and the solar wind electron interactions with the atmosphere and ionosphere of Mars. Unanswered questions remain regarding the origin of the low-energy electron fluxes in the vicinity of the Mars plasma boundary. Crider, in an analysis of Mars Global Surveyor Magnetometer/Electron Reflectometer measurements, has attributed the formation of the magnetic pile-up boundary to electron impact ionization of exospheric neutral species by solar wind electrons. However, the role of photoelectrons escaping from the lower ionosphere was not determined. In the proposed work, we will examine the role of solar wind and ionospheric photoelectrons in producing ionization in the upper ionosphere of Mars. Low-energy (< 4 keV) electrons will be modeled using the two-stream electron transport code of Link. The code models both external (solar wind) and internal (photoelectron) sources of ionization, and accounts for Auger electron production. The code will be used to analyze Mars Global Surveyor measurements of solar wind and photoelectrons down to altitudes below 200 km in the Mars ionosphere, in order to determine the relative roles of solar wind and escaping photoelectrons in maintaining plasma densities in the region of the Mars plasma boundary.
Mathematical models, and computer codes based on these models were developed which allow prediction of the product distribution in chemical reactors in which gaseous silicon compounds are converted to condensed phase silicon. The reactors to be modeled are flow reactors in which silane or one of the halogenated silanes is thermally decomposed or reacted with an alkali metal, H2 or H atoms. Because the product of interest is particulate silicon, processes which must be modeled, in addition to mixing and reaction of gas-phase reactants, include the nucleation and growth of condensed Si via coagulation, condensation, and heterogeneous reaction.
Modeling guidelines constitute one of the fundamental cornerstones for Model Based Development. Their relevance is essential when dealing with code generation in the safety-critical domain. This article presents the experience of a railway signaling systems manufacturer on this issue. Introduction of Model-Based Development (MBD) and code generation in the industrial safety-critical sector created a crucial paradigm shift in the development process of dependable systems. While traditional software development focuses on the code, with MBD practices the focus shifts to model abstractions. The change has fundamental implications for safety-critical systems, which still need to guarantee a high degree of confidence also at code level. Usage of the Simulink/Stateflow platform for modeling, which is a de facto standard in control software development, does not ensure by itself production of high-quality dependable code. This issue has been addressed by companies through the definition of modeling rules imposing restrictions on the usage of design tools components, in order to enable production of qualified code. The MAAB Control Algorithm Modeling Guidelines (MathWorks Automotive Advisory Board)[3] is a well established set of publicly available rules for modeling with Simulink/Stateflow. This set of recommendations has been developed by a group of OEMs and suppliers of the automotive sector with the objective of enforcing and easing the usage of the MathWorks tools within the automotive industry. The guidelines have been published in 2001 and afterwords revisited in 2007 in order to integrate some additional rules developed by the Japanese division of MAAB [5]. The scope of the current edition of the guidelines ranges from model maintainability and readability to code generation issues. The rules are conceived as a reference baseline and therefore they need to be tailored to comply with the characteristics of each industrial context. Customization of these recommendations has been performed for the automotive control systems domain in order to enforce code generation [7]. The MAAB guidelines have been found profitable also in the aerospace/avionics sector [1] and they have been adopted by the MathWorks Aerospace Leadership Council (MALC). General Electric Transportation Systems (GETS) is a well known railway signaling systems manufacturer leading in Automatic Train Protection (ATP) systems technology. Inside an effort of adopting formal methods within its own development process, GETS decided to introduce system modeling by means of the MathWorks tools [2], and in 2008 chose to move to code generation. This article reports the experience performed by GETS in developing its own modeling standard through customizing the MAAB rules for the railway signaling domain and shows the result of this experience with a successful product development story.
The Environment Power System Analysis Tool (EPSAT) is being developed to provide engineers with the ability to assess the effects of a broad range of environmental interactions on space power systems. A unique user-interface-data-dictionary code architecture oversees a collection of existing and future environmental modeling codes (e.g., neutral density) and physical interaction models (e.g., sheath ionization). The user-interface presents the engineer with tables, graphs, and plots which, under supervision of the data dictionary, are automatically updated in response to parameter change. EPSAT thus provides the engineer with a comprehensive and responsive environmental assessment tool and the scientist with a framework into which new environmental or physical models can be easily incorporated.
The PI (Cranmer) and Co-I (A. van Ballegooijen) made significant progress toward the goal of building a "unified model" of the dominant physical processes responsible for the acceleration of the solar wind. The approach outlined in the original proposal comprised two complementary pieces: (1) to further investigate individual physical processes under realistic coronal and solar wind conditions, and (2) to extract the dominant physical effects from simulations and apply them to a one-dimensional and time-independent model of plasma heating and acceleration. The accomplishments in the report period are thus divided into these two categories: 1a. Focused Study of Kinetic MHD Turbulence. We have developed a model of magnetohydrodynamic (MHD) turbulence in the extended solar corona that contains the effects of collisionless dissipation and anisotropic particle heating. A turbulent cascade is one possible way of generating small-scale fluctuations (easy to dissipate/heat) from a pre-existing population of low-frequency Alfven waves (difficult to dissipate/heat). We modeled the cascade as a combination of advection and diffusion in wavenumber space. The dominant spectral transfer occurs in the direction perpendicular to the background magnetic field. As expected from earlier models, this leads to a highly anisotropic fluctuation spectrum with a rapidly decaying tail in the parallel wavenumber direction. The wave power that decays to high enough frequencies to become ion cyclotron resonant depends on the relative strengths of advection and diffusion in the cascade. For the most realistic values of these parameters, though, there is insufficient power to heat protons and heavy ions. The dominant oblique waves undergo Landau damping, which implies strong parallel electron heating. We thus investigated the nonlinear evolution of the electron velocity distributions (VDFs) into parallel beams and discrete phase-space holes (similar to those seen in the terrestrial magnetosphere) which are an alternate means of heating protons via stochastic interactions similar to particle-particle collisions. 1b. Focused Study of the Multi-Mode Detailed Balance Formalism. The PI began to explore the feasibility of using the "weak turbulence," or detailed-balance theory of Tsytovich, Melrose, and others to encompass the relevant physics of the solar wind. This study did not go far, however, because if the "strong" MHD turbulence discussed above is a dominant player in the wind's acceleration region, this formalism is inherently not applicable to the corona. We will continue to study the various published approaches to the weak turbulence formalism, especially with an eye on ways to parameterize nonlinear wave reflection rates. 2. Building the Unified Model Code Architecture. We have begun developing the computational model of a time-steady open flux tube in the extended corona. The model will be "unified" in the sense that it will include (simultaneously for the first time) as many of the various proposed physical processes as possible, all on equal footing. To retain this generality, we have formulated the problem in two interconnected parts: a completely kinetic model for the particles, using the Monte Carlo approach, and a finite-difference approach for the self-consistent fluctuation spectra. The two codes are run sequentially and iteratively until complete consistency is achieved. The current version of the Monte Carlo code incorporates gravity, the zero-current electric field, magnetic mirroring, and collisions. The fluctuation code incorporates WKJ3 wave action conservation and the cascade/dissipation processes discussed above. The codes are being run for various test problems with known solutions. Planned additions to the codes include prescriptions for nonlinear wave steepening, kinetic velocity-space diffusion, and multi-mode coupling (including reflection and refraction).
Blading models added to TURBAN, an axial turbine meanline design code, and CSPAN, an axial compressor spanline design code, are presented. Calculations of blade chord length, blade axial length, and number of blades for each blade row are now provided by each code. In order to accomplish this, computational models for axial solidity and stagger angle were added to TURBAN. For a compressor, where solidity and chord were already available from input, the modeling additions to CSPAN were for stagger angle and for the blade angles to include incidence and deviation. All added computations utilized models already available in the literature.
The U.S. Department of Energy (DOE) Building Energy Codes Program (BECP) periodically evaluates state and national impacts associated with energy codes in residential and commercial buildings. Pacific Northwest National Laboratory (PNNL), funded by DOE, assessed the prospective impacts of national model building energy codes from 2010 through 2040. The model codes of interest in this report are the International Energy Conservation Code (IECC) for residential buildings and ASHRAE Standard 90.1 for commercial buildings. This report describes the methodology used for the assessment and presents the impacts in terms of energy savings, consumer cost savings, and reduced CO 2 emissions at the state and national levels.