Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dependability modeling”

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 1,585 records · Page 88

A remapping scheme for history-dependent material state variables

Here, we present an interpolation scheme for use with material models involving history-dependent state variables. The interpolation scheme is based on the solution of a constrained optimization problem, to identify the smoothest monotone curve which produces an admissible solution. Unlike classical schemes used for transport equations, the proposed methodology accounts for the interdependence of the history variables and produces admissible interpolations of the material's state variables. After introducing a general solution scheme for the problem, a simplified and more computationally efficient method is derived. Finally, as a validation problem for the proposed methods, we consider the advection of internal state variables belonging to material models with highly nonlinear yield surfaces.

42 ENGINEERING↗

Gaussian Mixture Model Solvers for the Boltzmann Equation

This report documents our experience constructing a numerical method for the collisional Boltzmann equation that is capable of accurately capturing the collisionless through strongly collisional limits. We explore three different functional representations and present a detailed account of a numerical method based on a spatially dependent Gaussian mixture model (GMM). The Kullback-Leibler divergence is used as a closeness measure and various expectation maximization (EM) solution algorithms are implemented to find a compact representation in velocity space for distribution functions that exhibit significant non-Maxwellian character. We discuss issues that appear with this representation over a range of Knudsen numbers for a prototypical test problem and demonstrate that the strongly collisional limit recovers a solution to Euler's equations. Looking forward, this approach is broadly applicable to the non-relativistic and relativistic collisional Vlasov equations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Sensitivity Study of Multiscale and Phenomenological Elasto-Viscoplastic Grade 91 Material Models for Component-Scale Response

Many advanced nuclear reactor concepts currently being developed are targeting higher operating temperatures relative to the current fleet of light water nuclear reactors, for efficiency gains and other operational considerations. The design of high temperature structural components with reliable long-term operational performance will depend on material models that accurately capture the inelastic deformation mechanisms active in these environments. In this work, we perform a detailed parameter sensitivity analysis of two unified elasto-viscoplastic Grade 91 material models capable of capturing long term high temperature creep deformation. The first model is a phenomelogical material model from the Nuclear Engineering Material Library (NEML) developed at Argonne National Lab. The NEML model parameters and their uncertainty were fit to a range of Grade 91 experimental data using Bayesian Markov Chain Monte Carlo analysis. The second model is a LAROMance data-driven surrogate material model developed at Los Alamos National Lab. The LAROMance model is fit to a large database of responses produced by a mechanistic crystal plasticity based polycrystal model. Parameters for the LAROMance surrogate material model reflect the pedigree of the Grade 91 microstructure. Both material models have been integrated into the Grizzly code, based on the open-source MOOSE multiphysics simulation framework, to simulate both the progression of aging mechanisms and the effects of that aging on nuclear power plant structures. Grizzly is used analyze a three-dimensional Grade 91 piping system to compare the long-term inelastic response predicted by these two fundamentally different models and assess the sensitivity of the material model input parameters on this quantity of interest.

42 ENGINEERING↗

Advancing Multi-Hazard Risk and Safety Considerations for Aging Nuclear Facilities (Revision 1)

