SEARCH · Engineering Papers

Results for “lower and upper parameter bounds”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 19 records

All sources NASA NTRS DOE OSTI

DOE OSTI2022

Enhanced Multi-Dimensional Inversion Through Target-Specific Inversion Parameter Bounds With an Application to Crosswell Electromagnetic for Sequestration Monitoring

In geophysical inversions, lower and upper model parameter bounds are a means of solution stabilization. Further, constraints that intend to let only geologically plausible inverse solutions pass are amenable to lower and upper bounds. Reliable prior information is paramount to construct such bound constraints. It is common practice to narrow and widen bound intervals for regions of, respectively, more and less certain prior information. Contrary to this practice, we experiment with widened bound intervals in zones that are poorly resolved by a given survey configuration but where prior information would suggest structural anomalies of interest. The purpose of enlarged parameter bounds that correlate spatially with predefined targets is to let the inversion explore a larger solution space, thus increasing the potential to resolve otherwise hidden anomalies. Application of the method is based on a carbon-sequestration baseline (pre-injection) crosswell electromagnetic (EM) field survey at the Containment and Monitoring Institute Field Research Station (Alberta, Canada), where impeded measurements led to generally reduced sensitivities for the interwell region. Synthetic-data proofs of concept use augmented bounds designed to boost the resolution of artificial plume targets, indicating an enhanced illumination compared to constant bounds. Comparative field data inversions with spatially variable bounds constructed from prior resistivity and velocity information highlight non-horizontal baseline structures.

3-D electrical resistivity imaging↗

DOE OSTI2022

Improved geophysical monitoring of carbon sequestration through parameter linkage to reservoir modeling

Predictive reservoir modeling, even if present in the form of only basic hydrogeological model assumptions, is expected to accompany the majority of carbon capture and sequestration monitoring activities. It thus represents a source of prior information about the migration of injected fluids that can benefit geophysical survey planning and ensuing monitoring. Constraining the imaging of geophysical monitoring data with reservoir modeling is preferable over standalone geophysical imaging because of additional complementary hydrogeological information. However, fully coupled hydrogeophysical data inversion for flow-modeling parameters that control saturation predictions is an involved process. Within the context of three-dimensional electromagnetic (EM) inversion of data from borehole-to-surface layouts, we employ a "poor people's" alternative. The approach constrains geophysical inversion parameters through saturation predictions. The coupling is realized through spatially variable lower and upper parameter bounds that scale with gas saturation magnitudes, the latter provided by reservoir modeling. In this work, enhancement of three-dimensional time-lapse plume EM imaging is demonstrated for simulated sequestration into a depleted gas reservoir.

58 GEOSCIENCES↗

DOE OSTI2024

Stochastic finite volume method for uncertainty quantification of transient flow in gas pipeline networks

We develop a weakly intrusive framework to simulate the propagation of uncertainty in solutions of generic hyperbolic partial differential equation systems on graph-connected domains with nodal coupling and boundary conditions. The method is based on the Stochastic Finite Volume (SFV) approach and can be applied for uncertainty quantification (UQ) of the dynamical state of fluid flow over actuated transport networks. The numerical scheme has specific advantages for modeling intertemporal uncertainty in time-varying boundary parameters, which cannot be characterized by strict upper and lower (interval) bounds. We describe the scheme for a single pipe, and then formulate the controlled junction Riemann problem (JRP) that enables the extension to general network structures. In conclusion, we demonstrate the method's capabilities and performance characteristics using a standard benchmark test network.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Entropic order parameters in weakly coupled gauge theories

The entropic order parameters measure in a universal geometric way the statistics of non-local operators responsible for generalized symmetries. In this article, we compute entropic order parameters in weakly coupled gauge theories. To perform this computation, the natural route of evaluating expectation values of physical (smeared) non-local operators is prevented by known difficulties in constructing suitable smeared Wilson loops. We circumvent this problem by studying the smeared non-local class operators in the enlarged non-gauge invariant Hilbert space. This provides a generic approach for smeared operators in gauge theories and explicit formulas at weak coupling. In this approach, the Wilson and ’t Hooft loops are labeled by the full weight and co-weight lattices respectively. We study generic Lie groups and discuss couplings with matter fields. Smeared magnetic operators, as opposed to the usual infinitely thin ones, have expectation values that approach one at weak coupling. The corresponding entropic order parameter saturates to its maximum topological value, except for an exponentially small correction, which we compute. On the other hand, smeared ’t Hooft loops and their entropic disorder parameter are exponentially small. We verify that both behaviors match the certainty relation for the relative entropies. In particular, we find upper and lower bounds (that differ by a factor of 2) for the exact coefficient of the linear perimeter law for thin loops at weak coupling. This coefficient is unphysical/non-universal for line operators. We end with some comments regarding the RG flows of entropic parameters through perturbative beta functions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2023

