Search NASA⌕ Search

SEARCH · Search NASA

Results for “Flow rate ratio”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Comparative investigation of Ga- and In-CHA in the non-oxidative ethane dehydrogenation reaction

Ga- and In-exchanged chabazite (CHA) zeolites with same Si/Al and metal/Al ratios were prepared via the incipient wetness impregnation method, were characterized using N 2 adsorption, electron microscopy, temperature-programed reactions and were evaluated for the ethane dehydrogenation reaction using flow microreactors. Ga-CHA has higher reaction rates and a lower activation energy of 107 kJ/mol than In-CHA (E a = 175 kJ/mol). Rietveld refinement of the X-ray powder diffraction pattern shows that the In + cation is predominantly located above the 6-ring of the CHA cage. It is proposed that the reaction proceeds through the alkyl mechanism based on stability of alkyl hydride intermediates as determined using DFT calculations. The oxidative addition of ethane to the metal shows much lower Gibbs free energy for Ga-CHA (+27.95 kJ/mol) vs In-CHA (+124.85 kJ/mol). Finally, these results indicate that oxidative addition may be the rate-limiting step of ethane dehydrogenation in these materials.

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

Creation of a Mixed-Mode Fracture Network at Mesoscale Through Hydraulic Fracturing and Shear Stimulation

Enhanced Geothermal Systems could provide a substantial contribution to the global energy demand if their implementation could overcome inherent challenges. Examples are insufficient created permeability, early thermal breakthrough, and unacceptable induced seismicity. In this study we report on the seismic response of a mesoscale hydraulic fracturing experiment performed at 1.5-km depth at the Sanford Underground Research Facility. We have measured the seismic activity by utilizing a 100-kHz, continuous seismic monitoring system deployed in six 60-m length monitoring boreholes surrounding the experimental domain in 3-D. The achieved location uncertainty was on the order of 1 m and limited by the signal-to-noise ratio of detected events. These uncertainties were corroborated by detections of fracture intersections at the monitoring boreholes. Three intervals of the dedicated injection borehole were hydraulically stimulated by water injection at pressures up to 33 MPa and flow rates up to 5 L/min. We located 1,933 seismic events during several injection periods. The recorded seismicity delineates a complex fracture network comprised of multistrand hydraulic fractures and shear-reactivated, preexisting planes of weakness that grew unilaterally from the point of initiation. We find that heterogeneity of stress dictates the seismic outcome of hydraulic stimulations, even when relying on theoretically well-behaved hydraulic fractures. Once hydraulic fractures intersected boreholes, the boreholes acted as a pressure relief and fracture propagation ceased. In order to create an efficient subsurface heat exchanger, production boreholes should not be drilled before the end of hydraulic stimulations.

58 GEOSCIENCES↗
DOE OSTI2024

Understanding and modeling current and future coastal wetland methane dynamics (Final Report)

The coastal terrestrial-aquatic interface (TAI) is a highly dynamic component of the Earth system that plays a critical role in biogeochemical cycling. Due to its dynamic nature, the processes that regulate decomposition and methane (CH 4 ) emissions are of greater significance at the TAI than in upland systems. Despite this, we have limited mechanistic understanding of how climate stressors interact to regulate the electron acceptors and donors that determine decomposition pathways within TAIs, including the generation of hot spots and hot moments. Accurately modeling these processes is critical for incorporating the coastal TAI into Earth systems models, such as DOE’s E3SM. With previous DOE support, we adapted an aerobic terrestrial representation of decomposition using PFLOTRAN, a reactive flow and transport model, and added anerobic decomposition pathways, salinity, and oxygen (O 2 ) that fluctuates independently of water table level. However, because this model (PFLOTRAN TAI ) is based on decomposition rates and organic matter carbon to nitrogen ratios from terrestrial systems, its performance in TAI systems is currently limited by the lack of empirical data to properly parameterize variables. In addition, while PFLOTRAN TAI can simulate movement of O 2 into sediments, it is not currently capable of tracking the movement of CH 4 gas through plant tissues due to both current model structure and lack of available data. Our overall objective of the project was the increase our mechanistic understanding of CH 4 dynamics in response to environmental change, such that we can improve the representation of these dynamics in PFLOTRAN. We installed automated flux chambers in a new field-scale active soil warming experiment and set up marsh organs (mesocosms) to test effects of warming, flooding, and salinity. This resulted in a new dataset consisting of chamber-level CH 4 flux measurements across multiple sites, ecological conditions, and timeframes, as well as corresponding measurements on porewater chemistry, soil carbon quality, plant biomass, and redox reaction rates. Using these data, we improved PFLOTRAN TAI to more accurately model CH 4 dynamics and successfully tested our hypotheses. This grant contributed to the professional development of 3 postdocs, 4 undergraduate interns, 4 teacher interns, 17 technicians, and 15 participatory scientists. The automated chamber technology designed for this grant has also been shared with multiple new projects.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

