Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Slow control”

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

Kinetic Model for Moisture-Controlled CO 2 Sorption

The understanding of the sorption/desorption kinetics is essential for practical applications of moisture-controlled CO 2 sorption. We introduce an analytic model of the kinetics of moisture-controlled CO 2 sorption and its interpretation in two limiting cases. In one case, chemical reaction kinetics on pore surfaces dominates, in the other case, diffusive transport through the sorbent defines the kinetics. Here, we show that reaction kinetics, which is dominant in the first case, can be expressed as a linear combination of 1st and 2nd order kinetics in agreement with the static isotherm equation derived and validated in a previous paper. The interior transport kinetics can be described by non-linear diffusion equations. By combining all carbon species into a single equation, we can eliminate — in certain limits — the source terms associated with chemical reactions. In this case, the governing equation is ∂θ/∂t = –∇ · (–D eff ∇ θ ). For a sorbent in a form of a flat sheet or a membrane, one can maintain the same functional form of a diffusion equation by introducing a generalized effective diffusivity D M that combines contributions from both surface chemical reaction kinetics and interior diffusive transport kinetics. Experimental data of transient CO 2 flux in a preconditioned commercial anion exchange membrane fit well to the 1st order model as long as very dry states are avoided, validating the theory. The observed DM for a preconditioned commercial anion exchange membrane ranges from 6.6× 10 -14 to 7.1× 10 -14 m 2 s -1 at 35°C. These small values compared to typical ionic diffusivities imply a very slow kinetics, which will be the largest issue that needs to be addressed for practical application. The collected transient CO 2 flux data are used to predict the magnitude of a continuous CO 2 pumping flux in an active membrane that transports CO 2 against a CO 2 concentration gradient. The pumped CO 2 flux is supported by water flux due to a water concentration gradient.

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

Heterogeneous oxidation of amorphous organic aerosol surrogates by O 3 , NO 3 , and OH at typical tropospheric temperatures

Typical tropospheric temperatures render possible phase states of amorphous organic aerosol (OA) particles of solid, semisolid, and liquid. This will affect the multiphase oxidation kinetics involving the organic condensed-phase and gaseous oxidants and radicals. To quantify this effect, we determined the reactive uptake coefficients (γ) of O 3 , NO 3 , and OH by substrate films composed of single and binary OA surrogate species under dry conditions for temperatures from 213 to 313 K. A temperature-controlled coated-wall flow reactor coupled to a chemical ionization mass spectrometer was applied to determine γ with consideration of gas diffusion transport limitation and gas flow entrance effects, which can impact heterogeneous reaction kinetics. The phase state of the organic substrates was probed via the poke-flow technique, allowing the estimation of the substrates' glass transition temperatures. γ values for O 3 and OH uptake to a canola oil substrate, NO 3 uptake to a levoglucosan and a levoglucosan / xylitol substrate, and OH uptake to a glucose and glucose / 1,2,6-hexanetriol substrate have been determined as a function of temperature. We observed the greatest changes in γ with temperature for substrates that experienced the largest changes in viscosity as a result of a solid-to-liquid phase transition. Organic substrates that maintain a semisolid or solid phase state and as such a relatively higher viscosity do not display large variations in heterogeneous reactivity. From 213 to 293 K, γ values of O 3 with canola oil, of NO 3 with a levoglucosan / xylitol mixture, and of OH with a glucose / 1,2,6-hexanetriol mixture and canola oil, increase by about a factor of 34, 3, 2, and 5, respectively, due to a solid-to-liquid phase transition of the substrate. These results demonstrate that the surface and bulk lifetime of the OA surrogate species can significantly increase due to the slowed heterogeneous kinetics when OA species are solid or highly viscous in the middle and upper troposphere. This experimental study will further our understanding of the chemical evolution of OA particles with subsequent important consequences for source apportionment, air quality, and climate.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

TULIP: An RNA-seq-based Primary Tumor Type Prediction Tool Using Convolutional Neural Networks

Background: With cancer as one of the leading causes of death worldwide, accurate primary tumor type prediction is critical in identifying genetic factors that can inhibit or slow tumor progression. There have been efforts to categorize primary tumor types with gene expression data using machine learning, and more recently with deep learning, in the last several years. Methods In this paper, we developed four 1-dimensional (1D) Convolutional Neural Network (CNN) models to classify RNA-seq count data as one of 17 highly represented primary tumor types or 32 primary tumor types regardless of imbalanced representation. Additionally, we adapted the models to take as input either all Ensembl genes (60,483) or protein coding genes only (19,758). Unlike previous work, we avoided selection bias by not filtering genes based on expression values. RNA-seq count data expressed as FPKM-UQ of 9,025 and 10,940 samples from The Cancer Genome Atlas (TCGA) were downloaded from the Genomic Data Commons (GDC) corresponding to 17 and 32 primary tumor types respectively for training and validating the models. Results: All 4 1D-CNN models had an overall accuracy of 94.7% to 97.6% on the test dataset. Further evaluation indicates that the models with protein coding genes only as features performed with better accuracy compared to the models with all Ensembl genes for both 17 and 32 primary tumor types. For all models, the accuracy by primary tumor type was above 80% for most primary tumor types. Conclusions: We packaged all 4 models as a Python-based deep learning classification tool called TULIP (TUmor CLassIfication Predictor) for performing quality control on primary tumor samples and characterizing cancer samples of unknown tumor type. Further optimization of the models is needed to improve the accuracy of certain primary tumor types.

