Search NASA⌕ Search

SEARCH · Search NASA

Results for “high capture rate”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Membrane Development for CO2 Capture from Steel Manufacturing

The U.S. government is targeting a net-zero carbon-emission economy by 2050, offering an exciting opportunity for membrane-based CO2 capture from various industrial point sources. Given that industrial flue gas has low CO2 partial pressures and high volumetric flow rates, high-permeance membranes are needed to make membrane technology economically viable for large-scale deployment. Thin film composite (TFC) membranes are necessary for this implementation because they can provide high permeance by forming a thin selective layer on top of a porous support. This presentation will report the rational design and fabrication of NETL’s highly permeable non-aging TFCs achieved by: synthesizing a high-performance rubbery selective material; developing a high-porosity membrane support; optimizing coating methods to assemble the two materials into scalable membranes; and scaling up membrane supports and TFCs via a roll-to-roll process. The novel rubbery selective material shows mixed-gas CO2 permeability of 930 Barrer and CO2/N2 selectivity of 44, exceeding the 2008 Robeson upper bound. The resulting TFCs yield remarkably high CO2 permeance of 4,500 GPU and CO2/N2 selectivity of 34 at 23ºC. Moreover, the TFCs exhibit not only excellent performance stability (or non-aging behavior) for 1,000 hours in the lab, but also maintain their separation properties in a 700-hour field test at the U.S. DOE’s National Carbon Capture Center using real humid flue gas.

Tran, Thien↗
DOE OSTI2025

Advanced Computational Modeling of High-Level Waste Vitrification at the Hanford Site

The U.S. Department of Energy (DOE) has selected vitrification for stabilizing legacy tank waste at the Hanford site, where radioactive waste from plutonium production was historically stored in underground tanks. This waste will be separated into low-activity waste (LAW) and high-level waste (HLW) fractions and processed at the Waste Treatment and Immobilization Plant (WTP). At WTP, glass melters are used for the vitrification of radioactive tank waste, transforming it into a stable borosilicate glass form for safe long-term storage. The melter vessel is constructed from highly durable and heat-resistant materials, where the vitrification process occurs. The main regions that are modeled are the melt pool, plenum, cold cap, riser/discharge chamber, and surrounding structure with insulation layers. Forced convection induced by air bubblers at the base of the melter ensure uniform temperature distribution and provide heat to the cold cap layer. The cold cap is a region of reacting batch feed that floats on top of the molten glass and is where the batch-to-glass reactions occur. Joule heating provided by electrodes mounted along the vertical walls of the melter and immersed directly in the glass, generates the necessary heat for the net endothermic conversion processes that occur in the cold cap. The high temperatures, radioactivity, and opaque nature of the glass prevent direct observation inside the melters. Therefore, computational models are essential for providing insight into factors that affect melter throughput. Thermocouples in the plenum provide operators with plenum temperature measurements. Operational adjustments include bubbling rate, voltage supplied to the electrodes, feed adjustments, and glass removal rate. Different computational fluid dynamics (CFD) models have been developed, each serving a specific purpose. There are CFD models of different scale melters, as well as models that capture the two-phase flow interfaces of rising bubbles in the molten glass or models with a simplified molten glass region so that the surrounding structure and plenum can be feasibly incorporated. Pilot-scale melter models have been developed to serve as validation of the methods employed in the simulation of the full-scale WTP melters. Models incorporating resolved bubbling are used to develop momentum source terms to implement into a single phase, multi-region, steady-state flow model that is being validated by measured process parameters such as glass production rate, voltage, input power, plenum temperatures, etc. The resolved bubbling model uses the multiphase volume of fluid approach to model the system with a high-resolution interface capturing scheme to maintain sharp interfaces between the molten glass and the air phase. The suite of CFD models is continually being improved to incorporate more realistic physics and achieve faster turnaround time. For example, an incremental controller is implemented to automatically adjust electrode voltage within the simulation to a molten glass set point temperature of 1150°C. Newer models feature improved meshes to ensure conformal meshes between regions and eliminate unnecessary mesh refinement in areas that are not of interest (such as boundary layers in offgas ports). Instead of explicitly modeling the structural, refractory, and insulation layers of the melter, a thermal resistance approach is used with published correlations used for boundary conditions. The development of robust and efficient CFD models will be instrumental in enabling the WTP to successfully fulfill its mission of safely stabilizing legacy nuclear waste.

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

