Search NASA⌕ Search

SEARCH · Search NASA

Results for “diffusion limitation”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Experimental investigation of elemental and isotopic evaporation processes by laser heating in an aerodynamic levitation furnace

In this report we carried out evaporation experiments on a B-type calcium–aluminium-rich inclusion (CAI) melt in a gas-mixing aerodynamic levitation laser furnace, at 1873 K and an oxygen partial pressure of 10 -9.1 atm, for durations ranging from 60 to 600 s. Evaporation of SiO 2 and MgO follow the same trend as those observed in vacuum furnace experiments at the same temperature and starting composition, showing that their evaporation relative to one another from the melt is independent of pressure, oxygen fugacity, and hydrodynamical regime specific to the furnace. Isotopic ratios of Mg and Si in evaporation residues are used to derive fractionation factors of α 26/24 Mg vap-liq = 0.9906± 0.0004 and α 30/28 Si vap-lip = 0.9943±0.0003, which are both significantly closer to unity than those found for evaporation in a vacuum, which translates to less isotope fractionation. The residues are also less isotopically fractionated than expected for cases in which transport of the gas species away from the melt is diffusion-controlled at 1-atm. By analysing the flow regimes in our furnace, we find that advection by the levitating gas is the primary mode of mass transport away from the melt surface, as opposed to diffusion-limited transport in a vacuum or 1-atm tube furnace. A modified Hertz–Knudsen–Langmuir formulation accounts for this process, and shows that isotopic fractionation of both Si and Mg reflect a saturation factor (ratio of the pressure of the evaporating species to vapour saturation pressure) equal to 0.75. This is in perfect accord with recent measurements of Cu isotopic fractionation using a similar furnace. The fact that three elements (Mg, Si, Cu) with varying equilibrium vapour pressures, activity coefficients in the liquid, and diffusion coefficients in the gas have the same scaling behaviour to saturation pressure is a strong indication that the mechanism controlling evaporation is driven by the hydrodynamical regime imposed in the furnace. Therefore, this class of experiments can be used to constrain processes in which advection dominates over diffusion, such as (but not limited to) planetary ejecta, tektites, giant impacts, nebular condensation in a turbulent flow, or nuclear fallout material. Finally, the possibility to reach high temperatures (in excess of 3500 K) in this furnace allows it to be used to evaluate the activity coefficients of melt components in extreme conditions relevant to molten planetary interiors (i.e., magma oceans), with a specific focus on refractory elements.

58 GEOSCIENCES↗
DOE OSTI2021

Recombination rate analysis in long minority carrier lifetime mid-wave infrared InGaAs/InAsSb superlattices

Gallium is incorporated into the strain-balanced In(Ga)As/InAsSb superlattice system to achieve the same mid-wave infrared cutoff tunability as conventional Ga-free InAs/InAsSb type-II superlattices, but with an additional degree of design freedom to enable optimization of absorption and transport properties. Time-resolved photoluminescence measurements of InGaAs/InAsSb superlattice characterization- and doped device structures are reported from 77 to 300 K and compared to InAs/InAsSb. The low-injection photoluminescence decay yields the minority carrier lifetime, which is analyzed with a recombination rate model, enabling the determination of the temperature-dependent Shockley–Read–Hall, radiative, and Auger recombination lifetimes and extraction of defect energy levels and capture cross section defect concentration products. The Shockley–Read–Hall-limited lifetime of undoped InGaAs/InAsSb is marginally reduced from 2.3 to 1.4 μs due to the inclusion of Ga; however, given that Ga improves the vertical hole mobility by a factor of >10×, a diffusion-limited InGaAs/InAsSb superlattice nBn could expect a lower bound of 2.5× improvement in diffusion length with significant impact on photodetector quantum efficiency and radiation hardness. At temperatures below 120 K, the doped device structures are Shockley–Read–Hall limited at 0.5 μs, which shows promise for detector applications.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Moisture-driven CO 2 pump for direct air capture

