Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “constrained dynamics”

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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Comparison of Conventional and Constrained Variational Methods for Computing Large-Scale Budgets and Forcing Fields

Analyses of atmospheric heat and moisture budgets serve as an effective tool to study convective characteristics over a region and to provide large-scale forcing fields for various modeling applications. This paper examines two popular methods for computing large-scale atmospheric budgets: the conventional budget method (CBM) using objectively gridded analyses based primarily on radiosonde data and the constrained variational analysis (CVA) approach which supplements vertical profiles of atmospheric fields with measurements at the top of the atmosphere and at the surface to conserve mass, water, energy, and momentum. Successful budget computations are dependent on accurate sampling and analyses of the thermodynamic state of the atmosphere and the divergence field associated with convection and the large-scale circulation that influences it. Utilizing analyses generated from data taken during Dynamics of the Madden-Julian Oscillation (DYNAMO) field campaign conducted over the central Indian Ocean from October to December 2011, we evaluate the merits of these budget approaches and examine their limitations. While many of the shortcomings of the CBM, in particular effects of sampling errors in sounding data, are effectively minimized with CVA, accurate large-scale diagnostics in CVA are dependent on reliable background fields and rainfall constraints. Finally, for the DYNAMO analyses examined, the operational model fields used as the CVA background state provided wind fields that accurately resolved the vertical structure of convection in the vicinity of Gan Island. However, biases in the model thermodynamic fields were somewhat amplified in CVA resulting in a convective environment much weaker than observed.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Cosmological constraints from a joint DESI DR1 Full-Shape and DR2 BAO

We present a cosmological analysis combining full-shape (FS) clustering measurements from the Dark Energy Spectroscopic Instrument (DESI) DR1 with baryon acoustic oscillation (BAO) measurements from DESI DR2. To achieve a robust combination that accounts for the correlation between the two data releases, we employ the ShapeFit compression method and estimate the joint covariance using EZmocks. This compressed approach inherently mitigates the prior volume effects that have previously dominated Bayesian constraints from DESI data with minimal external priors. Consequently, we obtain — for the first time within a Bayesian framework — reliable DESI-only constraints on extensions to ΛCDM using only a Big Bang Nucleosynthesis prior on the baryon density and a wide prior on the spectral index. In flat ΛCDM, we find Ω m = 0.3035 ± 0.0085, h = 0.6876 ± 0.0059, and σ 8 = 0.822 ± 0.034. For the w 0 w a CDM dynamical dark energy model, we measure w 0 = -0.49 ± 0.25 and w a = -1.52 ± 0.77, improving constraints by ∼ 30% relative to the analogous DR1 measurement and reducing the discrepancy with ΛCDM to 1.4σ when compared to BAO only analyses. We also report competitive limits on the sum of neutrino masses and spatial curvature. This work demonstrates that the ShapeFit compression provides a prior-robust and computationally efficient pathway to constrain beyond-ΛCDM physics with large-scale structure.

baryon acoustic oscillations↗
DOE OSTI2024

Exploring the Frontiers of Energy Efficiency using Power Management at System Scale

In the face of surging power demands for exascale HPC systems, this work tackles the critical challenge of understanding the impact of software-driven power management techniques like Dynamic Voltage and Frequency Scaling (DVFS) and Power Capping. These techniques have been actively developed over the past few decades. By combining insights from GPU benchmarking to understand application power profiles, we present a telemetry data-driven approach for deriving energy savings projections. This approach has been demonstrably applied to the Frontier supercomputer at scale. Our findings based on three months of telemetry data indicate that, for certain resource-constrained jobs, significant energy savings (up to 8.5%) can be achieved without compromising performance. This translates to a substantial cost reduction, equivalent to 1438 MWh of energy saved. The key contribution of this work lies in the methodology for establishing an upper limit for these best-case scenarios and its successful application. This work enables HPC professionals to optimize the power-performance trade-off within constrained power budgets, not only for the exascale era but also beyond.

Karimi, Ahmad Maroof↗
DOE OSTI2021

Combining visible shadowgraph imaging with x-ray imaging for dynamic shock experiments

