Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Modelling Capabilities”

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 649 records · Page 36

Introducing the 9500 Node Distribution Test System to Support Advanced Power Applications: An Operations-Focused Approach

The 9500 Node Test System is a representative section of distribution power system model developed as a part of the GridAPPS-D™ project, an effort funded by DOE as a part of the Grid Modernization Lab Consortium (GMLC) program. The test system was developed to fulfill a growing need to represent the rapidly evolving state of electric distribution systems by combining elements of legacy infrastructure systems, modern feeder topologies, and an anticipated future with smart grid technologies. It also provides a network model capable of supporting the simulation of operational scenarios such as the ones in a utility distribution control center. This test system allows the evaluation of the performance of advanced power applications in real-time, such as one that simulates the operations of an Advanced Distribution Management Systems (ADMS), Distributed Energy Resource Management Systems (DERMS), etc. in a Distribution control center. This model is an extension of the widely used IEEE 8500 Node Test Feeder and is currently being validated to become an IEEE test case to help increase adoption and widespread usage among both academia and industry. It is a full-size model representative of a section of a utility’s distribution system with multiple feeders fed from different substations. The model includes multiple distribution circuits, a sub-transmission system, multiple substations, behind the meter customer rooftop photovoltaics (PV), and multiple utility-scale distributed energy resources. To enable accurate simulations of operational scenarios, the 9500 Node Test System is designed to support procedure-based operations, with the ability to realistically demonstrate switching operations, feeder reconfiguration, adjustment of volt-var control equipment, dispatch of distributed generation, and response to planned and unplanned outages. The 9500 Node Test System includes three radial distribution feeders with 12.3 MVA of average load, consisting of both medium voltage and low voltage equipment each supplied by a different distribution substation. The three distribution feeders are connected to each other through Normally-Open switches which can be closed when needed to simulate restoration scenarios due to a fault. One feeder represents today’s grid with low penetration of customer-side renewables. The second represents a potential future grid with microgrids and 100% renewable penetration. The third has no customer generation resources, a district steam plant, and a utility-scale solar farm. The three diverse circuits were created to allow the simulation of both today’s situation as well as potential future scenarios. All three feeders have customers connected by low-voltage secondary triplex lines. This test system meets all requirements outlined in the report for creating a simulation environment that would enable discussion between key technical stakeholders as well as having the potential to accelerate operational application development and their subsequent testing and integration. The new model is a possible representation of what we believe the distribution grid may look like in the future: a high penetration of renewables, reconfigurable radial and mesh topology, numerous DERs, islanded microgrids, and significantly increased data and measurement density. The system supports both the solution of existing and newer algorithms but also enables the evaluation of applications in a realistic operational environment defined by task oriented procedural steps that represent the interaction between the control center operator and field personnel.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Code Coverage Status of the ARC Code PERSENT

The Argonne Reactor Code (ARC) software system supports users in their fast reactor design goals by providing neutronic, thermal-hydraulic, and structural analysis capabilities. PERSENT fulfills the role of generating reactivity coefficients for a given time point of a REBUS calculation usable in a point kinetics based safety analysis capability. PERSENT also provides a sensitivity coefficient capability on eigenvalue, reactivity worth, and several other key coefficients that are used in the follow-on safety analysis. Given a co-variance matrix, PERSENT can carry out the uncertainty quantification to indicate the amount of error in the reactivity coefficients derived from the errors in the cross section measurements. With continued improvement of computational resources, many of the geometry modeling capabilities in DIF3D that were primarily used in low order schemes are not really needed anymore. Today, the diffusion and transport capabilities of DIF3D-VARIANT are primarily used in the reactor design process with some scattered usage of DIF3D-FD and DIF3D-Nodal. PERSENT is part of the ARC code system and is built around DIF3D-VARIANT and the flux solution it provides. The purpose of the present work is to identify a set of test problems for PERSENT and assess the code coverage of PERSENT for those test problems. PERSENT treats the DIF3D executable as an external executable and thus the code coverage considerations only need to focus on the PERSENT source code and only a fraction of the connected modules in the existing ARC software library. The goal is to document what parts of the existing PERSENT code are touched by the set of test problems and which are not. Because the verification work done on PERSENT was focused on the most common uses of PERSENT for fast reactor analysis, the code coverage assessment of those capabilities is the highest priority. This will ensure that nothing is being missed by the existing verification test problems that users of PERSENT rely upon. The code coverage analysis of PERSENT was performed with the Code Coverage Tool of the Intel Fortran compiler which requires modifications to the compilation of PERSENT. The detailed coverage tables are given for each submodule of PERSENT. Most of the uncovered parts/files could be easily ignored because they are either for error message and debugging output or not needed by PERSENT today. Only a few uncovered parts of PERSENT deserve extending the verification test suite.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A comprehensive Guyton model analysis of physiologic responses to preadapting the blood volume as a countermeasure to fluid shifts

