Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “sample complexity analysis”

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

Silver-mediated separations: A comprehensive review on advancements of argentation chromatography, facilitated transport membranes, and solid-phase extraction techniques and their applications

The use of silver(I) ions in chemical separations, also known as argentation separations, is a powerful approach for the selective separation and analysis of many natural and synthetic organic compounds. In this review, a comprehensive discussion of the most common argentation separation techniques, including argentation-liquid chromatography (Ag-LC), argentation-gas chromatography (Ag-GC), argentation-facilitated transport membranes (Ag-FTMs), and argentation-solid phase extraction (Ag-SPE) is provided. For each of these techniques, notable advancements, optimized separations, and innovative applications are discussed. The review begins with an explanation of the fundamental chemistry underlying argentation separations, mainly the reversible π-complexation between silver(I) ions and carbon-carbon double bonds. Within Ag-LC, the use of silver(I) ions in thin-layer chromatography, high-performance liquid chromatography, as well as preparative LC are explored. This discussion focuses on how silver(I) ions are employed in the stationary and mobile phase to separate unsaturated compounds. For Ag-GC and Ag-FTMs, different silver compounds and supporting media are discussed, often with relation to olefin-paraffin separations. Ag-SPE has been widely employed for the selective extraction of unsaturated compounds from complex matrices in sample preparation. This comprehensive review of Ag-LC, Ag-GC, Ag-FTMs, and Ag-SPE techniques emphasizes the immense potential of argentation separations in separations science and serves as a valuable resource for researchers seeking to learn, optimize, and utilize argentation separations.

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

A hybrid numerical and machine learning framework for evaluating the performance of a 780 cm 2 aqueous organic redox flow battery

Aqueous organic redox flow battery (AORFB) is a promising cost-competitive technology for large-scale energy storage. Among existing work, the dihydroxyphenazine (DHP)-based AORFB has demonstrated high energy density and low capacity degradation in 10 cm2 cells during lab tests. However, its commercial-scale performance in more complex environments remains unknown, posing a barrier for commercialization. To address this gap, this work presents a comprehensive performance evaluation of a 780 cm 2 DHP-based AORFB by combining physics-based numerical model, machine learning (ML)-based surrogate models, and ML-derived sensitivity quantification. Specifically, we first select 12 key battery parameters that include 10 physicochemical quantities and 2 operation quantities, then select 6 performance metrics that include energy efficiency (EE), discharging capacity, charging energy, and power losses due to concentration, activation, and ohmic over-potentials. With such selection, 12800 combinations of the 12 parameters are subsequently generated using the Latin Hypercube Sampling method. These combinations, together with 38 pre-defined State of Charge, are then integrated to a validated AORFB model developed in COMSOL to compute the performance metrics. With both input parameters and performance metrics, 60 deep neural network (DNN) surrogate models are then trained to approximate the relationship between the 10 physicochemical quantities and 6 performance metrics at each flow rate and current density. Sensitivity scores are then calculated based on the DNN models. Two additional sensitivity analysis tools, i.e., MARS, and SHAP, are also used to cross-validate the sensitivity scores from the DNN. The results demonstrate that 1) the standard potential ranks the first in controlling EE and charging energy, 2) the membrane conductivity is most critical for power loss and EE, and 3) specific area and reaction rate control activation power loss.

25 ENERGY STORAGE↗
DOE OSTI2024

Evaluating large scale aqueous organic redox flow battery performance with a hybrid numerical and machine learning framework

