Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reference models”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,261 records · Page 70

Methodology for Digital Image Correlation and Infrared Measurement of Melting Aluminum Bars

Ultimately, our experiment measures two quantities on an aluminum bar: motion (which modeling must predict) and temperature (which sets thermal boundary conditions). For motion, stereo DIC is a technique to use imaging data to provide displacements relative to a reference image down to 1/100th of a pixel. We use a calibrated infrared imaging method for accurate temperature measurements. We will be capturing simultaneous data and then registering temperature data in space to the same coordinate system as the displacement data. While we will later show that our experiments are repeatable, indicating that separate experiments for motion and temperature would provide similar data, the simultaneous and registered data removes test to test variability as a source of uncertainty for model calibration and reduces the number of time-consuming tests that must be performed.

36 MATERIALS SCIENCE↗

Weak Form Scientific Machine Learning: Test Function Construction for System Identification

Weak form Scientific Machine Learning (WSciML) is a recently developed framework for data-driven modeling and scientific discovery. It leverages the weak form of equation error residuals to provide enhanced noise robustness in system identification via convolving model equations with test functions, reformulating the problem to avoid direct differentiation of data. The performance, however, relies on wisely choosing a set of compactly supported test functions. In this work, we mathematically motivate a novel data-driven method for constructing Single-scale-Local reference functions for creating the set of test functions. Our approach numerically approximates the integration error introduced by the quadrature and identifies the support size for which the error is minimal, without requiring access to the model parameter values. Through numerical experiments across various models, noise levels, and temporal resolutions, we demonstrate that the selected supports consistently align with regions of minimal parameter estimation error. We also compare the proposed method against the strategy for constructing Multi-scale-Global (and orthogonal) test functions introduced in our prior work, demonstrating the improved computational efficiency.

FOS: Computer and information sciences↗

Canadian Oil Sands Extraction and Upgrading: A Synthesis of the Data on Energy Consumption, CO 2 Emissions, and Supply Costs

As Canadian crude bitumen production from oil sands has increased in recent decades, the nation’s oil and gas industry has become a significant contributor to national greenhouse gas emissions. Canada has developed carbon emission reduction targets to meet its Nationally Determined Contributions and Mid-Century Strategy goals. A detailed profile of energy consumption pathways in the oil sands industry is necessary to identify potential areas of improvement and to monitor progress toward meeting emissions reduction targets. Much of the existing literature for oil sands modeling provides input assumptions with different technological boundaries. For a set of oil sands extraction and upgrading technologies, this study first reviews the literature and then quantifies energy input requirements, CO 2 emissions, and operating costs for a set of consistent technological boundaries and energy units. Summary results refer to requirements and costs at the production facility, excluding transportation and blending costs. An energy system diagram of oil sands production that matches these boundaries is provided, which can be used by integrated assessment models, oil sands companies, and government ministries to evaluate the present and future energy consumption and emissions pathways of the oil sands industry.

02 PETROLEUM↗

Initial steady-state core simulation capability or thermal and pool-type molten salt reactors, coupling reactor physics, thermal-hydraulics, and evolving chemistry

