Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “INTERPOLATION”

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,423 records · Page 79
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Computational Screening of Cationic Zeolites for n -Butane/Methane Separations Using Quantitatively Accurate First-Principles-Derived Force Fields

We developed force fields for linear alkanes in Na and Ca-exchanged zeolites based on periodic first-principles calculations with coupled cluster corrections. These force fields were validated by comparing the simulated adsorption properties with extensive experimental data. As an example of using these force fields, we screened silica and aluminosilicate zeolites for separation of n-butane/methane mixtures in several different separation processes. A key step in this screening was the use of an interpolation method to create accurate zeolite structures that take account of the lattice constant changes for different Si/Al ratios for a given zeolite topology. We focused on zeolite structures that have known synthesis routes, so the materials selected from our calculations can be tested experimentally. Several promising materials that show good separation performance were selected for more detailed simulations.

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

Smooth Dispersion Is Physically Appropriate: Assessing and Amending the D4 Dispersion Model

The addition of dispersion corrections to density functionals is essential for accurate energy and geometry predictions. Among them, the D4 scheme is popular due to its low computational cost and high accuracy. However, due to its design, the D4 correction can occasionally lead to anomalies, such as unphysical curvature and bumps in the potential energy surface. We find these anomalies are common in the D4 model, although observable consequences are rarer than in the D3 model for reasons we explain. Nevertheless, we uncover instances of unphysical local minima and stationary points with the D4 scheme and propose two solutions that yield smoother dispersion energy as a function of nuclear position. One is trivial to implement, based on a smoother reparametrization of Gaussian weighting (D4S) to find the effective coordination number. The other replaces Gaussian weighting with soft linear interpolation (D4SL). These new approaches usually remove artificial extremum points, while maintaining accuracy.

Tkachenko, Nikolay V↗
DOE OSTI2026

Minimization of Disorder as a Key Design Principle for Natural Sizes of Light Harvesting 2 Complexes

The light harvesting 2 (LH2) complex of purple bacteria has excellent energy conversion efficiency. Clarifying the design principle behind such efficiency at the atomistic level is crucial for understanding its structure–function relationship and can be utilized for the design of artificial light harvesting systems. To this end, we conducted comprehensive computational investigation of the dynamical and statistical nature of electronic excited states of pigment molecules in a natural LH2 complex with 9-fold symmetry and its two non-natural in silico analogues with 6- and 12-fold symmetries. To ensure reliable and efficient all-atomistic molecular dynamics simulations, we combined a well established interpolation approach for the construction of the potential energy surface with a neural network machine learning approach. Outcomes of these calculations clarify that non-natural forms of LH2-type complexes have significantly larger quasistatic disorder than those for the natural one. In addition, non-natural systems have more disruptions of the hydrogen bonding, underscoring its crucial role for reducing the disorder. On the other hand, local environmental dynamics are relatively insensitive to the structural changes although there is moderate enhancement in the anharmonic or interatomic components for the synthetic ones. These findings based on all-atomistic simulations provide direct computational evidence that the structure and sizes of natural LH2 complexes are designed to minimize the energetic disorder. We analyze quantitative implications of these for the energy transferring capability of the LH2 complex.

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

A Preprocessing Tool for Enhanced Ion Mobility–Mass Spectrometry-Based Omics Workflows

The ability to improve the data quality of ion mobility–mass spectrometry (IM-MS) measurements is of great importance for enabling modular and efficient computational workflows and gaining better qualitative and quantitative insights from complex biological and environmental samples. We developed the PNNL PreProcessor, a standalone and user-friendly software housing various algorithmic implementations to generate new MS-files with enhanced signal quality and in the same instrument format. Different experimental approaches are supported for IM-MS based on Drift-Tube (DT) and Structures for Lossless Ion Manipulations (SLIM), including liquid chromatography (LC) and infusion analyses. The algorithms extend the dynamic range of the detection system, while reducing file sizes for faster and memory-efficient downstream processing. Specifically, multidimensional smoothing improves peak shapes of poorly defined low-abundance signals, and saturation repair reconstructs the intensity profile of high-abundance peaks from various analyte types. Further, other functionalities are data compression and interpolation, IM demultiplexing, noise filtering by low intensity threshold and spike removal, and exporting of acquisition metadata. Several advantages of the tool are illustrated, including an increase of 19.4% in lipid annotations and a two-times faster processing of LC-DT IM-MS data-independent acquisition spectra from a complex lipid extract of a standard human plasma sample. The software is freely available at https://omics.pnl.gov/software/pnnl-preprocessor.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

