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 811 records · Page 45

The Conforming Virtual Element Method for the convection-diffusion-reaction equation with variable coeffcients.

This document describes the conforming formulations for virtual element approximation of the convection-reaction-diffusion equation with variable coefficients. Emphasis is given to construction of the projection operators onto polynomial spaces of appropriate order. These projections make it possible the virtual formulation to achieve any order of accuracy. We present the construction of the internal and the external formulation. The difference between the two is in the way the projection operators act on the derivatives (laplacian, gradient) of the partial differential equation. For the diffusive regime we prove the well-posedness of the external formulation and we derive an estimate of the approximation error in the H 1 -norm. For the convection-dominated case, the streamline diffusion stabilization (aka SUPG) is also discussed.

97 MATHEMATICS AND COMPUTING↗

Embedding Equality Constraints of Optimization Problems into a Quantum Annealer

Quantum annealers such as D-Wave machines are designed to propose solutions for quadratic unconstrained binary optimization (QUBO) problems by mapping them onto the quantum processing unit, which tries to find a solution by measuring the parameters of a minimum-energy state of the quantum system. While many NP-hard problems can be easily formulated as binary quadratic optimization problems, such formulations almost always contain one or more constraints, which are not allowed in a QUBO. Embedding such constraints as quadratic penalties is the standard approach for addressing this issue, but it has drawbacks such as the introduction of large coefficients and using too many additional qubits. In this paper, we propose an alternative approach for implementing constraints based on a combinatorial design and solving mixed-integer linear programming (MILP) problems in order to find better embeddings of constraints of the type Σ x i = k for binary variables x i . Our approach is scalable to any number of variables and uses a linear number of ancillary variables for a fixed k.

97 MATHEMATICS AND COMPUTING↗

Near-infrared Spectral Characterization of Solar-type Stars in the Northern Hemisphere

Although solar-analog stars have been studied extensively over the past few decades, most of these studies have focused on visible wavelengths, especially those identifying solar-analog stars to be used as calibration tools for observations. As a result, there is a dearth of well-characterized solar analogs for observations in the near-infrared, a wavelength range important for studying solar system objects. We present 184 stars selected based on solar-like spectral type and V−J and V−K colors whose spectra we have observed in the 0.8–4.2 μm range for calibrating our asteroid observations. Each star has been classified into one of three ranks based on spectral resemblance to vetted solar analogs. Of our set of 184 stars, we report 145 as reliable solar-analog stars, 21 as solar analogs usable after spectral corrections with low-order polynomial fitting, and 18 as unsuitable for use as calibration standards owing to spectral shape, variability, or features at low to medium resolution. We conclude that all but five of our candidates are reliable solar analogs in the longer wavelength range from 2.5 to 4.2 μm. The average colors of the stars classified as reliable or usable solar analogs are V−J = 1.148, V−H = 1.418, and V−K = 1.491, with the entire set being distributed fairly uniformly in R.A. across the sky between −27° and +67° in decl.

79 ASTRONOMY AND ASTROPHYSICS↗

A multitree approach for global solution of ACOPF problems using piecewise outer approximations

Electricity markets rely on the rapid solution of the optimal power flow (OPF) problem to determine generator power levels and set nodal prices. Traditionally, the OPF problem has been formulated using linearized, approximate models, ignoring nonlinear alternating current (AC) physics. These approaches do not guarantee global optimality or even feasibility in the real ACOPF problem. We introduce an outer-approximation approach to solve the ACOPF problem to global optimality based on alternating solution of upper- and lower-bounding problems. The lower-bounding problem is a piecewise relaxation based on strong second-order cone relaxations of the ACOPF, and these piecewise relaxations are selectively refined at each major iteration through increased variable domain partitioning. Our approach is able to efficiently solve all but one of the test cases considered to an optimality gap below 0.1%. This approach opens the door for global solution of MINLP problems with AC power flow equations.

42 ENGINEERING↗

Overview of the SCEC/USGS Community Stress Drop Validation Study Using the 2019 Ridgecrest Earthquake Sequence

