Search NASA⌕ Search

SEARCH · Search NASA

Results for “Reaction diffusion systems”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
NASA NTRS1997

Enhancement of Mixing and Reaction in High-speed Combustor Flowfields

Research has been underway for a number of years, both in the United States and abroad, to develop advanced aerospace propulsion systems for use late in this century and beyond. One program is now underway at the NASA Langley Research Research has been underway for a number of years, both in the United States and abroad, to develop advanced aerospace propulsion systems for use late in this century and beyond. One program is now underway at the NASA Langley Research Center to develop a hydrogen-fueled supersonic combustion ramjet (scramjet) that is capable of propelling a vehicle at hypersonic speeds in the atmosphere. A part of that research has been directed toward the optimization of the scramjet combustor and, in particular, the efficiency of fuel-air mixing and reaction taking place in the engine. In the high-speed vehicle configurations currently being considered, achieving a high combustor efficiency becomes particularly difficult. With increasing combustor Mach number, the degree of fuel-air mixing that can be achieved through natural convective and diffusive processes is reduced leading to an overall decrease in combustion efficiency and thrust. Compressible shear layers and jets provide a good model for studying the physical processes occurring in high speed mixing and combustion in a scramjet. Mixing layers are characterized by large scale eddies that form due to the high shear that is present between the fuel and the air streams. These eddies entrain the fuel and air into the mixing region between the fluids leading to increased surface area and locally steep concentration gradients. Molecular diffusion then occurs across the strained interfaces. In an early study of high-speed mixing, Brown and Roshko show that the spreading rate of a supersonic mixing layer decreases with increasing Mach number, exhibiting a factor of three decrease in spread rate as compared with an incompressible mixing layer with the same density ratio. They conclude that the reduced spread rate is primarily due to compressibility. Papamoschou and Roshko also observe that the spreading rate of compressible mixing layers is significantly reduced over that of incompressible layers. To characterize the structure of the flow quantitatively, they define a convective Mach number. The reduction in mixing layer spreading rate is shown in these experiments to correlate well with increasing convective Mach number.

Drummond, J. Philip↗
NASA NTRSDate not supplied

From the Interstellar Medium to Ocean Worlds: new challenges for Laboratory Astrophysics

From the Interstellar Medium to Ocean Worlds: new challenges for Laboratory Astrophysics For the past 35+ years, laboratory experiments in Astrophysics have been focusing on identifying and understanding the molecular species and chemistry occurring in interstellar environments, such as the diffuse interstellar medium (ISM) and dense molecular clouds. A significant portion of these laboratory efforts concerns investigations of carbon-based molecules ranging in size from a single carbon atom (e.g., CO, CO2, H2CO) to large polycyclic aromatic hydrocarbon (PAH) molecules potentially consisting of more than one hundred carbon atoms as well as their reactions with H2O and other interstellar species. The laboratory facilities built to conduct these investigations are capable of monitoring the chemistry and measuring the spectra under a wide range of interstellar conditions (e.g., low temperatures and pressures). The past 15 years of space exploration revealed the presence of at least nine other Ocean Worlds, besides Earth, in our Solar System. These worlds include moons and dwarf planets, ice-covered surfaces, subsurface oceans, as well as geysers venting H2O and other compounds, including organics in at least one environment, into space. Although these worlds represent environments that are significantly different from those traditionally considered the realm of Laboratory Astrophysics, their facilities are well situated to provide insight into the chemistry occurring in and around these worlds as well as experimental data for the interpretation of observations from missions visiting these Ocean World. This presentation will discuss the areas where Laboratory Astrophysics can provide experimental data to help understand the chemistry, and mission data of Ocean Worlds, and the challenges research into Ocean Worlds presents to our current facilities and how they can be overcome. Lastly, this presentation will introduce the new ICEE (In-situ Carbon Evolution Experiments) facility at NASA Ames Research Center as an example of how a lab dedicated to Laboratory Astrophysics can adapt itself to also meet the needs for Ocean Worlds research.

Andrew Lige Mattioda↗
DOE OSTI2024

Effect of f -element complexation on the radiolysis of 2-ethylhexylphosphonic acid mono-2-ethylhexyl ester (HEH[EHP])

