Search NASA⌕ Search

SEARCH · Search NASA

Results for “variable order”

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 127 records · Page 7

A reforecasting-based dynamic reserve estimation for variable renewable generation and demand uncertainty

The installed capacity of renewables-based energy sources has been increasing in traditional power systems. In order to accommodate the increased variability and uncertainty associated with the deeper penetration of renewable sources like solar and wind, adjusted amounts of dynamic reserve are needed. Although probabilistic dynamic reserve estimation methods have been previously developed, most of them consider the uncertainty to be represented by parametric density functions that tend to perform poorly under extreme events and, moreover, neglect uncertainty introduced by the forecasting model itself. Toward addressing these limitations, this work presents, for the first time, a dynamic reserve estimation method for flexibility that incorporates nonparametric density estimation and a machine learning based reforecasting to provide a day-ahead prediction of the mean and spread of uncertainty around the base forecast. The prediction is, in turn, used to estimate the up and down reserve relative to the base forecast. Here, the present method takes various endogenous and exogenous features, including the calendar variables, as input to estimate the day-ahead reserve. Using a combination of reforecasting and dynamic reserve estimation techniques, the method is shown to adjust better to the dynamic nature of reserve requirements providing only what is needed to accommodate the expected deviations. Considering California Independent System Operator (CAISO) solar, wind and load data over an 18 month period, up to 67% reduction in the amount of reserve capacity needed for a one day reserve and reserve penalty for solar uncertainty is demonstrated. Additionally, the risk of reserve insufficiency in meeting the net demand is reduced by 20% with the proposed method.

14 SOLAR ENERGY↗

Downscaled hyper-resolution (400 m) gridded datasets of daily precipitation and temperature (2008–2019) for the East–Taylor subbasin (western United States)

Abstract. High-resolution gridded datasets of meteorological variables are needed in order to resolve fine-scale hydrological gradients in complex mountainous terrain. Across the United States, the highest available spatial resolution of gridded datasets of daily meteorological records is approximately 800 m. This work presents gridded datasets of daily precipitation and mean temperature for the East–Taylor subbasin (in the western United States) covering a 12-year period (2008–2019) at a high spatial resolution (400 m). The datasets are generated using a downscaling framework that uses data-driven models to learn relationships between climate variables and topography. We observe that downscaled datasets of precipitation and mean temperature exhibit smoother spatial gradients (while preserving the spatial variability) when compared to their coarser counterparts. Additionally, we also observe that when downscaled datasets are upscaled to the original resolution (800 m), the mean residual error is almost zero, ensuring no bias when compared with the original data. Furthermore, the downscaled datasets are observed to be linearly related to elevation, which is consistent with the methodology underlying the original 800 m product. Finally, we validate the spatial patterns exhibited by downscaled datasets via an example use case that models lidar-derived estimates of snowpack. The presented dataset constitutes a valuable resource to resolve fine-scale hydrological gradients in the mountainous terrain of the East–Taylor subbasin, which is an important study area in the context of water security for the southwestern United States and Mexico. The dataset is publicly available at https://doi.org/10.15485/1822259 (Mital et al., 2021).

54 ENVIRONMENTAL SCIENCES↗

SCEC/USGS Community Stress-Drop Validation Study: How Spectral Fitting Approaches Influence Measured Source Parameters

