SEARCH · Search NASA

Results for “model determination”

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.

SourceResearch connectionFrom yearThrough year

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

At least 181 records · Page 10

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Crystal plasticity modeling and analysis for the transition from intergranular to transgranular failure in nickel-based alloy Inconel 740H at elevated temperature

The precipitation-strengthened Nickel alloy Inconel® 740H® (IN740H) exhibits increased ductility at higher applied strain rates during quasi-static tensile tests at an elevated temperature of 760°C. The examination of fracture surfaces in this context reveals a noteworthy transition of underlying fracture mechanisms from transgranular to intergranular fracture as the applied strain rate decreases from 1×10 -3 / s to 0.83×10 -4 / s . To thoroughly understand the mechanical response of IN740H under these conditions, this study develops a crystal plasticity finite element (CPFE) model. Further, this model incorporates various deformation mechanisms including dislocation slips, climb, and grain boundary sliding, which are relevant to the test conditions. The model is calibrated using data from both tensile tests at different strain rates and creep tests across a broad stress range at 760°C, enabling the accurate determination of model parameters for each mechanism. Simulation results well captured the experimental observations of different failure modes. At higher strain rates, the model shows a dominance of dislocation slip leading to heterogeneous plastic deformation and formation of transgranular shear bands causing the failure, while at lower strain rates, an increased activity of grain boundary sliding causes grain boundaries crack leading to intergranular failure.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

MODELING AND METHOD FOR OPTICAL PERFORMANCE OF A RECEIVER COATING AND ITS DEGRADATION DUE TO OPERATION

The optical properties of the coating used on a concentrated solar power (CSP)’s receiver directly impact the performance of the power plant. Receiver coatings are designed to absorb as much heat as possible and transfer it to the circulating fluid. Over time, the absorptivity of the receiver coating degrades due to environmental and operational conditions. This project aims to (1) create a test apparatus that allows coatings to undergo accelerated lifetime testing under conditions that mimic what the coating will experience on-tower; and (2) develop a mathematical model that based on the test data predicts absorptivity over time. To create the mathematical model, samples of BrightSource Energy’s Gen1 and Gen2 coatings were tested first looking at individual failure modes (such as high temperature cycling) and then for combined failure modes in the newly developed apparatus. In addition to taking optical measurements, the samples were destroyed to allow for metalogical evaluation to understand what was happening on a microscopic level. Pre-existing models from the literature that describe accelerated testing as well as individual phenomena observed in the coating were then combined to create the model. Additionally, the Gen2 and Gen3 coatings are being testing at the Plataforma Solar de Almería (PSA) facility in Almeria, Spain to confirm that the apparatus test results are similar to those received on-sun. (PSA’s facility closely mimics on-tower conditions). The goal is to enable industry players to test coatings using the new apparatus, gather absorptivity measurements (without needing to destroy samples), and input the information into the model to determine how the coating will perform over time, thus allowing the optical performance of coatings to be compared even before they are used on-tower. The model will be publicly available by the end of the project (currently projected as December 2024).

14 SOLAR ENERGY↗

DOE OSTI2023

Simulating organic aerosol in Delhi with WRF-Chem using the volatility-basis-set approach: exploring model uncertainty with a Gaussian process emulator

