Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “modeling and performance analysis”

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 217 records · Page 12

Validation of MFUEL Metal Fuel Performance Models of SAS4A/SASSYS-1

The fuel characterization models of SAS4A/SASSYS-1 (SAS) have recently been extended to include a new U-Pu-Zr metal fuel model, MFUEL. MFUEL is equipped with mechanistic physics based models to predict the pre-transient characterization and transient response of metal fuel, with emphasis on fuel melting, cladding failure, and the metal fuel’s impact on core reactivity. The MFUEL model will be available in the full version of SAS4A/SASSYS-1 5.7, which is scheduled to be released in June 2023. Fast reactor fuel pins that operated in EBR-II and FFTF with low smear density U-Zr and U-Pu-Zr metal fuels and irradiation resistant ferritic-martensitic cladding showed significant advantages in achieving high burnups and assuring inherent safety characteristics during anticipated transients, design basis events, and beyond design basis events. For safety analysis, a fuel performance model must be able to predict (1) Fuel pin mechanics and compositional and dimensional changes, (2) Clad failure, and (3) Fuel pin thermal resistance. Achieving these high level goals accurately is strongly related to the model performance of individual physical processes taking place within a fuel pin during its lifetime. Metal fuels typically operate above the mid-point temperature of melting during normal, as well as off-normal, conditions. At these elevated temperatures, the availability of thermal activation provides a driving force for various diffusional processes leading to complex phase transformations, micro-structure evolution, significant amounts of fuel swelling, interconnected porosity formation, excessive amounts of fission gas release, and fuel clad chemical interactions. Clad failure in fast reactors primarily occurs as a result of creep rupture augmented by clad wastage formation. The reaction is driven by thermal creep induced dislocation motion, grain boundary cavity nucleation, growth and breakup of grain boundaries. The high level complexity and limited available data requires introducing physics-based modeling approaches to gain extrapolation ability and sensitivity with respect to various conditions. The objective of this report is to perform validation of MFUEL using the experimental data for (1) Normal operation EBR-II fuel behavior, (2) Normal operation PHENIX fuel behavior, (3) HT9 Pressure tube ramp-and-hold creep rupture tests, (4) Whole Pin Furnace (WPF) creep strain, creep rupture and eutectic tests, (5) Fuel Behavior Test Apparatus (FBTA) eutectic tests, and (6) TREAT M5-7 OverPower tests up to clad failure. Section-2 includes a brief description of the MFUEL models. A detailed description of the MFUEL physics-based, semi-empirical models will be presented in the SAS V 5.7 theory manual. Section-3, 4, and 5 describes the validation effort for the pre-transient irradiation, furnace transients, and TREAT M-Series transients, respectively.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Modeling Joule Heating Effect on Lunar O2 Generation via Electrolytic Reduction.

Kennedy Space Center is leading research work on lunar O2 generation via electrolytic reduction of regolith; the metal oxide present in the regolith is dissociated in oxygen anions and metal cations leading to the generation of gaseous oxygen at the anode and liquid metal at the cathode. Electrical resistance of molten regolith is high, leading to heating of the melt when electrical current is applied between the electrodes (Joule heating). The authors have developed a 3D model using a rigorous approach for two coupled physics (thermal and electrical potential) to not only study the effect of Joule heating on temperature distribution throughout the molten regolith but also to evaluate and optimize the design of the electrolytic cells. This paper presents the results of the thermal analysis performed on the model and used to validate the design of the electrolytic cell.

Dominquez, Jesus↗

Methodology for Examining the Operator and the System Concurrently: Pilot Interaction with Automation