Image processing methods for Rayleigh scattering measurements of diesel spray mixing at high repetition rate

This work describes the diagnostic implementation and image processing methods to quantitatively measure diesel spray mixing injected into a high-pressure, high-temperature environment. We used a high-repetition-rate pulse-burst laser developed in-house, a high-speed CMOS camera, and optimized the optical configuration to capture Rayleigh scattering images of the vaporized fuel jets inside a constant volume chamber. The experimental installation was modified to reduce reflections and flare levels to maximize the images’ signal-to-noise ratios by anti-reflection coatings on windows and surfaces, as well as series of optical baffles. Because of the specificities of the high-speed system, several image processing techniques had to be developed and implemented to provide quantitative fuel concentration measurements. We find that these methods involve various correction procedures such as camera linearity, laser intensity fluctuation, dynamic background flare, as well as beam-steering effects. Image inpainting was also applied to correct the Rayleigh scattering signal from large scatterers (e.g. particulates). The experiments demonstrate that applying planar laser Rayleigh scattering at high repetition rate to quantitatively resolve the mixing of fuel and ambient gases in diesel jets is challenging, but possible. The thorough analysis of the experimental uncertainty and comparisons to past data prove that such measurements can be accurate, whilst providing valuable information about the mixing processes of high-pressure diesel jets.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Physics augmented machine learning discovery of composition-dependent constitutive laws for 3D printed digital materials

Multi-material 3D printing, particularly through polymer jetting, enables the fabrication of digital materials by mixing distinct photopolymers at the micron scale within a single build to create a composite with tunable mechanical properties. Here, this work presents an integrated experimental and computational investigation into the composition-dependent mechanical behavior of 3D printed digital materials. We experimentally characterize five formulations, combining soft and rigid UV-cured polymers under uniaxial tension and torsion across three strain and twist rates. The results reveal nonlinear and rate-dependent responses that strongly depend on composition. To model this behavior, we develop a physics-augmented neural network (PANN) that combines a partially input convex neural network (pICNN) for learning the composition-dependent hyperelastic strain energy function with a quasi-linear viscoelastic (QLV) formulation for time-dependent response. The pICNN ensures convexity with respect to strain invariants while allowing non-convex dependence on composition. To enhance interpretability, we apply $L_0$ sparsification. For the time-dependent response, we introduce a multilayer perceptron (MLP) to predict viscoelastic relaxation parameters from composition. The proposed model accurately captures the nonlinear, rate-dependent behavior of 3D printed digital materials in both uniaxial tension and torsion, achieving high predictive accuracy for interpolated material compositions. This approach provides a scalable framework for automated, composition-aware constitutive model discovery for multi-material 3D printing.

Constitutive modeling↗
DOE OSTI2021

Visualization and modeling of inhibition of IL-1β and TNF-α mRNA transcription at the single-cell level

Abstract IL-1β and TNF-α are canonical immune response mediators that play key regulatory roles in a wide range of inflammatory responses to both chronic and acute conditions. Here we employ an automated microscopy platform for the analysis of messenger RNA (mRNA) expression of IL-1β and TNF-α at the single-cell level. The amount of IL-1β and TNF-α mRNA expressed in a human monocytic leukemia cell line (THP-1) is visualized and counted using single-molecule fluorescent in-situ hybridization (smFISH) following exposure of the cells to lipopolysaccharide (LPS), an outer-membrane component of Gram-negative bacteria. We show that the small molecule inhibitors MG132 (a 26S proteasome inhibitor used to block NF-κB signaling) and U0126 (a MAPK Kinase inhibitor used to block CCAAT-enhancer-binding proteins C/EBP) successfully block IL-1β and TNF-α mRNA expression. Based upon this single-cell mRNA expression data, we screened 36 different mathematical models of gene expression, and found two similar models that capture the effects by which the drugs U0126 and MG132 affect the rates at which the genes transition into highly activated states. When their parameters were informed by the action of each drug independently, both models were able to predict the effects of the combined drug treatment. From our data and models, we postulate that IL-1β is activated by both NF-κB and C/EBP, while TNF-α is predominantly activated by NF-κB. Our combined single-cell experimental and modeling efforts show the interconnection between these two genes and demonstrates how the single-cell responses, including the distribution shapes, mean expression, and kinetics of gene expression, change with inhibition.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