The nature and origin of organic aerosol in the atmosphere remain unclear. The gas–particle partitioning of semi-volatile organic compounds (SVOCs) that constitute primary organic aerosols (POAs) and the multigenerational chemical aging of SVOCs are particularly poorly understood. The volatility basis set (VBS) approach, implemented in air quality models such as WRF-Chem (Weather Research and Forecasting model with Chemistry), can be a useful tool to describe emissions of POA and its chemical evolution. However, the evaluation of model uncertainty and the optimal model parameterization may be expensive to probe using only WRF-Chem simulations. Gaussian process emulators, trained on simulations from relatively few WRF-Chem simulations, are capable of reproducing model results and estimating the sources of model uncertainty within a defined range of model parameters. In this study, a WRF-Chem VBS parameterization is proposed; we then generate a perturbed parameter ensemble of 111 model runs, perturbing 10 parameters of the WRF-Chem model relating to organic aerosol emissions and the VBS oxidation reactions. This allowed us to cover the model's uncertainty space and to compare outputs from each run to aerosol mass spectrometer observations of organic aerosol concentrations and O:C ratios measured in New Delhi, India. The simulations spanned the organic aerosol concentrations measured with the aerosol mass spectrometer (AMS). However, they also highlighted potential structural errors in the model that may be related to unsuitable diurnal cycles in the emissions and/or failure to adequately represent the dynamics of the planetary boundary layer. While the structural errors prevented us from clearly identifying an optimized VBS approach in WRF-Chem, we were able to apply the emulator in the following two periods: the full period (1–29 May) and a subperiod period of 14:00–16:00 h LT (local time) on 1–29 May. The combination of emulator analysis and model evaluation metrics allowed us to identify plausible parameter combinations for the analyzed periods. We demonstrate that the methodology presented in this study can be used to determine the model uncertainty and to identify the appropriate parameter combination for the VBS approach and hence to provide valuable information to improve our understanding of OA production.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Bayesian analysis for estimating statistical parameter distributions of elasto-viscoplastic material models

High temperature design methods rely on constitutive models for inelastic deformation and failure typically calibrated against the mean of experimental data without considering the associated scatter. Variability may arise from the experimental data acquisition process, from heat-to-heat material property variations, or both and need to be accurately captured to predict parameter bounds leading to efficient component design. Applying the Bayesian Markov Chain Monte Carlo (MCMC) method to produce statistical models capturing the underlying uncertainty in the experimental data is an area of ongoing research interest. This work varies aspects of the Bayesian MCMC method and explores their effect on the posterior parameter distributions for a uniaxial elasto-viscoplastic damage model using synthetically generated reference data. From our analysis with the uniaxial inelastic model we determine that an informed prior distribution including different types of test conditions results in more accurate posterior parameter distributions. The parameter posterior distributions, however, do not improve when increasing the number of similar experimental data. Additionally, changing the amount of scatter in the data affects the quality of the posterior distributions, especially for the less sensitive model parameters. Moreover, we perform a sensitivity study of the model parameters against the likelihood function prior to the Bayesian analysis. The results of the sensitivity analysis help to determine the reliability of the posterior distributions and reduce the dimensionality of the problem by fixing the insensitive parameters. The comprehensive study described in this work demonstrates how to efficiently apply the Bayesian MCMC methodology to capture parameter uncertainties in high temperature inelastic material models. Quantifying these uncertainties in inelastic models will improve high temperature engineering design practices and lead to safer, more effective component designs.

42 ENGINEERING↗

DOE OSTI2024

A RELAP5-3D Model of the Lobo Lead Loop

The goal of our research is to build upon the capability of RELAP5-3D to model molten lead systems. Molten lead has several potential uses in future advanced reactors like the lead fast reactor or fusion reactors that utilize dual-coolant lead lithium blankets. This potential for use in future generations of reactors highlights the necessity to develop molten lead models to ensure that they can accurately predict the thermohydraulic behavior. We have developed a RELAP5-3D model of the Lobo Lead Loop facility, located at the University of New Mexico, to verify the accuracy of RELAP5-3D via comparison to existing computational fluid dynamics results and analytical calculations. It was found that RELAP5-3D accurately calculated radiative heat transfer (within < 1%) when compared to theoretical calculations. In addition, pressure drop calculations done in RELAP5-3D demonstrated reasonable agreement within 20 kPa, mostly within ~7-15%, when compared to the computational fluid dynamics model of the facility developed by the University of New Mexico, and captured the dependence of pressure drop on flow velocity accurately. Finally, a hypothetical loss of flow transient was imposed on the RELAP5-3D model to determine the feasibility of performing a similar experiment with the Lobo Lead Loop. It was found that such an experiment could be possible as the RELAP5-3D model indicated that the temperatures of the fluid would not exceed limiting temperatures of the structure (1658 K) nor the maximum temperature of the electromagnetic pump inlet (823 K). Although there is not experimental data to begin validation, the model will be readily available for future validation studies when the experimental data is generated, especially as the model continues to evolve over time. Furthermore, the results so far demonstrate a promising first step in the verification/validation of the RELAP5-3D model of the Lobo Lead Loop.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2020