Complex system description is problematic when considering operator task activities interacting with system dynamics. Engineering languages have matured sufficiently to allow machine system description at various levels of depth and breadth but without operator synergy. Concurrently, Task Analysis methods have evolved along diverse lines enabling a description of the operator in the system from various paradigms but not describing the system. A void exists when attempting to view the system and the operator in the same plane. We propose a methodology employing descriptive languages from different domains viewed in a single dimension. Finite Automata (FA) languages describe the machine system in the proposed approach. Operator task specifications, a form of task analysis output, examine the operator activities within the system. Operator task specifications were then selected for discrete task activities and overlaid on the system description to examine operator inputs and subsequently view system responses. Unexpected (surprise) and undesirable system behavior was expected to emerge from this analysis. In this paper we shall first describe the methodology and show how the two perspectives, machine model and operator task specifications are integrated. Following we describe the process of doing such analysis using an example from cockpit automation. The methodology was employed in the analysis of a new function that was added to an existing automatic flight control system. We begin by defining a flight scenario involving all aspects of pilot interaction with the new function. Then, we proceed to develop a basic model of the machine behavior, in the context of pilot actions. Finally we superimpose the operator task specification on the machine model and perform the analysis. The proposed methodology may have broad appeal to system designers and human factors specialists. A common language for engineers of diverse domains is a strong point of this approach. Systems engineers may not fully understand operational considerations and human limitations, and human factors specialists may not be exposed to the full extent of system behavior. We hope the proposed methodology is adopted by both disciplines and in the process each assimilates a common language to address potential shortcomings of either approach separately.

Austin, David↗

Enhanced Contact Graph Routing (ECGR) MACHETE Simulation Model

Contact Graph Routing (CGR) for Delay/Disruption Tolerant Networking (DTN) space-based networks makes use of the predictable nature of node contacts to make real-time routing decisions given unpredictable traffic patterns. The contact graph will have been disseminated to all nodes before the start of route computation. CGR was designed for space-based networking environments where future contact plans are known or are independently computable (e.g., using known orbital dynamics). For each data item (known as a bundle in DTN), a node independently performs route selection by examining possible paths to the destination. Route computation could conceivably run thousands of times a second, so computational load is important. This work refers to the simulation software model of Enhanced Contact Graph Routing (ECGR) for DTN Bundle Protocol in JPL's MACHETE simulation tool. The simulation model was used for performance analysis of CGR and led to several performance enhancements. The simulation model was used to demonstrate the improvements of ECGR over CGR as well as other routing methods in space network scenarios. ECGR moved to using earliest arrival time because it is a global monotonically increasing metric that guarantees the safety properties needed for the solution's correctness since route re-computation occurs at each node to accommodate unpredicted changes (e.g., traffic pattern, link quality). Furthermore, using earliest arrival time enabled the use of the standard Dijkstra algorithm for path selection. The Dijkstra algorithm for path selection has a well-known inexpensive computational cost. These enhancements have been integrated into the open source CGR implementation. The ECGR model is also useful for route metric experimentation and comparisons with other DTN routing protocols particularly when combined with MACHETE's space networking models and Delay Tolerant Link State Routing (DTLSR) model.

Segui, John S.↗

Fracture mechanics analysis for various fiber/matrix interface loadings

Fiber/matrix (F/M) cracking was analyzed to provide better understanding and guidance in developing F/M interface fracture toughness tests. Two configurations, corresponding to F/M cracking at a broken fiber and at the free edge, were investigated. The effects of mechanical loading, thermal cooldown, and friction were investigated. Each configuration was analyzed for two loadings: longitudinal and normal to the fiber. A nonlinear finite element analysis was performed to model friction and slip at the F/M interface. A new procedure for fitting a square-root singularity to calculated stresses was developed to determine stress intensity factors (K sub I and K sub II) for a bimaterial interface crack. For the case of F/M cracking at a broken fiber with longitudinal loading, crack tip conditions were strongly influenced by interface friction. As a result, an F/M interface toughness test based on this case was not recommended because nonlinear data analysis methods would be required. For the free edge crack configuration, both mechanical and thermal loading caused crack opening, thereby avoiding frictional effects. A F/M interface toughness test based on this configuration would provide data for K(sub I)/K(sub II) ratios of about 0.7 and 1.6 for fiber and radial normal loading, respectively. However, thermal effects must be accounted for in the data analysis.

Naik, R. A.↗

Fracture mechanics analysis for various fiber/matrix interface loadings