Nuclear-Level Effective Theory of 𝜇 → 𝑒 Conversion

The Mu2e and COMET 𝜇 → 𝑒 conversion experiments are expected to significantly advance limits on new sources of charged lepton flavor violation. Almost all theoretical work in the field has focused on just two operators. However, general symmetry arguments lead to a 𝜇 → 𝑒 conversion rate with six response functions, each of which, in principle, is observable by varying nuclear properties of targets. We construct a nucleon-level nonrelativistic effective theory (NRET) to clarify the microscopic origin of these response functions and to relate rate measurements in different targets. This exercise identifies three operators and their small parameters that control the NRET operator expansion. We note inconsistencies in past treatments of these parameters. The NRET is technically challenging, involving 16 operators, several distorted electron partial waves, bound muon upper and lower components, and an exclusive nuclear matrix element. We introduce a trick for treating the electron Coulomb effects accurately, which enables us to include all of these effects while producing transition densities whose one-body matrix elements can be evaluated analytically, greatly simplifying the nuclear physics. We derive bounds on operator coefficients from existing and anticipated 𝜇 → 𝑒 conversion experiments. We discuss how similar NRET formulations have impacted dark matter phenomenology, noting that the tools this community has developed could be adapted for charged lepton flavor violation studies.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2025

Shear strength of a refractory high entropy alloy MoNbTaVW under high pressure

Radial X-ray diffraction (R-XRD) was performed in situ using a Panoramic Diamond Anvil Cell on the refractory high entropy alloy MoNbTaVW. The lattice parameters were determined through a Le Bail fit using a body centered cubic (BCC) lattice of symmetry $Im$$\overline{3}$$m$ (international space group number 229). Upper and lower bounds to the shear strength were determined up to 80 GPa nonhydrostatically using copper as a pressure standard. The equation of state was derived at the magic angle ψ = 54.7° and yielded a bulk modulus of K 0 = 220.8 ± 1.65 GPa. The experimental lattice parameters and bulk modulus match closely with corresponding density functional theory (DFT) calculations. The BCC phase remains stable up to the highest pressure of 80 GPa studied and is shown to be elastically anisotropic. The shear strength was found to saturate around 70 GPa with a value of τ = 1.75 GPa, and the shear moduli are presented in different limits of iso-strain and iso-stress.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Multi-parametric analysis for mixed integer linear programming: An application to transmission upgrade and congestion management

Upgrading the capacity of existing transmission lines is essential for meeting the growing energy demands, facilitating the integration of renewable energy, and ensuring the security of the transmission system. This study focuses on the selection of lines whose capacities and by how much should be expanded from the perspective of the Independent System Operators (ISOs) to minimize the total system cost. We employ advanced multi-parametric programming and an enhanced branch-and-bound algorithm to address complex mixed-integer linear programming (MILP) problems, considering multi-period time constraints and physical limitations of generators and transmission lines. To characterize the various decisions in transmission expansion, we model the increased capacity of existing lines as parameters within a specified range. This study first relaxes the binary variables to continuous variables and applies the Lagrange method and Karush-Kuhn-Tucker (KKT) conditions to obtain optimal solutions and identify critical regions associated with active and inactive constraints. Moreover, we extend the traditional branch-and-bound (B&B) method by determining the problem’s upper and lower bounds at each node of the B&B decision tree, helping to manage computational challenges in large-scale MILP problems. Here, we compare the difference between the upper and lower bounds to obtain an approximate optimal solution within the decision-makers’ tolerable error range. In addition, the first derivative of the objective function on the parameters of each line is used to inform the selection of lines for easing congestion and maximizing social welfare. Finally, the capacity upgrades are selected by weighing the reductions in system costs against the expense of upgrading line capacities. The findings are supported by numerical simulations and provide transmission-line planners with decision-making guidance.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2021

A prediction interval method for uncertainty quantification of regression models

This paper considers calculation of prediction intervals (PIs) by neural networks (NNs) for quantifying uncertainty in regression tasks, so as to provide fast, accurate and robust emulators to accelerate scientific simulations. We propose a novel method to learn lower and upper bounds of the PI using independent NNs without defining an exclusive loss. Our method requires no distributional assumption, does not introduce extra hyper-parameters, and can effectively identify out-of-distribution samples and quantify their uncertainty. We demonstrate advantages of our method using a benchmark problem and two real-world scientific applications.

Zhang, Pei↗