Aqueous organic redox flow battery (AORFB) is a promising cost-competitive technology for large-scale energy storage. Among existing work, the dihydroxyphenazine (DHP)-based AORFB has demonstrated high energy density and low-capacity degradation in 10 cm$^2$ cells during lab tests. However, its commercial-scale performance in more complex environments remains unknown, posing a barrier to commercialization. To address this gap, this work presents a comprehensive performance evaluation of a 780 cm$^2$ DHP-based AORFB by combining a physics-based numerical model, machine learning (ML)-based surrogate models, and ML-derived sensitivity quantification. Specifically, we first select 12 key battery parameters that include 10 physicochemical and 2 operation quantities, then select 6 performance metrics that include energy efficiency (EE), discharging capacity, charging energy, and power losses due to concentration, activation, and ohmic over-potentials. With such selection, 12800 combinations of the 12 parameters are subsequently generated using the Latin Hypercube Sampling method. Such combinations, together with 38 pre-defined State of Charge, are then integrated to a validated AORFB model developed in COMSOL to compute the performance metrics. With both input parameters and performance metrics, 60 deep neural network (DNN) surrogate models are then trained to approximate the relationship between the 10 physicochemical quantities and 6 performance metrics at each flow rate and current density. Sensitivity scores are then calculated based on the DNN models. Two additional sensitivity analysis tools, i.e., MARS, and SHAP, are also used to cross-validate the sensitivity scores from the DNN. The results demonstrate that 1) the standard potential ranks first in controlling EE and charging energy, 2) the membrane conductivity is most critical for power loss and EE, and 3) specific area and reaction rate control activation power loss.

25 ENERGY STORAGE↗
DOE OSTI2015

Towards intelligent diagnostic system employing integration of mathematical and engineering model

The development of medical diagnostic system has been one of the main research fields during years. The goal of the medical diagnostic system is to place a nosological system that could ease the diagnostic evaluation normally performed by scientists and doctors. Efficient diagnostic evaluation is essentials and requires broad knowledge in order to improve conventional diagnostic system. Several approaches on developing the medical diagnostic system have been designed and tested since the earliest 60s. Attempts on improving their performance have been made which utilizes the fields of artificial intelligence, statistical analyses, mathematical model and engineering theories. With the availability of the microcomputer and software development as well as the promising aforementioned fields, medical diagnostic prototypes could be developed. In general, the medical diagnostic system consists of several stages, namely the 1) data acquisition, 2) feature extraction, 3) feature selection, and 4) classifications stages. Data acquisition stage plays an important role in converting the inputs measured from the real world physical conditions to the digital numeric values that can be manipulated by the computer system. One of the common medical inputs could be medical microscopic images, radiographic images, magnetic resonance image (MRI) as well as medical signals such as electrocardiogram (ECG) and electroencephalogram (EEG). Normally, the scientist or doctors have to deal with myriad of data and redundant to be processed. In order to reduce the complexity of the diagnosis process, only the significant features of the raw data such as peak value of the ECG signal or size of lesion in the mammogram images will be extracted and considered in the subsequent stages. Mathematical models and statistical analyses will be performed to select the most significant features to be classified. The statistical analyses such as principal component analysis and discriminant analysis as well as mathematical model of clustering technique have been widely used in developing the medical diagnostic systems. The selected features will be classified using mathematical models that embedded engineering theory such as artificial intelligence, support vector machine, neural network and fuzzy-neuro system. These classifiers will provide the diagnostic results without human intervention. Among many publishable researches, several prototypes have been developed namely NeuralPap, Neural Mammo, and Cervix Kit. The former system (NeuralPap) is an automatic intelligent diagnostic system for classifying and distinguishing between the normal and cervical cancerous cells. Meanwhile, the Cervix Kit is a portable Field-programmable gate array (FPGA)-based cervical diagnostic kit that could automatically diagnose the cancerous cell based on the images obtained during sampling test. Besides the cervical diagnostic system, the Neural Mammo system is developed to specifically aid the diagnosis of breast cancer using a fine needle aspiration image.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2017

Internal Rotation in the Globular Cluster M53

We present an analysis of the internal bulk rotation in the metal-poor globular cluster (GC) NGC 5024 (M53) using radial velocities (RVs) of individual cluster members. We use RV measurements from a previous abundance study of M53 done using the Hydra multi-object spectrograph on the WIYN 3.5 m telescope. The Hydra sample greatly increases the number of RVs available in the central regions of the cluster where the internal rotation is the strongest. The sample of cluster members is further increased through two previous kinematic studies of M53. The combined total sample contains 245 cluster members. With our sample, we are able to create a velocity dispersion profile of the cluster and derive a central velocity dispersion σ{sub 0}=4.0±0.3 km s{sup −1}; we find that M53 inner regions are characterized by a peak amplitude of rotation equal to 1.4±0.1 km s{sup −1} corresponding to a relatively high value of the ratio of the rotation speed to central velocity dispersion (V{sub rot}/σ{sub 0}=0.35±0.04). Our data also reveal a radial variation in the orientation of the projected rotation axis suggesting complex internal kinematics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Particle-image Velocimetry Sensitivity through Automatic Differentiation