The Guyton model of fluid, electrolyte, and circulatory regulation is an extensive mathematical model capable of simulating a variety of experimental conditions. It has been modified for use at NASA to simulate head-down tilt, a frequently used analog of weightlessness. Weightlessness causes a headward shift of body fluids that is believed to expand central blood volume, triggering a series of physiologic responses resulting in large losses of body fluids. We used the modified Guyton model to test the hypothesis that preadaptation of the blood volume before weightless exposure could counteract the central volume expansion caused by fluid shifts, and thereby attenuate the circulatory and renal responses that result in body fluid losses. Simulation results show that circulatory preadaptation, by a procedure resembling blood donation immediately before head-down bedrest, is effective in damping the physiologic responses to fluid shifts and reducing body fluid losses. After 10 hours of head-down tilt, preadaptation also produces higher blood volume, extracellular volume, and total body water for 20 to 30 days of bedrest, compared with non-preadapted control. These results indicate that circulatory preadaptation before current Space Shuttle missions may be beneficial for the maintenance of reentry and postflight orthostatic tolerance in astronauts. This paper presents a comprehensive examination of the simulation results pertaining to changes in relevant physiologic variables produced by blood volume reduction before a prolonged head-down tilt. The objectives were to study and develop the countermeasure theoretically, to aid in planning experimental studies of the countermeasure, and to identify potentially disadvantageous physiologic responses that may be caused by the countermeasure.

NASA Center JSC↗

Simulation of the Fast Reactor Fuel Assembly Duct-Bowing Reactivity Effect Using Monte Carlo Neutron Transport and Finite Element Analysis

This paper discusses a new method of simulating the fuel assembly duct-bowing reactivity coefficient for EBR-II run 138B. Quantification of the fuel assembly duct-bowing reactivity effect in liquid metal–cooled fast reactors has been a persistent problem since they were first designed and operated. Simulation of the duct-bowing reactivity effect is difficult because the level of detail required to simulate the effect has exceeded most modeling capabilities. The new method outlined in this paper utilizes the finite element analysis code ANSYS to analyze the thermal and structural components. Here, the displacement of the fuel assembly duct due to thermal expansion and mechanical interaction was calculated by ANSYS using recorded EBR-II run 138B temperature and power boundary value data. The displacement values were incorporated into to a Monte Carlo model of EBR-II run 138B and keff was calculated. Multiple Monte Carlo calculations were performed with duct displacement values corresponding to different reactor temperatures. Using the calculated keff values associated with the different duct displacement results allowed calculation of the duct-bowing reactivity coefficient. The duct-bowing reactivity coefficient was calculated to be –14.5 × 10 –4 $/°C/ ± 4.4%.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Computational methods and software systems for dynamics and control of large space structures