Direct Air Capture, the chemical separation of carbon dioxide (CO 2 ) from air, is considered a necessary approach to generate negative carbon emissions and in turn limit global warming, though energetic and monetary costs must continue to decline to meet the needed gigaton scale. The promise of energy-efficient and continuous carbon dioxide membrane separation has only recently been investigated with the development of facilitated chemical transport and novel electrochemically driven methods. This work demonstrates the potential for continuous pumping of CO 2 from air using commercial anion exchange membranes driven solely by a water vapor gradient. A reactive transport model that couples ionic transport with moisture sensitive bicarbonate chemistry in charged polymers, also termed moisture swing sorption, was developed alongside experimentally measured equilibrium and transport properties of CO 2 and water vapor in Fumasep FAA-3 to determine the mass transport limits of such a system. The results show that the commercial membrane is kinetically limited by the moisture sensitive chemistry, enabling a CO 2 current of 1.1 μmol/m 2 • s. In diffusion limited materials, a CO 2 pumping flux can reach 11 μmol/m 2 • s in and peak performance is expected with a dry air relative humidity around 50% pumping into a saturated water vapor permeate. In all configurations, water loss becomes the major cost in such a separation and must be managed through optimal polymer design or water recovery mechanisms. The work demonstrates the applicability of a new driving force, the evaporation of water, to meet the energetic demand of CO 2 separation from air.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Revisiting diffusion-based moisture functions: why do they fail?

Microbial activities in soil are controlled by the balance between physical and biological mechanisms that are currently implied by soil moisture-respiration relationships, but not explicitly represented. Diffusion based moisture functions focus on physical controls, while substrate kinetics based models primarily address biological drivers. How to harmonize these two types of models for robust predictions of soil respiration and other processes has not been addressed. Here we derived a generalizable analytical solution of diffusion-limited nonlinear microbial uptake kinetics, in which microbial metabolic rates were governed by both physical and biological drivers via substrate diffusion and microbial affinity. We proposed a soil trait concept to integrate both physical and biological perspectives into one unified function, and tested it using both laboratory and field based observations from a variety of soils. This approach provides a tractable and scalable approach that incorporates soil heterogeneity and functional complexity in macroscopic model descriptions.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Hydrogen and deuterium permeation in Hastelloy N

Hastelloy N was chosen as the fluoride salt-contacting structural material for the Molten Salt Reactor Experiment due to its excellent compatibility with the fuel salt FLiBe. FLiBe is currently investigated for several advanced fusion and fission reactor concepts where tritium generation in the FLiBe is anticipated. Knowledge of hydrogen transport properties through Hastelloy N is important to understand how tritium would permeate through this material and result in an unintentional release. In this study, the hydrogen and deuterium permeability, diffusivity, and solubility were measured from 500 to 700 °C at a primary-side pressure of 10 kPa in a well-characterized sample of Hastelloy N. The prepared polycrystalline Hastelloy N had C and O impurities present on the surface. These impurities were investigated using Auger Emission Spectroscopy and Ar depth profiling. Here, the adventitious C was removed upon the first Ar sputter cycle whereas O persisted deeper into the sample. For permeation experiments, applied deuterium pressures ranged from 13 Pa to 100 kPa and deuterium transport remained in the diffusion-limited regime (J ∝ P 0.5 ) throughout the pressure range examined. Two methods are employed to measure the effective hydrogen and deuterium diffusivity: rise and decline. The decline method produced improved statistical model fits for calculating the effective diffusion coefficient compared to the rise method. The resultant transport properties compared well to published values for other nickel alloys.

08 HYDROGEN↗
DOE OSTI2021

Resonant amplification of enzymatic chemical oscillations by oscillating flow