The influence of temperature history and flow mixing on the vapor-phase speciation of uranium oxide nanoparticles

The transport of radioactive material following a nuclear event depends heavily on fallout characteristics (e.g., composition, shape, and size), which in turn are controlled by the conditions under which the fallout forms. To this end, we use a novel plasma flow reactor (PFR) to explore how conditions including temperature history and gas mixing influence the chemical speciation and physical characteristics of uranium oxide nanoparticles. The PFR consists of a glass tube attached to an inductively coupled plasma (ICP) torch via an adaptor piece (ring flow injector) located downstream of the plasma-generating RF coil. We increase the gas flow rates from normal ICP operating conditions to higher values both upstream at the RF coil and downstream at the ring flow injector to create four distinct temperature profiles and two distinct gas flow mixing conditions (normal and higher Ar atomic ratios). Ex situ transmission electron microscopy (TEM) analysis is used to determine the crystal structure, morphology, and size distribution of the synthesized particles as a function of temperature and mixing. The particles display a striking dependence on these conditions in terms of the resulting speciation of uranium oxide: lower temperature (higher upstream) and slower cooling (normal downstream with normal Ar atomic ratio mixing) favor the formation of α-UO 3 particles, while faster cooling (higher downstream with higher Ar atomic ratio mixing) and higher temperature (normal upstream) favor the formation of UO 2 particles. The extent of gas mixing dictates particle aggregation: better mixing encourages the particles to aggregate more extensively. Here these results demonstrate the sensitivity of uranium to initial local conditions and can be used to improve our understanding of the kinetics that drive uranium speciation and ultimately inform fallout transport models.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Assessing hybrid vigour using the thermal sensitivity of physiological trade-offs in tiger salamanders

Hybridization between species affects biodiversity and population sustainability in numerous ways, many of which depend on the fitness of the hybrid relative to the parental species. Hybrids can exhibit fitter phenotypes compared to the parental lineages, and this ‘hybrid vigour’ can then lead to the extinction of one or both parental lines. In this study, we analysed the relationship between water loss and gas exchange to compare physiological performance among three tiger salamander genotypes—the native California tiger salamander (CTS), the invasive barred tiger salamanders (BTS) and CTS × BTS hybrids across multiple temperatures (13.5°C, 20.5°C and 23.5°C). We developed a new index of performance, the water-gas exchange ratio (WGER), which we define as the ratio of gas exchange to evaporative water loss (μL V O 2 /μL H 2 O). The ratio describes the ability of an organism to support energetically costly activities with high levels of gas exchange while simultaneously limiting water loss to lower desiccation risk. We used flow through respirometry to measure the thermal sensitivity of metabolic rate and resistance to water loss of each salamander genotype to compare indices of physiological performance. We found that temperature had a significant effect on metabolic rate and resistance to water loss, with both traits increasing as temperatures warmed. Across genotypes, we found that hybrids have a higher WGER than the native CTS, owing to a higher metabolic rate despite having a lower resistance to water loss. These results provide a greater insight into the physiological mechanisms driving hybrid vigour and offer a potential explanation for the rapid spread of salamander hybrids. More broadly, our introduction of the WGER may allow for species- or lineage-wide comparisons of physiological performance across changing environmental conditions, highlighting the insight that can be gleaned from multitrait analysis of organism performance.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Indium as a surfactant: Effects on growth morphology and background impurity in GaN films grown by ammonia-assisted molecular beam epitaxy