This report presents the development and validation of an initial steady-state multiphysics capability for molten salt reactors (MSRs) under the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program in Fiscal Year 2025. The framework integrates neutronics, thermal-hydraulics, species transport, and thermochemistry to capture the coupled dynamics of liquid-fueled systems. Implementation and testing were performed on two representative designs: the Molten Salt Reactor Experiment (MSRE), a thermal-spectrum, channeled-flow reactor, and the Lotus Molten Salt Reactor (L-MSR), a fast-spectrum, pool-type reactor. The modeling suite employs Griffin for reactor physics and depletion, Pronghorn and SAM for thermal-hydraulics, Thermochimica for chemistry, and Saline for thermophysical properties, with benchmarking and validation carried out against historical MSRE data, experimental flow-loop measurements, and reference depletion calculations from Monte Carlo codes. The framework demonstrated the ability to reproduce key reactor behaviors including temperature feedback, reactivity losses, delayed neutron precursor transport, xenon poisoning, and redox potential evolution. The results confirm the feasibility and accuracy of the coupled models in predicting steady-state and selected transient MSR behaviors. This latter ones are used in this report as a proxy indicating that the steady-state models from which the transient starts are accurate. For MSRE, validation showed good agreement with pump start-up and natural circulation tests, while for the L-MSR, benchmarking confirmed hydraulic calibration and consistency of neutronics–thermal coupling. The tools also provided new insights into species transport, noble metal deposition, and salt solidification dynamics. On the Xenon transport front, the code is validated against the steady state Xenon poisoining measurement and showed good agreement with the experimental value. Identified areas for future work include advanced void transport modeling, three-dimensional simulations, improved alloy corrosion models, and tighter integration with high-fidelity Monte Carlo codes. These developments provide a foundation for high-fidelity MSR simulations that can support reactor design optimization, safety assessments, and long-term operational strategies.

42 - ENGINEERING↗

Electrochemical Modeling and Experimental Verification of Lithiation Gradients in Oxide Cathodes of Lithium-Ion Cells

Lithiated nickel-cobalt-manganese oxides, such as NCM523, are used in the positive electrode (cathode) of Li-ion cells. Using operando X-ray diffraction profilometry, lithiation gradients in the cathode matrix can be observed and quantified by expansion into Legendre polynomials with time-dependent weights. These weights (referred to as gradients) increase in magnitude when electric current flows through the cell, decrease during potentiostatic hold and finally relax to zero when the current is interrupted during open circuit rest. Both physics-based electrochemical models and operando X-ray experiments suggest that the time constants for gradient growth and abatement are primarily determined by ionic diffusion in the oxide particles, which in turn depends on their lithium content. In contrast, the magnitude of gradients depends mainly on the applied current. The X-ray profilometry provides a way of directly probing the formation and disappearance of Li gradients across the cathode during fast cycling, which can help to diagnose the effects of material degradation in the cells.

25 ENERGY STORAGE↗

Conservation Voltage Reduction (CVR) via Two-Timescale Control in Unbalanced Power Distribution Systems

Voltage control devices are employed in power distribution systems to reduce the power consumption by operating the system closer to the lower acceptable voltage limits; this technique is called conservation voltage reduction (CVR). The different modes of operation for system’s legacy devices (with binary control) and new devices (e.g. smart inverters with continuous control) coupled with variable photovoltaic (PV) generation results in voltage fluctuations which makes it challenging to achieve CVR objective. This paper presents a two-timescale control of feeder’s voltage control devices to achieve CVR that includes (1)a centralized controller operating in a slower time-scale to coordinate voltage control devices across the feeder and (2) local controllers operating in a faster timescale to mitigate voltage fluctuations due to PV variability. The centralized controller utilizes a three-phase optimal power flow model to obtain the decision variables for both legacy devices and smart inverters. The local controllers operate smart inverters to minimize voltage fluctuations and restore nodal voltages to their reference values by adjusting the reactive power support. The proposed approach is validated using the IEEE-123 bus (medium-size) and R3-12.47-2 (large-size) feeders. It is demonstrated that the proposed approach is effective in achieving the CVR objective for unbalanced distribution systems.

optimal control, distributed generation, fluctuati↗

Allometric relationships and trade‐offs in 11 common M editerranean‐climate grasses