Primary $\gamma$-ray intensities and $\gamma$-strength functions from discrete two-step $\gamma$-ray cascades in radiative proton-capture experiments

Background: Reaction rates of radiative capture reactions can play a crucial role in the nucleosynthesis of heavy nuclei in explosive stellar environments. These reaction rates depend strongly on $\gamma$-ray decay widths in the reaction products, which are, for nonresonant capture reactions at high excitation energies, derived from the $\gamma$-ray strength function and the nuclear level density. Recently, the ratio method was applied to primary $\gamma$ rays observed from ($\textit{d, p}$) reactions and nuclear resonance fluorescence measurements to extract the dipole strength in atomic nuclei and to test the generalized Brink-Axel hypothesis. Purpose: The purpose of this work is to apply the ratio method to primary $\gamma$-ray intensities of the 63,65 Cu($\textit{p},\gamma$) reactions to extract $\gamma$-ray strength information on the nuclei 64,66 Zn. Here, the impact of spin distribution, total $\gamma$-ray decay widths, level densities, and width fluctuations on the application of the ratio method will be discussed. Additionally, by comparing the relative $\gamma$-ray strength at different excitation energies, conclusions on the validity of the generalized Brink-Axel hypothesis can be made. Method: The radiative proton capture reaction measurements have been performed at the HORUS $\gamma$-ray spectrometer of the University of Cologne at one excitation energy for each reaction. Primary $\gamma$-ray intensities have been determined by normalizing secondary $\gamma$-ray transitions in two-step cascades using their absolute branching ratio. The ratio method was applied to the measured primary $\gamma$-ray intensities as well as to previous measurements by Erlandsson et al. at different excitation energies. Results: The relative strength function curve for 64 Zn from our measurement shows no significant deviation from the previous measurement at a different excitation energy. The same is true for 66 Zn where both measurements were at almost the same excitation energy. Absolute $\gamma$-strength function values have been obtained by normalizing the relative curves to quasiparticle random phase approximation calculations because of the absence of experimental data in the respective energy region. Conclusion: The generalized Brink-Axel hypothesis, i.e., the independence of the strength function on the excitation energy, seems to hold in the studied energy region and nuclei. The method to obtain primary $\gamma$-ray intensities from two-step cascade spectra was shown to be a valuable and sensitive tool although its uncertainties are connected to the knowledge of the low-energy level scheme of the investigated nucleus. The scaling in the ratio method should be taken with care, because the relative strength is not a simple sum of $f_{E1}$ and $f_{M1}$ but a somewhat complex linear combination dependent on the excitation energy of the nucleus.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Electrocatalysis in CO 2 -Binding Organic Liquids with an Iron Porphyrin

Direct electrochemical upgrading of CO 2 in capture media is an attractive approach to carbon capture that can bypass the energy requirement for the thermal release of pure CO 2 . Here we investigate the electrocatalytic activity of iron(tetraphenylporphyrin) in the presence of organic solvents that convert into ionic liquids upon exposure to CO 2 . Four different solvent systems were tested, all of which capture CO 2 in the form of an alkyl carbonate (or carbamate) anion and an acidic ammonium cation. The electrocatalytic selectivity exhibited a strong dependence on the acidity of the capture medium, with the most basic solvent affording a high selectivity for production of CO instead of H 2 . Experimental and computational studies support a canonical mechanism in which the catalyst reacts with free CO 2 in solution, as opposed to a reaction with the alkyl carbonate that is present in high concentration. Kinetic analysis indicates that the rate-limiting step is changed from C–O protonolysis in traditional solvents to the binding of CO 2 in the capture media. Quantitative 13 C– 13 C EXSY revealed that the dissociation of the alkyl carbonate into free, solvated CO 2 is very rapid (~15 s –1 ) compared to the interconversion of HCO 3 – /CO 2 in aqueous solution. These results underscore the need to understand the mechanism and kinetics for both the release of captured CO 2 and its electrocatalytic conversion.