Jones, Sara↗
DOE OSTI2024

Co-, Ni- and Fe-rich grain-boundary phases enhance creep resistance in θ'-strengthened Al-Cu alloys

Here, microstructural evolution and creep response were investigated in the cast Al-5.0Cu-0.3Mn-0.2Zr (wt.%) alloy with and without addition of slow-diffusing, intermetallic-forming elements Fe, Ni, or Co. Baseline Al-5.0Cu-0.3Mn-0.2Zr alloy exhibits high creep resistance at 300 °C, up to ~75 MPa, which is attributed to high-aspect-ratio, intragranular θ'-Al 2 Cu precipitates that effectively suppress dislocation climb. However, θ'-Al 2 Cu precipitate-free zones form along grain boundaries upon heat treatment, whose extent is amplified during subsequent creep. Such weak regions experience dislocation creep, leading to strain localization and acceleration of grain-boundary sliding. Adding Ni and Co, individually or in combination, leads to the formation of grain-boundary precipitates (Al 9 Co 2 , Al 3 Ni 2 ) which are resistant to coarsening, thus suppressing the formation of θ'-Al 2 Cu precipitate-free zones. This microstructure provides high creep resistance at stresses up to 75–80 MPa, with strain rates much lower than in unmodified Al-5.0Cu-0.3Mn-0.2Zr. Adding Fe, which results in extensive decoration of grain boundaries with coarsening-resistant Al 7 Cu 2 Fe, and then increasing the Cu content to compensate for the Cu loss to this new phase, leads to a new Al-7.4Cu-1.6Fe-0.3Mn-0.2Zr alloy with creep resistance at 300 °C that surpasses known cast aluminum alloys. Adding Fe to improve the creep resistance of Al-Cu alloys is both cost-effective and sustainable. Our findings offer guidelines applicable to various alloy systems on controlling the evolution of precipitate-free zones and its ensuing effects on creep deformation.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Understanding Biases in Sample Preparation Techniques for Coupled Scanning Electron Microscopy and MAMA PuO 2 Morphological Analysis

In this project, the scanning electron microscopy (SEM) sampling method used during the statistical design study (SDS) was investigated to determine if any sampling biases were present in the analyzed data. Using standard particle size distribution powders from the National Institute of Standards and Technology (NIST 1984 standard reference material) with the origin wet dispersion method, it was determined that a bias to smaller particles was present. This was supported by theoretical calculations using Stokes’ law to determine the settling rate of spherical particles of roughly the same size and mass as those found in the SDS. Based on the theoretical calculations, it was determined that the settling rate for each of the 76 powder sets in the SDS could be unique based on specific particle shape and mass distributions, making a universal correction factor/formula not applicable. Therefore, priority shifted to developing an improved wet dispersion method that significantly reduced the particle settling rate for all particle size and shapes. This was achieved by replacing the original solvent (isopropyl alcohol) with a heavy liquid (lithium heteropolytungstates), which dramatically slowed the settling rate and allowed for the capture of a suitable homogeneous aliquot. SEM imaging and Morphological Analysis for Material Attribution (MAMA) software analysis were conducted on the NIST standard, and the SEM/MAMA data were compared to data captured by a dynamic image analysis particle size analyzer. The resulting data confirmed that the new wet dispersion method does indeed deliver an improved representative aliquot to the SEM stub. For instance, in the NIST certificate, the average particle size is ~17.1 µm ± 2.2 µm with a normal distribution. The initial wet dispersion method resulted in a drastically reduced average particle size of 6.1 µm in addition to a non-representative heavy bi-modal distribution whereas the improved LST wet dispersion method resulting in an average particle size that was much closer to the NIST certificate (12.7 µm) with a similar normal distribution. Although the improved method was still short of the NIST certificate average, atomic force microscopy analysis determined that the resulting ~20-25% reduction in size was due to particles sinking into the carbon sticky tape used for SEM imaging. It is believed that that this bias can be calibrated in a much more predicable manner than the original settling rate bias. In addition, the matching normal distribution curves between the NIST certificate and the heavy liquid method indicate a much-improved representative aliquot has been sampled and imaged. A surrogate CeO 2 powder was used to reflect PuO 2 more accurately and to aid in implementing radiological controls and shielding. The resulting data sets from the SEM/MAMA method and the particle size analyzer give almost identical average particle sizes and particle distribution statistics. Future work will re-analyze several select runs from the SDS to determine if morphological signatures can be found with the improved sampling method.

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

Initial Assessment of TeF 6 Adsorption

