Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Effective diffusivity”

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,261 records · Page 70
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Cross-Code verification and sensitivity analysis to effectively model the electrothermal instability

This manuscript presents verification cases that are developed to study the electrothermal instability (ETI). Specific verification cases are included to ensure that the unit physics components necessary to model the ETI are accurate, providing a path for fluid-based codes to effectively simulate ETI in the linear and nonlinear growth regimes. Two software frameworks with different algorithmic approaches are compared for accuracy in their ability to simulate diffusion of a magnetic field, linear growth of the ETI, and a fully nonlinear ETI evolution. The nonlinear ETI simulations show early time agreement, with some differences emerging, as noted in the wavenumber spectrum, late into the nonlinear development of ETI. Here, a sensitivity study explores the role of equation-of-state (EOS), vacuum density, and vacuum resistivity. EOS and vacuum resistivity are found to be the most critical factors in the modeling of nonlinear ETI development.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

How to Get the Best Gas Separation Membranes from State-of-the-Art Glassy Polymers

Herein we focus on fundamental polymer science factors in membrane-based gas separation to open another chapter on this important topic. Realistically thinking about nano and angstrom scale matrix attributes that are translatable into higher performance thin glassy polymer asymmetric membrane layers in practical large-scale devices is crucial. Such thinking should be guided by expertise in both polymer and membrane communities to advance the state of the art. Polymer structures, membrane morphologies and effects of operating conditions are discussed here, using specific examples to illustrate key issues. Rubbery polymers are not the focus of this discussion, since they lack diffusive discrimination for size-similar penetrants, which glassy polymers provide and make them the focus here. Four scalable subtopics guide necessary thinking: (i) plasticization; (ii) anti-plasticization†2; (iii) dual-mode transport involving saturation of unrelaxed free volume, and (iv) non-uniform free volume and stress profiles in thin skin glassy polymer asymmetric or composite membranes. Furthermore, using these subtopics in the context of macromolecular science we discuss issues needed to expand the state-of-the-art in gas separation membranes.

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

Influence of Thermal Treatment on Preceramic Polymer Grafted Nanoparticle Network Formation: Implications for Thermal Protection Systems and Aerospace Propulsion Components

Preceramic polymer grafted nanoparticles (PCP GNPs) offer potential advantages in the production of polymer-derived ceramic composites compared to neat preceramic polymers (PCPs), such as reduced thermal shrinkage and increased char yield. In comparison to traditional composites (particles + free polymer), PCPs avoid issues of compatibility, agglomeration, and phase separation during processing and provide routes to control nanoparticle arrangement. Prior literature has established that the conversion efficiency of PCPs and PCP GNPs to ceramics and ceramic composites is linked to the thermal pretreatment of these precursor polymers (e.g., curing). In this paper, we investigate the material transformations that occur when exposing PCP GNPs to a low-temperature thermal treatment (≤250 °C) prior to pyrolysis. The PCP GNPs explored in this study possessed silica nanoparticle cores but were distinct in their PCP coronas, having poly(1,1-dimethylpropylsilane) (allyl-GNP)- or poly(1,1-dimethylbenzylethylsilane) (styryl-GNP)-grafted polymers. After undergoing low-temperature thermal treatment, these PCP GNPs exhibited increased char yield at 800 °C; however, unlike commercial PCPs (e.g., allylhydridopolycarbosilane (SMP-10)), obvious cross-linkable sites are not present in the PCP corona structure. Differential scanning calorimetry, oscillatory rheology, and X-ray photon correlation spectroscopy were utilized to elucidate the thermally induced material changes in the allyl- and styryl-GNPs. Thermally induced changes to the structure and behavior of the PCP GNPs were determined to be linked to the chemical structure of the grafted polycarbosilane chains. Styryl-GNPs experienced a curing event between 180 and 250 °C, effecting the formation of a network, which contributed to char yield improvement (at 800 °C). The nanoscale dynamics of this material also shows a transition from diffuse behavior to ballistic behavior because of the cross-linking event. Conversely, the allyl-GNPs did not exhibit a curing event under heat treatment but did show an increase in thermally induced physical cross-linking. The allyl-GNP physical network is stronger after the first cycle of thermal treatment, and it is considered that this contributed to the improvement in char yield post thermal treatment. The advancements in the understanding of the cross-linking behavior of these hybrid materials are expected to advance the application of these materials to turbine engine, advanced friction, and heat-shielding components.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Evaluating the diffusion of Kr in UO 2 and ADOPT TM using time-of-flight elastic recoil detection analysis (ToF-erda)

