Search NASA⌕ Search

SEARCH · Search NASA

Results for “analysis of kinetic data”

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 289 records · Page 16

Kinetics of the Carbothermic Reduction of Manganese Oxide from Slag

Experiments were performed using a range of test conditions to elucidate the rate controlling step during the reaction of liquid iron-carbon droplets and slags containing manganese oxide. Four conditions were tested in the system: initial MnO content in the slag (5, 10, and 15 wt pct), initial carbon content of the metal (1, 2.5, 4.3 wt pct), initial droplet mass (0.5, 1.0, and 1.5 g), and reaction temperature (1823 K [1550 °C], 1873 K [1600 °C], and 1923 K [1650 °C]). Data were collected using the Constant Volume Pressure Increase (CVPI) technique which tracked the continuous pressure increase in the sealed furnace over time. Samples were quenched at the end of each experiment and chemistry was checked using LECO Carbon Analysis and ICP (Inductively Coupled Plasma) for manganese. The rate of reaction can be broken into a faster initial period related to internal CO formation, and a slower second reaction controlled by a complex mechanism involving transport of oxygen from slag to metal via CO{sub 2} and decomposition of the CO{sub 2} at the gas–metal interface.

36 MATERIALS SCIENCE↗

Validation Assessment of Turbulent Reacting Flow Model Using The Area-Validation Metric on Medium-Scale Methanol Pool Fire Results

Accident analysis and ensuring power plant safety are pivotal in the nuclear energy sector. Significant strides have been achieved over the past few decades regarding fire protection and safety, primarily centered on design and regulatory compliance. Yet, after the Fukushima accident a decade ago, the imperative to enhance measures against fire, internal flooding, and power loss has intensified. Hence, a comprehensive, multilayered protection strategy against severe accidents is needed. Consequently, gaining a deeper insight into pool fires and their behavior through extensive validated data can greatly aid in improving these measures using advanced validation techniques. A model validation study was performed at Sandia National Laboratories in which a 30-cm diameter methanol pool fire was modeled using the SIERRA/Fuego turbulent reacting flow code. This validation study used a standard validation experiment to compare model results against, and conclusions have been published. The fire was modeled with a Large Eddy Simulation (LES) turbulence model with subgrid turbulent kinetic energy closure. Combustion was modeled using a strained laminar flamelet library approach. Radiative heat transfer was accounted for with a model utilizing the gray-gas approximation. In the present study, additional validation analysis is performed using the area validation metric (AVM). These activities are done on multiple datasets involving different variables and temporal/spatial ranges and intervals. In conclusion, the results provide insight into the use of the area validation metric on such temporally varying datasets and the importance of physics-aware use of the metric for proper analysis.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

On the Rate Constant for NH 2 +HO 2 and Third-Body Collision Efficiencies for NH 2 +H(+M) and NH 2 +NH 2 (+M)

In low-temperature flash photolysis of NH 3 /O 2 /N 2 mixtures, the NH 2 consumption rate and the product distribution is controlled by the reactions NH 2 + HO 2 → products (R1), NH 2 + H (+M) → NH 3 (+M) (R2), and NH 2 + NH 2 (+M) → N 2 H 4 (+M) (R3). In the present work, published flash photolysis experiments by, among others, Cheskis and co-workers, are re-interpreted using recent direct measurements of NH 2 + H (+N 2 ) and NH 2 + NH 2 (+N 2 ) from Altinay and Macdonald. To facilitate analysis of the FP data, relative third-body collision efficiencies compared to N 2 for R2 and R3 were calculated for O 2 and NH 3 as well as for other selected molecules. We report results were in good agreement with the limited experimental data. Based on reported NH 2 decay rates in flash photolysis of NH 3 /O 2 /N 2 , a rate constant for NH 2 + HO 2 → NH 3 + O 2 (R1a) of $k_{1\text{a}}$ = 1.5(±0.5) × 10 14 cm 3 mol –1 s –1 at 295 K was derived. This value is higher than earlier determinations based on the FP results but in good agreement with recent theoretical work. Kinetic modeling of reported N 2 O yields indicates that NH 2 + HO 2 → H 2 NO + O (R1c) is competing with R1a, but perturbation experiments with addition of CH4 indicate that it is not a dominating channel. Measured HNO profiles indicate that this component is formed directly by NH 2 + HO 2 → HNO + H 2 O (R1b), but theoretical work indicates that R1b is only a minor channel. Based on this analysis, we estimate $k_{1\text{c}}$ = 2.5 × 10 13 cm 3 mol –1 s –1 and $k_{1\text{b}}$ = 2.5 × 10 12 cm 3 mol –1 s –1 at 295 K, with significant uncertainty margins.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Characterizing Defect Dynamics in Silicon Carbide Using Symmetry-Adapted Collective Variables and Machine Learning Interatomic Potentials

