Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Morris method”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

A new multi-model absolute difference-based sensitivity (MMADS) analysis method to screen non-influential processes under process model and parametric uncertainty

Process-based models have been widely used for hydrologic modeling, and it is a common practice to use sensitivity analysis methods for excluding non-influential hydrologic processes from further investigation and/or model improvement. This study develops a new method called multi-model absolute difference-based sensitivity (MMADS) analysis method to screen non-influential system processes and parameters. MMADS is conceptually similar to the Morris method for addressing parametric uncertainty, but has a unique feature to address both process model uncertainty (i.e., a process may be represented by multiple process models) and process model parameter uncertainty (i.e., parameters associated with a process model are random). MMADS first evaluates absolute differences of a quantity of interest (i.e., a system model output) by varying process models and/or process model parameter values, and then calculates the mean and variance of the differences for investigating process influence. The mean measures overall influence of the process on the quantity of interest, and the variance estimates influence of nonlinear effects of the process and/or its interactions with other processes. MMADS is an extension of the Morris method from a parameter space to a joint parameter-model space for explicitly addressing both process model uncertainty and model parameter uncertainty. The performance of MMADS is evaluated by using two numerical experiments. One experiment is based on Sobol’s G*-function with ten product elements, and has analytical solutions of the MMADS mean and variance of absolute differences. The other experiment is for groundwater flow modeling which considers three processes (i.e., recharge, geology, and snowmelt) that interact with each other. Finally, results indicate that MMADS is computationally efficient and can identify non-influential processes of complex hydrological systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Identifying high-impact and high-uncertainty parameters in MiniFuel model predictions

The MiniFuel irradiation platform at Oak Ridge National Laboratory's High Flux Isotope Reactor (HFIR) is a flexible, high-throughput separate effects test capability. Finite element thermal models are relied upon to design MiniFuel experiments and to achieve experimental objectives. Recent reports show good agreement in the model prediction of target fuel temperatures, but as the capability of the experiments is extended to higher temperatures, the uncertainty in the model predictions must be quantified. To that end, high-impact, high-uncertainty parameters that contribute the most uncertainty to the model are identified. The uncertainty quantification was accomplished through a series of screening and sensitivity analyses. The first analysis utilizes the method of Morris to perform a computationally efficient preliminary screening that considers uncertainty in a large number of the model inputs. The most important parameters identified in the Morris screening study were then considered in a Sobol sensitivity analysis that more robustly ranks and quantifies the uncertainty associated with each parameter. From these analyses, it was determined that thermal contact conductance between components is the parameter that contributes the highest uncertainty. The estimated uncertainty of the MiniFuel model fuel temperature predictions is ±80 °C in the removable beryllium and ±40 °C in the vertical experiment facilities. In conclusion, the framework established by the series of sensitivity analyses presented herein could easily be adapted to fit the needs of accelerated fuel qualification processes.

Fuel, Irradiation↗
DOE OSTI2026

Physiological Controls on Carbon Fluxes and Biomass Production in Miscanthus: Insights From a Process‐Based Agroecosystem Model

Biomass crops serve as essential feedstocks for renewable energy and bioproducts and play a critical role in achieving lower emissions in the transportation sector. However, dedicated perennial biomass crops such as Miscanthus × giganteus (Miscanthus) remain underrepresented in process-based agroecosystem models, limiting robust evaluation of their economic and environmental performance. In this study, we developed a data-constrained representation of the sterile triploid Miscanthus (IL clone) within the process-based model ecosys, integrating global sensitivity analysis, ensemble simulation, and parameter calibration. Planting, harvesting, and fertilization practices consistent with field management were incorporated, and phenology was constrained using PhenoCam-derived Green Chromatic Coordinate (GCC) data. Using the Morris global sensitivity analysis method, we identified 11 key physiological parameters governing plant carbon, water, and nutrient relations, particularly processes associated with CO 2 assimilation. We then conducted ensemble simulations by perturbing these parameters and calibrated the model against eddy covariance fluxes and field-measured biomass. Building on the calibrated operating state, parameter-response analyses show that different photosynthetic processes influence productivity in different ways. Protein allocation determines whether productivity increases toward a higher level, whereas electron transport capacity controls additional gains once protein allocation approaches saturation. These findings demonstrate that parameter importance depends on physiological context and on which photosynthetic processes remain limiting. Calibration and validation against observations show that ecosys can reliably reproduce carbon and water fluxes, as well as both above- and belowground biomass, with post-calibration GPP R 2 improving from 0.67 to 0.95 during the calibration period and remaining high during validation (R 2 = 0.95). These results provide a mechanistic foundation for regional simulations and sustainable bioenergy assessments.