This project will demonstrate a multi-hazard time-dependent probabilistic risk assessment (PRA) approach for nuclear facilities considering aging-related deterioration of structures. A generic pressurized water reactor (PWR) reactor subjected to seismic mainshock-aftershock sequences considering the aging of the containment structure will be used as a case study to demonstrate the multi-hazard PRA approach. Using advanced modeling and simulation, seismic mainshock-aftershock fragility functions will be simulated for the containment structure considering aging effects. A multi-hazard PRA model for a generic PWR reactor will be built to quantify the multi-hazard core damage frequency (CDF) and large early release frequency (LERF) with explicit time-dependent modeling of event sequences. To date, the cascading impacts of multi-hazards are not adequately accounted for in the PRA models for nuclear facilities. In addition, for both initial and periodic evaluation of facilities to withstand natural phenomena hazards (NPH), deterioration of the SSCs due to aging and other effects is not adequately considered. By advancing multi-hazard considerations accounting for aging effects, this project will contribute to improved understanding of the safety of aging nuclear facilities. Project outcomes such as the multi-hazard CDF and LERF will also allow facility owners to optimize upgrade and retrofit protocols. This project is divided into two components: (1) advanced modeling and simulations; and (2) multi-hazard PRA. For the first component, Idaho National Laboratory’s (INL) Multi-hazard Analysis for STOchastic time-DOmaiN phenomena (MASTODON) and BlackBear codes will be used to simulate the seismic response and damage of a containment structure under cascading mainshock-aftershock sequences with aging effects. Uncertainties related to the seismic inputs, material parameters, and environmental factors such as temperature and humidity will be identified and propagated to the fragility functions. These fragility functions, which consider aging effects, will be time dependent. The second component will take the fragility information from the first component to build a multi-hazard time-dependent PRA model starting from a generic PWR model to evaluate the multi-hazard CDF and LERF and characterize the associated consequences. For this component, OpenPRA Web Application and SAPHIRE code will be used. Events such as the Fukushima Daiichi accident have highlighted the importance of considering the cascading impacts of multi-hazards for PRA. Moreover, many of the reactors in the current nuclear fleet in the United States (US) are already operating well beyond their initially planned design life, and applications for further extensions to operating licenses are being considered. Therefore, considering multi-hazard effects and aging deterioration in the NPH risk assessment process will contribute to the safety of both existing and future nuclear facilities. Outcomes of this project will thus directly benefit the standards DOE-STD-1020-16 and DOE-HDBK-1224-18.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Analytical Modeling of the High Strain Rate Deformation of Polymer Matrix Composites

The results presented here are part of an ongoing research program to develop strain rate dependent deformation and failure models for the analysis of polymer matrix composites subject to high strain rate impact loads. State variable constitutive equations originally developed for metals have been modified in order to model the nonlinear, strain rate dependent deformation of polymeric matrix materials. To account for the effects of hydrostatic stresses, which are significant in polymers, the classical 5 plasticity theory definitions of effective stress and effective plastic strain are modified by applying variations of the Drucker-Prager yield criterion. To verify the revised formulation, the shear and tensile deformation of a representative toughened epoxy is analyzed across a wide range of strain rates (from quasi-static to high strain rates) and the results are compared to experimentally obtained values. For the analyzed polymers, both the tensile and shear stress-strain curves computed using the analytical model correlate well with values obtained through experimental tests. The polymer constitutive equations are implemented within a strength of materials based micromechanics method to predict the nonlinear, strain rate dependent deformation of polymer matrix composites. In the micromechanics, the unit cell is divided up into a number of independently analyzed slices, and laminate theory is then applied to obtain the effective deformation of the unit cell. The composite mechanics are verified by analyzing the deformation of a representative polymer matrix composite (composed using the representative polymer analyzed for the correlation of the polymer constitutive equations) for several fiber orientation angles across a variety of strain rates. The computed values compare favorably to experimentally obtained results.

Goldberg, Robert K.↗

OCHRE: The Object-oriented, Controllable, High-resolution Residential Energy Model for Dynamic Integration Studies

Electrification and the growth of distributed energy resources (DERs), including flexible loads, are changing the energy landscape of electric distribution systems and creating new challenges and opportunities for electric utilities. Changes in demand profiles require improvements in distribution system load models, which have not historically accounted for device controllability or impacts on customer comfort. Although building modeling research has focused on these features, there is a need to incorporate them into distribution load models that include DERs and can be used to study grid-interactive buildings. In this paper, we present the Object-oriented, Controllable, High-resolution Residential Energy (OCHRE) model. OCHRE is a controllable thermal-electric residential energy model that captures building thermal dynamics, integrates grid-dependent electrical behavior, contains models for common DERs and end-use loads, and simulates at a time resolution down to 1 minute. It includes models for space heaters, air conditioners, water heaters, electric vehicles, photovoltaics, and batteries that are externally controllable and integrated in a co-simulation framework. Using a proposed zero energy ready community in Colorado, we co-simulate a distribution grid and 498 all-electric homes with a diverse set of efficiency levels and equipment properties. We show that controllable devices can reduce peak demand within a neighborhood by up to 73% during a critical peak period without sacrificing occupant comfort. We also demonstrate the importance of modeling load diversity at a high time resolution when quantifying power and voltage fluctuations across a distribution system.

24 POWER TRANSMISSION AND DISTRIBUTION↗

A physics-informed and hierarchically regularized data-driven model for predicting fluid flow through porous media

