Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Polynomial models”

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 253 records · Page 14

Learning infinite-horizon average-reward restless multi-action bandits via index awareness

We consider the online restless bandits with average-reward and multiple actions, where the state of each arm evolves according to a Markov decision process (MDP), and the reward of pulling an arm depends on both the current state of the corresponding MDP and the action taken. Since finding the optimal control is typically intractable for restless bandits, existing learning algorithms are often computationally expensive or with a regret bound that is exponential in the number of arms and states. In this paper, we advocate \textit{index-aware reinforcement learning} (RL) solutions to design RL algorithms operating on a much smaller dimensional subspace by exploiting the inherent structure in restless bandits. Specifically, we first propose novel index policies to address dimensionality concerns, which are provably optimal. We then leverage the indices to develop two low-complexity index-aware RL algorithms, namely, (i) GM-R2MAB, which has access to a generative model; and (ii) UC-R2MAB, which learns the model using an upper confidence style online exploitation method. We prove that both algorithms achieve a sub-linear regret that is only polynomial in the number of arms and states. A key differentiator between our algorithms and existing ones stems from the fact that our RL algorithms contain a novel exploitation that leverages our proposed provably optimal index policies for decision-makings.

Xiong, Guojun↗

Parameterization of cloud transmittance for expeditious assessment and forecasting of all-sky DNI

Radiative transfer models require vast computing resources to solve cloud transmittance and reflectance from the radiative transfer equation. As a result, models offering precise simulation in operations often acquire individual cloud transmittance or reflectance from a lookup table precomputed for practicable scenarios. To further expedite the computation of global horizontal irradiance and to reduce the storage requirements, the Fast All-sky Radiation Model for Solar applications (FARMS) parameterized the lookup table using elementary functions with specified coefficients. This study extends FARMS direct normal irradiance (DNI) computation by utilizing hyperbolic tangent functions and various polynomial functions to parameterize the cloud transmittance for scattered solar radiation in the circumsolar region. The parameterization is implemented in FARMS with DNI (FARMS-DNI) and accounts for the circumsolar radiation when assessing or forecasting DNI. The evaluation, with long-term observations at the National Renewable Energy Laboratory's, Solar Radiation Research Laboratory, and the Atmospheric Radiation Measurement, Southern Great Plains, Central Facility, shows that the parameterized DNIs are virtually identical with those computed by coupling FARMS-DNI to a lookup table of cloud transmittance. This parameterization has diverse applications in radiative transfer models and numerical weather prediction models used to assess or forecast direct solar radiation.

14 SOLAR ENERGY↗

Pseudospectral particle-in-cell formulation with arbitrary charge and current-density time dependencies for the modeling of relativistic plasmas

This paper introduces a formulation of the particle-in-cell (PIC) method for the modeling of relativistic plasmas, that leverages the ability of the pseudospectral analytical time-domain solver (PSATD) to handle arbitrary time dependencies of the charge and current densities during one PIC cycle (applied to second-order polynomial dependencies here). Here, the formulation is applied to a modified set of Maxwell's equations that was proposed earlier in the context of divergence cleaning, and to recently proposed extensions of the PSATD-PIC algorithm. Detailed analysis and testings revealed that, under some condition, the formulation can expand the range of numerical parameters under which PIC simulations are stable and accurate when modeling relativistic plasmas such as, e.g., plasma-based particle accelerators.

43 PARTICLE ACCELERATORS↗

WRF-simulated low-level jets over Iowa: characterization and sensitivity studies

Abstract. Output from 6 months of high-resolution simulations with the Weather Research and Forecasting (WRF) model are analyzed to characterize local low-level jets (LLJs) over Iowa for winter and spring in the contemporary climate. Low-level jets affect rotor plane aerodynamic loading, turbine structural loading and turbine performance, and thus accurate characterization and identification are pertinent. Analyses using a detection algorithm wherein the wind speed above and below the jet maximum must be below 80 % of the jet wind speed within a vertical window of approximately 20–530 m a.g.l. (above ground level) indicate the presence of an LLJ in at least one of the 14 700 4 km×4 km grid cells over Iowa on 98 % of nights. Nocturnal LLJs are most frequently associated with stable stratification and low turbulent kinetic energy (TKE) and hence are more frequent during the winter months. The spatiotemporal mean LLJ maximum (jet core) wind speed is 9.55 m s−1, and the mean height is 182 m. Locations of high LLJ frequency and duration across the state are seasonally varying, with a mean duration of 3.5 h. The highest frequency occurs in the topographically complex northwest of the state in winter and in the flatter northeast of the state in spring. Sensitivity of LLJ characteristics to the (i) LLJ definition and (ii) vertical resolution at which the WRF output is sampled is examined. LLJ definitions commonly used in the literature are considered in the first sensitivity analysis. These sensitivity analyses indicate that LLJ characteristics are highly variable with definition. Use of different definitions identifies both different frequencies of LLJs and different LLJ events. Further, when the model output is down-sampled to lower vertical resolution, the mean jet core wind speed height decreases, but spatial distributions of regions of high frequency and duration are conserved. Implementation of a polynomial interpolation to extrapolate down-sampled output to full-resolution results in reduced sensitivity of LLJ characteristics to down-sampling.