ecosys↗
DOE OSTI2026

Data for Physiological Controls on Carbon Fluxes and Biomass Production in Miscanthus: Insights From a Process- Based Agroecosystem Model

Biomass crops serve as essential feedstocks for renewable energy and bioproducts and play a critical role in achieving lower emissions in the transportation sector. However, dedicated perennial biomass crops such as Miscanthus × giganteus (Miscanthus) remain underrepresented in process- based agroecosystem models, limiting robust evaluation of their economic and environmental performance. In this study, we developed a data- constrained representation of the sterile triploid Miscanthus (IL clone) within the process- based model ecosys, integrating global sensitivity analysis, ensemble simulation, and parameter calibration. Planting, harvesting, and fertilization practices consistent with field management were incorporated, and phenology was constrained using PhenoCam- derived Green Chromatic Coordinate (GCC) data. Using the Morris global sensitivity analysis method, we identified 11 key physiological parameters governing plant carbon, water, and nutrient relations, particularly processes associated with CO2 assimilation. We then conducted ensemble simulations by perturbing these parameters and calibrated the model against eddy covariance fluxes and field- measured biomass. Building on the calibrated operating state, parameter- response analyses show that different photosynthetic processes influence productivity in different ways. Protein allocation determines whether productivity increases toward a higher level, whereas electron transport capacity controls additional gains once protein allocation approaches saturation. These findings demonstrate that parameter importance depends on physiological context and on which photosynthetic processes remain limiting. Calibration and validation against observations show that ecosys can reliably reproduce carbon and water fluxes, as well as both above- and belowground biomass, with post- calibration GPP R2 improving from 0.67 to 0.95 during the calibration period and remaining high during validation (R2 = 0.95). These results provide a mechanistic foundation for regional simulations and sustainable bioenergy assessments.

Carbon↗
DOE OSTI2022

Assessing the Sensitivity of the Tropical Cyclone Boundary Layer to the Parameterization of Momentum Flux in the Community Earth System Model

Recent studies have demonstrated that high-resolution (~25 km) Earth System Models (ESMs) have the potential to skillfully predict tropical cyclone (TC) occurrence and intensity. However, biases in ESM TCs still exist, largely due to the need to parameterize processes such as boundary layer (PBL) turbulence. Building on past studies, we hypothesize that the depiction of the TC PBL in ESMs is sensitive to the configuration of the PBL parameterization scheme, and that the targeted perturbation of tunable parameters can reduce biases. The Morris one-at-a-time (MOAT) method is implemented to assess the sensitivity of the TC PBL to tunable parameters in the PBL scheme in an idealized configuration of the Community Atmosphere Model, version 6 (CAM6). The MOAT method objectively identifies several parameters in an experimental version of the Cloud Layers Unified by Binormals (CLUBB) scheme that appreciably influence the structure of the TC PBL. We then perturb the parameters identified by the MOAT method within a suite of CAM6 ensemble simulations and find a reduction in model biases compared to observations and a high-resolution, cloud-resolving model. Importantly, we demonstrate that the high-sensitivity parameters are tied to PBL processes that reduce turbulent mixing and effective eddy diffusivity, and that in CAM6 these parameters alter the TC PBL in a manner consistent with past modeling studies. In this way, we provide an initial identification of process-based input parameters that, when altered, have the potential to improve TC predictions by ESMs.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Sensitive parameter identification and uncertainty quantification for the stability of pipeline conveying fluid