A commonly encountered difficulty in modern experimental fluid dynamics experiments is the need for full-flow visualization and measurement; traditional techniques such as hot-wire probes and laser Doppler velocimetry cannot provide such information without difficult scanning procedures. Particle-image velocimetry (PIV) is a measurement technique which can overcome this spatial limitation to provide real-time three-dimensional fluid velocity vectors for an entire plane within the flow, as opposed to a single point that is provided by other measurement techniques. To accomplish this, a PIV system seeds the fluid flow with tracer particles; these must have nearly the same density as that of the fluid and a high drag coefficient to ensure that they follow the fluid flow without deviation. A laser is sent through optics to expand its beam into a sheet which illuminates the plane of interest in the fluid. The particles reflect more light than the fluid, and their image is recorded using a high speed camera. Due to limitations in bandwidth, video cannot be taken fast enough to be useful, and instead pairs of images are captured with a known separation in time. The illuminated particles will move between the two images taken by the camera, and their relative motion can be used to compute the fluid velocity at that point. When a single camera is used, only two-dimensional velocity vectors can be computed, but with the addition of another camera this can be improved to full three-dimensional vectors. The flow's time evolution can be measured by taking series of image pairs and processing each pair to get the velocity field at the time between each image pair. Only the two-dimensional version will be considered in this work. To compute the velocity vectors of the flow, the images are broken up into interrogation regions, usually in a uniform grid. In the most basic PIV algorithm, corresponding regions in the image pairs are cross-correlated to compute the most likely displacement from one image to the next. This is used since the cross correlation will achieve a maximum value when the two interrogation regions line up the best. The location of the peak in the correlation map and scaling information for the optical setup is used to compute the particle displacement from one image to the next. This process is similar to holding two semi-transparent images and moving one until they match as well as possible. Sub-pixel displacements between images can be obtained by fitting a Gaussian function to the correlation peak. Knowing the displacement vector for each set of interrogation regions in the image pair, it is a simple matter to divide by the time gap between the two images to compute the fluid velocity at the point in time and space. Since there are many interrogation regions in each image pair, the velocity field for a sheet is sampled according to the arrangement of the interrogation regions. The sensitivity of the final computed vector to the input images is difficult to compute, and has led most considerations of uncertainty in the PIV field to apply Monte Carlo methods and correlations to quantify the systematic and random errors which result from application of various PIV algorithms. In this work, automatic differentiation (AD) has been used to trace the dependence of velocity vectors computed using the correlation-peak PIV algorithm back to the input values given to an artificial image generator. These initial results form the foundation for continued work in the application of AD to algorithms in PIV. They demonstrate that AD methods can be successfully applied to PIV algorithms without undue complexity. The derivatives of each individual vector with regard to the coefficients of the velocity field used to create the artificial images can be computed and the statistics of them considered collectively, or correlated between each other to quantify the effectiveness of a PIV algorithm. With the initial implementation complete, adding extra algorithms is relatively simple and considerable data can be generated for analysis. It is hoped that application of these methods can lead to improved PIV algorithms which demonstrate better sensitivity to small variations in particle displacement. While these methods can be used as an analysis tool for refinement of new algorithms, computation of actual resultant vector sensitivities to input images' individual pixels is the next logical step of this work. Such information will be analyzed with experimental data to produce a better understanding of how output vectors' uncertainty relates to the actual image pairs, and lead to automatic bad vector removal. (authors)

42 ENGINEERING↗
DOE OSTI2015

Mineralogy, provenance, and diagenesis of a potassic basaltic sandstone on Mars: CheMin X-ray diffraction of the Windjana sample (Kimberley area, Gale Crater)