A systematic study of the impact on the chemical reactivity of the oxidising n-dodecane radical cation (RH˙ + ) with f-element complexed 2-ethylhexylphosphonic acid mono-2-ethylhexyl ester (HEH[EHP]) has been undertaken utilizing time-resolved electron pulse radiolysis/transient absorption spectroscopy and high-level quantum mechanical calculations. Lanthanide ion complexed species, [Ln((HEH[EHP]) 2 ) 3 ], exhibited vastly increased reactivity (over 10× faster) in comparison to the non-complexed ligand in n-dodecane solvent, whose rate coefficient was k = (4.66 ± 0.22) × 10 9 M -1 s -1 . Similar reactivity enhancement was also observed for the corresponding americium ion complex, k = (5.58 ± 0.30) × 10 10 M -1 s -1 . The vastly increased reactivity of these f-element complexes was not due to simple increased diffusion-control of these reactions; rather, enhanced hole transfer mechanisms for the complexes were calculated to become energetically more favourable. Interestingly, the observed reactivity trend with lanthanide ion size was not linear; instead, the rate coefficients showed an initial increase (Lu to Yb) followed by a decrease (Tm to Ho), followed by another increase (Dy to La). This behaviour was excellently predicted by the calculated reaction volumes of these complexes. In conclusion, complementary cobalt-60 gamma irradiations for select lanthanide complexes demonstrated that the measured kinetic differences translated to increased ligand degradation at steady-state timescales, affording ~38% increase in ligand loss of a 1:1 [La((HEH[EHP]) 2 ) 3 ]: HEH[EHP] ratio system.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

DNS of Flame Stabilization Dynamics of a Swirl-Stabilized Spray Burner Using Sustainable Aviation Fuels

Concern for emission reduction has motivated the development of new cost-effective alternative sustainable aviation fuels (SAFs). Drop-in SAFs with blends of the certified and the alternative fuels are beneficial because they do not require engine modifications for use in current aviation engines. Combustion characteristics are of particular concern when comparing a new fuel with Jet-A. Lean blow-off and flame stabilization have been identified to be governed by the fuel propensity to autoignition, i.e. derived cetane number. However, fundamental investigation is required to identify the impact of derived cetane number on the combustion mode of flame stabilization in a realistic combustor. The combustion mode has also a direct impact on turbulent combustion modeling closures. DNS simulations are performed in the low-Mach solver of the Pele Suite called PeleLMeX. Lagrangian multi-phase modeling is used to capture the liquid spray injection of Jet-A (reference fuel) and C1 as a representative of a low cetane number SAF. Adaptive Mesh Refinement (AMR) is used to enable a more efficient simulation of a more realistic domain size and embedded boundary treatment is used to model a bluff-body geometry. Local extinction and edge flame propagation were observed for both fuels. The edge flame propagation mode was quantified in terms of a Damkohler number defined as the ratio between progress variable reaction rate and its diffusive flux. Initial analysis suggests that a mixed mode combustion occurs for the edge flame propagation, with flame propagation assisted by ignition and autoignition co-existing for both Jet-A and C1 flames. The analysis of the extinction region shows a larger progress variable for Jet-A which can further increase the local displacement speed and present a faster reignition of the stoichiometric mixture fraction.

adaptive mesh refinement↗
DOE OSTI2025

Carbon mineralization pathways in interfacial adsorbed water nanofilms

Carbon mineralization in humidified carbon dioxide offers a promising route to mitigate anthropogenic emissions in a world stressed by water security. Despite its technological importance, our understanding of carbonation in water-poor environments lags, as traditional dissolution-precipitation pathways struggle to explain the adsorbed water nanofilm-mediated reactivity. Here, we utilize in operando X-ray diffraction (XRD) and advanced molecular simulations to investigate nanoconfined reactions driving forsterite carbonation, the magnesium-rich olivine. By examining magnesium ion dissolution and transport in atomistic simulations of the forsterite-water-carbon dioxide interface and comparing these with the in operando XRD activation energies, we identify both processes as rate-limiting at saturation. Our simulations reveal a mechanistic view of interfacial carbonation, where dissolution and precipitation are mediated by anomalous quasi two-dimensional diffusion. The transport process involves intermittent diffusive hopping in the desorbed state, separated by crawling events that are spatially short but temporally long. This understanding transcends carbon mineralization, with implications for understanding the transport of contaminants in geosystems, the design of multifunctional materials, water desalination, and molecular recognition systems.

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