Interrogating ejecta particles launched from target materials that are undergoing dynamic shock can be done with both x-ray imaging and visible shadowgraph imaging. Our dynamic testing must be done inside a containment vessel with limited access ports available. We designed an imaging system to relay both types of imaging systems through a single port using the same optical relay and then splitting the images onto three separate high-speed imaging cameras outside the containment vessel. X-ray imaging provides ejecta density measurements. Shadowgraph imaging that is done at two wavelengths (blue and red) constrains ejecta particle size distributions and provides areal density measurements of the ejecta cloud. The ejecta particles are positioned 225 mm before the x-ray scintillators; this arrangement permits a folded mirror system to allow the shadowgraph data to bypass the x-ray scintillators. This configuration results in spatial separations between the intermediate image planes of the x-ray and shadowgraph images along the optical axis. At the position of the x-ray intermediate image plane, mirrors are positioned such that the shadowgraph images are kicked out and their images are sent on to different cameras. Positioning of the large doublet relay lenses keeps shrapnel from impacting the vessel containment windows.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Knowledge-guided machine learning can improve carbon cycle quantification in agroecosystems

Abstract Accurate and cost-effective quantification of the carbon cycle for agroecosystems at decision-relevant scales is critical to mitigating climate change and ensuring sustainable food production. However, conventional process-based or data-driven modeling approaches alone have large prediction uncertainties due to the complex biogeochemical processes to model and the lack of observations to constrain many key state and flux variables. Here we propose a Knowledge-Guided Machine Learning (KGML) framework that addresses the above challenges by integrating knowledge embedded in a process-based model, high-resolution remote sensing observations, and machine learning (ML) techniques. Using the U.S. Corn Belt as a testbed, we demonstrate that KGML can outperform conventional process-based and black-box ML models in quantifying carbon cycle dynamics. Our high-resolution approach quantitatively reveals 86% more spatial detail of soil organic carbon changes than conventional coarse-resolution approaches. Moreover, we outline a protocol for improving KGML via various paths, which can be generalized to develop hybrid models to better predict complex earth system dynamics.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Determining the diffusivity for light quarks from experiment

Charge balance functions reflect the evolution of charged pair correlations throughout the stages of pair production, dynamical diffusion, and hadronization in heavy-ion collisions. In this work, microscopic modeling of these correlations in the full collision volume shows that the balance functions are sensitive to the diffusivity of light quarks when studied as functions of relative azimuthal angle. By restricting our analysis to K + K – and $p\bar{p}$ pairs, we find that the diffusivity of light quarks, a fundamental property not currently well understood, can be constrained by experimental measurement.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Machine learning for nonintrusive model order reduction of the parametric inviscid transonic flow past an airfoil

Fluid flow in the transonic regime finds relevance in aerospace engineering, particularly in the design of commercial air transportation vehicles. Computational fluid dynamics models of transonic flow for aerospace applications are computationally expensive to solve because of the high degrees of freedom as well as the coupled nature of the conservation laws. While these issues pose a bottleneck for the use of such models in aerospace design, computational costs can be significantly minimized by constructing special, structure-preserving surrogate models called reduced-order models. In this work, we propose a machine learning method to construct reduced-order models via deep neural networks and we demonstrate its ability to preserve accuracy with a significantly lower computational cost. In addition, our machine learning methodology is physics-informed and constrained through the utilization of an interpretable encoding by way of proper orthogonal decomposition. Application to the inviscid transonic flow past the RAE2822 airfoil under varying freestream Mach numbers and angles of attack, as well as airfoil shape parameters with a deforming mesh, shows that the proposed approach adapts to high-dimensional parameter variation well. Notably, the proposed framework precludes the knowledge of numerical operators utilized in the data generation phase, thereby demonstrating its potential utility in the fast exploration of design space for diverse engineering applications. Comparison against a projection-based nonintrusive model order reduction method demonstrates that the proposed approach produces comparable accuracy and yet is orders of magnitude computationally cheap to evaluate, despite being agnostic to the physics of the problem.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

A kinetic model for simulating non-equilibrium mass transport in oxides applied to hematite growth under irradiation

We propose an approach to simulating the dynamic evolution and transport of charged point defects within and through the oxide scales that form during corrosion. The method follows the cluster dynamics formalism widely adopted for radiation damage in solids, which can apply in both quasi-static and far from equilibrium conditions, such as in radiation environments. By treating each charged defect as a cluster of an atomic defect and associated unit charges, the proposed model flexibly allows charge state transitions and shifts in the Fermi level by absorption and emission of charge carriers from point defects in a reaction network with rates constrained by local equilibrium. Applying this model to hematite predicts changes in self-diffusion and oxidation kinetics in irradiation environments, surprisingly showing reduced oxidation rates in many conditions, despite enhanced self-diffusion. In conclusion, the origin of this effect lies in a change in Fermi level induced by excess vacancies formed under irradiation, which in turn suppresses the transport of cation interstitials which facilitate hematite growth.