Using theory and simulation, we analyzed the resonant amplification of chemical oscillations that occur due to externally imposed oscillatory fluid flows. The chemical reactions are promoted by two enzyme-coated patches located sequentially on the inner surface of a pipe that transports the enclosed chemical solution. In the case of diffusion-limited systems, the period of oscillations in chemical reaction networks is determined by the rate of the chemical transport, which is diffusive in nature and, therefore, can be effectively accelerated by the imposed fluid flows. In this work, we first identify the natural frequencies of the chemical oscillations in the unperturbed reaction–diffusion system and, then, use the frequencies as a forcing input to drive the system to resonance. We demonstrate that flow-induced resonance can be used to amplify the amplitude of the chemical oscillations and to synchronize their frequency to the external forcing. In particular, we show that even 10% perturbations in the flow velocities can double the amplitude of the resulting chemical oscillations. Particularly, effective control can be achieved for the two-step chemical reactions where during the first half-period, the fluid flow accelerates the chemical flux toward the second catalytic patch, while during the second half-period, the flow amplifies the flux to the first patch. The results can provide design rules for regulating the dynamics of coupled reaction–diffusion processes and can facilitate the development of chemical reaction networks that act as chemical clocks.

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

Evolution of intermetallic phases in an Al–Si–Ti alloy during solution treatment

A cast Al–Si–Ti alloy was solution treated at 540 °C for different periods between 0 and 72 h to understand the evolution of intermetallic phases. Only an (Al,Si) 3 Ti intermetallic phase with a low Si content was found in the as-cast alloy. The (Al,Si) 3 Ti particles were converted partly into a lamellar structure, a eutectoid phase, consisting of a Si-rich phase and an Al phase during solution treatment. The amount of the lamellar structure increased with the solution treatment time, but the composition of either constitute was kept almost unchanged regardless of the solution treatment times. The lamellar Si-rich phase is (Al,Si) 2 Ti (τ 2 ) with a TiSi 2 (C49-type) structure based on thermodynamic calculations and high-resolution TEM analyses. FCC Al phase is the product residing between the τ 2 lamellae after the completion of transformation from (Al,Si) 3 Ti to τ 2 phase. A near-rational orientation relationship (OR) between the Al and τ 2 phases is determined as Al [110]//τ 2 [100], Al ()//τ 2 [060]. The phase transformation from (Al,Si) 3 Ti to τ 2 being a result of the diffusion of Ti and Si within the original (Al,Si) 3 Ti particulates as well as the Si diffusion from the Al matrix during solution treatment is proposed. The formation of the lamellar structure in the microstructure is attributed mainly to the limited diffusivity of Ti element.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Elemental zoning enhances mass transport in zeolite catalysts for methanol to hydrocarbons

Mass transport limitations in zeolite catalysts pose major hurdles for their optimal performance in diverse chemical reactions. Most approaches to reduce these restrictions focus on the synthesis of either hierarchical or nanosized zeolites. Here we demonstrate that the existence of a siliceous, catalytically inactive exterior rim on ZSM-5 particles dramatically reduces the diffusion limitations, which leads to an enhanced catalyst lifetime for the methanol-to-hydrocarbon reaction. In this study, our findings reveal that binary inorganic and organic structure-directing agents enable a one-pot synthesis of Si-zoned ZSM-5 catalysts with diffusion properties that are characteristic of particles with a much smaller size. Operando ultraviolet–visible light diffuse reflectance spectroscopy reveals a marked reduction in external coking among Si-zoned samples. Molecular dynamics simulations to assess the diffusion of methanol and benzene in siliceous pores and in those with Brønsted acids reveal substantially reduced transport limitations in zoned regions, consistent with the improved catalyst activity of Si-zoned zeolites relative to that of ZSM-5 with a homogeneous acid-site distribution.

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

Fundamentals of Au(111) Surface Dynamics: Coarsening of Two-Dimensional Au Islands