A First-Principles Approach to Modeling Surface Site Stabilities on Multimetallic Catalysts

The study of multimetallic alloys and the multitude of possible surface compositions have sparked a tremendous interest in engineering low-cost materials with high activity and selectivity in heterogeneous catalysis. Multimetallic systems provide complementary functionalities and an unprecedented tunability when designing catalyst formulations. However, due to their immense structural and compositional complexity, the investigation and identification of an optimal catalyst is a tedious and time-consuming process, both experimentally and theoretically. Therefore, theoretical design principles are highly desirable to accelerate the screening of catalyst structures across the vast compositional space. In this paper, we introduce a simple and general model for predicting the site stability of multimetallic surfaces and nanoparticles, which is based on physical principles. The model requires only a small set of density functional theory (DFT) calculations of metal atom binding energies on monometallic and dilute alloy surface slabs to optimize the parameters in the simple model. The resulting model allows for the quantification of the stability of any particular atom site in any conceivable chemical environment across a wide range of morphologies, sizes, and arrangements by interpolating the derived parameters from a monometallic system to a completely diluted alloyed system. Herein, we demonstrate the robustness of the model across an extensive data set of transition metal alloy surfaces and 147-atoms cuboctahedral nanoparticles (NPs) composed of IrRhRu and PtPdRu. In conclusion, our approach yields mean absolute errors of ≈0.15 (IrRhRu), 0.20 (PtPdRu), 0.19 (IrRhRu NP), and 0.26 (PtPdRu NP) eV relative to site binding energies calculated using DFT.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Synergistic Interactions During Co-Hydrothermal Liquefaction of Food Waste and Biomass Model Compounds for Increased Sustainable Aviation Fuel Production

Hydrothermal liquefaction (HTL) of lignocellulosic biomass is plagued with low biocrude yields owing to the tendency of highly reactive oxygenated intermediates to condense to form biochars. By contrast, the high protein content in food waste is comprised of substantial nitrogen species, which are known to interact strongly with oxygenates through Maillard, amide, and peptide bond formation reactions. Co-feeding food waste and lignocellulose opens new reaction pathways for biocrude formation but is currently poorly understood. This work evaluated the molecular level interactions between food waste and lignocellulose model compounds and the corresponding effect on product yields and quality. Food waste–cellulose and food waste–xylan feedstock blends achieved maximum biocrude carbon yield improvements of 12.2% and 10.1%, respectively, relative to a simple linear model that interpolates between the yields of the pure feedstocks. Increases in biocrude yield were balanced by corresponding decreases in char yield, indicating synergistic interactions between the feeds during HTL. Biocrude volatility analysis revealed that increased biocrude yield preferentially benefitted the jet fuel fraction, which comprised up to 22.6% of the total carbon yield for food waste–cellulose blends. Biocrude and char were analyzed using GC–MS and FT-IR spectroscopy to investigate the source of synergistic trends and provide greater mechanistic understanding. Key cofeeding effects included the promotion of retro-aldol condensation reactions and trans-esterification of fatty acids, sequestering carbon in the biocrude phase via the inhibition of char formation while increasing biocrude volatility toward jet fuel-range compounds. These results indicate the potential for judicious selection of HTL cofeeds to increase both biocrude yield and selectivity to desired fuel precursors, including sustainable aviation fuel.

Maillard reaction↗
DOE OSTI2026