Spectral source parameters used to estimate an earthquake’s stress drop (⁠Δσ⁠) can vary significantly across measurement approaches. The Statewide California Earthquake Center/U.S. Geological Survey Community Stress‐Drop Validation Study was initiated to compare source parameter estimates, focusing initially on a dataset from the 2019 Ridgecrest earthquake sequence. As part of that validation effort, here we focus on one potential source of uncertainty: whether spectral fitting approaches alone, applied to a common set of spectra from the 2019 Ridgecrest sequence result in different source parameter estimates. By using a common set of benchmark spectra analyzed across a consistent frequency band of 1–40 Hz, we eliminate many sources of variability. A subgroup of validation study participants volunteered to estimate the low‐frequency displacement (⁠Ω 0 ⁠) and corner frequency (⁠ƒ c ⁠) by fitting a smooth function to benchmark displacement spectra. Participants used linear‐ or log‐sampled spectra, assumed a Brune or Boatwright spectral model, and applied different misfit criteria. We compare 17 approaches used to estimate ⁠Ω 0 ⁠, ƒ c ⁠, and Δσ for 54 earthquake spectra. Our results reveal that 35% of events have Δσ estimates within a factor of two, whereas others exhibit variations exceeding an order of magnitude. The variability in and can largely be attributed to whether a spectrum is consistent with the smooth function of an idealized simple crack model. The trade‐off between Ω 0 and ƒ c may be more pronounced when using linearly sampled spectra, as higher frequency spectral bumps control the fits. As expected, methods that assumed a Boatwright model tended to have lower Ω 0 and somewhat higher ƒ c compared to those assuming a Brune model, although resulting Δσ estimates are similar. Finally, when compared to the overall validation study results, the fitting approach alone may account for between 5% and 90% (25% on average) of the total variability in spectral Δσ⁠.

58 GEOSCIENCES↗

Detection of variable VHE γ -ray emission from the extra-galactic γ -ray binary LMC P3

Context. Recently, the high-energy (HE, 0.1–100 GeV) γ-ray emission from the object LMC P3 in the Large Magellanic Cloud (LMC) has been discovered to be modulated with a 10.3-day period, making it the first extra-galactic γ-ray binary. Aim. This work aims at the detection of very-high-energy (VHE, >100 GeV) γ-ray emission and the search for modulation of the VHE signal with the orbital period of the binary system. Methods. LMC P3 has been observed with the High Energy Stereoscopic System (H.E.S.S.); the acceptance-corrected exposure time is 100 h. The data set has been folded with the known orbital period of the system in order to test for variability of the emission. Results. VHE γ-ray emission is detected with a statistical significance of 6.4 σ. The data clearly show variability which is phase-locked to the orbital period of the system. Periodicity cannot be deduced from the H.E.S.S. data set alone. The orbit-averaged luminosity in the 1–10 TeV energy range is (1.4 ± 0.2) × 1035 erg s -1 . A luminosity of (5 ± 1) × 1035 erg s -1 is reached during 20% of the orbit. HE and VHE γ-ray emissions are anti-correlated. In conclucion, LMC P3 is the most luminous γ-ray binary known so far.

79 ASTRONOMY AND ASTROPHYSICS↗

A Combined Machine Learning and High-Energy X-ray Diffraction Approach to Understanding Liquid and Amorphous Metal Oxides

Determining the structure-property relations of liquid and amorphous metal oxides is challenging, due to their variable short-range order and polyhedral connectivity. To predict chemically realistic structures, we have developed a Machine Learned, Gaussian Approximation Potential (GAP) for HfO 2 , with a focus on enhanced sampling of the training database and accurate density functional theory calculations. By using training datasets for the GAP model at the level of Density Functional Theory-Strongly Constrained and Appropriately Normed (DFT-SCAN) level of theory, our results show that the topology of both the low viscosity liquid and the amorphous form are dominated by edge-shared chains and small corner-shared rings of polyhedra. This topology is shown to be consistent with the structure of other liquid and amorphous transition metal oxides of variable ion size, such as TiO 2 and ZrO 2 . Current limitations of the ML-GAP modeling method for obtaining glass structures and future perspectives are also discussed.

36 MATERIALS SCIENCE↗

An optimization-based approach to parameter learning for fractional type nonlocal models

Nonlocal operators of fractional type are a popular modeling choice for applications that do not adhere to classical diffusive behavior; however, one major challenge in nonlocal simulations is the selection of model parameters. In this study we propose an optimization-based approach to parameter identification for fractional models with an optional truncation radius. We formulate the inference problem as an optimal control problem where the objective is to minimize the discrepancy between observed data and an approximate solution of the model, and the control variables are the fractional order and the truncation length. For the numerical solution of the minimization problem we propose a gradient-based approach, where we enhance the numerical performance by an approximation of the bilinear form of the state equation and its derivative with respect to the fractional order. Several numerical tests in one and two dimensions illustrate the theoretical results and show the robustness and applicability of our method.

97 MATHEMATICS AND COMPUTING↗

