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 343 records · Page 19

Nonlinear dynamics of quadratic gravity in spherical symmetry

We present the first numerically stable nonlinear evolution for the leading-order gravitational effective field theory (quadratic gravity) in the spherically-symmetric sector. The formulation relies on (i) the harmonic gauge to cast the evolution system into quasilinear form (ii) the Cartoon method to reduce to spherical symmetry in keeping with the harmonic gauge, and (iii) order reduction to first order (in time) by means of introducing auxiliary variables. The well posedness of the respective initial-value problem is numerically confirmed by evolving randomly perturbed flat-space and black-hole initial data. In conclusion, our study serves as a proof-of-principle for the possibility of stable numerical evolution in the presence of higher derivatives.

79 ASTRONOMY AND ASTROPHYSICS↗

JGCRI/hydro-emulator

The hydrological emulator (HE) was built to mimic complex global hydrological models (GHMs) at a range of spatial scales. The HE was written in MATLAB and utilizes a genetic algorithm to calibrate the a, b, c, d, m parameters from the ABCD runoff algorithm to observed runoff from GHMs when forced by the same climate data. The HE employs two methods to choose from when processing basin data: lumped (by basin) and distributed (gridded). Both methods can be evaluated using the built-in Kling-Gupta efficiency as a measure of model performance when compared to the outputs from the more complex GHM. Testing shows that the HE is 7 orders of magnitude faster then the widely used Variable Infiltration Capacity (VIC) GHM.

Vernon, Chris↗

Uncertainty quantification for incident helium flux in plasma-exposed tungsten

In this work, the surface response of a tungsten plasma-facing component was simulated by a cluster-dynamics code, Xolotl, with a focus on quantifying the impact of uncertainty in one of the input parameters to Xolotl, namely, the incident helium flux. The simulated conditions involve a tungsten surface exposed to 100 eV helium ion implantations with a flux of either 4 × 10 22 or 4 × 10 25 He m -2 s -1 . Furthermore, two sources were used to describe the implanted helium depth distribution in tungsten, either molecular dynamics (MD) or a binary collision approximation code, the stopping and range of ions in matter (SRIM). The aim of this work is to evaluate and examine uncertain predictions on the helium retention based on these two different modeling methodologies that either neglect electronic energy loss or the crystalline structure of the solid, respectively. An embedded model-form error approach was pursued here in order to arrive at predictions that account for variability due to the two different data sources, and the impact of this model-form uncertainty in incident helium flux on Xolotl output was presented for the two implantation fluxes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Parallel scaling analysis for explicit solid dynamics in ALEGRA

Weak scaling studies were performed for the explicit solid dynamics component of the ALEGRA code on two Cray supercomputer platforms during the period 2012-2015, involving a production-oriented hypervelocity impact problem. Results from these studies are presented, with analysis of the performance, scaling, and throughput of the code on these machines. The analysis demonstrates logarithmic scaling of the average CPU time per cycle up to core counts on the order of 10,000. At higher core counts, variable performance is observed, with significant upward excursions in compute time from the logarithmic trend. However, for core counts less than 10,000, the results show a 3 × improvement in simulation throughput, and a 2 × improvement in logarithmic scaling. This improvement is linked to improved memory performance on the Cray platforms, and to significant improvements made over this period to the data layout used by ALEGRA.

97 MATHEMATICS AND COMPUTING↗

Results of the Photometric LSST Astronomical Time-series Classification Challenge (PLAsTiCC)

Abstract Next-generation surveys like the Legacy Survey of Space and Time (LSST) on the Vera C. Rubin Observatory (Rubin) will generate orders of magnitude more discoveries of transients and variable stars than previous surveys. To prepare for this data deluge, we developed the Photometric LSST Astronomical Time-series Classification Challenge (PLAsTiCC), a competition that aimed to catalyze the development of robust classifiers under LSST-like conditions of a nonrepresentative training set for a large photometric test set of imbalanced classes. Over 1000 teams participated in PLAsTiCC, which was hosted in the Kaggle data science competition platform between 2018 September 28 and 2018 December 17, ultimately identifying three winners in 2019 February. Participants produced classifiers employing a diverse set of machine-learning techniques including hybrid combinations and ensemble averages of a range of approaches, among them boosted decision trees, neural networks, and multilayer perceptrons. The strong performance of the top three classifiers on Type Ia supernovae and kilonovae represent a major improvement over the current state of the art within astronomy. This paper summarizes the most promising methods and evaluates their results in detail, highlighting future directions both for classifier development and simulation needs for a next-generation PLAsTiCC data set.

