Search NASA⌕ Search

SEARCH · Search NASA

Results for “Kinetic Modelling”

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 973 records · Page 54
All sourcesNASA NTRSDOE OSTI
DOE OSTI2018

Gas Phase Chemical Evolution of Uranium, Aluminum, and Iron Oxides

Abstract We use a recently developed plasma-flow reactor to experimentally investigate the formation of oxide nanoparticles from gas phase metal atoms during oxidation, homogeneous nucleation, condensation, and agglomeration processes. Gas phase uranium, aluminum, and iron atoms were cooled from 5000 K to 1000 K over short-time scales (∆t < 30 ms) at atmospheric pressures in the presence of excess oxygen. In-situ emission spectroscopy is used to measure the variation in monoxide/atomic emission intensity ratios as a function of temperature and oxygen fugacity. Condensed oxide nanoparticles are collected inside the reactor for ex-situ analyses using scanning and transmission electron microscopy (SEM, TEM) to determine their structural compositions and sizes. A chemical kinetics model is also developed to describe the gas phase reactions of iron and aluminum metals. The resulting sizes and forms of the crystalline nanoparticles (FeO-wustite, eta-Al 2 O 3 , UO 2 , and alpha-UO 3 ) depend on the thermodynamic properties, kinetically-limited gas phase chemical reactions, and local redox conditions. This work shows the nucleation and growth of metal oxide particles in rapidly-cooling gas is closely coupled to the kinetically-controlled chemical pathways for vapor-phase oxide formation.

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

Mechanism, thermochemistry, and kinetics of the reversible reactions: C 2 H 3 + H 2 ⇌ C 2 H 4 + H ⇌ C 2 H 5

High-level coupled cluster theory, in conjunction with Active Thermochemical Tables (ATcT) and E,J-resolved master equation calculations, was used in a study of the title reactions, which play an important role in the combustion of hydrocarbons. In the set of radical/radical reactions leading to soot formation in flames, the addition of H-atoms to alkenes is likely a common reaction, triggering the isomerization of complex hydrocarbons to aromatics. The heats of formation of C 2 H 3 , C 2 H 4 , and C 2 H 5 are established to be 301.26 ± 0.30 at 0 K (297.22 ± 0.30 at 298 K), 60.89 ± 0.11 (52.38 ± 0.11), and 131.38 ± 0.22 (120.63 ± 0.22) kJ mol -1 , respectively. The calculated rate constants from first principles agree well with experiments where they are available. Under conditions typical of high temperature combustion – where experimental work is very challenging with a consequent dearth of accurate data –here we provide high-level theoretical results for kinetic modeling.

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

Octane Requirements of Lean Mixed-Mode Combustion in a Direct-Injection Spark-Ignition Engine

Here, this study investigates the octane requirements of a hybrid flame propagation and controlled autoignition mode referred to as mixed-mode combustion (MMC), which allows for strong control over combustion parameters via a spark-initiated deflagration phase. Due to the throughput limitations associated with both experiments and 3-D computational fluid dynamics calculations, a hybrid 0-D and 1-D modeling methodology was developed, supported by experimental validation data. This modeling approach relied on 1-D, two-zone engine simulations to predict bulk in-cylinder thermodynamic conditions over a range of engine speeds, compression ratios, intake pressures, trapped residual levels, fueling rates, and spark timings. Those predictions were then transferred to a 0-D chemical kinetic model, which was used to evaluate the autoignition behavior of fuels when subjected to temperature–pressure trajectories of interest. Finally, the predicted autoignition phasings were screened relative to the progress of the modeled deflagration-based combustion in order to determine if an operating condition was feasible or infeasible due to knock or stability limits. The combined modeling and experimental results reveal that MMC has an octane requirement similar to modern stoichiometric spark-ignition engines in that fuels with high research octane number (RON) and high octane sensitivity (S) enable higher loads. Experimental trends with varying RON and S were well predicted by the model for 1000 and 1400 rpm, confirming its utility in identifying the compatibility of a fuel’s autoignition behavior with an engine configuration and operating strategy. However, the model was not effective in predicting (nor designed to predict) operability limits due to cycle-to-cycle variations, which experimentally inhibited operation of some fuels at 2000 rpm. Putting the operable limits and efficiency from MMC in the context of a state-of-the-art engine, the MMC showed superior efficiencies over the range investigated, demonstrating the potential to further improve fuel economy.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2019