Carbon Dioxide↗
DOE OSTI2023

Deformation, dislocation evolution and the non-Schmid effect in body-centered-cubic single- and polycrystal tantalum

A physically-informed continuum crystal plasticity model is presented to elucidate deformation mechanisms, dislocation evolution and the non-Schmid effect in body-centered-cubic (bcc) tantalum widely used as a key structural material for mechanical and thermal extremes. We show the unified structural modeling framework informed by mesoscopic dislocation dynamics simulations is capable of capturing salient features of the large inelastic behavior of tantalum at quasi-static (10 -3 s -1 ) to extreme strain rates (5000 s -1 ) and at low (77 K) to high temperatures (873 K) at both single- and polycrystal levels. Here we also present predictive capabilities of the model for microstructural evolution in the material. To this end, we investigate the effects of dislocation interactions on slip activities, instability and the non-Schmid behavior at the single crystal level. Furthermore, ex situ measurements on crystallographic texture evolution and dislocation density growth are carried out for polycrystal tantalum specimens at increasing strains. Numerical simulation results also support that the modeling framework is capable of capturing the main features of the polycrystal behavior over a wide range of strains, strain rates and temperatures. The theoretical, experimental and numerical results at both single- and polycrystal levels provide critical insight into the underlying physical pictures for micro- and macroscopic responses and their relations in this important class of refractory bcc materials undergoing large inelastic deformations.

36 MATERIALS SCIENCE↗
DOE OSTI2025

The Co-Firing of Pine Biomass and Waste Coal in 100 and 600 MW Power Plants: A Sustainable Approach to Reduce GHG Emissions

Climate change is a global issue that has gained much attention recently. Co-firing biomass with coal/waste coal reduces the electricity sector’s GHG emissions sustainably. This study uses commercial software to model waste coal and biomass co-firing in 100 MW and 600 MW power plants. The objective is to assess the effects of fluid types (subcritical and supercritical), plant capacities (100 MW and 600 MW), boiler types (pulverized coal and circulating fluidized bed boilers), biomass and waste coal co-firing ratios (0:100, 20:80, 40:60, 60:40, 80:20, and 100:0), and carbon capture and storage efficiencies (0%, 90%, 95%, and 97%) on performance parameters such as net plant efficiency, heat rate, net plant CO2 and SO2, and particulate matter emissions. The feedstocks selected for this investigation include anthracite waste coal and loblolly pine biomass. As the biomass fraction increases from 0% to 100%, co-fired power plants net efficiency increases by 3–8%. Supercritical plants had a 6% higher net plant efficiency than the subcritical plants. The study found that the biomass’s high heating value decreased the fuel flow rate and reduced plant CO2 emissions by 10–16%. With 100% biomass power plant feed and 90% carbon capture and storage efficiency, CO2 emissions drop by 83% and SO2 and PM emissions drop to zero.

Bhoi, Prakashbhai R. (ORCID:0000000333104888)↗
NASA NTRS1993

(abstract) Experimental and Modeling Studies of the Exchange Current at the Alkali Beta'-Alumina/Porous Electrode/Alkali Metal Vapor Three Phase Boundary

The microscopic mechanism of the alkali ion-electron recombination reaction at the three phase boundary zone formed by a porous metal electrode in the alkali vapor on the surface of an alkali beta'-alumina solid electrolyte (BASE) ceramic has been studied by comparison of the expected rates for the three simplest reaction mechanisms with known temperature dependent rate data; and the physical parameters of typical porous metal electrode/BASE/alkali metal vapor reaction zones. The three simplest reactions are tunneling of electrons from the alkali coated electrode to a surface bound alkali metal ion; emission of an electron from the electrode with subsequent capture by a surface bound alkali metal ion; and thermal emission of an alkali cation from the BASE and its capture on the porous metal electrode surface where it may recombine with an electron. Only the first reaction adequately accounts for both the high observed rate and its temperature dependence. New results include crude modeling of simple, one step, three phase, solid/solid/gas electrochemical reaction.

