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 505 records · Page 28
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Reactive particle-tracking solutions to a benchmark problem on heavy metal cycling in lake sediments

Geochemical systems are known to exhibit highly variable spatiotemporal behavior. This may be observed both in non-smooth concentration curves in space for a single sampling time and also in variability between samples taken from the same location at different times. However, most models that are designed to simulate these systems provide only single-solution smooth curves and fail to capture the noise and variability seen in the data. We apply a recently developed reactive particle-tracking method to a system that displays highly complex geochemical behavior. When the method is made to most closely resemble a corresponding Eulerian method, in its unperturbed form, we see near-exact match between solutions of the two models. More importantly, we consider two approaches for perturbing the model and find that the spatially-perturbed condition is able to capture a greater degree of the variability present in the data. This method of perturbation is a task to which particle methods are uniquely suited and Eulerian models are not well-suited. Additionally, because of the nature of the algorithm, noisy spatial gradients can be highly resolved by a large number of mobile particles, and this incurs negligible computational cost, as compared to expensive chemistry calculations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Learning in Modal Space: Solving Time-Dependent Stochastic PDEs Using Physics-Informed Neural Networks

One of the open problems in scientific computing is the long-time integration of nonlinear stochastic partial differential equations (SPDEs), especially with arbitrary initial data. We address this problem by taking advantage of recent advances in scientific machine learning and the spectral dynamically orthogonal (DO) and borthogonal (BO) methods for representing stochastic processes. The recently introduced DO/BO methods reduce the SPDE to solving a system of deterministic PDEs and a system of stochastic ordinary differential equations. Specifically, we propose two new physics-informed neural networks (PINNs) for solving time-dependent SPDEs, namely the neural network (NN)-DO/BO methods. The proposed methods incorporate the DO/BO constraints into the loss function (along with the modal decomposition of the SPDE) with an implicit form instead of generating explicit expressions for the temporal derivatives of the DO/BO modes. Hence, the NN-DO/BO methods can overcome some of the drawbacks of the original DO/BO methods. For example, we do not need the assumption that the covariance matrix of the random coefficients is invertible as in the original DO method, and we can remove the assumption of no eigenvalue crossing as in the original BO method. Moreover, the NN-DO/BO methods can be used to solve time-dependent stochastic inverse problems with the same formulation and same computational complexity as for forward problems. Furthermore, we demonstrate the capability of the proposed methods via several numerical examples, namely: (1) A linear stochastic advection equation with deterministic initial condition: we obtain good results with the proposed methods, while the original DO/BO methods cannot be applied directly in this case. (2) Long-time integration of the stochastic Burgers' equation: we show the good performance of NN-DO/BO methods, especially the effectiveness of the NN-BO approach for such problems with many eigenvalue crossings during the whole time evolution, while the original BO method fails. (3) Nonlinear reaction diffusion equation: we consider both the forward problem and the inverse problems, including very noisy initial point values, to investigate the flexibility of the NN-DO/BO methods in handling inverse and mixed type problems. Taken together, these simulation results demonstrate that the NN-DO/BO methods can be employed to effectively quantify uncertainty propagation in a wide range of physical problems, but future work should address the efficiency issue of PINNs for forward problems.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1996

Combustion Processes in Hybrid Rocket Engines

In recent years, there has been a resurgence of interest in the development of hybrid rocket engines for advanced launch vehicle applications. Hybrid propulsion systems use a solid fuel such as hydroxyl-terminated polybutadiene (HTPB) along with a gaseous/liquid oxidizer. The performance of hybrid combustors depends on the convective and radiative heat fluxes to the fuel surface, the rate of pyrolysis in the solid phase, and the turbulent combustion processes in the gaseous phases. These processes in combination specify the regression rates of the fuel surface and thereby the utilization efficiency of the fuel. In this paper, we employ computational fluid dynamics (CFD) techniques in order to gain a quantitative understanding of the physical trends in hybrid rocket combustors. The computational modeling is tailored to ongoing experiments at Penn State that employ a two dimensional slab burner configuration. The coordinated computational/experimental effort enables model validation while providing an understanding of the experimental observations. Computations to date have included the full length geometry with and with the aft nozzle section as well as shorter length domains for extensive parametric characterization. HTPB is sed as the fuel with 1,3 butadiene being taken as the gaseous product of the pyrolysis. Pure gaseous oxygen is taken as the oxidizer. The fuel regression rate is specified using an Arrhenius rate reaction, which the fuel surface temperature is given by an energy balance involving gas-phase convection and radiation as well as thermal conduction in the solid-phase. For the gas-phase combustion, a two step global reaction is used. The standard kappa - epsilon model is used for turbulence closure. Radiation is presently treated using a simple diffusion approximation which is valid for large optical path lengths, representative of radiation from soot particles. Computational results are obtained to determine the trends in the fuel burning or regression rates as a function of the head-end oxidizer mass flux, G=rho(e)U(e), and the chamber pressure. Furthermore, computation of the full slab burner configuration has also been obtained for various stages of the burn. Comparisons with available experimental data from small scale tests conducted by General Dynamics-Thiokol-Rocketdyne suggest reasonable agreement in the predicted regression rates. Future work will include: (1) a model for soot generation in the flame for more quantitative radiative transfer modelling, (2) a parametric study of combustion efficiency, and (3) transient calculations to help determine the possible mechanisms responsible for combustion instability in hybrid rocket motors.