Abstract Biomass allocation in plants is the foundation for understanding dynamics in ecosystem carbon balance, species competition, and plant–environment interactions. However, existing work on plant allometry has mainly focused on trees, with fewer studies having developed allometric equations for grasses. Grasses with different life histories can vary in their carbon investment by prioritizing the growth of specific organs to survive, outcompete co‐occurring plants, and ensure population persistence. Further, because grasses are important fuels for wildfire, the lack of grass allocation data adds uncertainty to process‐based models that relate plant physiology to wildfire dynamics. To fill this gap, we conducted a greenhouse experiment with 11 common California grasses varying in photosynthetic pathway and growth form. We measured plant sizes and harvested above‐ and belowground biomass throughout the life cycle of annual species, while for the establishment stage of perennial grasses to quantify allometric relationships for leaf, stem, and root biomass, as well as plant height and canopy area. We used basal diameter as a reference measure of plant size. Overall, basal diameter is the best predictor for leaf and stem biomass, height, and canopy area. Including height as another predictor can improve model accuracy in predicting leaf and stem biomass and canopy area. Fine root biomass is a function of leaf biomass alone. Species vary in their allometric relationships, with most variation occurring for plant height, canopy area, and stem biomass. We further explored potential trade‐offs in biomass allocation across species between leaf and fine root, leaf and stem, and allocation to reproduction. Consistent with our expectation, we found that fast‐growing plants allocated a greater fraction to reproduction. Additionally, plant height and specific leaf area negatively influenced the leaf‐to‐stem ratio. However, contrary to our hypothesis, there were no differences in root‐to‐leaf ratio between perennial and annual or C 4 and C 3 plants. Our study provides species‐specific and functional‐type‐specific allometry equations for both above‐ and belowground organs of 11 common California grass species, enabling nondestructive biomass assessment in California grasslands. These allometric relationships and trade‐offs in carbon allocation across species can improve ecosystem model predictions of grassland species interactions and environmental responses through differences in morphology.

54 ENVIRONMENTAL SCIENCES↗

Integrating chromatin conformation information in a self-supervised learning model improves metagenome binning

Metagenome binning is a key step, downstream of metagenome assembly, to group scaffolds by their genome of origin. Although accurate binning has been achieved on datasets containing multiple samples from the same community, the completeness of binning is often low in datasets with a small number of samples due to a lack of robust species co-abundance information. In this study, we exploited the chromatin conformation information obtained from Hi-C sequencing and developed a new reference-independent algorithm, Metagenome Binning with Abundance and Tetra-nucleotide frequencies—Long Range (metaBAT-LR), to improve the binning completeness of these datasets. This self-supervised algorithm builds a model from a set of high-quality genome bins to predict scaffold pairs that are likely to be derived from the same genome. Then, it applies these predictions to merge incomplete genome bins, as well as recruit unbinned scaffolds. We validated metaBAT-LR’s ability to bin-merge and recruit scaffolds on both synthetic and real-world metagenome datasets of varying complexity. Benchmarking against similar software tools suggests that metaBAT-LR uncovers unique bins that were missed by all other methods.

59 BASIC BIOLOGICAL SCIENCES↗

Future impacts of ozone driven damages on agricultural systems

Current ozone (O 3 ) concentration levels entail significant damages in crop yields around the world. The reaction of the emitted precursors (mostly methane and nitrogen oxides) with solar radiation contribute to O 3 levels that exceed established thresholds for crop damage. This paper shows current and projected (up to 2080) relative yield losses (RYLs) driven by O 3 exposure for different crops and the associated economic damages applying dynamic crop production and prices that are calculated per region and period. We adjust future crop yields in the Global Change Assessment Model (GCAM) to reflect the RYLs and analyze the effects on agricultural markets. We find that the changes (generally reductions) in O 3 precursor emissions in a reference scenario would reduce the agricultural damages, compared to present, for most of the regions, with a few exceptions including India, where higher future O3 concentrations have large negative impacts on crop yields. The annual economic impact of O3 driven losses from 2010-2080 are, in billion US dollars at 2015 prices ($B), 5.0-6.0, 9.8-18.8, 6.7-10.6 and 10.4-12.5 for corn, soybeans, rice and wheat, respectively, with the large losses for wheat and soybeans driven by their comparatively high responses to O 3 . When O 3 effects are explicitly modelled as exogenous yield shocks in future periods, there is a direct impact in future agricultural markets. Therefore, the aggregated net present value (NPV) of crop production would be reduced around by $90.8B at a global level. However, these changes are not distributed evenly across regions, and the net present market value of the crops would increase by up to $118.2B (India) or decrease by up to $59.2B (China).

