Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “trending”

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 865 records · Page 48
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Microbial community diversity changes during voltage reversal repair in a 12-unit microbial fuel cell

Microbial fuel cell stacks (MFC-Stack) are often confronted with voltage reversals, likely due to an interplay between microbial community dynamics and insufficient electric circuit balancing. Herein, we provide new insight into voltage reversals by examining the microbiomes of twelve MFC units of a 12-liter Pilot-MFC-Stack during repair. Different biofilm repair methods (self-healing, electrostimulation, and re-acclimatization upon cross-inoculation) were used to evaluate the microbial community response. In addition, MFC-Stack simulation was performed based on Kirchhoff’s Second Law to predict values for source potentials and post-evaluate internal resistances. Analysis of the 16S rRNA amplicon sequencing data suggests that the biofilm repair methods could slowly heal damaged biofilms. Notably, severely voltage reversed MFC units had low electrogen relative abundances (18%) and positive anode potentials, while strong bioanodes and contained more than 50% electrogens and had negative anode potentials. Between-community analyses (beta diversity ordination and multinomial regression) of the voltage reversed MFC units revealed differences among biofilms in contrast to healthy/strong MFC units. Permutational multivariate analysis of variance (PERMANOVA) confirmed that reversed biofilms were, indeed, significantly (p < 0.05) different from stronger ones. Overall, these analyses demonstrated the utility of combining electrotechnical and microbial community analyses, especially beta diversity ordination and multinomial regression, to understand problematic MFC units and the potential success of a biofilm repair method. Finally, thicker biofilms were usually healthier and stronger, although thickness was no guarantee for proper structure and power function as all factors were interdependent. There was an evolutionary trend that strong anodes became stronger/healthier and others weaker. This spontaneous trend has to be considered to avoid irreversible voltage reversals and to repair electrogenic biofilms in an MFC-Stack.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2026

Harnessing citizen science to contextualize adaptation mechanism discovery

Species occupying broad geographic regions have evolved multiple mechanisms to regulate phenological characteristics, enabling adaptations to diverse native habitats. By developing computer vision AI to process citizen science observations across native habitats over North America, we uncovered a consistent latitudinal trend of earlier flowering at higher latitudes in warm-season perennial grasses. To explore the underlying mechanisms of adaptation, we conducted common garden experiments with one species (switchgrass) and discovered the opposite latitudinal flowering-time trend. Integration of differential plasticity of GI-Hd1-FTL1 haplotypes of flowering time regulatory genes, haplotype range, and local environmental profiles found that observations from native habitats capture only part of the genotype-environment-phenotype spectrum established in common garden experiments, therefore reconciling the discrepancy. Two mechanisms emerged as key forces shaping current haplotype ranges and influencing future shifts. Our study highlights the power of combining citizen science observations with designed experiments to uncover mechanisms of adaptation across spatiotemporal scales.

FTL1↗
DOE OSTI2022

Chemical and isotopic evolution of flowback fluids from the Utica Gas Shale Play, Eastern Ohio USA

For this work, hydraulic fracturing flowback fluids were collected from two Utica/Point Pleasant well pads in eastern Ohio. One site was in the wet gas zone (UPPW), while the site ~50 km south consisted of four wells on the same pad in the dry gas zone (UPPS). Samples of input fluids also were collected before and during hydraulic fracturing. Flowback fluids are Na-Ca-Cl brines with total dissolved salt (TDS) concentrations that increase over several months from ~100 to 200 g/L. The slightly higher TDS of the dry gas fluids are in part due to recycled flowback used as input fluids, and in part due to lower volume of water used in hydraulic fracturing. Concentrations of most major ions (Ca 2+ , Mg 2+ , Na + , Sr2 + , Fe, Mn, Cl – , Br – ) are similar for the five wells sampled, although small but systematic changes occur in the major element ratios over time. Most notably, an increase in the Sr/Cl ratio corresponds to an increase in the 87 Sr/ 86 Sr ratio in the fluids, suggesting interactions with a more radiogenic Sr source in the subsurface. Dissolved Ba concentrations and Ra activities were different between the two sites, reflecting a high SO 4 2– fluid used at the UPPW site, and water-rock reactions occurring during hydraulic fracturing at the UPPS site. Water oxygen (δ 18 O) and hydrogen (δD) isotopes for input fluids used in the UPPW4 well and fresh water used for the UPPS wells fall on the Global Meteoric waterline (GMWL). Flowback fluids from both sites are relatively enriched in δ 18 O and δD compared to the input, but do not appear to follow a simple mixing trend, suggesting reaction and isotopic exchange with carbonates and fractionation due to imbibition in the rock. Cl isotopes, δ 37 Cl in the FP fluids varied from ~ –0.43 to +0.13‰, the largest variation was observed in the earlier stages of flowback, while in the later stages δ 37 Cl exhibited a small but systematic increase over time, suggesting diffusion control of isotopic composition. Some of the trace species measured (dissolved Fe, Mn, thiosulfate and organic acids) do not follow the same trends as the major ions, suggesting contributions from input fluids, or microbially mediated reactions are exerting control on their concentrations. Input water chemistry exerts an important control on the Sr and Ba concentrations in flowback water. High SO 4 2– in the input fluids used for hydraulically fracturing the wet gas well leads to precipitation of barite-celestite in flowback fluids.

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