Venkateswaran,S.↗
DOE OSTI2020

A mass-transfer particle-tracking method for simulating transport with discontinuous diffusion coefficients

The problem of a spatially discontinuous diffusion coefficient (D(x)) is one that may be encountered in hydrogeologic systems due to natural geological features or as a consequence of numerical discretization of flow properties. To date, mass-transfer particle-tracking (MTPT) methods, a family of Lagrangian methods in which diffusion is jointly simulated by random walk and diffusive mass transfers, have been unable to solve this problem. This manuscript presents a new mass-transfer (MT) algorithm that enables MTPT methods to accurately solve the problem of discontinuous D(x). To achieve this, we derive a semi-analytical solution to the discontinuous D(x) problem by employing a predictor-corrector approach, and we use this semi-analytical solution as the weighting function in a reformulated MT algorithm. Furthermore, this semi-analytical solution is generalized for cases with multiple 1D interfaces as well as for 2D cases, including a 2 × 2 tiling of 4 subdomains that corresponds to a numerically-generated diffusion field. The solutions generated by this new mass-transfer algorithm closely agree with an analytical 1D solution or, in more complicated cases, trusted numerical results, demonstrating the success of our proposed approach.

42 ENGINEERING↗
DOE OSTI2021

Mixing Cell Model: A One-Dimensional Numerical Model for Assessment of Water Flow and Contaminant Transport in the Unsaturated Zone

This report describes the Mixing Cell Model code, a one dimensional model for water flow and solute transport in the unsaturated zone under steady state or transient flow conditions. The model is based on the principles and assumptions underlying mixing-cell model formulations. The unsaturated zone is discretized into a series of independent mixing cells. Each cell may have unique hydrologic, lithologic, and sorptive properties. Ordinary differential equations describe the material (water and solute) balance within each cell. Water-flow equations are derived from the continuity equation, assuming that unit gradient conditions exist at all times in each cell. Pressure gradients are considered implicitly through model discretization. Unsaturated hydraulic conductivity and moisture contents are determined by the material specific moisture-characteristic curves. Solute-transport processes include explicit treatment of advective processes, first order chain decay, and linear sorption reactions. Dispersion is addressed through implicit and explicit dispersion. Implicit dispersion is an inherent feature of all mixing-cell models and originates from the formulation of the problem in terms of mass balance around fully mixed volume elements. Expressions are provided that relate implicit dispersion to the physical dispersion of the system. Two FORTRAN codes were developed to solve the water flow and solute-transport equations: (1) the Mixing Cell Model for Flow (MCMF) solves transient water-flow problems and (2) the Mixing Cell Model for Transport (MCMT) solves the solute-transport problem. The transient water-flow problem is typically solved first by estimating the water flux through each cell in the model domain as a function of time using the MCMF code. These data are stored in either ASCII or binary files that are later read by the solute transport code (MCMT). Code output includes solute pore water concentrations, water and solute inventories in each cell and at each specified output time, and water and solute fluxes through each cell and specified output time. Computer run times for coupled transient water flow and solute transport were typically several seconds on a 2 GHz Intel Pentium IV desktop computer. The model was benchmarked against analytical solutions and finite element approximations to the partial differential equations (PDE) describing unsaturated flow and transport. Differences between the maximum solute flux estimated by the mixing cell model and the PDE models were typically less than two percent. This revision includes an option for a fixed concentration lower boundary condition for diffusive fluxes for versions 020321 and later.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Reaction Diffusion Modelling of 3D Pillar Electrodes in Single-Catalyst CO 2 Reduction Cascades

Effective electrochemical CO 2 reduction to liquid fuels requires that the local catalytic environment facilitates the desired reactivity, yet a microscopic understanding of this environment is difficult to achieve from experiment alone. In this work, a 3D reaction-diffusion model was developed to explore the effects of electrode surface area and local geometry on the performance of a heterogeneous catalyst that performs a two-step CO 2 reduction cascade reaction to CO and then CH 3 OH under aqueous conditions. Kinetic parameters for the model were inspired by experimental results using a cobalt phthalocyanine (CoPc) catalyst. Three-dimensional architectures composed of arrays of square pillars with varying dimensions and either smooth or periodically modulated surfaces were tested, revealing the extent to which geometry modulates the performance of the cascade reactions. Although structural variations modulate local concentration gradients, we find that electrochemically active surface area predominantly governs the overall cascade reaction. Moreover, the results suggest that supersaturation of CO, with concentrations up to ten-fold higher than the equilibrium solubility limit, might be critical for more efficient conversion to CH 3 OH. For any given geometry, the spatially averaged ratio of [CO] to [CO 2 ] is dictated by the electrochemically active surface area and determines the yield of CH 3 OH. For a fixed surface area, geometries that spatially confine the electrolyte yield moderate local [CO] to [CO 2 ] ratios within small volumes. In contrast, less confining geometries result in a broader distribution of local ratios spread over larger volumes, with both configurations yielding the same spatially averaged [CO] to [CO 2 ] ratio. These insights provide valuable design principles—highlighting the critical importance of surface area and possibly CO supersaturation—for engineering advanced electrode architectures that leverage intermediate trapping and CO supersaturation to enhance overall performance in tandem CO 2 reduction systems.

