SEARCH · Engineering Papers

Results for “empirical validation”

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.

SourceResearch connectionFrom yearThrough year

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

At least 433 records · Page 24

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Revisiting Wagner’s Criteria to Predict Establishment and Retention of Alumina Scales on Ternary NiCrAl Alloys

Despite the advent of numerical modeling approaches and high-performance computing infrastructure, the design and development of corrosion-resistant high temperature alloys (> 500 °C) continue to be largely empirical and typically involve extensive experimentation. This is mainly due to the lack of a single unified physics-based model that can address the impact of multiple competing factors such as time, environment, alloy composition, microstructure, and geometry. The classical Wagner’s criteria have been foundational to estimate the minimum concentrations required of an oxide-forming element to establish and sustain a protective oxide scale. However, the formulation is primarily limited to lower-order alloy systems (binary alloys) and ignores the time dependence of subsurface compositional changes in the alloy. The lack of key data on the temperature and composition dependence of the solubility and transport of oxidants in multicomponent-multiphase alloys further exacerbates the problem. In the present work, a few of these limitations were addressed using a flux-based approach (FLAP) which tracks the spatiotemporal evolution of the fundamental flux balance between oxygen and the oxide-forming elements to enable the prediction of the formation of an external alumina scale in ternary NiCrAl alloys. In conclusion, the modeling results were validated with the literature findings and additional experimental work conducted in the present work.

36 MATERIALS SCIENCE↗

DOE OSTI2022

A Universal Refrigerant Charge Fault Detection and Diagnostics Method Based on Pump Down Operation

The performance of the heat pump system varies greatly depending on the refrigerant charge amount. Improving the refrigerant charge fault detection and diagnostics (FDD) method of vapor compression systems have the potential for increasing energy efficiency and reducing service cost. Previous studies to predict refrigerant charge amount are mostly empirical methods which require significant amount of experimental data for high accuracy. The primary goal of this research is to develop a universal charge fault detection method which requires only a few experimental data with high prediction accuracy.Currently, pump down operations are typical practices by HVAC technicians when they need to open the refrigerant circuit to make a repairment. In addition, compressors have a built-in low-pressure cut-off protection function, and the compressor performance maps are commonly available from manufacturers. The proposed method innovatively utilizes the typical pump down operation, the compressor low-pressure cut-off protection, and the compressor performance map. It does not require any geometry information of heat exchangers, refrigerant lines, or charge buffers.The new charge prediction method is firstly formulated through theoretical analysis, then verified and calibrated by a quasi-steady-state simulation of the pump down process for a residential heat pump system. The quasi steady-state simulation uses an HVAC system simulation framework driven by DOE/ORNL Heat Pump Design Model (HPDM). Preliminary experiment validations with heat pump refrigerant leakage tests demonstrate the deviation of the proposed charge prediction method compared with measurement is within 8%. This technology makes refrigerant charge amount available at the technician’s fingertips and leads to shorter maintenance time and fewer site visits.

Li, Zhenning↗

DOE OSTI2020

Medium Energy Electron Flux in Earth's Outer Radiation Belt (MERLIN): A Machine Learning Model

The radiation belts of the Earth, filled with energetic electrons, comprise complex and dynamic systems that pose a significant threat to satellite operation. While various models of electron flux both for low and relativistic energies have been developed, the behavior of medium energy (120–600 keV) electrons, especially in the MEO region, remains poorly quantified. At these energies, electrons are driven by both convective and diffusive transport, and their prediction usually requires sophisticated 4D modeling codes. In this paper, we present an alternative approach using the Light Gradient Boosting (LightGBM) machine learning algorithm. The Medium Energy electRon fLux In Earth's outer radiatioN belt (MERLIN) model takes as input the satellite position, a combination of geomagnetic indices and solar wind parameters including the time history of velocity, and does not use persistence. MERLIN is trained on >15 years of the GPS electron flux data and tested on more than 1.5 years of measurements. Tenfold cross validation yields that the model predicts the MEO radiation environment well, both in terms of dynamics and amplitudes o f flux. Evaluation on the test set shows high correlation between the predicted and observed electron flux (0.8) and low values of absolute error. The MERLIN model can have wide space weather applications, providing information for the scientific community in the form of radiation belts reconstructions, as well as industry for satellite mission design, nowcast of the MEO environment, and surface charging analysis.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2023