Research and development that supports the management of off-gases from nuclear fuel reprocessing has historically been focused on the off-gas streams that arise from aqueous reprocessing technology. However, as Gen-IV reactor development pathways move toward deployment, alternative spent nuclear fuel (SNF) processing and disposition pathways have been considered more actively. This work is focused upon aspects of fluoride volatility (FV) processing. The versatility of this method for deployment against multiple types of spent fuel encourages the continued advancement of both the primary separations processes and the secondary processes, including waste treatment, material control and accountability, and engineering designs. Recent work noted that TeF 6 , the most highly volatile fluorinated compound produced during FV processing, did not have a clear abatement technology recommended in the literature. Thus, this work performed scoping tests on activated alumina and copper shot to assess whether they could be used to remove TeF 6 from gas streams that bear F 2 . Previous work on this topic was not well described in the literature and was not performed with excess F 2 in the stream as would be typical of spent fuel processing via FV. To support an understanding of the concentration of TeF 6 contacting the adsorbent beds, a series of preliminary testing identified the TeF 6 production rate and the equilibrium concentration of TeF 6 in the gas stream contacting the adsorbent. Excess F 2 was determined to not affect the ability of activated alumina to remove TeF 6 from the gas stream quickly and completely. A determination of whether excess F 2 affected the distribution depth of Te in the sorbent bed is still pending analysis of the used sorbent. Literature suggests that when activated alumina is near saturation, TeF 6 could migrate from the sorbent bed. Future testing should investigate this possibility. Copper metal did not adsorb TeF 6 in the presence of F 2 across the sorbent temperature range of 50 to 335°C. F 2 was fully removed by the copper bed. Although preliminary thermodynamics would indicate adsorption to be energetically favorable, other factors that impact adsorption (e.g., slow kinetics, excess fluorine on the copper surface, an unfavorable transition state) are likely preventing the adsorption of TeF 6 at an easily measurable rate. The work also investigated TeF 6 production from multiple forms of Te in the temperature range of 100 to 250°C. No previous study had assessed the initial reaction rates for this process. A carefully designed study allowed determination of the activation energy for the production of TeF 6 from Te metal. The substantial amount of data collected during this study merits analysis beyond what is described here. A more in-depth kinetic analysis will be pursued. Additional analytical results will provide the ability to benchmark adsorption coefficients for TeF 6 , understand the distribution of TeF 6 within the alumina bed, and better understand the effect of F 2 partial pressure on Te fluorination. The data from this report, as supplemented by these additional analyses, will be submitted to a peer-reviewed journal.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Probing Carbon Mineralization Mechanisms in Pore and Bulk Fluids by Harnessing Architected Calcium Silicates

The ability to synthesize materials with well-controlled pore structures gives us unprecedented control over probing fluid interactions with reactive interfaces and advancing calibrated insights into coupled chemo-morphological interactions. One of the primary challenges in developing crystalline silicate materials lies in achieving ordered pore structures. Existing approaches of producing amorphous mesoporous metal silicates via sol–gel methods and heat-treatment of these materials to produce crystalline phases cause the pore structures in the amorphous phases to collapse. To overcome this challenge, carbon coating of amorphous mesoporous calcium silicate particles is carried out to retain the pore structure, while the material is heated to produce crystalline calcium silicate with calcium sulfate inclusions. The pore diameter in these materials is about 3.9 nm, with a surface area and a pore volume of 28.75 m 2 /g and 0.092 cm 2 /g, respectively. The mechanisms of carbon mineralization are investigated by reacting architected calcium silicate with 1 M Na 2 CO 3 and monitoring the evolution in the structural phases using operando wide-angle X-ray scattering (WAXS) measurements. Formation of stable calcium carbonate polymorph or calcite and metastable calcium carbonate polymorph or vaterite in pore and bulk fluids, respectively, resulting from the reaction between Na 2 CO 3 and CaSiO 3 , are noted. The mechanisms associated with carbon mineralization are delineated using ReaxFF molecular dynamics (MD) simulations. The surface dissolution reaction is initiated by 2H + ions that replace a Ca 2+ ion in the Ca–silicate matrix. Ca 2+ ions in the solution initially react with water to form calcium hydroxide and eventually form calcium (bi)carbonate. A slow and gradual increase in the formation of sodium silicate in the solution resulting from the reactions of silicic acid or the silicon dioxide reaction with sodium hydroxide is noted. When carbon mineralization occurs in environments bearing interfacial fluids, calcite is the dominant calcium carbonate polymorph, as determined using experiments with pore fluids and molecular-scale simulations. In conclusion, these studies provide fundamental insights into the mechanisms underlying the carbon mineralization of calcium silicate informed by experiments and molecular-scale simulations.

Calcium↗
DOE OSTI2023

Plasma-assisted deflagration to detonation transition in a microchannel with fast-frame imaging and hybrid fs/ps coherent anti-Stokes Raman scattering measurements