In this study, several uncertainty quantification and sensitivity analysis methods are used to determine the most sensitive geometric and material input parameters of a cantilevered pipeline conveying fluid when uncertainty is introduced to the system at the onset of instability. The full nonlinear equations of motion are modeled using the extended Hamilton’s principle and then discretized using Galerkin’s method. A parametric study is first performed, and the Morris elementary effects are calculated to obtain a preliminary understanding of how the onset speed changes when each parameter is introduced to a ± 5% uncertainty. Then, four different input uncertainty distributions, mainly, uniform and Gaussian distribution, are chosen to investigate how input distributions affect uncertainty in the output. A convergence analysis is used to determine the number of samples needed to maintain simulation accuracy while saving the most computational time. Then, Monte Carlo simulations are run, and the output distributions for each input distribution at ± 1%, ± 3% and ± 5% input uncertainty range are found and discussed. Additionally, the Pearson correlation coefficients are evaluated for different uncertainty ranges. A final Monte Carlo study is performed in which single parameters are held constant while all others still have uncertainty. Overall, the flow speed at the onset of instability is the most sensitive to changes in the outer diameter of the pipe.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Comparison of planetary boundary layer height from ceilometer with ARM radiosonde data

Abstract. Ceilometer measurements of aerosol backscatter profiles have been widely used to provide continuous planetary boundary layer height (PBLHT) estimations. To investigate the robustness of ceilometer-estimated PBLHT under different atmospheric conditions, we compared ceilometer- and radiosonde-estimated PBLHTs using multiple years of U.S. Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) ceilometer and balloon-borne sounding data at ARM fixed-location atmospheric observatories and from ARM mobile facilities deployed around the world for various field campaigns. These observatories cover from the tropics to the polar regions and over both ocean and land surfaces. Statistical comparisons of ceilometer-estimated PBLHTs from the Vaisala CL31 ceilometer data with radiosonde-estimated PBLHTs from the ARM PBLHT-SONDE Value-added Product (VAP) are performed under different atmospheric conditions including stable and unstable atmospheric boundary layer, low-level cloud-free conditions, and cloudy conditions at these ARM observatories. Under unstable conditions, good comparisons are found between ceilometer- and radiosonde-estimated PBLHTs at ARM low- and mid-latitude land observatories. However, it is still challenging to obtain reliable PBLHT estimations over ocean surfaces even using radiosonde data. Under stable conditions, ceilometer- and radiosonde-estimated PBLHTs have weak correlations. We compare different PBLHT estimations utilizing the Heffter, the Liu–Liang, and the bulk Richardson number methods applied to radiosonde data with ceilometer-estimated PBLHT. We find that ceilometer-estimated PBLHT compares better with the Liu–Liang method under unstable conditions and compares better with the bulk Richardson number method under stable conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Synthesis and thermodynamics of uranium-incorporated α-Fe 2 O 3 nanoparticles

Hematite nanoparticles were synthesized with U(VI) in circumneutral water through a coprecipitation and hydrothermal treatment process. XRD, TEM, and EXAFS analyses reveal that uranium may aggregate along grain boundaries and occupy Fe sites within hematite. The described synthesis method produces crystalline, single-phase iron oxide nanoparticles absent of surface-bound uranyl complexes. EXAFS data were comparable to spectra from existing studies whose syntheses were more representative of naturally occurring, extended aging processes. Herein this work provides and validates an accelerated method of synthesizing uranium-immobilized iron oxide nanoparticles for further mechanistic studies. High temperature oxide melt solution calorimetry measurements were performed to calculate the thermodynamic stability of uranium-incorporated iron oxide nanoparticles. Increasing uranium content within hematite resulted in more positive formation enthalpies. Standard formation enthalpies of U x Fe 2–2x O 3 were as high as 76.88 ± 2.83 kJ/mol relative to their binary oxides, or -764.04 ± 3.74 kJ/mol relative to their constituent elements, at x = 0.037. Data on the thermodynamic stability of uranium retention pathways may assist in predicting waste uranyl remobilization, as well as in developing more effective methods to retain uranium captured from aqueous environments.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Development and experimental verification of an active noise cancellation (ANC) method for non-invasive sensing