An adaptive sampling augmented Lagrangian method for stochastic optimization with deterministic constraints

The primary goal of this paper is to provide an efficient solution algorithm based on the augmented Lagrangian framework for optimization problems with a stochastic objective function and deterministic constraints. Our main contribution is combining the augmented Lagrangian framework with adaptive sampling, resulting in an efficient optimization methodology validated with practical examples. To achieve the presented efficiency, here we consider inexact solutions for the augmented Lagrangian subproblems, and through an adaptive sampling mechanism, we control the variance in the gradient estimates. Furthermore, we analyze the theoretical performance of the proposed scheme by showing equivalence to a gradient descent algorithm on a Moreau envelope function, and we prove sublinear convergence for convex objectives and linear convergence for strongly convex objectives with affine equality constraints. The worst-case sample complexity of the resulting algorithm, for an arbitrary choice of penalty parameter in the augmented Lagrangian function, is $\mathscr{O}$(ϵ -3-δ ) , where ϵ > 0 is the expected error of the solution and δ > 0 is a user-defined parameter. If the penalty parameter is chosen to be $\mathscr{O}$(ϵ -1 ), we demonstrate that the result can be improved to $\mathscr{O}$(ϵ -2 ) , which is competitive with the other methods employed in the literature. Moreover, if the objective function is strongly convex with affine equality constraints, we obtain $\mathscr{O}$(ϵ -1 log(1/ϵ)) complexity. Finally, we empirically verify the performance of our adaptive sampling augmented Lagrangian framework in machine learning optimization and engineering design problems, including topology optimization of a heat sink with environmental uncertainty.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Machine learning prediction on the fractional free volume of polymer membranes

Fractional free volume (FFV) characterizes the microstructural level features of polymers and affects their properties including thermal, mechanical, and separation performance. Experimental measurements and theoretical analyses have been used to quantify the FFV of polymers, but challenges remain because of their limitations. Experimental measurements are laborious and based on semi empirical equations, while Bondi’s group contribution theory involves ambiguities like the determination of van der Waals volume and the choice of factor values in the theoretical equation. To efficiently evaluate the FFV of polymers, this study utilizes high-throughput molecular dynamics (MD) simulations to build a large dataset regarding polymer’s FFV. Based on this large dataset, we further build machine learning (ML) models to establish the composition-structure relation. Inspired by group contribution theory which correlates polymer’s functional groups to FFV, our ML models correlate polymer’s substructures or physico-chemical indexes to FFV. Here, our study first benchmarks the MD simulation protocol to obtain reliable FFV of polymers and then carries out high-throughput MD simulations for more than 6,500 homopolymers and 1,400 polyamides. Such a large and diverse dataset makes the well-trained ML models more generalizable, compared with the group contribution theory. The efficiency of a feed forward neural network model is further demonstrated by applying it to a hypothetical polyimide dataset of more than 8 million chemical structures. The predicted FFVs of hypothetical polyimides are further validated by MD simulations. The obtained FFVs of the 8 million polymers, plus their previously reported gas separation performances, demonstrate the promising capability of ML virtual screening for the discovery of polymer membranes with exceptional permeability/selectivity.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Correlation-aware binning for small-angle neutron scattering via Gaussian-process inference