The deployment, assembly and mission-oriented maneuvering of space structures in orbit will trigger large motions of flexible, truss-type structures. In addition, the presence of on-board controls both for attitude stabilization and specified vibration tolerance requirements may further complicate the dynamic behavior of the orbiting structures. Because of safety and cost considerations, the dynamic response of the combined structural and control systems must be predicted reliably. This need can only be met through the development of reliable and efficient simulation capabilities, since there is general agreement that on-orbit experiments should be limited because of cost, time and facility constraints. The long-term objective of this research effort is to develop a next-generation computer simulator for the dynamics and control of large space structures. The simulator will be based on integrating four research thrusts: a new multibody dynamics formulation methodology, modeling capabilities in long/slender truss-beam components with realistic joints, efficient computational procedures that can be implemented either in sequential or concurrent computers, and prototype simulation modules that can be easily processed into a modern large-scale engineering software system such as the NASA/Computational Structural Mechanics (CSM) testbed.

Park, K. C.↗

GoldSim E-Area Low-Level Waste Lacility Vadose Zone Model Benchmarking

This report documents the development and benchmarking of the E-Area Low-Level Waste Facility (ELLWF) GoldSim Engineered and Slit Trench (ET and ST) vadose zone models. Subsequent activities (beyond the scope of this effort) will couple the GoldSim-based vadose zone models to the GoldSim aquifer model, with the intent of performing stochastic analyses that couple all pertinent fate and transport processes from the ground surface up to the 100-m Point of Assessment (POA). The benchmarking presented within this report focuses on the Slit Trench 06 (ST06) bydrostratigraphic location that resides within the center set of slit trenches. These GoldSim-based vadose zone models are generic models capable of being benchmarked to any of the other trench locations within E-Area.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

The importance of freeze–thaw cycles for lateral tracer transport in ice-wedge polygons

Abstract. A significant portion of the Arctic coastal plain is classified as polygonal tundra and plays a vital role in soil carbon cycling. Recent research suggests that lateral transport of dissolved carbon could exceed vertical carbon releases to the atmosphere. However, the details of lateral subsurface flow in polygonal tundra have not been well studied. We incorporated a subsurface transport process into an existing state-of-the-art hydrothermal model. The model captures the physical effects of freeze–thaw cycles on lateral flow in polygonal tundra. The new modeling capability enables non-reactive tracer movement within subsurface. We utilized this new capability to investigate the impact of freeze–thaw cycles on lateral flow in the polygonal tundra. Our study indicates the important role of freeze–thaw cycles and the freeze-up effect in lateral tracer transport, suggesting that dissolved species could be transported from the middle of the polygon to the sides within a couple of thaw seasons. Introducing lateral carbon transport into the climate models could substantially reduce the uncertainty associated with the impact of thawing permafrost.

54 ENVIRONMENTAL SCIENCES↗

Accelerating FEM-Based Corrosion Predictions Using Machine Learning

Atmospheric corrosion of metallic parts is a widespread materials degradation phenomena that is challenging to predict given its dependence on many factors (e.g. environmental, physiochemical, and part geometry). For materials with long expected service lives, accurately predicting the degree to which corrosion will degrade part performance is especially difficult due to the stochastic nature of corrosion damage spread across years or decades of service. The Finite Element Method (FEM) is a computational technique capable of providing accurate estimates of corrosion rate by numerically solving complex differential Eqs. characterizing this phenomena. Nevertheless, given the iterative nature of FEM and the computational expense required to solve these complex equations, FEM is ill-equipped for an efficient exploration of the design space to identify factors that accelerate or deter corrosion, despite its accuracy. In this work, a machine learning based surrogate model capable of providing accurate predictions of corrosion with significant computational savings is introduced. Specifically, this work leverages AdaBoosted Decision trees to provide an accurate estimate of corrosion current per width given different values of temperature, water layer thickness, molarity of the solution, and the length of the cathode for a galvanic couple of aluminum and stainless steel.

36 MATERIALS SCIENCE↗

Improving Mars-GRAM: Increasing the Accuracy of Sensitivity Studies at Large Optical Depths