Au(111) surfaces play a central role in many applications, yet studies of fundamental aspects of their dynamics are limited. Thus, using scanning tunneling microscopy (STM) at 300 K, we analyze the coarsening of first-layer 2D Au islands directly on the Au(111) substrate and also of second-layer 2D Au islands. Specifically, we monitor the decay of Au first-layer islands with areas of about 100–500 nm 2 in the vicinity of larger islands or extended step edges over a period of approximately 40 h, the relevant time scale for this process. Experimentally observed behavior is captured by analytical theory for terrace-diffusion-limited decay incorporating DFT results for the Au terrace diffusion barrier and the adatom formation energy.Lastly, experimental observations of second-layer island decay are also compared with appropriate analytical theory and stochastic simulations, thereby determining the effective Ehrlich–Schwoebel barrier for Au on Au(111).

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

Dynamic and static corrosion of chromium-alumina refractory in simulated nuclear waste glass

Monofrax® K-3, a chromium–alumina refractory, is widely used in nuclear waste glass melters due to its high durability. However, aggressive vitrification conditions still lead to corrosion that limits melter lifespan. Here, this study investigates two key corrosion modes, subsurface and melt-line corrosion, under static and dynamic conditions, using a custom setup enabling precise control of melt velocity and temperature. Experimental results, interpreted using diffusion-limited dissolution models, show that (i) subsurface corrosion increases by approximately 10% per mm/s increase in melt velocity, while melt-line corrosion remains virtually unaffected, (ii) corrosion rate dependence on temperature is inversely proportional to melt viscosity, and (iii) glass composition affects both corrosion modes similarly, based on Cr 2 O 3 , Al 2 O 3 , and alkali contents. The experimental data also allowed the determination of diffusion coefficients of major dissolving species, which were found to be consistent between the melt-line and subsurface corrosion models, indicating their suitability for future CFD studies.

12 - MGMT OF RADIOACTIVE AND NON-RADIOACTIVE WASTE↗
DOE OSTI2022

Mass Transport Limitations and Kinetic Consequences of Corn Stover Deacetylation

Alkaline pretreatment of herbaceous feedstocks such as corn stover prior to mechanical refining and enzymatic saccharification improves downstream sugar yields by removing acetyl moieties from hemicellulose. However, the relationship between transport phenomena and deacetylation kinetics is virtually unknown for such feedstocks and this pretreatment process. Here, we report the development of an experimentally validated reaction–diffusion model for the deacetylation of corn stover. A tissue-specific transport model is used to estimate transport-independent kinetic rate constants for the reactive extraction of acetate, hemicellulose and lignin from corn stover under representative alkaline conditions (5–7 g L -1 NaOH, 10 wt% solids loadings) and at low to mild temperatures (4–70°C) selected to attenuate individual component extraction rates under differential kinetic regimes. The underlying transport model is based on microstructural characteristics of corn stover derived from statistically meaningful geometric particle and pore measurements. These physical descriptors are incorporated into distinct particle models of the three major anatomical fractions (cobs, husks and stalks) alongside an unsorted, aggregate corn stover particle, capturing average Feret lengths of 917–1239 μm and length-to-width aspect ratios of 1.8–2.9 for this highly heterogeneous feedstock. Individual reaction–diffusion models and their resulting particle model ensembles are used to validate and predict anatomically-specific and bulk feedstock performance under kinetic-controlled vs. diffusion-controlled regimes. In general, deacetylation kinetics and mass transfer processes are predicted to compete on similar time and length scales, emphasizing the significance of intraparticle transport phenomena. Critically, we predict that typical corn stover particles as small as ~2.3 mm in length are entirely diffusion-limited for acetate extraction, with experimental effectiveness factors calculated to be 0.50 for such processes. Debilitatingly low effectiveness factors of 0.021–0.054 are uncovered for cobs—implying that intraparticle mass transfer resistances may impair observable kinetic measurements of this anatomical fraction by up to 98%. These first-reported quantitative maps of reaction vs. diffusion control link fundamental insights into corn stover anatomy, biopolymer composition, practical size reduction thresholds and their kinetic consequences. These results offer a guidepost for industrial deacetylation reactor design, scale-up and feedstock selection, further establishing deacetylation as a viable biorefinery pretreatment for the conversion of lignocellulosics into value-added fuels and chemicals.