Hourly PM 2.5 Estimates across California from 2018 to 2023

This study presents a new data set of hourly PM 2.5 concentrations across California from 2018 to 2023 at a three-kilometer resolution. This data set was developed by assimilating observations from PurpleAir and the U.S. EPA Air Quality System monitors into wildfire smoke forecasts from the High-Resolution Rapid Refresh Smoke (HRRR-Smoke) model using the Gridpoint Statistical Interpolation (GSI) three-dimensional variational data assimilation framework. Archived forecasts of modeled wildfire smoke PM 2.5 from HRRR-Smoke create the background field for assimilation, which is then corrected using surface observations of total PM 2.5 . The resulting reanalysis from GSI provides an estimate of total PM 2.5 that minimizes error from both the observational and the model data. Validation results indicate strong performance, with monthly R 2 values ranging from 0.73 to 0.91 across the six-year data set, comparable to other PM 2.5 data sets. Case studies are presented for three major fire events, the 2018 Camp Fire, 2019 Kincade Fire, and 2020 Lightning Complex Fires to demonstrate the data set’s fidelity in resolving plume dynamics and local exposure patterns. Root-mean-squared error averaged over each month scales with average PM 2.5 concentrations, resulting in a low error under typical conditions but higher absolute errors during extreme smoke events. This is the first long-term, hourly PM 2.5 data set of its kind for California and enables the generation of subdaily exposure metrics, such as peak hourly concentrations, exceedance durations, and time-of-day exposure peaks. The novelty and strong validation of this data set make it a compelling resource for future studies on the impact and significance of subdaily PM 2.5 exposure.

PM2.5↗
DOE OSTI2017

Implications of grain size variation in magnetic field alignment of block copolymer blends

Recent experiments have highlighted the intrinsic magnetic anisotropy in coil–coil diblock copolymers, specifically in poly(styrene- block -4-vinylpyridine) (PS- b -P4VP), that enables magnetic field alignment at field strengths of a few tesla. We consider here the alignment response of two low molecular weight (MW) lamallae-forming PS- b -P4VP systems. Cooling across the disorder–order transition temperature (T odt ) results in strong alignment for the higher MW sample (5.5K), whereas little alignment is discernible for the lower MW system (3.6K). This disparity under otherwise identical conditions of field strength and cooling rate suggests that different average grain sizes are produced during slow cooling of these materials, with larger grains formed in the higher MW material. Blending the block copolymers results in homogeneous samples which display T odt , d-spacings, and grain sizes that are intermediate between the two neat diblocks. Similarly, the alignment quality displays a smooth variation with the concentration of the higher MW diblock in the blends, and the size of grains likewise interpolates between limits set by the neat diblocks, with a factor of 3.5× difference in the grain size observed in high vs low MW neat diblocks. Finally, these results highlight the importance of grain growth kinetics in dictating the field response in block copolymers and suggests an unconventional route for the manipulation of such kinetics.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2014

Hybrid Quantum Mechanics/Molecular Mechanics Solvation Scheme for Computing Free Energies of Reactions at Metal–Water Interfaces

We report the development of a quantum mechanics/molecular mechanics free energy perturbation (QM/MM-FEP) method for modeling chemical reactions at metal-water interfaces. This novel solvation scheme combines planewave density function theory (DFT), periodic electrostatic embedded cluster method (PEECM) calculations using Gaussian-type orbitals, and classical molecular dynamics (MD) simulations to obtain a free energy description of a complex metal-water system. We derive a potential of mean force (PMF) of the reaction system within the QM/MM framework. A fixed-size, finite ensemble of MM conformations is used to permit precise evaluation of the PMF of QM coordinates and its gradient defined within this ensemble. Local conformations of adsorbed reaction moieties are optimized using sequential MD-sampling and QM-optimization steps. An approximate reaction coordinate is constructed using a number of interpolated states and the free energy difference between adjacent states is calculated using the QM/MM-FEP method. By avoiding on-the-fly QM calculations and by circumventing the challenges associated with statistical averaging during MD sampling, a computational speedup of multiple orders of magnitude is realized. The method is systematically validated against the results of ab initio QM calculations and demonstrated for C–C cleavage in double-dehydrogenated ethylene glycol on a Pt (111) model surface.