Binning in small-angle neutron scattering (SANS) is typically performed empirically, with fixed parameters chosen for convenience rather than statistical optimality. Such practices often fail to balance statistical precision and spatial resolution, leading to inconsistencies across instruments and datasets. Here we establish a correlation-aware framework that determines the optimal bin width from first principles by extending the classical Freedman–Diaconis (FD) rule to account for inter-bin correlations with a Gaussian process. In this formulation, the scattering intensity is treated as a smooth stochastic field whose statistical coherence is described by a covariance matrix. Analytical expressions of errors derived from this model yield closed-form criteria that separate the total deviation into contributions from counting noise, aliasing distortion and curvature-dependent correlation effects. Expressed in reduced variables, the resulting dimensionless error surface reveals a continuous transition from the uncorrelated FD regime to the correlation-dominated limit, providing a unified description of noise suppression and resolution control. Because the formulation depends only on the profile characteristics of scattering intensity I(Q), specifically its average intensity and first- and second-order derivatives, it applies generally to any SANS measurement regardless of sample, instrument or geometry. Experimental validation using small- and ultra-small-angle neutron scattering data confirms the predicted scaling behavior, demonstrating that correlation-aware inference systematically reduces mean-squared error and enables information-efficient reproducible data reduction across materials and instruments.

Tung, Chi-Huan [ORNL] (ORCID:0000000221972074)↗

DOE OSTI2020

Theoretical Understanding of the Linear Relationship between Convective Updrafts and Cloud-Base Height for Shallow Cumulus Clouds. Part II: Continental Conditions

This is the Part II of a two-part study that seeks a theoretical understanding of an empirical relationship for shallow cumulus clouds: subcloud updraft velocity covaries linearly with the cloud-base height. This work focuses on continental cumulus clouds that are more strongly forced by surface fluxes and more deviated from equilibrium than those over oceans (Part I). We use a simple analytical model for shallow cumulus that is well tested against a high-resolution (25 m in the horizontal) large-eddy simulation model. Consistent with a conventional idea, we find that surface Bowen ratio is the key variable that regulates the covariability of both parameters: under the same solar insolation, a drier surface allows for stronger buoyancy flux, triggering stronger convection that deepens the subcloud layer. We find that the slope of the Bowen-ratio-regulated relationship between the two parameters (defined as λ) is dependent on both the local time and the stability of the lower free atmosphere. The value of λ decreases with time exponentially from sunrise to early afternoon and linearly from early afternoon to sunset. The value of λ is larger in a more stable atmosphere. In addition, continental λ in the early afternoon more than doubles the oceanic λ. Validation of the theoretical results against ground observations over the Southern Great Plains shows a reasonable agreement. Physical mechanisms underlying the findings are explained from the perspective of different time scales at which updrafts and cloud-base height respond to a surface flux forcing.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2025

A statistical and simulation-informed model for estimating permeability from pore size distribution in saturated geomaterials

Accurate permeability estimation is essential across subsurface engineering applications but remains challenging due to the complex pore structures of natural geomaterials. Traditional empirical methods and simplified theoretical models often inadequately capture the role of pore size distribution and connectivity. Here, this study develops a statistical and simulation-informed permeability model that collapses pore-scale complexity into a compact scaling of the form k = αϕμ d 2 , where ϕ is porosity, μ d is mean pore size, and α is a weakly varying coefficient. By combining pore network simulations with statistical analysis of unimodal and bimodal pore size distributions, we identify three key findings: (i) permeability is much more sensitive to mean pore size than to porosity; (ii) across extensive datasets, the ratio σ d /μ d (standard deviation to mean) clusters around a characteristic value ∼0.4, allowing the effects of the full pore size distribution to be represented by μ d and a narrowly varying α ≈ 0.05; and (iii) for bimodal systems, there exists a critical fraction of small pores ∼0.78 above which flow becomes small-pore dominated, enabling the definition of an effective flow-controlling pore population and facilitating simplified permeability estimation for such systems. The resulting model, which requires only porosity and a representative mean pore size as inputs, is validated against comprehensive experimental datasets (>1700 samples) spanning diverse soils and rocks and achieves good predictive accuracy. Overall, this work provides a physically grounded yet practically simple permeability estimator suitable for subsurface engineering, environmental protection, and resource management applications.

Permeability↗

DOE OSTI2020

Reynolds-Averaged Turbulence Modeling Using Deep Learning with Local Flow Features: An Empirical Approach