A combination of 300 keV 84 Kr ion implantation and Time-of-Flight Elastic Recoil Detection Analysis is utilized to investigate the diffusion of Kr in UO 2 and ADOPT TM fuels. Composition depth-profiles on the nanometer scale were obtained, both for as-implanted samples and after annealing at 800°C for 1 hour. Observed drifts in the 84 Kr profiles could be associated with short-range diffusion mechanisms. The approach employed here provides the possibility to make direct comparisons with atomistic scale modelling data, and can be of service as a separate effect test in line with the Accelerated Fuel Qualification initiative.

ADOPT UO2↗
DOE OSTI2023

On the minimum transport required to passively suppress runaway electrons in SPARC disruptions

In Izzo et al (2022 Nucl. Fusion 62 096029), state-of-the-art modeling of thermal and current quench (CQ) magnetohydrodynamics (MHD) coupled with a self-consistent evolution of runaway electron (RE) generation and transport showed that a non-axisymmetric (n = 1) in-vessel coil could passively prevent RE beam formation during disruptions in SPARC, a compact high-field tokamak projected to achieve a fusion gain Q > 2 in DT plasmas. However, such suppression requires finite transport of REs within magnetic islands and re-healed flux surfaces; conservatively assuming zero transport in these regions leads to an upper bound of RE current ~1MA compared to ~8.7MA of pre-disruption plasma current. Further investigation finds that core-localized electrons, within r/a < 0.4 and with kinetic energies ~0.2–15 MeV, contribute most to the RE plateau formation. Yet only a relatively small amount of transport, i.e. a diffusion coefficient ~18 m 2 s -1 , is needed in the core to fully mitigate these REs. Properly accounting for (a) the CQ electric field's effect on RE transport in islands and (b) the contribution of significant RE currents to disruption MHD may help achieve this.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Evaluation of the Effect of Burnup Acceleration on UO 2 Microstructure Evolution

Accelerated fuel qualification (AFQ) is a methodology by which new nuclear fuels are developed in an accelerated time frame compared with historical fuel qualification approaches. AFQ generally relies on high-fidelity physics-based modeling and simulation tools to adequately describe fuel performance as well as on revolutionary methods to accelerate burnup accumulation and collect relevant data more quickly. This report summarizes the use of advanced fuel modeling and simulation tools to evaluate microstructures from commercially irradiated fuel and microstructures from proposed MiniFuel irradiations, in which burnup accumulation is accelerated while prototypic temperature conditions are maintained. In this milestone, we used the mesoscale fuel performance code MARMOT to model the evolution of irradiated UO 2 microstructures and their potential restructuring at high burnup. The simulation conditions were informed by BISON models of both commercially irradiated fuel and MiniFuel. A first set of simulations investigated the recrystallization behavior of fully dense microstructures and showed full recrystallization at burnups as low as 52 MWd/kgU at 950°C. However, these simulations did not account for the presence of fission gas bubbles (predicted by BISON). Therefore, a second set of simulations including fission gas bubbles was performed and indicated that at the lowest temperature considered (650°C), the porous UO 2 microstructures have the highest total Gibbs free energies and are likely to recrystallize earlier than higher temperature cases (800 and 950°C), which agrees with high-burnup fuel characterization data. The results also showed that at lower temperature (650°C), the total free energies of the PWR fuel and MiniFuel microstructures are not significantly different. However, at the highest temperature (950°C), MiniFuel microstructures have a lower free energy than that of the PWR fuel microstructure. The competing effects between the temperature-dependent grain nucleation rate and the reduction of the free energy of the microstructure at higher temperature as a result of diffusion indicated that restructuring may occur at even higher temperatures than those considered in this study. In addition to the microstructure evolution modeling efforts, the burnup gradient across a single fuel specimen was also considered. This evaluation was for the VXF-15 position of the High Flux Isotope Reactor (HFIR) using the code suite HFIRCON, which was developed to automate the workflow for evaluating targets and fuel as they are irradiated in HFIR. The burnup gradient evaluation showed a dependence on both the axial and radial locations within the specimen, with a maximum difference of 1.7 between the inner and outermost radial layers. This relationship was further supported by considering the fission product speciation with respect to location within the specimen, which showed a higher concentration of 239 Pu, 240 Pu, and 241 Pu on the outer radial locations of the specimen than the center. The findings of the burnup and speciation evaluation show that some amount of self-shielding is occurring in the specimen when irradiated in the high-flux environment of HFIR; however, this impact is more pronounced for natural uranium when compared to 6% enrichment due to the higher ratio of 238 U in the specimen. Further analyses are required to understand the sensitivity of this gradient to spatial mesh and enrichment of the specimen.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Carbonate Management to Enable Energy- and Carbon-Efficient CO 2 Electrolysis (Final Technical Report)