Silicon carbide (SiC) divacancies are attractive candidates for spin-defect qubits possessing long coherence times and optical addressability. The high activation barriers associated with SiC defect formation and motion pose challenges for their study by first-principles molecular dynamics. In this work, we develop and deploy machine learning interatomic potentials (MLIPs) to accelerate defect dynamics simulations while retaining ab initio accuracy. We employ an active learning strategy comprising symmetry-adapted collective variable discovery and enhanced sampling to compile configurationally diverse training data, calculation of energies and forces using density functional theory (DFT), and training of an E(3)-equivariant MLIP based on the Allegro model. Here, the trained MLIP reproduces DFT-level accuracy in defect transition activation free energy barriers, enables the efficient and stable simulation of multidefect 216-atom supercells, and permits an analysis of the temperature dependence of defect thermodynamic stability and formation/annihilation kinetics to propose an optimal annealing temperature to maximally stabilize VV divacancies.

Computer simulations↗

In Situ Characterization of Surface Recombination in p-Si/SiO x Based Photoelectrochemical Cells

Time-resolved infrared (TRIR) and electrochemical impedance spectroscopy (EIS) were utilized to quantify surface states present at silicon oxide (SiO x ) protected crystalline p-Si electrolyte interfaces. The primary goal was to identify p-Si/SiO x photoelectrodes with both low surface recombination rates and efficient multi-electron transfer to an acceptor present in the external electrolyte. Three SiO x layers were investigated: native oxide (nOx), chemical oxide (cOx), and rapid thermal annealed (RTA) thermal oxide (tOx). Comparative study with [Ru(bpy) 3 ](PF 6 ) 2 as the electron acceptor indicated that tOx was most optimal with a small effective recombination rate, multi-electron transfer capability, and photovoltage of 500 ± 50 mV. A secondary goal was to analyze the surface recombination rates with the Shockley–Read–Hall (SRH) kinetic model. Two surface states were identified from this analysis, one closer to the CB edge (V t,1 ) and the other near the midgap (V t,2 ). EIS and SRH analyses revealed that a forming gas (5% H 2 /N 2 ) anneal (FGA) decreased surface recombination for tOx and nOx through a lower density of surface states. In the case of tOx, the infrared data indicated that V t,2 was completely removed. Here, the energetic positions of the band edges were correlated with the surface state density; low densities corresponded to more favorable potentials for inversion layer formation, which is expected to be most optimal for photocatalysis. Collectively this study indicates that the free carrier dynamics provided by TRIR represent a powerful in situ probe of the band edge and the surface state energetics in silicon based photoelectrochemical cells.

Electrochemical Impedance Spectroscopy↗

Equilibrium reconstruction of DIII-D plasmas using predictive modeling of the pressure profile

New workflows have been developed for predictive modeling of magnetohydrodynamic (MHD) equilibrium in tokamak plasmas. The goal of this work is to predict the MHD equilibrium in tokamak discharges without having measurements of the kinetic profiles. The workflows include a cold start tool, which constructs all the profiles and power flows needed by transport codes; a Grad–Shafranov equilibrium solver; and various codes for the sources and sinks. For validation purposes, a database of DIII-D tokamak discharges has been constructed that is comprised of scans in the plasma current, toroidal magnetic field, and triangularity. Initial efforts focused on developing a workflow utilizing an empirically derived pressure model tuned to DIII-D discharges with monotonic safety factor profiles. This workflow shows good agreement with experimental kinetic equilibrium calculations, but is limited in that it is a single fluid (equal ion and electron temperatures) model and lacks H-mode pedestal predictions. The best agreement with the H-mode database is obtained using a theory-based workflow utilizing pressure profile predictions from a coupled TGLF turbulent transport and EPED pedestal models together with external magnetics and Motional Stark Effect (MSE) data to construct the equilibrium. Here, we obtain an average root mean square error of 5.1% in the safety factor profile when comparing the predicted and experimental kinetic equilibrium. We also find good agreement with the plasma stored energy, internal inductance, and pressure profiles. Including MSE data in the theory-based workflow results in noticeably improved agreement with the q-profiles in high triangularity discharges in comparison with the results obtained with magnetic data only. The predictive equilibrium workflow is expected to have wide applications in experimental planning, between-shot analysis, and reactor studies.