Extensively utilized for numerous mission applications, the Mars Global Reference Atmospheric Model (Mars-GRAM) is an engineering-level atmospheric model. In a Monte-Carlo mode, Mars-GRAM's perturbation modeling capability is used to perform high fidelity engineering end-to-end simulations for entry, descent, and landing (EDL). Mars-GRAM has been found to be inexact when used during the Mars Science Laboratory (MSL) site selection process for sensitivity studies for MapYear=0 and large optical depth values such as tau=3. Mars-GRAM is based on the NASA Ames Mars General Circulation Model (MGCM) from the surface to 80 km altitude. Mars-GRAM with the MapYear parameter set to 0 utilizes results from a MGCM run with a fixed value of tau=3 at all locations for the entire year. Imprecise atmospheric density and pressure at all altitudes is a consequence of this use of MGCM with tau=3. Density factor values have been determined for tau=0.3, 1 and 3 as a preliminary fix to this pressure-density problem. These factors adjust the input values of MGCM MapYear 0 pressure and density to achieve a better match of Mars-GRAM MapYear 0 with Thermal Emission Spectrometer (TES) observations for MapYears 1 and 2 at comparable dust loading. These density factors are fixed values for all latitudes and Ls and are included in Mars-GRAM Release 1.3. Work currently being done, to derive better multipliers by including variations with latitude and/or Ls by comparison of MapYear 0 output directly against TES limb data, will be highlighted in the presentation. The TES limb data utilized in this process has been validated by a comparison study between Mars atmospheric density estimates from Mars-GRAM and measurements by Mars Global Surveyor (MGS). This comparison study was undertaken for locations on Mars of varying latitudes, Ls, and LTST. The more precise density factors will be included in Mars-GRAM 2005 Release 1.4 and thus improve the results of future sensitivity studies done for large optical depths.

Justh, Hilary L.↗

An Internal Thermal Environment Model of an Aluminized Solid Rocket Motor with Experimental Validation

Due to the severity of the internal solid rocket motor (SRM) environment, very few direct measurements of that environment exist; therefore, the appearance of such data provides a unique opportunity to assess current thermal/fluid modeling capabilities. As part of a previous study of SRM internal insulation performance, the internal thermal environment of a laboratory-scale SRM featuring aluminized propellant was characterized with two types of custom heat-flux calorimeters: one that measured the total heat flux to a graphite slab within the SRM chamber and another that measured the thermal radiation flux. Therefore, in the current study, a thermal/fluid model of this lab-scale SRM was constructed using ANSYS Fluent to predict not only the flow field structure within the SRM and the convective heat transfer to the interior walls, but also the resulting dispersion of alumina droplets and the radiative heat transfer to the interior walls. The dispersion of alumina droplets within the SRM chamber was determined by employing the Lagrangian discrete phase model that was fully coupled to the Eulerian gas-phase flow. The P1-approximation was engaged to model the radiative heat transfer through the SRM chamber where the radiative contributions of the gas phase were ignored and the aggregate radiative properties of the alumina dispersion were computed from the radiative properties of its individual constituent droplets, which were sourced from literature. The convective and radiative heat fluxes computed from the thermal/fluid model were then compared with those measured in the lab-scale SRM test firings and the modeling approach evaluated.

Martin, Heath T.↗

Development, Verification and Use of Gust Modeling in the NASA Computational Fluid Dynamics Code FUN3D

This paper presents the implementation of gust modeling capability in the CFD code FUN3D. The gust capability is verified by computing the response of an airfoil to a sharp edged gust. This result is compared with the theoretical result. The present simulations will be compared with other CFD gust simulations. This paper also serves as a users manual for FUN3D gust analyses using a variety of gust profiles. Finally, the development of an Auto-Regressive Moving-Average (ARMA) reduced order gust model using a gust with a Gaussian profile in the FUN3D code is presented. ARMA simulated results of a sequence of one-minus-cosine gusts is shown to compare well with the same gust profile computed with FUN3D. Proper Orthogonal Decomposition (POD) is combined with the ARMA modeling technique to predict the time varying pressure coefficient increment distribution due to a novel gust profile. The aeroelastic response of a pitch/plunge airfoil to a gust environment is computed with a reduced order model, and compared with a direct simulation of the system in the FUN3D code. The two results are found to agree very well.