Modeling of the Coral Microbiome: the Influence of Temperature and Microbial Network

Host-associated microbial communities are shaped by extrinsic and intrinsic factors to the holobiont organism. Environmental factors and microbe-microbe interactions act simultaneously on the microbial community structure, making the microbiome dynamics challenging to predict. The coral microbiome is essential to the health of coral reefs and sensitive to environmental changes. Here, we develop a dynamic model to determine the microbial community structure associated with the surface mucus layer (SML) of corals using temperature as an extrinsic factor and microbial network as an intrinsic factor. The model was validated by comparing the predicted relative abundances of microbial taxa to the relative abundances of microbial taxa from the sample data. The SML microbiome from Pseudodiploria strigosa was collected across reef zones in Bermuda, where inner and outer reefs are exposed to distinct thermal profiles. A shotgun metagenomics approach was used to describe the taxonomic composition and the microbial network of the coral SML microbiome. By simulating the annual temperature fluctuations at each reef zone, the model output is statistically identical to the observed data. The model was further applied to six scenarios that combined different profiles of temperature and microbial network to investigate the influence of each of these two factors on the model accuracy. The SML microbiome was best predicted by model scenarios with the temperature profile that was closest to the local thermal environment, regardless of the microbial network profile. Our model shows that the SML microbiome of P. strigosa in Bermuda is primarily structured by seasonal fluctuations in temperature at a reef scale, while the microbial network is a secondary driver. Coral microbiome dysbiosis (i.e., shifts in the microbial community structure or complete loss of microbial symbionts) caused by environmental changes is a key player in the decline of coral health worldwide. Multiple factors in the water column and the surrounding biological community influence the dynamics of the coral microbiome. However, by including only temperature as an external factor, our model proved to be successful in describing the microbial community associated with the surface mucus layer (SML) of the coral P. strigosa. The dynamic model developed and validated in this study is a potential tool to predict the coral microbiome under different temperature conditions.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2021

Managing Solar Photovoltaic Integration in the Western United States Appendix: Reference and High Solar Photovoltaic Scenarios for Three Regions [Slides]

This slide deck is an appendix to a paper series that examines potential challenges related to planning future power systems with higher solar photovoltaic (PV) penetrations. The series uses the western U.S. power system for these investigations because it is a region the authors and their colleagues have already extensively studied. We are therefore well-suited to analyze even higher PV penetrations and then examine the results in multiple models to determine whether our current approaches are missing key details that only emerge at higher PV penetrations. This deck details the systems underlying those analyses and how they were modeled using the Resource Planning Model (RPM), a capacity expansion modeling tool. We examine both Western Interconnection-wide and regional results for three regions in the Western U.S. with significantly different existing power systems and connections to neighboring regions; this provides a more balanced picture as to how power systems with high PV penetration might emerge in different contexts and what the resulting grid challenges, if any, might be.

14 SOLAR ENERGY↗

DOE OSTI2025

DEVELOPMENT AND APPLICATION OF RISK ANALYSIS TOOLKIT FOR PLANT RESOURCE OPTIMIZATION

This paper presents the development of methods and tools that are being designed to optimize plant operations (e.g., maintenance/replacement schedules and optimal maintenance postures for plant components) in a manner that is more cost effective than current approaches and makes better use of available component health and cost data. These methods include both data- and model-based optimization methods. Model-based optimization methods directly include reliability and cost models to determine an optimal plant operational strategy. We consider gradient-based and evolutionary (based on genetic algorithms) optimization methods. The second class of methods target more specific use cases (e.g., project schedule optimization) and are not based on reliability models directly, but they require specific component reliability and cost data. This class of methods is based on variants of the knapsack problem with an aim to determine an optimal project schedule that maximizes the overall NPV. This paper also presents multi-objective methods designed to identify an optimal maintenance posture based on a Pareto frontier analysis. Rather than dictating the “right” tradeoff (i.e., identify the absolute best posture), we show how it is possible to perform a trade space exploration approach (i.e., identify value and costs of several postures and let the analysis account for desired value and cost metrics). This is performed by identifying maintenance postures that maximize value (e.g., system availability) and minimize operational costs, i.e., the Pareto frontier in a value-cost trade space. For all these methods we present detailed applicative examples that show their validity from a decision-making perspective.