We present initial findings from the ongoing Community Stress Drop Validation Study to compare spectral stress-drop estimates for earthquakes in the 2019 Ridgecrest, California, sequence. This study uses a unified dataset to independently estimate earthquake source parameters through various methods. Stress drop, which denotes the change in average shear stress along a fault during earthquake rupture, is a critical parameter in earthquake science, impacting ground motion, rupture simulation, and source physics. Spectral stress drop is commonly derived by fitting the amplitude-spectrum shape, but estimates can vary substantially across studies for individual earthquakes. Sponsored jointly by the U.S. Geological Survey and the Statewide (previously, Southern) California Earthquake Center our community study aims to elucidate sources of variability and uncertainty in earthquake spectral stress-drop estimates through quantitative comparison of submitted results from independent analyses. The dataset includes nearly 13,000 earthquakes ranging from M 1 to 7 during a two-week period of the 2019 Ridgecrest sequence, recorded within a 1° radius. Here, in this article, we report on 56 unique submissions received from 20 different groups, detailing spectral corner frequencies (or source durations), moment magnitudes, and estimated spectral stress drops. Methods employed encompass spectral ratio analysis, spectral decomposition and inversion, finite-fault modeling, ground-motion-based approaches, and combined methods. Initial analysis reveals significant scatter across submitted spectral stress drops spanning over six orders of magnitude. However, we can identify between-method trends and offsets within the data to mitigate this variability. Averaging submissions for a prioritized subset of 56 events shows reduced variability of spectral stress drop, indicating overall consistency in recovered spectral stress-drop values.

58 GEOSCIENCES↗

Investigation of pulsed mode operation with the frequency tuned CAPRICE ECRIS

In order to increase the intensity of the highly charged ions produced by the Electron Cyclotron Resonance Ion Sources (ECRISs), techniques like the frequency tuning and the afterglow mode have been developed and in this paper the effect on the ion production is shown for the first time when combining both techniques. Recent experimental results proved that the tuning of the operating frequency of the ECRIS is a promising technique to achieve higher ion currents of higher charge states. On the other hand, it is well known that the afterglow mode of the ECRIS operation can provide more intense pulsed ion beams in comparison with the continuous wave (cw) operation. These two techniques can be combined by pulsing the variable frequency signal driving the traveling wave tube amplifier which provides the high microwave power to the ECRIS. In order to analyze the effect of these two combined techniques on the ion source performance, several experiments were carried out on the pulsed frequency tuned CAPRICE (Compacte source A Plusiers Résonances Ionisantes Cyclotron Electroniques)-type ECRIS. Different waveforms and pulse lengths have been investigated under different settings of the ion source. The results of the pulsed mode have been compared with those of cw operation.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Predicting carbon dioxide and energy fluxes across global FLUXNET sites with regression algorithms

Spatio-temporal fields of land–atmosphere fluxes derived from data-driven models can complement simulations by process-based land surface models. While a number of strategies for empirical models with eddy-covariance flux data have been applied, a systematic intercomparison of these methods has been missing so far. In this study, we performed a cross-validation experiment for predicting carbon dioxide, latent heat, sensible heat and net radiation fluxes across different ecosystem types with 11 machine learning (ML) methods from four different classes (kernel methods, neural networks, tree methods, and regression splines). We applied two complementary setups: (1) 8-day average fluxes based on remotely sensed data and (2) daily mean fluxes based on meteorological data and a mean seasonal cycle of remotely sensed variables. The patterns of predictions from different ML and experimental setups were highly consistent. There were systematic differences in performance among the fluxes, with the following ascending order: net ecosystem exchange ( R 2 < 0.5), ecosystem respiration ( R 2 > 0.6), gross primary production ( R 2 > 0.7), latent heat ( R 2 > 0.7), sensible heat ( R 2 > 0.7), and net radiation ( R 2 > 0.8). The ML methods predicted the across-site variability and the mean seasonal cycle of the observed fluxes very well ( R 2 > 0.7), while the 8-day deviations from the mean seasonal cycle were not well predicted ( R 2 < 0.5). Fluxes were better predicted at forested and temperate climate sites than at sites in extreme climates or less represented by training data (e.g., the tropics). Finally, the evaluated large ensemble of ML-based models will be the basis of new global flux products.