36 MATERIALS SCIENCE↗
DOE OSTI2023

FPGA Acceleration of GCN in Light of the Symmetry of Graph Adjacency Matrix

Graph Convolutional Neural Networks (GCNs) are widely used to process large-scale graph data. Different from deep neural networks (DNNs), GCNs are sparse, irregular, and unstructured, posing unique challenges to hardware acceleration with regular processing elements (PEs). In particular, the adjacency matrix of a GCN is extremely sparse, leading to frequent but irregular memory access, low spatial/temporal data locality and poor data reuse. Furthermore, a realistic graph usually consists of unstructured data (e.g., unbalanced distributions), creating significantly different processing times and imbalanced workload for each node in GCN acceleration. To overcome these challenges, we propose an end-to-end hardware-software co-design to accelerate GCNs on resource-constrained FPGAs with the features including: (1) A custom dataflow that leverages symmetry along the diagonal of the adjacency matrix to accelerate feature aggregation for undirected graphs. We utilize either the upper or the lower triangular matrix of the adjacency matrix to perform aggregation in GCN to improve data reuse. (2) Unified compute cores for both aggregation and transform phases, with full support to the symmetry-based dataflow. These cores can be dynamically reconfigured to the systolic mode for transformation or as individual accumulators for aggregation in GCN processing. (3) Preprocessing of the graph in software to rearrange the edges and features to match the custom dataflow. This step improves the regularity in memory access and data reuse in the aggregation phase. Moreover, we quantize the GCN precision from FP32 to INT8 to reduce the memory footprint without losing the inference accuracy. We implement our accelerator design in Intel Stratix10 MX FPGA board with HBM2, and demonstrate 1.3x-110.5x improvement in end-to-end GCN latency as compared to the state-of the-art FPGA implementations, on the graph datasets of Cora, Pubmed, Citeseer and Reddit.

Nair, Gopikrishnan R.↗
DOE OSTI2024

MRS Modelling in OpenFAST

Multi-rotor wind turbines, where more than one rotor is connected to a single support structure, are increasing in popularity due to their balance of power generation with component size, weight, and cost. Multi-rotor turbines are also common in marine energy, where deployment sites tend to be depth-constrained, limiting the maximum size of a rotor. To support the industry and research communities in designing and modeling multi-rotor turbines, the open-source turbine modeling tool OpenFAST, developed and supported by the National Renewable Energy Laboratory, is being extended to model multi-rotor wind and marine turbines. This presentation gives an overview of the current state of OpenFAST, planned changes to enable coupled aero-hydro-servo-elastic dynamic modeling of multi-rotor turbines, and examples of the types of turbines that OpenFAST will support.

marine energy↗
DOE OSTI2023

Comparison of two pressure–temperature equilibration methods

We compare and contrast the traditionally used method of solving the pressure–temperature equilibration problem in hydrodynamics, where specific internal energy and density are considered independent variables, with a different method where pressure and temperature are independent variables. With the goal of examining the robustness of the two methods as the number of components increases, we examine 2-, 4-, 6-, and 8-component systems. After equilibrating more than 10 4 initial conditions for each system using both methods, we demonstrate that the latter method constrains the search space by lowering its dimensionality and forces a better initial guess, resulting in a higher probability of convergence to solution with fewer, cheaper iterations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Plastic homogeneity in nanoscale heterostructured binary and multicomponent metallic eutectics: An overview

Heterostructured materials comprised of relatively soft/hard disparate phases typically exhibit composite strengthening but lack plastic deformability at ambient temperatures. However, heterostructured systems comprised of nanoscale phases can simultaneously enhance yield strength and strain hardening, thereby promoting uniform distribution of plastic flow. In this review, the atomic-scale deformation mechanisms in model systems of eutectic alloys, Al-Al 2 Cu and Al-Si, refined to nanoscales via laser rapid solidification are discussed, and compared with literature on multi-component (high entropy) eutectics such as Ni-Al-Fe-based with Cr and/or Co additions. The nano-lamellar Al-Al 2 Cu structures exhibit unit defect mechanisms not reported in monolithic Al 2 Cu intermetallic: localized shear on {0 1 1} and shear-induced faults on {1 2 1} planes, constrained by closely-spaced dislocation arrays in Al confined by Al/Al 2 Cu interfaces. The unexpected plasticity mechanisms are enabled by slip continuity in nanoscale Al-Al 2 Cu eutectics associated with the orientation relationship and interface habit planes. In nano-fibrous Al-Si eutectic, tensile ductility at strength approaching 600 MPa is observed resulting from dislocation plasticity in the nano-Al channels and cracking in Si nanofibers. Molecular dynamics simulations show that Al dislocations easily cross-slip (screw) or climb (edge) along Al-Si interfaces, making slip transmission difficult. The propagation of nano-cracks is suppressed by surrounding strain hardening Al, retaining good ductility of the sample, in spite of lack of direct slip transmission. Finally, the critical unit mechanisms of slip transmission and interface-enabled plasticity observed in nanoscale eutectic binary systems can also explain the strength-ductility relationship in multi-component eutectics and homogeneously distributed plastic flow with increasing microstructural heterogeneity.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Nonlinearity of the post-spinel transition and its expression in slabs and plumes worldwide