97 - MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Comparing Regional Energy Consumption for Direct Drone and Truck Deliveries

Drone delivery, once thought of as fictitious, is becoming a reality with the efforts of both forward-looking enterprises and supportive government policies. This emerging mode of e-commerce delivery raises many concerns. One important concern is the energy efficiency of direct delivery drones compared with conventional delivery trucks at a regional systems level. Here, in this study, we develop and apply methods to quantify the regional energy impacts of drone delivery, then we assess these impacts and compare them with the impacts of truck delivery. To study this problem, we develop an optimization model that determines an optimal set of fulfillment centers (FCs) with variable service capacities that allow drones to make direct e-commerce deliveries. We adopt two drone delivery energy estimation models from the literature and use them as inputs to demonstrate the potential range of energy needs. We also develop another optimization model to account for the energy consumption of diesel trucks (DTs) and battery electric vehicles (BEVs). We test the models using validated simulation data for the Chicago metropolitan area in the U.S. to quantify the energy implications of these three delivery modes. For drone delivery, we further extend our analyses by considering the impact of wind speed and flight patterns. Our results show that direct delivery drones require 15.8% more energy than BEVs on an average windy day, and they need 15% more energy than DTs on a very windy day. We provide essential parameter values for reproducibility and list relevant open problems.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2022

Hanford Waste Treatment Plant Effluent Management Facility Stack Effluent Monitoring: Sampling Probe Location Qualification Evaluation

The Hanford Tank Waste Treatment and Immobilization Plant Effluent Management facility (EMF) stack monitor location was qualified using the LV-S1 scale model stack as a baseline, augmented by the LB-S1 and LV-S2 scale model stacks to address the Direct Feed Low Activity Waste Effluent Management Facility Vessel Vent Process (DVP) injection into the main Active Confinement Ventilation (ACV) system duct. As required by the American National Standards Institute/Health Physics Society (ANSI/HPS) N13.1-1999 standard, the scale models and its sampling locations were geometrically similar to the actual stack and the Reynolds numbers for both the actual and model stacks were >10,000. An additional criterion is that the product of the hydraulic diameter and mean velocity (DV) of the full-scale stack must be between 1/6 DV and 6 DV of the scale model stack tests. Verification tests of the EMF stack was performed at normal operating conditions. The minimum 1/6 DV value, along with the maximum 6 DV value from the scale model testing determines the range of stack flow rates for which the full-scale stack may be operated and remain in compliance with the stack verification criterion. Based on the LV-S1 scale model test DV values, the corresponding stack flow rates for the EMF stack are as listed in Table S1. Table S1. Effluent Management Facility Stack Qualified Flow Range. Stack Parameter EM-1 Minimum Qualified Stack Flow (scfm) 781 Maximum Qualified Stack Flow (scfm) 53,432 The remaining criteria for the stack verification to be considered valid involve the flow angle and velocity uniformity results. First, the flow angle at the full-scale stack must be =20°. Second, the velocity uniformity at the full-scale stack must be =20% coefficient of variance (COV). Finally, the velocity uniformity results for the actual and scale model stack tests must agree within 5% COV. These criteria were met through the full-scale stack test at the EMF. Flow angle results were less than 5°; all flow angle results were within the =20° criterion. The velocity uniformity results for each test condition ranged between 2.2% COV and 4.3% COV, which were all within the range of the target % COV values from the scale model tests on the LV-S1, LB-S1, and LV S2 scale models. Based on these stack verification test results, the EMF filtered exhaust stack sampling location meets the qualification criteria provided in the ANSI/HPS N13.1-1999 standard for all fan operating configurations. This includes each combination of ACV fan with the DVP exhausters. Further changes to the system configuration or operating conditions that are outside the qualified flow rates described in this report may require additional tests or analyses to determine compliance with the standard.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