The Windjana drill sample, a sandstone of the Dillinger member (Kimberley formation, Gale Crater, Mars), was analyzed by CheMin X-ray diffraction (XRD) in the MSL Curiosity rover. From Rietveld refinements of its XRD pattern, Windjana contains the following: sanidine (21% weight, ~Or 95 ); augite (20%); magnetite (12%); pigeonite; olivine; plagioclase; amorphous and smectitic material (~25%); and percent levels of others including ilmenite, fluorapatite, and bassanite. From mass balance on the Alpha Proton X-ray Spectrometer (APXS) chemical analysis, the amorphous material is Fe rich with nearly no other cations—like ferrihydrite. The Windjana sample shows little alteration and was likely cemented by its magnetite and ferrihydrite. From ChemCam Laser-Induced Breakdown Spectrometer (LIBS) chemical analyses, Windjana is representative of the Dillinger and Mount Remarkable members of the Kimberley formation. LIBS data suggest that the Kimberley sediments include at least three chemical components. The most K-rich targets have 5.6% K 2 O, ~1.8 times that of Windjana, implying a sediment component with >40% sanidine, e.g., a trachyte. A second component is rich in mafic minerals, with little feldspar (like a shergottite). A third component is richer in plagioclase and in Na 2 O, and is likely to be basaltic. The K-rich sediment component is consistent with APXS and ChemCam observations of K-rich rocks elsewhere in Gale Crater. The source of this sediment component was likely volcanic. Finally, the presence of sediment from many igneous sources, in concert with Curiosity's identifications of other igneous materials (e.g., mugearite), implies that the northern rim of Gale Crater exposes a diverse igneous complex, at least as diverse as that found in similar-age terranes on Earth.

58 GEOSCIENCES↗
DOE OSTI2021

Manipulating meso-scale solvent structure from Pd nanoparticle deposits in deep eutectic solvents

Deep Eutectic Solvents (DESs) are complex solutions that present unique challenges compared to traditional electrolytes. Unlike most electrolytes and ionic liquids, DESs have delicate hydrogen bond networks that are responsible for their highly sensitive compositional dependence on the melting point. Prior work has demonstrated unique nanoscale structure both experimentally and theoretically that brings both challenges and opportunities to their adoption in traditional electrochemical processes. In this study, we use in-situ sample-rotated ultra-small angle X-ray scattering to resolve the near-interface solvent structure after electrodepositing Pd nanoparticles onto a glassy carbon electrode in choline chloride:urea and choline chloride:ethylene glycol DESs. Our results indicate that solvent structure can be observed on the meso-scale in the choline chloride:urea and choline chloride:ethylene glycol systems. Importantly, this extended solvent structure peaks at -0.5 V (vs. Ag/AgCl) and decreases slightly with higher overpotentials during Pd electrodeposition. Experimentally, the nature of this structure is more pronounced in the ethylene glycol system, as evidenced by both the X-ray scattering as well as the electrochemical impedance spectroscopy (EIS). Here, molecular dynamics simulations and dipolar orientation analysis reveal that chloride delocalization near the Pd interface and long-range interactions between the choline and each hydrogen bond donor are very different and qualitatively consistent with the experimental data. These results show how the long-range solvent-deposit interactions can be tuned by changing the hydrogen bond donor in the DES and the applied potential.

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

Modeling thin layers of analytes on substrates for spectral analysis: use of solid/liquid n and k values to model reflectance spectra

Since solids are only sometimes seen en masse in a pure bulk form, and for liquids other than water almost never, a capability to model reflectance spectra from analytes deposited on various substrates would be highly advantageous. If available, the real, n ( ν ) , and imaginary, k ( ν ) , components of the complex refractive index, n n + i k , can be used to simulate infrared spectra, accounting for reflection, refraction, and absorption phenomena as a function of wavelength. We focus on using the Pacific Northwest National Laboratory (PNNL) derived n / k vectors for solid and liquid analytes deposited as thin layers on different types of substrates including conductors, such as aluminum, and inorganic dielectrics, such as glass. The model is an adaptation of the Monte Carlo ray trace modeling program, TracePro, extended through the use of its macrolanguage. The model is tested using thin films of organic liquids including silicone oil and no. 2 diesel fuel, as well as organic solids such as caffeine and acetaminophen on aluminum and glass. The predicted spectra for the solid films were compared to experimental hemispherical reflectance data measured using a Fourier transform spectrometer with an integrating sphere. The thickness of the calculated layer is a parameter for predicting the (transflectance) spectra and is obtained using the areal density measured from gravimetric methods to generate the thin-layer samples. Comparison of the calculated spectra with experimental hemispherical reflectance data shows excellent agreement, indicating promise for the use of measured n / k data to synthesize reference spectral data. In particular, accounting for the inhomogeneity of the deposits greatly improved the match with experimental data. Finally, the theoretical modeling shows that for thicker layers (ca. 20 to 100 μm) of typical organics possessing moderately strong k values, the longwave infrared features are often saturated and better spectral contrast is obtained from the overtone/combination bands in the shortwave infrared.