Ultrasonic wave propagation may be used in non-invasive sensing to identify the properties of a fluid in a container. An ultrasonic excitation can be coupled into the liquid or the gas through the walls of the container and detected by a receiver on the other side. However, in addition to this fluid path, another signal will also propagate through the structure and is referred as noise in this study. For most materials, this structural noise will be of higher amplitude compared to the fluid signal due to the large acoustic impedance difference between the solid container and the fluid. In this paper, an active noise cancellation (ANC) method is developed and demonstrated on an experimental setup. It was shown that this method cancels more than 90% of the noise signal in a given time window of interest. A sequential version of this cancellation method is also proposed, and it was shown that its performance improves as the number of canceling excitations is increased. The performance was found to be independent of the location of canceling transducers and position of the silence window but dependent on the width of the silence window. As an example, the ANC method was applied to noise cancellation of a signal travelling through water in a stainless-steel container. The signal-to-noise ratio (SNR) of the water signal in the window of interest was increased by 148% from 13.98 dB to 34.68 dB by using the proposed ANC method.

42 ENGINEERING↗
DOE OSTI2024

Cambium 2023 Scenario Descriptions and Documentation

The National Renewable Energy Laboratory's (NREL's) Cambium data sets are annually released sets of simulated hourly emission, cost, and operational data for a range of modeled futures of the U.S. electric sector with metrics designed to be useful for long-term decision- making. The 2023 Cambium data set is the fourth annual release. The data sets are a companion product to NREL's Standard Scenarios, which are likewise released annually and are a set of projections of how the U.S. electric sector could evolve across a suite of different potential futures, but covering more scenarios with less temporal granularity (Gagnon et al. 2024). Information about Cambium and related publications can be found at https://www.nrel.gov/analysis/cambium.html, and the Cambium data sets can be viewed and downloaded at https://scenarioviewer.nrel.gov/. In this documentation, we describe Cambium 2023's scenarios, define the metrics, and document the Cambium-specific methods for calculating those metrics.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

A deep learning-accelerated data assimilation and forecasting workflow for commercial-scale geologic carbon storage

Fast assimilation of monitoring data to update forecasts of pressure buildup and carbon dioxide (CO 2 ) plume migration under geologic uncertainties is a challenging problem in geologic carbon storage. The high computational cost of data assimilation with a high-dimensional parameter space impedes fast decision-making for commercial-scale reservoir management. We propose to leverage physical understandings of porous medium flow behavior with deep learning techniques to develop a fast data assimilation-reservoir response forecasting workflow. Applying an Ensemble Smoother Multiple Data Assimilation (ES-MDA) framework, the workflow updates geologic properties and predicts reservoir performance with quantified uncertainty from pressure history and CO 2 plumes interpreted through seismic inversion. As the most computationally expensive component in such a workflow is reservoir simulation, we developed surrogate models to predict dynamic pressure and CO 2 plume extents under multi-well injection. The surrogate models employ deep convolutional neural networks, specifically, a wide residual network and a residual U-Net. The workflow is validated against a flat threedimensional reservoir model representative of a clastic shelf depositional environment. Intelligent treatments are applied to bridge between quantities in a true-3D reservoir model and those in a single-layer reservoir model. The workflow can complete history matching and reservoir forecasting with uncertainty quantification in less than one hour on a mainstream personal workstation.