Development of fiber orientation in injection molding: Comparison of glass fiber, carbon fiber and their hybrid composites

Fiber orientation distribution (FOD) in injection-molded panels with respect to distance from the gate was analyzed using X-ray computed tomography (X-CT) for glass fiber (GF), carbon fiber (CF), and hybrid CF/GF (CGF) reinforced nylon 66. To understand the reason behind the FOD with different fiber types, computational fluid dynamics (CFD) and rheology were performed. Samples were extracted at three locations: near the gate, center, and opposite end. Thickness of the layers of typical skin-shell-core type FOD varied with fiber type and location. GF achieved flow direction alignment (in shell) earlier than viscous CF and CGF near the gate, whereas CF showed the highest flow-direction alignment at the center due to shear induced orientation. At the opposite end, GF experienced more backflow than others indicating faster mold filling owing to its lower viscosity. Hybrid CGF exhibited GF-dominated center and CF-dominated end region. The numerical model used to obtain FOD and rheological predictions for the CF and GF composites served to corroborate the trends observed in the experimental trials. The FOD responses across fiber types and location were reflected in their longitudinal and transverse properties. Only GF showed higher longitudinal modulus over transverse modulus near the gate attributed to rapid alignment, whereas CF and CGF exhibited opposite trend. However, fountain flow enhanced the longitudinal modulus over transverse modulus with the distance for all, particularly for CF. This study offers insights into mold filling behavior of different fibers which are critical in optimizing injection molding conditions for tailored final properties.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Fuel molecular structure effect on autoignition of highly branched iso-alkanes at low-to-intermediate temperatures: Iso-octane versus iso-dodecane

Highly branched iso-alkanes are an important class of hydrocarbons found in conventional petroleum-derived and alternative renewable fuels used for combustion applications. Recognizing that chemical kinetics for most of these iso-alkanes, especially at low-to-intermediate temperatures, has not been well studied, an experimental and modeling investigation of two selected iso-alkanes, iso-octane (2,2,4-trimethylpentane, iC8) and iso-dodecane (2,2,4,6,6-pentamethylheptane, iC12), is conducted to understand the fuel molecular structure effect on their autoignition characteristics. Using a rapid compression machine (RCM), the ignition responses of iC8 and iC12 at varying pressures, temperatures, and equivalence ratios are characterized and compared. The newly-acquired experimental ignition delay times have been compared with the literature RCM and shock tube data, demonstrating the complementary nature of the current dataset. Further comparison of the experimental pressure traces and ignition delay times illustrates the reactivity crossover between iC8 and iC12. Namely, there exists a temperature window in the negative temperature coefficient regime within which iC12 is less reactive than iC8, but iC12 becomes more reactive outside this temperature window. Furthermore, a chemical kinetic model of iso-alkanes including both iC8 and iC12 is developed. Simulated results using this model are then compared to the experimental data obtained in this study and available in the literature, showing its ability to predict the experimental trends. Chemical kinetic analyses have also been conducted to identify the important reaction pathways controlling autoignition at varying conditions, and to elucidate the underlying mechanism leading to different reactivity trends between iC8 and iC12.

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

An experimental and kinetic modeling study of NO x sensitization on methane autoignition and oxidation