Reynolds-Averaged Navier-Stoke (RANS) models offer an alternative avenue in predicting flow characteristics when the corresponding experiments are difficult to achieve due to geometry complexity, limited budget, or knowledge. RANS models require the knowledge of subgrid scale physics to solve conservation equations for mass, energy, and momentum. Mechanistic turbulence models, such as k-epsilon, are generally evaluated and calibrated for specific flow conditions with various degrees of uncertainty. These models have limited capability to assimilate a substantial amount of data due to model form constraints. Meanwhile, deep learning (DL) has been proven to be universal approximators with the potential to assimilate available, relevant, and adequately evaluated data. Moreover, deep neural networks (DNNs) can create surrogate models without knowing function forms. Such a data-driven approach can be used in updating fluid models based on observations as opposed to hard-wiring models with precalibrated correlations. The paper presents progress in applying DNNs to model Reynolds stress using two machine learning (ML) frameworks. A novel flow feature coverage mapping is proposed to quantify the physics coverage of DL-based closures. It can be used to examine the sufficiency of training data and input flow features for data-driven turbulence models. The case of a backward-facing step is formulated to demonstrate that not only can DNNs discover underlying correlation behind fluid data but also they can be implemented in RANS to predict flow characteristics without numerical stability issues. Finally, the presented research is a crucial stepping-stone toward the data-driven turbulence modeling, which potentially benefits the design of data-driven experiments that can be used to validate fluid models with ML-based fluid closures.

42 ENGINEERING↗

DOE OSTI2020

Evolving Trust in Robots: Specification Through Sequential and Comparative Meta-Analyses

Objective The objectives of this meta-analysis are to explore the presently available empirical findings on the antecedents of trust in robots and use this information to expand upon a previous meta-analytic review of the area. Background Human–robot interaction (HRI) represents an increasingly important dimension of our everyday existence. Currently, the most important element of these interactions is proposed to be whether the human trusts the robot or not. We have identified three overarching categories that exert effects on the expression of trust. These consist of factors associated with (a) the human, (b) the robot, and (c) the context in which any specific HRI event occurs. Method The current body of literature was examined and all qualifying articles pertaining to trust in robots were included in the meta-analysis. A previous meta-analysis on HRI trust was used as the basis for this extended, updated, and evolving analysis. Results Multiple additional factors, which have now been demonstrated to significantly influence trust, were identified. The present results, expressed as points of difference and points of commonality between the current and previous analyses, are identified, explained, and cast in the setting of the emerging wave of HRI. Conclusion The present meta-analysis expands upon previous work and validates the overarching categories of trust antecedent (human-related, robot-related, and contextual), as well as identifying the significant individual precursors to trust within each category. A new and updated model of these complex interactions is offered. Application The identified trust factors can be used in order to promote appropriate levels of trust in robots.

Behavioral Sciences↗

DOE OSTI2023

Enhancing the accuracy of density functional tight binding models through ChIMES many-body interaction potentials

Semi-empirical quantum models such as Density Functional Tight Binding (DFTB) are attractive methods for obtaining quantum simulation data at longer time and length scales than possible with standard approaches. However, application of these models can require lengthy effort due to the lack of a systematic approach for their development. In this work, we discuss the use of the Chebyshev Interaction Model for Efficient Simulation (ChIMES) to create rapidly parameterized DFTB models, which exhibit strong transferability due to the inclusion of many-body interactions that might otherwise be inaccurate. We apply our modeling approach to silicon polymorphs and review previous work on titanium hydride. We also review the creation of a general purpose DFTB/ChIMES model for organic molecules and compounds that approaches hybrid functional and coupled cluster accuracy with two orders of magnitude fewer parameters than similar neural network approaches. In all cases, DFTB/ChIMES yields similar accuracy to the underlying quantum method with orders of magnitude improvement in computational cost. In conclusion, our developments provide a way to create computationally efficient and highly accurate simulations over varying extreme thermodynamic conditions, where physical and chemical properties can be difficult to interrogate directly, and there is historically a significant reliance on theoretical approaches for interpretation and validation of experimental results.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Determination of confinement regime boundaries via separatrix parameters on Alcator C-Mod based on a model for interchange-drift-Alfvén turbulence