The rapid growth and plummeting cost of solar energy have spurred growing interest in using CO 2 electrolysis to produce chemicals and fuels as an alternative to conventional petrochemical processes. High-temperature (>800 °C) solid oxide electrolyzers that convert CO 2 into CO and O 2 have recently become commercially available. Low-temperature electrolysis cells offer the prospect of more convenient and flexible operation, which is critical for utilizing intermittent solar energy, and provide access to more valuable C 2+ products such as ethylene, ethanol, and propanol. Over the past 10 years, research in this area has yielded substantial progress in both fundamental understanding of the requisite electrocatalytic reactions and design of prototype devices. Leveraging insights from fuel cells and membrane water electrolyzers, researchers have developed electrolysis cells with gas diffusion electrodes (GDE) that have demonstrated high CO 2 electrolysis current densities (>100 mA cm –2 ) as well as promising selectivity and stability. Despite these advances, the energy efficiency (electrical energy-to-product) and carbon efficiency (CO 2 -to-product) of low-temperature CO 2 electrolysis remain far too low for large-scale deployment. A preponderance of evidence indicates that the principal source of efficiency losses is the rapid and thermodynamically favorable reaction of CO 2 with hydroxide (OH – ) to form carbonate (CO 3 2– ). Carbonate formation imposes steady-state electrolysis conditions that result in large voltage and CO 2 losses for all known (photo)electrochemical CO 2 cells. While much current research remains focused on CO 2 reduction catalyst design, this largely overlooked CO 3 2– problem presents a fundamental scientific barrier to creating a viable electrochemical option for converting solar energy into chemicals and fuels. The project pursues an integrated, multi-PI research effort that establishes a fundamental science of CO 3 2– management. PI Kanan and Co-PI Mani’s contribution to the project is to evaluate strategies to mitigate the CO 3 2– problem by changing the properties of the electrolyte and the environment in which CO 2 reduction catalysis takes place. Experimental studies showed that electrolytes composed of a high concentration of both CO 3 2– and HCO 3 – , which serve as moderately alkaline buffers, improved the cell voltage by compared to all-HCO 3 – electrolytes, but these buffered systems still show substantial CO 2 uptake that reduces pH over time. Computational studies developed a homogenized model of a CO 2 reduction catalyst layer that permits a low-cost exploration of the high-dimensional parameter space associated with catalyst layers on gas diffusion electrodes. The model was validated by accurately reproducing experimental data for the related but simpler reaction of CO reduction and then used to probe the effects of catalyst layer architecture on CO 2 reduction. Minimizing the size of catalyst and hydrophobic domains in the catalyst layer is predicted to mitigate CO 3 2– formation and thereby enable prolonged operation at elevated pH. In support of the CO 2 electrolysis studies, a new method for rapidly prototyping electrochemical cells was developed and validated. The method uses a combination of 3D printing and electroless plating to generate conductive cell components for evaluating new cell designs. The carbonate problem encompasses mass transport processes and acid-base reactions that are relevant to many other electrochemical systems. Investigation of strategies to address the carbonate problem led to an additional line of inquiry into the physicochemical phenomena that determine the efficiency of electrochemical acid-base production, which has numerous applications in the broader field of carbon management. New strategies for using the supporting electrolyte to inhibit H + /OH – recombination in electrochemical acid-base production were evaluated, leading to the development of a novel acid-base producing system that eliminates the need for ion exchange membranes and exhibits promising efficiency and current densities for scalable applications.

25 ENERGY STORAGE↗
DOE OSTI2020

Impact of Laboratory-Accelerated Aging Methods to Study Alkali–Silica Reaction and Reinforcement Corrosion on the Properties of Concrete