An experimental and kinetic modeling study of the influence of NO x (i.e. NO 2 , NO and N 2 O) addition on the ignition behavior of methane/‘air’ mixtures is performed. Ignition delay time measurements are taken in a rapid compression machine (RCM) and in a shock tube (ST) at temperatures and pressures ranging from 900–1500 K and 1.5–3.0 MPa, respectively for equivalence ratios of 0.5–2.0 in ‘air’. The conditions chosen are relevant to spark ignition and homogeneous charge compression ignition engine operating conditions where exhaust gas recirculation can potentially add NO x to the premixed charge. The RCM measurements show that the addition of 200 ppm NO 2 to the stoichiometric CH 4 /oxidizer mixture results in a factor of three increase in reactivity compared to the baseline case without NOx for temperatures in the range 600–1000 K. However, adding up to 1000 ppm N 2 O does not show any appreciable effect on the measurements. The promoting effect of NO 2 was found to increase with temperature in the range 950–1150 K, while the sensitization effect decreases at higher pressures. The experimental results measured are simulated using NUIGMech1.2 comprising an updated NO x sub-chemistry in this work. A kinetic analysis indicates that the competition between the reactions ĊH 3 + NO 2 ↔ CH 3 Ȯ + NO and ĊH 3 + NO 2 (+M) ↔ CH 3 NO 2 (+M), the former being a propagation reaction and the latter being a termination reaction governs NO x sensitization on CH 4 ignition. Recent calculations by Matsugi and Shiina (A. Matsugi, H. Shiina, J. Phys. Chem. A. 121 (2017) 4218–4224) for the nitromethane formation reaction CH 3 + NO 2 (+M) ↔ CH 3 NO 2 (+M), together with the recently calculated rate constants for HONO/HNO 2 reactions significantly improve ignition delay time predictions in the temperature range 600–1000 K. Furthermore, the experiments with NO addition reveal a non-monotonous sensitization impact on CH 4 ignition at lower temperatures with NO initially acting as an inhibitor at low NO concentrations and then as a promoter as NO concentrations increase in the mixture. This non-monotonous trend is attributed to the role of the chain-termination reaction ĊH 3 + NO 2 (+M) ↔ CH 3 NO 2 (+M) and the impact of NO on the transition to the chain-branching steps CH 2 O + HȮ 2 ↔ HĊO + H 2 O 2 , H 2 O 2 (+M) ↔ ȮH + ȮH (+M), HĊO ↔ CO + Ḣ followed by CO + O 2 ↔ CO 2 + Ö and Ḣ + O 2 ↔ Ö + ȮH. NUIGMech1.2 is systematically validated against the new ignition delay measurements taken here together with species measurements and high temperature ignition delay time data available in the literature for CH 4 /oxidizer mixtures diluted with NO 2 /N 2 O/NO and is observed to accurately capture the sensitization trends.

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

Isomeric effects on the reactivity of branched alkenes: An experimental and kinetic modeling study of methylbutenes

Here, a detailed experimental study of the low-to-intermediate temperature combustion of methylbutene isomers, i.e., branched C 5 alkenes, has been undertaken with multiple experimental facilities. Ignition delay times were measured at equivalence ratios 0.5–2.0, 685–1020 K and up to 45 bar condition from two rapid compression machines and showed slight deviation from an Arrhenius behavior for all three isomers, while their reactivity order differs as temperature changes. Sampled intermediates formed during the oxidation process of mixtures at 900–1150 K and 0.82 bar from a flow reactor and at 730 K and 20 bar from a rapid compression machine were analyzed using gas chromatography techniques. Trends in the formation and consumption of sampled intermediates were modeled using a kinetic model developed in this work for all three isomers. Rate of production and sensitivity analyses emphasize the role of double bond-specific reactions governing the global reactivity of these fuels. Additional studies of the addition reactions of HO 2 radicals to the double bond and to allylic radicals may improve the model performance.

2-Methyl-1-butene↗
DOE OSTI2026

Carbon fiber design improvements based on economic models