Phase transitions in the mantle control its internal dynamics and structure. The post-spinel transition marks the upper–lower mantle boundary, where ringwoodite dissociates into bridgmanite plus ferropericlase, and its Clapeyron slope regulates mantle flow across it. This interaction has previously been assumed to have no lateral spatial variations, based on the assumption of a linear post-spinel boundary in pressure and temperature. Here we present laser-heated diamond anvil cell experiments with synchrotron X-ray diffraction to better constrain this boundary, especially at higher temperatures. Combining our data with results from the literature, and using a global analysis based on machine learning, we find a pronounced nonlinearity in the post-spinel boundary, with its slope ranging from –4 MPa/K at 2100 K, to –2 MPa/K at 1950 K, and to 0 MPa/K at 1600 K. Changes in temperature over time and space can therefore cause the post-spinel transition to have variable effects on mantle convection and the movement of subducting slabs and upwelling plumes.

58 GEOSCIENCES↗
DOE OSTI2025

Thermal-Fluid and Thermal-Structural Response of the T-Tube Modular Divertor to Spatiotemporally Varying Heat Loads

Tungsten (W) is the leading candidate for divertor target plates because of its high melting point (>3000°C), thermal conductivity, and ultimate tensile stress. While W and its alloys are the only solid materials that can survive the high heat fluxes incident on the divertor, W’s low-ductility high ductile-to-brittle transition temperature of ~600°C and relatively low recrystallization temperature (RT) of ~1300°C pose structural (among other) challenges. The objective of this work is to estimate the thermal-fluid and thermal-structural performance of the helium (He)-cooled T-tube divertor, which was originally developed by the Advanced Reactor Innovation and Evaluation Study (ARIES) using numerical simulations. Here, predictions of temperature distributions across the plasma-facing structural component and surface pressures from computational fluid dynamics simulations are used to determine stress distributions using commercial structural finite element modeling software over a range of fusion-relevant conditions. The maximum allowable incident heat fluxes are determined based on the temperature limits imposed by the ITER elastic Structural Design Criteria for In-vessel Components (SDC-IC) and the maximum RT over a range of He mass flow rates and presented in the form of performance design charts. Our recent work found that thermal- structural criteria accounting for the low ductility of W in a finger-type modular divertor constrain the maximum incident heat fluxes to values well below the ITER specifications, and those based on considering only the RT demonstrate that integrated thermal-fluid and elastic structural performance evaluation are required for accurate assessment of divertor performance. This novel analysis of the T-tube considers how nonuniform and transient incident heat fluxes affect its thermal-fluid and thermal-structural performance, as well as the effect of volumetric heating, which can be as great as 27% of the power incident on the divertor surface. The W tile of the T-tube, with its relatively large plasma-facing area of ~15 cm 2 , will likely experience significant spatial variations in incident heat flux. This work therefore assesses whether steady-state incident heat flux profiles with a peak of 10 MW/m 2 and maximum heat flux gradients of 200 MW/m 2 per m exceed the structural limits imposed by the ITER elastic SDC-IC and the maximum RT over a range of fusion-relevant conditions. The effect of transient heat fluxes typical of plasma detachment and reattachment from the target plate due, for example, to gas injection are also evaluated

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Persistence of soil organic carbon caused by functional complexity

Soil organic carbon management has the potential to aid climate change mitigation through drawdown of atmospheric carbon dioxide. To be effective, such management must account for processes influencing carbon storage and re-emission at different space and time scales. Achieving this requires a conceptual advance in our understanding to link carbon dynamics from the scales at which processes occur to the scales at which decisions are made. Here, we propose that soil carbon persistence can be understood through the lens of decomposers as a result of functional complexity derived from the interplay between spatial and temporal variation of molecular diversity and composition. For example, co-location alone can determine whether a molecule is decomposed, with rapid changes in moisture leading to transport of organic matter and constraining the fitness of the microbial community, while greater molecular diversity may increase the metabolic demand of, and thus potentially limit, decomposition. This conceptual shift accounts for emergent behaviour of the microbial community and would enable soil carbon changes to be predicted without invoking recalcitrant carbon forms that have not been observed experimentally. Finally, functional complexity as a driver of soil carbon persistence suggests soil management should be based on constant care rather than one-time action to lock away carbon in soils.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Entanglement Structure of Non-Gaussian States and How to Measure It