This paper presents a new deep learning data-driven model for predicting structure dependent pore-fluid velocity fields in rock. The model is based on a Convolutional Auto-Encoder (CAE) artificial neural network capable of learning from image data generated by direct numerical simulations of fluid flow through pore-structures, such as by Lattice Boltzmann or molecular dynamics methods. The main novelty of the model in comparison to previous CAE-based data-driven approaches consists of three parts. The first is a methodology for decomposing the full-domain of the porous media into sub-regions, or “sub-domains”, in order to reduce the overall size of the CAE, batch process the sub-domains in parallel, and enable the CAE to learn local and generalizable nonlinear mappings of pore-fluid velocities. The second consists of embedding the finite difference solutions of the incompressible Navier-Stokes and continuity equations into convolutional layers prior to the CAE in order to provide the CAE with knowledge of fluid dynamics physics (PhyFlow). The third main novelty is that the training of the CAE is regularized with a hierarchical loss function that encourages the learning of fluid flow patterns (in a way similar to ranked modes in principal component analysis), ranking from most to least important. This is shown to increase the stability in learning, reduce over-fitting, and promote interpretability of the CAE neural network layers (HierCAE). The comprehensive new data-driven model, which we call the PhyFlow-HierCAE model, is shown to exhibit improved accuracy and generalizability of flow field predictions over conventional CAE models, attributable to the embedded physical knowledge and the hierarchical regularization, as well as realize orders of magnitude speed-ups in computation times as a surrogate for the direct numerical simulations. Examples of training and forward predictions on unseen pore-structures are provided and evaluated for data from Lattice Boltzmann and molecular dynamics simulations of pore-fluid flow. The model is shown to be a fast and accurate emulator (or “surrogate”) for predicting effective permeability of unseen pore-structures based on learning from relatively small direct numerical simulation datasets.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

A Dynamic Risk Framework for the Physical Security of Nuclear Power Plants

This paper describes ongoing work within the Light Water Reactor Sustainability pathway at Idaho National Laboratory (INL) to optimize the security and cost of nuclear power plants. It introduces the dynamic risk assessment tool developed at INL, Event Modeling Risk Assessment using Linked Diagrams (EMRALD). EMRALD is leveraged to optimize the security posture of a nuclear power plant by integrating force-on-force (FOF) simulations and operator mitigation actions, including dynamic and flexible coping strategies (FLEX). To illustrate the methodology, four attack scenarios are modeled in a commercially available FOF simulation tool using a hypothetical nuclear power plant facility. The simulation results provide valuable insights into possible attack outcomes, as well as the probabilistic risk of a core damage event given these outcomes. Safety mitigation procedures are modeled in EMRALD dependent on the attack outcomes by considering human operator uncertainties. The results demonstrate that the number of armed responders can be optimized, while still maintaining the same protection level as the initial security posture. The proposed modeling and simulation framework of integrating FLEX equipment with FOF models enables the nuclear power plants to credit FLEX portable equipment in the plant security posture, resulting in an efficient and optimized physical security system.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Assessing the impact of valence sd neutrons and protons on fusion

Experimental near-barrier fusion cross sections for 17 F + 12 C are compared to the fusion excitation functions for 16,17,18 O, 19 F, and 20 Ne ions on a carbon target. Normalized or reduced fusion cross sections are utilized in order to compare the different systems and account for the differing static size of the incident ions as well as changes in fusion barrier. Fusion excitation functions for the case of the mirror nuclei 17 F and 17 O with a single sd nucleon are compared. The 17 F data are also juxtaposed with nuclei involving multiple sd nucleons. Trends of the fusion cross section above the barrier beyond the expected systematic behavior are observed. These trends are interpreted as the interplay of the sd protons and neutrons. Furthermore, the experimental data are also compared to a widely used analytic model of near-barrier fusion, a time-dependent Hartree-Fock model, and coupled channels calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Performance analysis of wideband data and television channels

Several aspects are discussed of space shuttle communications, including the return link (shuttle-to-ground) relayed through a satellite repeater (TDRS). The repeater exhibits nonlinear amplification and an amplitude-dependent phase shift. Models were developed for various link configurations, and computer simulation programs based on these models are described. Certain analytical results on system performance were also obtained. For the system parameters assumed, the results indicate approximately 1 db degradation relative to a link employing a linear repeater. While this degradation is dependent upon the repeater, filter bandwidths, and modulation parameters used, the programs can accommodate changes to any of these quantities. Thus the programs can be applied to determine the performance with any given set of parameters, or used as an aid in link design.