A coupled multipoint stress–multipoint flux mixed finite element method for the Biot system of poroelasticity

In this work, we present a mixed finite element method for a five-field formulation of the Biot system of poroelasticity that reduces to a cell-centered pressure–displacement system on simplicial and quadrilateral grids. A mixed stress–displacement–rotation formulation for elasticity with weak stress symmetry is coupled with a mixed velocity–pressure Darcy formulation. The spatial discretization is based on combining the multipoint stress mixed finite element (MSMFE) method for elasticity and the multipoint flux mixed finite element (MFMFE) method for Darcy flow. It uses the lowest order Brezzi–Douglas–Marini mixed finite element spaces for the poroelastic stress and Darcy velocity, piecewise constant displacement and pressure, and continuous piecewise linear or bilinear rotation. A vertex quadrature rule is applied to the velocity, stress, and stress–rotation bilinear forms, which block-diagonalizes the corresponding matrices and allows for local velocity, stress, and rotation elimination. This leads to a cell-centered positive-definite system for pressure and displacement at each time step. We perform error analysis for the semidiscrete and fully discrete formulations, establishing first order convergence for all variables in their natural norms. The numerical tests confirm the theoretical convergence rates and illustrate the locking-free property of the method.

42 ENGINEERING↗

Adjoint DSMC for nonlinear spatially-homogeneous Boltzmann equation with a general collision model

We derive an adjoint method for the Direct Simulation Monte Carlo (DSMC) method for the spatially homogeneous Boltzmann equation with a general collision law. This generalizes our previous results in Caflisch et al., which was restricted to the case of Maxwell molecules, for which the collision rate is constant. The main difficulty in generalizing the previous results is that a rejection sampling step is required in the DSMC algorithm in order to handle the variable collision rate. We find a new term corresponding to the so-called score function in the adjoint equation and a new adjoint Jacobian matrix capturing the dependence of the collision parameter on the velocities. The new formula works for a much more general class of collision models.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Radon kinetics in a natural indoor radon chamber

An unusual 180 m 3 storage room in the basement of a two-story laboratory building is unventilated, and separated from occupiable rooms by double steel doors. The space completely borders on soil through the concrete floor and two of its concrete walls. The room also contains a separate inner chamber with 1 m thick concrete walls designed to damp vibrations in the room above it. The space boasts a relatively high radon level, 1083 Bq m –3 , which varies with local outdoor environmental conditions. Measurements were made of radon concentrations at various locations and heights within the facility. More than a year of continuous radon concentration data corresponding to a single location are also available, along with measurements of indoor and outdoor pressure, temperature, and humidity. Data were also collected with as many as five fans placed in different locations and cycled on for variable time periods. First order linear kinetic models were created to explain the observed approaches to steady state due to changing conditions and wash-out resulting from intentional ventilation. Results demonstrate a good fit between changes in the radon concentration level and the developed compartmental models. Furthermore, no significant differences were observed between radon concentration at different locations or heights in the chamber.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

From calibration to parameter learning: Harnessing the scaling effects of big data in geoscientific modeling

Abstract The behaviors and skills of models in many geosciences (e.g., hydrology and ecosystem sciences) strongly depend on spatially-varying parameters that need calibration. A well-calibrated model can reasonably propagate information from observations to unobserved variables via model physics, but traditional calibration is highly inefficient and results in non-unique solutions. Here we propose a novel differentiable parameter learning (dPL) framework that efficiently learns a global mapping between inputs (and optionally responses) and parameters. Crucially, dPL exhibits beneficial scaling curves not previously demonstrated to geoscientists: as training data increases, dPL achieves better performance, more physical coherence, and better generalizability (across space and uncalibrated variables), all with orders-of-magnitude lower computational cost. We demonstrate examples that learned from soil moisture and streamflow, where dPL drastically outperformed existing evolutionary and regionalization methods, or required only ~12.5% of the training data to achieve similar performance. The generic scheme promotes the integration of deep learning and process-based models, without mandating reimplementation.

54 ENVIRONMENTAL SCIENCES↗

The role of diffusion-driven pure climb creep on the rheology of bridgmanite under lower mantle conditions