We report on the improvement of the surface morphology of c-plane GaN films grown at high growth rates (~1 µm/h) using ammonia molecular beam epitaxy through a series of growth optimizations as well as the introduction of indium as a surfactant. The indium surfactant was expected to help with the adatom mobility and, thus, provide smoother growth surfaces. Through a combination of varying V/III ratios, In flux, and growth temperatures, an optimal condition for surface morphology, characterized by atomic force microscopy, was achieved. At higher Ga fluxes for fast growth rates (~1 µm/h and beam equivalent pressures of ~5 × 10 –7 Torr), higher ammonia flows were necessary to preserve the surface morphology. In addition, indium was an effective surfactant—reducing the roughness and improving the overall surface morphology. However, excessive indium causes the surface morphology to degrade, potentially due to the enhancement of the Ga desorption from the surface as a result of the reaction of indium with ammonia for high indium fluxes. The indium surfactant also resulted in a reduction of background Si impurity concentrations in the film. These effects allow for the growth of thick drift layers with low background dopant concentrations for vertical GaN power devices.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Universal behavior in fragmenting brittle, isotropic solids across material properties

A bonded particle model is used to explore how variations in the material properties of brittle, isotropic solids affect critical behavior in fragmentation. To control material properties, a model is proposed which includes breakable two- and three-body particle interactions to calibrate elastic moduli and mode I and mode II fracture toughnesses. In the quasistatic limit, fragmentation leads to a power-law distribution of grain sizes which is truncated at a maximum grain mass that grows as a nontrivial power of system size. In the high-rate limit, truncation occurs at a mass that decreases as a power of increasing rate. A scaling description is used to characterize this behavior by collapsing the mean-square grain mass across rates and system sizes. Consistent scaling persists across all material properties studied, although there are differences in the evolution of grain size distributions with strain as the initial number of grains at fracture and their subsequent rate of production depend on Poisson's ratio. Importantly, this evolving granular structure is found to induce a unique rheology where the ratio of the shear stress to pressure, an internal friction coefficient, decays approximately as the logarithm of increasing strain rate. The stress ratio also decreases at all rates with increasing strain as fragmentation progresses and depends on elastic properties of the solid.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Refining a Novel Process Monitoring Method to Safeguard Continuously Cycling Designs Using Isotopic Ratios

In advanced reactor (AR) designs, a common feature is continuous chemical processing and circulation of the nuclear material. This work bridges a significant measurement gap in safeguarding reactors with circulating fuel or continuous refueling by leveraging and building on the isotope ratio method first developed by our team under an FY21 Advanced Reactors International Safeguards Engagement (ARISE) project (Uribe et al. 2021). In circulating fuel designs, the radioisotope inventory changes from traditional effects (e.g., radioactive decay, fission) but also includes material transport due to pressure and temperature gradients. Such designs may also require regular or continuous additions or removals during operation, which significantly increases the rate of inventory change compared to traditional pressurized water reactor (PWR)s. Thus, directly tracking the nuclear inventory is ineffective since the isotopes are continuously added and removed. The isotope ratio method instead focuses on detecting changes to the input and output flows of radioisotopes. Previous work showed that for well-chosen pairs of isotopes, the isotopic ratio provides a sensitive and lasting indicator of deviation from normal conditions (e.g., startup, shutdown, diversion). The isotope ratio method is a process monitoring method with potential for application in for forward-looking approaches to International Atomic Energy Agency (IAEA) safeguards. The original process monitoring method was developed for a specific case—the decay tank of a thorium-fueled molten salt breeder reactor. In this expanded work, we explored other types of reactors and processes with nonstationary (e.g., flowing) nuclear material, which are difficult to safeguard with traditional methods because of the transient nature of the systems. The goal of the present work is to generalize the isotope ratio method for use in processes with continuously flowing nuclear material. All continuous processes have an average time for isotopes to be replaced in the system. The isotope ratio method works by choosing radioisotopes with half-lives both above and below the processing time. The present work seeks to explore which isotopes are suitable for the method by simulating the nuclear inventory, radioisotope emissions, and detector responses for several classes of advanced reactors. While the method can in principle be applied to other processes (e.g., enrichment or reprocessing facilities), the present work limits scope to ARs with continuously flowing fuel. Section 2 details the mathematics supporting the isotope ratio method, and Section 3 introduces the representative ARs selected for this work. Section 4 discusses how each reactor was analyzed, and Section 5 showcases the results for each representative reactor. Finally, Section 6 provides concluding remarks and suggests pathways for further analysis.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2026