Our study examines kinetic enhancement by nanosecond dielectric barrier discharge (ns-DBD) plasma on fuellean dimethyl ether (DME), oxygen (O 2 ), and argon (Ar) premixtures during deflagration to detonation transition (DDT) experiments in a microchannel. Non-equilibrium plasma produces active species and radicals and creates fast and slow heating of a mixture to promote ignition due to electronic and vibrational excitation. Experiments have been conducted to examine the influence of the plasma discharge on the premixture and on the resultant deflagration to detonation transition (DDT) onset time and distance through the use of high speed imaging and one-dimensional, two-beam, femtosecond/picosecond, coherent anti-Stokes Raman scattering (CARS). A highspeed camera is used to trace the time histories of flame front position and velocity and to identify the dynamics and onset of DDT. The results show that plasma discharge can nonlinearly affect the onset time and distance of DDT. It is shown that a small number of plasma discharge pulses prior to ignition result in reduced DDT onset time and distance by 60% and 40%, respectively, when compared to the results without pre-excitation by ns discharges. The results also show that an increase of plasma discharge pulses results in an extended DDT onset time and distance of 224% and 94%, respectively. Time history of the deflagration wave speed of DME and the analysis of ignition timescale under the choking condition of the deflagration front suggest low temperature ignition may play a role for DME near the isobaric choking condition of the burned gas and the DDT. Plasmainduced conversion of the reactive mixture was assessed via the O 2 to CO 2 ratio as measured through fs/ps CARS during the DBD discharges. The present experiments demonstrate the ability of non-equilibrium plasma to alter the chemistry of DME/O 2 /Ar premixtures in order to control DDT for applications in advanced propulsion engines.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Switching-cycle-based Startup for Grid-connected Inverters

Conventional inverter startups, or grid synchronization, are hindered by slow dynamics and inrush current issues, which impede the integration of more renewable energy resources into the power grid. This paper overcomes the barriers by introducing a novel switching cycle-based startup approach for grid-connected inverters, eliminating the need for voltage sensors and phase-locked loops (PLLs). The proposed method surpasses the bandwidth constraints of traditional PLL-based techniques and achieves grid synchronization in just two switching cycles. Relying solely on sensed three-phase currents, synchronization is accomplished without inrush currents or the need for additional hardware or increased software complexity, resulting in a faster inverter startup. The paper also presents a sensitivity analysis to discuss the impact of critical parameters on the method's accuracy. As a result, the effectiveness of the proposed technique is demonstrated through simulation and experimental results.

14 SOLAR ENERGY↗
DOE OSTI2024

Ambient energy for buildings: Beyond energy efficiency

The following Key Messages comprise the salient findings of this study: 1. Ambient energy (from sun, air, ground, and sky) can heat and cool buildings; provide hot water, ventilation and daylighting; dry clothes; and cook food. These services account for about three-quarters of building energy consumption and a third of total US demand. Biophilic design (direct and indirect connections with nature) is an intrinsic adjunct to ambient energy systems, and improves wellness and human performance. 2. The current strategy of electrification and energy efficiency for buildings will not meet our climate goals, because the transition to an all-renewable electric grid is too slow. Widespread adoption of ambient energy is needed. Solar-heated buildings also flatten the seasonal demand for electricity compared to all-electric buildings, reducing required production capacity and long-term energy storage. In addition, ambient-conditioned buildings improve resilience by remaining livable during power outages. 3. National policies, incentives, and marketing should be enacted to promote ambient energy use. Federal administrative priorities should reflect the importance of ambient energy for buildings. Use of ambient energy should be encouraged through existing and new building codes and standards. 4. Ambient energy system design tools are needed for architects, engineers, builders, building scientists, realtors, appraisers, and consumers. PVWatts is used over 100 million times per year for photovoltaic system design. A similar, simple, and accessible tool for ambient design is crucial. 5. Training on ambient energy is needed throughout secondary, post-secondary, and continuing education for workforce development. Currently, only about 10% of colleges teach courses on passive heating and cooling systems. 6. Ambient-conditioned buildings should be demonstrated in all US climate zones. Performance should be monitored and reported, with quantitative case studies made widely available. 7. While current technology is sufficient to build high-performance ambient buildings now, research is needed to develop new technologies to harness ambient energy more effectively and more economically. Such advancements will facilitate adoption of ambient energy technologies in a wider range of buildings, including retrofits. Examples include windows with much lower thermal losses, use of the building shell as thermal storage, alternative light-weight thermal storage systems, sky radiation cooling systems, automated controls for solar gains and passive cooling, and ground coupling.

14 SOLAR ENERGY↗
DOE OSTI2025

Staircases of passive and active scalar concentration in cellular flow