COMSOL↗
DOE OSTI2021

Elucidating Mass Transport Regimes in Gas Diffusion Electrodes for CO 2 Electroreduction

Gas diffusion electrodes (GDEs) have shown promising performance for the electrochemical reduction of CO 2 (CO 2 R). In this study, a resolved, pore scale model of electrochemical reduction of CO 2 within a liquid-filled catalyst layer is developed. Three CO 2 mass transport regimes are identified in which the CO 2 penetration depth is controlled by CO 2 consumption in the electrolyte, CO 2 conversion along the solid-electrolyte double-phase boundaries (DPBs), and CO 2 conversion concentrated around the gas–solid–electrolyte triple-phase boundaries (TPBs). While it is possible for CO 2 R to be localized around the TPBs, in systems with submicron pore radii operating at <1 A cm –2 CO 2 R will be distributed across the DPBs within the catalyst layer. This validates the assumption of pore-scale uniformity implicit in popular, volume-averaged GDE models. Overall, the CO 2 conversion efficiency depends strongly on the governing mass transport regime, and operational-phase diagrams are constructed to guide the catalyst layer design.

10 SYNTHETIC FUELS↗
DOE OSTI2023

AN ACID BAKING APPROACH TO ENHANCE RARE EARTH ELEMENT RECOVERY FROM BITUMINOUS COAL SOURCES

Rare earth elements (REE) are a group of 17 elements typically classified as light and heavy rare earth elements, which play a crucial role in developing the latest technologies for energy, defense, and medical sectors. Even though REEs have been found in more than 200 minerals, only bastnaesite, monazite and xenotime are commercially exploited for REE extraction. However, the recent exponential increase in REE demand has spurred countries such as the United States into research for the extraction of REEs from secondary sources such as coal, acid mine drainage, and coal ash. Several coal sources (e.g., Fire Clay seam coal) across the United States have been identified to contain elevated concentrations of rare earth elements (>600 ppm), and various researchers have investigated the feasibility of both physical and hydrometallurgical extraction techniques for rare earth concentration and subsequent extraction. However, both physical and direct leaching were concluded to be inefficient for RE beneficiation and extraction due to low recoveries. Alternatively, thermal treatment provides promising means for RE recovery from bituminous coal sources. The positive impact of thermal treatment/calcination was established to be due to the decarbonization and dehydroxylation of the clays, which released entrapped rare earth elements within the dominant minerals and converted them into an acid-soluble form. Nonetheless, the improvement in recovery was limited to the light REEs (LREEs) with an insignificant increase in the heavy REEs (HREE). It was demonstrated that the light and heavy REEs in the material were associated with difficult-to-leach minerals such as monazite, xenotime, and zircon, which were not decomposed by simple calcination due to their high thermal stability. Hence, roasting in the presence of chemicals was necessary to ensure the decomposition of those REE-containing minerals. As such, this study was focused on the acid baking treatment of bituminous coals with an aim to enhance REE recovery, especially HREEs. Based on the presence of REE minerals like monazite and xenotime, three pre-leach treatment methods, i.e., 1) roasting, 2) direct acid baking, and 3) acid baking after roasting were investigated. Roasting tests at 600 ⁰C revealed that the recovery of light REEs (LREEs) was enhanced while the recovery of HREEs remained relatively unaffected. LREE and HREE recovery values of 38.3% and 21.3%, respectively, were achieved using a 50 g/L sulfuric acid solution at 5% solid concentration and a solution temperature of 75 ⁰C for 2 hours. Comparatively, direct acid baking at 250 °C provided substantial increased LREE and HREE recovery values to approximately 49.4% and 53.0%, respectively, using an equivalent acid dosage. Recoveries were maximized to 77.0% and 79.6% for LREE and HREE, respectively, by roasting followed by acid baking. Similar results were obtained from the treatment of a second bituminous coal source. Due to strong correlations between REE and Al recovery values, tests were performed on kaolinite and illite, which were prominent clay minerals within the source coals. These experiments revealed that the REE recovery improvements were likely a result of dehydroxylation of clays and subsequent release and decomposition of REE-bearing minerals such as monazite, xenotime and zircon. Subsequently, a parametric study was conducted to identify the impact of acid baking parameters on rare earth element recovery. The factors investigated using a three-level statistical experimental program were acid baking time, acid solution concentration, baking temperature, and acid solution-to-solids ratio which were found to significantly impact REE and contaminant element (Al, Fe, and Ca) recovery. An increase in baking temperature up to around 250 ⁰C improved the light and heavy REE recovery values by more than 50 absolute percentage points relative to performances achieved when direct leaching. As aforementioned, acid baking was needed to both decompose the clay minerals and liberate the REE minerals, which allowed access for the acid to solubilize the REEs. Acid concentration of the solution used for acid baking was studied as a means of minimizing the amount of acid needed to achieve a target REE recovery. However, thermo-gravimetric and differential scanning calorimetry analysis (TGA-DSC) of sulfuric acid under oxidizing atmosphere revealed that the addition of water decreased the evaporation temperature, which explains the lower REE recovery values obtained when using acid concentrations less than 100%. Using pure sulfuric acid at an acid-to-solid ratio of 0.8:1 resulted in recovery values of around 70% for both LREEs and HREEs. The decomposition reaction time was relatively quick with 65% of the TREEs recovered within the first 10 minutes. Following the identification of optimum operating conditions through the parametric study, a systematic leaching study was carried out to examine the impact of leaching parameters, such as solid-to-liquid (S/L) ratio, temperature, and time, on REE recovery using acid baking conditions of 1:1 acid to solids ratio at 250 °C for 30 minutes. The solid-to-liquid ratio was varied from 1-20% by weight at 25 °C, 50 °C, and 75 °C solution temperatures. The results indicated that reaction time and solution temperature considerably impacted the recovery of heavy and light REEs. Interestingly, LREE recovery reduced from 68% at 5% S/L to 58% at 20%S/L, whereas HREE recovery of 78% remained unaffected. The decrease in the LREE recovery was determined to be due to La and Ce precipitation, likely through isomorphic substitution with calcium in gypsum. The kinetic data indicated that 67% LREE and 77% HREE recovery could be obtained within the first 15 minutes of the reaction, suggesting fast reaction kinetics. Furthermore, raising the solution temperature from 25 °C to 75 °C increased the LREE and HREE recovery from 60% and 32% to 67% and 79%, respectively. The kinetic modeling results demonstrated that the rate-limiting step in the LREE dissolution was diffusion and chemical reaction, whereas the HREE extraction was controlled by only chemical reaction. The leaching study concluded that using 20% S/L at 75 °C for 15 minutes maximized LREE and HREE recovery. The lab-scale precipitation study showed that Fe and Al in solution could be removed at pH 4.5 followed by REE precipitation at pH 6.0 using 6 mol/L NaOH. Finally, the bench-scale data was used to develop a process flowsheet for REE recovery from low-grade bituminous coal sources using acid baking. Finally, based on the proposed flowsheet, a concentrated RE-cake obtained through selective precipitation at pH 6.5 was re-leached using HCl at pH 1.5. The resultant leachate was used to identify the impact of various operating parameters on REE recovery and purity with an aim to maximize REE precipitation efficiency while minimizing the oxalic acid dosage. The operating parameters for this investigation were oxalic acid dosage, iron (III) contamination, solution pH and temperature. The resultant model suggested that oxalic acid dosage and reaction pH are the most significant factors for the REE precipitation efficiency, followed by the interaction of oxalic dosage and Fe concentration. Test results indicated that increasing the oxalic acid concentration from 0g/L to 80g/L improved the REE precipitation efficiency from approximately 4.2% to 95.0%. Furthermore, raising the solution pH from 0.5 to 2.5 considerably enhanced the precipitation efficiency from 0.0% to 98.9%. A solution temperature elevation decreased REE recovery, which indicated an exothermic reaction between REEs and oxalate anions. Finally, a high level of Fe contamination adversely impacted REE precipitation efficiency. The speciation analysis revealed that the dominant iron species in the solution system were Fe-(C₂O₄)₃³⁻, Fe-(C₂O₄)²⁻, and Fe-(C₂O₄)⁺, which consumed the majority of the oxalate anions