Mechanistic Modeling of TEG Dehydrator Emissions in Oil and Gas Industry

This work presents a mechanistic modeling approach for simulating methane emissions from triethylene glycol (TEG) dehydrators used in oil & gas (O&G) operations. The model was developed as a modular component of the Mechanistic Air Emissions Simulator (MAES) tool, incorporating species-specific absorption and emission dynamics through two-level, second-order polynomial regression (PR) models trained on ProMax simulation data: (1) species-level regression models that track the transfer rates of individual gas species within the dehydrator unit streams, and (2) outlet flow stream regression models that predict the fraction of inlet gas distributed among the outlet streams of the dehydrator unit. These behaviors were characterized over a range of glycol circulation ratios, wet gas pressures, and temperatures. The model was validated using root mean square error (RMSE) analysis. The species-level PR achieved low root mean square error (RMSE) values (<0.03) for light hydrocarbon species across all dehydrator components, ranging from 0.0009 for methane to 0.029 for normal pentane. Similarly, the outlet-level PR yielded RMSE values below 0.002 for the dry gas fraction, 0.001 for the flash tank fraction, and 0.002 for the still vent fraction, demonstrating strong agreement between predicted and reference ProMax values. When deployed at field facilities, the model significantly improved MAES-simulated dehydrator emissions, revealing that gas-assisted glycol pump emissions are the dominant contributors to both dehydrator-level and site-level methane emissions under uncontrolled conditions. Further analysis of the 154 dehydrator units reported by operators under the AMI 2024 project showed that 54 units (31%) used gas-driven glycol pumps, of which 6 units (11%) operated with uncontrolled flash tanks, and 22 units (40.7%) were identified as potentially oversized. Of the six dehydrator units with uncontrolled gas-assisted pumps, pump emissions accounted for 90.25% of total dehydrator emissions and 63.10% of total site-level emissions. These findings highlight substantial opportunities for emissions mitigation through equipment upgrades.

MAES↗
DOE OSTI2025

Sedimentation and shear-induced dynamics of spheroids in fluids with spatial viscosity variations

A generalized reciprocal theorem is used to relate the force and torque induced on a particle in an inertia-less fluid with small variation in viscosity to integrals involving Stokes flow fields and the spatial dependence of viscosity. These resistivity expressions are analytically evaluated using spheroidal harmonics and then used to obtain the mobility of the spheroid during sedimentation, and in linear flows, of a fluid with linear viscosity stratification. The coupling between the rotational and translational motion induced by stratification rotates the spheroid’s centerline, creating a variety of rotational and translational dynamics dependent upon the particle’s aspect ratio, κ, and the component of the stratification unit vector in the gravity direction, d g . Spheroids with 0.55 ⪅ κ ⪅ 2.0 exhibit the largest variety of settling behaviors. Interestingly, this range covers most microplastics and typical microorganisms. One of the modes include a stable orientation dependent only on κ and d g , but independent of initial orientation, thus allowing for the potential control of settling angles and sedimentation rates. In a simple shear flow, cross-streamline migration occurs due to the stratification-induced force generated on the particle. Similarly, a particle no longer stays at the stagnation point of a uniaxial extensional flow. While fully analytical results are obtained for spheroids, numerical simulations provide a source of validation. These simulations also provide additional insights into the stratification-induced force- and torque-producing mechanisms through the stratification-induced stress, which is not accessed in the reciprocal theorem-based analytical calculations.

Geophysical and Geological Flows: Stratified flows↗
DOE OSTI2020

Adaptive simulations enable computational design of electron beam processing of nanomaterials with supersonic micro-jet precursor