42 ENGINEERING↗
DOE OSTI2020

Retrospective on Recent DOE-Funded Studies Concerning the Extraction of Rare Earth Elements & Lithium from Geothermal Brines (Final Report)

Rare earth elements (REE) and lithium are non-toxic metals that are considered critical materials due to their use in electronics, magnets, batteries, and a wide variety of industrial processes important for the economy and military preparedness. Demand for REE and lithium is increasing and these critical materials are imported, so identifying and exploiting domestic sources of REE and lithium is a national priority. The U.S. Department of Energy (DOE) Geothermal Technologies Office (GTO) has been in the forefront of sponsoring research investigating the potential recovery of REE, lithium, and other critical minerals from geothermal brines. It has been proposed that the future of geothermal energy should include “hybrid systems” that combine electricity generation with other revenue-generating activities, such as recovery of valuable and critical minerals, including REE and lithium. Two recent GTO funding opportunities have focused on the recovery of REE and other valuable minerals from geothermal brines. The research supported by the GTO’s mineral recovery program is focused on three areas: resource characterization, technology for the extraction of REE, and technology for the extraction of lithium (Tables 1 and 2). This report is a retrospective study examining the outcome of GTO’s two recent mineral recovery programs (DE-FOA-0001016 in FY 2014 and DE-FOA-0001376 in FY 2016). In this report, the knowledge, technology, and techniques that were developed by researchers funded by GTO are summarized and discussed. Four projects were funded to assess the concentrations and amounts of REE found in geothermal brines and oil field produced waters. The GTO-funded studies compiled publically available data on REE concentrations from brines and produced water from all over the USA. In addition, new samples were collected and characterized from major geothermal and hydrocarbon basins in the Western USA. The studies examined the relationship between lithology and REE concentrations and developed models examining the influence of geology on REE concentrations in produced brines. It was determined that REE are frequently found at higher concentrations in oil field produced water than geothermal brines, but that some geothermal areas had significant REE resources. Significant reservoirs of REE were identified in the Western USA. In some cases, concentrations of REE were more than 1000 times the concentrations found in seawater. Collectively, these studies represent a comprehensive picture of REE resources associated with geothermal and hydrocarbon systems in the USA. The studies did not examine lithium resources, but in some cases, lithium concentration data was collected. Data from these studies are housed in the Geothermal Data Repository (GDR) and represent a significant information resource and it is recommended that these data be further analyzed in a future study. Eight projects were funded to develop new technology for REE extraction from geothermal fluids. These projects investigated sorption as an approach for removal and recovery of REE from geothermal brines. The projects investigated cutting-edge technology for selective sorption of ions from complex solutions, including the application of metal-organic frameworks and biosorbent proteins. The REE sorption studies tested different combinations of metal-binding ligands and solid supports. The most promising metal-binding ligands for REE included phosphonic acid, thiol, and carboxylic acid functional groups. Ligands were attached or incorporated into a wide variety of solid supports. In most cases, attachment was via covalent bonding to organic resins, polymers, or silica-based supports. Most of the REE projects were conducted at a low technology readiness level (TRL) and showed promise, but direct comparison between technologies was not possible based on the available information. It is recommended that testing and reporting be standardized to the extent possible to facilitate comparisons between technologies. Two projects were directed at novel lithium extraction technology. Both projects investigated the use of inorganic sorbents, including manganese oxides. One study also examined the use of metal- ion imprinted polymers as selective ion-exchange resins for the separation of lithium and manganese from brines. Both approaches showed promise for the selective extraction of lithium from brines, including potentially geothermal brines. Results from these GTO studies indicated that selective REE and lithium extraction is possible, but interference from co-occurring solutes, such as calcium, magnesium, or heavy metals, will interfere with process efficiency and negatively impact process economics. Techno-economic analysis conducted as part of the resource and technology studies suggest extraction of REE from geothermal brines is unlikely to be economically viable, especially since non-geothermal produced waters frequently have higher REE concentrations. It is recommended that benchmarks for techno-economic analysis be established to the extent possible for future studies, to facilitate direct comparison of various technologies. Based on the collective results of this program, it appears that hybrid geothermal power would benefit more from recovery of lithium and other metals, rather than REE. It is recommended that future studies be conducted at a higher- TRL and that sorbents be tested against actual geothermal fluid samples. Prior higher-TRL efforts to extract metals from geothermal brines should be further evaluated for lessons learned.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Effect of D2O on growth properties and chemical structure of annual ryegrass (Lolium multiflorum)