This study focuses on two separate investigations of the main aging mechanisms: alkali–silica reactivity (ASR) and the corrosion of reinforcing steel (rebar) concrete, both of which may result in a premature failure to meet the serviceability or strength requirements of a concrete structure. However, these processes occur very slowly, spanning decades. The impact of direct chemical additives to fresh concrete to accelerate ASR and the corrosion of reinforcing steel on the fresh and hardened properties of the ensuing material are investigated to inform the potential use of chemicals in large-scale studies. The deterioration of reinforced concrete (RC) is determined by means of expansion, cracking, bulk diffusivity and surface resistivity measurements, and compressive, split tensile and flexural strength tests. The results indicate that the addition of sodium hydroxide and calcium chloride can effectively accelerate the crack formation and propagation in concrete due to ASR and the corrosion of rebar, respectively. The ASR-induced cracks maintained a constant crack width from 0.05 mm to 0.1 mm over the measurement period regardless of the intensity of aging acceleration. Adding 4% chloride by weight of cement for accelerating rebar corrosion resulted in an average crack that was 82% larger than in the case of ASR accelerated with the addition of sodium hydroxide. The addition of alkali resulted in an increase in early-age (7-day) strength. At a total alkali loading of 2.98 kg/m3, 3.84 kg/m3 and 5.57 kg/m3, the 28-day compressive strength of concrete decreased by 3%, 10% and 24%, respectively. Similarly, a higher early-age strength and a lower later-age strength was observed for the concrete in the presence of corrosive calcium chloride. The results from this research are expected to inform future studies on the long-term performance of RC structures under accelerated ASR and corrosion.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Parallelized domain decomposition for multi-dimensional Lagrangian random walk mass-transfer particle tracking schemes

Lagrangian particle tracking schemes allow a wide range of flow and transport processes to be simulated accurately, but a major challenge is numerically implementing the inter-particle interactions in an efficient manner. This article develops a multi-dimensional, parallelized domain decomposition (DDC) strategy for mass-transfer particle tracking (MTPT) methods in which particles exchange mass dynamically. We show that this can be efficiently parallelized by employing large numbers of CPU cores to accelerate run times. In order to validate the approach and our theoretical predictions we focus our efforts on a well-known benchmark problem with pure diffusion, where analytical solutions in any number of dimensions are well established. In this work, we investigate different procedures for “tiling” the domain in two and three dimensions (2-D and 3-D), as this type of formal DDC construction is currently limited to 1-D. An optimal tiling is prescribed based on physical problem parameters and the number of available CPU cores, as each tiling provides distinct results in both accuracy and run time. We further extend the most efficient technique to 3-D for comparison, leading to an analytical discussion of the effect of dimensionality on strategies for implementing DDC schemes. Increasing computational resources (cores) within the DDC method produces a trade-off between inter-node communication and on-node work. For an optimally subdivided diffusion problem, the 2-D parallelized algorithm achieves nearly perfect linear speedup in comparison with the serial run-up to around 2700 cores, reducing a 5 h simulation to 8 s, while the 3-D algorithm maintains appreciable speedup up to 1700 cores.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Quantitative nanoscale tracking of oxygen vacancy diffusion inside single ceria grains by in situ transmission electron microscopy

Oxygen vacancy formation and migration in ceria is critical to its electrochemical and catalytic properties in systems for chemical and energy transformation, but its quantification is rather challenging especially at atomic-scale because of disordered distribution. Here we report a rational approach to track oxygen vacancy diffusion in single grains of pure and Sm-doped ceria at -20 °C to 160 °C using in situ (scanning) transmission electron microscopy ((S)TEM). To create a gradient in oxygen vacancy concentration, a small region (~30 nm in diameter) inside a ceria grain is reduced to the C-type CeO 1.68 phase by the ionization or radiolysis effect of a high-energy electron beam. The evolution in oxygen vacancy concentration is then mapped through lattice expansion measurement using scanning nano-beam diffraction or 4D STEM at a spatial resolution better than 2 nm; this allows direct determination of local oxygen vacancy diffusion coefficients in a very small domain inside pure and Sm-doped ceria at different temperatures. Further, the activation energies for oxygen transport are determined to be 0.59, 0.66, 1.12, and 1.27 eV for pure CeO 2 , Ce 0.94 Sm 0.06 O 1.97 , Ce 0.89 Sm 0.11 O 1.945 , and Ce 0.8 Sm 0.2 O 1.9 , respectively, implying that activation energy increases due to impurity scattering. The results are qualitatively supported by density functional theory (DFT) calculations. In addition, our in situ TEM investigation reveals that dislocations impede oxygen vacancy diffusion by absorbing oxygen vacancies from the surrounding areas and pinning them locally. With more oxygen vacancies absorbed, dislocations show extended strain fields with local tensile zone sandwiched between the compressed ones. Therefore, dislocation density should be reduced in order to minimize the resistance to oxygen vacancy diffusion at low temperatures.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Tritium Transport Phenomena in Molten-Salt Reactors: Molten Salt Tritium Transport Experiment Design