conversion↗
DOE OSTI2023

A Theoretical Study of NH 2 Radical Reactions with Propane and Its Kinetic Implications in NH 3 -Propane Blends’ Oxidation

The reaction of NH 2 radicals with C 3 H 8 is crucial for understanding the combustion behavior of NH 3 /C 3 H 8 blends. In this study, we investigated the temperature dependence of the rate coefficients for the hydrogen abstraction reactions of C 3 H 8 by NH 2 radicals using high-level theoretical approaches. The potential energy surface was constructed at the CCSD(T)/cc-pV(T, Q)//M06-2X/aug-cc-pVTZ level of theory, and the rate coefficients were computed using conventional transition state theory, incorporating the corrections for quantum tunneling and hindered internal rotors (HIR). The computed rate coefficients showed a strong curvature in the Arrhenius behavior, capturing the experimental literature data well at low temperatures. However, at T > 1500 K, the theory severely overpredicted the experimental data. The available theoretical studies did not align with the experiment at high temperatures, and the possible reasons for this discrepancy are discussed. At 300 K, the reaction of NH 2 with C 3 H 8 predominantly occurs at the secondary C-H site, which accounts for approximately 95% of the total reaction flux. However, the hydrogen abstraction reaction at the primary C-H site becomes the dominant reaction above 1700 K. A composite kinetic model was built, which incorporated the computed rate coefficients for NH 2 + C 3 H 8 reactions. The importance of NH 2 + C 3 H 8 reactions in predicting the combustion behavior of NH 3 /C 3 H 8 blends was demonstrated by kinetic modeling.

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

Spatial dynamics of lithiation and lithium plating during high-rate operation of graphite electrodes

The principal inhibitor of fast charging lithium ion cells is the graphite negative electrode, where favorable conditions for lithium plating occur at high charge rates, causing accelerated degradation and safety concerns. The local response of graphite, both at the electrode and particle level, when exposed to fast charging conditions of around 6C is not well understood. Consequently, the conditions that lead to the onset of lithium plating, as well as the local dynamics of lithium plating and stripping, have also remained elusive. Here, we use high-speed (100 Hz) pencil-beam X-ray diffraction to repeatedly raster along the depth of a 101 μm thick graphite electrode in 3 μm steps during fast (up to 6C) charge and discharge conditions. Consecutive depth profiles from separator to current collector were each captured in 0.5 seconds, giving an unprecedented spatial and temporal description of the state of the electrode and graphite's staging dynamics during high rate conditions. The electrode is preferentially activated near the separator, and the non-uniformity increases with rate and is influenced by free-energy barriers between graphite's lithiation stages. The onset of lithium plating and stripping was quantified, occurring only within the first 15 μm from the separator. The presence of lithium plating changed the behavior of the underlying graphite, such as causing co-existence of LiC 6 and graphite in the fully discharged state. Finally, the staging behavior of graphite at different rates was quantified, revealing a high dependency on rate and drastic hysteresis between lithiation and delithiation.

25 ENERGY STORAGE↗
DOE OSTI2023

Microstructure dependent shear localization in 316L SS

Shear localization is characterized by the concentration of plastic deformation within narrow bands of material and constitutes the precursor to failure in materials subjected to high strain rates. These bands have been observed experimentally, but classical continuum material models have not been able to fully capture their development. The aim of this research is to investigate the role of microstructural differences in the dynamic shear response of wrought and additively manufactured (AM) 316L SS. Compact Forced Simple Shear specimens for Digital Image Correlation (CFSS-DIC) of wrought and additively manufactured 316L SS were tested on a Split Hopkinson Pressure Bar (SHPB) in conjunction with infrared thermometry measurements. Results indicate that the microstructures of wrought and AM 316L SS exhibit differences that affect the shear response of this materials. In conclusion, CFSS-DIC samples tested on the SHPB did not exhibit recrystallization inside the shear region, fact confirmed by thermometry measurements.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Effect of Approach Speed and Electrode Geometry on Electrostatic Discharges Off Floating Dielectrics