This paper develops a unified model for staircase formation in both passive and active scalar systems, building upon prior numerical studies by offering new heuristic and physical insights. While prior studies primarily reported numerical results, they did not explore the underlying unifying physics that governs both types of scalar transport; this work addresses that gap by identifying shared mechanisms across both cases. Results of studies of passive and active scalar staircase formation in cellular flows are presented. Staircase formation in cellular flows occurs due to the interplay of fast mixing within cells and slow transport across the inter-cell boundary. The cell boundary emerges as a de facto transport barrier. Special attention is focused on the effects of cellular fluctuations and noise upon staircase structure. A forced, fluctuating vortex array model is used to drive the underlying flow structure. Cellular Peclet number and staircase profile curvature are identified as figures-of-merit to quantify the resiliency of layering. These are related to simple, multi-scatterer scalar random walk models. Results for Peclet number and curvature scaling with flow excitation are presented. We also study staircases of magnetic potential evolving in two-dimensional magnetohydrodynamics as examples of layering of active scalar concentration. Formation of magnetic potential staircases is indeed observed. Flux expulsion inhibits the intercellular transport of magnetic potential and strengthens staircase barriers. Magnetic staircases can be supported against resistive decay by magnetic potential noise forcing. Implications for staircase formation in magnetic confinement experiments are discussed.

Control theory↗
DOE OSTI2020

GaAsP/Si Tandem Solar Cells: Pathway to Low-Cost, High-Efficiency Photovoltaics

Si is the dominant PV technology, now and for the foreseeable future, due to its extensive manufacturing infrastructure, supply chain, feedstock availability, and highly optimized degree of fabrication processes, which altogether has produced an economic scenario where PV electricity generation is often cheaper than conventional fossil based generation. In many places, the overarching goal of grid parity has been achieved, but further improvement in performance-cost metrics are still needed to sustain the continued LCOE reductions needed to not only compete with conventional generation, but displace it on a global scale; a matter of critical importance if we stand any hope of slowing climate change. Nevertheless, single-junction Si PV is already nearing its physical limit, both in performance and cost, and is thus cannot meet these long-term goals alone. To this end, we are working on the development of monolithic III-V/Si tandem solar cells, which improve upon the performance of pure Si by providing enhanced utilization (reduced thermalization) of high-energy photons. This architecture nominally combines the substantial existing knowledge base, manufacturing infrastructure, and low cost of Si PV with the high efficiencies afforded by the well-established multijunction approach — the only proven way to break the single-junction limit. Although the metal-halide perovskite/Si tandem architecture has garnered substantial attention in recent years, serious questions regarding reliability and service lifetime remain, whereas III-V PV has a proven track record, including in the harsh concentrator and space environments. Additionally, there are multiple fabrication approaches to producing III-V/Si tandem cells, but we are focused on monolithic epitaxial integration as it is the most likely to yield the lowest ultimate LCOE in a fully mature, scaled technology. In this work we have produced multiple generations of GaAsP/Si tandem solar cells, demonstrating a more than 10% absolute AM1.5G efficiency improvement within the time frame of the project, including two verified world records. We have done this using industry-standard fabrication methods, showing that this platform can ultimately be manufactured at scale using existing or only slightly upgraded Si and III-V tooling. Our scientific and engineering advances across a range of fundamental and applied areas – III-V/Si heteroepitaxial integration, defect control in metamorphic III-V epitaxy, fundamental materials-oriented solar cell design and modeling methodology, and more – have created clear pathways for continued advances toward the goal of >30% AM1.5G cell efficiency (and >25% module) and will serve to inform the broader research community for well beyond this immediate application. Techno-economic modeling indicates that our approach can indeed meet SunShot/SETO LCOE targets, but as with any “post-Si” technology there are difficult, but not insurmountable barriers, requiring continued focused research and development efforts.

14 SOLAR ENERGY↗
DOE OSTI2025

Plasma assisted NH 3 /H 2 /air ignition in nanosecond discharges with non-equilibrium energy transfer

Ammonia (NH 3 ), with its high energy density and easiness to store and transport as a hydrogen carrier, has become a promising alternative green fuel. However, its adoption in power generation is hindered by challenges such as low burning velocity, slow low-temperature oxidation, high NO x emissions, and ignition difficulty. Here, this work computationally investigates the effects of non-equilibrium energy transfer by nanosecond discharges on NH 3 ignition and flame propagation in an NH 3 /H 2 /air flow at 700 K and 1 atm. The simulation results demonstrate that NH 3 /air mixtures require a large ignition energy due to their large critical ignition radius. It is shown that adding 30 % hydrogen significantly reduces the critical ignition radius and minimum ignition energy. Two-dimensional modeling further shows a non-monotonic dependence of ignition kernel volume on the applied voltage and reduced electric field. The optimum ignition enhancement occurs at 200 Td where the generation of electronically excited species and radicals including N 2 (B), O( 1 D) and OH becomes most efficient. Higher voltages divert electron energy toward ionization, which makes it less effective for NH 3 ignition. The study also identifies an optimal electrode gap size for a given pulse energy. Smaller gap sizes increase deposited energy density, raising temperature and radical concentrations. However, excessive reduction of the gap distance reduces flame propagation speed due to the flame stretch effect in rich mixtures with the effective Lewis number greater than unity. A nonlinear relationship between pulse repetition frequency and ignition kernel volume is observed in a nanosecond pulsed high frequency discharge (NPHFD). An optimal frequency range of 200 kHz to 2 MHz is found when two pulses are used. In addition, an optimal number of pulses exists for each pulse repetition frequency, with higher frequencies requiring more pulses to maximize the overlap region. These findings provide critical insights on developing controlled plasma discharge techniques for efficient NH 3 ignition in reactive flows within internal combustion engines and gas turbines.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Dioxygen Activation Kinetics over Distinct Cu Site Types in Cu-Chabazite Zeolites