54 ENVIRONMENTAL SCIENCES↗

BEYONDPLANCK IV. On end-to-end simulations in CMB analysis — Bayesian versus frequentist statistics

End-to-end simulations play a key role in the analysis of any high-sensitivity cosmic microwave background (CMB) experiment, providing high-fidelity systematic error propagation capabilities that are unmatched by any other means. In this paper, we address an important issue regarding such simulations, namely, how to define the inputs in terms of sky model and instrument parameters. These may either be taken as a constrained realization derived from the data or as a random realization independent from the data. We refer to these as posterior and prior simulations, respectively. We show that the two options lead to significantly different correlation structures, as prior simulations (contrary to posterior simulations) effectively include cosmic variance, but they exclude realization-specific correlations from non-linear degeneracies. Consequently, they quantify fundamentally different types of uncertainties. We argue that as a result, they also have different and complementary scientific uses, even if this dichotomy is not absolute. In particular, posterior simulations are in general more convenient for parameter estimation studies, while prior simulations are generally more convenient for model testing. Before BEYONDPLANCK, most pipelines used a mix of constrained and random inputs and applied the same hybrid simulations for all applications, even though the statistical justification for this is not always evident. BEYONDPLANCK represents the first end-to-end CMB simulation framework that is able to generate both types of simulations and these new capabilities have brought this topic to the forefront. The BEYONDPLANCK posterior simulations and their uses are described extensively in a suite of companion papers. In this work, we consider one important applications of the corresponding prior simulations, namely, code validation. Specifically, we generated a set of one-year LFI 30 GHz prior simulations with known inputs and we used these to validate the core low-level BEYONDPLANCK algorithms dealing with gain estimation, correlated noise estimation, and mapmaking.

79 ASTRONOMY AND ASTROPHYSICS↗

Deep learning for time series forecasting: a survey of recent advances

Time series forecasting plays a critical role in numerous real-world applications, such as finance, healthcare, transportation, and scientific computing. In recent years, deep learning has become a powerful tool for modeling complex temporal patterns and improving forecasting accuracy. This survey provides an overview of recent deep learning approaches for time series forecasting, involving various architectures including RNNs, CNNs, GNNs, transformers, large language models, MLP-based models, and diffusion models. We first identify key challenges in the field, such as temporal dependency, efficiency, and cross-variable dependency, which drive the development of forecasting techniques. Then, the general advantages and limitations of each architecture are discussed to contextualize their adaptation in time series forecasting. Furthermore, we highlight promising design trends like multi-scale modeling, decomposition, and frequency-domain techniques, which are shaping the future of the field. This paper serves as a compact reference for researchers and practitioners seeking to understand the current landscape and future trajectory of deep learning in time series forecasting.

97 MATHEMATICS AND COMPUTING↗

Contributions From Cloud Morphological Changes to the Interannual Shortwave Cloud Feedback Based on MODIS and ISCCP Satellite Observations

The surface temperature-mediated change in cloud properties, referred to as the cloud feedback, continues to dominate the uncertainty in climate projections. A larger number of contemporary global climate models (GCMs) project a higher degree of warming than the previous generation of GCMs. This greater projected warming has been attributed to a less negative cloud feedback in the Southern Ocean. Here, we apply a novel “double decomposition method” that employs the “cloud radiative kernel” and “cloud regime” concepts, to two data sets of satellite observations to decompose the interannual cloud feedback into contributions arising from changes within and shifts between cloud morphologies. Our results show that contributions from the latter to the cloud feedback are large for certain regimes. We then focus on interpreting how both changes within and between cloud morphologies impact the shortwave cloud optical depth feedback over the Southern Ocean in light of additional observations. Results from the former cloud morphological changes reveal the importance of the wind response to warming increases low- and mid-level cloud optical thickness in the same region. Results from the latter cloud morphological changes reveal that a general shift from thick storm-track clouds to thinner oceanic low-level clouds contributes to a positive feedback over the Southern Ocean that is offset by shifts from thinner broken clouds to thicker mid- and low-level clouds. Our novel analysis can be applied to evaluate GCMs and potentially diagnose shortcomings pertaining to their physical parameterizations of particular cloud morphologies.