rare earth elements, acid baking, high-temperature↗
NASA NTRS1993

Influence of test configuration on the combustion characteristics of polymers as ignition sources

The experimental evaluation of polymers as ignition sources for metals was accomplished at the NASA White Sands Test Facility (WSTF) using a standard promoted combustion test. These tests involve the transient burning of materials in high-pressure oxygen environments. They have provided data from which design decisions can be made; data include video recordings of ignition and non-ignition for specific combinations of metals and polymers. Other tests provide the measured compositions of combustion products for polymers at select burn times and an empirical basis for estimating burn rates. With the current test configuration, the detailed analysis of test results requires modeling a three-dimensional, transient convection process involving fluid motion, thermal conduction and convection, the diffusion of chemical species, and the erosion of sample surface. At the high pressure extremes, it even requires the analysis of turbulent, transient convection where the physics of the problem are not well known and the computation requirements are not practical at this time. An alternative test configuration that can be analyzed with a relatively-simple convection model was developed during the summer period. The principal change constitutes replacing a large-diameter polymer disk at the end of the metal test rod with coaxial polymer cylinders that have a diameter nearer to that of the metal rod. The experimental objective is to assess the importance of test geometries on the promotion of metal ignition by testing with different lengths of the polymer and, with an extended effort, to analyze the surface combustion in the redesigned promoted combustion tests through analytical modeling of the process. The analysis shall use the results of cone-calorimeter tests of the polymer material to model primary chemical reactions and, with proper design of the promoted combustion test, modeling of the convection process could be conveniently limited to a quasi-steady boundary layer analysis where the economical solution of parabolic equations is involved. The products for the summer period are: (1) a conceptual-level redesign of the test apparatus, and (2) the development and use of an approximate integral boundary layer analysis to demonstrate the influence of geometry changes prior to testing. A computer code STAN5, an accurate numerical boundary layer model whose earlier versions were developed for the NASA Lewis Research Center by the Fellow, also was installed and validated on the WSTF and New Mexico State University computer systems as a starting point in the development of a more detailed fluid mechanics and combustion model.

Julien, Howard L.↗
NASA NTRS2026

Development of a One-Domain Volume-Averaged Navier–Stokes Solver