In present paper, we report the production and detailed structural analysis of deuterium-enriched rye grass (Lolium multiflorum) for neutron scattering experiments. An efficient method to produce deuterated biomass was developed by designing hydroponic perfusion chambers. In preliminary studies, the partial deuterated rye samples were grown in increasing levels of D2O to study the seed germination and the level of deuterium incorporation as a function of D2O concentration. Solution NMR method indicated 36.9 % deuterium incorporation in 50 % D2O grown annual rye samples and further significant increase in the deuterium incorporation level was observed by germinating the rye seedlings in H2O and growing in 50 % D2O inside the perfusion chambers. Moreover, in an effort to compare the substrate characteristics related to enzymatic hydrolysis on deuterated and protiated version of biomass, annual rye grown in 50 % D2O was selected for detailed biomass characterization studies. The compositional analyses, degree of polymerization and cellulose crystallinity were compared with its protiated control. The cellulose molecular weight indicated slight variation with deuteration; however, hemicellulose molecular weights and cellulose crystallinity remain unaffected with the deuteration. Besides the minor differences in biomass components, the development of deuterated biomass for neutron scattering application is essential to understand the complex biomass conversion processes.

deuteration↗
DOE OSTI2014

Quantitative Analysis of Human Salivary Gland-Derived Intact Proteome Using Top-Down Mass Spectrometry

There are several notable challenges inherent to fully characterizing the entirety of the human saliva proteome using bottom-up approaches, including polymorphic isoforms, post-translational modifications, unique splice variants, deletions, and truncations. To address these challenges, we have developed a top-down based liquid chromatography-mass spectrometry (LC-MS) approach, which cataloged 20 major human salivary proteins with a total of 83 proteoforms, containing a broad range of post-translational modifications. Among these proteins, several previously reported disease biomarker proteins were identified at the intact protein level, such as beta-2 microglobulin (B2M). In addition, intact glycosylated proteoforms of several saliva proteins were also characterized, including intact N-glycosylated protein prolactin inducible protein (PIP) and O-glycosylated acidic protein rich protein (aPRP). These characterized proteoforms constitute an intact saliva proteoform database, which was used for quantitative comparison of intact salivary proteoforms among six healthy individuals. Human parotid (PS) and submandibular/sublingual gland (SMSL) secretion samples (2 μg of protein each) from six healthy individuals were compared using RPLC coupled with the 12T FTICR mass spectrometer. Significantly different protein and PTM patterns were resolved with high reproducibility between PS and SMSL glands. The results from this study provide further insight into the potential mechanisms of PTM pathways in oral glandular secretion, expanding our knowledge of this complex yet easily accessible fluid. Intact protein LC-MS approach presented herein can potentially be applied for rapid and accurate identification of biomarkers from only a few microliters of human glandular saliva.

Environmental Molecular Sciences Laboratory↗
DOE OSTI2018

Thermodynamics of Reactions Affected by Medium Reorganization