Circular fiber geometries are predominant in commercial carbon fiber material systems, but the use of this fiber shape has numerous limitations. Circular geometries have continuous symmetry, which is helpful for various processing considerations, but also have the largest possible maximum diffusion thickness for a given fiber area. This characteristic means that circular carbon fibers always have the highest material processing cost and lowest production throughput compared to any other fiber shape with the same area and tow count. To quantify material cost and other benefits for non-circular carbon fiber geometries, process models for polyacrylonitrile based carbon fiber production are developed in relationship to the carbon fiber shape, size, and tow count. For a given fiber shape, precursor production costs are shown to favor maximizing fiber size while conversion costs are minimized by the smallest fiber size. These competing cost trends result in a numerically optimal fiber size for a given shape and tow count, while both cost components are decreased by increasing tow count. Cost–performance tradeoffs for three lobe fiber geometries are studied by supplementing the cost trends with a numerical failure model to predict compressive strength for discrete shape variants. The shape selection is shown to be more sensitive to variations in cost than compressive strength while suboptimal shape designs can improve manufacturing robustness and achievable fiber volume fractions. Finally, an optimal three lobe carbon fiber is identified that balances the set of considerations while reducing costs and embodied energy and increasing production throughput compared to a commercial carbon fiber.

Carbon fiber↗
DOE OSTI2026

Assessment of the hydromechanical higher-order MPM for the simulation of geotechnical problems

The Material Point Method (MPM) has been increasingly used to simulate large strain deformations. Linear interpolation functions are commonly used to perform the spatial integration. It is well-known that the discontinuities in the interpolation function derivatives induce shock-like artifacts known as ‘cell-crossing’ error. These errors compound with volumetric locking errors when used with hydromechanical formulations for porous media, where different velocity fields are used for each phase. The capabilities of higher-order MPM frameworks have not been explored for real-scale geotechnical problems. As such, this paper aims to assess, validate, and further discuss a higher-order B-spline MPM (BS-MPM) framework. First, the BS-MPM framework is verified against the large-strain oedometer consolidation problem. Second, the framework is validated against a real-scale slope failure experiment triggered by pore water pressure recharge. Landslide features that are captured using the higher-order framework are specifically highlighted, and results (e.g., pore water pressure and deformation) are validated with field measurements. A generally convergent numerical solution is observed when using cubic interpolation functions. Third, a footing penetration problem is simulated using the multi-patch BS-MPM. Trends are examined with respect to penetration velocity and variation in hydraulic conductivity. The BS-MPM framework ultimately presents a stabilized numerical solution that captures plausible hydromechanical interaction trends important in geotechnical engineering applications.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Advancing equitable value chains for the global hydrogen economy

Hydrogen is a rapidly growing focus for countries seeking to develop green industries, but there are many questions about how the nascent global hydrogen economy will develop, and what this implies for equitable sharing of benefits and burdens between nations. In this perspective we summarize emerging trends in national hydrogen strategies and develop recommendations for researchers and policymakers to center equity in hydrogen development. This will require integrating innovation and development perspectives on international technology transfer, developing more detailed representation of hydrogen trade in systems models, building equity considerations into national and international planning processes, and establishing robust technology transfer efforts. In conclusion, policymakers will also need to grapple with the difficulties of verifying life cycle emissions of hydrogen if hydrogen trade emerges as a significant trend, potentially requiring new methods of emissions accounting and trade reforms that prioritize international equity.

08 HYDROGEN↗
DOE OSTI2021

Miscellaneous Electric Loads: Characterization and Energy Savings Potential

Over time, miscellaneous electric loads (MELs) are expected to increase both in magnitude and share of residential and commercial building energy consumption. This trend is most apparent in North America, but it is also occurring in Japan and Europe. However, the contribution of MELs to building energy use is not currently well understood, both because the products in this category are transforming rapidly and the definition and classification of MELs is ambiguous. This study estimated the national energy consumption of 36 MELs using best-available data and found them to comprise 12% of delivered electricity to the U.S. residential and commercial building sectors. If 26 of these MELs were replaced with the most energy-efficient product models available on the market, their energy consumption could be halved to 6% of delivered electricity. National energy models will better account for building energy consumption by incorporating the MELs data collected and analyzed for this study, leading to improved policy decisions.

Miscellaneous electric loads, Plug loads, Taxonomy↗
DOE OSTI2025

The short-term comprehensive impact of the phase-out of global coal combustion on air pollution and climate change