Modelling electron transfer in photosystem I: limits and perspectives

Uncovering the parameters underlying the electron transfer (ET) in photosynthetic reaction centres is of importance for understanding the molecular mechanisms underpinning their functionality. The reductive nature of most cofactors involved in photosynthetic ET makes the direct estimation of their properties difficult. Photosystem I (PSI) operates in a highly reducing regime, making the assessment of cofactor properties even more difficult. Kinetic modelling coupled to a non-adiabatic description of ET is a useful approach in overcoming this hindrance. Here we review the theory and modelling approaches that have been used in assessing parameters associated with ET reactions in PSI, with particular attention to ET reactions involving the phylloquinones and the iron–sulphur clusters. In most modelling studies, the goal is to estimate the driving force of ET, which is usually associated with the cofactor midpoint potentials. The driving force is sensitive to many factors, which define the ET rate, i.e. the reorganisation energy, the coupling with nuclear modes and the electronic matrix elements, which are explored and discussed here. The importance of an inclusive modelling of both forward and reverse ET processes is discussed and highlighted. It is shown that although estimates are indeed sensitive to the exact parameter sets employed in the modelling, a general consensus is still attained, pointing to a scenario where $ΔG^0_{A_{tA→Fx}}/ΔG^0_{A_{tB→Fx}}$ is weakly endergonic/exergonic, respectively. It is emphasised that to further refine those estimates, it will require a joint effort between computational modelling and more wide-ranging experimental studies.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

A deep learning approach to programmable RNA switches

Engineered RNA elements are programmable tools capable of detecting small molecules, proteins, and nucleic acids. Predicting the behavior of these synthetic biology components remains a challenge, a situation that could be addressed through enhanced pattern recognition from deep learning. Here, we investigate Deep Neural Networks (DNN) to predict toehold switch function as a canonical riboswitch model in synthetic biology. To facilitate DNN training, we synthesize and characterize in vivo a dataset of 91,534 toehold switches spanning 23 viral genomes and 906 human transcription factors. DNNs trained on nucleotide sequences outperform (R 2 = 0.43–0.70) previous state-of-the-art thermodynamic and kinetic models (R 2 = 0.04–0.15) and allow for human-understandable attention-visualizations (VIS4Map) to identify success and failure modes. This work shows that deep learning approaches can be used for functionality predictions and insight generation in RNA synthetic biology.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2019

Formation of solitary zonal structures via the modulational instability of drift waves