Cu-exchanged zeolites activate dioxygen to form active sites for partial methane oxidation (PMO), nitrogen oxide decomposition, and carbon monoxide oxidation. Apparent rates of O 2 activation depend both on the intrinsic kinetics of distinct Cu site types and the distributions of such sites within a given zeolite, which depend on the density and arrangement of the framework Al atoms. In this work, we use hydrothermal synthesis methods to control the arrangement of framework Al sites in chabazite (CHA) zeolites and, in turn, the distinct Cu site types formed. Time-resolved in situ Raman spectroscopy reveals the kinetics of O 2 adsorption and activation within these well-defined Cu-CHA materials and the concomitant structural evolution of copper–oxygen (Cu x O y ) complexes, which are interpreted alongside Cu(I) oxidation kinetics extracted from in situ X-ray absorption spectroscopy (XAS). Raman spectra of several plausible Cu x O y species simulated using density functional theory suggest that experimental spectra (λ ex = 532 nm) capture the formation of mono(μ-oxo)dicopper species (ZCuOCuZ). Transient experiments show that the timescales required to form Cu x O y structures that no longer change in Raman spectra correspond to the durations of oxidative treatments that maximize CH 3 OH yields in stoichiometric PMO cycles (approximately 2 h). Yet, these periods extend well beyond the timescales for the complete conversion of the initial Cu(I) intermediates to their Cu(II) states (<0.3 h, reflected in X-ray absorption near edge spectroscopy spectra), which demonstrates that Cu x O y complexes continue to evolve structurally following rapid Cu(I) oxidation. The dependence of ZCuOCuZ formation rates on O 2 pressure, H 2 O pressure, and temperature is consistent with a mechanism in which ZCuOH reduces to form ZCu + sites that bind molecular oxygen and form ZCu-O 2 intermediates. Subsequent reaction with proximate ZCu + forms bridging peroxo dicopper complexes that cleave O–O bonds to form ZCuOCuZ in steps facilitated by water. These data and interpretations provide evidence for the chemical processes that link rapid and kinetically irrelevant Cu oxidation steps (frequently probed by XAS and UV–vis spectroscopy) to the relatively slow genesis of reactive Cu complexes that form CH 3 OH during PMO. In doing so, we reveal previously unrecognized complexities in the processes by which Cu ions in zeolites activate O 2 to form active Cu x O y complexes, which underscore the insight afforded by judicious combinations of experimental and theoretical techniques.

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

RELAP-7 Application and Enhancement for FLEX Strategies and ATF Behavior under Extended Loss of AC Power Conditions

This report summarizes the results of a three-year research project sponsored by the U.S. Department of Energy (DOE) Nuclear Energy University Program (NEUP) to enhance and apply the RELAP-7 code by adding and improving several important components (e.g., a mechanistic Reactor Core Isolation Cooling (RCIC) system model) for thermal hydraulic studies of LWRs under ELAP conditions and evaluating the time available for transition to portable FLEX equipment. The project team included University of Massachusetts–Lowell, The Ohio State University, Texas A&M University, Idaho National Laboratory and Oak Ridge National Laboratory. In the Fukushima accident, it was found that the RCIC system played a crucial role in delaying core meltdown by almost three days in Fukushima Daiichi Unit 2, because of self-regulated operation of the steam driven RCIC turbine-pump injection system. Steam flow in the convergent-divergent nozzles of the RCIC Terry turbine is two-phase non-equilibrium transonic flow with homogenous nucleation condensation. To more accurately predict the dynamic process and behavior of the transonic compressible steam flow, a one-dimensional transient two-phase analytical model is presented. A simplified four-fluid model was employed in the present work with the consideration of four separate fluid fields: vapor, liquid film, entrained droplets and condensed droplets. The mass, momentum and energy interactions between the fluids were considered and modeled. An extended seven-equation non-equilibrium critical flow model was developed to obtain the critical pressure and velocities of each phase at the nozzle throat. To predict the wetness in the divergent section, a mechanistic nucleation condensation model was integrated in the nozzle analysis model, considering the generation and consequent growth of droplets. The governing differential equations on a staggered grid were discretized using the second-order Lax-Wendroff scheme with a flux limiter, and the Semi-Implicit Method for Pressure-Linked Equation (SIMPLE) algorithm was employed to solve the discrete linear system. To demonstrate the predictability and reliability of the physical models and the numerical method proposed in the present work, three representative nozzles were modeled and simulated. The results show good agreement with the available experimental data, even for condensation shock. Then, the 1D nozzle model was employed to obtain nozzle flow tables of the Terry turbine nozzle for different working pressures which can cover the operation pressure range of the RCIC system. A mechanistic RCIC turbine-pump system model was developed and implemented in the system code TRACE to simulate dynamic responses of the RCIC system under Beyond Design Basis Accident (BDBA) conditions. The turbine-pump governing equations are based on the control volume approach of the angular momentum balance. The physics based mechanistic RCIC model was developed using the TRACE control system components (i.e., signal variables, control blocks, and tables), and incorporated into a TRACE boiling water reactor (BWR) model. The TRACE model in this report has a detailed nodalization of the reactor pressure vessel (RPV), and all of the major flow paths and system components, including the safety relief valves (SRVs) and the containment suppression pool and drywell. Based on the nozzle flow tables generated from the 1D nozzle model developed, the turbine drive torque can be calculated from table lookup. Since the detailed specifications of the RCIC pump are unavailable, the homologous curves for a Bingham pump were used in the current pump component. A station black-out (SBO) accident test problem was selected to demonstrate the TRACE RCIC model. The short-term SBO simulations were performed for two cladding materials: Zircaloy and FeCrAl, to demonstrate the effect of the accident tolerant fuel cladding on fuel heat-up under BDBA accident conditions. The wetwell plays a vital safety role in SBO and other BWR accident scenarios in that it can reduce containment pressure and supply additional core make-up water. The suppression pool temperature distribution has a very large impact on both RPV and containment pressure. Thus, another novel contribution of the project comes mainly from an improved, systems-level wetwell model which can capture buoyancy-induced thermal stratification effects due to steam injection and condensation. A two-zone stratified wetwell model has been implemented in RELAP-7 and some results from that model are presented. This wetwell model is capable of simulating thermal stratification due to a low steam mass injection rate. With a low mass flow rate, the model assumes that all the steam condenses within the pipe and the resulting plume can be approximated with a purely buoyant, heat-source driven model. The wetwell model developed with these assumptions is adequate to simulate slow transients such as extended SBO transients.

