Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spatial constraints”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

ClimSim: A large multi-scale dataset for hybrid physics-ML climate emulation

Modern climate projections lack adequate spatial and temporal resolution due to computational constraints. A consequence is inaccurate and imprecise predictions of critical processes such as storms. Hybrid methods that combine physics with machine learning (ML) have introduced a new generation of higher fidelity climate simulators that can sidestep Moore’s Law by outsourcing compute-hungry, short, high-resolution simulations to ML emulators. However, this hybrid ML-physics simulation approach requires domain-specific treatment and has been inaccessible to MLexperts because of lack of training data and relevant, easy-to-use workflows. Wepresent ClimSim, the largest-ever dataset designed for hybrid ML-physics research. It comprises multi-scale climate simulations, developed by a consortium of climate scientists and ML researchers. It consists of 5.7 billion pairs of multivariate input and output vectors that isolate the influence of locally-nested, high-resolution, high-fidelity physics on a host climate simulator’s macro-scale physical state. The dataset is global in coverage, spans multiple years at high sampling frequency, and is designed such that resulting emulators are compatible with downstream coupling into operational climate simulators. We implement a range of deterministic and stochastic regression baselines to highlight the ML challenges and their scoring. The data (https://huggingface.co/datasets/LEAP/ClimSim_high-res2) and code(https://leap-stc.github.io/ClimSim)arereleasedopenlytosupport the development of hybrid ML-physics and high-fidelity climate simulations for the benefit of science and society.

artificial intelligence, machine learning↗
DOE OSTI2021

Estimating basis functions in massive fields under the spatial mixed effects model

Abstract Spatial prediction is commonly achieved under the assumption of a Gaussian random field by obtaining maximum likelihood estimates of parameters, and then using the kriging equations to arrive at predicted values. For massive datasets, fixed rank kriging using the expectation–maximization algorithm for estimation has been proposed as an alternative to the usual but computationally prohibitive kriging method. The method reduces computation cost of estimation by redefining the spatial process as a linear combination of basis functions and spatial random effects. A disadvantage of this method is that it imposes constraints on the relationship between the observed locations and the knots. We develop an alternative method that utilizes the spatial mixed effects model, but allows for additional flexibility by estimating the range of the spatial dependence between the observations and the knots via an alternating expectation conditional maximization algorithm. Experiments show that our methodology improves estimation without sacrificing prediction accuracy while also minimizing the additional computational burden of extra parameter estimation. The methodology is applied to a temperature dataset archived by the United States National Climate Data Center, with improved results over previous methodology.

Pazdernik, Karl↗
DOE OSTI2022

Microscale spatial distribution and soil organic matter persistence in top and subsoil

The spatial distribution of organic substrates and microscale soil heterogeneity significantly influence organic matter (OM) persistence as constraints on OM accessibility to microorganisms. However, it is unclear how changes in OM spatial heterogeneity driven by factors such as soil depth affect the relative importance of substrate spatial distribution on OM persistence. This work evaluated the decomposition and persistence of 13 C and 15 N labeled water-extractable OM inputs over 50 days as either hotspot (i.e., pelleted in 1–2 mm-size pieces) or distributed (i.e., added as OM < 0.07 μm suspended in water) forms in topsoil (0–0.2 m) and subsoil (0.8–0.9 m) samples of an Andisol. We observed greater persistence of added C in the subsoil with distributed OM inputs relative to hotspot OM, indicated by a 17% reduction in cumulative mineralization of the added C and a 10% higher conversion to mineral-associated OM. A lower substrate availability potentially reduced mineralization due to OM dispersion throughout the soil. NanoSIMS (nanoscale secondary ion mass spectrometry) analysis identified organo-mineral associations on cross-sectioned aggregate interiors in the subsoil. On the other hand, in the topsoil, we did not observe significant differences in the persistence of OM, suggesting that the large amounts of particulate OM already present in the soil outweighed the influence of added OM spatial distribution. Here, we demonstrated under laboratory conditions that the spatial distribution of fresh OM input alone significantly affected the decomposition and persistence of OM inputs in the subsoil. On the other hand, spatial distribution seems to play a lower role in topsoils rich in particulate OM. The divergence in the influence of OM spatial distribution between the top and subsoil is likely driven by differences in soil mineralogy and OM composition.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Ensemble flow reconstruction in the atmospheric boundary layer from spatially limited measurements through latent diffusion models

Due to costs and practical constraints, field campaigns in the atmospheric boundary layer typically only measure a fraction of the atmospheric volume of interest. Machine learning techniques have previously successfully reconstructed unobserved regions of flow in canonical fluid mechanics problems and two-dimensional geophysical flows, but these techniques have not yet been demonstrated in the three-dimensional atmospheric boundary layer. Here, we conduct a numerical analogue of a field campaign with spatially limited measurements using large-eddy simulation. We pose flow reconstruction as an inpainting problem, and reconstruct realistic samples of turbulent, three-dimensional flow with the use of a latent diffusion model. The diffusion model generates physically plausible turbulent structures on larger spatial scales, even when input observations cover less than 1% of the volume. Through a combination of qualitative visualization and quantitative assessment, we demonstrate that the diffusion model generates meaningfully diverse samples when conditioned on just one observation. These samples successfully serve as initial conditions for a large-eddy simulation code. We find that diffusion models show promise and potential for other applications for other turbulent flow reconstruction problems.

17 WIND ENERGY↗
DOE OSTI2025

Dark Energy Survey Year 3 Results: Cosmological Constraints from Cluster Abundances, Weak Lensing, and Galaxy Clustering

Galaxy clusters provide a unique probe of the late-time cosmic structure and serve as a powerful independent test of the $\Lambda$CDM model. This work presents the first set of cosmological constraints derived with ~16,000 optically selected redMaPPer clusters across nearly 5,000 $\rm{deg}^2$ using DES Year 3 data sets. Our analysis leverages a consistent modeling framework for galaxy cluster cosmology and DES-Y3 joint analyses of galaxy clustering and weak lensing (3x2pt), ensuring direct comparability with the DES-Y3 3x2pt analysis. We obtain constraints of $S_8 = 0.864 \pm 0.035$ and $\Omega_{\rm{m}} = 0.265^{+0.019}_{-0.031}$ from the cluster-based data vector. We find that cluster constraints and 3x2pt constraints are consistent under the $\Lambda$CDM model with a Posterior Predictive Distribution (PPD) value of $0.53$. The consistency between clusters and 3x2pt provides a stringent test of $\Lambda$CDM across different mass and spatial scales. Jointly analyzing clusters with 3x2pt further improves cosmological constraints, yielding $S_8 = 0.811^{+0.022}_{-0.020}$ and $\Omega_{\rm{m}} = 0.294^{+0.022}_{-0.033}$, a $24\%$ improvement in the $\Omega_{\rm{m}}-S_8$ figure-of-merit over 3x2pt alone. Moreover, we find no significant deviation from the Planck CMB constraints with a probability to exceed (PTE) value of $0.6$, significantly reducing previous $S_8$ tension claims. Finally, combining DES 3x2pt, DES clusters, and Planck CMB places an upper limit on the sum of neutrino masses of $\sum m_\nu < 0.26$ eV at 95% confidence under the $\Lambda$CDM model. These results establish optically selected clusters as a key cosmological probe and pave the way for cluster-based analyses in upcoming Stage-IV surveys such as LSST, Euclid, and Roman.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Self‐Consistent Convolutional Density Functional Approximations: Application to Adsorption at Metal Surfaces

The exchange-correlation (XC) functional in density functional theory is used to approximate multi-electron interactions. A plethora of different functionals are available, but nearly all are based on the hierarchy of inputs commonly referred to as “Jacob's ladder.” This paper introduces an approach to construct XC functionals with inputs from convolutions of arbitrary kernels with the electron density, providing a route to move beyond Jacob's ladder. We derive the variational derivative of these functionals, showing consistency with the generalized gradient approximation (GGA), and provide equations for variational derivatives based on multipole features from convolutional kernels. A proof-of-concept functional, PBEq, which generalizes the PBEα framework with mathematical equation being a spatially-resolved function of the monopole of the electron density, is presented and implemented. It allows a single functional to use different GGAs at different spatial points in a system, while obeying PBE constraints. Analysis of the results underlines the importance of error cancellation and the XC potential in data-driven functional design. After testing on small molecules, bulk metals, and surface catalysts, the results indicate that this approach is a promising route to simultaneously optimize multiple properties of interest.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

The Milky Way, coming into focus: Precision astrometry probes its evolution and its dark matter

The growing trove of precision astrometric observations from the Gaia space telescope and other surveys is revealing the structure and dynamics of the Milky Way in ever more exquisite detail. We summarize the current status of our understanding of the structure and the characteristics of the Milky Way, and we review the emerging picture: the Milky Way is evolving through interactions with the massive satellite galaxies that stud its volume, with evidence pointing to a cataclysmic past. It is also woven with stellar streams, and observations of streams, satellites, and field stars offer new constraints on its dark matter, both on its spatial distribution and its fundamental nature. The recent years have brought much focus to the study of dwarf galaxies found within our Galaxy’s halo and their internal matter distributions. In this review, we focus on the predictions of the cold dark matter paradigm at small mass scales through precision astrometric measurements, and we summarize the modern consensus on the extent to which small-scale probes are consistent with this paradigm. We note the discovery prospects of these studies, and also how they intertwine with probes of the dynamics and evolution of the Milky Way in various and distinct ways.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Magnetic structure of Fe 16 N 2 determined by polarized neutron diffraction on thin-film samples

Strain-free Fe 16 N 2 thin films were obtained on MgO (001) single-crystal substrates with Cr seed layers. Additionally, x-ray diffraction and transmission electron microscopy confirmed that the film is textured in all three spatial directions with crystallographic orientation determined by epitaxial constraints from the substrate. Average grain size was estimated to be in the range of tens of nanometers. Ten pieces of 40-nm-thick Fe 16 N 2 thin-film samples were stacked together for half-polarized neutron diffraction study at room temperature. A total of six reflections were usable, from which the magnetic structure of Fe 16 N 2 was determined and compared to existing models. Contrary to most electronic structure calculations, the present results suggest that the magnetic moment of Fe atoms at the 8h site is large and larger than that of the 4d site. Our results support the giant magnetization of Fe 16 N 2 .

36 MATERIALS SCIENCE↗
DOE OSTI2025

Serpentine Magnet Designs for the Interaction Region of the Electron-Ion Collider (EIC)

The Electron-Ion Collider (EIC), hosted by Brookhaven National Laboratory, is designed to deliver a peak luminosity of 1 × 10 34 cm −2 sec −1 . The interaction region (IR) of the EIC imposes several constraints in terms of field quality, aperture, and spatial layout, which necessitates the development of several unique superconducting serpentine direct wind magnets. These magnets are constructed using either a single strand or a small-diameter 6-around-1 NbTi cable, presenting unique challenges for design and optimization. This paper introduces a new computational code specifically developed to streamline and integrate the design process for these magnets, enabling faster design iterations while addressing their complex requirements. Here, in this paper, we first introduce the code, which builds on established electromagnetic fundamentals. The code incorporates tools for optimizing winding patterns and for correcting magnetic multipoles; additionally, it interfaces with established magnet design software. We also present the design of several serpentine magnets for the EIC IR, demonstrating the code’s capability to deliver precise and efficient solutions. These designs highlight the code’s ability to accelerate the development cycle, ensuring the serpentine magnets meet the demanding specifications of the EIC project.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Standardizing Dainotti-correlated gamma-ray bursts, and using them with standardized Amati-correlated gamma-ray bursts to constrain cosmological model parameters

ABSTRACT We show that each of the three Dainotti-correlated gamma-ray burst (GRB) data sets recently compiled by Wang et al. and Hu et al., that together probe the redshift range 0.35 ≤ z ≤ 5.91, obey cosmological-model-independent Dainotti correlations and so are standardizable. We use these GRB data in conjunction with the best currently available Amati-correlated GRB data, that probe 0.3399 ≤ z ≤ 8.2, to constrain cosmological model parameters. The resulting cosmological constraints are weak, providing lower limits on the non-relativistic matter density parameter, mildly favouring non-zero spatial curvature, and largely consistent with currently accelerated cosmological expansion as well as with constraints determined from better-established data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Dark energy survey year 3 results: Cosmological constraints from cluster abundances, weak lensing, and galaxy clustering

Galaxy clusters provide a unique probe of the late-time cosmic structure and serve as a powerful independent test of the Λ⁢CDM model. This work presents the first set of cosmological constraints derived with ∼16,000 optically selected redMaPPer clusters across nearly 5000 deg 2 using DES year 3 datasets. Our analysis leverages a consistent modeling framework for galaxy cluster cosmology and DES-Y3 joint analyses of galaxy clustering and weak lensing (3 × 2⁢pt), ensuring direct comparability with the DES-Y3 3 × 2⁢pt analysis. Here, we obtain constraints of 𝑆 8 = 0.864 ± 0.035 and Ω m = $0.26⁢5^{+0.019}_{−0.031}$ from the cluster-based data vector. We find that cluster constraints and 3 × 2⁢pt constraints are consistent under the Λ⁢CDM model with a posterior predictive distribution (PPD) value of 0.53. The consistency between clusters and 3 × 2⁢pt provides a stringent test of Λ⁢CDM across different mass and spatial scales. Jointly analyzing clusters with 3 × 2⁢pt further improves cosmological constraints, yielding 𝑆 8 = $0.81⁢1^{+0.022}_{−0.020}$ and Ω m = $0.29⁢4^{+0.022}_{−0.033}$, a 24% improvement in the Ω m − 𝑆 8 figure of merit over 3 × 2⁢pt alone. Moreover, we find no significant deviation from the Planck CMB constraints with a probability to exceed (PTE) value of 0.6, significantly reducing previous 𝑆 8 tension claims. Finally, combining DES 3 × 2⁢pt, DES clusters, and Planck CMB places an upper limit on the sum of neutrino masses of ∑𝑚 𝜈 < 0.26 eV at 95% confidence under the Λ⁢CDM model. These results establish optically selected clusters as a key cosmological probe and pave the way for cluster-based analyses in upcoming stage-IV surveys such as LSST, Euclid, and Roman.

Abbott, T. M. C. [NSF’s National Optical-Infrared ↗
DOE OSTI2024

Continuously bounds-preserving discontinuous Galerkin methods for hyperbolic conservation laws

For finite element approximations of transport phenomena, it is often necessary to apply a form of limiting to ensure that the discrete solution remains well-behaved and satisfies physical constraints. However, these limiting procedures are typically performed at discrete nodal locations, which is not sufficient to ensure the robustness of the scheme when the solution must be evaluated at arbitrary locations (e.g., for adaptive mesh refinement, remapping in arbitrary Lagrangian–Eulerian solvers, overset meshes, etc.). In this work, a novel limiting approach for discontinuous Galerkin methods is presented which ensures that the solution is continuously bounds-preserving (i.e., across the entire solution polynomial) for any arbitrary choice of basis, approximation order, and mesh element type. Through a modified formulation for the constraint functionals, the proposed approach requires only the solution of a single spatial scalar minimization problem per element for which a highly efficient numerical optimization procedure is presented. Here, the efficacy of this approach is shown in numerical experiments by enforcing continuous constraints in high-order unstructured discontinuous Galerkin discretizations of hyperbolic conservation laws, ranging from scalar transport with maximum principle preserving constraints to compressible gas dynamics with positivity-preserving constraints.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Performance of explicit and IMEX MRI multirate methods on complex reactive flow problems within modern parallel adaptive structured grid frameworks

Large-scale multiphysics simulations are computationally challenging due to the coupling of multiple processes with widely disparate time scales. The advent of exascale computing systems exacerbates these challenges since these systems enable ever-increasing size and complexity. In recent years, there has been renewed interest in developing multirate methods as a means to handle the large range of time scales, as these methods may afford greater accuracy and efficiency than more traditional approaches of using implicit-explicit (IMEX) and low-order operator splitting schemes. However, to date there have been few performance studies that compare different classes of multirate integrators on complex application problems. In this work, we study the performance of several newly developed multirate infinitesimal (MRI) methods, implemented in the SUNDIALS solver package, on two reacting flow model problems built on structured mesh frameworks. The first model revisits prior work on a compressible reacting flow problem with complex chemistry that is implemented using BoxLib but where we now include comparisons between a new explicit MRI scheme with the multirate spectral deferred correction (SDC) methods in the original paper. The second problem uses the same complex chemistry as the first problem, combined with a simplified flow model, but runs at a large spatial scale where explicit methods become infeasible due to stability constraints. Two recently developed IMEX MRI multirate methods are tested. These methods rely on advanced features of the AMReX framework on which the model is built, such as multilevel grids and multilevel preconditioners. The results from these two problems show that MRI multirate methods can offer significant performance benefits on complex multiphysics application problems and that these methods may be combined with advanced spatial discretization to compound the advantages of both.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

A scalable planning framework of energy storage systems under frequency dynamics constraints

As the penetration of renewables increases in power systems, the declining system inertia can cause frequency stability issues. Battery energy storage systems (BESSs) respond fast and therefore can relieve the low inertia difficulty but need to be appropriately sized considering the associated cost. This paper presents a novel stochastic optimization model for economically planning BESS capacity while considering the spatial–temporal correlation of wind generation and generator outages under frequency stability constraints, which include the rate-of-change of frequency (RoCoF), frequency nadir (FN), and quasi-steady-state (QSS) frequency. A set of new FN constraints that can be easily linearized is developed. To account for renewable uncertainties, a realistic uncertainty modeling approach, Random Field, is adopted to generate wind generation scenarios by considering both spatial and temporal evolutions of wind speed profiles. The ESS sizing is formulated as a mixed-integer linear programming problem and solved by using a scalable decomposition-and-coordination approach, Surrogate Absolute Value Lagrangian Relaxation (SAVLR). To further improve the scalability and reduce computational burdens, a rolling-horizon-based update is developed and incorporated into SAVLR for providing a practical solution to the long-term planning of very large-scale power systems. Finally, a modified IEEE 118-bus system and the Polish system are used to validate the effectiveness and scalability of the model and solution methodology.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2020

Cosmological constraints from higher redshift gamma-ray burst, H ii starburst galaxy, and quasar (and other) data

ABSTRACT We use higher redshift gamma-ray burst (GRB), H ii starburst galaxy (H iiG), and quasar angular size (QSO-AS) measurements to constrain six spatially flat and non-flat cosmological models. These three sets of cosmological constraints are mutually consistent. Cosmological constraints from a joint analysis of these data sets are largely consistent with currently accelerating cosmological expansion and with cosmological constraints derived from a combined analysis of Hubble parameter (H(z)) and baryon acoustic oscillation (BAO, with Planck-determined baryonic matter density) measurements. A joint analysis of the H(z) + BAO + QSO-AS + H iiG + GRB data provides fairly model-independent determinations of the non-relativistic matter density parameter $\Omega _{\rm m_0}=0.313\pm 0.013$ and the Hubble constant $H_0=69.3\pm 1.2\, \rm {km \, s^{-1} \, Mpc^{-1}}$. These data are consistent with the dark energy being a cosmological constant and with spatial hypersurfaces being flat, but they do not rule out mild dark energy dynamics or a little spatial curvature. We also investigate the effect of including quasar flux measurements in the mix and find no novel conclusions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Magnetic nulls in interacting dipolar fields

The prominence of nulls in reconnection theory is due to the expected singular current density and the indeterminacy of field lines at a magnetic null. Electron inertia changes the implications of both features. Magnetic field lines are distinguishable only when their distance of closest approach exceeds a distance Δ d . Electron inertia ensures Δ d ≳ c/ω pe . The lines that lie within a magnetic flux tube of radius Δ d at the place where the field strength B is strongest are fundamentally indistinguishable. If the tube, somewhere along its length, encloses a point where B = 0 vanishes, then distinguishable lines come no closer to the null than ≈ (a 2 c/ω pe ) 1/3 , where a is a characteristic spatial scale of the magnetic field. Here, the behaviour of the magnetic field lines in the presence of nulls is studied for a dipole embedded in a spatially constant magnetic field. In addition to the implications of distinguishability, a constraint on the current density at a null is obtained, and the time required for thin current sheets to arise is derived.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Quantum space, quantum time, and relativistic quantum mechanics

We treat space and time as bona fide quantum degrees of freedom on an equal footing in Hilbert space. Motivated by considerations in quantum gravity, we focus on a paradigm dealing with linear, first-order Hamiltonian and momentum constraints that lead to emergent features of temporal and spatial translations. Unlike the conventional treatment, we show that Klein-Gordon and Dirac equations in relativistic quantum mechanics can be unified in our paradigm by applying relativistic dispersion relations to eigenvalues rather than treating them as operator-valued equations. With time and space being treated on an equal footing in Hilbert space, we show symmetry transformations to be implemented by unitary basis changes in Hilbert space, giving them a stronger quantum mechanical footing. Global symmetries, such as Lorentz transformations, modify the decomposition of Hilbert space; and local symmetries, such as U(1) gauge symmetry are diagonal in coordinate basis and do not alter the decomposition of Hilbert space. Here, we briefly discuss extensions of this paradigm to quantum field theory and quantum gravity.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Climate vs Energy Security: Quantifying the Trade-offs of BECCS Deployment and Overcoming Opportunity Costs on Set-Aside Land

Bioenergy with carbon capture and storage (BECCS) sits at the nexus of the climate and energy security. We evaluated trade-offs between scenarios that support climate stabilization (negative emissions and net climate benefit) or energy security (ethanol production). Our spatially explicit model indicates that the foregone climate benefit from abandoned cropland (opportunity cost) increased carbon emissions per unit of energy produced by 14–36%, making geologic carbon capture and storage necessary to achieve negative emissions from any given energy crop. The toll of opportunity costs on the climate benefit of BECCS from set-aside land was offset through the spatial allocation of crops based on their individual biophysical constraints. Dedicated energy crops consistently outperformed mixed grasslands. We estimate that BECCS allocation to land enrolled in the Conservation Reserve Program (CRP) could capture up to 9 Tg C year –1 from the atmosphere, deliver up to 16 Tg CE year –1 in emissions savings, and meet up to 10% of the US energy statutory targets, but contributions varied substantially as the priority shifted from climate stabilization to energy provision. Furthermore, our results indicate a significant potential to integrate energy security targets into sustainable pathways to climate stabilization but underpin the trade-offs of divergent policy-driven agendas.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