The viscosity of Earth’s lower mantle is poorly constrained due to the lack of knowledge on some fundamental variables that affect the deformation behaviour of its main mineral phases. This study focuses on bridgmanite, the main lower mantle constituent, and assesses its rheology by developing an approach based on mineral physics. Following and revising the recent advances in this field, pure climb creep controlled by diffusion is identified as the key mechanism driving deformation in bridgmanite. The strain rates of this phase under lower mantle pressures, temperatures and stresses are thus calculated by constraining diffusion and implementing a creep theoretical model. The viscosity of MgSiO3 bridgmanite resulting from pure climb creep is consequently evaluated and compared with the viscosity profiles available from the literature. We show that the inferred variability of viscosity in these profiles can be fully accounted for with the chosen variables of our calculation, e.g., diffusion coefficients, vacancy concentrations and applied stresses. A refinement of these variables is advocated in order to further constrain viscosity and match the observables.

58 GEOSCIENCES↗

Recovery of Damaged Information and the Out-of-Time-Ordered Correlators

The evolution with a complex Hamiltonian generally leads to information scrambling. A time-reversed dynamics unwinds this scrambling and thus leads to the original information recovery. In this work, we show that if the scrambled information is, in addition, partially damaged by a local measurement, then such a damage can still be treated by application of the time-reversed protocol. This information recovery is described by the long-time saturation value of a certain out-of-time-ordered correlator of local variables. We also propose a simple test that distinguishes between quantum and reversible classical chaotic information scrambling.

97 MATHEMATICS AND COMPUTING↗

Remote Distributed Vibration Sensing Through Opaque Media Using Permanent Magnets

Vibration sensing is critical for a variety of applications from structural fatigue monitoring to understanding the modes of airplane wings. In particular, remote sensing techniques are needed for measuring the vibrations of multiple points simultaneously, assessing vibrations inside opaque metal vessels, and sensing through smoke clouds and other optically challenging environments. Here, in this paper, we propose a method which measures high-frequency displacements remotely using changes in the magnetic field generated by permanent magnets. We leverage the unique nature of vibration tracking and use a calibrated local model technique developed specifically to improve the frequency-domain estimation accuracy. The results show that two-dimensional local models surpass the dipole model in tracking high-frequency motions. A theoretical basis for understanding the effects of electronic noise and error due to correlated variables is generated in order to predict the performance of experiments prior to implementation. Simultaneous measurements of up to three independent vibrating components are shown. The relative accuracy of the magnet-based displacement tracking with respect to the video tracking ranges from 40 to 190 μm when the maximum displacements approach ±5 mm and when sensor-to-magnet distances vary from 25 to 36 mm. Finally, vibration sensing inside an opaque metal vessel and mode shape changes due to damage on an aluminum beam are also studied using the wireless permanent-magnet vibration sensing scheme.

42 ENGINEERING↗

Assessment of Aquifer Baseline Conditions at the Idaho National Laboratory Remote-Handled Low-Level Waste Disposal Facility

This document provides an analysis of radiologic concentrations in the aquifer and vadose zone near Idaho National Laboratory’s (INL’s) Remote-Handled Low-Level Waste (RH-LLW) disposal facility. This document identifies pre-existing sources of radiologic contaminants that could be transported to the aquifer near the RH-LLW disposal facility based on releases from nearby facility operations. It contains a summary of radionuclide concentration data in the upstream perched water and in the aquifer upstream and downstream of the RH LLW disposal facility in order to establish known variability from these pre existing radiologic sources. Using this historical data allows meaningful aquifer monitoring targets to be established for the RH-LLW disposal facility monitoring system (i.e., vadose zone wells and aquifer wells) to allow potential releases from the RH-LLW disposal facility to be distinguished from the pre existing sources of contaminants that are likely to be detected in the monitoring system.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

An optimization-based approach to parameter learning for fractional type nonlocal models