With the damaging effects of electrostatic discharges (ESDs), it is important to study them within various parameter regimes relevant to real-world scenarios. One such scenario studied here is a floating dielectric with no nearby ground plane with a grounded object approaching at a high rate of speed. One may encounter this scenario when moving objects by hand. Discharge current and radiated fields are captured for discharges drawn from polytetrafluoroethylene (PTFE) and poly methyl methacrylate (PMMA) with approach speeds ranging from 20 to 100 mm/s. Previous studies have shown that at lower speeds in metal-to-metal or metal-to-ground backed dielectric geometries, the peak current and discharge distance decrease with increasing speed. However, in the speed range studied here, an increase in distance and current are observed for increasing speed. Namely, the highest speeds coincide with the highest peak currents and discharge distances of approximately 800 mA and up to 24 mm. Further, with no grounded backing, as opposed to setups in other ESD studies, the electric field between the dielectric and an approaching electrode has a more uniform distribution with highest fields at the electrode, which is elucidated to be the driving factor in the differences revealed in this study. Two spheres of differing diameters and a conical electrode are used to draw discharges off the dielectric surfaces. Captured radiated fields via B-dot sensor, plasma imagery via intesified charge-coupled device (ICCD), and mapping of surface potential reveal that a majority of the energy lost during a discharge goes to gas losses and radiated fields.

36 MATERIALS SCIENCE↗
DOE OSTI2022

NUMERICAL PREDICTIONS OF MEAN PERFORMANCE AND DYNAMIC BEHAVIOR OF A 10 MWe SCO2 COMPRESSOR WITH TEST DATA VALIDATION

High-fidelity aerodynamic analysis has been demonstrated for mean performance and unsteady dynamics in a sCO2 compressor designed by Hanwha Power Systems Americas for a 10 MWe Concentrating Solar Power plant. Simulations were performed with CRUNCH CFD® software tool that was matured to accurately model near critical real fluid effects in sCO2. Pre-test predictions for mean performance were validated with test data collected later. Performance predictions were accurate and captured sensitivity of the efficiency to inlet temperature of CO2 as well as steep drop-off at high flow rates due to condensation in the inlet throat. Detailed analysis was performed to understand the source of these performance losses at near critical conditions. Unsteady dynamic effects in the compressor at off-design conditions were also identified and quantified. In particular a system wide “condensation surge” condition was detected at high flow coefficients that results in large amplitude pulsations with accompanying mass flow fluctuations at low frequencies and has potential to cause damage in closed loop systems.

hosangadi, ashvin↗
DOE OSTI2024

Canopy cover and microtopography control precipitation-enhanced thaw of ecosystem-protected permafrost

Abstract Northern high-latitudes are projected to get warmer and wetter, which will affect rates of permafrost thaw and mechanisms by which thaw occurs. To better understand the impact of rain, as well as other factors such as snow depth, canopy cover, and microtopography, we instrumented a degrading permafrost plateau in south-central Alaska with high-resolution soil temperature sensors. The site contains ecosystem-protected permafrost, which persists in unfavorable climates due to favorable ecologic conditions. Our study (2020–2022) captured three of the snowiest years and three of the four wettest years since the site was first studied in 2015. Average thaw rates along an across-site transect increased nine-fold from 6 ± 5 cm yr −1 (2015–2020) to 56 ± 12 cm yr −1 (2020–2022). This thaw was not uniform. Hummock locations, residing on topographic high points with relatively dense canopy, experienced only 8 ± 9 cm yr −1 of thaw, on average. Hollows, topographic low points with low canopy cover, and transition locations, which had canopy cover and elevation between hummocks and hollows, thawed 44 ± 6 cm yr −1 and 39 ± 13 cm yr −1 , respectively. Mechanisms of thaw differed between these locations. Hollows had high warm-season soil moisture, which increased thermal conductivity, and deep cold-season snow coverage, which insulated soil. Transition locations thawed primarily due to thermal energy transported through subsurface taliks during individual rain events. Most increases in depth to permafrost occurred below the ∼45 cm thickness seasonally frozen layer, and therefore, expanded existing site taliks. Results highlight the importance of canopy cover and microtopography in controlling soil thermal inputs, the ability of subsurface runoff from individual rain events to trigger warming and thaw, and the acceleration of thaw caused by consecutive wet and snowy years. As northern high-latitudes become warmer and wetter, and weather events become more extreme, the importance of these controls on soil warming and thaw is likely to increase.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Insights into mixing of non-isothermal multi-polymer melts for complex plastics recycling