The separatrix operational space (SepOS) model (Eich and Manz 2021 Nucl. Fusion 61 086017) is shown to predict the L–H transition, the L-mode density limit, and the ideal magnetohydrodynamic ballooning limit in terms of separatrix parameters for a wide range of Alcator C-Mod plasmas. The model is tested using Thomson scattering measurements across a wide range of operating conditions on C-Mod, spanning $\overline{n}_{\mathrm{e}}$ = 0.3–5.5 $\times 10^{20}\,$m −3 , $B_{\mathrm{t}} = 2.5$–8.0 T, and $B_{\mathrm{p}}$= 0.1–1.2 T. An empirical regression for the electron pressure gradient scale length, $\lambda_{{p}_{\mathrm{e}}}$, against a turbulence control parameter, $\alpha_{\mathrm{t}}$, and the poloidal fluid gyroradius, $\rho_{\mathrm{s,p}}$, is constructed for H-modes and found to require positive exponents for both regression parameters, indicating turbulence widening of near-scrape-off layer widths at high $\alpha_{\mathrm{t}}$ and an inverse scaling with $B_{\mathrm{p}}$, consistent with results on ASDEX Upgrade. The SepOS model is also tested in the unfavorable drift direction and found to apply well to all three boundaries, including the L–H transition as long as a correction to the Reynolds energy transfer term, $\alpha_\mathrm{RS} \lt 1$ is applied. I-modes typically exist in the unfavorable drift direction for values of $\alpha_{\mathrm{t}} \lesssim 0.3$. Finally, an experiment studying the transition between the Type-I ELMy and EDA H-mode is analyzed using the same framework. It is found that a recently identified boundary $\alpha_{\mathrm{t}} = 0.55$ at excludes most EDA H-modes but that the balance of wavenumbers responsible for the L-mode density limit, namely $k_\mathrm{EM} = k_\mathrm{RBM}$, may better describe the transition on C-Mod. The ensemble of boundaries validated and explored is then applied to project regime access and limit avoidance for the SPARC primary reference discharge parameters.

ELM suppression↗

DOE OSTI2026

A Nonergodic Ground-Motion Model for the San Francisco Bay Area for Small-Magnitude Earthquakes

ABSTRACT Recently, generative models have become a computationally efficient alternative to physics-based numerical simulations of ground motions. Neural networks can learn from existing ground-motion data to generate unobserved ground-motion data at new source and site locations. A key challenge with generative models is ensuring that predicted ground motions remain within a physically realistic range. For this purpose, we developed an empirical, nonergodic ground-motion model (GMM) for small-magnitude earthquakes in the San Francisco Bay area based on about 5000 recordings per component for Mw ≤ 4 earthquakes. The nonergodic GMM predicts spatially varying median source, site, and path effects for both the Fourier amplitude spectrum (FAS) and the Fourier phase derivative (a proxy for duration), as well as the corresponding epistemic uncertainty for each term. For FAS, our model shows above-average source and site effects in the western part of the region and below-average effects in the eastern part, with regional effects exhibiting larger spatial correlation lengths with increasing frequency. For duration, the source term is negligible for small-magnitude earthquakes, and the site term leads to site-specific variations up to 5 s. Path effects for FAS and duration depend on the source–site pair and are extrapolated spatially using recent methods for path-effect modeling. The aleatory variability of the within-site within-path residuals is similar to the variability found in previous studies for other regions. The nonergodic model provides two key contributions: first, median adjustment terms that are transferable to larger magnitude earthquakes, further reducing aleatory variability in probabilistic seismic hazard analysis; second, region-specific criteria for validating machine learning-based ground-motion generators to evaluate whether synthetic ground motions exhibit physically realistic source, site, and path effects.

Lacour, Maxime↗

DOE OSTI2026

Development of MOSCATO: A CFD-Level Electrochemistry and Corrosion Simulator for Molten Salt Systems

