Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “constrained regions”

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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Memory-efficient emulation of physical tabular data using quadtree decomposition

Computationally expensive functions are sometimes replaced in simulations with an emulator that approx-imates the true function (e.g., equations of state, wavelength-dependent opacity, or composition-dependent materials properties). For functions that have a constrained domain of interest, this can be done by discretizing the domain and performing a local interpolation on the tabulated function values of each local domain. For these so-called tabular data methods, the method of discretizing the domain and mapping the input space to each subdomain can drastically influence the memory and computational costs of the emulator. This is especially true for functions that vary drastically in different regions. We present a method for domain discretization and mapping that utilizes quadtrees, which results in significant reductions in the size of the emulator with minimal increases to computational costs or loss of global accuracy. We apply our method to the electron-positron Helmholtz free energy equation of state and show over an order of magnitude reduction in memory costs for reasonable levels of numerical accuracy.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Regional Moment Tensor Inversion Using Rotational Observations

There are benefits from the addition of rotational motions to translational displacements for moment tensor (MT) inversions. The rotational radiation pattern is orthogonal to the shear radiation pattern, thus incorporating rotations is equivalent to gaining another observation point on the focal sphere. We demonstrated this by simulating curl and displacement wavefields for a regional distance station. Thus, one 6-C station (3-Component translational + 3-Component rotational) gathers the same information on radiation pattern as two 3-C stations at 90° azimuth from one another along the focal plane axis, which is sometimes difficult to obtain when restricted to surface sensors. We added rotational Green's functions to a regional MT inversion scheme (long-period, time-domain, and linear inversion) by computing spatial gradients from f-κ reflectivity synthetics. For rotational data, we used Array Derived Rotations (ADR) from Piñon Flats Observatory Array and Golay array deployed during IRIS Community Wavefield Demonstration Experiment. The hope is to ultimately use compact and field deployable broadband rotational seismometers instead of ADRs. Rotational motions were predicted from well-constrained deviatoric MT solutions of nine earthquakes recorded by the arrays and other seismic networks. Additionally, we formed three station sparse datasets with (3-C and 6-C) and without rotational ground motions (3-C only) for eight earthquakes testing the benefits of including rotational waveforms in MT inversions when they are equally weighted with translational waveforms. Adding rotational motions improved the double-couple components and reduced compensated-linear-vector-dipole and isotropic components in full-MT solutions using sparse datasets with poor coverage.

58 GEOSCIENCES↗
DOE OSTI2024

Quantifying Variability and Controls of Riverine Dissolved Organic Carbon Exported to Arctic Coastal Margins of North America (Final Report)

This project involved implementation and application of a coupled permafrost hydrology and dissolved organic carbon process models to investigate how spatial and seasonal variations in terrestrial hydrology and soil freeze/thaw dynamics influence the mobilization, loading, and export of organic carbon to the stream network for selected arctic basins across northern Alaska and northwest Canada. The model simulations were constrained by detailed observations of in- stream chemistry, soil active layer profile moisture and temperature dynamics, streamflow, soil carbon inventories and satellite microwave remote sensing based assessments of surface soil freeze-thaw dynamics. We developed and applied the numerical modeling and data analysis, incorporating observed data for calibration and validation, in order to investigate the terrestrial hydrology, permafrost dynamics, and associated DOC production and loading to rivers across a region encompassing watersheds draining to the coast. The project produced six publications and three datasets archived in public repositories.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Bayesian Exploration of Phenomenological EoS of Neutron/Hybrid Stars with Recent Observations

The description of the stellar interior of compact stars remains as a big challenge for the nuclear astrophysics community. The consolidated knowledge is restricted to density regions around the saturation of hadronic matter ρ 0 =2.8 × 10 14 g cm -3 , regimes where our nuclear models are successfully applied. As one moves towards higher densities and extreme conditions up to the quark/gluons deconfinement, little can be said about the microphysics of the equation of state (EoS). Here, we employ a Markov Chain Monte Carlo (MCMC) strategy to access the variability at high density regions of polytropic piecewise models for neutron star (NS) EoS or possible hybrid stars, i.e., a NS with a small quark-matter core. With a fixed description of the hadronic matter for low density, below the nuclear saturation density, we explore a variety of models for the high density regimes leading to stellar masses near to 2.5 M ⊙ , in accordance with the observations of massive pulsars. The models are constrained, including the observation of the merger of neutrons stars from VIRGO-LIGO and with the pulsar observed by NICER. In addition, we also discuss the possibility of the use of a Bayesian power regression model with heteroscedastic error. The set of EoS from the Laser Interferometer Gravitational-Wave Observatory (LIGO) was used as input and treated as the data set for the testing case.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Understanding the Charging Flexibility of Shared Automated Electric Vehicle Fleets