The interaction between a high-enthalpy flow and a thermal protection material is inherently multiscale and multiphysics. In conventional aerothermal analyses, the external flow and material response are generally modeled using separate computational domains coupled through boundary conditions at the material surface. Although this approach has supported many practical applications, it requires assumptions about the location and behavior of the interface and may become difficult to apply when material decomposition, internal reactions, and surface recession substantially alter the porous structure. This report presents the development of a one-domain formulation in which the free-fluid and porous-material regions are represented within a single computational domain. The formulation is based on the volume-averaged Navier–Stokes (VANS) equations, derived from the governing equations for reacting, compressible flow and condensed material. Volume averaging transfers the influence of the unresolved material microstructure to the macroscale equations through effective transport properties, interfacial source terms, and dispersion fluxes. Particular attention is given to regions in which porosity and permeability vary rapidly, including the diffuse transition between a porous material and the surrounding fluid. The resulting equations are implemented in the Porous-material Analysis Toolbox based on OpenFOAM (PATO). The report describes the pressure–velocity coupling strategy used by the solver, examines spatial filtering techniques for deriving effective properties, and evaluates the influence of a smoothly varying interface permeability. Numerical demonstrations include canonical porous-flow configurations, a flow-tube configuration representative of FiberForm® permeability experiments, and the oxidation of a porous carbon material. The purpose of this work is to establish a mathematical and computational foundation for a unified treatment of flow and thermal protection material response. The present formulation is intended to support the progressive inclusion of additional physical processes, including multicomponent transport, finite-rate gas–surface chemistry, pyrolysis, internal oxidation, and material recession. It also provides a framework for connecting pore-scale simulations and microstructural characterization with macroscale aerothermal-response calculations. This report is intended for researchers and engineers working in computational fluid dynamics, porous-media transport, material response, and thermal protection system modeling. It documents both the theoretical development and the initial numerical assessment of the one-domain approach, while identifying the closure of effective and dispersion terms as an important subject for continued investigation.

Ablation↗
DOE OSTI2022

Electrochemical Ammonia Compression

An electrochemical (EC) compressor is a solid-state compression device. For decades, researchers have been studying EC compressors for applications for a variety of energy systems. As the world transitions away from fossil fuel energy sources, EC compression emerges as a promising technique for energy storage, specifically energy stored in the form of pressurized ammonia. Ammonia is also a commonly used refrigerant. The present studies examine the viability of EC compression for ammonia storage and refrigeration. EC ammonia compression increases the concentration of an ammonia-hydrogen mixture via the input of electrical energy and a series of chemical reactions. A polymer electrolyte membrane separates the low-concentration side from the high-concentration side of the device. On the low-pressure (anode) side hydrogen atoms oxidize and react with ammonia molecules, forming positively charge ammonium ions. The ions traverse the membrane electrolytically. The ammonium ions are reduced and revert back to hydrogen and ammonia upon reaching the high-pressure (cathode). An external circuit provides the current needed to sustain the reactions. In this project, we studied the performance of the ammonia EC compressor under a variety of different conditions. We replicated preliminary data and using a small cell with 5 cm2 of active area. We demonstrated the operation of larger cells with 100 cm2 of active area. Further, we used a commercially available hydrogen EC compressor stack to analyze the scaled-up compressor performance. While previous experiments examined only the transient EC compressor performance, we developed test facilities that allowed the compressor to reach steady state. We analyzed the effects of pressure and current on the EC compressor performance. We measured the flow rates of gas leaving the compressors and analyzed the composition using gas chromatography. We tested methods of separating ammonia from the effluent vapor, which contained hydrogen and water vapor. We found that back diffusion adversely affected the performance, especially when we maintained high pressure lifts.

25 ENERGY STORAGE↗
DOE OSTI2025

Visualizing Size-Dependent Dynamics of CeO 2-δ {100}-Supported CoO x Nanoparticles Under CO 2 Hydrogenation Conditions

Carbon dioxide is a major greenhouse gas. In order to optimize processes focused on its chemical valorization, one needs detailed information about the effects of CO 2 and/or CO 2 /H 2 mixtures on the structure and morphology of metal/oxide catalysts. Here, in this study, the evolution of a catalyst with cobalt supported on CeO 2 -cube nanostructures under CO 2 hydrogenation conditions was investigated by using a set of in situ characterization techniques (X-ray absorption fine structure, X-ray diffraction, diffuse reflectance infrared Fourier transform spectroscopy, and environmental transmission electron microscopy (TEM)). The {100} facets of the ceria support displayed an unexpectedly high stability due to strong interactions with the aggregates of cobalt oxide. A significant influence of interfacial bonding between CoO x and CeO 2-δ {100} is evident through a clear preference in the orientation of CoO x nanoparticles (NPs) with respect to the substrate. For initially reduced Co/CeO 2 -cube nanostructures, a kinetically controlled oxidation of cobalt upon the introduction of CO 2 was observed during the early stages of CO 2 hydrogenation. Environmental TEM revealed the size-dependent morphological behavior of cobalt oxide NPs due to strong interactions with the CeO 2 {100} surface. When the environment was switched from H 2 to a mixture of H 2 and CO 2 at 250 °C, small CoO x NPs (in the largest dimension < 2.5 nm) rapidly transformed from a pyramidal three-dimensional (3D) form to a planar, monatomic layer attached to the concurrently oxidized CeO 2-δ {100} surface. This maximizes the number of sites available for the binding of CO 2 or reaction intermediates. The shape transformation reflected the oxophilic character of cobalt and strong metal–support interactions. The removal of CO 2 from the gas phase led to a reduction of the cobalt oxide NPs by hydrogen and a reversible two-dimensional → 3D transformation. In contrast, no significant morphological changes, apart from further oxidation, were observed for big CoO x NPs (in the largest dimension > 3 nm). These trends are not seen for nanoparticles of noble metals. The observed morphological and structural changes in the small CoO x NPs affected the stability of reaction intermediates and modified the selectivity of the CoO x /CeO 2 catalyst system for methane production.