58 GEOSCIENCES↗

Neuromorphic intermediate representation: A unified instruction set for interoperable brain-inspired computing

Abstract Spiking neural networks and neuromorphic hardware platforms that simulate neuronal dynamics are getting wide attention and are being applied to many relevant problems using Machine Learning. Despite a well-established mathematical foundation for neural dynamics, there exists numerous software and hardware solutions and stacks whose variability makes it difficult to reproduce findings. Here, we establish a common reference frame for computations in digital neuromorphic systems, titled Neuromorphic Intermediate Representation (NIR). NIR defines a set of computational and composable model primitives as hybrid systems combining continuous-time dynamics and discrete events. By abstracting away assumptions around discretization and hardware constraints, NIR faithfully captures the computational model, while bridging differences between the evaluated implementation and the underlying mathematical formalism. NIR supports an unprecedented number of neuromorphic systems, which we demonstrate by reproducing three spiking neural network models of different complexity across 7 neuromorphic simulators and 4 digital hardware platforms. NIR decouples the development of neuromorphic hardware and software, enabling interoperability between platforms and improving accessibility to multiple neuromorphic technologies. We believe that NIR is a key next step in brain-inspired hardware-software co-evolution, enabling research towards the implementation of energy efficient computational principles of nervous systems. NIR is available atneuroir.org

Science & Technology - Other Topics↗

Time-temperature-stress superposition of PBX 9502 compressive creep data for lifetime predictions

The plastic-bonded explosive (PBX) 9502 is a viscoelastic, high-solids loaded polymer bound composite comprised of 95 weight% (wt%) TATB explosive crystals and 5 wt% FK-800 polymer binder. The uniaxial quasi-static mechanical properties have been studied and characterized extensively over 25+ years, including creep. Creep is the strain evolution over time as a result of an asymmetric load typically at temperatures above ambient. Creep measurements to failure can take excessively long periods of time depending on various factors including applied stress, temperature, and loading rate that render a single test to failure impossible. Using the principles of time-temperature-stress superposition (TTSSP), short-term creep tests at various applied stresses and temperatures can be used for long-term creep predictions. The TTSSP model shifts short-term test data using two factors, one for temperature and stress, creating a master curve at a reference test condition. The master curve can be shifted to predict the creep response at any desired test condition. PBX 9502 compression specimens were tested at three temperatures and various applied stresses. This work discusses the TTSSP analysis of PBX 9502 compressive creep data to create a single comprehensive master curve.

36 MATERIALS SCIENCE↗

Unraveling H2 chemisorption and physisorption on metal decorated graphene using quantum Monte Carlo

Molecular hydrogen has the potential to significantly reduce the use of carbon dioxide emitting energy processes. However, hydrogen gas storage is a major bottleneck for its large-scale use as current storage methods are energy intensive. Among different storage methods, physisorbing molecular hydrogen at ambient pressure and temperatures is a promising alternative—particularly in light of the advancements in tunable lightweight nanomaterials and high throughput screening methods. Nonetheless, understanding hydrogen adsorption in well-defined nanomaterials remains experimentally challenging and reference information is scarce despite the proliferation of works predicting hydrogen adsorption. We focus on Li, Na, Ca, and K, decorated graphene sheets as substrates for molecular hydrogen adsorption, and compute the most accurate adsorption energies available to date using quantum diffusion Monte Carlo (DMC). Building on our previous insights at the density functional theory (DFT) level, we find that a weak covalent chemisorption of molecular hydrogen, known as Kubas interaction, is feasible on Ca decorated graphene according to DMC, in agreement with DFT. This finding is in contrast to previous DMC predictions of the 4H2/Ca+ gas cluster (without graphene) where chemisorption is not favored. However, we find that the adsorption energy of hydrogen on metal decorated graphene according to a widely used DFT method is not fully consistent with DMC. The reference adsorption energies reported herein can be used to find better work-horse methods for application in large-scale modeling of hydrogen adsorption. Furthermore, the implications of this work affect strategies for finding suitable hydrogen storage materials and high-throughput methods.