DOE OSTI2023

Lepto-axiogenesis and the scale of supersymmetry

If the Peccei-Quinn field containing the QCD axion undergoes rotations in the early universe, the dimension-five operator responsible for neutrino masses can generate a lepton asymmetry that ultimately gives rise to the observed baryon asymmetry of the Universe. This lepto-axiogenesis scenario requires a flat potential for the radial direction of the Peccei-Quinn field, naturally realized in supersymmetric models. We carefully compute the efficiency of this mechanism for the Dine-Fischler-Srednicki-Zhitnitsky (DFSZ) and Kim-Shifman-Vainshtein-Zakharov (KSVZ) axion models and place lower bounds on the masses of scalar superpartners required to reproduce the observed baryon asymmetry. For the KSVZ model, we find an efficiency for generation of the asymmetry six times larger than the previously extant computation after including scattering channels involving superpartners. In this case, the superpartner scale should be above 30 TeV for a domain wall number of one; the lower bound weakens for larger domain wall numbers. We find that the superpartner mass scale may also be as low as ~ 30 TeV for the DFSZ model. In all cases, the lower bound on the superpartner masses is inversely proportional to the sum of the squares of the neutrino masses and so can strengthen as the upper bound on the neutrino mass improves. We identify the parameter space where the axion rotation can simultaneously produce axion dark matter via kinetic misalignment; in this case it is possible to put an upper bound of order PeV on the masses of scalar superpartners.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

On the Impact of Bus Dwelling on Macroscopic Fundamental Diagrams

Network macroscopic fundamental diagrams (MFDs) have recently been shown to exist in real-world urban traffic networks. When present, MFDs can be used to model traffic dynamics within an urban network by dividing the network into a set of spatially compact homogeneous regions and tracking the average level of congestion in each region. Existing analytical methods to estimate MFD mostly focus on the behavior of a single type of vehicle and do not capture the patterns of mixed traffic (e.g., cars and buses). The existence of buses matters since a bus will block the movements of other vehicles when it dwells at the bus stop. This paper proposes an analytical method to estimate the impact of bus dwelling on a network’s MFD based on the network’s geometric features, traffic control strategies, and bus operation parameters, and validates the performance of the proposed method using simulations based on microscopic traffic models. Comparisons of the analytical and simulation results show that the proposed analytical method can generally provide a good estimate of the lower bound and upper bound of the network’s MFD.

Xu, Guanhao↗

DOE OSTI2026

Testing Scale-Dependent Modified Gravity with DESI DR1

The Dark Energy Spectroscopic Instrument (DESI) provides an unprecedented opportunity to test deviations from general relativity (GR) that introduce a new physical scale within its redshift range. Using the connection between a Yukawa-like potential and the Hu-Sawicki $f(R)$ model, we place strong constraints on the range of a hypothetical fifth force mediated by a massive scalar field. We analyze the power spectrum measurements from DESI Data Release 1 using a baseline EFT model that employs the fkpt approach for the loop integrals. We find no evidence for deviations from GR and obtain the constraint $\log_{10} |f_{R_0}| < -4.59$ (95% C.L.). This corresponds to an upper bound at redshift zero on the scale at which corrections to GR become important, $λ< 17.81$ Mpc, or equivalently, a lower bound on the mass of the additional gravitational mediator of $m_ϕ> 3.60 \times 10^{-31}$ eV. We find that the modified gravity parameter $f_{R_0}$ is largely orthogonal to the cosmological parameters in the model, such that no additional projection effects relative to the GR case are introduced in this Full-Shape analysis. Furthermore, a second modified gravity parameter, the power index $n$, which modulates the time-variation of the associated mass, is found to be consistent with previous analyses that fixed it to unity. Adding DESI BAO data or other cosmological probes does not significantly change these results. The conclusions remain similar if the background evolution is described by evolving dark energy instead of a cosmological constant. Additionally, we test the robustness of the baseline model by varying the maximum wavenumber used in the Full-Shape analysis and analyzing the DESI targets separately. Finally, we analyze the degeneracies between the modified-gravity parameters and the sum of neutrino masses.