Geist, J. M.↗

The SYSGEN user package

The user documentation of the SYSGEN model and its links with other simulations is described. The SYSGEN is a production costing and reliability model of electric utility systems. Hydroelectric, storage, and time dependent generating units are modeled in addition to conventional generating plants. Input variables, modeling options, output variables, and reports formats are explained. SYSGEN also can be run interactively by using a program called FEPS (Front End Program for SYSGEN). A format for SYSGEN input variables which is designed for use with FEPS is presented.

Carlson, C. R.↗

Hot-spot tectonics of Eistla Regio, Venus: Results from Magellan images and Pioneer Venus gravity

Eistla Regio (ER) is a broad, low, discontinuous topographic rise striking roughly EW at low northern latitudes of Venus. Some 2000 x 7000 km in dimensions, it is the third largest rise in planform on Venus after Aphrodite Terra and Beta Phoebe Regiones. These rises are the key physiographic elements in a hot spot model of global tectonics including transient plume behavior. Since ER is the first such rise viewed by Magellan and the latitude is very favorable for Pioneer Venus gravity studies, some of the predictions of a time dependent hot spot model are tested. Western ER is defined as the rise including Gula and Sif Mons and central ER as that including Sappho Patera. Superior conjunction prevented Magellan from returning data on eastern ER (Pavlova) during the first mapping cycle. It is concluded that the western and central portions of ER, while part of the same broad topographic rise and tectonic framework, have distinctly different surface ages and gravity signatures. The western rise, including Gula and Sif Mons, is the expression of deep seated uplift with volcanism limited to the individual large shields. The eastern portion has been widely resurfaced more recently by thermal anomalies in the mantle.

Grimm, Robert E.↗

The galactic distribution of aliphatic hydrocarbons in the diffuse interstellar medium

The infrared absorption feature near 2950(exp -1) (3.4 micron), characteristic of dust in the diffuse interstellar medium (ISM), is attributed to C-H stretching vibrations of aliphatic hydrocarbons. We show here that the strength of the band does not scale linearly with visual extinction everywhere, but instead increases more rapidly for objects near the center of the Galaxy, a behavior that parallels that of the Si-O stretching band due to silicate materials in the diffuse ISM. This implies that the grains responsible for the diffuse medium aliphatic C-H and silicate Si-O stretching bands are different from those responsible for much of the observed visual extinction. It also suggests that the distribution of the carbonaceous component of the diffuse ISM is not uniform throughout the Galaxy, but instead may increase in density toward the center of the Galaxy. The similar behavior of the C-H and Si-O stretching bands suggests that these two components may be coupled, perhaps in the form of silicate-core, organic-mantle grains. Several possible models of the distribution of this material are presented and it is demonstrated that the inner parts of the Galaxy has a carrier density that is 5 to 35 times higher than in the local ISM. Depending on the model used, the density of aliphatic material in the local ISM is found to be about 1 to 2 -CH3 groups m(exp -3) and about 2 to 5 -CH2- groups m(exp -3). These densities are consistent with the strengths of the 2955 and 2925 cm(exp -1) (3.4 micron) band being described by the relations A(sub nu)/tau(sub 2955 cm(exp -1)) = 270 +/- 40 and A(sub nu)/tau(sub 2925 cm(exp -1)) = 250 +/- 40 in the local diffuse ISM.

Sandford, Scott A.↗

The Galactic Distribution of Aliphatic Hydrocarbons in the Diffuse Interstellar Medium

The infrared absorption feature near 2950/ cm (3.4 microns), characteristic of dust in the diffuse interstellar medium (ISM), is attributed to C-H stretching vibrations of aliphatic hydrocarbons. We show here that the strength of the band does not scale linearly with visual extinction everywhere, but instead increases more rapidly for objects near the center of the Galaxy, a behavior that parallels that of the Si-O stretching band due to silicate materials in the diffuse ISM. This implies that the grains responsible for the diffuse medium aliphatic C-H and silicate Si-O stretching bands are different from those responsible for much of the observed visual extinction. It also suggests that the distribution of the carbonaceous component of the diffuse ISM is not uniform throughout the Galaxy, but instead may increase in density toward the center of the Galaxy. The similar behavior of the C-H and Si-O stretching bands suggests that these two components may be coupled, perhaps in the form of silicate-core, organic-mantle grains. Several possible models of the distribution of this material are presented and it is demonstrated that the inner parts of the Galaxy has a carrier density that is 5 to 35 times higher than in the local ISM. Depending on the model used, the density of aliphatic material in the local ISM is found to be about 1 to 2 -CH3 groups /cc and about 2 to 5 -CH2- groups/cc. These densities are consistent with the strengths of the 2955 and 2925/ cm (3.38 and 3.42 microns) subfeatures (due to -CH3 and -CH2- groups, respectively) within the overall 2950/ cm (3.4 microns) band being described by the relations A(sub upsilon)/tau(sub 2925/cm) = 270 +/- 40 and A(sub upsilon)/tau(sub 2925/cm) = 250 +/- 40 in the local diffuse ISM.