Focused Electron Beam Induced Processing (FEBIP) is a powerful tool for the “direct-write” of nanomaterials with the possibility of atomistic control on suspended 2D material substrates. FEBIP capabilities have been significantly expanded by using a localized jet-based delivery of precursors, and especially when a thermally energized supersonic micro-jet enhances the delivery of mass flux along with controlling the far-from-equilibrium thermodynamic state of adsorbed adatoms. The possibilities of growing nanomaterials with “dialed-in” composition with ultra-high growth rates and aspect ratios up to 100:1 have been demonstrated using the supersonic micro-jet-FEBIP. Bringing this scientific discovery to the level of maturity required for practical applications in additive nanomanufacturing requires simulation tools that are capable of capturing the complex flow physics of micro-jet-substrate interactions that bridge a wide range of flow regimes from the high-density gas micro-jet expanding into the vacuum environment of FEBIP. To address this significant computational challenge, a new approach has been developed and described in this work for a multiscale adaptive DSMC (Direct Simulation Monte Carlo) algorithm to predict the micro-jet gas dynamics in an FEBIP environment, spanning the full range of flow regimes from low Knudsen (Kn) number O(0.01) continuum flow to high Kn of O(10) for the molecular flow within a unified computational framework. Here, the fundamental principles of the adaptive DSMC algorithm are described, its viability as a DSMC technique is demonstrated, and computational improvements for the benchmark cases are discussed in the context of advancing 3D nanofabrication with FEBIP. Ultimately, combining the first principle simulations via adaptive DSMC with the complementary experimental data will enable the creation of the powerful CAD tools for in silico design and optimal operation of micro-jet-FEBIP.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Optimizing the fuel efficiency of an opposed piston engine for electric power generation

This paper investigates the optimal crankshaft motion for an opposed piston (OP) engine in a novel hybrid architecture to maximize fuel efficiency. The OP engine was selected for this work due to its inherent thermodynamic benefits and the balanced nature of the engine which can achieve downsizing through reducing the number of cylinders rather than the individual cylinder volume. The typical geartrain required on an OP engine was exchanged for two electric motors, reducing friction loss and decoupling the crankshafts. Using the motors to control the crankshaft motion profiles, this architecture introduces capabilities to dynamically vary compression ratio, combustion volume, and scavenging dynamics. To leverage these opportunities, an optimization scheme was developed utilizing nonlinear optimization of a 0-D model to compute the crankshaft motion profile that maximizes the work generated by the system. This optimization was then iteratively coupled with a high fidelity model which supplies the cylinder flow boundary conditions. This iterative approach reduces the model complexity used in the optimal control problem (OCP) while capturing the gas exchange dynamics critical to the 2-stroke cycle of the OP engine. By using the rate of change of motor torque as the input to the OCP, the torque fluctuation in a single cycle can be limited to ensure tracking feasibility. The results show crankshaft velocity slows during the compression stroke and conversely accelerates during the expansion stroke, reducing the peak motor torque required for control and thus reducing the motor losses. The extended residence time at top dead center, however, leads to an increase in heat transfer, illustrating the trade-off between the work extraction efficiency and the indicated engine efficiency.

Engineering↗
DOE OSTI2020

Effect of inorganic-to-organic mass ratio on the heterogeneous OH reaction rates of erythritol: implications for atmospheric chemical stability of 2-methyltetrols

The 2-methyltetrols have been widely chosen as chemical tracers for isoprene-derived secondary organic aerosols. While they are often assumed to be relatively unreactive, a laboratory study reported that pure erythritol particles (an analog of 2-methyltetrols) can be heterogeneously oxidized by gas-phase OH radicals at a significant rate. This might question the efficacy of these compounds as tracers in aerosol source-apportionment studies. Additional uncertainty could arise as organic compounds and inorganic salts often coexist in atmospheric particles. To gain more insights into the chemical stability of 2-methyltetrols in atmospheric particles, this study investigates the heterogeneous OH oxidation of pure erythritol particles and particles containing erythritol and ammonium sulfate (AS) at different dry inorganic-to-organic mass ratios (IOR) in an aerosol flow tube reactor at a high relative humidity of 85 %. The same reaction products are formed upon heterogenous OH oxidation of erythritol and erythritol–AS particles, suggesting that the reaction pathways are not strongly affected by the presence and amount of AS. On the other hand, the effective OH uptake coefficient, $γ_{eff}$, is found to decrease by about a factor of ~20 from 0.45±0.025 to 0.02±0.001 when the relative abundance of AS increases and the IOR increases from 0.0 to 5.0. One likely explanation is the presence of dissolved ions slows down the reaction rates by decreasing the surface concentration of erythritol and reducing the frequency of collision between erythritol and gas-phase OH radicals at the particle surface. Hence, the heterogeneous OH reactivity of erythritol and likely 2-methyltetrols in atmospheric particles would be slower than previously thought when the salts are present. Given 2-methyltetrols often coexist with a significant amount of AS in many environments, where ambient IOR can vary from ~1.89 to ~250, our kinetic data would suggest that 2-methyltetrols in atmospheric particles are likely chemically stable against heterogeneous OH oxidation under humid conditions.

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