Rapidly growing capabilities of quantum simulators to probe quantum many-body phenomena require new methods to characterize increasingly complex states. Here, we present a protocol that constrains quantum states using experimentally measured correlation functions. This method enables measurement of a quantum state’s entanglement structure, opening a new route to study entanglement-related phenomena. Our approach extends Gaussian state parameterizations by systematically incorporating higher-order correlations. We show the protocol’s usefulness in conjunction with current and forthcoming experimental capabilities, focusing on weakly interacting fermions as a proof of concept. Here, the lowest nontrivial expansion quantitatively predicts early time thermalization dynamics, including signaling the onset of quantum chaos indicated by the entanglement Hamiltonian.

Fermi gases↗
DOE OSTI2022

The MASSIVE Survey. XVI. The Stellar Initial Mass Function in the Center of MASSIVE Early-type Galaxies

The stellar initial mass function (IMF) is a fundamental property in the measurement of stellar masses and galaxy star formation histories. In this work, we focus on the most massive galaxies in the nearby universe $\mathrm{log}({M}_{\star }/\,{M}_{\odot })\gt 11.2$. We obtain high-quality Magellan/LDSS-3 long-slit spectroscopy with a wide wavelength coverage of 0.4–1.01 μm for 41 early-type galaxies (ETGs) in the MASSIVE survey and derive high signal-to-noise spectra within an aperture of R e /8. Using detailed stellar synthesis models, we constrain the elemental abundances and stellar IMF of each galaxy through full spectral modeling. All the ETGs in our sample have an IMF that is steeper than a Milky Way (Kroupa) IMF. The best-fit IMF mismatch parameter, α IMF = (M/L)/(M/L) MW , ranges from 1.1 to 3.1, with an average of $\langle$α IMF $\rangle$ = 1.84, suggesting that on average, the IMF is more bottom heavy than Salpeter. Comparing the estimated stellar masses with the dynamical masses, we find that most galaxies have stellar masses that are smaller than their dynamical masses within the 1σ uncertainty. We complement our sample with lower-mass galaxies from the literature and confirm that $\mathrm{log}({\alpha }_{\mathrm{IMF}})$ is positively correlated with $\mathrm{log}(\sigma )$, $\mathrm{log}({M}_{\star })$, and $\mathrm{log}({M}_{\mathrm{dyn}})$. From the combined sample, we show that the IMF in the centers of more massive ETGs is more bottom heavy. In addition, we find that $\mathrm{log}({\alpha }_{\mathrm{IMF}})$ is positively correlated with both [Mg/Fe] and the estimated total metallicity [Z/H]. We find suggestive evidence that the effective stellar surface density Σ Kroupa might be responsible for the variation of α IMF . We conclude that σ, [Mg/Fe], and [Z/H] are the primary drivers of the global stellar IMF variation.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Galaxy Clustering in the Mira-Titan Universe. I. Emulators for the Redshift Space Galaxy Correlation Function and Galaxy–Galaxy Lensing

Abstract We construct accurate emulators for the projected and redshift space galaxy correlation functions and excess surface density as measured by galaxy–galaxy lensing, based on halo occupation distribution modeling. Using the complete Mira-Titan suite of 111 N -body simulations, our emulators vary over eight cosmological parameters and include the effects of neutrino mass and dynamical dark energy. We demonstrate that our emulators are sufficiently accurate for the analysis of the Baryon Oscillation Spectroscopic Survey DR12 CMASS galaxy sample over the range 0.5 ≤ r ≤ 50 h −1 Mpc. Furthermore, we show that our emulators are capable of recovering unbiased cosmological constraints from realistic mock catalogs over the same range. Our mock catalog tests show the efficacy of combining small-scale galaxy–galaxy lensing with redshift space clustering and that we can constrain the growth rate and σ 8 to 7% and 4.5%, respectively, for a CMASS-like sample using only the measurements covered by our emulator. With the inclusion of a cosmic microwave background prior on H 0 , this reduces to a 2% measurement of the growth rate.

79 ASTRONOMY AND ASTROPHYSICS↗