Chemical calculations↗
DOE OSTI2015

Charge Hopping Dynamics along a Disordered Chain in Quantum Environments: Comparative Study of Different Rate Kernels

This work presents a computational study of charge hopping dynamics along a one-dimensional chain with Gaussian site energy disorder and linearly coupled quantum bath. Time dependent square displacements are calculated directly from numerical solutions of Pauli master equations, for five different rate kernels: exact Fermi golden rule (FGR) rate expression, stationary phase interpolation (SPI) approximation, semiclassical (SC) approximation, classical Marcus rate expression, and Miller-Abrahams expression. All results demonstrate diffusive behavior in the steady state limit. The results based on the FGR rate expression show that the charge transport in the quantum bath can be much more sensitive to the disorder than the prediction from the classical Marcus expression. While the SPI approximation captures this general trend reasonably well, the SC approximation tends to be unreliable at both quantitative and qualitative levels and becomes even worse than the classical Marcus expression under certain conditions. Furthermore, these results offer useful guidance in the choice of approximate rate kernels for larger-scale simulations and also demonstrate significant but fragile positive effects of quantum environments on the charge hopping dynamics.

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

Assessing the Skill of the Improved Treatment of Riverine Freshwater in the Community Earth System Model (CESM) Relative to a New Salinity Climatology

Abstract Recent studies have explored the sensitivity of global ocean model simulations to the treatment of riverine freshwater and the representation of estuarine processes via an estuary box model applied within Community Earth System Model (CESM). This study builds on these efforts by assessing the model skill score relative to a new salinity climatology. The new climatology averages the original observational data of the World Ocean Database directly onto the CESM ocean component tracer grid cells without spatial interpolation, smoothing, or other gap‐filling techniques to mitigate coastal ocean salinity bias present in the World Ocean Atlas. The mean square error for coastal upper ocean salinity relative to climatology is reduced by up to 14%, and the mean square error of near‐surface salinity stratification is reduced by up to 28% near major river mouths in the simulations with improved treatments of river runoff. The improvement in upper ocean bulk salinity is attributed primarily to focusing runoff as point sources thereby avoiding the artificial horizontal spreading of the control run and to applying a locally varying instead of a global constant reference salinity for riverine virtual salt fluxes. The improvements in near‐surface salinity stratification are primarily attributed to adding parameterized estuarine mixing with the estuary box model. Salinity and salinity stratification skill improvements are achieved not just near large rivers but also along the global coast and skill improvements extend far offshore. Despite these improvements, many other sources of model‐climatology mismatch in coastal salinity and stratification remain and merit further attention.

58 GEOSCIENCES↗
DOE OSTI2020

Simulating and Evaluating Atmospheric River-Induced Precipitation Extremes Along the U.S. Pacific Coast: Case Studies From 1980–2017

Atmospheric rivers (ARs) are responsible for a majority of extreme precipitation and floodevents along the U.S. West Coast. To better understand the present-day characteristics of AR-relatedprecipitation extremes, a selection of nine most intense historical AR events during 1980–2017 is simulatedusing a dynamical downscaling modeling framework based on the Weather Research and ForecastingModel. We find that the chosen framework and Weather Research and Forecasting Model configurationreproduces both large-scale atmospheric features—including parent synoptic-scale cyclones—as well asthe filamentary corridors of integrated vapor transport associated with the ARs themselves. The accuracyof simulated extreme precipitation maxima, relative to in situ and interpolated gridded observations,improves notably with increasing model resolution, with improvements as large as 40–60% for fine scale(3 km) relative to coarse-scale (27 km) simulations. A separate set of simulations using smoothedtopography suggests that much of these gains stem from the improved representation of complex terrain.Additionally, using the 12 December 1995 storm in Northern California as an example, we demonstratethat only the highest-resolution simulations resolve important fine-scale features—such as localizedorographically forced vertical motion and powerful near hurricane-force boundary layer winds. Given thedemonstrated ability of a targeted dynamical downscaling framework to capture both local extremeprecipitation and key fine-scale characteristics of the most intense ARs in the historical record, we arguethat such a configuration may be highly conducive to understanding AR-related extremes and associatedchanges in a warming climate

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