Reactive Transport Modeling of Hydrogen Production from Serpentinization of Olivine

Hydrogen production from serpentinization of ultramafic rocks represents a promising natural pathway for generating carbon-free energy, yet its kinetics and controlling factors remain incompletely understood. A key challenge in advancing serpentinization research lies in the heterogeneity of porosity and permeability in rocks, which leads to nonuniform fluid velocity fields, as well as uncertainties in estimating reactive surface area and identifying appropriate mineral reaction equilibria. Additional complexities arise from the role of dissolved SiO 2 , Fe 2+ /Fe 3+ partitioning, and the limited effect of pH variations within the strongly alkaline regime on hydrogen yields. These challenges hinder straightforward extrapolation from laboratory tests to practical applications of hydrogen production from natural rocks. Here, in this work, we address these questions using a simulation-based reactive transport modeling framework calibrated against controlled laboratory experiments reported elsewhere. The model couples geochemical kinetics, multiphase flow, and mineralogical feedbacks, enabling systematic evaluation of how surface area, dissolved silica concentration, Fe redox state, temperature, and pressure govern serpentinization and H2 generation. We find that surface area exerts the strongest control on reaction rates and hydrogen yields, while Fe 2+ /Fe 3+ ratios act as secondary modulators. Elevated dissolved silica concentrations suppress hydrogen production but accelerate serpentine precipitation, whereas increasing pH beyond 12 within the strongly alkaline regime produces only marginal gains. Finally, we demonstrate that integrating targeted experiments with calibrated simulations offers a powerful and efficient approach for predicting hydrogen yields and assessing parameter trade-offs in industrial-scale applications. This integration can substantially reduce the experimental burden while improving predictive capability, thereby enhancing both the mechanistic understanding and the practical feasibility of hydrogen production from serpentinization.

08 HYDROGEN↗
DOE OSTI2025

Utilization of Existing Pipelines in Hydrogen Transport: Literature Review Report

This report critically reviews the flow behavior of hydrogen-natural gas (H 2 -NG) mixtures in pipelines and examines the critical factors of hydrogen integration into existing natural gas infrastructure. It addresses the choking behavior characterized by velocity increase and pressure drop, as well as the effects of flow restrictions and pressure losses during hydrogen transport. Computational and analytical models are used to investigate these effects, and their effects on thermodynamic properties and system performance are evaluated. The study also reviews the energy efficiency and flow dynamics of hydrogen and methane-hydrogen mixtures and optimizes the hydrogen flow rate. In addition, the effects of these mixtures on the flow characteristics are discussed in detail, with special emphasis on the compressibility factor (z factor) and fluid properties based on equations of state for hydrogen-natural gas mixtures. The study also analyzes the mixture ratios and highlights the thermophysical properties, flow dynamics, and hydrogen-blended natural gas application potential. These investigations assess flow stability, material interactions, and operational feasibility of transporting hydrogen mixtures through natural gas pipelines, which contribute to developing sustainable and efficient energy systems.

08 HYDROGEN↗
DOE OSTI2025

Chemical Reactor Network Modeling of Ammonia Rich-Quench-Lean Combustion Using a Partially Stirred Reactor Approach