36 MATERIALS SCIENCE↗
NASA NTRS2003

Self-Propagating Frontal Polymerization in Water at Ambient Pressure

Advances in polymer chemistry have led to the development of monomers and initiation agents that enable propagating free-radical polymerization fronts to exist. These fronts are driven by the exothermicity of the polymerization reaction and the transport of heat from the polymerized product to the reactant monomer/solvent/initiator solution. The thermal energy transported to the reactant solution causes the initiator to decompose, yielding free radicals, which start the free radical polymerization process as discussed in recent reviews. The use of polymerization processes based on propagating fronts has numerous applications. Perhaps the most important of these is that it enables rapid curing of polymers without external heating since the polymerization process itself provides the high temperatures necessary to initiate and sustain polymerization. This process also enables more uniform curing of arbitrarily thick samples since it does not rely on heat transfer from an external source, which will necessarily cause the temperature history of the sample to vary with distance from the surface according to a diffusion-like process. Frontal polymerization also enables filling and sealing of structures having cavities of arbitrary shape without having to externally heat the structure. Water at atmospheric pressure is most convenient solvent to employ and the most important for practical applications (because of the cost and environmental issues associated with DMSO and other solvents). Nevertheless, to our knowledge, steady, self-propagating polymerization fronts have not been reported in water at atmospheric pressure. Currently, polymerization fronts require a high boiling point solvent (either water at high pressures or an alternative solvent such as dimethyl sulfoxide (DMSO) (boiling point 189 C at atmospheric pressure.) Early work on frontal polymerization, employed pressures up to 5000 atm in order to avoid boiling of the monomer/solvent/initiator solution. High boiling point solutions are needed because in order to produce a propagating front, a high front temperature is needed to produce sufficiently rapid decomposition of the free radical initiator and subsequent free radical polymerization and heat release at a rate faster than heat losses remove thermal energy from the system. (While the conduction heat loss rate increases linearly with temperature, the free radical initiator decomposition is a high activation energy process whose rate increases much more rapidly than linearly with temperature, thus as the temperature decreases, the ratio of heat loss to heat generation increases, eventually leading to extinction of the front if the temperature is too low.) In order to obtain atmospheric pressure frontal polymerization in water, it is necessary to identify a monomer/initiator combination that is water soluble and will not extinguish even when the peak temperature (T*) is less than 100 C. In this work acrylic acid (AA) was chosen as the monomer because is it one of the most reactive monomers and can polymerize readily at low temperatures even without initiators. Ammonium persulfate (AP) was chosen as the initiator because it decomposes readily at low temperatures, produces relatively few bubbles and is commercially available. The propagation rates and extinction conditions of the fronts are studied for a range of AA and AP concentrations. Small amounts of fumed silica powder (Cab-o-sil, Cabot Corporation) were added to the solutions to inhibit buoyancy induced convection in the solutions; future studies will investigate the effects of buoyant convection within the solutions.

Olten, Nesrin↗
DOE OSTI2021

Influence of microstructure on replacement and porosity generation during experimental dolomitization of limestones

Replacement reactions commonly alter the multiscale pore structures of rocks during fluid-rock interactions. Analysis of these processes in various model fluid-rock systems during controlled laboratory experiments provides insights into the origins of microstructures found in natural materials. This study focused on understanding the effects of initial starting material permeability and resultant differences in transport pathways on porosity and mineralogical changes during limestone dolomitization.A series of replacement experiments (32–317 days in duration) have been conducted in which 1.59 cm (5/8 in.) diameter cores of two different limestones were reacted with saturated MgCl2 solutions at 200 °C. The Texas Cream (Austin Chalk) is a high-porosity, high-permeability limestone, whereas both the porosity and permeability of the Carthage Marble (Burlington Limestone) are relatively low. Altered limestones were imaged using scanning electron microscopy with energy-dispersive X-ray spectroscopy (SEM-EDX), Time-of-Flight Secondary Ion Mass Spectrometry (ToF–SIMS) and electron microprobe analysis (EMPA). A representative grain boundary of the low-porosity limestone was targeted for a focused ion beam (FIB) lift-out and characterized using transmission electron microscopy (TEM). These results were coupled with analyses of radial changes in the porosity distribution of the core derived from X-ray and neutron small- and ultra-small angle scattering ((U)SANS/(U)SAXS).The high-porosity/permeability limestone showed a four times faster bulk replacement rate than the lower-porosity/permeability material, and a different mechanism of porosity development. For the low-porosity limestone, a two-stage replacement occurred, with the reacted region of the core consisting of an inner rim in which the limestone was replaced by two calcite-dolomite solid solutions, and an outer rim in which the dolomite was replaced by magnesite. Elongated pores formed along grain boundaries at the initial limestone/dolomite reaction interface, and additional nanometer-scale porosity was formed at the secondary magnesite replacement rim. Grain boundaries were identified as preferential pathways for transport leading to dolomitization and a grain boundary diffusion rate was calculated based on microstructural characterization. In contrast, replacement in the high-porosity limestone was accompanied by porosity generation through replacement of individual grains by dolomite throughout the sample and, in longer runs, magnesite in outer parts of the core. These observations emphasize that both the mechanisms of the replacement reaction and the microstructure and chemistry of the replaced product are contingent on the initial structure of the starting material.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Modeling Interfaces to Support Low-Level Waste Disposal System Performance Assessments - 20366

In low-level waste (LLW) disposal sites, interfaces between cementitious materials, used as waste forms and/or engineered barriers, and the surrounding soil or backfill material are often encountered. Reactions across these interfaces may lead to chemical and structural alteration of the cementitious and backfill materials that may ultimately affect long-term performance. For example, the ingress of carbon dioxide from soil gas into a waste tank concrete shell at the Hanford Site or Savannah River Site (SRS) or a low-level waste disposal vault at SRS leads to carbonation of cement hydration products (e.g. portlandite, calcium-silicate hydrate or C-S-H, and ettringite). The result of carbonation is a decrease in pH in the cement paste portion of the tank wall or vault concrete structure that could ultimately lead to de-passivation of the embedded structural steel and, eventually, to cracking. Cracking can lead to increased ingress of water into the structure, corresponding increased release of constituents of concern, and to subsequent increased transport of these constituents to the surrounding environment. The mobility of trace constituents (e.g., radionuclides of concern) in waste forms may increase in response to changes in pH, pore structure, and mineralogical gradients. Depletion of portlandite results in subsequent decalcification of C-S-H that leads to changes in the cement strength and may lead to structural failure. At the boundary between waste forms and concrete vaults, the migration of sulfate ions from the salt waste form (e.g., the saltstone waste form used for LLW disposal at the Savannah River Site in Saltstone Disposal Units, or SDUs) into the barrier or vault concrete has been predicted to result in expansive mineral phase formation (or 'sulfate attack'). Expansive mineral formation could result in cracking and potential loss of structural integrity in the vault concrete. In the performance assessment (PA) for the SRS SDUs, prediction of both the carbonation and sulfate attack fronts are critical to assessing the long-term performance of these LLW disposal vaults. A general purpose geochemical reactive transport model has been developed in LeachXS/ORCHESTRA to evaluate typical interfaces for LLW disposal environments, including SDUs at SRS. The model may be used to predict the long-term performance of interfaces between cementitious materials and saltstone or backfill materials with respect to primary phases by considering mass transport, cement chemistry, geochemical speciation, and multi-ionic diffusion over interfaces between different materials. The present study compares the performance and interactions between vault concrete and saltstone waste form for a representative SDU scenario. The materials were characterized using pH-dependence leaching test (EPA 1313) and semi-dynamic transport tests (EPA 1315) with test results used to parameterize the interface model. The development of a credible inter model provides an important basis for supporting future PAs where carbonation and sulfate attack are important aging and degradation mechanisms. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS2002

Association between simian virus 40 and non-Hodgkin lymphoma

BACKGROUND: Non-Hodgkin lymphoma has increased in frequency over the past 30 years, and is a common cancer in HIV-1-infected patients. Although no definite risk factors have emerged, a viral cause has been postulated. Polyomaviruses are known to infect human beings and to induce tumours in laboratory animals. We aimed to identify which one of the three polyomaviruses able to infect human beings (simian virus 40 [SV40], JC virus, and BK virus) was associated with non-Hodgkin lymphoma. METHODS: We analysed systemic non-Hodgkin lymphoma from 76 HIV-1-infected and 78 HIV-1-uninfected patients, and non-malignant lymphoid samples from 79 HIV-1-positive and 107 HIV-1-negative patients without tumours; 54 colon and breast carcinoma samples served as cancer controls. We used PCR followed by Southern blot hybridisation and DNA sequence analysis to detect DNAs of polyomaviruses and herpesviruses. FINDINGS: Polyomavirus T antigen sequences, all of which were SV40-specific, were detected in 64 (42%) of 154 non-Hodgkin lymphomas, none of 186 non-malignant lymphoid samples, and none of 54 control cancers. This difference was similar for HIV-1-infected patients and HIV-1-uninfected patients alike. Few tumours were positive for both SV40 and Epstein-Barr virus. Human herpesvirus type 8 was not detected. SV40 sequences were found most frequently in diffuse large B-cell and follicular-type lymphomas. INTERPRETATION: SV40 is significantly associated with some types of non-Hodgkin lymphoma. These results add lymphomas to the types of human cancers associated with SV40.

Non-NASA Center↗
DOE OSTI2025

MOSCATO Development and Integration in Fiscal Year 2025: Implementation of Multiphase, Multiphysics Modeling Capabilities for Molten Salt Systems

MOSCATO (Molten Salt Chemistry and Transport) is a multiphysics code that provides high-fidelity, coupled simulations of fluid flow, heat transfer, mass transfer, chemistry, electrochemical phenomena, and alloy corrosion for molten salt systems. In FY25, significant developments were made to the code package, enhancing its capabilities for modeling all relevant phenomena within operating moltens salt reactors (MSRs). The developments and activities in FY25 included: 1. Implementation of Level-Set methods to enable modeling of single-bubble behavior in molten salts. In FY25, the Level-Set two-phase flow modeling implementation was improved to simulate single bubble behavior with molten salt media. The large density and viscosity ratios between typical gases and molten salt liquids present challenges for these types of numerical solvers. With enhancements to the pressure projection method, MOSCATO’s Level-Set solver was able to be successfully validated to experiments related to helium bubble rise in stagnant molten salt. The simulated bubble rising velocity showed reasonable good agreement with experimental measurements. The bubble shape and dynamics were also visually compared with experimental snapshots, demonstrating a good qualitative match. 2. Generation of mass transfer correlations for multiphase flow systems. To enable calculations of the tritium transport across the interface between gas bubbles and salt, we modeled high- Schmidt-number mass transfer around a sphere across a broad range of Reynolds numbers. The mesh near the sphere surface was highly refined to resolve steep concentration gradients caused by the low diffusion coefficient. Literature-based mass transfer correlations were compared with the numerical results, and modifications were proposed to improve agreement, particularly at higher Schmidt numbers. These mass transfer correlations were subsequently provided to other national laboratories to help enable high quality mass transfer simulations using lower-order solvers under development within the NEAMS program. 3. Preliminary implementation of a bubbly flow solver. To model bubbly flow in molten salt, we implemented a bubbly flow solver for void fractions less than 5%. To do so, an algebraic relative velocity model that assumes small bubbles with rapid momentum equilibration was added to MOSCATO to compute bubble velocities. Preliminary comparisons with experimental data showed reasonable agreement, and further improvements are underway. 4. Generation of mass transfer correlations for MSRE subchannel The Molten-Salt Reactor Experiment (MSRE) was a landmark historical project that demonstrated the feasibility of molten-salt reactor technology. The MSRE campaign also generated a significant body of experimental data and reports that continue to support molten-salt–related research. In this report, the MSRE core subchannel was used as the reference geometry for a mass transfer study performed with MOSCATO. The geometry and computational mesh were adapted from a previous study, providing adequate resolution for the relatively low Reynolds number in this case. Additional mesh refinement was applied to reach higher Schmidt numbers, enabling the derivation of a reliable mass-transfer correlation for the present scenario. 5. Simulations of oxygen ingressions into molten salt. In the previous fiscal year, we initiated a study to simulate oxygen ingression in stagnant salt. As oxygen enters the salt through its surface, it reacts with Ce 3+ to form solid CeO 2 and other reaction products. To more fully capture the complex diffusion-convection-reaction mechanisms, capabilities for modeling natural convection in the salt vessel were added. These were needed as the flow of the ingressed gas induced flow in the salt caused by surface shear and non-isothermal effects. With these updated physics in place, we were able to successfully reproduce the experimental results for the rate of change of CeCl 3 concentrations versus time. 6. Flow corrosion model validation. In FY24, MOSCATO’s corrosion model was validated against static corrosion experiments. In FY25, this work was extended to a flow corrosion experiment, where FLiNaK salt was driven by natural convection, with initial salt impurities to initiate corrosion. Despite uncertainties in parameters such as elemental diffusion coefficients in the alloy and unknown H + concentrations, the simulations achieved good agreement with experimental results, especially in predicting sample mass losses.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Enhanced Electrocatalytic and Cathode‐Electrolyte Interfacial Properties With a Pr‐Based Simple Perovskite/Ruddlesden‐Popper Nanocomposite Cathode in Protonic Ceramic Fuel Cells

The sluggish kinetics and poor stability of the oxygen reduction reaction (ORR) remain the primary bottleneck for achieving high performance in protonic ceramic fuel cells (PCFCs) at intermediate temperatures (400–650°C). In this work, a Pr-based nanocomposite cathode comprised of simple perovskite phase (PrNi 0.7 Co 0.3 O 3-δ ) and Ruddlesden-Popper phase (Co-doped Pr 4 Ni 3 O 10+δ ) is developed. Although PrNi 0.7 Co 0.3 O 3-δ solely stands as a good cathode with facile proton transfer, combining the superior catalytic activity against oxygen on the Ruddlesden-Popper phase boosts the ORR performance further. The designed nanocomposite cathode outperforms the simple perovskite cathode, attributed to enhanced oxygen absorption and surface diffusion with the Ruddlesden-Popper phase. A precursor-based cathode deposition technique is also developed to achieve cathode grain sizes of ∼100 nm. A single cell with the nanocomposite cathode delivers a peak power density of 1.38 W cm −2 at 650°C, among the highest in reported PCFCs with Pr-based cathodes, with a small degradation rate of 0.145 mV h −1 during 250 h stability test. Further investigation of cathode-electrolyte interface revealed interfacial PrO 2 phase formation, promoted by abundant Pr 6 O 11 in the nanocomposite precursor powder, thereby improving both ohmic resistance and stability. These findings highlight the effectiveness of the nanocomposite cathode and underscore its advantages on interfacial properties.

08 - HYDROGEN↗