Interconnected cathode-electrolyte double-layer enabling continuous Li-ion conduction throughout solid-state Li-S battery

All-solid-state lithium (Li) batteries with high energy density are a promising solution for the next-generation energy storage systems in large-scale devices. To simultaneously overcome the challenges of poor ionic conduction of solid electrolytes and shuttling of active materials, here we introduce a functional electrolyte-cathode bilayer framework with interconnected LLAZO channels from the electrolyte into the cathode for advanced solid-state Li-S batteries. Differing from the traditional solid-state batteries with separated layer compositions, the introduced bilayer framework provides ultrafast and continuous ion/electron conduction. Instead of transferring Li+ across the polymer and garnet phases which involve huge interfacial resistance, Li+ is directly conducted through the LLAZO channels created continuously from the cathode layer to the solid electrolyte layer, significantly shortening the diffusion distance and facilitating the redox reaction of sulfur and sulfides. A stable cycle life is demonstrated in the prototype Li-S solid-state batteries assembled with the introduced LLAZO-LLAZO@CNF interconnected bilayer framework. High capacity is obtained at room temperature, indicating the superior electrochemical properties of the bilayer framework that result from the unique design of the interconnected LLAZO garnet phase.

25 ENERGY STORAGE↗
NASA NTRS2001

Frontal Polymerization in Microgravity: Bubble Behavior and Convection on the KC-135 Aircraft

Frontal polymerization is a mode of converting monomer into polymer via a localized exothermic reaction zone that propagates through the coupling of thermal diffusion and Arrhenius reaction kinetics. Frontal polymerization was discovered in Russia by Chechilo and Enikolopyan in 1972. The macrokinetics and dynamics of frontal polymerization have been examined in detail and applications for materials synthesis considered. Large temperature and concentration gradients that occur in the front lead to large density gradients. A schematic is presented for a liquid monomer, usually a monoacrylate, being converted to a liquid (thermoplastic) polymer. The velocity can be controlled by the initiator concentration but is on the order of a cm/min. If the liquid monomer is multifunctional, then a solid (thermoset) polymer is formed. Convection can occur with all types of monomers if the front propagates up a tube. Bowden et al. studied liquid/solid systems. McCaughey et al. studied liquid polymer systems. Descending fronts in thermoplastic systems are also susceptible to the Rayleigh-Taylor instability.

Pojman, John A.↗
DOE OSTI2022

Developing novel electrodes with ultralow catalyst loading for high-efficiency hydrogen production in proton exchange membrane electrolyzer cells

Hydrogen plays more crucial roles for decarbonizing the planets and meeting the climate challenges because of its high energy density and zero-emission. It can be produced with proton exchange membrane electrolyzer cells (PEMECs) driven by sustainable and renewable energy resources. Although PEMECs have a number of advantages, including high purity production, quick response, and the ability to operate at high pressure facilitating the gas delivering, their performance and cost greatly hinder their commercial-scale applications. To achieve high-efficiency and cost-reduced hydrogen production in PEMECs, we proposed thin engineered liquid/gas diffusion layers (LGDLs) and associated electrodes, i.e., catalyst-coated LGDLs (CCLGDLs), over conventional porous transport layers (PTLs) and catalyst-coated membranes (CCMs). The research approaches in this project are based on material synthesis, in-situ and ex-situ characterizations, component design and treatment, numerical modeling, and cost analysis. The thin and tunable LGDLs (TT-LGDLs) and CCLGDLs were successfully developed with great performance improvement as demonstrated in lab-scale, bench-scale, and system-scale electrolyzer tests. The electrode thickness was reduced from 370 µm to less than 100 µm with simplified fabrication processes. With the catalytically enhanced Ir-based catalyst coating, the as-developed CCLGDLs with a catalyst loading of 0.34 mg Ir /cm 2 achieved a cell performance of 1.77 V at 2 A cm -2 , exhibiting the catalyst mass activity enhanced by >20 times with significant catalyst saving over conventional catalyst cell design. In-situ PEMEC characterizations, including the current distribution mapping and high-speed and multiscale visualizations, were conducted for a deeper understanding of mass transport and electrochemical reactions within an electrolyzer with LGDLs and CCLGDLs. A 2D cell model was developed and validated for the enhanced performance on TT-LGDL through reducing ohmic losses due to nonuniform hydration and water transport. Further, the cost analysis results have shown a path to move beyond equivalency and surpass costs associated with the project baseline. In this project, the design and fabrication of TT-LGDLs and CCLGDLs will contribute to the performance enhancement, manufacturing simplification, and cost reduction for PEMECs and other energy conversion devices, thus shortening their pathways towards commercialization. This project also provides a good foundation for furthering the in-situ reaction interface research.