54 ENVIRONMENTAL SCIENCES↗

Effects of Cr substitution on structural and magnetic properties in La{sub 0.7}Pr{sub 0.3}Fe{sub 11.4}Si{sub 1.6} compound

In an effort to explore the effect of substitution Fe by Cr in NaZn{sub 13}-type La{sub 0.7}Pr{sub 0.3}Fe{sub 11.4−x}Cr{sub x}Si{sub 1.6} (x = 0, 0.06, 0.12, 0.26, and 0.34) compounds, the structure and magnetic properties have been investigated by high intensity of x-ray and neutron diffraction, scanning electron microscopy, specific heat, and magnetization measurement. It has been found that a substitution of Cr for Fe in this compounds leads to decrease in the lattice parameter a at room temperature but variation on Curie temperature (T{sub C}). While the first order nature of magnetic phase transition around T{sub C} does not change with increasing Cr content up to x = 0.34. High intensity x-ray and neutron diffraction study at variable temperatures for highest Cr concentration x = 0.34 confirmed the presence of strong magneto-volume effect around T{sub C} and indicated the direct evident of coexistence across magnetic transition as characteristic of first order nature. The values of −ΔS{sub M} around T{sub C} decrease from 17 J kg{sup −1} K{sup −1} for x = 0 to 12 J kg{sup −1} K{sup −1} for x = 0.06 and then increases with further increasing Cr content up to 17.5 J kg{sup −1} K{sup −1} for x = 0.34 under a change of 0–5 T magnetic field. Similar behavior on relative cooling power which is decrease from 390 J kg{sup −1} for x = 0 to 365 J kg{sup −1} for x = 0.06 at the beginning and then increases up to 400 J kg{sup −1} for x = 0.34. From the point of this view with the highest Cr concentration (x = 0.34) exhibits favourable material candidate for magnetic refrigerator application therefore should inspire further study concerning on higher Cr concentration in this compound.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

An updated Lagrangian discontinuous Galerkin hydrodynamic method for gas dynamics

Here, we present a new Lagrangian discontinuous Galerkin (DG) hydrodynamic method for gas dynamics. The new method evolves conserved unknowns in the current configuration, which obviates the Jacobi matrix that maps the element in a reference coordinate system or the initial coordinate system to the current configuration. The density, momentum, and total energy (ρ, ρu, E ) are approximated with conservative higher-order Taylor expansions over the element and are limited toward a piecewise constant field near discontinuities using a limiter. Two new limiting methods are presented for enforcing the bounds on the primitive variables of density, velocity, and specific internal energy (ρ, u, e). The nodal velocity, and the corresponding forces, are calculated by solving an approximate Riemann problem at the element nodes. An explicit second-order method is used to temporally advance the solution. This new Lagrangian DG hydrodynamic method conserves mass, momentum, and total energy. 1D Cartesian coordinates test problem results are presented to demonstrate the accuracy and convergence order of the new DG method with the new limiters.

97 MATHEMATICS AND COMPUTING↗

Assessment of NHTSA’s Report “Relationships Between Fatality Risk, Mass, and Footprint in Model Year 2004-2011 Passenger Cars and LTVs” (LBNL Phase 1)