With the continuous intensification of global warming, the reduction and ultimate phase-out of coal combustion is an inevitable trend in the future global energy transformation. This study comprehensively analyzed the impact of phasing out coal combustion on global emissions and concentrations of air pollutants, radiative fluxes, meteorology and climate using Community Earth System Model 2 (CESM2). The results indicate that after the global phase-out of coal combustion, there is a marked decrease in the concentrations of sulfur dioxide (SO 2 ), nitrogen oxides (NO x ) and fine particulate matter (PM 2.5 ), with some regions experiencing a reduction of exceeding 50%. There is no significant change in global ozone (O 3 ) concentration. There are decreasing AOD and positive radiative fluxes globally in the short term, though the cloud contributes minor negative radiative fluxes. The global air temperature may increase by approximately (0.02 ± 0.15) °C on average with regional and seasonal variations, and the precipitation may potentially increase by approximately (2.7 ± 40.6) mm yr −1 globally and over 20% in equatorial regions in the short term. But combined with the decreasing trend of cloud water content in the Northern Hemisphere, it indicates a potential increase in the extremity of precipitation events. In conclusion, this study provides references for global control of air pollution, mitigation strategies of climate change, and transformation of energy structures under the objective of “carbon neutrality”, such as focusing on the negative climate impacts of exacerbating regional warming and increasing extreme precipitation resulting from the rapid reduction of aerosols in the short term.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2023

The effect of terrain-influenced winds on fire spread in QUIC-Fire

Here in this manuscript, we describe the implementation of the terrain-following version of QUIC-URB into QUIC-Fire and a demonstration of the impacts of terrain-influenced winds on QUIC-Fire-simulated fire spread. No changes to the underlying QUIC-Fire fire spread algorithm were made other than what was required to correctly account for the inclusion of terrain. This paper summarizes simulations used to understand how the QUIC-URB terrain-influenced winds affect upslope fire behavior without additional changes to the QUIC-Fire fire spread algorithm. Previously published FIRETEC results are compared to simulation results from the modified QUIC-Fire (incorporating terrain-influenced winds) that use the same topographies and fuels. QUIC-Fire results showed overall similar behaviors in terms of how the topographies affected fire shapes and trends in spread rates. Due to the terrain-following version of QUIC-URB being unable to generate flow separations at the crest of hills, fire spread rates in these regions across all non-flat topographies were over-predicted when compared to FIRETEC. Lateral fire growth showed similar trends with FIRETEC between topographies but did not capture the increase in spread due to a diagonal interface between grassland and forested fuel regions in the test domain. These simulations suggest possible refinements that are necessary to improve QUIC-Fire and thus guide ongoing efforts related to: how flame tilt angle is accounted for, the incorporation of non-local drag effects, and the inclusion of the wake-eddy parameterizations that are used in QUIC-URB.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Micromovements and discomfort associated with flight mission with helmet operation tasks with different levels of cognitive workload

When performing stationary tasks under elevated cognitive workload, individuals must perform continual muscle contractions to maintain stability of the body, resulting in fatigue of the postural muscles. When the muscles perform these contractions in a prolonged manner, the body potentially responds through small changes in body movements—micromovements that may lead to discomfort. The study purpose was to evaluate impact of cognitive load on micromovements. The micromovements were measured during three different cognitive workloads; low, medium, and high. The NASA-TLX score was used to evaluate the perceived mental workload and discomfort was assessed by visual analog scale. In total, 60 subjects (30 males and 30 females) were recruited and performed cognitive tasks that simulated flight operations such as changing the radio frequency based on air traffic control messages, balancing the fuel levels in simulated fuel tanks, and aiming a reticle in a designated moving target using the cyclic control. Cognitive load was defined by the frequency of events. Micromovements were defined by changes in the center of pressure (COP) of the seat pan and COP standard deviation. It was found that the high cognitive workloads had the highest NASA-TLX scores including mental demands, temporal demands, and effort. The neck area had the highest overall levels of discomfort followed by upper back. The highest standard deviation for COP shift and number of micromovements occurred for medium cognitive workloads. In conclusion, while there were some interesting trends, few trends reached a statistical significance due to high variability among subjects for the outcome variables.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Degradation and Modeling of Large-Format Commercial Lithium-Ion Cells as a Function of Chemistry, Design, and Aging Conditions