17 WIND ENERGY↗

Axion domain walls, small instantons, and non-invertible symmetry breaking

Non-invertible global symmetry often predicts degeneracy in axion potentials and carries important information about the global form of the gauge group. When these symmetries are spontaneously broken they can lead to the formation of stable axion domain wall networks which support topological degrees of freedom on their worldvolume. Such non-invertible symmetries can be broken by embedding into appropriate larger UV gauge groups where small instanton contributions lift the vacuum degeneracy, and provide a possible solution to the domain wall problem. We explain these ideas in simple illustrative examples and then apply them to the Standard Model, whose gauge algebra and matter content are consistent with several possible global structures. Each possible global structure leads to different selection rules on the axion couplings, and various UV completions of the Standard Model lead to more specific relations. As a proof of principle, we also present an example of a UV embedding of the Standard Model which can solve the axion domain wall problem. The formation and annihilation of the long-lived axion domain walls can lead to observables, such as gravitational wave signals. Observing such signals, in combination with the axion coupling measurements, can provide valuable insight into the global structure of the Standard Model, as well as its UV completion.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

On amplitudes and field redefinitions

We derive an off-shell recursion relation for correlators that holds at all loop orders. This allows us to prove how generalized amplitudes transform under generic field redefinitions, starting from an assumed behavior of the one-particle-irreducible effective action. The form of the recursion relation resembles the operation of raising the rank of a tensor by acting with a covariant derivative. This inspires a geometric interpretation, whose features and flaws we investigate.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

tesuract v.1.0

SAND2021-14370 O Tesuract (tensor surrogate approximation and computation) provides Python tools to build machine learning regressors that include polynomial chaos expansions. The software also offers methods for creation regressors for tensor outputs such as spatially varying fields and using a reduced order modeling methodology. The library is built on and compatible with the scikit-learn API, which allows integration with one of the most ubiquitous machine learning Python libraries out there. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Chowdhary, Kamaljit↗

New Determinations of the UV Luminosity Functions from z ~ 9 to 2 Show a Remarkable Consistency with Halo Growth and a Constant Star Formation Efficiency

Here we provide the most comprehensive determinations of the rest-frame UV luminosity function (LF) available to date with the Hubble Space Telescope (HST) at z ~ 2–9. Essentially all of the noncluster extragalactic legacy fields are utilized, including the Hubble Ultra Deep Field, the Hubble Frontier Fields parallel fields, and all five CANDELS fields, for a total survey area of 1136 arcmin{sup 2}. Our determinations include galaxies at z ~ 2–3 leveraging the deep HDUV, UVUDF, and ERS WFC3/UVIS observations available over an ~150 arcmin{sup 2} area in the GOODS-North and GOODS-South regions. All together, our collective samples include >24,000 sources, >2.3× larger than previous selections with HST. We identify 5766, 6332, 7240, 3449, 1066, 601, 246, and 33 sources at z ~ 2, 3, 4, 5, 6, 7, 8, and 9, respectively. Combining our results with an earlier z ~ 10 LF determination by Oesch et al., we quantify the evolution of the UV LF. Our results indicate that there is (1) a smooth flattening of the faint-end slope α from α ~ -2.4 at z ~ 10 to α ~ -1.5 at z ~ 2, (2) minimal evolution in the characteristic luminosity M* at z ≥ 2.5, and (3) a monotonic increase in the normalization log{sub 10}ϕ* from z ~ 10 to 2, which can be well described by a simple second-order polynomial, consistent with an “accelerated” evolution scenario. We find that each of these trends (from z ~ 10 to 2.5 at least) can be readily explained on the basis of the evolution of the halo mass function and a simple constant star formation efficiency model.

47 OTHER INSTRUMENTATION↗

UQ Toolkit v 2.0