In its 2012 report NHTSA simulated the effect four fleetwide mass reduction scenarios would have on the change in annual fatalities. NHTSA estimated that the most aggressive of these scenarios (reducing mass 5.2% in heavier light trucks and 2.6% in all other vehicles types except lighter cars) would result in a small reduction in societal fatalities. LBNL replicated the methodology NHTSA used to simulate six mass reduction scenarios, including the mass reductions recommended in the 2015 NRC committee report, and estimated in 2021 and 2025 by EPA in the TAR, using the updated data through 2012. The analysis indicates that the estimated x change in fatalities under each scenario based on the updated analysis is comparable to that in the 2012 analysis, but less beneficial or more detrimental than that in the 2016 analysis. For example, an across the board 100-lb reduction in mass would result in an estimated 157 additional annual fatalities based on the 2012 analysis, but would result in only an estimated 91 additional annual fatalities based on the 2016 analysis, and an additional 87 fatalities based on the current analysis. The mass reductions recommended by the 2015 NRC committee report would result in a 224 increase in annual fatalities in the 2012 analysis, a 344 decrease in annual fatalities in the 2016 analysis, and a 141 increase in fatalities in the current analysis. The mass reductions EPA estimated for 2025 in the TAR would result in a 203 decrease in fatalities based on the 2016 analysis, but an increase of 39 fatalities based on the current analysis. These results support NHTSA’s conclusion from its 2012 study that, when footprint is held fixed, “no judicious combination of mass reductions in the various classes of vehicles results in a statistically significant fatality increase and many potential combinations are safety-neutral as point estimates.”Like the previous NHTSA studies, this updated report concludes that the estimated effect of mass reduction while maintaining footprint on societal U.S. fatality risk is small, and not statistically significant at the 95% or 90% confidence level for all vehicle types based on the jack-knife method NHTSA used. This report also finds that the estimated effects of other control variables, such as vehicle type, specific safety technologies, and crash conditions such as whether the crash occurred at night, in a rural county, or on a high-speed road, on risk are much larger, in some cases two orders of magnitude larger, than the estimated effect of mass or footprint reduction on risk. Finally, this report shows that after accounting for the many vehicle, driver, and crash variables NHTSA used in its regression analyses, there remains a wide variation in risk by vehicle make and model, and this variation is unrelated to vehicle mass. Although the purpose of the NHTSA and LBNL reports is to estimate the effect of vehicle mass reduction on societal risk, this is not exactly what the regression models are estimating. Rather, they are estimating the recent historical relationship between mass and risk, after accounting for most measurable differences between vehicles, drivers, and crash times and locations. In essence, the regression models are comparing the risk of a 2600-lb Dodge Neon with that of a 2500-lb Honda Civic, after attempting to account for all other differences between the two vehicles. The models are not estimating the effect of literally removing 100 pounds from the Neon, leaving everything else unchanged. In addition, the analyses are based on the relationship of vehicle mass and footprint on risk for recent vehicle designs (model year 2004 to 2011). These relationships may or may not continue into the future as manufacturers utilize new vehicle designs and incorporate new technologies, such as more extensive use of strong lightweight materials and specific safety technologies. Therefore, throughout this report we use the phrase “the estimated effect of mass (or footprint) reduction on risk” as shorthand for “the estimated change in risk as a function of its relationship to mass (or footprint) for vehicle models of recent design.”

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Conformational space exploration of cryo-EM structures by variability refinement

Cryo-EM observation of biological samples enables visualization of sample heterogeneity, in the form of discrete states that are separable, or continuous heterogeneity as a result of local protein motion before flash freezing. Variability analysis of this continuous heterogeneity describes the variance between a particle stack and a volume, and results in a map series describing the various steps undertaken by the sample in the particle stack. While this observation is absolutely stunning, it is very hard to pinpoint structural details to elements of the maps. Here, in order to bridge the gap between observation and explanation, we designed a tool that refines an ensemble of structures into all the maps from variability analysis. Using this bundle of structures, it is easy to spot variable parts of the structure, as well as the parts that are not moving. Comparison with molecular dynamics simulations highlights the fact that the movements follow the same directions, albeit with different amplitudes. Ligand can also be investigated using this method. Variability refinement is available in the Phenix software suite, accessible under the program name phenix.varref.

59 BASIC BIOLOGICAL SCIENCES↗

Degenerate higher derivative theories beyond Horndeski: evading the Ostrogradski instability

Theories with higher order time derivatives generically suffer from ghost-like instabilities, known as Ostrogradski instabilities. This fate can be avoided by considering ''degenerate'' Lagrangians, whose kinetic matrix cannot be inverted, thus leading to constraints between canonical variables and a reduced number of physical degrees of freedom. In this work, we derive in a systematic way the degeneracy conditions for scalar-tensor theories that depend quadratically on second order derivatives of a scalar field. We thus obtain a classification of all degenerate theories within this class of scalar-tensor theories. The quartic Horndeski Lagrangian and its extension beyond Horndeski belong to these degenerate cases. We also identify new families of scalar-tensor theories with the property that they are degenerate despite the nondegeneracy of the purely scalar part of their Lagrangian.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Application of the Aquifer Impact Model to support decisions at a CO 2 sequestration site