Fiber/matrix (F/M) cracking was analyzed to provide better understanding and guidance in developing F/M interface fracture toughness tests. Two configurations, corresponding to F/M cracking at a broken fiber and at the free edge, were investigated. The effects of mechanical loading, thermal cooldown, and friction were investigated. Each configuration was analyzed for two loadings: longitudinal and normal to the fiber. A nonlinear finite element analysis was performed to model friction and slip at the F/M interface. A new procedure for fitting a square-root singularity to calculated stresses was developed to determine stress intensity factors (K sub I and K sub II) for a bimaterial interface crack. For the case of F/M cracking at a broken fiber with longitudinal loading, crack tip conditions were strongly influenced by interface friction. As a result, an F/M interface toughness test based on this case was not recommended because nonlinear data analysis methods would be required. For the free edge crack configuration, both mechanical and thermal loading caused crack opening, thereby avoiding frictional effects. An F/M interface toughness test based on this configuration would provide data for K(sub I)/K(sub II) ratios of about 0.7 and 1.6 for fiber and radial normal loading, respectively. However, thermal effects must be accounted for in the data analysis.

Naik, R. A.↗

Fracture mechanics analysis for various fiber/matrix interface loadings

Fiber/matrix (F/M) cracking was analyzed to provide better understanding and guidance in developing F/M interface fracture toughness tests. Two configurations, corresponding to F/M cracking at a broken fiber and at the free edge, were investigated. The effects of mechanical loading, thermal cooldown, and friction were investigated. Each configuration was analyzed for two loadings: longitudinal and normal to the fiber. A nonlinear finite element analysis was performed to model friction and slip at the F/M interface. A new procedure for fitting a square-root singularity to calculated stresses was developed to determine stress intensity factors (K sub I and K sub II) for a bimaterial interface crack. For the case of F/M cracking at a broken fiber with longitudinal loading, crack tip conditions were strongly influenced by interface friction. As a result, an F/M interface toughness test based on this case was not recommended because nonlinear data analysis methods would be required. For the free edge crack configuration, both mechanical and thermal loading caused crack opening, theory avoiding fractional effects. A F/M interface toughness test based on this configuration would provide data for K(sub I/K(sub II) ratios of about 0.7 and 1.6 for fiber and radial normal loading, respectively. However, thermal effects must be accounted for in the data analysis.

Naik, Rajiv A.↗

International Space Station Model Correlation Analysis

This paper summarizes the on-orbit structural dynamic data and the related modal analysis, model validation and correlation performed for the International Space Station (ISS) configuration ISS Stage ULF7, 2015 Dedicated Thruster Firing (DTF). The objective of this analysis is to validate and correlate the analytical models used to calculate the ISS internal dynamic loads and compare the 2015 DTF with previous tests. During the ISS configurations under consideration, on-orbit dynamic measurements were collected using the three main ISS instrumentation systems; Internal Wireless Instrumentation System (IWIS), External Wireless Instrumentation System (EWIS) and the Structural Dynamic Measurement System (SDMS). The measurements were recorded during several nominal on-orbit DTF tests on August 18, 2015. Experimental modal analyses were performed on the measured data to extract modal parameters including frequency, damping, and mode shape information. Correlation and comparisons between test and analytical frequencies and mode shapes were performed to assess the accuracy of the analytical models for the configurations under consideration. These mode shapes were also compared to earlier tests. Based on the frequency comparisons, the accuracy of the mathematical models is assessed and model refinement recommendations are given. In particular, results of the first fundamental mode will be discussed, nonlinear results will be shown, and accelerometer placement will be assessed.

Modal Analysis↗

Proportional and scale change models to project failures of mechanical components with applications to space station

In this paper we develop the mathematical theory of proportional and scale change models to perform reliability analysis. The results obtained will be applied for the Reaction Control System (RCS) thruster valves on an orbiter. With the advent of extended EVA's associated with PROX OPS (ISSA & MIR), and docking, the loss of a thruster valve now takes on an expanded safety significance. Previous studies assume a homogeneous population of components with each component having the same failure rate. However, as various components experience different stresses and are exposed to different environments, their failure rates change with time. In this paper we model the reliability of a thruster valves by treating these valves as a censored repairable system. The model for each valve will take the form of a nonhomogeneous process with the intensity function that is either treated as a proportional hazard model, or a scale change random effects hazard model. Each component has an associated z, an independent realization of the random variable Z from a distribution G(z). This unobserved quantity z can be used to describe heterogeneity systematically. For various models methods for estimating the model parameters using censored data will be developed. Available field data (from previously flown flights) is from non-renewable systems. The estimated failure rate using such data will need to be modified for renewable systems such as thruster valve.

Taneja, Vidya S.↗

PERTS: A Prototyping Environment for Real-Time Systems

We discuss an ongoing project to build a Prototyping Environment for Real-Time Systems, called PERTS. PERTS is a unique prototyping environment in that it has (1) tools and performance models for the analysis and evaluation of real-time prototype systems, (2) building blocks for flexible real-time programs and the support system software, (3) basic building blocks of distributed and intelligent real time applications, and (4) an execution environment. PERTS will make the recent and future theoretical advances in real-time system design and engineering readily usable to practitioners. In particular, it will provide an environment for the use and evaluation of new design approaches, for experimentation with alternative system building blocks and for the analysis and performance profiling of prototype real-time systems.

Liu, Jane W. S.↗

Modeling interconnections of safety and financial performance of nuclear power plants, part 3: Spatiotemporal probabilistic physics-of-failure analysis and its connection to safety and financial performance

Here, this paper is a byproduct of a line of research by the authors to analyze interrelationships of safety and financial performance of nuclear power plants (NPPs). The result of this line of research is summarized in three parts: Part 1 covers a categorical review of relevant literature and the theoretical bases that support the methodological developments in Part 2. Part 2 introduces an Integrated Enterprise Risk Management (I-ERM) methodological framework to quantify the interconnections of safety and financial performance with a focus on operation and maintenance (O&M) of NPPs. Part 2 has also demonstrated the applicability and values of the I-ERM methodology through an NPP case study. This paper is Part 3, where detailed development and implementation of one of the I-ERM modules, i.e., probabilistic physics-of-failure (PPoF) analysis, and its connection with safety and financial performance is reported. In this article, the physical failure modeling for hardware components is advanced by incorporating finite element analysis (FEA) into PPoF analysis and coupling the FEA-based PPoF with the maintenance performance through a renewal process model. This article covers two scientific contributions: (i) first-of-its-kind incorporation of FEA into the PPoF model of thermal fatigue for NPP components; and (ii) advancing the interface between the PPoF analysis and the renewal process model in order to deal with spatiotemporal FEA outputs and to efficiently estimate the physical transition rates even when the PPoF outputs are dominated by success data. Through the incorporation of FEA, the resolution of the PPoF analysis is enhanced as spatiotemporal conditions such as stress and temperature can be considered explicitly instead of relying on simplified assumptions or analytical models with reduced spatiotemporal dimensions. To demonstrate an application of the FEA-based PPoF analysis and its coupling with maintenance through the renewal process model, a case study is conducted using excess letdown elbow piping in the chemical and volume control system of a Pressurized Water Reactor.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Visual HPC Workflows for the Analysis of System Dynamics Models

Visual analytics supported by high performance computing (HPC) accelerates and enhances the discovery, exploration, and analysis of causal patterns in complex system dynamics (SD) models. We present a suite of visualization-assisted ensemble-based techniques for hypothesis generation and testing, and for sensitivity analysis. By employing HPC to provide parallel, on-demand simulation of SD models, one can “steer” an ensemble of simulated scenarios in real time as one first formulates and then informally tests those hypotheses: this provides rapid feedback for analysts to refine their understanding of the causal relationships emergent from a model. Such understandings can be followed and augmented by rigorous application of statistical methods, namely global variance-based sensitivity analysis, Monte-Carlo filtering, adaptive regional sensitivity analysis, and self-organized maps: here timely computation relies on HPC, while effective presentation emphasizes high-dimensional multivariate data visualization. Immersive visualization in virtual 3D environments provides an excellent adjunct to the traditional 2D graphics typically used for SD models, as it generates an embodied understanding of model behavior and facilitates an active, collaborative critique of model structure and output. Finally, we summarize prospects for HPC-enabled visual analytics applied to SD modeling.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Techno-Economic Analysis of Synthetic Fuels Pathways Integrated with Light Water Reactors

The purpose of this work is to identify, model, perform techno-economic analysis, and compare two possible synfuel production routes utilizing CO 2 as the feedstock. Heat from an LWR nuclear plant is integrated wherever possible to positively affect the economics of the LWR by converting power to fuels during times of low grid electricity demand. Process and economic modeling for a conceptual synfuel production plant co-located (or in near proximity) with an LWR is presented, including the cost of CO 2 captured from an ethanol plant, compressed, and transported to the LWR hybrid plant, co-electrolysis of the CO 2 with water in a solid oxide electrolyzing cell (SOEC) system to produce syngas, and thermocatalytic conversion of the syngas to transportation fuel. The hybrid LWR/synfuels plant is assumed to be located within 50–150 miles of an ethanol plant (e.g., located in the midwest region of the United States). Performance and nth-plant economics for the co-electrolysis-based processes are evaluated and compared with biomass-gasification-based technology for the synfuel routes considered. Sensitivity analysis around the price of CO 2 and electricity, two of the major cost drivers, is presented for each case. Consideration of a carbon credit is also included in the sensitivity analysis.

10 SYNTHETIC FUELS↗

Uncertainty Quantification Framework for Predicting Material Response with Large Number of Parameters: Application to Creep Prediction in Ferritic-Martensitic Steels Using Combined Crystal Plasticity and Grain Boundary Models

This paper presents an uncertainty quantification (UQ) framework for the physics-based model prediction of material response with a large number of parameters. The application problem presented in this work is that of predicting creep in Grade 91 steel at 600°C. The material response is defined with a physically based microstructural model with constitutive equations emulating several observed phenomena in Grade 91 and embodied into an explicit geometry mesoscale finite element model for prior austenite grains and grain boundaries. Creep within the grains and in grain boundaries are represented by crystal plasticity for dislocation motion and a physics-based model for cavity growth and nucleation, respectively. The creep behavior of this material is influenced by several parameters, some of which have a wide range of variation based on experimental data. UQ combined with microstructural modeling can discover the core microstructural causes of experimental variability, leading to improved materials with lower variability in critical long-term material properties. In this study, we investigate the model's uncertainty to identify material properties that may be modified during production to increase creep life and analyze different components of the crystal plasticity model for improvements. For this purpose, a quantity of interest is defined as time to minimum creep rate, which correlates well to the creep failure of the material. A deep neural network model was trained and validated to be used as a surrogate for the finite element model. Then, a variance-based sensitivity analysis is performed on the surrogate model to find the Sobol indices of the input parameters in respect to the output quantity of interest. The Sobol indices are used to reduce the dimensionality of the model. Generalized polynomial chaos expansion is used on the reduced basis models to propagate the uncertainty from the input parameters to the quantity of interest using the deep neural network surrogate model. These results are benchmarked against uncertainty propagation using Monte Carlo simulations. In conclusion, the UQ performed through the reduced basis model captures almost all the uncertainty in the model with significantly fewer simulations, making it possible to perform the UQ directly via simulations with the finite element model rather than surrogate machine-learned models.

36 MATERIALS SCIENCE↗

IDAES-PSE 2.3.0 Release

The Institute for the Design of Advanced Energy Systems (IDAES) Integrated Platform is a versatile computational environment offering extensive process systems engineering (PSE) capabilities for optimizing the design and operation of complex, interacting technologies and systems. IDAES enables users to efficiently search vast, complex design spaces to discover the lowest cost, most environmentally sustainable solutions while supporting the full process modeling lifecycle, from conceptual design to dynamic optimization and control. The extensible, open platform empowers users to create models of novel processes and rapidly develop custom analyses, workflows, and end-user applications.. New Features and Models • New diagnostics toolboxes and examples o Tutorial for using the DiagnosticsToolbox o Methods to identify potential evaluation errors in models o SVDToolbox for performing singular value decomposition analysis on models to identify possible degeneracies and poor scaling o DegeneracyHunter for identifying irreducible degenerate sets in models • Model for solid-liquid separation which takes two inputs (solid and liquid streams) and produces three outlets (solids, liquid entrained with solids and pure liquid) • Example of temperature swing adsorption in models_extra Deprecation Warnings • With the update to Pyomo 6.7.0, the plate heat exchanger model has begun failing tests on some CI platforms. If this is not addressed by the February release, this model may be removed due to failing tests. See #1294. Offers to help identify the cause and fix this issue are welcome. Improvements and bug fixes • Fixed oversight which precluded using Mixer and Separator models when inherent reactions were present in property package • Added reporting methods to the MSContactor model • Minor corrections in some examples where values were being fixed outside of declared variable bounds

DiagnosticsToolbox,IDAES,PSE,Process Systems Engin↗

Quantitative evaluations of subtropical westerly jet simulations over East Asia based on multiple CMIP5 and CMIP6 GCMs

As a salient feature of the Asian monsoon system, the East Asian subtropical westerly jet (EASWJ) exerts significant impacts on weather and climate changes in China and even throughout East Asia. In this paper, we applied a new self-adaptive algorithm to detect the EASWJ, identify its boundaries, and then represent its characteristics by defining three indices: the intensity index, meridional displacement index, and width index. Compared to the reanalysis data, we carried out a comprehensive, objective, and quantitative EASWJ evaluation using historical experiments from multiple global climate models (GCMs) in the Coupled Model Intercomparison Project Phases 5 and 6 (CMIP5 and CMIP6). The results show that the multimodel ensemble mean (MME) of both CMIP5 and CMIP6 can simulate the characteristics of winter EASWJ well. While for the other three seasons, the MME of both phase models underestimate the 200-hPa zonal wind (U200) strength in the jet coverage area and overestimate the U200 outside the jet area, such simulation results weaken the meridional shear of the wind field. The EASWJ simulations from the CMIP5 GCMs had no consistent intensity or location characteristic tendencies, and most CMIP5 GCMs tended to simulate relatively wide-coverage jets. In contrast, most CMIP6 GCMs are inclined to simulate significantly weaker, wider, and more-northward jets. Compared to the predecessors in CMIP5, about half of CMIP6 GCMs significantly minimized the jet intensity bias, but remarkable errors were still observed in their jet location and coverage representations. Furthermore, the comparative analysis performed by classifying models based on their evaluated simulation results suggested that the simulated performance of the meridional temperature gradient was important for capturing the EASWJ characteristics. Further in-depth study of the causes of model differences is warranted to improve simulation results.

54 ENVIRONMENTAL SCIENCES↗

Photovoltaic performance models - A report card

Models for the analysis of photovoltaic (PV) systems' designs, implementation policies, and economic performance, have proliferated while keeping pace with rapid changes in basic PV technology and extensive empirical data compiled for such systems' performance. Attention is presently given to the results of a comparative assessment of ten well documented and widely used models, which range in complexity from first-order approximations of PV system performance to in-depth, circuit-level characterizations. The comparisons were made on the basis of the performance of their subsystem, as well as system, elements. The models fall into three categories in light of their degree of aggregation into subsystems: (1) simplified models for first-order calculation of system performance, with easily met input requirements but limited capability to address more than a small variety of design considerations; (2) models simulating PV systems in greater detail, encompassing types primarily intended for either concentrator-incorporating or flat plate collector PV systems; and (3) models not specifically designed for PV system performance modeling, but applicable to aspects of electrical system design. Models ignoring subsystem failure or degradation are noted to exclude operating and maintenance characteristics as well.

Smith, J. H.↗

Analysis of Cryogenic Propellant Liquefaction Rates in Cooled Constant-Wall-Temperature Tanks

NASA has been focused on developing technology that would allow the production of cryogenic propellants on the Lunar and Martian surfaces. Utilizing Lunar/Martian resources, the produced gaseous propellants must first be liquefied and stored prior to use on the Moon or Mars ascent vehicle. Liquefaction of cryogenic propellants is a necessary technical development to enable NASA’s future spaceflight goals. This paper presents an overview of a proposed model for a propellant liquefaction system, and the effect of tank wall temperatures on cryogenic propellant liquefaction rates. The propellant system was assumed to be a receiver of a gas from an In-Situ Resource Utilization (ISRU) harvester and condensed the gas by cooling the tank walls. First, the heat transfer principles pertinent to condensation are summarized, followed by an overview of the model which performed the analysis, including a summary of the Fortran algorithm implementation. The sensitivity of the condensation rate relative to varying tank wall temperatures is also discussed.

Anson Koch↗