08 HYDROGEN↗
DOE OSTI2025

Statistical distributions for transient transport

Here, this paper introduces the use of statistical distributions based on transport differential equations for clear distinction of transport modes within transient kinetic experiments. More specifically, novel techniques are developed for the transient data obtained through the Temporal Analysis of Products (TAP) reactor and are applicable to experiments where pulse response into a gas flow is used. The methodology allows distinguishing between two domains of diffusion transport in heterogeneous catalytic systems, i.e., Knudsen and non-Knudsen diffusion, using statistical fingerprints, and finding the transition domain. Two distribution parameters were obtained that directly result in coefficients that correspond to the concentration and the rate of transport. Using a linear relationship between the rate and concentration coefficients, Knudsen diffusion is revealed when the rate of transport is constant and non-Knudsen diffusion is confirmed when the rate of transport coefficient is a function of the concentration coefficient. As a result, accurate transport information can be extracted from experimental data even in the presence of comprising instrument drift or noise particularly when analyzing higher pressure pulse responses with complex transport. This enables more direct investigation of experiments influenced by gas-phase reactions.

TAP reactor↗
DOE OSTI2026

Thermodynamic Control of Interface Directs MnO 2 Nucleation Chemistry for Dense and Conformal Electrodeposition

Manganese dioxide (MnO 2 ) is widely recognized as a promising material for high-energy-density energy storage systems due to its broad applicability and facile electrodeposition. However, achieving uniform, thin, and high-mass-loading MnO 2 coatings on high-surface-area electrodes remains a significant challenge. Conventional electrodeposition methods typically yield nonuniform, thick layers with poor conductivity and limited material utilization, restricting their practical use. Here, we uncover a thermodynamically engineered vanadyl/pervanadyl (VO 2+ /VO 2 + ) interface that fundamentally reshapes MnO 2 electrodeposition chemistry, enabling highly uniform and dense coatings. Here, combining in situ AFM measurement, Classical Nucleation Theory, and Johnson–Mehl–Avrami–Kolmogorov modeling, we show that this interface reduces early-stage detectable MnO 2 island size by 35-fold and shifts the MnO 2 growth from diffusion-limited to reaction-limited progressive nucleation. This thermodynamically controlled interface yields highly dense and conformal MnO 2 films with record-high mass loading of 241 mg cm –2 (1607 mg cm –3 ) on 3D-printed graphene aerogels, without compromising porosity or inducing thickness gradient. As a prototype demonstration, the resulting MnO 2 electrodes deliver record-setting volumetric performance in both capacitors (106 F cm –3 ) and Zn//MnO 2 pouch cells (162 mAh cm –3 ). Beyond energy storage, our findings demonstrate the significance of thermodynamic interface control in MnO 2 nucleation chemistry for achieving dense and uniform coatings on various substrates, with implications for electrocatalysis, semiconductor processing, and advanced materials manufacturing.

Batteries↗
DOE OSTI2022

The growth and nanothermite reaction of 2Al/3NiO multilayer thin films