An Assessment of Global Positioning System Velocity Uncertainty in California

We analyze data from 580 continuous global positioning system (GPS) stations in California to quantify differences in published velocity estimates from five analysis centers. Horizontal and vertical rates for individual stations can differ up to 5 mm/yr, with systematic differences in some areas comparable to deformation rates. Published velocity uncertainties vary between analysis centers and are systematically underreported in the horizontal relative to empirical uncertainties calculated from the scatter of analysis center velocities. In the vertical, published velocity uncertainties are both over and underreported and vary more widely between centers. An interpolated ensemble vertical velocity field shows high-subsidence regions in the Central Valley and Salton Trough have the largest empirical uncertainties, while station density has a modest impact on uncertainties. Applications that rely on subcentimeter GPS accuracy should consider the possibility that formal errors published with velocity rate estimates understate true velocity uncertainties in both the horizontal and vertical.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

A 3D Full Stress Tensor Model for Oklahoma

Abstract The stress tensor is an important property for upper crustal studies such as those that involve pore fluids and earthquake hazards. At tectonic plate scale, plate boundary forces and mantle convection are the primary drivers of the stress field. In many local settings (10–100 s of km and <10 km depth) in tectonic plate interiors, we can simplify by assuming a constant background stress field that is perturbed by local heterogeneity in density and elasticity. Local stress orientation and sometimes magnitude can be estimated from earthquake and borehole‐based observations when available. Modeling of the local stress field often involves interpolating sparse observations. We present a new method to estimate the 3D stress field in the upper crust and demonstrate it for Oklahoma. We created a 3D material model by inverting multiple types of geophysical observations simultaneously. Integrating surface‐wave dispersion, local travel times and gravity observations produces a model of P‐wave velocity, S‐wave velocity, and density. The stress field can then be modeled using finite element simulations. The simulations are performed using our simplified view of the local stress field as the sum of a constant background stress field that is perturbed by local density and elasticity heterogeneity and gravitational body forces. An orientation of N82°E, for the maximum compressive tectonic force, best agrees with previously observed stress orientations and faulting types in Oklahoma. The gravitational contribution of the horizontal stress field has a magnitude comparable to the tectonic contribution for the upper 5 km of the subsurface.

58 GEOSCIENCES↗
DOE OSTI2020

Simulation of Continental Shallow Cumulus Populations Using an Observation-Constrained Cloud-System Resolving Model

Continental shallow cumulus (ShCu) clouds observed on 30 August 2016 during the Holistic Interactions of Shallow Clouds, Aerosols, and Land-Ecosystems (HI-SCALE) field campaign are simulated by using an observation-constrained cloud-system resolving model. On this day, ShCu forms over Oklahoma and southern Kansas and some of these clouds transition to deeper, precipitating convection during the afternoon. We apply a four-dimensional ensemble-variational (4DEnVar) hybrid technique in the Community Gridpoint Statistical Interpolation (GSI) system to assimilate operational data sets and unique boundary layer measurements including a Raman lidar, radar wind profilers, radiosondes, and surface stations collected by the U.S. Department of Energy's (DOE) Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) atmospheric observatory into the Weather Research and Forecasting (WRF) model to ascertain how improved environmental conditions can influence forecasts of ShCu populations and the transition to deeper convection. Independent observations from aircraft, satellite, as well as ARM's remote sensors are used to evaluate model performance in different aspects. Several model experiments are conducted to identify the impact of data assimilation (DA) on the prediction of clouds evolution. The analyses indicate that ShCu populations are more accurately reproduced after DA in terms of cloud initiation time and cloud base height, which can be attributed to an improved representation of the ambient meteorological conditions and the convective boundary layer. Extending the assimilation to 18 UTC (local noon) also improved the simulation of shallow-to-deep transitions of convective clouds.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Quantitative Precipitation Estimation of Extremes in CONUS With Radar Data