79 ASTRONOMY AND ASTROPHYSICS↗

Coupled Neutron and Gamma Heating Calculation Based on VARIANT Nodal Transport Method

Accurate determination of temperature distributions in fuel pins and duct walls is essential for fuel-pin performance and assembly bowing analyses. Currently, the multi-assembly, sub-channel analysis code SE2-ANL is used in thermalhydraulic analysis of sodium cooled fast reactors (SFR). SE2- ANL is interfaced with the DIF3D/GAMSOR/DIF3D procedure for coupled neutron and gamma heating calculations, developed at Argonne National Laboratory (ANL). In this procedure, the neutron flux distribution is first calculated using the finite-difference diffusion theory option of the DIF3D code. Using the neutron flux solution, the GAMSOR code computes the gamma source distribution for the subsequent gamma flux calculation with DIF3D. However, in the current heating calculation scheme, the neutron and gamma flux calculations are limited to six triangular meshes per hexagonal assembly. In addition, the pin power distributions are estimated by assuming that the power distribution within an assembly is separable in the radial and axial directions. In each assembly, the axial power profile is obtained by averaging the power densities of six triangular meshes at each axial mesh interval. The radial power profile is evaluated at a single axial mesh (e.g., at the core mid-plane) by approximating the neutron and gamma flux distributions by a linear or quadratic polynomial of two variables x and y, depending on the polynomial order specified by the user. The coefficients of the linear and quadratic polynomials are determined by the least squares fit of six and 24 triangular mesh fluxes, respectively. Furthermore, the total heat generation rate in an assembly duct wall is simply determined by a user-specified fraction of the total gamma heating within that assembly. With the advancement of computer technology, the nodal transport theory code VARIANT is routinely used in core design analyses of SFR instead of DIF3D. It is therefore necessary to develop a new heating calculation procedure based on VARIANT neutron and gamma transport calculations. In this study, a new, coupled neutron and gamma heating calculation scheme has been developed to determine the pin power distributions more accurately by eliminating the aforementioned limitations of the current procedure. For this, the GAMSOR code has been extended to generate the intra-nodal gamma source distribution to be consistent with the trial functions of VARIANT. In the new heating calculation procedure, the neutron flux distribution is determined first by solving an eigenvalue transport problem using VARIANT. The obtained neutron flux is then used in GAMSOR to calculate the gamma source distribution. Using this gamma source distribution, the gamma flux distribution is determined by solving a fixed-source transport problem using VARIANT again. Using the computed neutron and gamma flux distributions, the power distribution in each fuel pin and duct wall is determined using a new utility code CURVE, which has been developed to evaluate the VARIANT flux solution at each fuel pin position and along the duct mid-wall. For each pin segment contained in a VARIANT node, the axial power shape is represented as a quadratic profile. The performance of the new heating calculation scheme based on VARIANT transport solutions has been tested against the existing method of SE2-ANL, using the ABR-1000 metal core design. The gamma source distribution and the assembly, fuel pin and duct wall power distributions were calculated and compared. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Are Baby Boomers' Non-Work Trip-Making Behavior Different than Millennials? Lessons Learned from NHTS Data

This paper presents a comparison between Millennials' and Baby Boomers' non-work travel behaviors using data from the 2017 National Household Travel Survey. Bootstrapped segmented ordered logit models are employed to capture the variability in travel preferences and trip frequency across these generational groups, providing more robust insights into their non-work travel. Millennials, particularly those who work from home, are found to have a negative association with higher non-work trip frequency, whereas Baby Boomers have a positive association with higher non-work trip frequency. The model results show that female Millennials who are heads of households are more likely to make non-work trips but less likely when living in urban areas. Ride sharing among Baby Boomers shows a higher association with non-work travel compared to Millennials. These insights could have implications for travel demand management, as shifting travel patterns necessitate adjustments in infrastructure investments and management strategies to support effective long-term transportation planning.

Patwary, Latif [ORNL] (ORCID:0000000189174928)↗

Navigation and Meteorological Data from Multiple Sensors on Airborne Platform (NAVMET-AIR) Value-Added Product Report