Chemistry↗

Intrinsic Limits of Charge Carrier Mobilities in Layered Halide Perovskites

Layered halide perovskites have emerged as potential alternatives to three-dimensional (3D) halide perovskites due to their improved stability and larger material phase space, allowing fine tuning of structural, electronic, and optical properties. However, their charge carrier mobilities are significantly smaller than those of 3D halide perovskites, which has a considerable impact on their application in optoelectronic devices. Here, we employ state-of-the-art approaches to unveil the electron-phonon mechanisms responsible for the diminished transport properties of layered halide perovskites. Starting from a prototypical A M X 3 halide perovskite, we model the case of n = 1 and n = 2 layered structures and compare their electronic and transport properties to the 3D reference. The electronic and phononic properties are investigated within density functional theory (DFT) and density functional perturbation theory (DFPT), while transport properties are obtained via the Boltzmann transport equation. The vibrational modes contributing to charge carrier scattering are investigated and associated with polar-phonon scattering mechanisms arising from the long-range Fröhlich coupling and deformation-potential scattering processes. Our investigation reveals that the lower mobilities in layered systems primarily originate from the increased electronic density of states at the vicinity of the band edges, while the electron-phonon coupling strength remains similar. Such an increase is caused by the dimensionality reduction and the break in octahedra connectivity along the stacking direction. Our findings provide a fundamental understanding of the electron-phonon coupling mechanisms in layered perovskites and highlight the intrinsic limitations of the charge carrier transport in these materials. Published by the American Physical Society 2024

Cucco, Bruno (ORCID:0000000331564143)↗

Computational study on the impact of gasoline-ethanol blending on autoignition and soot/NO x emissions under low-load gasoline compression ignition conditions

Here, in the present work, computational fluid dynamics (CFD) simulations of a single-cylinder gasoline compression ignition (GCI) engine are performed to investigate the impact of gasoline-ethanol blending on autoignition, nitrogen oxide (NO x ), and soot emissions under low-load conditions. In order to represent the test gasoline (RD5-87), a four-component toluene primary reference fuel (TPRF)+ethanol (ETPRF) surrogate (with 10% ethanol by volume; E10) is employed. A three-dimensional (3D) engine CFD model employing finite-rate chemistry with a skeletal kinetic mechanism (including NO x sub-mechanism), adaptive mesh refinement (AMR), and hybrid method of moments (HMOM) is adopted to capture the in-cylinder combustion phenomena and soot/NO x emissions. The engine CFD model is validated against experimental data for three gasoline-ethanol blends: E10, E30 and E100, with varying ethanol content by volume. Model validation is carried out for a broad range of start-of-injection (SOI) timings (−21, −27, −36, and −45 crank angle degrees (°CA) after top-dead-center (aTDC)) with respect to in-cylinder pressure, heat release rate, combustion phasing, NO x and soot emissions. For relatively later injection timings (−21 and −27 °CA aTDC), E30 yields higher amount of soot than E10; while the trend reverses for early injection cases (−36 and −45 °CA aTDC ). On the other hand, E100 yields the lowest amount of soot among all fuels irrespective of SOI timing. Further, E10 shows a non-monotonic trend in soot emissions with SOI timing: SOI-36>SOI-45>SOI-21>SOI-27, while soot emissions from E30 exhibit monotonic decrease with advancing SOI timing. NO x emissions from various fuels follow a trend of E10>E30>E100. On the other hand, NO x emissions increase as SOI timing is advanced for all fuels, with an anomaly for E10 and E100 where NO x decreases when SOI is advanced beyond −36 °CA aTDC. Detailed analysis of the numerical results is performed to investigate the soot/NO x emission trends and elucidate the impact of chemical composition and physical properties on autoignition and emissions characteristics.

Computational fluid dynamics↗