Gonzalez, D. [Guanajuato U.; UNAM, CCF] (ORCID:000↗

DOE OSTI2023

Constraints on a split superconducting transition under uniaxial strain in Sr 2 RuO 4 from scanning SQUID microscopy

More than two decades after the discovery of superconductivity in Sr 2 RuO 4 , it is still unclear whether the order parameter has a single component or two degenerate components. For any two-component scenario, application of uniaxial strain is expected to lift the degeneracy, generating two distinct phase transitions. The presence of a second (lower-temperature) transition may be observable by probes that are sensitive to changes in the London penetration depth, λ, as a function of temperature, T. Furthermore, we use scanning SQUID microscopy combined with a uniaxial strain device to test for a second transition under strain. We only observe a single transition. Within the temperature range where a second transition has been suggested by μSR measurements, we further place a tight upper bound of less than 1% on the change in the zero temperature superfluid density n s ∝ λ –2 (0) due to a second transition, suggesting that such a transition does not occur. These results constrain theories of the order parameter in Sr 2 RuO 4 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

Mercury’s Chaotic Secular Evolution as a Subdiffusive Process

Abstract Mercury’s orbit can destabilize, generally resulting in a collision with either Venus or the Sun. Chaotic evolution can causeg 1 to decrease to the approximately constant value ofg 5 and create a resonance. Previous work has approximated the variation ing 1 as stochastic diffusion, which leads to a phenomological model that can reproduce the Mercury instability statistics of secular andN-body models on timescales longer than 10 Gyr. Here we show that the diffusive model significantly underpredicts the Mercury instability probability on timescales less than 5 Gyr, the remaining lifespan of the solar system. This is becauseg 1 exhibits larger variations on short timescales than the diffusive model would suggest. To better model the variations on short timescales, we build a new subdiffusive phenomological model forg 1 . Subdiffusion is similar to diffusion but exhibits larger displacements on short timescales and smaller displacements on long timescales. We choose model parameters based on the behavior of theg 1 trajectories in theN-body simulations, leading to a tuned model that can reproduce Mercury instability statistics from 1–40 Gyr. This work motivates fundamental questions in solar system dynamics: why does subdiffusion better approximate the variation ing 1 than standard diffusion? Why is there an upper bound ong 1 , but not a lower bound that would prevent it from reachingg 5 ?

Astronomy & Astrophysics↗

DOE OSTI2022

Variational quantum algorithm for estimating the quantum Fisher information

The quantum Fisher information (QFI) quantifies the ultimate precision of estimating a parameter from a quantum state and can be regarded as a reliability measure of a quantum system as a quantum sensor. However, estimation of the QFI for a mixed state is in general a computationally demanding task. In this paper we present a variational quantum algorithm called variational quantum Fisher information estimation (VQFIE) to address this task. By estimating lower and upper bounds on the QFI, based on bounding the fidelity, VQFIE outputs a range in which the actual QFI lies. This result can then be used to variationally prepare the state that maximizes the QFI, for the application of quantum sensing. In contrast to previous approaches, VQFIE does not require knowledge of the explicit form of the sensor dynamics. We simulate the algorithm for a magnetometry setup and demonstrate the tightening of our bounds as the state purity increases. For this example, we compare our bounds with literature bounds and show that our bounds are tighter.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Discontinuous Aligned Carbon Fiber Intermediates for Automotive and Related Applications

This work focused on preferentially aligning discontinuous carbon fibers in wet-laid or air-laid processes. It is well known that aligned fibers provides higher directional strength and stiffness. Discontinuous fibers further allow higher degree of draw and formability as the gaps in the fibers allow for higher material movement. The current processes are limited in their ability to align carbon fibers during processing. The aligned fibers have several benefits - (a) in applications where chopped fibers can replace continuous fibers for targeted strength and stiffness metrics, but at a substantially reduced cost; (b) they can tolerate deeper draws than continuous fiber composites in thermo-stamping and compression molding processes; (c) they can be tailored for pultrusion and unidirectional applications. Although pultrusion is primarily a process that adopts continuous fibers, stitch bonded entangled discontinuous fibers can provide unique intermediates. This is analogous to natural coir fibers which get aligned and entangled to produce ropes/rods for example, (d) they can be processed in cross-ply and multi-directional formats, like composite laminates. In this work Neenah Paper partnered with IACMI, UT and ORNL to evaluate structure-process-property relationships with Zoltek carbon fiber. A few process parameters such as machine speed, weight basis, fiber length, effect of fiber sizing, direction of mat lay-up etc. were investigated. The produced mats were converted to thermoplastic composite laminates using polyamide 6 (PA6, nylon) resin. The specific objective of this project is to produce a wet-laid nonwoven carbon fiber mat with a high degree of unidirectional fiber alignment, using discontinuous carbon fibers. The report provides details about the processing, characterization, and lower-upper bound properties.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Environmental Influences on Deep Convective Upscale Growth Rate in Central Argentina From a Convection‐Permitting Simulation

This study uses a convection‐permitting model simulation to describe the environmental conditions under which convective upscale growth occurs in central Argentina, particularly examining environmental parameters when deep convection initially forms that could differentiate the rate of initial upscale growth. Simulated mesoscale convective systems (MCSs) are separated into slow and rapid growth by the rate of spatial growth from convection initiation until reaching the MCS scale. A low‐level jet (LLJ) is found more frequently near the deep convection that experiences rapid growth to an MCS, but its presence alone is not predictive of rapid growth. Using spatially‐averaged parameters, we find that rapid growth to MCSs also occurs in environments that are significantly more thermodynamically favorable with greater low‐level moisture and instability. Fewer significant differences are found in the kinematic environment with only the 0–2 km vertical wind shear magnitude being significantly larger for rapid growth MCSs compared to slow growth MCSs, potentially related to LLJs often peaking near this height. When focusing only on MCSs with the slowest and fastest growth rates, elevated‐layer shear is significantly smaller for very rapid growth MCSs, suggesting elevated‐layer shear may help discriminate between the upper and lower bounds of growth rate. Finally, when upscale growth occurs near the Sierras de Córdoba (SDC) with a LLJ present, rapid growth is also supported by favorable wind shear orientation. However, this does not hold for upscale growth occurring away from the SDC, highlighting the importance of interpreting shear direction relative to the orientation of features initiating deep convection.

Sasaki, Clayton R. S. [Univ. of Washington, Seattl↗

DOE OSTI2024

Search for photons above 10 18 eV by simultaneously measuring the atmospheric depth and the muon content of air showers at the Pierre Auger Observatory

The Pierre Auger Observatory is the most sensitive instrument to detect photons with energies above 1 0 17 eV . It measures extensive air showers generated by ultrahigh energy cosmic rays using a hybrid technique that exploits the combination of a fluorescence detector with a ground array of particle detectors. The signatures of a photon-induced air shower are a larger atmospheric depth of the shower maximum ( X max ) and a steeper lateral distribution function, along with a lower number of muons with respect to the bulk of hadron-induced cascades. In this work, a new analysis technique in the energy interval between 1 and 30 EeV ( 1 EeV = 1 0 18 eV ) has been developed by combining the fluorescence detector-based measurement of X max with the specific features of the surface detector signal through a parameter related to the air shower muon content, derived from the universality of the air shower development. No evidence of a statistically significant signal due to photon primaries was found using data collected in about 12 years of operation. Thus, upper bounds to the integral photon flux have been set using a detailed calculation of the detector exposure, in combination with a data-driven background estimation. The derived 95% confidence level upper limits are 0.0403, 0.01113, 0.0035, 0.0023, and 0.0021 km − 2 sr − 1 yr − 1 above 1, 2, 3, 5, and 10 EeV, respectively, leading to the most stringent upper limits on the photon flux in the EeV range. Compared with past results, the upper limits were improved by about 40% for the lowest energy threshold and by a factor 3 above 3 EeV, where no candidates were found and the expected background is negligible. The presented limits can be used to probe the assumptions on chemical composition of ultrahigh energy cosmic rays and allow for the constraint of the mass and lifetime phase space of super-heavy dark matter particles. Published by the American Physical Society 2024

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Provably accurate simulation of gauge theories and bosonic systems

Quantum many-body systems involving bosonic modes or gauge fields have infinite-dimensional local Hilbert spaces which must be truncated to perform simulations of real-time dynamics on classical or quantum computers. To analyze the truncation error, we develop methods for bounding the rate of growth of local quantum numbers such as the occupation number of a mode at a lattice site, or the electric field at a lattice link. Our approach applies to various models of bosons interacting with spins or fermions, and also to both abelian and non-abelian gauge theories. We show that if states in these models are truncated by imposing an upper limit Λ on each local quantum number, and if the initial state has low local quantum numbers, then an error at most ϵ can be achieved by choosing Λ to scale polylogarithmically with ϵ − 1 , an exponential improvement over previous bounds based on energy conservation. For the Hubbard-Holstein model, we numerically compute a bound on Λ that achieves accuracy ϵ , obtaining significantly improved estimates in various parameter regimes. We also establish a criterion for truncating the Hamiltonian with a provable guarantee on the accuracy of time evolution. Building on that result, we formulate quantum algorithms for dynamical simulation of lattice gauge theories and of models with bosonic modes; the gate complexity depends almost linearly on spacetime volume in the former case, and almost quadratically on time in the latter case. We establish a lower bound showing that there are systems involving bosons for which this quadratic scaling with time cannot be improved. By applying our result on the truncation error in time evolution, we also prove that spectrally isolated energy eigenstates can be approximated with accuracy ϵ by truncating local quantum numbers at Λ = polylog ( ϵ − 1 ) .

Tong, Yu↗