The Navigation and Meteorological Data from Multiple Sensors on Airborne Platform (NAVMET-AIR) value-added product encompasses the aafnaviwg data set in the ARM Data Center (ADC), known as the IWG file, named for the Inter-Agency Working Group for Airborne Data and Telemetry Systems (IWGADTS). Its purpose is to produce a suite of tools to promote standardization of instrument interface, data format, and data processing. This data set’s contents are navigational and meteorological state variables at 1 Hz. It also contains higher-order data flags created by ARM Aerial Facility (AAF) scientists to aid analysis such as periods when the aircraft is flying level, operating in maneuvers, or flying through cloud. Due to the nature of research flights, the payload of the aircraft consists of many duplicate and overlapping instruments to ensure, in case of instrument failure, there is a backup to record the data. Measurements from multiple instruments are consolidated into a single file. Each variable of this data set is carefully chosen from the onboard instrumentation and quality checked by AAF scientists to create this wholistic and most accurate data set for airborne research. This document is intended to be a hub towards the individual “read me” files produced for each campaign. For reference about what instrument was used for calculations on specific days, the read me files are located in the ARM IOP archive at the addresses listed, or access the IOP (intensive operational period) database at https://www.archive.arm.gov/.

54 ENVIRONMENTAL SCIENCES↗

Extreme Radio Flares and Associated X-Ray Variability from Young Stellar Objects in the Orion Nebula Cluster

Young stellar objects are known to exhibit strong radio variability on timescales of weeks to months, and a few reports have documented extreme radio flares with at least an order of magnitude change in flux density on timescales of hours to days. However, there have been few constraints on the occurrence rate of such radio flares or on the correlation with pre-main sequence X-ray flares, although such correlations are known for the Sun and nearby active stars. Here we report simultaneous deep VLA radio and Chandra X-ray observations of the Orion Nebula Cluster, targeting hundreds of sources to look for the occurrence rate of extreme radio variability and potential correlation with the most extreme X-ray variability. We identify 13 radio sources with extreme radio variability, with some showing an order of magnitude change in flux density in less than 30 minutes. All of these sources show X-ray emission and variability, but we find clear correlations with extreme radio flaring only on timescales <1 hr. Strong X-ray variability does not predict the extreme radio sources and vice versa. Radio flares thus provide us with a new perspective on high-energy processes in YSOs and the irradiation of their protoplanetary disks. Finally, our results highlight implications for interferometric imaging of sources violating the constant-sky assumption.

79 ASTRONOMY AND ASTROPHYSICS↗

Impact of Weak Agostic Interactions in Nickel Electrocatalysts for Hydrogen Oxidation