{\rtf1\ansi\ansicpg1252\cocoartf1561\cocoasubrtf600{\fonttbl\f0\fswiss\fcharset0 Helvetica;}{\colortbl;\red255\green255\blue255;\red0\green0\blue0;}{\*\expandedcolortbl;;\cssrgb\c0\c0\c0;}\margl1440\margr1440\vieww10800\viewh8400\viewkind0\pard\tx887\tx1775\tx2662\tx3550\tx4438\tx5325\tx6213\tx7101\tx7988\tx8876\tx9764\tx10651\tx11539\tx12427\tx13314\tx14202\tx15090\tx15977\tx16865\tx17753\tx18640\tx19528\tx20416\tx21303\tx22191\tx23079\tx23966\tx24854\tx25742\tx26629\tx27517\tx28405\tx29292\tx30180\tx31067\tx31955\tx32843\tx33730\tx34618\tx35506\tx36393\tx37281\tx38169\tx39056\tx39944\tx40832\tx41719\tx42607\tx43495\tx44382\tx45270\tx46158\tx47045\tx47933\tx48821\tx49708\tx50596\tx51484\tx52371\tx53259\tx54147\tx55034\tx55922\tx56810\tx57697\tx58585\tx59472\tx60360\tx61248\tx62135\tx63023\tx63911\tx64798\tx65686\tx66574\tx67461\tx68349\tx69237\tx70124\tx71012\tx71900\tx72787\tx73675\tx74563\tx75450\tx76338\tx77226\tx78113\tx79001\tx79889\tx80776\tx81664\tx82552\tx83439\tx84327\tx85215\tx86102\tx86990\tx87877\tx88765\slleading20\pardirnatural\partightenfactor0\f0\fs38 \cf2 The dynamics of the radial envelope of a weak coherent drift wave is approximately governed by a nonlinear Schr\'f6dinger equation, which emerges as a limit of the modified Hasegawa\'97Mima equation. The nonlinear Schr\'f6dinger equation has well-known soliton solutions, and its modulational instability can naturally generate solitary structures. In this paper, we demonstrate that this simple model can adequately describe the formation of solitary zonal structures in the modified Hasegawa\'97Mima equation, but only when the amplitude of the coherent drift wave is relatively small. At larger amplitudes, the modulational instability produces stationary zonal structures instead. Furthermore, we find that incoherent drift waves with beam-like spectra can also be modulationally unstable to the formation of solitary or stationary zonal structures, depending on the beam intensity. Notably, we show that these drift waves can be modeled as quantumlike particles (\'93driftons\'94) within a recently developed phase-space (Wigner\'97Moyal) formulation, which intuitively depicts the solitary zonal structures as quasi-monochromatic drifton condensates. Quantumlike effects, such as diffraction, are essential to these condensates; hence, the latter cannot be described by wave-kinetic models that are based on the ray approximation.\}

Solitons↗
NASA NTRS2003

Production, Outflow Velocity, and Radial Distribution of H2O and OH in the Coma of Comet C/1995 O1 [Hale-Bopp] from Wide Field Imaging of OH

Observations of OH are a useful proxy of the water production rate (Q(sub H2O)) and outflow velocity (V(sub out)) in comets. From wide field images taken on 03/28/1997 and 04/08/1997 that capture the entire scale length of the OH coma of comet C/1995 O1 (Hale-Bopp), we obtain Q(sub H2O) from the model-independent method of aperture summation. With an adaptive ring summation algorithm, we extract the radial brightness distribution of OH 0-0 band emission out to cometocentric distances of up to 10(exp 6) km, both as azimuthal averages and in quadrants covering different position angles relative to the comet-Sun line. These profiles are fit using both fixed and variable velocity 2-component spherical expansion models to estimate V(sub OH) with increasing distance from the nucleus. The OH coma of Hale-Bopp was more spatially extended than previous comets, and this extension is best matched by a variable acceleration of H2O and OH that acted across the entire coma, but was strongest within 1-2 x 10(exp 4) km from the nucleus. Our models indicate that V(sub OH) at the edge of our detectable field of view (10(exp 6) km) was approx. 2-3 times greater in Hale-Bopp than for a 1P/Halley-class comet at 1 AU, which is consistent with the results of more sophisticated gas-kinetic models, extrapolation from previous observations of OH in comets with Q(sub H2O) greater than 10(exp 29)/s , and direct radio measurements of the outer coma Hale-Bopp OH velocity. The most probable source of this acceleration is thermalization of the excess energy of dissociation of H2O and OH over an extended collisional coma. When the coma is broken down by quadrants in position angle, we find an azimuthal asymmetry in the radial distribution that is characterized by an increase in the spatial extent of OH in the region between the orbit-trailing and anti-sunward directions. Model fits specific to this area and comparison with radio OH measurements suggest greater acceleration here, with V(sub OH) approx. 1.5 times greater at a 10(exp 6) km cometocentric distance than elsewhere in the coma. We discuss several mechanisms that may have acted within the coma to produce the observed effect.

Harris, Walter M.↗
DOE OSTI2022

BETO 2021 Peer Review - 3.4.2.302 - Process Scale Up for Production Environments

This project supports the Department of Energy's Bioenergy Technology Office mission of transitioning bioenergy technologies to market by de-risking integration, developing scaling relations and modeling that de-risks scale up, and providing a feedback loop with lower TRL projects to understand process and technology fundamentals. Two objectives were pursued during this merit cycle. First, the Thermal and Catalytic Process Development Unit (TCPDU), a 0.5 ton/day pilot plant, was used to validate kinetic models developed by the Consortium for Computational Physics and Chemistry (CCPC). Validating the models accurately predict product yields and composition is a first step toward de-risking a common industry failure - taking too large of steps between scales. In addition, these models help with design and troubleshooting of new unit operations for the TCPDU. The second objective was to conduct the FY22 verification campaign around ex situ CFP. The primary challenge of this effort was the coordination of multiple national labs and projects. This was overcome by a dedicated leadership role and effective communication strategy between the projects. Currently the project is undergoing a pivot based on a stagegate decision to not conduct the verification campaign in the TCPDU and it is focused on closing out the CFP technology at a smaller scale by the end of the current fiscal year.

BIOMASS FUELS↗
DOE OSTI2018

Developing the test conditions for the first fueled transients in the TREAT reactor

The restart of the Transient Reactor Test (TREAT) facility at Idaho National Laboratory, re-establishes a national in-pile testing capability for nuclear fuels in accident scenarios. The purpose of the first tests in TREAT are to study the response of nuclear fuel durring a reactivity initiated accident (RIA). The initial objective is to demonstrate the ability to create a desired enthalpy increase in a representative fuel test specimen over a desired time step and to better understand the uncertainty. The five test rods initially available for these irradiations are UO2-Zircaloy with standard Pressurized Water Reactor (PWR) geometry. The rods are 15 cm long with 9 cm of active fuel. The initial testing environment is dry helium. BISON fuel performance calculations were run to model the heat transfer during postulated transients. Based on those results energy deposition targets were selected around the range of observable phenomena in the testing environment. Reactor / experiment power coupling factors, step insertion size and transient clipping response times were then calculated using results from a RELAP point kinetics model. Calculated energy deposition targets are 150, 350, 650, 960, and 1300 J/gUO2. The energy coupling factor for the experiments was calculated to be ~1.5 J/gUO2-MJ. Two step insertion sizes were selected: a 0.6% ?k/k step for the 150 J/g transient and a 1.8% ?k/k step for the other transients. Clipping response times ranged from 1.3 seconds to unclipped transients. The experiments will be irradiated in the TREAT facility in the fall of 2018. The first transient will be the 150 J/g pulse with a 0.6% ?k/k step that is unclipped. Follow on tests may modify the reactor operating parameters based on results from previous transients. Any preliminary data available on the results of the transients at the time of the NUMAT Conference will be presented.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS2022

An ODE-Fitting Approach to Estimate Critical Aircraft Performance Parameters for Trajectory Prediction

Ground-based decision support tools (DST) in air traffic management (ATM) typically perform trajectory prediction based on aircraft performance model (APM) parameters, but some or all of these parameters might not be readily available. In particular, the three critical parameters required for trajectory prediction are the thrust setting, drag coefficients, and takeoff weight of the aircraft. Unfortunately, these parameters are coupled and appear together in physics-based kinetic models. Past approaches utilize data from a specific phase of flight (climb, level flight or descent), where one or more of the parameters are assumed to be known and estimate the remaining unknown parameters. This approach introduces bias/errors and also does not extend to scenarios where all of the above parameters are not known with sufficient accuracy. This paper is the first of its kind to propose a generalized framework for simultaneous estimation of all three critical APM parameters (thrust, drag, and mass). The proposed approach utilizes data from both the climb and descent phases and fits the ordinary differential equation for altitude in each phase using historical trajectory data available from radar tracks or ADS-B. The approach yields a set of optimized APM parameters that are best suited to fit each historical flight record. The methodology is applied on on sample flights from three different aircraft types, and the results demonstrate low fit error and consequently will yield a high level of prediction accuracy.

trajectory prediction, aircraft performance model,↗
NASA NTRSDate not supplied

An ODE-Fitting Approach to Estimate Critical Aircraft Performance Parameters for Trajectory Prediction

Ground-based decision support tools (DST) in air traffic management (ATM) typically perform trajectory prediction based on aircraft performance model (APM) parameters, but some or all of these parameters might not be readily available. In particular, the three critical parameters required for trajectory prediction are the thrust setting, drag coefficients, and takeoff weight of the aircraft. Unfortunately, these parameters are coupled and appear together in physics-based kinetic models. Past approaches utilize data from a specific phase of flight (climb, level flight or descent), where one or more of the parameters are assumed to be known and estimate the remaining unknown parameters. This approach introduces bias/errors and also does not extend to scenarios where all of the above parameters are not known with sufficient accuracy. This paper is the first of its kind to propose a generalized framework for simultaneous estimation of all three critical APM parameters (thrust, drag, and mass). The proposed approach utilizes data from both the climb and descent phases and fits the ordinary differential equation for altitude in each phase using historical trajectory data available from radar tracks or ADS-B. The approach yields a set of optimized APM parameters that are best suited to fit each historical flight record. The methodology is applied on on sample flights from three different aircraft types, and the results demonstrate low fit error and consequently will yield a high level of prediction accuracy.

trajectory prediction↗
DOE OSTI2024

A Detailed Reaction Mechanism for Thiosulfate Oxidation by Ozone in Aqueous Environments

The ozone oxidation, or ozonation, of thiosulfate is an important reaction for wastewater processing, where it is used for remediation of mining effluents, and for studying aerosol chemistry, where its fast reaction rate makes it an excellent model reaction. Although thiosulfate ozonation has been studied since the 1950s, challenges remain in developing a realistic reaction mechanism that can satisfactorily account for all observed products with a sequence of elementary reaction steps. Here, we present novel measurements using trapped microdroplets to study the pH-dependent thiosulfate ozonation kinetics. We detect known products and intermediates, including SO32-, SO42-, S3O62-, and S4O62-, establishing agreement with the literature. However, we identify S2O42- as a new reaction intermediate and find that the currently accepted mechanism does not directly explain observed pH effects. Thus, we develop a new mechanism, which incorporates S2O42- as an intermediate and uses elementary steps to explain the pH dependence of thiosulfate ozonation. The proposed mechanism is tested using a kinetic model benchmarked to the experiments presented here, then compared to literature data. We demonstrate good agreement between the proposed thiosulfate ozonation mechanism and experiments, suggesting that the insights in this paper can be leveraged in wastewater treatment and in understanding potential climate impacts.

Deal, Alexandra M↗
DOE OSTI2024

Direct Evidence for a Sequential Electron Transfer–Proton Transfer Mechanism in the PCET Reduction of a Metal Hydroxide Catalyst

Here, the proton-coupled electron transfer (PCET) mechanism for the reaction M ox -OH + e − + H + → M red -OH 2 was de-termined through kinetic resolution of the independent electron transfer (ET) and proton transfer (PT) steps. The reaction of interest was triggered by visible light excitation of [Ru II (tpy)(bpy′)H 2 O] 2+ , Ru II -OH 2 , where tpy is 2,2′:6′,2″-terpyridine and bpy′ is 4,4′-diaminopropylsilatrane-2,2′-bipyridine, anchored to In 2 O 3 :Sn (ITO) thin films in aqueous solutions. Interfacial kinetics for the PCET reduction reaction were quantified by nanosecond transient absorption spectroscopy as a function of solution pH and applied potential. Data acquired from pH = 5-10 revealed a stepwise electron-transfer proton-transfer (ET-PT) mechanism, while kinetic measurements made below the pK a (Ru III -OH/OH 2 ) = 1.3 were used to study the analogous interfacial reaction where electron transfer was the only mechanistic step. Analysis of this data with a recently reported multi-channel kinetic model was used to construct a PCET zone diagram and supported the assignment of an ET-PT mechanism from pH = 5-10. Ultimately, this study represents a unique example amongst M ox -OH/M red -OH 2 reactivity where the protona-tion and oxidation state of the intermediate was kinetically and spectrally resolved to firmly establish the PCET mechanism.

Charge transfer↗
DOE OSTI2021

Magnetic Energy Release, Plasma Dynamics, and Particle Acceleration in Relativistic Turbulent Magnetic Reconnection

In strongly magnetized astrophysical plasma systems, magnetic reconnection is believed to be the primary process during which explosive energy release and particle acceleration occur, leading to significant high-energy emission. Past years have witnessed active development of kinetic modeling of relativistic magnetic reconnection, supporting this magnetically dominated scenario. A much less explored issue in studies of relativistic reconnection is the consequence of three-dimensional dynamics, where turbulent structures are naturally generated as various types of instabilities develop. This paper presents a series of three-dimensional, fully kinetic simulations of relativistic turbulent magnetic reconnection (RTMR) in positron–electron plasmas with system domains much larger than kinetic scales. Our simulations start from a force-free current sheet with several different modes of long-wavelength magnetic field perturbations, which drive additional turbulence in the reconnection region. Because of this, the current layer breaks up and the reconnection region quickly evolves into a turbulent layer filled with coherent structures such as flux ropes and current sheets. Additionally, we find that plasma dynamics in RTMR is vastly different from its 2D counterpart in many aspects. The flux ropes evolve rapidly after their generation, and can be completely disrupted by the secondary kink instability. This turbulent evolution leads to superdiffusive behavior of magnetic field lines as seen in MHD studies of turbulent reconnection. Meanwhile, nonthermal particle acceleration and the timescale for energy release can be very fast and do not depend strongly on the turbulence amplitude. Additionally, the main acceleration mechanism is a Fermi-like acceleration process supported by the motional electric field, whereas the nonideal electric field acceleration plays a subdominant role. We also discuss possible observational implications of three-dimensional RTMR in high-energy astrophysics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Methanol Decomposition on Copper Surfaces under Ambient Conditions: Mechanism, Surface Kinetics, and Structure Sensitivity

Here, we study the adsorption of methanol vapor under ambient pressure and temperature conditions on low-index Cu surfaces using surface-sensitive infrared (IR) and X-ray spectroscopy techniques. The first step of methanol decomposition, i.e. , breaking of the O—H bond to form surface-bound methoxy, readily occurs under ambient conditions. Time-lapse IR spectra clearly indicate a gradually decreasing methoxy coverage, which does not obey well established kinetic models. We rationalize the initial temperature-independent, high, nonequilibrium coverage of methoxy by a H-bonded methanol assembly in the precursor state. A temperature-dependent equilibrium coverage is achieved as the excess methoxy is eliminated gradually via further dehydrogenation to CO that desorbs to the gas phase. The kinetics of this process displays a significant structure sensitivity with considerably faster kinetics on the Cu(110) surface compared to Cu(111) and Cu(100) surfaces.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2019

Formation of solitary zonal structures via the modulational instability of drift waves

The dynamics of the radial envelope of a weak coherent drift-wave (DW) is approximately governed by a nonlinear Schrödinger equation (NLSE), which emerges as a limit of the modified Hasegawa–Mima equation (mHME). The NLSE has well-known soliton solutions, and its modulational instability (MI) can naturally generate solitary structures. Here, we demonstrate that this simple model can adequately describe the formation of solitary zonal structures in the mHME, but only when the amplitude of the coherent DW is relatively small. At larger amplitudes, the MI produces stationary zonal structures instead. Furthermore, we find that incoherent DWs with beam-like spectra can also be modulationally unstable to the formation of solitary or stationary zonal structures, depending on the beam intensity. Notably, we show that these DWs can be modeled as quantum-like particles ('driftons') within a recently developed phase-space (Wigner–Moyal) formulation, which intuitively depicts the solitary zonal structures as quasi-monochromatic drifton condensates. Quantum-like effects, such as diffraction, are essential to these condensates; thus, the latter cannot be described by wave-kinetic models that are based on the ray approximation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Integrating Carbon Capture, Utilization, & Sequestration into Chemical Pulp Mills

The U.S. pulp and paper industry presents a unique and largely untapped opportunity for large- scale carbon dioxide removal (CDR). Unlike most industrial sectors, pulp mills rely heavily on biomass, meaning that much of their carbon emissions originate from atmospheric CO₂ that was recently captured by plants. If this biogenic CO₂ can be captured and permanently stored, pulp mills can be transformed from carbon emitters into net carbon removal facilities. This project was motivated by that opportunity and aimed to develop and evaluate integrated, low-cost strategies for capturing, utilizing, and sequestering CO₂ within existing chemical pulping operations. The scope of this work focused on four complementary innovations designed to integrate seamlessly into kraft pulp mill infrastructure: (1) in situ CO₂ capture within the recovery cycle, (2) oxy-fuel retrofitting of the rotary lime kiln to produce a high-purity CO₂ stream, (3) ex situ CO₂ capture and mineralization using pulp mill residues (dregs, grits, and lime mud), and (4) beneficial reuse of these residues as mineral carbonate fertilizers. The project combined process modeling, laboratory experimentation, life cycle assessment (LCA), and field trials to evaluate the technical feasibility, economic viability, and environmental impact of these approaches. The results demonstrate that pulp mills can serve as effective platforms for carbon removal when equipped with integrated carbon capture systems. Process modeling showed that combining sodium spiking with oxy-fuel calcination significantly enhances CO₂ capture efficiency while reducing costs by up to 31% compared to conventional configurations. Experimental work further revealed that calcination behavior in high-CO₂ environments differs substantially from traditional systems, leading to the development of a new kinetic model that predicts reaction rates under these conditions. This model provides essential design guidance for next-generation decarbonized lime kilns. In parallel, the project demonstrated that alkaline mineral residues generated during pulping operations can be repurposed as a sustainable alternative to agricultural lime. Across a wide range of soils in the southeastern United States, these materials performed equivalently to commercial lime in adjusting soil pH while offering lower greenhouse gas emissions and reduced cost. Field and greenhouse studies confirmed that crop and tree growth responses were comparable, supporting their viability as a drop-in replacement. This co-product pathway provides a practical utilization strategy that offsets costs and improves overall system economics. A major contribution of this project is the first comprehensive life cycle assessment of carbon removal in pulp and paper systems across multiple system boundaries. Results show that retrofitted mills can achieve carbon removal efficiencies ranging from 12% to 92%, depending on how the system is defined. This finding highlights a critical issue in carbon accounting: reported performance is highly sensitive to methodological choices. By explicitly quantifying these differences, this work provides valuable guidance for policymakers, carbon registries, and project developers working to standardize carbon removal metrics. From a commercialization perspective, the technologies investigated in this project are well- aligned with existing industrial infrastructure, minimizing the need for entirely new facilities. 3 DE-EE0009413 Industry engagement throughout the project—including collaboration with pulp and paper companies, equipment manufacturers, and carbon removal developers—has accelerated the transition from research to deployment. Notably, a commercial developer is actively pursuing carbon capture projects at pulp mills in the southeastern United States and has cited this research as a contributing foundation. The emergence of voluntary carbon markets and long-term offtake agreements further strengthens the business case for implementation. The broader public benefits of this work are significant. By enabling large-scale carbon removal using existing industrial systems, this approach offers a near-term pathway to reduce atmospheric CO₂ concentrations while supporting domestic manufacturing and rural economies. The reuse of industrial residues as fertilizers reduces reliance on mined materials, lowers costs for farmers, and decreases environmental impacts associated with conventional lime production. In addition, the project has supported workforce development by training graduate students and researchers in carbon capture technologies, helping to build capacity in a critical area of national interest. In conclusion, this project demonstrates that integrated carbon capture, utilization, and sequestration in pulp mills is both technically feasible and economically promising. By combining process innovation, experimental validation, and systems-level analysis, the work advances the understanding of how biomass-based industries can contribute to climate mitigation. The findings provide a strong foundation for commercial deployment and offer a scalable solution for transforming a major U.S. industry into a source of durable carbon removal.

09 BIOMASS FUELS↗
NASA NTRS1992

Modeling of hypersonic reacting flows

The effect of nonequilibrium thermochemical processes on the aerodynamic and heat transfer characteristics of a hypersonic vehicle is illustrated with examples. It is shown that the conventional method of predicting chemical reactions always predict that the flow is closer to equilibrium than it actually is, leading to incorrect predictions of aerodynamic characteristics of a vehicle. A method is presented for predicting chemical processes using a two-temperature kinetic model, which is developed on the basis of combining conservation equations for vibrational energy and for electron electronic energy, which are also derived. The validity of the two-temperature model is proven using two types of experimental data: radiation and shock shapes.

Park, Chul↗