Kinsey, J. E. (ORCID:0000000193347473)↗

A Predictive Transient Model of the TREAT-SIRIUS Experiments

To support experiments modeling and transient analysis of the NASA-sponsored SIRIUS experiments in the Transient Reactor Test Facility (TREAT), an innovative predictive transient model has been developed to simulate various transient experiments. The developed model is based on a new approach that utilizes steady-state Monte Carlo calculations along with a surrogate model based on polynomial regression to determine the core reactivity. The reactivity is supplied to an exact point kinetics model that determines the reactor power, and it is coupled to an adiabatic feedback model to determine the fuel average temperature. In this work, the newly developed model is introduced along with validation test results considering SIRIUS-1 experiments at different power levels. The initial test results show a very good agreement with experimental data at different power levels, which indicates that the model can be used to predict the power of other transient tests.

42 ENGINEERING↗

Tuning Copolymer Microstructure Using Ring-Opening Cross-Metathesis Polymerization

The capability of ring-opening cross-metathesis (RO/CM) polymerization to produce alternating copolymers was studied. By treating commercial polybutadiene (PB) with bulky oxanorbornene monomers and Ru-based olefin metathesis catalysts, alternating copolymers were produced under mild conditions with high sequence fidelities. Here, we found that alternating copolymers could be produced starting from a variety of butadiene sources including PB, cyclooctadiene (COD), or t,t,t-1,5,9-cyclododecatriene (CDT), highlighting for the first time the kinetic pathway independence of this process. Kinetic copolymerization analysis of an oxanorbonene monomer with CDT revealed that much higher monomer conversions were obtained compared with the analogous homopolymerizations and showed evidence of alternating monomer incorporation. Copolymerization of these monomers also enabled good control when targeting different molecular weights. Copolymer thermal analysis revealed a strong correlation between thermal behavior and alternating sequence fidelity, providing a second lever beyond composition to tune thermal behavior. These data demonstrate that a broad variety of polymer microstructures can be accessed via RO/CM polymerization and highlight the potential of CDT in alternating copolymer synthesis.