Nanothermite NiO–Al is a promising material system for low gas emission heat sources; still, its reactive properties are highly dependent on material processing conditions. In the current study, sputter deposition is used to fabricate highly controlled nanolaminates comprised of alternating NiO and Al layers. Films having an overall stoichiometry of 2Al to 3NiO were produced with different bilayer thicknesses to investigate how ignition and self-sustained, high temperature reactions vary with changes to nanometer-scale periodicity and preheat conditions. Ignition studies were carried out with both hot plate and laser irradiation and compared to slow heating studies in hot-stage x-ray diffraction. Ignition behavior has bilayer thickness and heating rate dependencies. The 2Al/3NiO with λ ≤ 300 nm ignited via solid/solid diffusion mixing (activation energy, E a = 49 ± 3 kJ/mole). Multilayers having λ ≥ 500 nm required a more favorable mixing kinetics of solid/liquid dissolution into molten Al (E a = 30 ± 4 kJ/mole). This solid/liquid dissolution E a is a factor of 5 lower than that of the previously reported powder compacts due to the elimination of a passivating Al oxide layer present on the powder. The reactant mixing mechanism between 300 and 500 nm bilayer thicknesses was dependent on the ignition source's heating rate. The self-propagating reaction velocities of 2Al/3NiO multilayers varied from 0.4 to 2.5 m/s. Pre-heating nanolaminates to temperatures below the onset reaction temperatures associated with forming intermediate nickel aluminides at multilayer interfaces led to increased propagation velocities, whereas pre-heating samples above the onset temperatures inhibited subsequent attempts at laser ignition.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Local and Global Sensitivity Analysis of a Reactive Transport Model Simulating Floodplain Redox Cycling

Reactive transport models (RTMs) are essential tools that simulate the coupling of advective, diffusive, and reactive processes in the subsurface, but their complexity makes them difficult to understand, develop and improve without accompanying statistical analyses. Although global sensitivity analysis (SA) can address these issues, the computational cost associated with most global SA techniques limits their use with RTMs. In this study, we apply distance-based generalized sensitivity analysis (DGSA), a novel and computationally efficient method of global SA, to a floodplain-scale RTM and compare DGSA results to those from local SA. Our test case focuses on the impact of 17 uncertain environmental parameters on spatially and temporally variable redox conditions within a floodplain aquifer. The input parameters considered include flow and diffusion rates, geochemical reaction rates, and the spatial distribution of sediment facies. Sensitivity was evaluated for three distinct components of the model response, encompassing both multidimensional and categorical output. Parameter rankings differ between local SA and DGSA, due to nonlinear effects of individual parameters and interaction effects between parameters. DGSA results show that fluid residence time, which is controlled by aquifer permeability, generally exerts a stronger control on redox conditions than do geochemical reaction rates. Sensitivity indices also demonstrate that sulfate reduction is key for establishing and maintaining reducing conditions throughout the aquifer. Furthermore, these results provide insights into the key drivers of heterogeneous redox processes within floodplain aquifers, as well as the main sources of uncertainty when modeling complex subsurface systems.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2001

The sugar model: catalytic flow reactor dynamics of pyruvaldehyde synthesis from triose catalyzed by poly-l-lysine contained in a dialyzer

The formation of pyruvaldehyde from triose sugars was catalyzed by poly-l-lysine contained in a small dialyzer with a 100 molecular weight cut off (100 MWCO) suspended in a much larger triose substrate reservoir at pH 5.5 and 40 degrees C. The polylysine confined in the dialyzer functioned as a catalytic flow reactor that constantly brought in triose from the substrate reservoir by diffusion to offset the drop in triose concentration within the reactor caused by its conversion to pyruvaldehyde. The catalytic polylysine solution (400 mM, 0.35 mL) within the dialyzer generated pyruvaldehyde with a synthetic intensity (rate/volume) that was 3400 times greater than that of the triose substrate solution (12 mM, 120 mL) outside the dialyzer. Under the given conditions the final yield of pyruvaldehyde was greater than twice the weight of the polylysine catalyst. During the reaction the polylysine catalyst was poisoned presumably by reaction of its amino groups with aldehyde reactants and products. Similar results were obtained using a dialyzer with a 500 MWCO. The dialyzer method of catalyst containment was selected because it provides a simple and easily manipulated experimental system for studying the dynamics and evolutionary development of confined autocatalytic processes related to the origin of life under anaerobic conditions.

Biogenesis↗
DOE OSTI2024