Demand for large-format (>10 Ah) lithium-ion batteries has increased substantially in recent years, due to the growth of both electric vehicle and stationary energy storage markets. The economics of these applications is sensitive to the lifetime of the batteries, and end-of-life can either be due to energy or power limitations. Despite this, there is little information from cell manufacturers on the sensitivity of cell degradation to environmental conditions or battery use. This work reports accelerated aging test data from four commercial large-format lithium-ion batteries from three manufacturers, with varying design (thickness, casings, ...), chemistry (lithium-iron-phosphate (LFP) or lithium-nickel-manganese-cobalt-oxide positive electrodes (NMC), with graphite (Gr) negative electrodes), and capacity (50 to 250 Amp hours). The tested LFP|Gr cell is found to be relatively insensitive to cycling conditions like temperature or voltage window, while NMC|Gr cells have varying sensitivity. Degradation trends are further investigated by training predictive models: simple polynomial trend lines, a semi-empirical reduced-order model, and an empirical reduced-order model identified using machine-learning based on symbolic regression. Calendar and cycle life are simulated over a variety of conditions to directly compare the various batteries. Cell size and thickness are found to substantially impact sensitivity to temperature during cycle aging, while electrode chemistry impacts depth-of-discharge sensitivity. Real-world battery lifetime is evaluated by simulating residential energy storage and commercial frequency containment reserve systems in several U.S. climate regions. Predicted lifetime across cell types varies from 7 years to 20+ years, though all cells are predicted to have at least 10 year life in certain conditions.

battery lifetime↗
DOE OSTI2023

Range-wide population assessments for subalpine fir indicate widespread disturbance-driven decline

Subalpine forests in western North America are threatened by rapid climate change, increased activity by endemic and exotic insects and diseases, and changing wildfire regimes. The interactive effects of these stressors have resulted in pronounced population declines in many subalpine tree species; however, a systematic assessment of the status and trends of subalpine forests is lacking. Subalpine fir (Abies lasiocarpa) is a widespread species across the western United States, with documented population declines in many parts of its distribution. Here we use subalpine fir as an initial leverage point to build a more complete understanding of subalpine forest baseline conditions and responses to environmental change. Specifically, we leverage the USDA Forest Service Forest Inventory and Analysis (FIA) database to (1) ask how subalpine fir populations are changing across the species’ distribution in the western US, (2) assess the drivers of recent subalpine fir population trends, and (3) explore whether those changes imply generalized species-wide and/or system-wide decline. We found that subalpine fir abundance and basal area are declining concurrently across ~ 62% of the species’ distribution, and increasing across ~ 19%. Range-wide, we estimated 25.02 ± 2.74 % subalpine fir mortality between 2000 and 2009 and 2010–2019 FIA inventory periods, with higher mortality concentrated in the eastern Oregon Cascades, central Idaho, and parts of southern Colorado. High regeneration density did not predict positive population trajectories, which were instead associated with higher rates of adult recruitment. While the importance of different mortality agents varied substantially between ecoregions, 83.4% of total range-wide mortality was related to fire or biological disturbance. Declining subalpine fir basal area coincided with declines in the basal area of other co-occurring tree species in 39% of subalpine forest area, and with increases in conspecific basal area in 22% of forest area. Fire disturbance was the single largest cause of subalpine fir mortality; however, even where subalpine fir fire mortality was high, mortality among other species was primarily caused by insects. In conclusion, our results suggest that subalpine fir declines across large portions of the western United States are driven by forest disturbance, and that declines in subalpine fir populations may be indicative of negative change in subalpine forest systems broadly.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Numerical analysis of soot emissions from gasoline-ethanol and gasoline-butanol blends under gasoline compression ignition conditions