For both coolant and fueled variants of molten salt reactors (MSRs), the corrosion of structural materials is a significant challenge. The corrosion stems from chemical and electrochemical reactions initiated by fissile material, fission products, and impurities in the salt. Lower-fidelity models rely on empirical correlations for mass transfer, simplified lumped temperature profiles, and similar assumptions. They do not capture detailed spatial variations in complex geometries, creating the need for high-fidelity modeling to bridge this gap.As we approach the demonstration and possible deployment of MSRs in this decade, the development of a high-fidelity, high-performance simulator becomes imperative. To simulate the complex electrochemical environment and corrosion within molten salt systems, we have developed the Molten Salt Chemistry And TranspOrt (MOSCATO) code. This endeavor is comprised of three essential components. First, mass transfer equations are coupled with the Navier-Stokes equations in order to account for the transport of species in the salt. Second, the diffusion of alloy constituents, such as Cr, Fe, Ni, etc. is simulated within the structural metals. Third, the alloy and salt domains are coupled to account for the heterogeneous chemical and electrochemical reactions that occur at the salt-alloy interface.MOSCATO manages all three components within the framework of the highly scalable, open-source spectral element method computational fluid dynamics code Nek5000/NekRS. This integration enables MOSCATO to harness the immense computational power of modern high-performance computing resources, ensuring both high fidelity and computational speed.In addition to code development, we have initiated a comprehensive verification and validation campaign, utilizing data from diverse sources. First, MOSCATO's electrochemical solver was verified with reference numerical data. Then validation occurred against experiments: one of a thermal galvanic cell and the other for corrosion in flowing molten salt of FLiNaK (LiF-NaF-KF). This campaign verified and validated MOSCATO as a reliable tool for simulating electrochemical environments and corrosion in molten salt systems.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

DOE OSTI2023

Building a DFT+U machine learning interatomic potential for uranium dioxide

Despite uranium dioxide (UO 2 ) being a widely used nuclear fuel, fuel performance models rely extensively on empirical correlations of material behavior, leveraging the historical operating experience of UO 2 . Mechanistic models that consider an atomistic understanding of the processes governing fuel performance (such as fission gas release and creep) will enable a better description of fuel behavior under non-prototypical conditions such as in new reactor concepts or for modified UO 2 fuel compositions. To this end, molecular dynamics simulation is a powerful tool for rapidly predicting physical properties of proposed fuel candidates. However, the reliability of these simulations depends largely on the accuracy of the atomic forces. Traditionally, these forces are computed using either a classical force field (FF) or density functional theory (DFT). While DFT is relatively accurate, the computational cost is burdensome, especially for f-electron elements, such as actinides. By contrast, classical FFs are computationally efficient but are less accurate. For these reasons, we report a new accurate machine learning interatomic potential (MLIP) for UO 2 that provides high-fidelity reproduction of DFT forces at a similar low cost to classical FFs. We employ an active learning approach that autonomously augments the DFT training data set to iteratively refine the MLIP. To further improve the quality of our predictions, we utilize transfer learning to retrain our MLIP to higher-accuracy DFT+U data. We validate our MLIPs by comparing predicted physical properties (e.g., thermal expansion and elastic properties) with those from existing classical FFs and DFT/DFT+U calculations, as well as with experimental data when available.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

DOE OSTI2023

Comparing simulated demand flexibility against actual performance in commercial office buildings

Commercial building energy benchmarking has been used as a mechanism to evaluate energy use of a single building over time, relative to other similar buildings, or to simulations of a reference building conforming to various energy standards. Lack of empirical demand flexibility data and consistent flexibility metrics has limited the ability to compare demand flexibility performance with estimated demand flexibility in buildings. In this study, we collected demand response performance data for a total of 831 demand response events from 192 sites as a first step to build such a building demand flexibility dataset, and propose a standard core data schema to consolidate field data from different sources. We also performed parametric simulations of a control strategy called “global temperature adjustment” using commercial office prototype building models. We then compared the simulated demand flexibility performance against the actual data for offices with global temperature adjustment strategy implemented. During demand response events with an average outside air temperature of 34 °C (range 23 °C–42 °C), the measured demand decrease intensity of the demand flexibility metrics were 6.1 watts per square meter (W/m 2 ), 10.0 W/m 2 , 11.1 W/m 2 , 7.1 W/m 2 , and 4.7 W/m 2 for small, small–medium, medium, medium–large, and large office buildings, respectively. Compared to the measured data in medium- and large-size buildings, the simulated demand decrease intensity was 0.7 W/m 2 (17%) lower on average. The discrepancy between simulated and measured peak demand intensities fell within one standard deviation of the mean measured data. Here, the comparison results validate the credibility of simulations in capturing real building data for assessing the technical potential of building demand flexibility.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