Degradation of Poly- and Perfluoroalkyl Substances (PFAS) in Water via High Power, Energy-Efficient Electron Beam Accelerator

The goal of the 2-year workplan was to see if electron beam (EB) could be used to break down a sub-set of the larger chemical family of per and polyfluoroalkylated substances (PFAS) in an energy efficient and economical manner when compared to conventional water treatment technologies. Year one (Y1) work focused on sample EB treatment work in the Fermi National Accelerator Laboratory’s (FNALs) Accelerator Applications Demonstration and Development (A2D2) EB accelerator. While there are reportedly thousands of types of PFAS, for the point of most of the work herein, a small subset was examined, typically perfluorooctane sulfonate (PFOS) and perfluorooctanoate (PFOA). PFOA and PFOS are two of the most well studied PFAS and are studied for baseline evaluations and are considered most useful. The work from Y1 provided information about the optimal operating parameters and additives to use when treating PFOS and PFOA via EB. The data were then used to see where in a water treatment system an EB accelerator would be best suited to treat PFAS. A conventional water treatment technology, GAC, was then compared to e-beam treatment technology with respect to energy and costs for treatment. In year two (Y2), several conventional e-beam accelerator designs, and FNAL’s developmental compact SRF accelerator design, were evaluated for their suitability in PFAS treatment, from an energy efficiency and cost standpoint. Several EB parameters were evaluated and optimized for the removal of PFOA and PFOS from water at normal pressure and temperature, measured as total PFAS removal. Under the optimized test conditions both PFOA showed complete destruction to inorganic fluoride, and PFOS to inorganic fluoride and sulfate, with mass balance. The effect on PFAS removal relative to solution pH, total EB dose, EB dose rate, dissolved oxygen concentration (DO), temperature, and initial PFAS concentration were evaluated. In general, PFOA was easier to destroy than PFOS. Degradation products, typically observed under less-than-optimal EB conditions, provided insight to degradation mechanisms. Products were identified to rule out possible deleterious biproduct formation. The water radiolysis radical reaction kinetics with PFOS and PFOA were not dependent on the initial concentration over 5-orders of magnitude from 2 μg/L to 20 mg/L. This is thought to be because there was an overabundance of the reactive water radiolysis radicals relative to PFAS molecules and largely attributed to aqueous electrons. The reaction rates appeared to be diffusion limited. Testing at higher concentrations (100-200 mg/L) showed a decrease in removal efficiency, suggesting alternative kinetics, possibly second order rates, at higher concentrations. In all, we successfully defined a set of optimal EB parameters to treat PFOA and PFOS at concentrations of 20 mg/L in water with destruction efficiencies near 100%. We further tested the optimized EB parameters with other types of PFAS, including shorter and longer fluorocarbon chain homologs of PFOA and PFOS, and PFAS with alternative functional groups such as sulfonamides. Based on our results EB can be optimized as an effective destructive technology for removing PFAS from water. The conditions optimized for PFOA and PFOS were less effective with ultra-short fluorocarbon compounds like TFMS, PFES, PFPS and PFBS, and likely require re-optimization of parameters to them. In all, it was determined that from a cost and energy efficiency standpoint, EB would be best applied to waste streams with relatively high concentrations of PFOS and PFOA and is not as cost effective as GAC treatment for removing low concentrations of PFAS from water. Higher concentrations of PFAS can be found in the wastewater of conventional treatment processes such as RO and IE and therefore EB may be used to supplement such treatment technologies. Some real-world IE regeneration wash water and RO reject water containing higher concentrations of PFAS and obtained from pilot scale industrial wastewater treatment system at a fluorochemical manufacturing facility, showed that EB could remove PFAS from such types of wastewaters. The IE regenerant wash water appeared to be the most efficient of the two types of wastewaters tested. However, some further optimization of the EB parameters for the specific PFAS types present in those wastewaters may be required. Also, the effects of co-present TOC and mineral salts should be considered during such optimization efforts. From the experimental Y1 results it was seen that the aqueous electron drives degradation of the PFAS. In a hypothetical water treatment skid using EB for PFAS destruction the parameters of the system should be optimized to promote aqueous electron production. Before EB treatment, the PFAS should be preconcentrated when possible, the pH should be raised to pH 10 or higher to enhance aqueous electron production, and the water should be nitrogen purged to remove dissolved oxygen to minimize aqueous electron scavenging. An excel spreadsheet was created that calculates optimal conditions based on inlet PFAS concentration and desired outlet concentration, by optimizing the accelerator power, dose rate, water treatment rate, pH and dissolved oxygen levels to reach the desired endpoint. Given this information on accelerator operating conditions five different EB accelerator systems were compared. One EB system was a continuous-wave, linear superconducting accelerator being designed at Fermilab. Three other EB systems (IMPELA at 5% and 25% duty factor and the ILU-14) were normal conducting pulsed linear accelerators. The fifth system was an IBA Rhodotron which is a normal conducting, circular, continuous-wave accelerator. The accelerator efficiency (% of the incoming power that is used in water treatment) was the dominating factor in accelerator choice. The radio frequency (RF) power supply and the accelerator design (superconducting versus warm technology) drive the accelerator efficiency. The IBA Rhodotron was seen to be the most energy efficient commercially available technology with a wall-plug (total) power efficiency of 43% at 400 kW. The Fermilab design, with a prototype for a different application currently being fabricated, was the most energy efficient at 55% when driven by a Klystron RF power supply and as high as 77% when powered by a magnetron. As the Fermilab design was the most energy efficient by approximately 10-30%, further design work was done on the accelerator and beam delivery system specific to the destruction of PFAS in water. The Fermilab design is unique from industrial accelerators in that is superconducting. Superconducting technology allows for the acceleration of electrons without losses. The accelerator must be cooled to below the point where it is superconducting and is operated around 4 degrees Kelvin. The bulk of the design work for the accelerator is on making the accelerator as energy efficient as possible so that it does not require liquid helium and can be cooled with conduction cooling via cryocoolers. Final design work resulted in an EB accelerator that would operate at minimally 200 kW and 10 MeV. Prototype construction would cost $\$ $7.8 million dollars when driven by a Klystron power supply. A second version of the same accelerator would cost $\$ $5.5 million dollars when driven by a magnetron that is still under development. The commercially available 300 kW IBA Rhodotron cost was estimated at approximately $\$ $9 million. While it is hard to directly compare, an operational GAC system used by 3M for groundwater treatment capital cost (2022 dollars) was estimated to cost $\$ $3.3 million. While the capital expense of the EB accelerator systems was higher than GAC, the accelerator EB treatment would result in destruction of the PFAS and not just sequestration of PFAS to form a new waste stream that requires further treatment or disposal. The operating cost to destroy the PFAS via 400 kw EB system was less than $\$ $1000/kg of PFAS destroyed when treating at a 20 mg/L PFAS concentration, compared to GAC with operating costs that calculated at $\$ $27,530 per kg of PFAS sequestered when treating 100 μg/L PFOA and PFOS combined concentration.

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