The Uncertainty Quantification (UQ) Toolkit is a software library for the characterizaton and propagation of uncertainties in computational models. For the characterization of uncertainties, Bayesian inference tools are provided to infer uncertain model parameters, as well as Bayesian compressive sensing methods for discovering sparse representations of high-dimensional input-output response surfaces, and also Karhunen-Loève expansions for representing stochastic processes. Uncertain parameters are treated as random variables and represented with Polynomial Chaos expansions (PCEs). The library implements several spectral basis function types (e.g. Hermite basis functions in terms of Gaussian random variables or Legendre basis functions in terms of uniform random variables) that can be used to represent random variables with PCEs. For propagation of uncertainty, tools are provided to propagate PCEs that describe the input uncertainty through the computational model using either intrusive methods (Galerkin projection of equations onto basis functions) or non-intrusive methods (perform deterministic operation at sampled values of the random values and project the obtained results onto basis functions).

Safta, Cosmin↗

A classical, non-singular, bouncing universe

We present a model for a classical, non-singular bouncing cosmology without violation of the null energy condition (NEC). The field content is General Relativity plus a real scalar field with a canonical kinetic term and only renormalizable, polynomial-type self-interactions for the scalar field in the Jordan frame. The universe begins vacuum-energy dominated and is contracting at t=-∞. We consider a closed universe with a positive spatial curvature, which is responsible for the universe bouncing without any NEC violation. An Rφ 2 coupling between the Ricci scalar and the scalar field drives the scalar field from the initial false vacuum to the true vacuum during the bounce. The model is sub-Planckian throughout its evolution and every dimensionful parameter is below the effective-field-theory scale M P , so we expect no ghost-type or tachyonic instabilities. This model solves the horizon problem and extends co-moving particle geodesics to past infinity, resulting in a geodesically complete universe without singularities. We solve the Friedman equations and the scalar-field equation of motion numerically, and analytically under certain approximations.

Astronomy & Astrophysics↗

A Generalized Copula-Polynomial Chaos Expansion for Probabilistic Power Flow Considering Nonlinear Correlations of PV Injections

This paper develops a generalized Copula-polynomial chaos expansion (PCE) framework for power system probabilistic power flow that can handle both linear and nonlinear correlations of uncertain power injections, such as wind and PVs. A data-driven Copula statistical model is used to capture the correlations of uncertain power injections. This allows us to resort to the Rosenblatt transformation to transform correlated variables into independent ones while preserving the dependence structure. This paves the way of leveraging the PCE for surrogate modeling and uncertainty quantification of power flow results, i.e., achieving the probabilistic distributions of power flows. Simulations carried out on the IEEE 57-bus system show that the proposed framework can get much more accurate results than other alternatives with different linear and nonlinear power injection correlations.

copula↗

Feature selection with distance correlation

Choosing which properties of the data to use as input to multivariate decision algorithms—also known as feature selection—is an important step in solving any problem with machine learning. While there is a clear trend towards training sophisticated deep networks on large numbers of relatively unprocessed inputs (so-called automated feature engineering), for many tasks in physics, sets of theoretically well-motivated and well-understood features already exist. Working with such features can bring many benefits, including greater interpretability, reduced training and run time, and enhanced stability and robustness. We develop a new feature selection method based on distance correlation, and demonstrate its effectiveness on the tasks of boosted top- and W -tagging. Using our method to select features from a set of over 7,000 energy flow polynomials, we show that we can match the performance of much deeper architectures, by using only ten features and two orders-of-magnitude fewer model parameters. Published by the American Physical Society 2024

Astronomy & Astrophysics↗

High-Order Wall-Modeled Large-Eddy Simulation of High-Lift Configuration

This paper presents the assessment of several recent enhancements for a high-order wall-modeled large-eddy simulation (WMLES) approach and demonstrates order independence with a fixed data exchange location in the wall model. The two enhancements include the use of isotropic tetrahedral elements to improve accuracy and an explicit subgrid-scale model, the Vreman model, to improve accuracy and robustness. The [Formula: see text] study focused on the high-lift Common Research Model (HL-CRM) at the angle of attack of 19.57 deg, a benchmark problem from the 4th AIAA High-Lift Prediction Workshop. Solution polynomial orders of [Formula: see text], and 5 were used in the study. The study demonstrated [Formula: see text] independence in integrated forces, pitch moment, velocity profile in the wall-normal direction, and surface flow topology. It also showed that a [Formula: see text] order of at least 3 ([Formula: see text]) was needed to correctly predict the external inviscid flow and the surface flow topology. Thereafter, [Formula: see text] simulations over several other angles of attack demonstrated that the high-order WMLES approach can correctly predict the maximum lift and flow separation regions for HL-CRM with about 40 million degrees of freedom (DOF) compared to at least 250 million DOF required by second-order methods.

Engineering↗

Fast-forwarding quantum evolution