Foster, Jeffrey C. [Oak Ridge National Laboratory ↗

HERO CarbonSAFE Phase 2 Project in the Columbia River Basalt Group

The Hermiston, Oregon Basalt CarbonSAFE Phase II project (HERO CarbonSAFE) seeks to accelerate the deployment of commercial carbon dioxide (CO2) storage projects in basaltic rocks. Basalt CO2 storage has several advantages to conventional saline storage reservoirs including 1. The potential for rapid mineralization of CO2, 2. Associated decreases in pressure and CO2 migration risks, 3. Reduced long-term monitoring requirements with respect to plume tracking, 4. Widespread geographic distribution and, 5. Large storage potential due to thickness, porosity, and CO2 interactions with basalt. And for locations such as the Pacific Northwest, Hawaii, Iceland, India and Japan, basalts may offer the only economically feasible option for local CO2 storage. However, there are limited field-scale assessments of CO2 storage in basalt, and current carbon capture utilization and storage (CCUS) permitting and regulatory frameworks were developed for conventional saline reservoirs. HERO CarbonSAFE is designed to address research gaps and uncertainties associated with basalt storage. Specifically, the project will assess the feasibility of CO2 injection in the deep layered basalts, long-term storage (mineralization), practical approaches for large-scale implementation (50+ million metric tons of CO2 over 30 years), lithology-specific risks, and the technoeconomic potential for CO2 storage in basalts. The HERO CarbonSAFE project will assess feasibility of developing a commercial-scale (50+ million metric tons of CO2) geological storage complex within the Columbia River Basalt Group (CRBG), a layered continental flood basalt complex that underlies Calpine’s natural gas-fired Hermiston Power Project (HPP) in Hermiston, OR (Figure 1). Under this 2-year CarbonSAFE Phase II project, the HERO team will conduct a data acquisition campaign that includes drilling a stratigraphic well to a total depth of ~1,500 m into the thick layered basalts proximal to HPP. A comprehensive well logging and hydrologic testing program will be augmented with new core collected from flow zones and sealing units, and comprehensive laboratory testing to help refine the kinetic rates of mineralization. The newly acquired information will be integrated with existing data from regional wells to correlate basalt injection zone properties to develop storage hub/commercial-scale models. Using these models, the project team will evaluate injection scenarios to define the technical and economic potential for storing a minimum of 50 million metric tons of CO2 over a 30-year period, along with a robust sensitivity analysis on key parameters governing reservoir viability for sustainable injection over a commercial project lifetime. Specific technical objectives of HERO are: (1) assessing the reservoir response of a series of stacked layered reservoir flowtop sequences occurring in this area of the CRBG to commercial-scale injection volumes; (2) extending prior efforts by the project team to characterize the deep layered basalts encountered in regional studies, to leverage prior investments by U.S. Department of Energy’s (DOE) Carbon Storage program; (3) leveraging DOE’s mineralization characterization efforts to advance model parametrization for commercial scale injection of CO2 in basalts; (4) conducting risk assessments associated with scaling up to commercial storage hub injection goals, while validating DOE’s National Risk Assessment Partnership (NRAP) tools, to identify potential constraints that would prevent the CRBG from serving as a commercial-scale storage complex; (5) developing mitigation plans to address identified risks; (6) developing a commercial-scale injection and monitoring, verification and accounting (MVA) strategy; (7) utilizing computational models to define and minimize, if possible, the Area of Review (AoR) under Class VI regulations; and (8) developing a robust CO2 management strategy for CRBG that also considers a regional source/sink approach that is responsive to stakeholder needs and industrial demand. Specific institutional objectives are: (1) identifying and developing plans to mitigate the nontechnical challenges associated with the build-out of a commercial-scale storage complex within the CRBG with integrated CO2 sources; (2) implementing the community outreach plan; (3) conducting regulatory research, including a survey of issues related to pore space ownership, MVA and long-term assurance of mineralization-based storage, to support an eventual application for a UIC Class VI permit; (4) advancing the project’s plan for CO2 liability management; and (5) continuing to refine and update the project’s economic model. The final objective is the preparation of a comprehensive Site Characterization Plan that draws upon the technical and institutional feasibility assessments to prepare the project for future commercialization efforts.

58 GEOSCIENCES↗

A near magnetic-to-kinetic energy equipartition flare from the relativistic jet in AO 0235 + 164 during 2013–2019

ABSTRACT We present the multiwavelength flaring activity of the blazar AO 0235 + 164 during its recent active period from 2013 to 2019. From a discrete correlation function analysis, we find a significant ($\gt 95~{{\ \rm per\ cent}}$) correlation between radio and gamma-ray light curves with flares at longer wavelengths following flares at shorter wavelengths. We identify a new jet component in 43 GHz Very Long Baseline Array data that was ejected from the radio core on MJD $57246^{+26}_{-30}$ (2015 August 12), during the peak of the 2015 radio flare. From the analysis of the jet component, we derived a Doppler factor of δvar = 28.5 ± 8.4, a bulk Lorentz factor of $\Gamma =16.8^{+3.6}_{-3.1}$, and an intrinsic viewing angle of $\theta _{\rm v}=1.42^{+1.07}_{-0.52}\textrm {~degrees}$. Investigation of the quasi-simultaneous radio data revealed a partially absorbed spectrum with the turnover frequency varying in the range of 10−70 GHz and the peak flux density varying in the range of 0.7−4 Jy. We find the synchrotron self-absorption magnetic field strength to be $B_{\rm SSA}=15.3^{+12.6}_{-14.0}\textrm {~mG}$ at the peak of the 2015 radio flare, which is comparable to the equipartition magnetic field strength of $B_{\rm EQ}=43.6^{+10.6}_{-10.4}\textrm {~mG}$ calculated for the same epoch. Additional analysis of the radio emission region in the relativistic jet of AO 0235 + 164 suggests that it did not significantly deviate from equipartition during its recent flaring activity.

Astronomy & Astrophysics↗

Kinetic Study of Polyvinyl Chloride Pyrolysis with Characterization of Dehydrochlorinated PVC

In this paper, we study the kinetics of polyvinyl chloride (PVC) decomposition using a combination of experimental and computational approaches. Here, we develop a simplified kinetic model that contains only two steps: dehydrochlorination of PVC and further decomposition of the PVC residue. The model is consistent with density functional theory (DFT) calculations and experimental data. Dehydrochlorination is an autocatalytic reaction that starts with a tertiary chloride (Cl) and generates hydrogen chloride (HCl) and benzene as the main products. Benzene and HCl formation rates showed similar trends, indicating that HCl likely catalyzes a homolytic carbon-carbon (C-C) bond cleavage, which gives rise to benzene and an aliphatic fragment. We characterized the structure of dehydrochlorinated PVC (PVC residue) by using thermal gravimetric analysis (TGA), Fourier-transform infrared spectroscopy (FTIR), and nuclear magnetic resonance spectroscopy (NMR). FTIR and NMR results indicate that the PVC residue contains 20% quaternary carbon content, indicating a high concentration of cross-linked molecules. We predict that the most probable structure in the cross-linked centers of the PVC residue is cyclohexadiene, which is supported by DFT calculations, FTIR, and NMR.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Forward and inverse predictive transient models of TREAT using surrogate reactivity models

In this work, we present two novel approaches for predicting the power evolution and control rod height of the Transient Reactor Test Facility (TREAT) to support experiment modeling; specifically, transient analysis of the NASA-sponsored Sirius series of experiments. These approaches utilize steady-state Monte Carlo model, point kinetics model, and surrogate models to predict power evolution and control rod axial position during transient experiments. Both approaches were tested and validated against several Sirius experiments that were performed in the TREAT facility at different power levels. The validation test results show very good agreement with the experimental data, and the models were able to accurately predict the power evolution and the axial control rod position with an average error within 3.0%. Finally, this indicates that these approaches will help the reactor engineering team of the TREAT facility in preparing and predicting the power and temperature of the experiment.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Differentiating Double-Layer, Psuedocapacitance, and Battery-like Mechanisms by Analyzing Impedance Measurements in Three Dimensions

Electrochemical energy storage arises from processes that are broadly categorized as capacitive, pseudocapacitive, or battery-like. Advanced charge-storing materials that are designed to deliver high capacity at a high rate often exhibit a multiplicity of such mechanisms, which complicates the understanding of their charge-storage behavior. In this paper, we apply a “3D Bode analysis” technique to identify key descriptors for fast Li-ion storage processes, where AC impedance data, such as the real capacitance ($C'$) or phase angle ($Φ$), are represented versus the frequency ($f$) and a third independent variable, the applied DC cell voltage. For double-layer processes, a near-constant $C'$ or $Φ$ is supported across the entire voltage range, and the decrease in these values shows a near-linear decrease at higher $f$. For pseudocapacitance, an increase in $C'$ is delivered, accompanied by high $C'$ retention at higher $f$ compared to double-layer processes. Interestingly, the lower $Φ$ values, where $C'$ is highest, suggest that this is a key descriptor for pseudocapacitance, where high-rate charge storage is still facilitated within a kinetically limited regime. For battery-like processes, a high $C'$ is only observed at the voltage at which the material stores charge, while outside that voltage, $C'$ is negligible. The three-dimensional (3D) Bode analysis allows charge-storage dynamics to be mapped out in great detail with more delineation between mechanisms compared to the more frequently deployed kinetic analyses derived from cyclic voltammetry.

25 ENERGY STORAGE↗

Resolving Discrepancies Between Theory and Experiment for the NCN + H Reaction

Flame modeling studies have highlighted the role of branching in the kinetics for the prompt NO switch reaction, NCN + H. In this reaction, there is strong competition between the CH + N2 (R1) and HCN + N (R2) product channels. Increased branching towards reaction R2 promotes the subsequent formation of Fenimore or Prompt NO. Recent direct shock tube studies on this complex reaction conclude that branching to reaction R2 is much more favorable than the predicted branching from theoretical studies. The experimentally predicted prompt NO switch temperature (TS), i.e., the temperature at which k2/(k1+k2) = 0.5, is at 1670 K, in contrast to theory that predicts Ts > 3200 K. The present theory/modeling work was initiated to highlight potential causes for this significant discrepancy between experiment and theory. Our analysis of the simulations reported in the experimental studies indicate that including a fall-off representation for C2H5 radical dissociation (with C2H5 produced from C2H5I dissociation and acting as an H-atom source) has a noticeable influence on the simulations for NCN decay and HCN formation in the shock tube studies. With the inclusion of a theory-based pressure fall-off representation for C2H5 dissociation, we show that this radical persists for much longer timescales even at the high-temperatures in the shock tube studies. This persistence then facilitates the rapid conversion of reactive H-atoms to CH3 radicals at the shock tube conditions via the H + C2H5 → CH3 + CH3 reaction. An updated theoretical analysis for this classic addition-elimination reaction is also provided in this work. The rapid formation of CH3 radicals necessitates the inclusion of additional reactions to adequately simulate the NCN and HCN data. Our simulations conclude that, with the inclusion of these additional reactions, the NCN and HCN profiles can be reasonably well simulated by the most recent theoretical predictions.

Ethyl Radicals↗

Cloud based tool for analysis of chemical kinetic mechanisms

The size and complexity of chemical kinetic models used to simulate combustion of hydrocarbon fuels continue to increase with our improved understanding of the underlying physical mechanisms, along with the improved ability to utilize such models as the result of increased computational power and efficiency. As mechanisms grow beyond thousands of species and tens of thousands of reactions, it becomes increasingly difficult to manually check for errors and physical inconsistencies. In this work, we present several automated methods to check for such issues in the specification of chemical reaction models. First, we demonstrate how discontinuities in thermodynamic data can cause simulation difficulties manifested as long wall-clock times and convergence failures. To correct this type of problem, we describe an automated method for refitting thermodynamic parameters. We also outline several methods to check the timescales of reaction rate coefficients to ensure physical consistency. All the methods are made available through a web-based tool (https://combustiontools.llnl.gov) intended to aid mechanism developers and users to improve the accuracy and performance of fuel models used by the combustion community.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Protein folding from heterogeneous unfolded state revealed by time-resolved X-ray solution scattering

One of the most challenging tasks in biological science is to understand how a protein folds. In theoretical studies, the hypothesis adopting a funnel-like free-energy landscape has been recognized as a prominent scheme for explaining protein folding in views of both internal energy and conformational heterogeneity of a protein. Despite numerous experimental efforts, however, comprehensively studying protein folding with respect to its global conformational changes in conjunction with the heterogeneity has been elusive. Here we investigate the redox-coupled folding dynamics of equine heart cytochrome c (cyt-c) induced by external electron injection by using time-resolved X-ray solution scattering. A systematic kinetic analysis unveils a kinetic model for its folding with a stretched exponential behavior during the transition toward the folded state. With the aid of the ensemble optimization method combined with molecular dynamics simulations, we found that during the folding the heterogeneously populated ensemble of the unfolded state is converted to a narrowly populated ensemble of folded conformations. These observations obtained from the kinetic and the structural analyses of X-ray scattering data reveal that the folding dynamics of cyt-c accompanies many parallel pathways associated with the heterogeneously populated ensemble of unfolded conformations, resulting in the stretched exponential kinetics at room temperature. This finding provides direct evidence with a view to microscopic protein conformations that the cyt-c folding initiates from a highly heterogeneous unfolded state, passes through still diverse intermediate structures, and reaches structural homogeneity by arriving at the folded state.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Direct Air Capture Using Trapped Small Amines in Hierarchical Nanoporous Capsules on Porous Electrospun Fibers (Final Technical Report)

This report summarizes the carbon capture research and development conducted by The State University of New York at Buffalo (UB) and GTI Energy (GTI) for award “DE-FE0031969: Direct Air Capture Using Trapped Small Amines in Hierarchical Nanoporous Capsules on Porous Electrospun Fibers” sponsored by the U.S. Department of Energy (DOE). The objective of this project is to develop an innovative sorbent structure of trapped small amines in HNC embedded in PEF for DAC. This involves tailoring both sorbent and PEF materials to achieve a compact system for DAC with high capacity for CO 2 at concentrations typically available in air and at near ambient conditions. An innovative sorbent structure of trapped small amines in hierarchical nanoporous capsules (HNC) embedded in porous electrospun fibers (PEF) was developed for direct air capture (DAC). This involves tailoring both sorbent and PEF materials to achieve a compact system for DAC with high capacity for CO 2 at concentrations typically available in air and at near ambient conditions. An interfacial polymerization process was developed, which utilized loaded amines inside mesoporous silica and trimesoyl chloride (TMC) dissolved in organic solvents as the precursors, to generate a polyamide (PA) coating layer on mesoporous silica and thus trap amines. Reaction conditions, including TMC concentration, organic solvents, reaction time, etc., for interfacial polymerization were optimized to effectively trap loaded amines, and cyclic heating-cooling operation was conducted to evaluate the coating quality. Larger pore volume mesoporous silica was also synthesized to increase amine loading and thus increase CO 2 capacity. The optimized sorbent material exhibited CO 2 capacity as high as 4.88 mmol/g under humid DAC conditions and negligible loss (<1%) during 10 cyclic heating-cooling operations. The optimized PA-coated sorbent also showed fast adsorption and desorption kinetics, with <20% t1/2 increase compared to uncoated sorbent. PEF fabrication conditions, including organic solvents for dissolving core and shell polymers, voltage, distance from the nozzle to the collection panel, etc. were adjusted to better incorporate HNC. After incorporating the optimized sorbent material into PEF, the structured sorbent had a CO 2 capacity of approximately 4.0 mmol/g under humid DAC conditions, with capacity loss of 0.17% per cycle and t1/2 increase less than 10%. A techno-economic analysis (TEA) for the process design for a DAC system based on our developed sorbent structure of trapped small amines in HNC embedded in PEF was conducted. The process design included process description and major equipment sizing and energy and mass balances in addition to scale-up research results and estimated capture cost. Aspen Adsorption Simulator was used to fit the experimentally measured breakthrough curves and extract equilibrium and kinetic data of the optimized sorbent. Our results indicated that for a DAC plant with CO 2 productivity of 3,000 tonne/year, the levelized cost of CO 2 capture was $\$$612/tonne, with the largest contribution of 44.33% from the fixed operation cost. Increasing CO 2 productivity, while maintaining similar fixed operation cost, is expected to significantly reduce the CO 2 capture cost. A sensitivity study was also conducted to understand the influence of total plant cost, sorbent cost, CO 2 concentration in the feed, sorbent mat lifetime, sorbent regeneration electricity, and adsorption blower pressure drop on the levelized cost of CO 2 capture, revealing a capture cost range of $\$$520-870/tonne.

36 MATERIALS SCIENCE↗

HFIR LEU High Density Silicide Dispersion Optimized Design Neutronics Analyses with PHAME

A high-fidelity neutronics model of the Oak Ridge National Laboratory High Flux Isotope Reactor (HFIR) with the low-enriched uranium (LEU) high-density silicide dispersion Optimized fuel design was updated and analyzed to generate reactor physics-based metrics to support follow-on thermal hydraulic and transient analyses of this design. The Python HFIR Analysis and Measurement Engine (PHAME) was also updated to enhance the automation capabilities of the framework developed and maintained to perform these reactor physics modeling and simulation efforts. The automated framework significantly increases the efficiency and reproducibility to design and thoroughly analyzes HFIR LEU core designs, changes, and uncertainties. Reactor physics metrics evaluated include but are not limited to fuel depletion, cycle length, fission rate density distributions, axial power peaking factors, kinetics data, reactivity coefficients, control element worths, heat deposition rates, and decay heat. These neutronics results provide essential input to follow-on steady state thermal, thermal hydraulic and reactor transient analyses, which are subject of other reports. The Optimized design operates at 95 MW to maintain HFIR’s current highly enriched uranium core performance level at 85 MW.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