The combined anticipated trends of vehicle sharing (ride-hailing), automated control, and powertrain electrification are poised to disrupt the current paradigm of predominately owner-driven gasoline vehicles with low levels of utilization. Shared, automated, electric vehicle (SAEV) fleets offer the potential for lower cost and emissions and have garnered significant interest among the research community. While promising, unmanaged operation of these fleets may lead to unintended negative consequences. One potentially unintended consequence is a high quantity of SAEVs charging during peak demand hours on the electric grid, potentially increasing the required generation capacity. This research explores the flexibility associated with charging loads demanded by SAEV fleets in response to servicing personal mobility travel demands. Travel demand is synthesized in four major United States metropolitan areas: Detroit, MI; Austin, TX; Washington, DC; and Miami, FL. In each of these four cities, SAEV simulations are performed using local projected electricity prices from the Regional Energy Deployment System (ReEDS) for a handful of supply side scenarios. A coordinated charging model is introduced that seeks to reduce fleet charging costs in response to time-varying electricity prices and increasing shares of personal trip demand served (between 1% and 25% of all metro trips served by the SAEV fleet). Simulation results are presented across various scenarios assuming fleetwide coordination to minimize charging energy costs while constrained by offering adequate mobility service to fleet customers. The results indicate that the SAEV charging load is highly flexible; energy costs were shown to reduce between 13% and 46% across a wide range of simulated scenarios. In addition, these savings were realized without detrimentally impacting the fleet’s ability to service trips.

autonomous vehicle↗
DOE OSTI2025

Uncertainty quantification and optimization of precipitating hydrometeor parameters for winter precipitation in a cloud microphysics scheme

The precipitating hydrometeor parameters used in cloud microphysics schemes carry inherent uncertainties. The quantification of these uncertainties, together with parameter optimization, can significantly improve precipitation forecasts. This study investigates the effects of 13 parameters in the Weather Research and Forecasting (WRF) Double-Moment 6-class (WDM6) microphysics scheme, which define the hydrometeor characteristics such as fall velocity–diameter and mass–diameter relationships, as well as the shape parameter of the drop size distribution for precipitating particles such as rain, snow, and graupel on simulated winter precipitation. A comparison between the model's pre-defined parameters and observations from the International Collaborative Experiments for the PyeongChang 2018 Olympic and Paralympic winter games (ICE-POP 2018) field campaign reveals that the fall velocity–diameter relationship for rain, the mass–diameter relationships for snow and graupel, and the shape parameters for all precipitating particles in the WDM6 scheme deviate from the median values observed by the two-dimensional video disdrometer (2DVD). To quantify parameter sensitivities, a perturbed parameter ensemble (PPE) of 256 simulations was conducted within parameter ranges constrained by 2DVD observations for three winter precipitation cases. Bayesian optimization was then applied to identify parameter sets that minimized the root mean square error (RMSE) for each case, achieving reductions of up to 30.2 %. These results demonstrate that ensemble-based uncertainty quantification and parameter optimization can help identify key parameters and provide a pathway to improving precipitation simulation performance. In addition, measurement sites can be strategically selected based on regions that show high sensitivity to variations in hydrometeor characteristic parameters.

Bayesian optimization↗
DOE OSTI2021

Evidence of hydrogen-helium immiscibility at Jupiter-interior conditions

The phase behaviour of warm dense hydrogen–helium (H–He) mixtures affects our understanding of the evolution of Jupiter and Saturn and their interior structures. For example, precipitation of He from a H–He atmosphere at about 1–10 megabar and a few thousand kelvin has been invoked to explain both the excess luminosity of Saturn, and the depletion of He and neon (Ne) in Jupiter’s atmosphere as observed by the Galileo probe. But despite its importance, H–He phase behaviour under relevant planetary conditions remains poorly constrained because it is challenging to determine computationally and because the extremes of temperature and pressure are difficult to reach experimentally. Here we report that appropriate temperatures and pressures can be reached through laser-driven shock compression of H 2 –He samples that have been pre-compressed in diamond-anvil cells. This allows us to probe the properties of H–He mixtures under Jovian interior conditions, revealing a region of immiscibility along the Hugoniot. A clear discontinuous change in sample reflectivity indicates that this region ends above 150 gigapascals at 10,200 kelvin and that a more subtle reflectivity change occurs above 93 gigapascals at 4,700 kelvin. Considering pressure–temperature profiles for Jupiter, these experimental immiscibility constraints for a near-protosolar mixture suggest that H–He phase separation affects a large fraction—we estimate about 15 per cent of the radius—of Jupiter’s interior. Furthermore, this finding provides microphysical support for Jupiter models that invoke a layered interior to explain Juno and Galileo spacecraft observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