When constructing an accurate, continental, in-situ-based, kilometer-scale, long-term record of the precipitation field and its spatiotemporal changes remains a significant challenge. Here, we determine the extreme-value behavior of the NEXRAD Stage IV radar-based quantitative precipitation estimate. We find that the climatology of 5-year daily return values in the contiguous United States East of the Rocky Mountains shows only slight variability on spatial scales smaller than (Formula presented.) ~100 km. In light of this finding, we test whether rain-gauge-only daily precipitation data sets can produce accurate extreme-value behavior at spatial scales finer than the spacing between gauges. We find that the 5-year daily return values are accurate at locations far from rain gauges only if the interpolation between gauges is carried out appropriately for extremes. Precipitation statistics derived from in-situ rain gauge data are therefore of sufficient spatial resolution to faithfully capture daily extremes over much of the eastern United States.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Distributed Fiber Optics Measurements of Rock Deformation and Failure in Triaxial Tests

The implementation of distributed strain measurement methods in triaxial and uniaxial tests have demonstrated the development of strain localization, even at early stages. This implies that single point measurement methods are location dependent. The use of distributed methods is required not only to improve the interpreted constitutive parameters obtained from triaxial tests, but also to understand the implications of strain localization in the failure process. Here, this work uses optical fibers in triaxial tests. The developed distributed measurement method was implemented on granite, gabbro, and sandstone samples and tested under different confining pressures, reaching 200 MPa on the granite sample. Using a temporal resolution of 0.25 Hz and 5 mm of spatial resolution, the strain evolution at over 300 locations at the sample surface was measured during testing. When compared to point sensing methods, the use of optical fiber greatly increases the number of measurements at the surface of the sample, and their interpolation provides the entire deformation of the sample surface. In all the tests performed, strain localizations were revealed before failure. A three-dimensional interpretation of a test, combining an optical scan of the sample and the distributed measurements, show good correlation between the fractures and the strain localization.

58 GEOSCIENCES↗
DOE OSTI2022

Equation‐Free Surrogate Modeling of Geophysical Flows at the Intersection of Machine Learning and Data Assimilation

Abstract There is a growing interest in developing data‐driven reduced‐order models for atmospheric and oceanic flows that are trained on data obtained either from high‐resolution simulations or satellite observations. The data‐driven models are non‐intrusive in nature and offer significant computational savings compared to large‐scale numerical models. These low‐dimensional models can be utilized to reduce the computational burden of generating forecasts and estimating model uncertainty without losing the key information needed for data assimilation (DA) to produce accurate state estimates. This paper aims at exploring an equation‐free surrogate modeling approach at the intersection of machine learning and DA in Earth system modeling. With this objective, we introduce an end‐to‐end non‐intrusive reduced‐order modeling (NIROM) framework equipped with contributions in modal decomposition, time series prediction, optimal sensor placement, and sequential DA. Specifically, we use proper orthogonal decomposition (POD) to identify the dominant structures of the flow, and a long short‐term memory network to model the dynamics of the POD modes. The NIROM is integrated within the deterministic ensemble Kalman filter (DEnKF) to incorporate sparse and noisy observations at optimal sensor locations obtained through QR pivoting. The feasibility and the benefit of the proposed framework are demonstrated for the NOAA Optimum Interpolation Sea Surface Temperature (SST) V2 data set. Our results indicate that the NIROM is stable for long‐term forecasting and can model dynamics of SST with a reasonable level of accuracy. Furthermore, the prediction accuracy of the NIROM gets improved by almost one order of magnitude by the DEnKF algorithm.

Pawar, Suraj↗