Tritium is produced from neutron interactions with both lithium and beryllium. Large quantities of tritium are generated in Molten Salt Reactors (MSRs) which use LiF/BeF2 (FLiBe) as the fuel salt. Tritium is unique among the radionuclide hazards as it readily permeates through metal structural materials at high temperatures. All metal surfaces are potential release paths for tritium. For adequate safety analysis and eventual licensing of new reactors, predictive models for tritium transport and release from MSRs must be developed. These models must account for the multiple transport phenomena involved with tritium: fuel salt phase mass transport, dissociation/recombination reactions on metal surfaces, interstitial diffusion through the metal structure, and salt or gas phase mass transport in the downstream fluid. These models also must also be validated with representative experiments. Our previous report outlined tritium transport phenomena involved in MSRs, made suggestions on gaps in the transport dataset, and proposed an experimental test stand to test combined transport effects – tritium transport through pipe walls in a convective salt flow. In this report, we summarize an updated analysis framework for tritium transport in MSRs, report our results on hydrogen and deuterium permeation through Hastelloy N, and describe the final design of the Molten Salt Tritium Transport Experiment (MSTTE, pronounced “misty”). The MSRE provides the only wholistic experimental data set for tritium transport in MSRs and understanding the transport phenomena involved in the MSRE is crucial for future model development. One set of parameters in our analysis framework was unknown for the MSRE—surface reaction rates for tritium on Hastelloy N. This warranted our hydrogen and deuterium permeation campaign to assess the permeability, diffusivity, and solubility of hydrogen isotopes in clean Hastelloy N. Surface reaction rate constants were probed by low pressure measurements, however, no surface effects were observed in the limits of our permeation apparatus. Permeation experiments on oxidized Hastelloy N were not performed for this report but are planned in future work. The experimental test stand, MSTTE, measures combined transport properties of the salt-metal system. MSTTE is a forced convection FLiBe loop with custom designed test section to measure tritium transport through candidate structural materials. We use MSRE relevant dimensionless numbers to design and scale the test section. Hastelloy N is a candidate loop and test section material due to the relevance for the MSRE and related designs, however, other metals are being considered (e.g. 316H SS) which may better align with current vendor concepts.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2023

Phase-field modeling of solid-state metathesis reactions with the charge neutrality constraint

In this work, we present a phase-field model that captures the evolution of ionic concentrations and phase fractions during solid-state metathesis (SSM) reactions where diffusion limits the rate of transformation. The evolution of the mole fraction of each ion is obtained via governing equations that describe the reduction of a free energy, which includes an energy landscape with local minima located at compositions corresponding to stable products. We utilized two Lagrange multipliers to impose constraints of electroneutrality as well as on the sum of mole fractions, which were then eliminated to derive set of two partial differential equations that describe the dynamics of the mole fraction evolution. From these governing equations, the expressions for effective mobilities for the cations and the anions were obtained. We first study the effect of mobilities of ions on the reaction kinetics, using a simple model considering the ions with an identical absolute value of charge numbers. The simulation results show that the overall characteristic mobility, defined as the sum of the two effective ionic mobilities, provides an excellent measure of the rate at which reaction progresses and that the ratio of the effective mobilities of the anions and the cations signifies the manner by which the reaction progresses. We then generalize the model to consider ions with different charge numbers and tuned the mobility of ions based on their diffusion coefficients reported in the literature and experimental data from a thin-film experiment for the synthesis of FeS 2 to demonstrate the capability of the model to predict the phase evolution during SSM reactions. In particular, the simulation predicts nonplanar phase evolution, which is recently observed in thin-film reactions for the synthesis of FeS 2 via transmission electron microscopy. The approach can serve as a basis for models for phase transformations in other multiphase ionic mixtures, such as in all-solid-state batteries and in ionic liquids.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Diffusion of Brønsted acidic dopants in conjugated polymers