The Jefferson Lab Eta Factory (JEF) experiment

The Jefferson Lab Eta Factory (JEF) experiment is an upcoming experiment designed to run in Hall D at Jefferson Lab using an upgraded GlueX spectrometer to study various decays of the ? meson. The GlueX spectrometer consists of a ?2 Tesla solenoid magnet housing a liquid hydrogen target and drift chambers used for tracking charged particles and an array of lead glass blocks (the Forward Calorimeter (FCAL)) downstream of the magnet for detecting neutral particles. The decays of the ? meson can be used to measure the light quark mass ratio via the ???+???0 channel and allow access to higher-order terms in Chiral Perturbation Theory via the ???0?? channel. The decays can also be used to constrain new charge-conjugation violating/parity conserving (CVPC) reactions and to search for signatures of dark matter. In particular, the ???0?? channel can be used to search for lepto-phobic dark vector (B) bosons in the reaction ??B? (B??0?) or dark scalar (S) bosons in the reaction ???0S (S???). Studying the rare radiative decay channel ???0?? requires replacing the 4×4×45 cm3 lead glass blocks in the inner region of the FCAL with 2×2×20 cm3 lead tungstate crystals, which will provide twice better position and energy resolution. This talk will describe the JEF physics program and the upgrade to the FCAL.

Taylor, Simon↗
DOE OSTI2023

Constraints on the Gamma-Ray Emission from Small Solar System Bodies with the Fermi Large Area Telescope Data

All known small solar system bodies have diameters between a few meters and a few thousands of kilometers. Based on the collisional evolution of solar system bodies, a larger number of asteroids with diameters down to ~2 m is thought to exist. As all solar system bodies, small bodies can be passive sources of high-energy gamma-rays, produced by the interaction of energetic cosmic rays impinging on their surfaces. Since the majority of known asteroids are in orbits between Mars and Jupiter (in a region known as the Main Belt), we expect them to produce a diffuse emission close to the ecliptic plane. In this work, we have studied the gamma-ray emission coming from the ecliptic using the data collected by the Large Area Telescope (LAT) onboard the Fermi satellite. We have fit the results with simulations of the gamma-ray intensity at the source level (calculated with the software FLUKA) to constrain the small solar system bodies population. Finally, we have proposed a model describing the distribution of asteroid sizes and we have used the LAT data to constrain the gamma-ray emission expected from this model and, in turn, on the model itself.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Transforming US agriculture for carbon removal with enhanced weathering

Abstract Enhanced weathering (EW) with agriculture uses crushed silicate rocks to drive carbon dioxide removal (CDR) 1,2 . If widely adopted on farmlands, it could help achieve net-zero emissions by 2050 2–4 . Here we show, with a detailed US state-specific carbon cycle analysis constrained by resource provision, that EW deployed on agricultural land could sequester 0.16–0.30 GtCO 2 yr −1 by 2050, rising to 0.25–0.49 GtCO 2 yr −1 by 2070. Geochemical assessment of rivers and oceans suggests effective transport of dissolved products from EW from soils, offering CDR on intergenerational timescales. Our analysis further indicates that EW may temporarily help lower ground-level ozone and concentrations of secondary aerosols in agricultural regions. Geospatially mapped CDR costs show heterogeneity across the USA, reflecting a combination of cropland distance from basalt source regions, timing of EW deployment and evolving CDR rates. CDR costs are highest in the first two decades before declining to about US$100–150 tCO 2 −1 by 2050, including for states that contribute most to total national CDR. Although EW cannot be a substitute for emission reductions, our assessment strengthens the case for EW as an overlooked practical innovation for helping the USA meet net-zero 2050 goals 5,6 . Public awareness of EW and equity impacts of EW deployment across the USA require further exploration 7,8 and we note that mobilizing an EW industry at the necessary scale could take decades.

Science & Technology - Other Topics↗
DOE OSTI2022

Global wave energy resource classification system for regional energy planning and project development