In the present work, computational fluid dynamics (CFD) simulations of a single-cylinder gasoline compression ignition (GCI) engine were performed to investigate the impact of blending two biofuels, ethanol and n-butanol, with gasoline on the trade-off between combustion phasing and soot emissions under low load conditions. Here, in order to represent market gasoline (RD5-87), a four-component toluene primary reference fuel (TPRF) + ethanol (ETPRF) surrogate (with 20% ethanol by mole; E20) was formulated using a neural network based octane predictor such that the surrogate had the same ethanol content, Research Octane Number (RON) and Octane Sensitivity (S). In addition, a novel skeletal kinetic mechanism for ETPRF and TPRF + n-butanol (BTPRF) blends, incorporating polycyclic aromatic hydrocarbon (PAH) chemistry, was developed. A three-dimensional (3D) engine CFD formulation employing the skeletal mechanism, adaptive mesh refinement (AMR), finite-rate chemistry approach, and hybrid method of moments (HMOM) was adopted to capture the in-cylinder combustion phenomena and soot emissions. The engine CFD model was validated against RD5-87 experimental data for a broad range of start-of-injection (SOI) timings (-21/-27/-36/-45 crank angle degrees (CAD) after top-dead center (aTDC)), with respect to in-cylinder pressure, heat release rate, combustion phasing, and soot emissions. The closed-cycle simulation results were analyzed to elucidate the non-monotonic trend of soot emissions versus SOI timing: SOI-36 > SOI-45 > SOI-21 > SOI-27. Thereafter, the validated CFD model was employed to simulate the combustion of a gasoline-ethanol blend with 45% (by mole) ethanol (E45) and a gasoline-butanol blend with 45% (by mole) n-butanol (B45) under the same operating conditions to study the effects of fuel composition and SOI timing on combustion phasing and soot emissions. The sooting propensity followed the trend: B45 > E20 > E45 at all SOI timings. Overall, it was observed that the autoignition propensity was primarily related to fuel chemistry. On the other hand, sooting propensity showed strong coupling with both fuel chemistry and physical properties, with greater impact of fuel physical properties at advanced SOI timings.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2022

Experimental study of the partitioning of some platinum group elements (Pd and Ir) between orthopyroxene and silicate melt

Past experiments and observations on natural samples have largely focused on the roles of olivine and chromite in controlling the behaviour of the platinum-group elements (PGE) during melting and solidification, whereas other phases, such as pyroxene, have gone largely uncharacterized. Here, to address this, experiments have been done to measure the partitioning of Pd (with a subset of results for Ir), between orthopyroxene and silicate melt at 1340 °C, 0.1 MPa and log fO 2 of FMQ - 1 to FMQ + 6 (FMQ = Fayalite-Magnetite-Quartz). The X-ray Absorption Near-Edge Structure (XANES) was measured in a subset of experiment glasses. Glass concentrations of Pd (corrected to unit Pd activity) increase from ~6 to ~650 ug/g over the fO 2 range of experiments. The slope of the solubility-fO 2 relation is consistent with Pd 1+ as the dominant oxidation state, with evidence for Pd 0 and Pd 2+ at the lowest and highest experiment fO 2 , respectively. Consistent with this result, the XANES reveal spectral features similar to Pd 0 and Pd 2+ spectral reference materials (specRM) at the most reduced and oxidized synthesis conditions, respectively. Other lines of evidence require the presence of a third melt species, here interpreted to be Pd 1+ . Values of orthopyroxene/melt partition coefficients for Pd (D Pd Opx/melt ) are 0.0051 (+/-0.006) at log fO 2 < ΔFMQ + 3, increasing with fO 2 to a maximum of 0.013 at ~FMQ + 6. Sodium partition coefficients, expected to be similar to Pd, range from 0.0061 (+/-0.00061) at FMQ + 3, increase to 0.007–0.009 at higher fO 2 , but with no clear systematic trend. A value for D Ir Opx-melt of ~0.6 was measured at ~FMQ + 4, indicating significantly more compatible behaviour for Ir relative to Pd. Partitioning results are interpreted in the context of the Blundy-Wood elastic strain model in which the variation in partitioning is related to ionic radius mismatch to an optimal crystallographic site size. Based on the trend in ionic radius with oxidation state, the estimated ionic radius of Pd 1+ in octahedral coordination is similar to Na 1+ , and comparison to previous orthopyroxene-melt partitioning experiments suggests D Pd1+ opx/melt and DNaopx/melt should be nearly identical, consistent with the results of this study. The ionic radius of VI-fold Pd 2+ is close to the optimal M2 site size, so an increased proportion of this species with fO 2 accounts for the larger values of D Pd opx/melt at the highest fO 2 investigated. The much larger partition coefficient for Ir is consistent with the presence of Ir 2+ , whose estimated ionic radius is close to Fe 2+ and Mg 2+ , as well as predictions for the optimal M1 site size. With the assumption that D Pd opx/melt = D Na opx/melt , combined with a revised value for the Pd content of the primitive mantle, a melting model is presented that better reproduces the Pd concentration of high degree melts from sulfide-free mantle sources.

58 GEOSCIENCES↗