09 BIOMASS FUELS↗
DOE OSTI2021

Representing methane emissions from wet tropical forest soils using microbial functional groups constrained by soil diffusivity

Tropical ecosystems contribute significantly to global emissions of methane (CH 4 ), and landscape topography influences the rate of CH 4 emissions from wet tropical forest soils. However, extreme events such as drought can alter normal topographic patterns of emissions. Here we explain the dynamics of CH 4 emissions during normal and drought conditions across a catena in the Luquillo Experimental Forest, Puerto Rico. Valley soils served as the major source of CH 4 emissions in a normal precipitation year (2016), but drought recovery in 2015 resulted in dramatic pulses in CH 4 emissions from all topographic positions. Geochemical parameters including (i) dissolved organic carbon (C), acetate, and soil pH and (ii) hydrological parameters like soil moisture and oxygen (O 2 ) concentrations varied across the catena. During the drought, soil moisture decreased in the slope and ridge, and O 2 concentrations increased in the valley. We simulated the dynamics of CH 4 emissions with the Microbial Model for Methane Dynamics-Dual Arrhenius and Michaelis–Menten (M3D-DAMM), which couples a microbial functional group CH 4 model with a diffusivity module for solute and gas transport within soil microsites. Contrasting patterns of soil moisture, O 2 , acetate, and associated changes in soil pH with topography regulated simulated CH 4 emissions, but emissions were also altered by rate-limited diffusion in soil microsites. Changes in simulated available substrate for CH 4 production (acetate, CO 2 , and H 2 ) and oxidation (O 2 and CH 4 ) increased the predicted biomass of methanotrophs during the drought event and methanogens during drought recovery, which in turn affected net emissions of CH 4 . A variance-based sensitivity analysis suggested that parameters related to aceticlastic methanogenesis and methanotrophy were most critical to simulate net CH 4 emissions. This study enhanced the predictive capability for CH 4 emissions associated with complex topography and drought in wet tropical forest soils.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Modeling the Size Distribution and Chemical Composition of Secondary Organic Aerosols during the Reactive Uptake of Isoprene-Derived Epoxydiols under Low-Humidity Condition

Reactive uptake of isoprene epoxydiols (IEPOX), which are isoprene oxidation products, onto acidic sulfate aerosols is recognized to be an important mechanism for the formation of isoprene-derived secondary organic aerosol (SOA). While a mechanistic understanding of IEPOX-SOA formation exists, several processes affecting their formation remain uncertain. Evaluating mechanistic IEPOX-SOA models with controlled laboratory experiments under longer atmospherically relevant timescales is critical. Here, we implement our latest understanding of IEPOX-SOA formation within a box model to simulate the measured reactive uptake of IEPOX on polydisperse ammonium bisulfate seed aerosols within an environmental Teflon chamber. The model is evaluated with single-particle measurements of size distribution, volume, density, and composition of aerosols due to IEPOX-SOA formation at timescales of hours. Additionally, we find that the model can simulate the growth of particles due to IEPOX multiphase chemistry, as reflected in increases of the mean particle size and volume concentrations, and a shift of number size distribution to larger sizes. The model also predicts the observed evolution of particle number mean diameter and total volume concentrations at the end of the experiment. We show that in addition to self-limiting effects of IEPOX-SOA coatings, accounting for the molar balance between inorganic and organic sulfate and mass accommodation coefficient of IEPOX are important parameters governing modeling of IEPOX-SOA formation. Thus, models which do not account for the molar sulfate balance and/or diffusion limitations within IEPOX-SOA coatings are likely to predict too high IEPOX-SOA formation.

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