DOE OSTI2025

Using system identification in modeling the yaw response of tail fins for small wind turbines with bearing friction

Here, we describe three main developments of our previous study of the nonlinear yaw dynamics of the tail fins for small wind turbines [Khedr et al., J. Renewable Sustainable Energy 16, 053305 (2024)]. First, the model constants derived from archived computational and experimental studies are adjusted by employing system identification (SI) to maximize the model's agreement with wind tunnel tests. This adjustment was done for high wind speeds, where yaw bearing friction can be ignored. When starting a turbine at low wind speed, however, friction can become important. Our second development is to implement a model for the frictional resistive torque and use SI to maximize its accuracy. These developments used wind tunnel experiments on generic delta, elliptical, and rectangular planforms that were described by Khedr et al. [J. Renewable Sustainable Energy 16, 053305 (2024)]. Since the aerodynamic and friction models employ a large number of constants, we describe ways to constrain the values using linearized solutions of the response equations for small and large yaw angles. Third, we test the generality of the aerodynamic and friction modeling using the complex planform from a commercial small turbine, for which limited theoretical and computational guidance is available in selecting the model constants. The model implemented with SI is shown to provide an accurate description of the yaw response of the complex planform. Guidelines are given for the use of wind tunnel tests to determine the model constants for tail fins of any planform.

17 WIND ENERGY↗

DOE OSTI2024

Machine learning informed control systems for extrusion printing processes

Systems and methods for controlling a material extrusion device to extrude a filament of an ink are provided. An extrusion printing control system collects from one or more sensors measurements representing an internal state of material extrusion processing during extrusion of the filament. In addition, the system collects an image of the filament as the filament is extruded. The system applies a classifier to the collected image to generate an image-derived state characterizing the filament. Based on the internal state and the image-derived state, the system estimates a derived state using a model. The system determines control parameters using the model to achieve a desired quality of the filament by minimizing a cost function based on the internal state, the image-derived state, the derived state, and constraints of the material extrusion device. Finally, the system provides the control parameters to a controller of the material extrusion device.

Howell, Brian↗

DOE OSTI2026

Robust Automatic EXAFS First-Shell Fits

Extended X-ray absorption fine structure (EXAFS) is a widely used technique for atomic structure determination. Fourier transformation connects EXAFS in k space and R space. However, determining the appropriate k-range for the transformation can be challenging, but critical for the first-shell fit. In this study, we present an automatic method to determine the k-range using the Larch package and a Python program. The first step is to estimate spectral noise across a series of k-ranges with a fixed minimum value and identify the optimal maximum value in the k-range (k max ). The k max is determined by an empirical noise threshold that marks the point where the noise level in the Fourier transformed spectrum changes dramatically. Using the obtained k max value, the first shell is modeled to determine the minimum k value (k min ) by optimizing the background function through alignment of the spectrum with theory. The optimal k min corresponds to the point of the minimum R-factor, which quantifies the difference between the experimental and fitted spectrum. Our method was tested on various typical datasets and yielded suitable k-ranges for Fourier transformation and accurate first-shell fits. This approach helps avoid unreliable, irreproducible data analysis, especially for noisy data from diluted samples, and enables robust automatic first-shell EXAFS fitting.

EXAFS analysis↗

DOE OSTI2023

BISON analysis of FeCrAl and Zircaloy cladding deformation during simulated BWR cyclic dryout conditions

The BISON fuel performance analysis code has been used to study FeCrAl and Zircaloy cladding performance to better understand cladding deformation behavior in boiling water reactors during conditions caused by anticipated operational occurrences (AOOs) or by anticipated transients without scram (ATWS). BISON models of novel separate effects tests simulating ATWS events predicted mechanical responses to thermal cycling boundary conditions. Here, we assessed FeCrAl and Zircaloy deformation and burst models to determine those best suited for cyclic dryout event analysis. In addition, we compared diametral strain profiles of the post-experiment claddings to simulation results. The capability of the BISON code for modeling BWR AOO or ATWS events is promising.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Heat transfer through a box beam from an impinging hydrogen flame