Efforts to streamline and codify wave energy resource characterization and assessment for regional energy planning and wave energy converter (WEC) project development have motivated the recent development of resource classification systems. Given the unique interplay between WEC absorption and resource attributes, viz, available wave power frequency, directionality, and seasonality, various consensus resource classification metrics have been introduced. However, the main international standards body for the wave energy industry has not reached consensus on a wave energy resource classification system designed with clear goals to facilitate resource assessment, regional energy planning, project site selection, project feasibility studies, and selection of WEC concepts or archetypes that are most suitable for a given wave energy climate. In this work, a primary consideration of wave energy generation is the available energy that can be captured by WECs with different resonant frequency and directional bandwidths. Therefore, the proposed classification system considers combinations of three different wave power classifications: the total wave power, the frequency-constrained wave power, and the frequency-directionally constrained wave power. The dominant wave period bands containing the most wave power are sub-classification parameters that provide useful information for designing frequency and directionally constrained WECs. The bulk of the global wave energy resource is divided into just 22 resource classes representing distinct wave energy climates that could serve as a common language and reference framework for wave energy resource assessment if codified within international standards.

16 TIDAL AND WAVE POWER↗
DOE OSTI2024

Congo Basin Water Balance and Terrestrial Fluxes Inferred From Satellite Observations of the Isotopic Composition of Water Vapor

Large spatio-temporal gradients in the Congo basin vegetation and rainfall are observed. However, its water-balance (evapotranspiration minus precipitation, or ET - P) is typically measured at basin-scales, limited primarily by river-discharge data, spatial resolution of terrestrial water storage measurements, and poorly constrained ET. We use observations of the isotopic composition of water vapor to quantify the spatio-temporal variability of net surface water fluxes across the Congo Basin between 2003 and 2018. These data are calibrated at basin scale using satellite gravity and total Congo river discharge measurements and then used to estimate time-varying ET - P over four quadrants representing the Congo Basin, providing first estimates of this kind for the region. We find that the multi-year record, seasonality, and interannual variability of ET - P from both the isotopes and the gravity/river discharge based estimates are consistent. Additionally, we use precipitation and gravity-based estimates with our water vapor isotope-based ET - P to calculate time and space averaged ET and net river discharge within the Congo Basin. These quadrant-scale moisture flux estimates indicate (a) substantial recycling of moisture in the Congo Basin (temporally and spatially averaged ET/P > 70%), consistent with models and visible light-based ET estimates, and (b) net river outflow is largest in the Western Congo where there are more rivers and higher flow rates. Our results confirm the importance of ET in modulating the Congo water cycle relative to other water sources.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Splitting regions and shrinking islands from higher point constraints

We study constraints from higher-point amplitudes on 2 → 2 scattering in the context of effective field theory (EFT) using the perturbative numerical S-matrix bootstrap. Specifically, we investigate the class of weakly coupled EFTs with amplitudes that obey the hidden zero and split conditions that are known to hold both for Tr(Φ 3 ) theory and for certain string tree amplitudes, including at 4-point the beta function. Requiring the splitting condition for the 5-point amplitude not only fixes nearly all its contact terms, but it also imposes non-linear constraints among the 4-point EFT Wilson coefficients. When included in the bootstrap, the resulting allowed region consistent with positivity is no longer convex but is restricted to a smaller non-convex region — which has a sharp corner near the string beta function! Assuming the absence of an infinite spin tower at the mass gap, the allowed region bifurcates into a trivial region (with states only above a chosen cutoff) and an island that continues to shrink around the string as more constraints are included in the bootstrap. The numerics indicate that in the absence of single-mass infinite spin towers the string beta function is the unique 4-point amplitude compatible with hidden zero and the 5-point splitting constraints. The analysis provides a prototype example for how features of higher-point amplitudes constrain the bootstrap of 4-point amplitudes.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Eigenfrequency constrained topology optimization of finite strain hyperelastic structures

This paper incorporates hyperelastic materials, nonlinear kinematics, and preloads in eigenfrequency constrained density–based topology optimization. The formulation allows for initial finite deformations and subsequent small harmonic oscillations. The optimization problem is solved by the method of moving asymptotes, and the gradients are calculated using the adjoint method. Both simple and degenerate eigenfrequencies are considered in the sensitivity analysis. A well-posed topology optimization problem is formulated by filtering the volume fraction field. Numerical issues associated with excessive distortion and spurious eigenmodes in void regions are reduced by removing low volume fraction elements. The optimization objective is to maximize stiffness subject to a lower bound on the fundamental eigenfrequency. Numerical examples show that the eigenfrequencies drastically change with the load magnitude, and that the optimization is able to produce designs with the desired fundamental eigenfrequency.

42 ENGINEERING↗
DOE OSTI2023