Ammonia is a promising alternative to hydrogen with high energy density and favorable storage and transport characteristics. However, low flammability and a propensity for high nitrogen oxide (NO x ) emissions make direct utilization challenging. Recently, two-stage rich-quench-lean (RQL) combustion strategies have shown promise in achieving low NO x emissions with ammonia. In this approach, the rich stage serves to oxidize a portion of the fuel while thermally decomposing as much of the remaining ammonia as possible, generating hydrogen. In the second (lean) stage, air is rapidly introduced, burning out the hydrogen and residual ammonia. Two-stage RQL combustion of ammonia has been investigated in the open literature both experimentally and numerically. In general, idealized chemical reactor network (CRN) models predict NO x concentrations below those of 2D/3D computational fluid dynamics models and experiments. The primary drivers of these discrepancies may be largely attributed to finite rate mixing nonadiabatic operation. The typical CRN model is comprised of a perfectly-stirred-reactor (PSR), followed by a plug-flow-reactor (PFR), meant to represent the flame, and postflame zones, respectively. In the two-stage RQL approach two PSR-PFR networks are arranged sequentially, corresponding to the rich and lean stages, with secondary air injection in between. In the authors' past work, this arrangement has demonstrated the significant sensitivity of exit NO x to the rich stage equivalence ratio, while the amount of secondary air injection was shown to be less critical. In this paper, the CRN model is extended to (1) include the impacts of heat loss and (2) utilize a partially-stirred-reactor (PaSR) approach to study the impacts of mixing on emissions performance. Varying amounts of heat loss are applied to the rich relaxation zone to understand emissions performance and changes to optimization of equivalence ratio and residence time. Premixed and nonpremixed configurations are considered in the rich stage PaSR, with varying degrees of mixing intensity to study the interaction between mixing, transport, and kinetic timescales. Critically, the impact of mixing between hot products and secondary air injection is studied to understand practical injector needs. Results show unburnt ammonia leaving the rich stage as a primary contributor to NO x emissions – driven both by increased heat loss and reduced mixing rates. Furthermore, heat losses have been shown to create conditions that are conducive to increased N 2 O formation in the lean stage. In conclusion, the results of this study will be considered in the context of developing optimized two-stage RQL combustors for ammonia.

Combustion↗
DOE OSTI2020

Deposition Behavior of Supersaturated Silicic Acid on the Surfaces of Chlorite and Sericite around Geological Repository - 20138

Around the repository for HLW underground, cementitious materials used for the repository may complicate the migration of radionuclides in the host rock by the alteration of groundwater to a high alkaline condition. Specially, the alkaline groundwater dissolves silicate minerals of host rocks but the plenty of silicic acid dissolved from host rocks in groundwater will become supersaturated and deposit on the surface of flow-paths after decreasing the pH value to about 8 with dispersion and dilution. The deposition of silicic acid is presumed to recover and improve the barrier performance of host rocks for the radionuclide migration because of the narrowing of flow paths and the decrease of the permeability. In this study, the deposition behavior of supersaturated silicic acid was examined by the deposition experiments and the deposition rate constants, k [m/s], were evaluated from the experimental data. As the minerals of the surface of flow-paths, this study focused on chlorite and sericite, hydrothermally-altered minerals distributed to fill fractures of granite. In the experiments, deposition behavior of silicic acid on the powder of chlorite or sericite was observed after the pH of the solution was adjusted to 8 in order to make the silicic acid supersaturated. The solution of CaCl{sub 2} and NaCl were added into a silicic acid to simulate Ca{sup 2+} from the cementitious materials of the repository, Na{sup +} of the saline groundwater in Japan. From the results, the value of k was estimated to be about the order of 10{sup -10}∼10{sup -9} m/s for both chlorite and sericite and the value was almost the same as those for amorphous silica in the previous study. Therefore, the deposition of silicic acid may be assessed by similar apparent deposition rate constants even in actual underground containing various silicate minerals. Moreover, the Damkoehler number which is the dimensionless number expressed by the ratio of characteristic time of mass transfer and chemical reaction was calculated to be 10{sup 4} or more from the value of k. This value shows that the deposition reaction of silicic acid was estimated as much more predominant than the groundwater flow under the situation simulated in this study. These results suggest that the deposition of supersaturated silicic acid on the flow paths composed with various silicate minerals will be able to retard the migration of radionuclides because of narrowing micro flow-paths. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Technology Demonstration of a High-Pressure Swirl Oxy-Coal Combustor