Nonlocal operators of fractional type are a popular modeling choice for applications that do not adhere to classical diffusive behavior; however, one major challenge in nonlocal simulations is the selection of model parameters. In this work we propose an optimization-based approach to parameter identification for fractional models with an optional truncation radius. We formulate the inference problem as an optimal control problem where the objective is to minimize the discrepancy between observed data and an approximate solution of the model, and the control variables are the fractional order and the truncation length. For the numerical solution of the minimization problem we propose a gradient-based approach, where we enhance the numerical performance by an approximation of the bilinear form of the state equation and its derivative with respect to the fractional order. Several numerical tests in one and two dimensions illustrate the theoretical results and show the robustness and applicability of our method.

97 MATHEMATICS AND COMPUTING↗

In-Situ Monitoring Assisted Large-Scale Additive Manufacturing of Mild Steel and 316L Alloys for Nuclear Application

The Advanced Materials and Manufacturing Technologies (AMMT) Program is aimed at developing cross-cutting technologies in support of a broad range of nuclear reactor parts, and to maintain U.S. leadership in materials and manufacturing technologies for nuclear energy applications. The overarching vision of the AMMT program is to accelerate the development, qualification, demonstration and deployment of advanced materials and manufacturing technologies to enable reliable and economical nuclear energy. One of its three goals is to target big challenges and game-changing technologies, to realize the mission and vision of AMMT program. Based on this context, this multi-year work package focuses on understanding the current state of large-scale additive manufacturing (AM) technology for the deposition of 316L stainless steel materials for final components and mild steel for use in nuclear manufacturing processes. The targeted AM modality is directed energy deposition (DED), capable of fabricating components on the size scale of meters including valves, pumps, impellers. etc. that are challenging or difficult to source, especially when developing new systems or replacing obsolete components. Accordingly, the current writeup aims at providing a baseline literature survey on structure-property relationships in mild steel and 316L alloys. Also, information on preliminary trials to date involving these two alloys show tremendous potential of printing parts having complex geometry and thin- walled structures, such as nuclear valve and Hot isostatic Press (HIP) can, using wire based (Wire Arc Additive Manufacturing and Hybrid Additive Manufacturing) as well as blown powder DED machines. All of this is aimed towards (i) demonstrating the ability to fabricate large components for pressure boundary applications relevant to the nuclear community and nuclear manufacturing technology, and (ii) understanding the effect of different manufacturing technology on AM can production and post HIPed material for nuclear applications. Another target of this writeup is to compile various in-situ monitoring tools that have been incorporated for different DED AM modalities, in order to understand process variability during the entire fabrication process. This can be correlated with processing-structure-property response surfaces and would add confidence around process quality verification and ultimately component certification for nuclear applications.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Earth as an Exoplanet. I. Time Variable Thermal Emission Using Spatially Resolved Moderate Imaging Spectroradiometer Data

Among the more than 4000 exoplanets known today, some terrestrial planets have been detected in the so-called habitable zone of their host stars and their number is expected to increase in the near future, energizing a drive to understand and interpret the eagerly awaited wealth of data to identify signs of life beyond our solar system. So far, Earth remains the best and only example of a habitable (and inhabited) world. Although, it seems extremely unlikely that any other exoplanets will be true Earth twins, it is important to explore and understand the full range of spectral signatures and variability of Earth in order to inform the design of future instruments and missions, and understand their diagnostic power as well as potential limitations. In this work we use Earth observation data collected by the MODIS instrument aboard the Aqua satellite. The complete data set comprises 15 years of thermal emission observations in the 3.66–14.40 μm range for five different locations on Earth (Amazon Rainforest, Antarctica, Arctic, Indian Ocean, and the Sahara Desert). We then determine flux levels and variations as a function of wavelength and surface type (i.e., climate zone and surface thermal properties) and investigate whether periodic signals indicating Earth’s tilted rotation axis can be detected. Our findings suggest that (1) viewing geometry plays an important role when thermal emission data is analyzed as Earth’s spectrum varies by a factor of three and more depending on the dominant surface type underneath; (2) typically strong absorption bands from CO 2 (15 μm) and O 3 (9.65 μm) are significantly less pronounced and partially absent in data from the polar regions implying that estimating correct abundance levels for these molecules might be challenging in these cases; and (3) the time-resolved thermal emission spectrum encodes information about seasons/planetary obliquity, but the significance depends on the viewing geometry and spectral band considered.

79 ASTRONOMY AND ASTROPHYSICS↗