Many semiconductor devices (e.g., light emitting diodes and photovoltaics) utilize heterojunctions of doped and undoped layers or depend on gradients of electronic doping to control charge transport. Understanding of the formation and stability of gradients in doping requires an understanding of diffusion of dopants and the complex changes in polymer properties that arise during doping. Conjugated polymers can be electrically doped by strong acids, but the details of the reaction mechanism and subsequent stability are not understood. Here, we show a clear kinetic isotope effect in the doping of thin films of poly(3-hexylthiophene) (P3HT) by bis(trifluoromethane)sulfonimide (HTFSI) from solution indicating that this doping process is limited by proton transfer to the polymer. Complementary X-ray photoelectron spectroscopy and dynamic secondary ion mass spectrometry (DSIMS) depth profiling of dopant concentrations show definitive evidence of dopant enrichment at the P3HT surface. These surface-limited concentration profiles suggest that diffusivity of dopants vary inversely with dopant concentration due to doping-induced changes to the structure of the conjugated polymer.

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

Graphene oxide precursor effects on 3D-printed carbon scaffolds

Manganese oxide (MnO 2 ), an earth-abundant material, is a promising component for energy storage devices, with uses in both pseudocapacitors and batteries. However, high MnO 2 loading often leads to reduced performance due to poor ion diffusion. 3D printing, particularly using the direct ink writing (DIW) technique, offers a solution by enabling the fabrication of electrodes with hierarchical porous structures and open channels that enhance mass transport and ion diffusion. Previous work demonstrated that 3D-printed graphene aerogels with MnO 2 coatings exhibited excellent electrochemical performance, even with thick electrodes, due to their optimized structure. Building on this work, the current study investigates the performance differences between aerogels developed using graphene oxide (GO) and reduced graphene oxide (rGO) as carbon precursors. Both materials were incorporated into thixotropic inks, 3D-printed into lattice structures, and carbonized. Despite expected similarities between the final graphene aerogel, rGO-based aerogels exhibited superior areal capacitance, compared to GO-based aerogels. These differences are attributed to the lower oxygen content and defect density of rGO, which influence its interaction with cellulose viscosifiers in the ink formulation. Brunauer–Emmett–Teller (BET) surface area analysis revealed that rGO aerogels exhibit a larger surface area and mesoporous structure, further enhancing their performance. When coated with MnO 2 , rGO-based aerogels maintained their superior capacitive behavior over GO-based aerogels. This study highlights the effect of carbon precursor on the end performance of graphene aerogels.

Materials science↗
DOE OSTI2024

A Lake Biogeochemistry Model for Global Methane Emissions: Model Development, Site‐Level Validation, and Global Applicability

Abstract Lakes are important sentinels of climate change and may contribute over 30% of natural methane (CH 4 ) emissions; however, no earth system model (ESM) has represented lake CH 4 dynamics. To fill this gap, we refined a process‐based lake biogeochemical model to simulate global lake CH 4 emissions, including representation of lake bathymetry, oxic methane production (OMP), the effect of water level on ebullition, new non‐linear CH 4 oxidation kinetics, and the coupling of sediment carbon pools with in‐lake primary production and terrigenous carbon loadings. We compiled a lake CH 4 data set for model validation. The model shows promising performance in capturing the seasonal and inter‐annual variabilities of CH 4 emissions at 10 representative lakes for different lake types and the variations in mean annual CH 4 emissions among 106 lakes across the globe. The model reproduces the variations of the observed surface CH 4 diffusion and ebullition along the gradients of lake latitude, depth, and surface area. The results suggest that OMP could play an important role in surface CH 4 diffusion, and its relative importance is higher in less productive and/or deeper lakes. The model performance is improved for capturing CH 4 outgassing events in non‐floodplain lakes and the seasonal variability of CH 4 ebullition in floodplain lakes by representing the effect of water level on ebullition. The model can be integrated into ESMs to constrain global lake CH 4 emissions and climate‐CH 4 feedback.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

A direct numerical simulation of Jet A flame kernel quenching