An analysis was performed to determine whether a hydrogen jet flame impinging on a tunnel ceiling composed of multiple prestressed steel reinforced concrete box beams could result in permanent damage to the tunnel. The lower layer of the concrete box beam was modeled to determine whether heat reaches the steel reinforcing bars and whether spalling could occur. Heat transfer analysis shows that the temperature remains constant at the location of the steel reinforcing bars after 1.3 minutes of impingement and reaches a maximum of 130°C after 5 minutes. However, assuming a constant impingement for 5 minutes is an over estimation due the existing fire model which includes conservative assumptions. Explosive spalling may occur at a thin layer (~0.05 in. at 50 seconds, 0.1 in. at 5 minutes) at the bottom surface of the concrete box beam, but the steel reinforcing bars will not be exposed to the hydrogen flame.

42 ENGINEERING↗

DOE OSTI2023

Airport Infrastructure Planning Using Multi-Stage Stochastic Programming

The Athena project, funded by the Department of Energy, has worked to identify the critical infrastructure at Dallas Fort Worth (DFW) Airport which influences mobility between the airport and the surrounding city of Dallas. Using scalable methods that can leverage HPC resources we have developed a multi-stage stochastic infrastructure expansion model for determining parking and curb modifications to the DFW Airport over a 20-year horizon. Additionally, we have explored the impacts of congestion pricing in conjunction with infrastructure modifications. Our multi-stage stochastic model is implemented using the mpi-sppy software and solved in parallel using progressive hedging on the National Renewable Energy Laboratory's HPC system Eagle. In this talk we present results from solving this model at scale.

airport planning↗

DOE OSTI2021

Dynamics of laser-generated magnetic fields using long laser pulses

In this work, we report on the experimental investigation of magnetic field generation with a half-loop gold sheet coil driven by long-duration (10 ns) and high-power (0.5 TW) laser pulses. The amplitude of the magnetic field was characterized experimentally using proton deflectometry. The field rises rapidly in the first 1 ns of laser irradiation, and then increases slowly and continuously up to 10 ns during further laser irradiation. The transient dynamics of current shape were investigated with a two-dimensional (2D) numerical simulation that included Ohmic heating of the coil and the resultant change of electrical resistivity determined by the coil material temperature. The numerical simulations show rapid heating at the coil edges by current initially localized at the edges. This current density then diffuses to the central part of the sheet coil in a way that depends both on normal current diffusion as well as temporal changes of the coil resistance induced by the Ohmic heating. The measured temporal evolution of the magnetic field is compared with a model that determines a solution to the coil current and voltage that is consistent with a plasma diode model of the drive region and a 2D simulation of current diffusion and dynamic resistance due to Ohmic heating in the laser coil.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2022

Utilizing CellBox to Describe Varied Dynamical Systems

CellBox is a hybrid modeling framework that integrates machine-learning approaches with explicit mathematical modeling to determine the behavior of dynamical systems, intended for determining the complex interactions of molecule and phenotype within a cellular environment using perturbation/response-based data. The original work utilized analysis of a melanoma cell line, SK-Mel-133, which we sought to expand to and evaluate effectiveness in describing other nonlinear systems both within and outside of biology. We evaluated with two models. One model is based upon the free-fall of an object under a velocity-dependent drag force; this produced inherent perturbation/response data. The other model is based upon a form of mass-action molecular kinetics; for this, we ran the simulation to steady state in each species for each perturbation and select the steady state as input for CellBox. We determined that CellBox, outside of its original use case, can still strongly predict even non-biological dynamical systems in the small-perturbation regime, with r = 0.99977 to 0.85123, p < 0.05 for freefall. However, CellBox cannot accurately predict behavior within large-perturbation regimes, p >> 0.05 for freefall. In addition, we performed a naïve comparison of CellBox network structure to known network structure in the mass-action model to assess accuracy.

97 MATHEMATICS AND COMPUTING↗