42 ENGINEERING↗
DOE OSTI2025

Slow Strain Rate Testing of A537 Tank Wall Material

At Savannah River Site (SRS), High-Level Waste is stored in below-grade carbon steel tanks. This waste in part consists of sludge, salt cake, and/or supernate. Preparation of this waste for future processing involves dissolution of the salt cake layer. The salt dissolution process can create conditions that leave the carbon steel tanks susceptible to localized corrosion. The salt to be dissolved contains high concentrations of nitrate, that once released, create an environment that may be conducive to pitting corrosion and/or stress corrosion cracking (SCC) of carbon steel. The salt dissolution process also liberates interstitial liquid trapped between the salt crystals. This liquid is initially high in nitrite and hydroxide concentration. High pH and greater ratios of nitrite to nitrate act as inhibitors to minimize corrosion of carbon steel in high nitrate environments. However, as dissolution proceeds, the concentration of nitrate will increase, while the hydroxide and nitrite concentration of the interstitial liquid will deplete and become insufficient to prevent the onset of corrosion attack. Tank blending and the addition of inhibitors are used to ensure adequate concentrations of hydroxide and nitrite. However, this is not desirable during salt dissolution as it can reduce process efficiency and increase the amount of waste that needs processing. This testing program was designed to examine the risk of SCC associated with utilizing the pitting factor (PF) and nitrite/nitrate (NO 2 - /NO 3 - ) ratio limits for handling dissolved salt solutions at an elevated temperature in the carbon steel waste tanks. The previously identified limits are a PF of 1.2 and an NO 2 - /NO 3 - ratio of 0.15. The results indicate that as long as the NO 2 - /NO 3 - ratio exceeds 0.1 and the PF is above approximately 0.8, there is a discernible safety margin between the open circuit potential (OCP) and the critical cracking potential (CCP) observed during applied potential testing. However, this margin, defined by the difference between the OCP and CCP, is relatively narrow, ranging from 0.1 to 0.25 volts. This small margin raises concerns about potential shifts in OCP during waste retrieval operations, which could inadvertently increase the risk of SCC if the OCP approaches or exceeds the CCP. These results confirm that dissolved salt solutions provide a potent chemistry that, under certain conditions, makes carbon steel susceptible to SCC. The next question to consider is the influence these results have on decisions for storage and retrieval of waste from the tanks. For Type III/IIIA waste tanks, the risk of SCC remains very low. First, and most importantly, the post-weld stress relief of the tanks has reduced the residual stress near the welds. Thus, without the stress component, SCC risk is minimized. The material of construction (A537 Carbon steel) for the Type III/IIIA tanks is superior to the steel in its resistance to SCC than the steel that was utilized for the Type I, II, and IV tanks (A285 carbon steel). From a chemistry control standpoint for a Type III/IIIA tank directly involved with handling dissolved salt solutions, the PF and NO 2 - /NO 3 - ratio limits may be utilized wherein chemistry control provides an extra layer of defense against SCC. Chemistry control for a Type III/IIIA tank minimizes the risk for a tank that may receive the dissolved salt solution, particularly if that tank is a Type I, II, or IV waste tank. On the other hand, if the dissolved salt solution is handled by a Type I, II, or IV waste tank the risk of SCC is real. The potent chemistry, absence of stress relief, and inferior material result in a condition that is conducive to cracking. Efforts should be made to either avoid transferring waste that may not meet the PF and NO 2 - /NO 3 - ratio criteria to one of these tanks or if it is unavoidable, take measures to minimize the consequences of a leak. As shown by these tests, even if the PF and NO 2 - /NO 3 - ratio criteria are met, there is a risk that the tank potential may be disturbed in the positive direction and the risk of SCC increase.

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