Bartels, Robert E.↗

Strategies to Improve the Accuracy of Mars-GRAM Sensitivity Studies at Large Optical Depths

The Mars Global Reference Atmospheric Model (Mars-GRAM) is an engineering-level atmospheric model widely used for diverse mission applications. Mars-GRAM s perturbation modeling capability is commonly used, in a Monte-Carlo mode, to perform high fidelity engineering end-to-end simulations for entry, descent, and landing (EDL). It has been discovered during the Mars Science Laboratory (MSL) site selection process that Mars-GRAM when used for sensitivity studies for MapYear=0 and large optical depth values such as tau=3 is less than realistic. A comparison study between Mars atmospheric density estimates from Mars- GRAM and measurements by Mars Global Surveyor (MGS) has been undertaken for locations of varying latitudes, Ls, and LTST on Mars. The preliminary results from this study have validated the Thermal Emission Spectrometer (TES) limb data. From the surface to 80 km altitude, Mars- GRAM is based on the NASA Ames Mars General Circulation Model (MGCM). MGCM results that were used for Mars-GRAM with MapYear=0 were from a MGCM run with a fixed value of tau=3 for the entire year at all locations. Unrealistic energy absorption by uniform atmospheric dust leads to an unrealistic thermal energy balance on the polar caps. The outcome is an inaccurate cycle of condensation/sublimation of the polar caps and, as a consequence, an inaccurate cycle of total atmospheric mass and global-average surface pressure. Under an assumption of unchanged temperature profile and hydrostatic equilibrium, a given percentage change in surface pressure would produce a corresponding percentage change in density at all altitudes. Consequently, the final result of a change in surface pressure is an imprecise atmospheric density at all altitudes. To solve this pressure-density problem, a density factor value was determined for tau=.3, 1 and 3 that will adjust the input values of MGCM MapYear 0 pressure and density to achieve a better match of Mars-GRAM MapYear=0 with MapYears 1 and 2 MGCM output at comparable dust loading. Currently, these density factors are fixed values for all latitudes and Ls. Results will be presented of the work underway to derive better multipliers by including possible variation with latitude and/or Ls. This is achieved by comparison of Mars-GRAM MapYear=0 output with TES limb data. The addition of these density factors to Mars-GRAM will improve the results of the sensitivity studies done for large optical depths. Answers may also be provided to the issues raised in a recent study by Desai(2008). Desai has shown that the actual landing sites of Mars Pathfinder, the Mars Exploration Rovers and the Phoenix Mars Lander have been further downrange than predicted by models prior to landing. Desai s reconstruction of their entries into the Martian atmosphere showed that the models consistently predicted higher densities than those found upon EDL. The solution of this problem would be important to the Mars Program since future exploration of Mars by landers and rovers will require more accurate landing capabilities, especially for the proposed Mars Sample Return mission.

Justh, Hilary L.↗

Life prediction of aging aircraft wiring systems

The program goal is to develop a computerized life prediction model capable of identifying present aging progress and predicting end of life for aircraft wiring. A summary is given in viewgraph format of progress made on phase 1 objectives, which were to identify critical aircraft wiring problems; relate most common failures identified to the wire mechanism causing the failure; assess wiring requirments, materials, and stress environment for fighter aircraft; and demonstrate the feasibility of a time-temperature-environment model.

Slenski, George↗

Modeling Specular Exchange Between Concentric Cylinders in a Radiative Shielded Furnace