We present a thermodynamic analysis of the activation barrier for reactions which can be monitored through the difference in the energies of reactants and products defined as the reaction coordinate (electron and atom transfer, enzyme catalysis, etc). The free-energy surfaces along the reaction coordinate are separated into the enthalpy and entropy surfaces. For the Gaussian statistics of the reaction coordinate, the free-energy surfaces are parabolas, and the entropy surface is an inverted parabola. Its maximum coincides with the transition state for reactions with zero value of the reaction free energy. Maximum entropic depression of the activation barrier, anticipated by the concept of transition-state ensembles, can be achieved for such reactions. From Onsager’s reversibility, the entropy of equilibrium fluctuations encodes the entropic component of the activation barrier. The reorganization entropy thus becomes the critical parameter of the theory reducing the problem of activation entropy to the problem of reorganization entropy. Standard solvation theories do not provide reorganization entropy sufficient for the barrier depression. Complex media, characterized by many relaxation processes, need to be involved and proteins provide several routes for achieving large entropic effects related to incomplete (nonergodic) sampling of the complex energy landscape and active role of water in the reaction mechanism.

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

Parametric Sensitivity and Uncertainty Quantification in the Version 1 of E3SM Atmosphere Model Based on Short Perturbed Parameter Ensemble Simulations

Abstract The atmospheric component of Energy Exascale Earth System Model version 1 has included many new features in the physics parameterizations compared to its predecessors. Potential complex nonlinear interactions among the new features create a significant challenge for understanding the model behaviors and parameter tuning. Using the one‐at‐a‐time method, the benefit of tuning one parameter may offset the benefit of tuning another parameter, or improvement in one target variable may lead to degradation in another target variable. To better understand the Energy Exascale Earth System Model version 1 model behaviors and physics, we conducted a large number of short simulations (three days) in which 18 parameters carefully selected from parameterizations of deep convection, shallow convection, and cloud macrophysics and microphysics were perturbed simultaneously using the Latin hypercube sampling method. From the perturbed parameter ensemble simulations and use of different skill score functions, we identified the most sensitive parameters, quantified how the model responds to changes of the parameters for both global mean and spatial distribution, and estimated the maximum likelihood of model parameter space for a number of important fidelity metrics. Comparison of the parametric sensitivity using simulations of two different lengths suggests that perturbed parameter ensemble using short simulations has some bearing on understanding parametric sensitivity of longer simulations. Results from this analysis provide a more comprehensive picture of the Energy Exascale Earth System Model version 1 behavior. The difficulty in reducing biases in multiple variables simultaneously highlights the need of characterizing model structural uncertainty (so‐called embedded errors) to inform future development efforts.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Tritium Transport Modeling at the Pore Scale in Ceramic Breeder Materials Using TMAP8

Fusion reactors depend on the blanket material to breed and release tritium at the same rate or faster than it is consumed by the fusion reaction. Cellular ceramic breeders (CCBs) are dense materials that can maintain a high tritium breeding ratio while promoting tritium release because of highly connected pores. Assessing the tritium breeding capabilities of these materials requires a combination of extensive experimental and modeling efforts. In this work, we develop and calibrate a multiphysics model of tritium transport. This novel model accounts for ceramic and pore diffusion, trapping and detrapping, and several surface reactions at the pore surface. We perform a sensitivity analysis and calibrate the model by comparing its predictions against experimental measurements of deuterium absorption. The calibrated model is then used to model tritium absorption in samples with different pore microstructures to investigate the effect of pore interconnectivity on tritium absorption. The model is part of the development of the multiscale, multiphysics framework for tritium transport [i.e., the Tritium Migration Analysis Program (TMAP8)], which is itself built on top of the finite-element multiphysics framework multiphysics object-oriented simulation environment (MOOSE). This study demonstrates some of TMAP8’s capabilities and is the first step toward assessing the tritium breeding capabilities of ceramic breeder material designs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Measurement Error and Resolution in Quantitative Stable Isotope Probing: Implications for Experimental Design