This technical report presents the exploration of the design and prototyping of a High-Pressure Swirl Oxy-coal Combustor. Pressurized oxy-coal combustion systems have the potential to improve efficiency along with an increased carbon capture rate. Reduction of flue gas at higher pressure, smaller system size, and capital cost reductions render high-pressure oxy-coal systems particularly attractive as next-generation energy-producing systems. High-pressure oxy-coal combustion systems are a recent concept, and thus operability issues of combustor designs for such systems are not fully understood. Significant challenges exist to maintain oxy-coal combustion stability at elevated pressure and a high CO 2 diluent environment. Although a body of knowledge exists for high-pressure oxygen combustion in rocket engines (or similar applications), it is yet to be strategized how these fundamental concepts can be translated to low-temperature CO 2 diluent combustion regimes. The realization of the pressurized oxy-coal based systems requires combustor components to be designed and demonstrated for an operating pressure over 10 bar. However, pressurized oxy-coal combustor design information at this pressure range and scale relevant to validate those proposed systems is currently limited. Experimental data from MWth scale oxy-coal combustors are needed to identify the optimal trade-off between net efficiency and systems size. The proposed effort is aimed at demonstrating a 1 MWth down-fired swirl Oxy-Coal combustor and investigate the interrelation between combustor operating conditions (pressure; flame stability; flue gas recirculation ratio) and conversion efficiencies to minimize oxygen requirements. One of the key challenges is to configure burner design (i.e., swirl number and injector) and operating conditions for high-pressure oxy-coal combustion systems. These experiments differ from current systems partly due to the high theoretical flame temperature and related burner operability issues associated with oxy-combustion. An ASPEN PLUS® model study for 550 MWe TIPS and ENEL pressurized oxy-coal systems with CO 2 recirculation was performed to evaluate system design, subsystems sizing, and operating condition determination. The system analysis effort included TRL and technology gap determination of subsystems and critical components. This information was scaled to develop design requirements (design pressure and flue gas recirculation: RR Flue Gas = $\frac{m_{flue}}{m_{total}}$) for the 1 MWth combustor. The effects of a wide range of carbon dioxide recirculation ratios on the thermal efficiency of ENEL and TIPS cycles are studied. The pressure of 10 bar and 80 bar are used for ENEL and TIPS cycles, respectively. The thermal efficiency of ENEL is significantly higher than the efficiency of TIPS at a pressure of less than 10 bar. The insights from system analysis were then used to design a 1 MWth swirl oxy-coal combustor. Flame temperature analysis and material strength analysis was performed to determine the combustor thickness. The structural integrity of the combustor was validated by finite element analysis using Abacus® and Hypermesh®. Feasibility of igniters and secondary burners are investigated in successful high-pressure oxy-methane combustion. The secondary burners are designed in such a way that it can operate between 100 to 500 kW firing input. Three generations of the pintle injector were designed based on swirl numbers (S=0, 0.9, and 1.2). Key pintle injector parameters such as pintle size, pintle orifice size, spray pattern were investigated by cold flow tests. Information from these tests was used to modify injector design for smooth and successful operation. A 5 mm pintle orifice size was decided upon as the optimum size for oxy-coal operation for the combustor. Shadow sizing experiments were performed to identify the atomization rate of each injector. Different coal water slurry mixtures (30 – 50% coal by wt% in the mixture) at various total momentum ratios (TMR) were investigated for this purpose. These experiments provided decisive information to choose the best design of the injector. The injector with 1.2 swirl provided higher atomization in all cases than other designs. The mean equivalent droplet size of the jet was similar at different TMR and mixture ratios using this injector, thus making it suitable for use in most cases. Therefore, the 1.2 swirl-pintle injector was chosen for the shakedown test. The combustor and other sub-systems, including feed systems and control and data acquisition, have been manufactured, assembled, and integrated. The total system integration and installation began on July 1, 2020. The shake-down tests and initial operational capability demonstration are expected to be completed by September 30, 2020.

01 COAL, LIGNITE, AND PEAT↗