An electrochemical mesoscale tool for modeling the corrosion of structural alloys by molten salt

Understanding the impact of microstructure on corrosion rates can aid the development of corrosion-resistant alloys for molten salt reactors. Here in this work, we develop an electrochemical phase-field model for capturing the microstructure-dependent corrosion of structural alloys by molten salts. As a demonstration problem, we apply this model to capture the selective depletion of Cr from Ni-Cr grain boundaries during corrosion in molten FLiBe salt. We perform sensitivity analysis and model verification on 1D simulations to confirm that the model predicts diffusion-limited kinetics. The model is validated using 1D, 2D, and 3D simulations against experimental data for Ni-5Cr and Ni-20Cr corrosion in molten FLiBe. The 1D simulations predict the corrosion behavior with reasonable accuracy when using an effective diffusion coefficient that accurately represents the grain boundary diffusion. 2D simulations that represent the grain structure underpredict the corrosion. 3D simulations that represent the grain structure predict the corrosion with reasonable accuracy. The corrosion rate predicted by the 3D simulations is proportional to the average grain size at the alloy/salt interface.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Impact of Site Identity, Location, and Accessibility on Polyethylene Conversion Rates and Product Selectivities over Metal-Free MFI Zeolites

Brønsted (BAS), Lewis (LAS), and surface Brønsted (SBAS) acid sites have been investigated for polyethylene (PE) upcycling by zeolite catalysts, but there is no clear consensus regarding their catalytic roles, partly due to the complexity of the catalysts used and varying reaction conditions across studies. This work systematically determined how these sites impact PE conversion rates and product distributions by utilizing a suite of microporous MFI catalysts with varying Si/Al ratios, acid site densities, and inherent mesoporosities. PE conversion rates did not trend with total BAS or LAS densities due to a combination of internal mass transfer limitation and the apparent inability of LAS alone to cleave C–C bonds, but a strong, statistically significant correlation was present with respect to SBAS density and mesopore surface area, jointly, owing to accelerated polymer activation on external surfaces to smaller diffusion-limited chains. However, ingress of these SBAS-derived fragments ultimately remained rate limiting, as demonstrated by solid conversion rates that increased with mesopore surface area at similar SBAS density and likewise increased with SBAS density at similar mesopore surface area. In batch PE cracking reactions, light gaseous product selectivities were most sensitive to total BAS, with higher densities generally exhibiting higher selectivity to C 3 and linear C 4 –C 7 products and higher alkane/alkene product ratios, consistent with increased β-scission turnovers. Insights from this work help systematically clarify the roles of BAS, LAS, SBAS, and mesopores in PE cracking reactions and inform the development of tailored zeolite catalysts for efficient polyolefin upcycling.

accessibility↗
DOE OSTI2024

Host and nonhost bacteria support bacteriophage dissemination along mycelia and abiotic dispersal networks

Abstract Bacteriophages play a crucial role in shaping bacterial communities, yet the mechanisms by which nonmotile bacteriophages interact with their hosts remain poorly understood. This knowledge gap is especially pronounced in structured environments like soil, where spatial constraints and air-filled zones hinder aqueous diffusion. In soil, hyphae of filamentous microorganisms form a network of ‘fungal highways’ (FHs) that facilitate the dispersal of other microorganisms. We propose that FHs also promote bacteriophage dissemination. Viral particles can diffuse in liquid films surrounding hyphae or be transported by infectable (host) or uninfectable (nonhost) bacterial carriers coexisting on FH networks. To test this, two bacteriophages that infect Pseudomonas putida DSM291 (host) but not KT2440 (nonhost) were used. In the absence of carriers, bacteriophages showed limited diffusion on 3D-printed abiotic networks, but diffusion was significantly improved in Pythium ultimum-formed FHs when the number of connecting hyphae exceeded 20. Transport by both host and nonhost carriers enhanced bacteriophage dissemination. Host carriers were five times more effective in transporting bacteriophages, particularly in FHs with over 30 connecting hyphae. This study enhances our understanding of bacteriophage dissemination in nonsaturated environments like soils, highlighting the importance of biotic networks and bacterial hosts in facilitating this process.