Melting and defect transitions in FeO up to pressures of Earth’s core-mantle boundary

The high-pressure melting curve of FeO controls key aspects of Earth’s deep interior and the evolution of rocky planets more broadly. However, existing melting studies on wüstite were conducted across a limited pressure range and exhibit substantial disagreement. Here we use an in-situ dual-technique approach that combines a suite of >1000 x-ray diffraction and synchrotron Mössbauer measurements to report the melting curve for Fe 1-x O wüstite to pressures of Earth’s lowermost mantle. We further observe features in the data suggesting an order-disorder transition in the iron defect structure several hundred kelvin below melting. This solid-solid transition, suggested by decades of ambient pressure research, is detected across the full pressure range of the study (30 to 140 GPa). At 136 GPa, our results constrain a relatively high melting temperature of 4140 ± 110 K, which falls above recent temperature estimates for Earth’s present-day core-mantle boundary and supports the viability of solid FeO-rich structures at the roots of mantle plumes. The coincidence of the defect order-disorder transition with pressure-temperature conditions of Earth’s mantle base raises broad questions about its possible influence on key physical properties of the region, including rheology and conductivity.

58 GEOSCIENCES↗
DOE OSTI2026

Study of Nuclear Structure Functions at Jefferson Lab

This thesis presents a study of nuclear modification of quark distributions (the EMC effect) using inclusive electron-scattering data from Jefferson Lab Hall C, with emphasis on experiment E12-10-008 (XEM2). The analysis spans nuclei from light to heavy targets and combines measurements at multiple spectrometer settings to constrain EMC ratios over a broad range in Bjorken-x and Q2. A complete analysis framework was developed to extract charge-normalized and efficiency-corrected yields, including detector calibrations, beam-current calibration, density-loss corrections for cryogenic targets, background subtraction, radiative and Coulomb corrections, and systematic studies. Instrumental cross-checks, including HMS–SHMS comparisons and reconstruction studies, show that residual spectrometer differences are predominantly multiplicative in the kinematic region relevant to the EMC analysis. The extracted EMC ratios are broadly consistent with previous measurements while extending coverage across many nuclei. The EMC slope increases from light to heavy nuclei and shows signs of saturation at large A. After accounting for the dominant A dependence, the results are consistent with no isospin dependence, but also consistent with the predicted modification from the isovector model.

Sharda, Abhyuday [Univ. of Tennessee, Knoxville, T↗
DOE OSTI2021

First measurement of the |t|-dependence of coherent J/ψ photonuclear production

The first measurement of the cross section for coherent J/ψ photoproduction as a function of |t|, the square of the momentum transferred between the incoming and outgoing target nucleus, is presented. The data were measured with the ALICE detector in ultra-peripheral Pb–Pb collisions at a centre-of-mass energy per nucleon pair $\sqrt{s_{NN}}$=5.02TeV with the J/ψ produced in the central rapidity region |y|<0.8, which corresponds to the small Bjorken-x range (0.3-1.4)×10 -3 . The measured |t|-dependence is not described by computations based only on the Pb nuclear form factor, while the photonuclear cross section is better reproduced by models including shadowing according to the leading-twist approximation, or gluon-saturation effects from the impact-parameter dependent Balitsky–Kovchegov equation. These new results are therefore a valid tool to constrain the relevant model parameters and to investigate the transverse gluonic structure at very low Bjorken-x.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Geodetic Measurements of Slow-Slip Events Southeast of Parkfield, CA

Tremor and low-frequency earthquakes are presumed to be indicative of surrounding slow, aseismic slip that is often below geodetic detection thresholds. This study uses data from borehole seismometers and long-baseline laser strainmeters to observe both the seismic and geodetic signatures of episodic tremor and slip on the Parkfield region of the San Andreas Fault near Cholame, CA. We note the observed occurrence rates of both the tremors and co-located families of low-frequency earthquakes are not steady but instead exhibit quasiperiodic bursts of increased activity.We show that these periods of elevated seismic activity correlate with statistically significant stacked strain signals consisting of 44 slow-slip events. Modeled individual slow-slip events and their total summed moment, which are constrained by seismic signals and stacked strain, respectively, indicate that the individual moment magnitudes of these events range from $M_w$ 4.6–5.2.We find that the measured geodetic signal likely precedes the seismic signal by several hours, consistent with the aseismic slip preceding and driving the observed seismic tremor activity.We confirm that strike-slip faults, in addition to subduction zones, are capable of producing episodic tremor and slip.

58 GEOSCIENCES↗