Quantitative stable isotope probing (qSIP) estimates isotope tracer incorporation into DNA of individual microbes and can link microbial biodiversity and biogeochemistry in complex communities. As with any quantitative estimation technique, qSIP involves measurement error, and a fuller understanding of error, precision, and statistical power benefits qSIP experimental design and data interpretation. We used several qSIP data sets—from soil and seawater microbiomes—to evaluate how variance in isotope incorporation estimates depends on organism abundance and resolution of the density fractionation scheme. We assessed statistical power for replicated qSIP studies, plus sensitivity and specificity for unreplicated designs. As a taxon’s abundance increases, the variance of its weighted mean density declines. Nine fractions appear to be a reasonable trade-off between cost and precision for most qSIP applications. Increasing the number of density fractions beyond that reduces variance, although the magnitude of this benefit declines with additional fractions. Our analysis suggests that, if a taxon has an isotope enrichment of 10 atom% excess, there is a 60% chance that this will be detected as significantly different from zero (with alpha 0.1). With five replicates, isotope enrichment of 5 atom% could be detected with power (0.6) and alpha (0.1). Finally, we illustrate the importance of internal standards, which can help to calibrate per sample conversions of %GC to mean weighted density. These results should benefit researchers designing future SIP experiments and provide a useful reference for metagenomic SIP applications where both financial and computational limitations constrain experimental scope.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Non-destructive simulation of node defects in additively manufactured lattice structures

Additive Manufacturing (AM), commonly referred to as 3D printing, offers the ability to not only fabricate geometrically complex lattice structures but parts in which lattice topologies in-fill volumes bounded by complex surface geometries. However, current AM processes produce defects on the strut and node elements which make up the lattice structure. This creates an inherent difference between the as-designed and as-fabricated geometries, which negatively affects predictions (via numerical simulation) of the lattice’s mechanical performance. Although experimental and numerical analysis of an AM lattice’s bulk structure, unit cell and struts have been performed, there exists almost no research data on the mechanical response of the individual as-manufactured lattice node elements. Here we propose a methodology that, for the first time, allows non-destructive quantification of the mechanical response of node elements within an as-manufactured lattice structure. A custom-developed tool is used to extract and classify each individual node geometry from micro-computed tomography scans of an AM fabricated lattice. Voxel-based finite element meshes are generated for numerical simulation and the mechanical response distribution is compared to that of the idealised computer-aided design model. The method demonstrates compatibility with Uncertainty Quantification methods that provide opportunities for efficient prediction of a population of nodal responses from sampled data. Overall, the non-destructive and automated nature of the node extraction and response evaluation is promising for its application in qualification and certification of additively manufactured lattice structures.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Position-specific carbon isotope analysis of serine by gas chromatography/Orbitrap mass spectrometry, and an application to plant metabolism

Position-specific 13 C/ 12 C ratios within amino acids remain largely unexplored in environmental samples due to methodological limitations. We hypothesized that natural-abundance isotope patterns in serine may serve as a proxy for plant metabolic fluxes including photorespiration. Here we describe an Orbitrap method optimized for the position-specific carbon isotope analysis of serine to test our hypothesis and discuss the generalizability of this method to other amino acids. Position-specific carbon isotope ratios of serine were measured using a Thermo Scientific™ Q Exactive™ GC Orbitrap™. Amino acids were hydrolyzed from Arabidopsis biomass, purified from potential matrix interferences, and derivatized alongside standards. Derivatized serine (N,O-bis(trifluoroacetyl)methyl ester) was isolated using gas chromatography, trapped in a reservoir, and purged into the electron ionization source over tens of minutes, producing fragment ions containing different combinations of atoms from the serine-derivative molecule. The 13 C/ 12 C ratios of fragments with monoisotopic masses of 110.0217, 138.0166, and 165.0037 Da were monitored in the mass analyzer and used to calculate position-specific δ 13 C values relative to a working standard. This methodology constrains position-specific δ 13 C values for nanomole amounts of serine isolated from chemically complex mixtures. The δ 13 C values of fragment ions of serine were characterized with ≤1‰ precisions, leading to propagated standard errors of 0.7–5‰ for each carbon position. Position-specific δ 13 C values differed by up to ca 28 ± 5‰ between serine molecules hydrolyzed from plants grown under contrasting pCO 2 , selected to promote different fluxes through photosynthesis and photorespiration. The method was validated using pure serine standards characterized offline. Here this study presents the first Orbitrap-based measurements of natural-abundance, position-specific carbon isotope variation in an amino acid isolated from a biological matrix. We present a method for the precise characterization of isotope ratios in serine and propose applications probing metabolism in plants. We discuss the potential for extending these approaches to other amino acids, paving the way for novel applications.

59 BASIC BIOLOGICAL SCIENCES↗