The National Risk Assessment Partnership (NRAP) has developed a suite of tools to assess and manage risk at CO 2 sequestration sites. The NRAP tool suite includes the Aquifer Impact Model (AIM), which evaluates the potential for groundwater impacts from leaks of CO 2 and brine through abandoned wellbores. There are two aquifer reduced-order models (ROMs) included with the AIM tool, a confined alluvium aquifer, and an unconfined carbonate aquifer. The models accept aquifer parameters as a range of variable inputs so they may have broad applicability. The generic aquifer models may be used at the early stages of site selection, when site-specific data is not available. Guidelines have been developed for determining when the generic ROMs might be applicable to a new site. This work considers the application of the AIM to predicting the impact of CO 2 or brine leakage were it to occur at the Illinois Basin Decatur Project (IBDP). Results of the model sensitivity analysis can help guide characterization efforts; the hydraulic parameters and leakage source term magnitude are more sensitive than clay fraction or cation exchange capacity. Sand permeability was the only hydraulic parameter measured at the IBDP site. More information on the other hydraulic parameters could reduce uncertainty in risk estimates. Some non-adjustable parameters are significantly different for the ROM than for the observations at the IBDP site. The generic ROMs could be made more useful to a wider range of sites if the initial conditions and no-impact threshold values were adjustable parameters.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Tetrafluoroethane (R134a) hydrate formation within variable volume reactor accompanied by evaporation and condensation

Vast size hydrate formation reactors with fast conversion rate are required for the economic implementation of seawater desalination utilizing gas hydrate technology. The commercial target production rate is order of thousand tons of potable water per day per train. Various heat and mass transfer enhancement schemes including agitation, spraying, and bubbling have been examined to maximize the production capacities in scaled up design of hydrate formation reactors. The present experimental study focused on acquiring basic knowledge needed to design variable volume reactors to produce tetrafluoroethane hydrate slurry. Test vessel was composed of main cavity with fixed volume of 140 ml and auxiliary cavity with variable volume of 0 ∼ 64 ml. Temperatures at multiple locations within vessel and pressure were monitored while visual access was made through front window. Alternating evaporation and condensation induced by cyclic volume change provided agitation due to density differences among water and vapor, liquid and hydrate R134a as well as extended interface area, which improved hydrate formation kinetics coupled with latent heat release and absorption. Influences of coolant temperature, piston stroke/speed, and volume change period on hydrate formation kinetics were investigated. Suggestions of reactor design improvement for future experimental study are also made.

42 ENGINEERING↗

Deep drilling in the time domain with DECam – II: characterizing the light curves of candidates in the extragalactic fields

ABSTRACT In this second paper on the DECam deep-drilling field (DDF) program, we release 2020 optical $gri$-band light curves for transients and variables in the extragalactic COSMOS and ELAIS fields based on time series observations with a 3-d cadence from semester 2021A through 2023A. In order to demonstrate the wide variety of time domain events detected by the program and encourage others to use the data set, we characterize the sample by presenting a brief analysis of the light-curve parameters such as time span, amplitude, and peak brightness. We also present preliminary light-curve categorizations, and identify potential stellar variables, active galactic nuclei, tidal disruption events, supernovae (SNe) (such as Type Ia, Type IIP, superluminous, and gravitationally lensed SNe), and fast transients. Where relevant, the number of identified transients is compared to the predictions of the original proposal. We also discuss the challenges of analyzing DDF data in the context of the upcoming Vera C. Rubin Observatory and its Legacy Survey of Space and Time, which will include DDFs. Images from the Dark Energy Camera DDF program are available without proprietary period and the light curves presented in this work are publicly available for analysis.

Graham, Melissa L. (ORCID:0000000291543136)↗