Périat, Claire (ORCID:000000020825625X)↗
DOE OSTI2023

Re-evaluating cloud chamber constraints on depositional ice growth in cirrus clouds – Part 1: Model description and sensitivity tests

Abstract. Ice growth from vapor deposition is an important process for the evolution of cirrus clouds, but the physics of depositional ice growth at the low temperatures (<235 K) characteristic of the upper troposphere and lower stratosphere is not well understood. Surface attachment kinetics, generally parameterized as a deposition coefficient αD, control ice crystal habit and also may limit growth rates in certain cases, but significant discrepancies between experimental measurements have not been satisfactorily explained. Experiments on single ice crystals have previously indicated the deposition coefficient is a function of temperature and supersaturation, consistent with growth mechanisms controlled by the crystal's surface characteristics. Here we use observations from cloud chamber experiments in the Aerosol Interactions and Dynamics in the Atmosphere (AIDA) aerosol and cloud chamber to evaluate surface kinetic models in realistic cirrus conditions. These experiments have rapidly changing temperature, pressure, and ice supersaturation such that depositional ice growth may evolve from diffusion limited to surface kinetics limited over the course of a single experiment. In Part 1, we describe the adaptation of a Lagrangian parcel model with the Diffusion Surface Kinetics Ice Crystal Evolution (DiSKICE) model (Zhang and Harrington, 2014) to the AIDA chamber experiments. We compare the observed ice water content and saturation ratios to that derived under varying assumptions for ice surface growth mechanisms for experiments simulating ice clouds between 180 and 235 K and pressures between 150 and 300 hPa. We found that both heterogeneous and homogeneous nucleation experiments at higher temperatures (>205 K) could generally be modeled consistently with either a constant deposition coefficient or the DiSKICE model assuming growth on isometric crystals via abundant surface dislocations. Lower-temperature experiments showed more significant deviations from any depositional growth model, with different ice growth rates for heterogeneous and homogeneous nucleation experiments.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

3D simulations of spinlike flames in Co/Al multilayers with enhanced conduction losses

Reactive Co/Al multilayers are uniformly structured materials that may be ignited to produce rapid and localized heating. Prior studies varying the bilayer thickness (i.e., sum of two individual layers of Co and Al) have revealed different types of flame morphologies, including: (a) steady/planar, (b) wavy/periodic, and (c) transverse bands, originating in the flame front. These instabilities resemble the “spin waves” first observed in the early studies of solid combustion (i.e., Ti cylinder in a N 2 atmosphere), and are likewise thought to be due to the balance of heat released by reaction and heat conduction forward into the unreacted multilayer. However, the multilayer geometry and three-dimensional (3D) edge effects are relatively unexplored. In this work, a new diffusion-limited reaction model for Co/Al multilayers was implemented in large, novel 3D finite element analysis (FEA) simulations, in order to study the origins of these spinlike flames. This reaction model builds upon previous work by introducing three new phase-dependent property models for: (1) the diffusion coefficient, (2) anisotropic thermal conductivity tensor, and (3) bulk heat capacity, as well as one additional model for the bilayer-dependent heat of reaction. These novel 3D simulations are the first to predict both steady and unsteady flames in Co/Al multilayers. Moreover, two unsteady modes of flame propagation are identified, which depend on the enhanced conduction losses with slower flames, as well as flame propagation around notched edges. Finally, future work will consider the generality of the current modeling approach and also seek to define a more generalized set of stability criteria for additional multilayer systems.

42 ENGINEERING↗