We investigate the problem of fast-forwarding quantum evolution, whereby the dynamics of certain quantum systems can be simulated with gate complexity that is sublinear in the evolution time. We provide a definition of fast-forwarding that considers the model of quantum computation, the Hamiltonians that induce the evolution, and the properties of the initial states. Our definition accounts for any asymptotic complexity improvement of the general case and we use it to demonstrate fast-forwarding in several quantum systems. In particular, we show that some local spin systems whose Hamiltonians can be taken into block diagonal form using an efficient quantum circuit, such as those that are permutation-invariant, can be exponentially fast-forwarded. We also show that certain classes of positive semidefinite local spin systems, also known as frustration-free, can be polynomially fast-forwarded, provided the initial state is supported on a subspace of sufficiently low energies. Last, we show that all quadratic fermionic systems and number-conserving quadratic bosonic systems can be exponentially fast-forwarded in a model where quantum gates are exponentials of specific fermionic or bosonic operators, respectively. Our results extend the classes of physical Hamiltonians that were previously known to be fast-forwarded, while not necessarily requiring methods that diagonalize the Hamiltonians efficiently. We further develop a connection between fast-forwarding and precise energy measurements that also accounts for polynomial improvements.

97 MATHEMATICS AND COMPUTING↗

Anatomy of information scrambling and decoherence in the integrable Sachdev-Ye-Kitaev model

The growth of information scrambling, captured by out-of-time-order correlation functions (OTOCs), is a central indicator of the nature of many-body quantum dynamics. Here, we compute analytically the complete time dependence of the OTOC for an integrable Sachdev-Ye-Kitaev (SYK) model, 𝑁 Majoranas with random two-body interactions of infinite range, coupled to a Markovian bath at finite temperature. In the limit of no coupling to the bath, the time evolution of scrambling experiences different stages. For 𝑡≲√𝑁, after an initial polynomial growth, the OTOC approaches saturation in a power-law fashion with oscillations superimposed. At 𝑡∼√𝑁, the OTOC reverses trend and starts to decrease linearly in time. The reason for this linear decrease is that, despite being a subleading 1/𝑁 effect, the OTOC in this region is governed by the spectral form factor of the antisymmetric couplings of the SYK model. The linear decrease stops at 𝑡∼2⁢𝑁, the Heisenberg time, where saturation occurs. The effect of the environment is an overall exponential decay of the OTOC for times longer than the inverse of the coupling strength to the bath. Here, the oscillations at 𝑡≲√𝑁 indicate lack of thermalization—a desired feature for better performance of quantum information devices.

Open quantum systems & decoherence↗

Calibration, measurement, and characterization of soil moisture dynamics in a central Amazonian tropical forest

Soil moisture plays a key role in hydrological, biogeochemical, and energy budgets of terrestrial ecosystems. Accurate soil moisture measurements in remote ecosystems such as the Amazon are difficult and limited because of logistical constraints. Time domain reflectometry (TDR) sensors are widely used to monitor soil moisture and require calibration to convert the TDR's dielectric permittivity measurement (K a ) to volumetric water content (θ v ). In this study, our objectives were to develop a field-based calibration of TDR sensors in an old-growth upland forest in the central Amazon, to evaluate the performance of the calibration, and then to apply the calibration to determine the dynamics of soil moisture content within a 14.2 -m-deep vertical soil profile. Depth-specific TDR calibration using local soils in a controlled laboratory setting yielded a novel K a –θ v third-degree polynomial calibration. The sensors were later installed to their specific calibration depth in a 14.2-m pit. The widely used K a –θ v relationship (Topp model) underestimated the site-specific θ v by 22–42%, indicating significant error in the model when applied to these well-structured, clay-rich tropical forest soils. The calibrated wet- and dry-season θ v data showed a variety of depth and temporal variations highlighting the importance of soil textural differentiation, root uptake depths, as well as event to seasonal precipitation effects. Data such as these are greatly needed for improving our understanding of ecohydrological processes within tropical forests and for improving models of these systems in the face of changing environmental conditions.

54 ENVIRONMENTAL SCIENCES↗

Uncertainty propagation in feed-forward neural network models

We develop new uncertainty propagation methods for feed-forward neural network architectures with leaky ReLU activation functions subject to random perturbations in the input vectors. In particular, we derive analytical expressions for the probability density function (PDF) of the neural network output and its statistical moments as a function of the input uncertainty and the parameters of the network, i.e., weights and biases. A key finding is that an appropriate linearization of the leaky ReLU activation function yields accurate statistical results even for large perturbations in the input vectors. This can be attributed to the way information propagates through the network. We also propose new analytically tractable Gaussian copula surrogate models to approximate the full joint PDF of the neural network output. To validate our theoretical results, we conduct Monte Carlo simulations and a thorough error analysis on a multi-layer neural network representing a nonlinear integro-differential operator between two polynomial function spaces. Our findings demonstrate excellent agreement between the theoretical predictions and Monte Carlo simulations.

MLP networks↗