Resilient Cyber Infrastructure for the Minimum Wind Curtailment Remedial Control Scheme

Integration of Distributed Energy Resources (DER) has been significantly increasing in the electric power grid. One of the predominant renewable energy sources, wind energy, possess a significant challenge for grid operations, due to its intermittency, uncertainty and variability. In changing operating conditions with wind uncertainty, part of wind may need to be curtailed for the safe utilization of the transmission system. In order to minimize the wind power curtailment and to keep the power flow on a transmission line under limits, a cyber resilient Distributed Remedial Action Scheme (DRAS) is developed in this paper. The RAS is formulated as an optimization problem, considering the real power flow constraints of the transmission lines and given wind variability. In this paper, DRAS is developed to operate even with failures in the computing nodes to make it cyber-resilient. Distributed simplex algorithm is used for optimization in the distributed RAS. The DRAS is implemented in a decentralized platform called Resilient Information Architecture Platform for Smart Grid (RIAPS) to enable cyber resiliency of DRAS. The effectiveness of the proposed approach is validated through offline simulations of the New England 39 bus system and online simulations of IEEE 14 bus test system. In conclusion, a cyber physical test bed utilizing Real Time Digital Simulator (RTDS), Phasor Measurement Units (PMUs), and BeagleBones has been used for the online simulation and validation of the developed cyber-resilient RAS.

17 WIND ENERGY↗

The Challenge of Identifying Controls on Cloud Properties and Precipitation Onset for Cumulus Congestus Sampled During MC3E

Here, the controls on precipitation onset and the transition from shallow cumulus to congestus are explored using a suite of 16 large–eddy simulations based on the 25 May 2011 event from the Midlatitude Continental Convective Clouds Experiment (MC3E). The thermodynamic variables in the model are relaxed at various timescales to observationally constrained temperature and moisture profiles in order to better reproduce the observed behavior of precipitation onset and total precipitation. Three of the simulations stand out as best matching the precipitation observations and also perform well for independent comparisons of cloud fraction, precipitation area fraction, and evolution of cloud top occurrence. All three simulations exhibit a destabilization over time, which leads to a transition to deeper clouds, but the evolution of traditional stability metrics by themselves is not able to explain differences in the simulations. Conditionally sampled cloud properties (in particular, mean cloud buoyancy), however, do elicit differences among the simulations. The inability of environmental profiles alone to discern subtle differences among the simulations and the usefulness of conditionally sampled model quantities argue for hybrid observational/modeling approaches. These combined approaches enable a more complete physical understanding of cloud systems by combining observational sampling of time–varying three–dimensional meteorological quantities and cloud properties, along with detailed representation of cloud microphysical and dynamical processes from numerical models.

54 ENVIRONMENTAL SCIENCES↗

Enhanced data validation strategy of air quality monitoring network

Quick validation and detection of faults in measured air quality data is a crucial step towards achieving the objectives of air quality networks. Therefore, the objectives of this paper are threefold: (i) to develop a modeling technique that can be used to predict the normal behavior of air quality variables and help provide accurate reference for monitoring purposes; (ii) to develop fault detection method that can effectively and quickly detect any anomalies in measured air quality data. For this purpose, a new fault detection method that is based on the combination of generalized likelihood ratio test (GLRT) and exponentially weighted moving average (EWMA) will be developed. GLRT is a well-known statistical fault detection method that relies on maximizing the detection probability for a given false alarm rate. In this paper, we propose to develop GLRT-based EWMA fault detection method that will be able to detect the changes in the values of certain air quality variables; (iii) to develop fault isolation and identification method that allows defining the fault source(s) in order to properly apply appropriate corrective actions. In this paper, reconstruction approach that is based on Midpoint-Radii Principal Component Analysis (MRPCA) model will be developed to handle the types of data and models associated with air quality monitoring networks. All air quality modeling, fault detection, fault isolation and reconstruction methods developed in this paper will be validated using real air quality data (such as particulate matter, ozone, nitrogen and carbon oxides measurement).

54 ENVIRONMENTAL SCIENCES↗