Mechanisms of mechanochemical synthesis of cesium lead halides: pathways toward stabilization of α-CsPbI 3

Cesium lead iodide with cubic perovskite structure (α-CsPbI 3 ) is gaining significant interest in photovoltaic applications due to its excellent absorbance of the visible solar light and other attractive optoelectronic properties. However, the synthesis of stable α-CsPbI 3 poses a significant challenge. Mechanochemical synthesis is emerging as a suitable method for the preparation of cesium lead halides. This work investigates the ball milling-induced synthesis of cesium lead halides perovskite phase using halide mixing or doping approaches. The synthesis in the CsI + PbI 2 , CsBr + PbBr 2 , CsBr + PbI 2 , and CsI + PbI 2 + NdI 3 mixtures and halide exchange reactions in the CsPbBr 3 + 3KI and CsBr + PbBr 2 + 3KI systems are investigated to elucidate the mechanism of this process. Then, CsPb(I 1–x Br x ) 3 and CsPb( 1–y )Nd y I 3 materials with different x and y ratios are prepared, and their stability is probed in the air using light absorption spectroscopy. These results suggest that Nd doping is more efficient in the stabilization of the perovskite structure than partial replacement of iodine with bromine. Microstructure observations reveal the existence of two different product formation mechanisms depending on the mechanical properties of reactants. The results reveal that the milling temperature has a significant impact on the reaction kinetics. Here, the produced particles nucleate and grow at the reactant interface and retard the synthesis reaction by creating a diffusion barrier. Extended milling reduces the product particle size and creates fresh contact between reactants, thus facilitating reaction completion.