25 ENERGY STORAGE↗
DOE OSTI2022

Correction and calibration of atmospheric impact observations in GOES GLM data

The Earth's atmosphere is impacted daily by both meteoroids and artificial objects. Calibrated observations of the emitted light at sufficiently high sampling rates can enable or improve the estimation of impactor attributes such as size, cohesion, trajectory, and composition, but are difficult to obtain owing to the unpredictability, brevity, and high dynamic (brightness) range of impacts. Ground-based camera systems have successfully monitored small regions of the atmosphere at video frame rates and with limited radiometric capabilities, but most impacts occur over the 70% of the Earth's surface covered by water and are therefore missed by these networks. The Geostationary Lightning Mapper (GLM) instruments aboard Geostationary Operational Environmental Satellites 16 and 17 provide near-hemispherical coverage at 500 frames per second. These data have been shown to contain the signatures of many independently confirmed impacts, often from both viewing angles simultaneously, and constitute an observational resource that is currently unparalleled in the public domain. NASA's Asteroid Threat Assessment Project has implemented an automated impact detection pipeline that processes data from GLM daily. Given a detected impact, the GLM data contain a wealth of information for use in quantitative follow-up analyses. However, impact events differ from lightning in ways that violate key assumptions built into GLM's design. The result is that GLM's onboard processing introduces errors into pixel observations of impact events and the calibrated energies near the periphery of the detector may be substantially overestimated. We present methods for mitigating these and other issues to produce a data product more suitable for impact analyses than the existing GLM lightning product.

58 GEOSCIENCES↗
DOE OSTI2022

Modeling protected species distributions and habitats to inform siting and management of pioneering ocean industries: A case study for Gulf of Mexico aquaculture

Marine Spatial Planning (MSP) provides a process that uses spatial data and models to evaluate environmental, social, economic, cultural, and management trade-offs when siting (i.e., strategically locating) ocean industries. Aquaculture is the fastest-growing food sector in the world. The United States (U.S.) has substantial opportunity for offshore aquaculture development given the size of its exclusive economic zone, habitat diversity, and variety of candidate species for cultivation. However, promising aquaculture areas overlap many protected species habitats. Aquaculture siting surveys, construction, operations, and decommissioning can alter protected species habitat and behavior. Additionally, aquaculture-associated vessel activity, underwater noise, and physical interactions between protected species and farms can increase the risk of injury and mortality. In 2020, the U.S. Gulf of Mexico was identified as one of the first regions to be evaluated for offshore aquaculture opportunities as directed by a Presidential Executive Order. We developed a transparent and repeatable method to identify aquaculture opportunity areas (AOAs) with the least conflict with protected species. First, we developed a generalized scoring approach for protected species that captures their vulnerability to adverse effects from anthropogenic activities using conservation status and demographic information. Next, we applied this approach to data layers for eight species listed under the Endangered Species Act, including five species of sea turtles, Rice’s whale, smalltooth sawfish, and giant manta ray. Next, we evaluated four methods for mathematically combining scores (i.e., Arithmetic mean, Geometric mean, Product, Lowest Scoring layer) to generate a combined protected species data layer. The Product approach provided the most logical ordering of, and the greatest contrast in, site suitability scores. Finally, we integrated the combined protected species data layer into a multi-criteria decision-making modeling framework for MSP. This process identified AOAs with reduced potential for protected species conflict. These modeling methods are transferable to other regions, to other sensitive or protected species, and for spatial planning for other ocean-uses.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Accurate prediction of short-range order and its effect on thermodynamic, structural, and electronic properties of disordered alloys: exemplified in archetypical Cu 3 Au