The safe operation of aeronautical engines requires an understanding of flame ignition, propagation and extinction. In this study, direct numerical simulations are performed using a 29 species reduced chemical mechanism for jet fuel surrogate Jet A to understand the flame quenching process. Here, initially laminar spherical flames of varying sizes and equivalence ratios are subject to an identical periodic domain of decaying and isotropic high intensity turbulence with a turbulent Reynolds number of 2400. All cases become quenched, except for the larger kernel with lower Karlovitz number. An analysis of the flame structure shows broadened preheat zone, flame shortening on the product side, differential species diffusion and partial fuel pyrolysis in the fresh mixture. Two extinction mechanisms are identified arising from flame shortening and high flame stretch. Flame shortening occurs due to turbulence-chemistry interactions that resemble the flame–flame interaction in a laminar counterflow reactant-to-reactant configuration, which contorts and breaks up the ignition kernel. Flame stretch is a local effect that attenuates the heat release rate and causes the flame to retreat towards the product mixtures, similar to what has been observed for reactant-to-product laminar counterflow flames. Chemical explosive mode analysis was also performed to quantify the flame structure and local combustion mode. The diffusion–reaction balance in pinched-off flame islands favors extinction of these smaller structures, while auto-ignition modes are observed within the flame kernel after fresh mixture is engulfed and preheated in the product kernel. Statistics of the density-weighted displacement speed conditional on local combustion mode indicates strong correlation between the local extinction mode and negative displacement speed. The local balance between diffusion and reaction ultimately determines the propensity for local extinction in both laminar and turbulent flames, the extent of which has an impact on global flame propagation.

42 ENGINEERING↗
DOE OSTI2025

The Effect of Luminance Pattern on Nighttime Discomfort Glare Response (CRADA Final Report)

LED luminaires with apertures containing visually resolvable “bright” (i.e., high luminance) spots are commonly used to illuminate outdoor environments. However, the impact of the distance between LEDs, which determines spatial frequency at a given viewing distance, on perceived discomfort caused by these luminaires remains unclear. This study involved 29 participants who were each shown 68 stimuli varying in spatial frequency, intensity, diffusion level, and ambient lighting in a dark laboratory setting. Participants reported their experiences of discomfort from glare, ability to resolve individual LEDs, and afterimages. The findings revealed that increased spatial frequency heightened both discomfort and the likelihood of experiencing afterimages. It is hypothesized that the point spread function of the eye contributed to these effects, where sources began to be perceived as a single larger source with equal or greater intensity. These results suggest that LED luminaire designs should favor configurations with lower spatial frequencies to minimize discomfort. A simple quantity, direct illuminance at the eye, is recommended for use to predict discomfort from glare because it performed like other more complex models.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Numerical Study of PAHs and Soot Emissions from Gasoline–Methanol, Gasoline–Ethanol, and Gasoline– n -Butanol Blend Surrogates

Soot formation is an intricate phenomenon, and soot propensity of a fuel is interwoven with the fuel composition, physical and chemical properties, and combustion environment. Here, the present study examines the hypothesis that in addition to the chemical composition of the fuel, the sooting nature of the fuel is closely coupled with its chemical property known as octane sensitivity (S). With this motivation, the present study numerically investigates the effects of gasoline surrogate composition and its property, octane sensitivity (S), on polycyclic aromatic hydrocarbons (PAHs) and soot emissions. Four-component toluene primary reference fuel (TPRF)–alcohol blends, comprising iso-octane, n-heptane, toluene, and one of the three different alcohols- methanol, ethanol, and n-butanol, are used as gasoline surrogates. A total of 320 TPRF–alcohol mixtures, with S in the range of 1–10, are examined under laminar counterflow diffusion flame conditions. A detailed chemical mechanism coupled with a comprehensive soot model, which includes reactions for soot inception, surface growth, PAH condensation, and oxidation, is adopted. The analysis indicates that the toluene content in the fuel mixture has a prominent effect, while the alcohol content and octane sensitivity of the fuel have a weak correlation with the PAHs and soot. Thus, it is not clear if any of these three variables, namely, toluene content in the fuel, alcohol content in the fuel, and S, are individually sufficient to characterize the PAHs and soot across various blends. For this reason, a new variable (X CHO ) based on the elemental composition of the fuel mixture is identified and it is shown that X CHO along with S of the fuel characterize soot emissions satisfactorily. Further, a reaction path analysis indicates that the efficacy of alcohols in reducing soot emissions follows the order: methanol > ethanol > n-butanol.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