Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “total variation”

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 433 records · Page 24
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Electric and magnetic dipole strength in 112,114,116,118,120,124 Sn

Background: There is renewed interest in electric dipole strength distributions for a variety of reasons including the extraction of the dipole polarizability related to properties of the symmetry energy and a measure for the neutron skin thickness, understanding the structure of low-energy $\textit{E}$1 strength in nuclei with neutron excess, and establishing the systematics of the isovector giant dipole resonance (IVGDR). Inelastic proton scattering at energies of a few hundred MeV and very forward angles including 0° has been established as a tool for the study of electric and magnetic dipole strength distributions in nuclei. Purpose: The present work aims at a systematic investigation of the electric and magnetic dipole strength distributions in the chain of stable even-mass tin isotopes. Methods: Inelastic proton scattering experiments were performed at the Research Center for Nuclear Physics, Osaka, with a 295-MeV beam covering laboratory angles 0°–6° and excitation energies 6–22 MeV. Cross sections due to $\textit{E}$1 and $\textit{M}$1 excitations were extracted with a multipole decomposition analysis (MDA) and then converted to reduced transition probabilities with the “virtual photon method” for $\textit{E}$1 and the “unit cross section method” for $\textit{M}$1 excitations, respectively. Including a theory-aided correction for the high-excitation-energy region not covered experimentally, the electric dipole polarizability was determined from the $\textit{E}$1 strength distributions. Results: In this work, total photoabsorption cross sections derived from the $\textit{E}$1 and $\textit{M}$1 strength distributions show significant differences compared to those from previous ($\textit{γ, xn}$) experiments in the energy region of the IVGDR. The widths of the IVGDR deduced from the present data with a Lorentz parametrization show an approximately constant value of about 4.5 MeV in contrast to the large variations between isotopes observed in previous work. The IVGDR centroid energies are in good correspondence to expectations from empirical systematics of their mass dependence. Furthermore, a study of the dependence of the IVGDR energies on bulk matter properties is presented. The $\textit{E}$1 strengths below neutron threshold show fair agreement in this work with results from ($\textit{γ, γ'}$) experiments on 112,116,120,124 Sn in the energy region between 6 and 7 MeV, where also isoscalar $\textit{E}$1 strength was found for 124 Sn. At higher excitation energies, large differences are observed, pointing to a different nature of the excited states with small ground-state branching ratios. The isovector spin-$\textit{M}$1 strengths exhibit a broad distribution between 6 and 12 MeV in all studied nuclei. Conclusions: The present results contribute to the solution of a variety of nuclear structure problems including the systematics of the energy and width of the IVGDR, the structure of low-energy $\textit{E}$1 strength in nuclei, new constraints to energy density functionals (EDFs) aiming at a systematic description of the dipole polarizability across the nuclear chart, from which properties of the symmetry energy can be derived, and the systematics of the isovector spin-$\textit{M}$1 strength in heavy nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Quantifying sensitivity in numerical weather prediction-modeled offshore wind speeds through an ensemble modeling approach

A decade of research has shown that numerical weather prediction (NWP)-modeled wind speeds can be highly sensitive to the inputs and setups within the NWP model. For wind resource characterization applications, this sensitivity is often addressed by constructing a range of setups and selecting the one that best validates against observations. However, this approach is not possible in areas that lack high-quality hub height observations, especially offshore wind areas. In such cases, techniques to quantify and disseminate confidence in NWP-modeled wind speeds in the absence of observations are needed. We address this need in the present study and propose best practices for quantifying the spread in NWP-modeled wind speeds. We implement an ensemble approach in which we consider 24 different setups to the Weather Research and Forecasting (WRF) model. We construct the ensemble by considering variations in WRF version, WRF namelist, atmospheric forcing, and sea surface temperature (SST) forcing. Our analysis finds that the standard deviation produces more consistent estimates compared to the interquartile range and tends to be the more conservative estimator for ensemble variability. We further find that model spread increases closer to the surface and on shorter time scales. In conclusion, we explore methods to attribute total ensemble variability to the different ensemble components (e.g., atmospheric forcing and SST product) and find that contributions by components also vary depending on time scale. We anticipate that the methods and results presented in this paper will provide a reasonable foundation for further research into ensemble-based wind resource data sets.

17 WIND ENERGY↗
DOE OSTI2022

Holocene paleoclimate change in the western US: The importance of chronology in discerning patterns and drivers

We report sediment in lakes and meadows forms a powerful archive that can be used to reconstruct environmental change through time. Reconstructions of lake level, of chemical, biological, and hydrological conditions, and of surrounding vegetation provide detailed information about past climate conditions, both locally and regionally. Indeed, most of our current knowledge of centennial- to millennial-scale climate variability in the arid western United States, where information about past hydroclimate is particularly important, comes from such sediment-based reconstructions. The pressing need for robust, precise predictions of future conditions is a significant motivation for paleoclimate science, and current research questions frequently require Holocene reconstructions to be resolved at sub-centennial timescales. Increasingly, regional syntheses seek to identify synoptic-scale patterns similar to those defined from modern observations (seasonal, interannual, multi-decadal, etc.) or to compare with the output of climate model simulations. However, the age control on existing records, especially those more than about 20 years old, is often sufficient only for millennial-scale interpretation. Here we assess the age control for 84 published and unpublished records from lakes and meadows in the Great Basin, California, and desert southwest, and use Bayesian modeling to evaluate the 95% uncertainty ranges for the 42 best-dated records. In the Late Holocene, about half of the 42 records have <400-year mean uncertainty ranges; however, high-precision age control is especially critical for young records, used to develop an accurate understanding of a proxy’s response to known climate variations. In the Middle Holocene, records vary from 400 to >800-year mean uncertainty and records of the Early Holocene have 600- to >1400-year mean uncertainty ranges. We find that the largest control on modeled uncertainties is dating density, with at least 2 dates/kyr being optimal and suggest obtaining “range-finder” dates at the onset of a study to better predict the total number of dates needed for an adequate age model. Such a density avoids a commonly observed phenomenon of significant peaks in uncertainty arising in gaps between age control points. Analysis of the uncertainties associated with proxy shifts reveal that more than half are >400 years. Although such large uncertainties currently prevent sub-centennial interpretations in most cases, increased dating density, strategic use of limited funds (including budgeting for a 2 date/kyr minimum at the proposal stage), construction of age-depth models with Bayesian methods, and critical evaluation of chronological uncertainty will shed light on past climate variability at finer timescales, enhancing our understanding of global and regional drivers of western U.S. climate.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Stochastic evaluation of fourth-order many-body perturbation energies

A scalable, stochastic algorithm evaluating the fourth-order many-body perturbation (MP4) correction to energy is proposed. Three hundred Goldstone diagrams representing the MP4 correction are computer generated and then converted into algebraic formulas expressed in terms of Green’s functions in real space and imaginary time. They are evaluated by the direct (i.e., non-Markov, non-Metropolis) Monte Carlo (MC) integration accelerated by the redundant-walker and control-variate algorithms. The resulting MC-MP4 method is efficiently parallelized and is shown to display O(n 5.3 ) size-dependence of cost, which is nearly two ranks lower than the O(n 7 ) dependence of the deterministic MP4 algorithm. Furthermore, it evaluates the MP4/aug-cc-pVDZ energy for benzene, naphthalene, phenanthrene, and corannulene with the statistical uncertainty of 10 mE h (1.1% of the total basis-set correlation energy), 38 mE h (2.6%), 110 mE h (5.5%), and 280 mE h (9.0%), respectively, after about 10 9 MC steps.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2026

Exploring composition mixing in kilonova ejecta with ray-by-ray simulations

Binary neutron star merger (BNSM) ejecta are considered a primary repository of r-process nucleosynthesis and a source of the observed heavy-element abundances. We implement composition mixing into ray-by-ray radiation-hydrodynamic simulations of BNSM ejecta, coupled with an online nuclear network (NN). We model mixing via a gradient-based mixing approximation that evolves simultaneously with the hydrodynamics. Here, we find that mixing occurs in regions where the electron fraction changes rapidly. While mixing smooths composition gradients in transition regions, it has a negligible impact on the heavy-element yields. This is because the primary r-process site (the equatorial ejecta) is initially homogeneous in free neutrons, leaving no strong gradients for mixing to act upon. In each angular ray, the abundances of the most produced elements are robust under mixing, while the less abundant ones are more affected. The total global abundances change only slightly from mixing, since each angular ray contributes its most abundant elements. Furthermore, the predicted kilonova light curves show only minor reddening, with differences below the detectability of state-of-the-art telescopes. In general, we do not observe significant effects from mixing in the time span of the r-process. Consequently, mixing only leads to minor variations in abundances and light curves in ray-by-ray simulations.

Explosive burning↗
DOE OSTI2024

Impact of Pressure-Dependent Interfacial Tension and Contact Angle on Capillary Heterogeneity Trapping of CO2 in Storage Aquifers

Summary Carbon dioxide (CO2) capillary trapping increases the total amount of CO2 that can be effectively immobilized in storage aquifers. This trapping, manifesting itself as accumulated CO2 columns at a continuum scale, is because of capillary threshold effects that occur below low-permeability barriers. Considering that capillary pressure is dictated by heterogeneous pore throat size, the trapped CO2 column height and associated CO2 saturation will vary spatially within a storage aquifer. This variation will be influenced by two pressure-dependent interfacial parameters—CO2/brine interfacial tension (IFT) and CO2/brine/rock contact angle. Our objective is to understand how the pressure dependence of these two parameters affects the heterogeneity of capillary trapped CO2 at a continuum scale. Our conceptual model is a 1D two-zone system with the upper zone being a flow barrier (low permeability) and the lower zone being a flow path (high permeability). The inputs to this model include microfacies-dependent capillary pressure vs. saturation curves and permeability values. The input capillary pressure curves were collected in the literature that represents carbonate microfacies (e.g., dolograinstone) in a prevalent formation in the Permian Basin. We then used the Leverett j-function to scale the capillary pressure curve for the two zones that are assigned with the same or different microfacies. During scaling, we considered the influence of pressure on both the IFT and contact angle of CO2/brine/dolomite systems. We varied the zone permeability contrast ratio from 2 to 50. We then assumed capillary gravity equilibriums and calculated the CO2 saturation buildup corresponding to various trapped CO2 column heights. The CO2 saturation buildup is defined as the CO2 saturation in the lower layer minus that in the upper one. We found that the saturation buildup can be doubled when varying pressure in a storage aquifer, after considering pressure-dependent IFT and contact angles. Thus, assuming these two parameters to be constant across such aquifers would cause large errors in the quantification of capillary trapping of CO2. The whole study demonstrates the importance of considering pressure-dependent interfacial properties in predicting the vertical distribution of capillary trapped CO2. It has important implications in developing a better understanding of leakage risks and consequent storage safety.

Engineering↗
DOE OSTI2021

Using TRMM Latent Heat as a Source to Estimate Convection Induced Gravity Wave Momentum Flux in the Lower Stratosphere

Due to a lack of observations, it is a big challenge to quantify the momentum flux transported from the troposphere into the lower stratosphere by gravity waves generated by convection. This limits our understanding of the dynamics of general circulation in the stratosphere and the realism of simulations of the Quasi Biennial Oscillation. In the past decade, some general circulation models have linked the momentum flux to the latent heating from subgrid scale convective precipitation with some success. However, there is still a large uncertainty in the sources of subgrid scale convection in these models. This study applies the links between the momentum flux and convective latent heating directly to the convective precipitation derived from 16 years of TRMM precipitation radar observations in the tropics and subtropics. The total and directional momentum flux at 100 hPa are derived at individual convective pixels by using the TRMM Spectrum Latent Heating product and large-scale wind profiles from the ERA-Interim dataset. For the first time, we are able to estimate the geographical distribution of momentum flux at 100 hPa in tropics and subtropics from the observed convective sources. The diurnal, seasonal, and interannual variations of the derived momentum flux are presented and discussed. These results could provide a reference to validate the dynamic coupling between the tropical troposphere and the stratosphere in the general circulation models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Genome-wide approaches delineate the additive, epistatic, and pleiotropic nature of variants controlling fatty acid composition in peanut ( Arachis hypogaea L.)

Abstract The fatty acid composition of seed oil is a major determinant of the flavor, shelf-life, and nutritional quality of peanuts. Major QTLs controlling high oil content, high oleic content, and low linoleic content have been characterized in several seed oil crop species. Here, we employ genome-wide association approaches on a recently genotyped collection of 787 plant introduction accessions in the USDA peanut core collection, plus selected improved cultivars, to discover markers associated with the natural variation in fatty acid composition, and to explain the genetic control of fatty acid composition in seed oils. Overall, 251 single nucleotide polymorphisms (SNPs) had significant trait associations with the measured fatty acid components. Twelve SNPs were associated with two or three different traits. Of these loci with apparent pleiotropic effects, 10 were associated with both oleic (C18:1) and linoleic acid (C18:2) content at different positions in the genome. In all 10 cases, the favorable allele had an opposite effect—increasing and lowering the concentration, respectively, of oleic and linoleic acid. The other traits with pleiotropic variant control were palmitic (C16:0), behenic (C22:0), lignoceric (C24:0), gadoleic (C20:1), total saturated, and total unsaturated fatty acid content. One hundred (100) of the significantly associated SNPs were located within 1000 kbp of 55 genes with fatty acid biosynthesis functional annotations. These genes encoded, among others: ACCase carboxyl transferase subunits, and several fatty acid synthase II enzymes. With the exception of gadoleic (C20:1) and lignoceric (C24:0) acid content, which occur at relatively low abundance in cultivated peanuts, all traits had significant SNP interactions exceeding a stringent Bonferroni threshold (α = 1%). We detected 7682 pairwise SNP interactions affecting the relative abundance of fatty acid components in the seed oil. Of these, 627 SNP pairs had at least one SNP within 1000 kbp of a gene with fatty acid biosynthesis functional annotation. We evaluated 168 candidate genes underlying these SNP interactions. Functional enrichment and protein-to-protein interactions supported significant interactions (P-value < 1.0E−16) among the genes evaluated. These results show the complex nature of the biology and genes underlying the variation in seed oil fatty acid composition and contribute to an improved genotype-to-phenotype map for fatty acid variation in peanut seed oil.

Otyama, Paul I.↗
DOE OSTI2024

The role of Li doping in layered/layered Na x Li y Ni 0.4 Fe 0.2 Mn 0.4 O 2 intergrowth electrodes for sodium ion batteries

Here, the layered NaTMO 2 (TM = Ni, Fe, Mn) materials with the O3-type structure are attractive as positive electrodes for sodium ion batteries because of their high theoretical capacity. Additionally, Li doping in these materials has been shown to offer substantial enhancements to their electrochemical properties by promoting the formation of intergrowth structures, which are combinations of specific phases. However, the mechanism by which the intergrowth modifies the electrochemical properties is often unclear. Systematic variation of Li content in Na x Li y Ni 0.4 Fe 0.2 Mn 0.4 O 2 (NFM-Li y ) was conducted to identify its role in structural modification and electrochemical performance. Li contents of 0.15 and greater generate a layered/layered Na-O3/Li-O’3 intergrowth structure. 7 Li and 23 Na nuclear magnetic resonance and x-ray absorption spectroscopy identify that when the total solubility for alkali ions in the layered structure is exceeded, Li continues to form the Li-O’3 phase while the excess Na forms residual sodium compounds such as Na 2 O. Higher Li content is associated with improved capacity retention in the initial cycles from the superior stability of the mechanically linked NaO3/Li-O’3 structure that suppresses the P3 to OP2 phase transition during charge. However, high Li contents are associated with increased rates of parasitic side reactions that reduce long-term cycling stability. These side reactions are associated with the instability of the cathode-electrolyte interphase, which can be partially mitigated by atomic layer deposition (ALD) coating with alumina, which significantly enhances the capacity retention and Coulombic efficiency over many cycles. Overall, we find that the layered/layered Na-O3/Li-O’3 intergrowth structure is able to provide structural stability and suppress undesired phase transformations but is overwhelmed by the increased reactivity of the surface if not protected by surface coating.

25 ENERGY STORAGE↗
DOE OSTI2020

Nitrogen Fertilization Restructured Spatial Patterns of Soil Organic Carbon and Total Nitrogen in Switchgrass and Gamagrass Croplands in Tennessee USA

Nitrogen (N) fertilizers can potentially alter spatial distribution of soil organic carbon (SOC) and total nitrogen (TN) concentrations in croplands such as switchgrass (SG: Panicum virgatum L.) and gamagrass (GG: Tripsacum dactyloides L.), but it remains unclear whether these effects are the same between crops and under different rates of fertilization. 13 C and 15 N are two important proxy measures of soil biogeochemistry, but they were rarely examined as to their spatial distributions in soil. Based on a three-year long fertilization experiment in Middle Tennessee, USA, the top mineral horizon soils (0–15 cm) were collected using a spatially explicit design within two 15-m 2 plots under three fertilization treatments in SG and GG croplands. A total of 288 samples were collected based on 12 plots and 24 samples in each plot. The fertilization treatments were no N input (NN), low N input (LN: 84 kg N ha -1 in urea) and high N input (HN: 168 kg N ha -1 in urea). The SOC, TN, SOC/TN (C: N), δ 13 C and δ 15 N were quantified and their within-plot variations and spatial distributions were achieved via descriptive and geostatistical methods. Results showed that SG generally displayed 10~120% higher plot-level variations in all variables than GG, and the plot-level variations were 20~77% higher in NN plots than LN and HN plots in SG but they were comparable in unfertilized and fertilized plots in GG. Relative to NN, LN and HN showed more significant surface trends and spatial structures in SOC and TN in both croplands, and the fertilization effect appeared more pronounced in SG. Spatial patterns in C: N, δ 13 C and δ 15 N were comparable among different fertilization treatments in both croplands. The descending within-plot variations were also identified among variables (SOC > TN > δ 15 N > C: N > δ 13 C). This study demonstrated that N fertilizations generally reduced the plot-level variance and simultaneously re-established spatial structures of SOC and TN in bioenergy croplands, which little varied with fertilization rate but was more responsive in switchgrass cropland.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

CODEX weak lensing mass catalogue and implications on the mass–richness relation

The COnstrain Dark Energy with X-ray clusters (CODEX) sample contains the largest flux limited sample of X-ray clusters at 0.35 < z < 0.65. It was selected from ROSAT data in the 10 000 square degrees of overlap with BOSS, mapping a total number of 2770 high-z galaxy clusters. We present here the full results of the CFHT CODEX programme on cluster mass measurement, including a reanalysis of CFHTLS Wide data, with 25 individual lensing-constrained cluster masses. We employ lensfit shape measurement and perform a conservative colour–space selection and weighting of background galaxies. Using the combination of shape noise and an analytic covariance for intrinsic variations of cluster profiles at fixed mass due to large-scale structure, miscentring, and variations in concentration and ellipticity, we determine the likelihood of the observed shear signal as a function of true mass for each cluster. We combine 25 individual cluster mass likelihoods in a Bayesian hierarchical scheme with the inclusion of optical and X-ray selection functions to derive constraints on the slope α, normalization β, and scatter σln λ|μ of our richness–mass scaling relation model in log-space: ${\langle {\rm In}\,\, \lambda\!\!\mid\!\!\mu\rangle = \alpha\mu + \beta,} $ with μ = ln (M 200c /M piv ), and M piv = 10 14.81 M ⊙ . We find a slope $\alpha = 0.49^{+0.20}_{-0.15}$, normalization $\exp (\beta) = 84.0^{+9.2}_{-14.8}$, and $\sigma _{\ln \lambda | \mu } = 0.17^{+0.13}_{-0.09}$ using CFHT richness estimates. In comparison to other weak lensing richness–mass relations, we find the normalization of the richness statistically agreeing with the normalization of other scaling relations from a broad redshift range (0.0 < z < 0.65) and with different cluster selection (X-ray, Sunyaev–Zeldovich, and optical).

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

SDSS-IV MaNGA: Understanding Ionized Gas Turbulence Using Integral Field Spectroscopy of 4500 Star-forming Disk Galaxies

The Sloan Digital Sky Survey MaNGA program has now obtained integral field spectroscopy for over 10,000 galaxies in the nearby universe. We use the final MaNGA data release DR17 to study the correlation between ionized gas velocity dispersion and galactic star formation rate, finding a tight correlation in which σ Hα from galactic H ιι regions increases significantly from ~18–30 km s –1 , broadly in keeping with previous studies. In contrast, σ Hα from diffuse ionized gas increases more rapidly from 20–60 km s –1 . Using the statistical power of MaNGA, we investigate these correlations in greater detail using multiple emission lines and determine that the observed correlation of σ Hα with local star formation rate surface density is driven primarily by the global relation of increasing velocity dispersion at higher total star formation rate, as are apparent correlations with stellar mass. Assuming H ιι region models consistent with our finding that σ [O III] < σ Hα < σ [O I] , we estimate the velocity dispersion of the molecular gas in which the individual H ιι regions are embedded, finding values σ Mol = 5–30 km s –1 consistent with ALMA observations in a similar mass range. Finally, we use variations in the relation with inclination and disk azimuthal angle to constrain the velocity dispersion ellipsoid of the ionized gas σ z /σ r = 0.84 ± 0.03 and σ Φ /σ r = 0.91 ± 0.03, similar to that of young stars in the Galactic disk. Our results are most consistent with the theoretical models in which turbulence in modern galactic disks is driven primarily by star formation feedback.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2026

A variational mimetic finite difference method for elliptic interface problems on non-matching polytopal meshes with geometric interface inconsistencies

A new variational mimetic finite difference method for elliptic interface problems with perfect and imperfect thermal contacts on non-matching polytopal meshes with geometric interface inconsistencies is developed and analyzed theoretically and numerically. The method is defined on multiple non-matching submeshes with gaps and overlaps along their interfaces. The discrete equations are derived from a minimization problem for the augmented Dirichlet functional. For a perfect thermal contact, the functional uses a modified mimetic gradient with extended stencil which couples unknowns from both sides of an interface, as well as penalty terms to enforce weak continuity of temperature across the interface. The method leads to a symmetric positive definite matrix for any scaling of the penalty terms. For an imperfect thermal contact, the Dirichlet functional is supplemented with a quadratic jump term along the interface related to the interface thermal resistance. We prove that the method conserves the total heat flux across each interface. In conclusion, the obtained results are verified with numerical experiments showing convergence in the discrete L 2 and L ∞ norms.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

A careful examination of closure models in Euler–Lagrange Simulations of compressible multiphase flow in a planar shock particle curtain problem

In this work we present a comprehensive investigation of state-of-the-art closure models employed to represent interphase momentum, thermal, and work exchange between the gas and particulate phases for Euler–Lagrange (EL) simulations in shock-driven flows. A complete list of closures for the force, torque, heat transfer, and work exchange models is provided. In particular, the present work includes a stochastic closure for the particle-to-particle variation in the quasi-steady force and a deterministic closure for particle-to-particle variation in the added mass force in an EL framework. These variations arise due to the presence of neighboring particles and particle–particle interactions. To investigate the importance of each closure term, we carry out fully three-dimensional simulations for a planar shock propagating over a random bed of inert particles. The primary goal is to evaluate the role of each closure term on the gas dynamic features (such as transmitted and reflected shock locations) and particle curtain features (such as upstream and downstream curtain locations). To this end, thirteen cases are considered, with each case progressively including a closure model with the goal to identify and quantify its contribution to the simulated dynamics. We show that the volume fraction dependence of the mean force models plays an important role in generating wave-like instabilities that lead to concentration bands. In addition, fluctuations in quasi-steady and added mass forces primarily decrease the internal instabilities that tend to enhance local volume fraction variations. Particle rotation is primarily due to inter-particle collisions, is generally weak, and does not play an important role in the translational dynamics for the present configuration. Inter-phase heat transfer has a strong effect on gas phase temperature, slows down the transmitted and reflected shocks, and decreases the width of the curtain. Furthermore, the absence of a work-coupling model fails to conserve the total energy, greatly under-predicts the gas temperature which in turn affects the particle dynamics.

Compressible flow↗
DOE OSTI2024

Spectral induced polarization (SIP) measurements across a PFAS-contaminated source zone

There is a pressing need for the development of field-scale, in situ screening technologies for assessing variations in aqueous film forming foam (AFFF) concentrations in soils at former fire training and storage sites. Field-scale Spectral Induced Polarization (SIP) geophysical measurements were acquired on a transect crossing an AFFF source zone. Soil samples were acquired from ten locations and used to determine variations in poly- and per-fluoroalkyl substances (PFAS) concentrations in soils and soil texture. These samples were also used to create triplicate soil columns for laboratory-grade SIP measurements. Field and laboratory observations provide evidence that SIP measurements are sensitive to the concentration of AFFF constituents associated with the pore surface in soils. The phase of the SIP measurements on the laboratory samples was linearly correlated with total soil-sorbed PFAS concentration. The phase from the field SIP measurements was highest over the location of maximum PFAS concentration measured on the laboratory samples, although a significant correlation between field-measured phase and laboratory-measured total PFAS concentration was not established. The sensitivity of the SIP response to the removal of soil PFAS using a methanol wash procedure (total PFAS concentration drop of 366 ppb) adds evidence for the case for SIP characterization of AFFF source zones. The results of these studies suggest that SIP might be developed into a field-scale technology for rapid, indirect assessment of AFFF source zones. Such a technology could improve the effectiveness of AFFF source zone characterization at reduced costs.

58 GEOSCIENCES↗
DOE OSTI2023

Comparison of the Chemical Composition of Liquids from the Pyrolysis and Hydrothermal Liquefaction of Lignocellulosic Materials

Major differences in thermal stability and hydrotreatment behavior of HTL and pyrolysis oils have been reported in the literature. However, little is known about the variations in the chemical composition of these oils that could explain such differences. Two commercial wood pyrolysis oils (Pyrovac and BTG), and their water-soluble (WS) and water-insoluble fractions (WIS) were analyzed and compared with the aqueous (WS WD-57 ) and oily (WIS WS-57 ) fractions obtained from hydrothermal liquefaction (HTL) of Douglas-fir. The samples were characterized by GC/MS, Karl Fischer titration, carbonyl content, total acid number, elemental composition, calorific value, proximate analysis, Fourier Transform Infrared Spectroscopy (FTIR), Folin-Ciocalteu (FC), and UV fluorescence. All the fractions were also analyzed by Fourier Transform Ion Cyclotron Resonance Mass Spectroscopy (FT-ICR-MS) and by Electrospray Ionization (EI). The most prevalent class of compounds in the water insoluble phases were phenols derived from lignin. Water-soluble phases contain mostly the oxygenated compounds derived from cellulose and hemicellulose and were richer in carbonyl functional groups. The water content of the resulting aqueous phases were between: 65 (WS BTG ) and 96 (WS Pyrovac ) wt. %. The bio-oil from BTG has higher water content and lower HHV, compared to Pyrovac oil. The GC/MS results of BTG oil show the presence of a more prominent acetic acid peak and higher TAN number than the Pyrovac oil. The GC/MS of Pyrovac oil showed more obvious mono-phenol peaks. The quantification of this family by Folin-Ciocalteu method confirmed higher content of monophenolic compounds compared with the BTG oil. The lower thermal stability of pyrolysis oils compared with HTL biocrudes can be partially explained by the fact that pyrolysis oils (BTG and Pyrovac) contain carbohydrates while HTL biocrude (WIS WD-57 ) doesn’t. Thus, we decided to further investigate the chemical differences between the phenolic rich fractions insoluble in water and the holocellulose derived compounds soluble in water. Even after water extraction, the acid content of the water insoluble fraction from BTG (WISBTG) was higher than the acid content of the water insoluble fraction obtained by HTL (WISWSD-57). Likewise, the acid content of the aqueous phases derived from pyrolysis oils (WS Pyrovac , WS BTG ) was also higher than for the aqueous phase obtained by HTL (WSWSD-57). This result is in part due to the use of bases in the HTL process that neutralizes the acid formed in that process. Moreover, the starting feedstock may also influence the differences between the oils. Although, the UV-Fluorescence spectra, ICR-MS and the EI analyses showed some minor differences in the molecular weight and chemical make-up of the oligomers soluble and insoluble in water from pyrolysis and HTL; the differences observed were not large enough to justify the differences in behavior between these oils reported in the literature. Our results suggest that the differences observed between HTL biocrudes and pyrolysis oils are likely partially due to the presence of holocellulose derived products in the pyrolysis oils and higher acid contents.

09 BIOMASS FUELS↗
DOE OSTI2023

TRACER-Vertical Profile of Aerosol in Convective System (TRACER-VPACS) Field Campaign Report

We investigated physicochemical properties of aerosol at different altitudes during the U.S. Department of Energy (DOE)’s Atmospheric Radiation Measurement (ARM) user facility Tracking Aerosol Convection Interactions Experiment (TRACER). We used ARM’s tethered balloon system (TBS) and multi-modal offline analysis techniques available at DOE’s Environmental Molecular Sciences Laboratory (EMSL) at Pacific Northwest National Laboratory (PNNL) and at the Advanced Light Source (ALS) at the Lawrence Berkley National Laboratory (LBNL). The main scientific aim of this project was to understand the influence of the vertical profile of aerosol composition on the properties of deep convective clouds. The effect of aerosol properties on convective activity is still poorly understood and debated. For example, previous studies suggested that higher concentration of hygroscopic aerosol may lead to lower collection efficiency and eventually delay in onset of precipitation and total amount of precipitation. However, some observations suggested that increased aerosol concentration may result in vigorous convection and a higher amount of precipitation. In addition, the radiative forcing of the ambient aerosol is strongly influenced by the vertical distribution of aerosols and their chemical composition, and not well constrained in climate models. The estimation of the indirect effect of aerosol also has a significant uncertainty due to the spatial variation of the aerosol, especially at the cloud base. Detailed size-resolved chemical composition studies of aerosol and their microphysical properties, along with in situ measurements from ARM’s mobile facilities, help to improve our understanding of the role of aerosol in deep convective cloud formation over a broad range of environmental thermodynamic and aerosol conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

ELG spectroscopic systematics analysis of the DESI Data Release 1

Dark Energy Spectroscopic Instrument (DESI) uses more than 2.4 million Emission Line Galaxies (ELGs) for 3D large-scale structure (LSS) analyses in its Data Release 1 (DR1). Such large statistics enable thorough research on systematic uncertainties. In this study, we focus on spectroscopic systematics of ELGs. The redshift success rate (f goodz ) is the relative fraction of secure redshifts among all measurements. It depends on observing conditions, thus introduces non-cosmological variations to the LSS. We, therefore, develop the redshift failure weight (w zfail ) and a per-fibre correction (η zfail ) to mitigate these dependences. They have minor influences on the galaxy clustering. For ELGs with a secure redshift, there are two subtypes of systematics: 1) catastrophics (large) that only occur in a few samples; 2) redshift uncertainty (small) that exists for all samples. The catastrophics represent 0.26% of the total DR1 ELGs, composed of the confusion between [O ii] and sky residuals, double objects, total catastrophics and others. We simulate the realistic 0.26% catastrophics of DR1 ELGs, the hypothetical 1% catastrophics, and the truncation of the contaminated 1.31 < z < 1.33 in the A BACUS S UMMIT ELG mocks. Their P ℓ show non-negligible bias from the uncontaminated mocks. But their influences on the redshift space distortions (RSD) parameters are smaller than 0.2σ. The redshift uncertainty of DR1 ELGs is 8.5km s -1 with a Lorentzian profile. The code for implementing the catastrophics and redshift uncertainty on mocks can be found in https://github.com/Jiaxi-Yu/modelling_spectro_sys.

79 ASTRONOMY AND ASTROPHYSICS↗