Electronic-structure methods based on density-functional theory (DFT) were used to quantify the effect of chemical short-range order (SRO) on thermodynamic, structural, and electronic properties of archetypal face-centered-cubic (fcc) Cu3Au alloy. We showed that SRO can be tuned to alter bonding and lattice dynamics (i.e., phonons) and detail how these properties are changed with SRO. Thermodynamically favorable SRO significantly improved the phase stability of fcc Cu3Au from -0.0343 eV-atom -1 to –0.0682 eV-atom -1 . We used our DFT-based linear-response theory to predict SRO and its electronic origin, and accurately estimate the observed transition temperature, ordering instability (L1 2 ), and Warren-Cowley SRO parameters, in agreement with experiments. The accurate prediction of real-space SRO gives an edge over computationally and resource intensive approaches such as monte-carlo methods or experiments, which will enable large scale molecular dynamic simulations by providing supercells with optimized SRO. Here we also analyzed phonon dispersion and estimated the vibrational entropy change (from 9kB at 300 K to 6kB at 100 K) in fcc Cu3Au. We established from SRO analysis that exclusion of chemical interactions may lead to a skewed view of true properties in chemically complex alloys. The first-principles methods described in this work are generally applicable to any arbitrary solid-solution alloys, including multi-principal-element alloys, therefore, holds promise for designing technologically useful materials.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Rare isotope-containing diamond colour centres for fundamental symmetry tests

Detecting a non-zero electric dipole moment in a particle would unambiguously signify physics beyond the Standard Model. A potential pathway towards this is the detection of a nuclear Schiff moment, the magnitude of which is enhanced by the presence of nuclear octupole deformation. However, due to the low production rate of isotopes featuring such ‘pear-shaped’ nuclei, capturing, detecting and manipulating them efficiently is a crucial prerequisite. Incorporating them into synthetic diamond optical crystals can produce defects with defined, molecule-like structures and isolated electronic states within the diamond band gap, increasing capture efficiency, enabling repeated probing of even a single atom and producing narrow optical linewidths. In this study, we used density functional theory to investigate the formation, structure and electronic properties of crystal defects in diamond containing 229 ⁢Pa, a rare isotope that is predicted to have an exceptionally strong nuclear octupole deformation. In addition, we identified and studied stable lanthanide-containing defects with similar electronic structures as non-radioactive proxies to aid in experimental methods. Our findings hold promise for the existence of such defects and can contribute to the development of a quantum information processing-inspired toolbox of techniques for studying rare isotopes.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2023

Transforming microseismic clouds into near real-time visualization of the growing hydraulic fracture

SUMMARY Microseismic observations during unconventional reservoir stimulation are typically seen as a proxy for clusters of hydraulic fractures and the extent of the stimulated reservoir. Such straightforward interpretation is often misleading and fails to provide a physically reasonable image of the fracturing process. This paper demonstrates the application of a physics-based machine learning algorithm which enables a rapid and accurate fracture mapping from the microseismic data. Our training and validation data set relies on a history-matched geomechanical modelling workflow implemented in GEOS software for the Hydraulic Fracturing Test Site 1 (HFTS-1) project. For this study we augmented the simulated fracture growth through geostatistical modelling of induced seismicity, so that the synthetic microseismic catalogue matches the main statistical properties of the field observations. We formulated the problem of mapping the actual fracture in the clutter of events to parallel common video segmentation workflows: several past video frames (microseismic density snapshots) are passed through a deep convolutional network to classify whether a given voxel is associated with a fracture or intact rock. We found that for accurate fracture mapping, the network’s input and architecture must be augmented to incorporate the fluid injection parameters (pressure, rate, concentration of proppant, and location of the perforation within the cluster). The error rate for the network reached as little as 10 per cent of the fracture area, while a conventional microseismic interpretation approach yielded ∼300 per cent. Our approach also yields must faster predictions than conventional methods (minutes instead of weeks), and could enable engineers to make rapid decisions regarding engineering parameters (pumping rate, viscosity) in real time during stimulation.

58 GEOSCIENCES↗