Quantification of early nonpharmaceutical interventions aimed at slowing transmission of Coronavirus Disease 2019 in the Navajo Nation and surrounding states (Arizona, Colorado, New Mexico, and Utah)

During an early period of the Coronavirus Disease 2019 (COVID-19) pandemic, the Navajo Nation, much like New York City, experienced a relatively high rate of disease transmission. Yet, between January and October 2020, it experienced only a single period of growth in new COVID-19 cases, which ended when cases peaked in May 2020. The daily number of new cases slowly decayed in the summer of 2020 until late September 2020. In contrast, the surrounding states of Arizona, Colorado, New Mexico, and Utah all experienced at least two periods of growth in the same time frame, with second surges beginning in late May to early June. Here, we investigated these differences in disease transmission dynamics with the objective of quantifying the contributions of non-pharmaceutical interventions (NPIs) (e.g., behaviors that limit disease transmission). We considered a compartmental model accounting for distinct periods of NPIs to analyze the epidemic in each of the five regions. We used Bayesian inference to estimate region-specific model parameters from regional surveillance data (daily reports of new COVID-19 cases) and to quantify uncertainty in parameter estimates and model predictions. Our results suggest that NPIs in the Navajo Nation were sustained over the period of interest, whereas in the surrounding states, NPIs were relaxed, which allowed for subsequent surges in cases. Our region-specific model parameterizations allow us to quantify the impacts of NPIs on disease incidence in the regions of interest.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Verification of diesel spray ignition phenomenon in dual-fuel diesel-piloted premixed natural gas engine

Dual-fuel (DF) engines, in which premixed natural gas and air in an open-type combustion chamber is ignited by diesel-fuel pilot sprays, have been more popular for marine use than pre-chamber spark ignition (PCSI) engines because of their superior durability. However, control of ignition and combustion in DF engines is more difficult than in PCSI engines. In this context, here we focus on the ignition stability of n-heptane pilot-fuel jets injected into a compressed premixed charge of natural gas and air at low-load conditions. To aid understanding of the experimental data, chemical-kinetics simulations were carried out in a simplified engine-environment that provided insight into the chemical effects of methane (CH4) on pilot-fuel ignition. The simulations reveal that CH4 has an effect on both stages of n-heptane autoignition: the small, first-stage, cool-flame-type, low-temperature ignition (LTI) and the larger, second-stage, high-temperature ignition (HTI). As the ratio of pilot-fuel to CH4 entrained into the spray decreases, the initial oxidization of CH4 consumes the OH radicals produced by pilot-fuel decomposition during LTI, thereby inhibiting its progression to HTI. Using imaging diagnostics, the spatial and temporal progression of LTI and HTI in DF combustion are measured in a heavy-duty optical engine, and the imaging data are analyzed to understand the cause of severe fluctuations in ignition timing and combustion completeness at low-load conditions. Images of cool-flame and hydroxyl radical (OH*) chemiluminescence serve as indicators of LTI and HTI, respectively. The cycle-to-cycle and spatial variation in ignition extracted from the imaging data are used as key metrics of comparison. The imaging data indicate that the local concentration of the pilot-fuel and the richness of the surrounding natural-gas air mixture are important for LTI and HTI, but in different ways. In particular, higher injection pressures and shorter injection durations increase the mixing rate, leading to lower concentrations of pilot-fuel more quickly, which can inhibit HTI even as LTI remains relatively robust. Decreasing the injection pressure from 80 MPa to 40 MPa and increasing the injection duration from 500 µs to 760 µs maintained constant pilot-fuel mass, while promoting robust transition from LTI to HTI by effectively slowing the mixing rate. This allows enough residence time for the OH radicals, produced by the two-stage ignition chemistry of the pilot-fuel, to accelerate the transition from LTI to HTI before being consumed by CH4 oxidation. Thus from a practical perspective, for a premixed natural gas fuel–air equivalence-ratio, it is possible to improve the “stability” of the combustion process by solely manipulating the pilot-fuel injection parameters while maintaining constant mass of injected pilot-fuel. This allows for tailoring mixing trajectories to offset changes in fuel ignition chemistry, so as to promote a robust transition from LTI to HTI by changing the balance between the local concentration of the pilot-fuel and richness of the premixed natural gas and air. This could prove to be a valuable tool for combustion design to improve fuel efficiency or reduce noise or perhaps even reduce heat-transfer losses by locating early combustion away from in-cylinder walls.

42 ENGINEERING↗