DOE OSTI2022

A comprehensive experimental and kinetic modeling study of di-isobutylene isomers: Part 2

A wide variety of high temperature experimental data obtained in this study complement the data on the oxidation of the two di-isobutylene isomers presented in Part I and offers a basis for an extensive validation of the kinetic model developed in this study. Due to the increasing importance of unimolecular decomposition reactions in high-temperature combustion, we have investigated the di-isobutylene isomers in high dilution utilizing a pyrolysis microflow reactor and detected radical intermediates and stable products using vacuum ultraviolet (VUV) synchrotron radiation and photoelectron photoion coincidence (PEPICO) spectroscopy. Additional speciation data at oxidative conditions were also recorded utilizing a plug flow reactor at atmospheric pressure in the temperature range 725-1150 K at equivalence ratios of 1.0 and 3.0 and at residence times of 0.35 s and 0.22 s, respectively. Combustion products were analyzed using gas chromatography (GC) and mass spectrometry (MS). Ignition delay time measurements for di-isobutylene were performed at pressures of 15 and 30 bar at equivalence ratios of 0.5, 1.0, and 2.0 diluted in 'air' in the temperature range 900-1400 K using a high-pressure shock-tube facility. New measurements of the laminar burning velocities of di-isobutylene/air flames are also presented. The experiments were performed using the heat flux method at atmospheric pressure and initial temperatures of 298-358 K. Moreover, data consistency was assessed with the help of analysis of the temperature and pressure dependencies of laminar burning velocity measurements, which was interpreted using an empirical power-law expression. Electronic structure calculations were performed to compute the energy barriers to the formation of many of the product species formed. The predictions of the present mechanism were found to be in adequate agreement with the wide variety of experimental measurements performed.

09 BIOMASS FUELS↗

DOE OSTI2021

A New Field Instrument for Leaf Volatiles Reveals an Unexpected Vertical Profile of Isoprenoid Emission Capacities in a Tropical Forest

Both plant physiology and atmospheric chemistry are substantially altered by the emission of volatile isoprenoids (VI), such as isoprene and monoterpenes, from plant leaves. Yet, since gaining scientific attention in the 1950’s, empirical research on leaf VI has been largely confined to laboratory experiments and atmospheric observations. Here, we introduce a new field instrument designed to bridge the scales from leaf to atmosphere, by enabling precision VI detection in real time from plants in their natural ecological setting. With a field campaign in the Brazilian Amazon, we reveal an unexpected distribution of leaf emission capacities (EC) across the vertical axis of the forest canopy, with EC peaking in the mid-canopy instead of the sun-exposed canopy surface, and moderately high emissions occurring in understory specialist species. Compared to the simple interpretation that VI protect leaves from heat stress at the hot canopy surface, our results encourage a more nuanced view of the adaptive role of VI in plants. We infer that forest emissions to the atmosphere depend on the dynamic microenvironments imposed by canopy structure, and not simply on canopy surface conditions. We provide a new emissions inventory from 52 tropical tree species, revealing moderate consistency in EC within taxonomic groups. We highlight priorities in leaf volatiles research that require field-portable detection systems. Our self-contained, portable instrument provides real-time detection and live measurement feedback with precision and detection limits better than 0.5 nmol VI m -2 leaf s -1 . We call the instrument ‘PORCO’ based on the gas detection method: photoionization of organic compounds. We provide a thorough validation of PORCO and demonstrate its capacity to detect ecologically driven variation in leaf emission rates and thus accelerate a nascent field of science: the ecology and ecophysiology of plant volatiles.

54 ENVIRONMENTAL SCIENCES↗