36 MATERIALS SCIENCE↗
NASA NTRS1987

Potential surfaces for O atom-polymer reactions

Ab initio quantum chemistry methods are used to study the energetics of interactions of O atoms with organic compounds. Polyethylene (CH2)n has been chosen as the model system to study the interactions of O(3P) and O(1D) atoms with polymers. In particular, H abstraction is investigated and polyethylene is represented by a C3 (propane) oligomeric model. The gradient method, as implemented in the GRADSCF package of programs, is used to determine the geometries and energies of products and reactants. The saddle point, barrier geometry is determined by minimizing the squares of the gradients of the potential with respect to the internal coordinates. To correctly describe the change in bonding during the reaction at least a two configuration MCSCF (multiconfiguration self consistent field) or GVB (generalized valence bond) wave function has to be used. Basis sets include standard Pople and Dunning sets, however, increased with polarization functions and diffuse p functions on both the C and O atoms. The latter is important due to the O(-) character of the wave function at the saddle point and products. Normal modes and vibrational energy levels are given for the reactants, saddle points and products. Finally, quantitative energetics are obtained by implementing a small CAS (complete active space) approach followed by limited configuration interaction (CI) calculations. Comparisons are made with available experimental data.

Laskowski, B. C.↗
NASA NTRS1997

Correspondence between discrete and continuous models of excitable media: trigger waves

We present a theoretical framework for relating continuous partial differential equation (PDE) models of excitable media to discrete cellular automata (CA) models on a randomized lattice. These relations establish a quantitative link between the CA model and the specific physical system under study. We derive expressions for the CA model's plane wave speed, critical curvature, and effective diffusion constant in terms of the model's internal parameters (the interaction radius, excitation threshold, and time step). We then equate these expressions to the corresponding quantities obtained from solution of the PDEs (for a fixed excitability). This yields a set of coupled equations with a unique solution for the required CA parameter values. Here we restrict our analysis to "trigger" wave solutions obtained in the limiting case of a two-dimensional excitable medium with no recovery processes. We tested the correspondence between our CA model and two PDE models (the FitzHugh-Nagumo medium and a medium with a "sawtooth" nonlinear reaction source) and found good agreement with the numerical solutions of the PDEs. Our results suggest that the behavior of trigger waves is actually controlled by a small number of parameters.

Non-NASA Center↗
DOE OSTI2023

Using Homogenized Models to Explore the Effect of Fracture Densities on Weathering

Despite its importance, only a few researchers have incorporated the effects of fracturing into models of reactive transport for rock weathering. Here we explore 2D simulations that describe weathering under conditions of diffusive and advective transport within heterogeneous media consisting of rocky blocks and fractures. In our simulations, the Darcy velocities vary in space and time and depend on weathering processes within the rock matrix. We explore simulations with saturated and unsaturated flow for weathering bedrock that consists of blocks separated by inert or weathered material. The simulations show that a simplified homogenized model can approximate exact solutions for some of the simulated columns and hills and can allow exploration of coupling between flow and reaction in fractured rock. These hillslope simulations document that, even in the presence of 2D water flow, i) an increase in fracture density results in faster weathering advance rates; and ii) the water table locates deeper for a rock system that is weathered and fractured rather than weathered and unfractured. Some of these patterns have also been observed for natural systems. But these simulations also highlight how simplified models that do not use appropriate averaging of heterogeneities can be inaccurate in predicting weathering rate for natural systems. For example, if water flows both vertically and laterally through the vadose zone of a hill, then a prediction of the depth of regolith that is based on modeling strictly unidirectional downward infiltration will be unrealistically large. Likewise, if the fracture density observed near the land surface is used in a model to predict depth of weathering for a system where the fracture density decreases downward, the model will overestimate regolith depth. Learning how to develop accurately homogenized models could thus enable better conceptual models and predictions of weathering advance in natural systems.

58 GEOSCIENCES↗