Sandford, Scott A.↗

Regression Model Term Selection for the Analysis of Strain-Gage Balance Calibration Data

The paper discusses the selection of regression model terms for the analysis of wind tunnel strain-gage balance calibration data. Different function class combinations are presented that may be used to analyze calibration data using either a non-iterative or an iterative method. The role of the intercept term in a regression model of calibration data is reviewed. In addition, useful algorithms and metrics originating from linear algebra and statistics are recommended that will help an analyst (i) to identify and avoid both linear and near-linear dependencies between regression model terms and (ii) to make sure that the selected regression model of the calibration data uses only statistically significant terms. Three different tests are suggested that may be used to objectively assess the predictive capability of the final regression model of the calibration data. These tests use both the original data points and regression model independent confirmation points. Finally, data from a simplified manual calibration of the Ames MK40 balance is used to illustrate the application of some of the metrics and tests to a realistic calibration data set.

Ulbrich, Norbert Manfred↗

The Impact of All-Sky Hyperspectral Infrared Radiance Assimilation on the Simulation and Forecast of Hurricane Sally in GEOS

Hyperspectral infrared (IR) radiance observations have been one of the major data sources assimilated in the data assimilation system over the last 20 years. However, observations peaking in the lower and mid-troposphere are underutilized in clear-sky radiance data assimilation as the quality control procedure removes a significant portion of cloud-affected observations from the assimilation. To include these cloud-affected observations, one approach is to assimilate them directly under all-sky conditions. The framework for all-sky assimilation of hyperspectral IR radiance observations has been developed in GEOS. The preliminary evaluation of simulated cloud-affected IR observations has been conducted, along with the corresponding sensitivities with respect to all hydrometeors. The symmetric cloud effect has been used as a cloud proxy in the observation error modeling, providing a balanced representation that mitigates the discrepancies between observations and model simulations for IR all-sky assimilation. The observation errors vary for different cloud conditions and have been modeled as cloud amount dependent and inter-channel correlated. Preliminary evaluation of model simulations for IR all-sky assimilation indicates that the simulated brightness temperature tends to have broader structures and lacks small-scale details. Excessive clouds generated by the model are also observed. To address this issue, Hurricane Sally (2020), which formed near the Bahamas, is used as a case study. Several cloud overlap schemes and cloud lookup tables (LUTs), along with other important aspects of all-sky IR radiance assimilation, will be evaluated and tested, aiming to improve the impact on Hurricane Sally's simulation and forecast.

Wei Gu↗

A Dynamic Risk Framework for the Optimization of Physical Security Posture of Nuclear Power Plants

This paper describes an ongoing work within the Light Water Reactor Sustainability pathway at Idaho National Laboratory (INL) to optimize security and cost of nuclear power plants. It introduces the dynamic risk assessment tool developed at INL, Event Modeling Risk Assessment using Linked Diagrams (EMRALD). EMRALD was leveraged to optimize the security posture of a nuclear power plant by integrating force-on-force (FOF) simulations and operator mitigation actions including the dynamic and flexible coping strategies (FLEX). To illustrate the methodology, four attack scenarios were modeled in a commercially available FOF simulation tool using a hypothetical nuclear power plant facility. The simulation results provide valuable insights into possible attack outcomes, as well as the probabilistic risk of core damage event given these outcomes. Safety mitigation procedures were modeled in EMRALD dependent on the attack outcomes by considering human operator uncertainties. The results demonstrate that the number of armed responders can be optimized, while still maintaining the same protection level as the initial security posture. The proposed modeling and simulation framework of integrating FLEX equipment with FOF models enables the nuclear power plants to credit FLEX portable equipment in the plant security posture, resulting in an efficient and optimized physical security system.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