The objective of this research is to develop and validate mathematical models to characterize the thermal performance of a radiative shielded furnace, the University of Alabama in Huntsville (UAH) Isothermal Diffusion Oven. The mathematical models are validated against experimental data obtained from testing the breadboard oven in a terrestrial laboratory environment. It is anticipated that the validation will produce math models capable of predicting the thermal performance of the furnace over a wide range of operating conditions, including those for which no experimental data is available. Of particular interest is the furnace core temperature versus heater power parametric and the transient thermal response of the furnace. Application to a microgravity environment is not considered, although it is conjectured that the removal of any gravity dependent terms from the math models developed for the terrestrial application should yield adequate results in a microgravity environment. The UAH Isothermal Diffusion Oven is designed to provide a thermal environment that is conducive to measuring the diffusion of high temperature liquid metals. In addition to achieving the temperatures required to melt a sample placed within the furnace, reducing or eliminating convective motions within the melt is an important design consideration [1]. Both of these influences are reflected in the design of the furnace. Reducing unwanted heat losses from the furnace is achieved through the use of low conductivity materials and reflective shielding. As evidenced by the highly conductive copper core used to house the sample within the furnace, convective motions can be greatly suppressed by providing an essentially uniform thermal environment. An oven of this design could ultimately be utilized in a microgravity environment, presumably as a experiment payload. Such an application precipitates other design requirements that limit the resources available to the furnace such as power, mass, volume, and possibly even time. Through the experimental and numerical results obtained, the power requirements and thermal response time of the breadboard furnace are quantified.

Schunk, Richard Gregory↗

New Developments In NASA’s Entry Systems Modeling Project

This paper describes recent developments for modeling entry, descent, and landing (EDL)of spacecraft in support of NASA’s exploration missions. Mission-specific research and model development for entry systems occurs across the Agency (e.g., within flight programs like Artemis/Orion and Mars Sample Return), however the aim of this paper is to discuss the research conducted by NASA’s Entry Systems Modeling (ESM) Project, which serves as the Agency’s only effort dedicated to advancing modeling capabilities that cross-cut multiple technical disciplines, missions and destinations. The ESM portfolio is developed to address the specific needs expressed by a cross-section of NASA stakeholders, including flight and research projects, technical leadership, and subject matter experts. Technology development in ESM is organized and prioritized from a system-level perspective, resulting in four broad technical areas of investment: (1) Thermal Protection System (TPS) material modeling, (2) Shock layer kinetics and radiation, (3) Aerosciences, and (4) Guidance, navigation, and control. In addition to the core technical areas, special topics are rolled into the project portfolio as specific demands arise. Current special topics include TPS Certification by Analysis, improving understanding of woven TPS material performance; Hypersonic Wake Flows, assessing and improving predictions of base flows; and the MEDLI2 Deep Dive, furthering analysis of data obtained during Mars2020 mission’s entry and descent at Mars. Key results from each of these areas are presented in this paper, along with associated references to serve as a roadmap for other EDL researchers to access NASA’s publications.

Aaron M Brandis↗

Transition mixing study

A computer model capable of analyzing the flow field in the transition liner of small gas turbine engines is developed. A FORTRAN code has been assembled from existing codes and physical submodels and used to predict the flow in several test geometries which contain characteristics similar to transition liners, and for which experimental data was available. Comparisons between the predictions and measurements indicate that the code produces qualitative results but that the turbulence models, both K-E and algebraic Reynolds Stress, underestimate the cross-stream diffusion. The code has also been used to perform a numerical experiment to examine the effect of a variety of parameters on the mixing process in transition liners. Comparisons illustrate that geometries with significant curvature show a drift of the jet trajectory toward the convex wall and weaker wake region vortices and decreased penetration for jets located on the convex wall of the liner, when compared to jets located on concave walls. Also shown were the approximate equivalency of angled slots and round holes and a technique by which jet mixing correlations developed for rectangular channels can be used for can geometries.

Reynolds, R.↗