Catalytic recycling or upcycling of plastics is often limited not by catalyst performance, but by transport, arising from highly viscous, non-Newtonian polymer melts. In this work, we develop a reactor-scale framework that integrates rheological measurements, constitutive modeling, computational fluid dynamics (CFD), and experiments to quantify mixing, heat transfer, and dispersion in surrogate hydrocarbon melts representing mixed plastics systems. Temperature- and shear rate-dependent viscosity of low-density polyethylene (LDPE) and high-density polyethylene (HDPE) is measured to create two surrogate polymers (PLD and PHD) that capture the dominant shear-thinning flow behavior while neglecting strong elastic effects, enabling tractable simulation of non-isothermal, polymer-melt mixing using a Carreau-Arrhenius generalized Newtonian framework. Three-dimensional CFD simulations are employed to evaluate impeller performance in PLD using mixing time, cavern volume, thermal uniformity, and interfacial area for regimes in which viscoelastic effects are not dominant. We show that magnetic stir bars commonly used in lab-scale studies produce large thermal gradients (~60 °C) and poor mixing, even under idealized power delivery and polymer flow conditions. In contrast, close-clearance anchor impellers achieve near-isothermal operation, reduce mixing times by up to 5×, and provide >90% active circulation volume. We further demonstrate that, at low pseudo-Deborah number (De*), motor power requirements can be predicted directly from shear rate-dependent rheology using the Carreau-Arrhenius framework, enabling rational selection of operating conditions. Extension to surrogate immiscible multi-polymer systems based on PLD and PHD shows that interfacial area is highly sensitive to operating conditions and impeller design, with coaxial anchor-turbine configurations enhancing dispersion by up to 4 × .

Close-clearance impellers↗
NASA NTRS2022

Flash Drought Onset and Development Mechanisms Captured With Soil Moisture and Vegetation Data Assimilation

Flash droughts evolve and intensify rapidly under the influence of anomalous atmospheric conditions. In this study, we investigate the role of assimilating remotely sensed soil moisture (SM) and vegetation properties in capturing the evolution and impacts of two flash droughts in the Northern Great Plains. We find that during 2016 drought triggered by anomalously high temperatures and excessive evaporative demands, multivariate data assimilation (DA) of MODIS-derived leaf area index (LAI) and Soil Moisture Active Passive SM within Noah-Multiparameterization model helps capture elevated transpiration at onset. Assimilation of LAI particularly helped model the resulting rapid decline in SM during onset with as high as 10.0% steeper rate of decline compared to the simulation without any assimilation. Modeled-SM anomalies exhibit a 7.5% and 11.7% increase in similarity with Evaporative Stress Index (ESI) data and U.S. Drought Monitor (USDM) maps, respectively. In contrast, during 2017 flash drought driven by record-low precipitation during summers, SM assimilation resulted in largest rates of decline in rootzone SM, as large as 48.4% compared to results from no assimilation. Multivariate DA of SM and LAI results in 6.7% and 14.3% higher spatial similarity with ESI and USDM, respectively, and is necessary to model rapid intensification caused by anomalous precipitation deficits. This study elucidates the need to incorporate multiple observational constraints from remote sensing to effectively capture rapid onset rates, intensification, and severity of flash drought following different propagation mechanisms. This is fundamental for drought early detection to provide a wider window of response and implement efficient mitigation strategies.

Shahryar K. Ahmad↗