To understand how H2 binding and oxidation is influenced by [Ni(PR2NR'2)2]2+ PR2NR'2 catalysts with H2 binding energies close to thermoneutral, two [Ni(PPh2NR'2)2]2+ (R = Me or C14H29) complexes with phenyl substituents on phosphorous and varying alkyl chain lengths on the pendant amine were studied. In the solid state, [Ni(PPh2NMe2)2]2+ exhibits an anagostic interaction between the Ni(II) center and the a-CH3 of the pendant amine, and DFT and variable-temperature 31P NMR experiments suggest than the anagostic interaction persists in solution. The equilibrium constants for H2 addition to these complexes was measured by 31P NMR spectroscopy, affording free energies of H2 addition (?G°H2) of –0.8 kcal mol–1 in benzonitrile and –1.6 to –2.3 kcal mol–1 in THF. The anagostic interaction contributes to the low driving force for H2 binding by stabilizing the four-coordinate Ni(II) species prior to binding of H2. The pseudo-first order rate constants for H2 addition at 1 atm were measured by variable scan rate cyclic voltammetry, and were found to be similar for both complexes, less than 0.2 s–1 in benzonitrile and 3 –6 s–1 in THF. In the presence of exogenous base and H2 , turnover frequencies of electrocatalytic H2 oxidation were measured to be less than 0.2 s–1 in benzonitrile and 4 –9 s–1 in THF. These complexes are slower electrocatalysts for H2 oxidation than previously studied [Ni(PR2NR'2)2]2+ complexes due to a competition between H2 binding and formation of the anagostic interaction. However, the decrease in catalytic rate is accompanied by a beneficial 130 mV decrease in overpotential. This research was supported as part of the Center for Molecular Electrocatalysis, an Energy Frontier Research Center funded by the U.S. Department of Energy (DOE), Office of Science, Office of Basic Energy Sciences. Computational resources were provided at the National Energy Research Scientific Computing Center (NERSC) at Lawrence Berkeley National Laboratory. Mass spectrometry experiments were performed in the William R. Wiley Environmental Molecular Sciences Laboratory, a DOE national scientific user facility sponsored by the DOE’s Office of Biological and Environmental Research and located and the Pacific Northwest National Laboratory (PNNL). The authors thank Dr. Rosalie Chu for mass spectroscopy analysis. PNNL is operated by Battelle for DOE.

Klug, Christina↗

Evaluating fine-resolution, regional outputs of a variable resolution global climate model

Climate models have been used to study water resources and regional hydrologic responses to climate change, but climate model outputs must be downscaled to provide relevant regional data. However, the accuracy of this regional data is limited by uncertainties across and within downscaling methods, uncertainty across global outputs, and discontinuities at downscaled boundaries. A new alternative to traditional downscaling is a variable resolution model that incorporates fine-resolution regions directly into a coarse-resolution, global climate simulation in order to capture contiguous dynamics across resolution boundaries. In this study, we used the Variable-Resolution Community Earth System Model (VR-CESM) to generate one-eighth degree (14 km) fine-resolution outputs for the western U.S. and eastern China from 1970-2006. We focus our evaluation on precipitaiton, temperature, snow pack, solar radiation, and wind. We compare the model outputs with remote-sensing-based precipitation data, and both reanalysis and gridded weather station data for precipitation and temperature. VR-CESM generally has a cold bias in winter and a warm bias in summer in the western U.S., which compensate each other to reduce the annual bias. In eastern China, however, the sign of temperature biases are more consistent throughout the year with cold biases in the higher mountains and warm biases throughout most of the rest of the region. Precipitation biases are dependent upon reference data, and show slight overestimation in high mountain regions in both the U.S. and China with respect to gridded weather station data. Simulated snow cover in the western U.S. is reasonable compared to remote sensing data, but snow cover and snow water equivalent have larger biases when compared to reanalysis data. In eastern China there are widespread snow cover biases compared to remote sensing data. VR-CESM underestimates downward shortwave radiation to a greater degree in summer than in winter, and underestimates surface layer windspeed over mountains to a greater degree than in other areas. Comparison between VR-CESM and a coarser simulation (1-degree Beijing Climate Center model) shows reduced precipitation biases in the mountainous regions with finer resolution, indicating the value of variable-resolution modeling for reigonal studies.

Di Vittorio, Alan↗

Variable Transmission Line Ratings

This paper outlines the concept of variable-transmission facility ratings. This topic is of heightened importance following the recent Federal Energy Regulatory Commission (FERC) Order No. 881 which requires the implementation of ambient adjusted ratings (AARs) for FERC jurisdictional portions of the transmission system. Additionally, FERC issued an Advanced Notice of Proposed Rulemaking (ANOPR) implying a future order on Dynamic Line Rating (DLR) requirements. The various types of facility-rating methodologies will be discussed with a focus on how they can benefit the stability, capacity, and reliability of the grid.

24 - POWER TRANSMISSION AND DISTRIBUTION↗

Ship-based measurements of ice nuclei concentrations over the Arctic, Atlantic, Pacific and Southern oceans

Ambient concentrations of ice-forming particles measured during ship expeditions are collected and summarised with the aim of determining the spatial distribution and variability in ice nuclei in oceanic regions. The presented data from literature and previously unpublished data from over 23 months of ship-based measurements stretch from the Arctic to the Southern Ocean and include a circumnavigation of Antarctica. In comparison to continental observations, ship-based measurements of ambient ice nuclei show 1 to 2 orders of magnitude lower mean concentrations. To quantify the geographical variability in oceanic areas, the concentration range of potential ice nuclei in different climate zones is analysed by meridionally dividing the expedition tracks into tropical, temperate and polar climate zones. We find that concentrations of ice nuclei in these meridional zones follow temperature spectra with similar slopes but vary in absolute concentration. Typically, the frequency with which specific concentrations of ice nuclei are observed at a certain temperature follows a log-normal distribution. A consequence of the log-normal distribution is that the mean concentration is higher than the most frequently measured concentration. Finally, the potential contribution of ship exhaust to the measured ice nuclei concentration on board research vessels is analysed as function of temperature. We find a sharp onset of the influence at approximately -36 °C but none at warmer temperatures that could bias ship-based measurements.

54 ENVIRONMENTAL SCIENCES↗

Second-Order Uncertainty Analysis for Radiation Shielding Problems Using MCNP

Standard uncertainty propagation is based on a first-order Taylor expansion. This approximation is usually justified by assuming that 'uncertainties are very small numbers compared with the values of their corresponding quantities'. Furthermore, the second derivatives needed for a second-order analysis are not usually available, and a second-order analysis requires more knowledge regarding the distribution of the random variables. For radiation shielding problems, the MCNP6 Monte Carlo code computes first and second derivatives of a response with respect to nuclide density or reaction cross section. This paper explores the application of these derivatives to second-order uncertainty analysis.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A family of independent Variable Eddington Factor methods with efficient preconditioned iterative solvers

We present a family of discretizations for the Variable Eddington Factor (VEF) equations that have high-order accuracy on curved meshes and efficient preconditioned iterative solvers. The VEF discretizations are combined with the Discontinuous Galerkin transport discretization from to form effective high-order, linear transport methods. The VEF discretizations are derived by extending the unified analysis of Discontinuous Galerkin methods for elliptic problems presented by Arnold et al. to the VEF equations. This framework is used to define analogs of the interior penalty, second method of Bassi and Rebay, minimal dissipation local Discontinuous Galerkin, and continuous finite element methods. The analysis of subspace correction preconditioners, which use a continuous operator to iteratively precondition the discontinuous discretization, is extended to the case of the non-symmetric VEF system. Numerical results demonstrate that the VEF discretizations have arbitrary-order accuracy on curved meshes, preserve the thick diffusion limit, and are effective on a proxy problem from thermal radiative transfer in both outer transport iterations and inner preconditioned linear solver iterations. We demonstrate that the VEF solution converges to the S N transport solution as the mesh is refined on both problems with smooth and non-smooth behavior in angle. Parallel performance studies show that the interior penalty VEF discretization's linear solve weak scales out to 1024 processors and strong scales well on a single node. Particular attention is paid to the parallel performance of the VEF algorithm when used in combination with a parallel block Jacobi transport sweep.

97 MATHEMATICS AND COMPUTING↗

Transient Relativistic Fluid Dynamics in a General Hydrodynamic Frame

In this work, we propose a new theory of second-order viscous relativistic hydrodynamics which does not impose any frame conditions on the choice of the hydrodynamic variables. It differs from Mueller-Israel-Stewart theory by including additional transient degrees of freedom, and its first-order truncation reduces to Bemfica-Disconzi-Noronha-Kovtun theory. Conditions for causality and stability are explicitly given in the conformal regime. As an illustrative example, we consider Bjorken flow solutions to our equations and identify variables which make a hydrodynamic attractor manifest.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

$t\overline{t}b\overline{b}$ at NLO precision in a variable flavor number scheme

Top-quark pair production in association with two $b$-jets is computed at next-to-leading order QCD precision, including effects of the $b$-quark mass, and matched to a $t\overline{t}$+jets simulation in a variable flavor number scheme. The Monte Carlo realization of this method, called fusing, consistently embeds the four-flavor calculation in a particle-level event generator. As a first phenomenological application, we present observables relevant to the data-driven estimation of irreducible backgrounds to $t\overline{t}H$ -production.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Harnessing the Second-Order Metal−Insulator Transition for Neuromorphic Computing

Vanadium oxides are widely studied phase change materials for brain-inspired computing architectures. Systems like VO 2 and V 2 O 3 exhibit first-order metal−insulator transitions (MITs) with hysteresis and percolative switching, increasing stochasticity and device variability. Here, we focus on the less explored Magnéli phase V 4 O 7 , which undergoes a continuous, non-hysteretic, second-order MIT. This surprisingly enables highly reproducible volatile resistive switching in spiking-neuron-type devices. We synthesize V 4 O 7 films, characterize their structural and transport properties, and demonstrate voltage and current-driven threshold switching with electrothermal feedback. In a Pearson–Anson oscillator, V 4 O 7 devices produce stable, tunable spiking across 20–200 kHz, with consistent operation among multiple devices. We introduce a numerical analog leaky-integrate-and-fire (aLIF) model that captures waveform shapes and their dependence on load resistance, temperature, and voltage. Furthermore, these findings suggest that second-order MIT materials like V 4 O 7 are promising for deterministic, scalable spiking neuron arrays for neuromorphic computing.

V4O7↗