Search NASA⌕ Search

SEARCH · Search NASA

Results for “Mixing front”

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.

305 records · Page 17
All sourcesNASA NTRSDOE OSTI
NASA NTRS2015

Effects of Unsteady Flow Interactions on the Performance of a Highly-Loaded Transonic Compressor Stage

The primary focus of this paper is to investigate the loss sources in an advanced GE transonic compressor design with high reaction and high stage loading. This advanced compressor has been investigated both experimentally and analytically in the past. The measured compressor efficiency is significantly lower than the efficiency calculated with various existing tools based on RANS and URANS. The general understanding is that some important flow physics in this modern compressor design are not represented in the current tools. To pinpoint the source of the efficiency miss, an advanced test with detailed flow traverse was performed for the front one and a half stage at the NASA Glenn Research Center. In the present paper, a Large Eddy Simulation (LES) is employed to determine whether a higher-fidelity simulation can pick up any additional flow physics that can explain past efficiency miss with RANS and URANS. The results from the Large Eddy Simulation were compared with the NASA test results and the GE interpretation of the test data. LES calculates lower total pressure and higher total temperature on the pressure side of the stator, resulting in large loss generation on the pressure side of the stator. On the other hand, existing tools based on the RANS and URANS do not calculate this high total temperature and low total pressure on the pressure side of the stator. The calculated loss through the stator from LES seems to match the measured data and the GE data interpretation. Detailed examination of the unsteady flow field from LES indicates that the accumulation of high loss near the pressure side of the stator is due to the interaction of the rotor wake with the stator blade. The strong rotor wake interacts quite differently with the pressure side of the stator than with the suction side of the stator blade. The concave curvature on the pressure side of the stator blade increases the mixing of the rotor wake with the pressure side boundary layer significantly. On the other hand, the convex curvature on the suction side of the stator blade decreases the mixing and the suction side blade boundary layer remains thin. The jet velocity in the rotor wake in the stator frame seems to magnify the curvature effect in addition to inviscid redistribution of wake fluid toward the pressure side of the blade.

LES↗
NASA NTRS2006

Single- and Multiple-Objective Optimization with Differential Evolution and Neural Networks

Genetic and evolutionary algorithms have been applied to solve numerous problems in engineering design where they have been used primarily as optimization procedures. These methods have an advantage over conventional gradient-based search procedures became they are capable of finding global optima of multi-modal functions and searching design spaces with disjoint feasible regions. They are also robust in the presence of noisy data. Another desirable feature of these methods is that they can efficiently use distributed and parallel computing resources since multiple function evaluations (flow simulations in aerodynamics design) can be performed simultaneously and independently on ultiple processors. For these reasons genetic and evolutionary algorithms are being used more frequently in design optimization. Examples include airfoil and wing design and compressor and turbine airfoil design. They are also finding increasing use in multiple-objective and multidisciplinary optimization. This lecture will focus on an evolutionary method that is a relatively new member to the general class of evolutionary methods called differential evolution (DE). This method is easy to use and program and it requires relatively few user-specified constants. These constants are easily determined for a wide class of problems. Fine-tuning the constants will off course yield the solution to the optimization problem at hand more rapidly. DE can be efficiently implemented on parallel computers and can be used for continuous, discrete and mixed discrete/continuous optimization problems. It does not require the objective function to be continuous and is noise tolerant. DE and applications to single and multiple-objective optimization will be included in the presentation and lecture notes. A method for aerodynamic design optimization that is based on neural networks will also be included as a part of this lecture. The method offers advantages over traditional optimization methods. It is more flexible than other methods in dealing with design in the context of both steady and unsteady flows, partial and complete data sets, combined experimental and numerical data, inclusion of various constraints and rules of thumb, and other issues that characterize the aerodynamic design process. Neural networks provide a natural framework within which a succession of numerical solutions of increasing fidelity, incorporating more realistic flow physics, can be represented and utilized for optimization. Neural networks also offer an excellent framework for multiple-objective and multi-disciplinary design optimization. Simulation tools from various disciplines can be integrated within this framework and rapid trade-off studies involving one or many disciplines can be performed. The prospect of combining neural network based optimization methods and evolutionary algorithms to obtain a hybrid method with the best properties of both methods will be included in this presentation. Achieving solution diversity and accurate convergence to the exact Pareto front in multiple objective optimization usually requires a significant computational effort with evolutionary algorithms. In this lecture we will also explore the possibility of using neural networks to obtain estimates of the Pareto optimal front using non-dominated solutions generated by DE as training data. Neural network estimators have the potential advantage of reducing the number of function evaluations required to obtain solution accuracy and diversity, thus reducing cost to design.

Rai, Man Mohan↗
NASA NTRS1988

Cut-and-connect of two antiparallel vortex tubes

Motivated by an early conjecture that vortex cut-and-connect plays a key role in mixing and production of turbulence, helicity and aerodynamic noise, the cross-linking of two antiparallel viscous vortex tubes via direct numerical simulation is studied. The Navier-Stokes equations are solved by a dealiased pseudo-spectral method with 64 cubed grid points in a periodic domain for initial Reynolds numbers Re up to 1000. The vortex tubes are given an initial sinusoidal perturbation to induce a collision and keep the two tubes pressed against each other as annihilation continues. Cross-sectional and wire plots of various properties depict three stages of evolution: (1) Inviscid induction causing vortex cores to first approach and form a contact zone with a dipole cross-section, and then to flatten and stretch; (2) Vorticity annihilation in the contact zone accompanied by bridging between the two vortices at both ends of the contact zone due to a collection of cross-linked vortex lines, now orthogonal to the initial vortex tubes. The direction of dipole advection in the contact zone reverses; and (3) Threading of the remnants of the original vortices in between the bridges as they pull apart. The crucial stage 2 is shown to be a simple consequence of vorticity annihilation in the contact zone, link-up of the un-annihilated parts of vortex lines, and stretching and advection by the vortex tube swirl of the cross-linked lines, which accumulate at stagnation points in front of the annihilating vortex dipole. It is claimed that bridging is the essence of any vorticity cross-linking and that annihilation is sustained by stretching of the dipole by the bridges. Vortex reconnection details are found to be insensitive to asymmetry. Modeling of the reconnection process is briefly examined. The 3D spatial details of scalar transport (at unity Schmidt number), enstrophy production, dissipation and helicity are also examined.

Melander, Mogens V.↗
DOE OSTI2020

Mixed-Material Scintillators (MMSS Quarterly Report FY20Q2)

This project aims to invent, model, and prototype mixed-material scintillator systems (MMSSs), a new class of radiation detectors that use heterogenous internal structures of different scintillating materials to detect additional properties of radiation. These internal structures can be produced using additive manufacture (3D printing) of scintillator, a currently emerging application of additive manufacture technology. MMSSs combine the low cost and complexity of conventional scintillation detectors with capabilities currently only available in more expensive and complex detectors. By identifying promising MMSS designs, this project will enable a new class of detectors to meet DNN’s mission needs for SNM detection. This quarter, we delivered one the MMSS Particle-ID Modeling report, a major deliverable of this project. This report describes our results on simulations of the particle-ID (PID) class of detectors. We found that these detectors outperform competing approaches for neutron-gamma discrimination, neutron source pointing, and neutron spectroscopy. Further highlights of the report are described below. During this quarter, We alerted HQ to a change in schedule for the remaining simulations planned for this project: we plan to put off further simulation and reporting of the position-resolving (PR) detector class until July. This change in schedule has allowed us to go into additional detail in the PID work that is both very promising and motivates our upcoming proposal. The other active task in this project is prototyping 3D-printed scintillators with the aim of realizing the designs we’ve simulated. On that front, we had three focuses: measuring the light output, increasing the size of printed parts, and combining blue and green scintillators. We succeeded in measuring the light output, showing an output of 30% of commercial standard. This result leaves room for improvement but is within striking distance of where we need to be. We were on track to demonstrate a key requirement for increasing the size of parts and also to show combined green/blue prints before the laboratory moved to a Minimum Safe (MinSafe) Operations posture as a result of local Shelter In Place orders due to Covid-19.

42 ENGINEERING↗
DOE OSTI2022

ExaSGD: 2022 Kernel Thrust Activities

The Kernel Thrust milestone ADSE22-407 covers the development of device-capable optimization algorithms and solvers technologies required by the ExaSGD project’s software stack in order to solve security-constrained alternating current optimal power flow (SC-ACOPF) problems on emerging exascale architectures. To this extent, in FY22 the main objective of the Kernel Thrust was (i) provide sparse optimization solver that runs efficiently on hardware accelerator devices (i.e., NVIDIA and AMD GPUs) to perform intra-node computations, (ii) strengthen the reliability and increase the performance of the mixed-dense sparse (MDS) solver of HiOp for deployment on the FY22 target architectures, Summit and Crusher, and (iii) increase performance by improving the mathematical algorithm and refining the parallel MPI-based implementation of the coarse-grain parallel solver HiOp-PriDec for capabilities deployment on the FY22 target architectures, Summit and Crusher. This document presents the developments and contributions done by the Kernels Thrust Team in FY22 toward completion of the above-mentioned objectives. These contributions progressed along four main development (sub)thrusts: (1) Design and implementation of a sparse optimization solver for use on hardware accelerators; (2) Improvement of the mathematical algorithm and of the parallel implementation of HiOp-PriDec to ensure readiness and efficient coarse-grain parallelism for FY23 target exascale machine; and (3) Support Software and Application Development Thrusts of the exaSGD project in their deployment of the project’s software stack on AMD- and NVIDIA-based architectures. The development of the sparse optimization solver (thrust 1 above) was new in FY22 and resulted in a new sparse solver in HiOp (available as of version 0.6). The second development thrust was a continuation of the efforts from FY21 and improved the mathematical algorithm and the communication strategy of the HiOp-PriDec solver. The last developement thrust is a large collaborative effort. Namely, the project’s teams from multiple labs (LLNL, PNNL, ORNL, and NREL) performed large-scale demonstration of the ExaSGD software stack, namely the optimization solvers of HiOp interfaced with the modeling front-end ExaGO and the stochastic sampler PowerScenarios. These demonstration efforts solved large-scale instances of the SC-ACOPF challenge problem of medium network sizes (10, 000-bus system) and large number of contingencies on Summit (NVIDIA accelerators) and Crusher (AMD accelerators) systems at ORNL.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Direct observations and characterization of crack closure during microstructurally small fatigue crack growth via in-situ high-energy X-ray characterization

This study provides direct observation of the crack closure mechanism of a naturally occurring, tortuous, 3D microstructurally small fatigue crack (SFC) in additively manufactured Inconel 718. In-situ non-destructive characterization is performed using high-energy X-ray diffraction techniques to capture the evolution of the 3D microstructure and micromechanical response in the vicinity of the crack front. Based on the stress state of twelve grains analyzed at the crack tip, the crack closure events of the SFC front was found to be spatially heterogeneous with respect to loading progression governed by the local stress state of the grains (specifically the stress reversal from compression to tension). From this analysis, three grains that displayed different degrees of crack closure were further investigated, based on the orientation of the crack relative to the grains and the associated modality of crack growth. The stress normal to the crack plane and the associated degree of Mode I crack behavior were correlated with events of the crack opening earlier during the loading cycle. This was further corroborated by diffraction spot spreading analysis that quantified the crystallographic lattice distortion caused by the opening crack. Additionally, the detailed characterization of the opening behavior of the grains located at the crack tip and their associated states of stress elucidates the mechanism governing crack closure and will inform future modeling efforts of this phenomenon.

36 MATERIALS SCIENCE↗
DOE OSTI2023

A 32-channels readout ASIC for X-ray spectrometry and tracking in the GAPS experiment

Here this work describes the architecture and the experimental results from the characterization of a 32-channels mixed-signal Application-Specific Integrated Circuit (ASIC) developed for the readout of the lithium-drifted silicon, Si(Li), detectors of the General AntiParticle Spectrometer (GAPS) experiment dedicated to searching for dark matter. The instrument is designed for the identification of antiprotons, antideuterons and antihelium nuclei from cosmic rays during an Antarctic balloon mission scheduled for late 2024. A full custom integrated circuit, named SLIDER32 (32-channels Si-LI DEtector Readout) ASIC, has been produced in a commercial 180 nm CMOS technology. The ASIC is comprised of 32 low-noise analog readout channels featuring dynamic signal compression to comply with the wide input range, an 11-bit SAR ADC and a digital back-end section which is responsible for channel setting and for sending digital information to the data acquisition system (DAQ). The circuit design criteria and the experimental results are discussed in the paper.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

GRAPH — an readout ASIC for large MCP based detectors

We present a programmable 16 channel, mixed signal, low power readout ASIC, having the project historically named Gigasample Recorder of Analog waveforms from a PHotodetector (GRAPH). It is designed to read large aperture single photon imaging detectors using micro channel plates for charge multiplication, and measuring the detector's response on crossed strips anodes to extrapolate the incoming photon position. Each channel consists of a fast, low power and low noise charge sensitive amplifier, which provides a myriad of coarse and fine programmable options for gain and shaping settings. Further, the amplified signal is recorded using, to our knowledge novel, the Hybrid Universal sampLing Architecture (HULA) ADC. A kind of mixed signal double buffer memory, that enables concurrent waveform recording, and selected event digitized data extraction. The sampling frequency is freely adjustable between few kHz up to 125 MHz, while the chip's internal digital memory holds a history 2048 samples for each channel, with a digital headroom of 12 bits. An optimized region of interest sample-read algorithm allows to extract the information just around the event pulse peak, while selecting the next event, thus substantially reducing the operational dead time. The chip is designed in 130 nm TSMC CMOS technology, and its power consumption is around 47 mW per channel.

47 OTHER INSTRUMENTATION↗
NASA NTRS1992

The interactions of atmospheric cosmogenic radionuclides with spacecraft surfaces

The discovery of the cosmogenic radionuclide Be-7 on the front surface of the Long Duration Exposure Facility (LDEF) spacecraft has opened opportunities to investigate new phenomena in several disciplines of space science. The experiments performed for this work show that the Be-7 results only if the source of the isotope is the atmosphere through which the spacecraft passed. We should expect that the uptake of beryllium in such circumstances will depend on the chemical form of the Be and the chemical nature of the substrate. It was found that the observed concentration of Be-7 does, in fact, differ between metal surfaces and organic surfaces such as PTFE (teflon). It is noted, however, that: (1) organic surfaces, even PTFE, are etched by the atomic oxygen found under these orbital conditions, and (2) the relative velocity of the species is 8 km(exp -1)s relative to the surface and the interaction chemistry and physics may differ from the norm. The Be-7 is formed by spallation of O and N nuclei under cosmic ray proton bombardment. The principal source region is at altitudes of 12-15 km. While very small quantities are produced above 300 km, the amount measured on the LDEF was 3 to 4 orders of magnitude higher than expected from production at orbital altitude. The most reasonable explanation is that Be-7 is rapidly transported from low altitudes by some unknown mechanism. The process must take place on a time scale similar to the half-life of the isotope (53 days). Many other isotopes are produced by cosmic ray reactions, and some of these are suited to measurement by the extremely sensitive methods of accelerator mass spectrometry. A program was initiated to search for these isotopes and it is hoped that such studies will provide new methods for studying mixing in the upper atmosphere.

Gregory, J. C.↗
NASA NTRS1991

The interactions of atmospheric cosmogenic radionuclides with spacecraft surfaces

The discovery of the cosmogenic radionuclide Be-7 on the front surface (and the front surface only) of the Long Duration Exposure Facility (LDEF) spacecraft has opened opportunities to investigate new phenomena in several disciplines of space science. The experiments performed for this work show that the Be-7 results only if the source of the isotope is the atmosphere through which the spacecraft passed. We should expect that the uptake of beryllium in such circumstances will depend on the chemical form of the Be and the chemical nature of the substrate. It was found that the observed concentration of Be-7 does, in fact, differ between metal surfaces and organic surfaces such as PTFE (teflon). It is noted, however, that: (1) organic surfaces, even PTFE, are etched by the atomic oxygen found under these orbital conditions, and (2) the relative velocity of the species is 8 km(exp -1)s relative to the surface and the interaction chemistry and physics may differ from the norm. The Be-7 is formed by spallation of O and N nuclei under cosmic ray proton bombardment. The principal source region is at altitudes of 12-15 km. While very small quantities are produced above 300 km, the amount measured on the LDEF was 3 to 4 orders of magnitude higher than expected from production at orbital attitude. The most reasonable explanation is that Be-7 is rapidly transported from low altitudes by some unknown mechanism. The process must take place on a time scale similar to the half-life of the isotope (53 days). Many other isotopes are produced by cosmic ray reactions, and some of these are suited to measurement by the extremely sensitive methods of accelerator mass spectrometry. A program was initiated to search for these isotopes and it is hoped that such studies will provide new methods for studying mixing in the upper atmosphere.

Gregory, John C.↗
DOE OSTI2023

Front-End Engineering and Design: Project Tundra Carbon Capture System (Final Report)

As the next phase of Project Tundra, Minnkota Power Cooperative (MPC) completed a front-end engineering and design (FEED) study to install a post-combustion CO 2 capture system (CCS) at the MPC operated (Square Butte Electric Cooperative-owned) Milton R. Young Station (MRYS) Unit 2 (MRY2), a 477-MW power plant fueled by North Dakota lignite. The project team completed the FEED with Fluor’s Econamine FG Plus℠ (EFG+) technology and is taking the next steps leading up to start of construction. Team members included FEED technical lead Fluor Enterprises (Fluor); owner’s engineer and balance of plant (BOP) engineer Burns & McDonnell (BMcD); leading carbon capture consultants (David Greeson Consulting, Hunt International Energy Services); environmental consultants (AECOM, Agora Environmental), other BOP engineering consultants (Golder Associates, Nels); cost-share funding agency, the North Dakota Industrial Commission (NDIC); and the Energy & Environmental Research Center (EERC). Project Tundra’s goal is to implement carbon capture, utilization, and storage (CCUS) in North Dakota, while preserving the use of lignite. Future options could lead to revitalizing legacy oil fields and creating a new CO 2 enhanced oil recovery (EOR) industry. The topic of this report is the FEED study for the carbon dioxide capture portion of Project Tundra. Building on the findings of a pre-FEED study for MRY2, the key deliverables contained in this FEED study are: a) design, costing, and performance data needed to commence project financing activities; b) engineering and material balances required to file for all project permits; and c) a final project schedule. Based on the results of the previous pre-FEED study, MPC and its team evaluated two options to provide the large amount of steam needed to operate the CCS: 1) installation of new natural gas package boilers and a new pipeline to supply natural gas to the facility and 2) extraction of steam from the existing MRYS steam turbine generator. The FEED study was performed using the package boilers option, as it was deemed to have the lowest technical risk, overall cost, and cost uncertainty. A key feature of the design, however, was the mixing the natural gas boilers’ flue gas with the MRY2 flue gas. This increased the size of the CCS system. This design also enabled the tie-in of Unit 1 flue gas for times when Unit 2 is offline. The CCS was designed to capture a combined 12,978 short tons per day (STPD) of CO 2 from the flue gases produced by the existing MRY Unit 2 along with flue gas produced by the new boiler package within the CCS. The 12,978 STPD CO 2 capture is achieved by normally processing 100% of the total available MRY Unit 2 flue gas (considering full load operation) and 100% of the boiler package flue gas, then removing 90% of the CO 2 available from the processed flue gas streams.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2023

COMPUTATIONAL MODELING OF IGNITION AND PREMIXED FLAME PROPAGATION INITIATED BY A PRE-CHAMBER TURBULENT JET

Addressing the pressing need for reduced carbon emissions, Turbulent Jet Ignition (TJI) emerges as a promising technology for ultra-lean combustion, offering enhanced thermal efficiencies and minimized cyclic variability in spark-ignited engines. To facilitate rapid testing and integration of this technology, a robust computational modeling framework is crucial. This study delves into the predictive capabilities of computational models for main-chamber ignition and premixed flame propagation using a single-cycle TJI rig measured by Biswas et al. (Applied Thermal Engineering, vol 106, 2016). Employing an open-source compressible flow simulation solver with Large Eddy Simulation (LES) for turbulence modeling, the investigation integrates the conventional Laminar Finite Rate Chemistry (LFRC) model alongside the transported Probability Density Method (PDF) for turbulence-chemistry interaction. A fully-consistent Eulerian Monte-Carlo Fields (EMCF) method is utilized to approximate the transported PDF, while Interaction by Exchange with Mean is employed to close micro-mixing terms in stochastic differential equations. A reduced chemical reaction mechanism with 21 species and 84 reactions (DRM-19) is used for solving chemical kinetics, and a double Gaussian energy deposition model is used to approximate the spark ignition in the pre-chamber. An unstructured O-grid mesh with 0.3 million cells in the prechamber and 1 million cells in the main chamber is employed. Results are divided into two phases: pre-chamber initialization and full TJI simulations. Validation of the predicted pre-chamber flame propagation and the lean ignition in the main-chamber is carried out by using available experimental data. Under quiescent conditions, both the LFRC and transported PDF methods largely underestimate the flame speed and subsequent pressure growth in the pre-chamber. A linear momentum forcing technique is applied to investigate the impact of initial turbulence in the pre-chamber, demonstrating a notable influence on flame propagation. Fine-tuning of the forcing coefficient reproduces the sudden pressure growth observed in the experiment. The experimentally validated pre-chamber simulation serves as the initial condition for the full TJI simulations. It is found that the LFRC model fails to predict lean-ignition in the main-chamber, resulting in a misfiring event. Incorporation of turbulence-chemistry interaction using the transported PDF method substantially improves the prediction of the ignition event in the main-chamber, achieving fair qualitative agreement and quantitative validation of combustion parameters within ±10% of the reported experimental data. The rich simulation results consisting of a full set of statistical description of the thermo-chemical states enable us to gain deep insights into the ignition mechanisms in the main chamber, which is limited when done experimentally. A novel dual ignition phenomenon is revealed in the TJI rig for the first time. Initially, a primary ignition kernel is formed at a downstream location which eventually detaches from the main jet. As the jet momentum decreases, a secondary ignition event follows, this time at a more upstream location which eventually combines with the primary ignition kernel to form a single connected flame front. Investigation of these ignition sequences in chemical composition space reveal distinct differences between the two. The primary ignition event in the main-chamber is followed by a large concentration of active radicals from the pre-chamber jet, accelerating the chain-branching steps, characterizing what has been referred to as flame ignition. In contrast, the secondary ignition occurs in the absence of active radicals in the pre-chamber jet, hence characterized as jet ignition. Further analysis of the effect of pre-chamber jet characteristics on lean ignition in the main-chamber is conducted by setting up cases with different initial pressure ratios (por) between the two chambers, a non-dimensional parameter, ranging from 1.2 to 3.2. As the initial pressure ratio increases, jet momentum increases, with dual ignition observed in cases above por= 2.2. Case with por= 3.2 lead to misfiring. The effect of ignition sequence on global combustion characteristics of TJI is analyzed. Dual ignition events lead to non-monotonicity in combustion characteristics such as global reaction progress variable, flame penetration, and global heat release rate. In dual ignition events, although the rate of fuel consumption and global heat release rate is initially lower, the secondary ignition leads to a sudden increase in flame surface area, resulting in a sudden jump and promoting the overall performance of the TJI system.

42 ENGINEERING↗
NASA NTRS2004

Exhaust Nozzle for a Multitube Detonative Combustion Engine

An improved type of exhaust nozzle has been invented to help optimize the performances of multitube detonative combustion engines. The invention is applicable to both air-breathing and rocket engines used to propel some aircraft and spacecraft, respectively. In a detonative combustion engine, thrust is generated through the expulsion of combustion products from a detonation process in which combustion takes place in a reaction zone coupled to a shock wave. The combustion releases energy to sustain the shock wave, while the shock wave enhances the combustion in the reaction zone. The coupled shockwave/reaction zone, commonly referred to as a detonation, propagates through the reactants at very high speed . typically of the order of several thousands of feet per second (of the order of 1 km/s). The very high speed of the detonation forces combustion to occur very rapidly, thereby contributing to high thermodynamic efficiency. A detonative combustion engine of the type to which the present invention applies includes multiple parallel cylindrical combustion tubes, each closed at the front end and open at the rear end. Each tube is filled with a fuel/oxidizer mixture, and then a detonation wave is initiated at the closed end. The wave propagates rapidly through the fuel/oxidizer mixture, producing very high pressure due to the rapid combustion. The high pressure acting on the closed end of the tube contributes to forward thrust. When the detonation wave reaches the open end of the tube, it produces a blast wave, behind which the high-pressure combustion products are expelled from the tube. The process of filling each combustion tube with a detonable fuel/oxidizer mixture and then producing a detonation repeated rapidly to obtain repeated pulses of thrust. Moreover, the multiple combustion tubes are filled and fired in a repeating sequence. Hence, the pressure at the outlet of each combustion tube varies cyclically. A nozzle of the present invention channels the expansion of the pulsed combustion gases from the multiple combustion tubes into a common exhaust stream, in such a manner as to enhance performance in two ways: (1) It reduces the cyclic variations of pressure at the outlets of the combustion tubes so as to keep the pressure approximately constant near the optimum level needed for filling the tubes, regardless of atmospheric pressure at the altitude of operation; and (2) It maximizes the transfer of momentum from the exhaust gas to the engine, thereby maximizing thrust. The figure depicts a typical engine equipped with a nozzle according to the invention. The nozzle includes an interface section comprising multiple intake ports that couple the outlets of the combustion tubes to a common plenum. Proceeding from its upstream to its downstream end, the interface section tapers to a larger cross-sectional area for flow. This taper fosters expansion of the exhaust gases flowing from the outlets of the combustion tubes and contributes to the desired equalization of exhaust combustion pressure. The cross-sectional area for flow in the common plenum is greater than, or at least equal to, the combined cross-sectional flow areas of the combustor tubes. In the common plenum, the exhaust streams from the individual combustion tubes mix to form a single compound subsonic exhaust stream. Downstream of the common plenum is the throat that tapers to a smaller flow cross section. In this throat, the exhaust gases become compressed to form a compound sonic gas stream. Downstream of the throat is an expansion section, which typically has a bell or a conical shape. (The expansion section can be truncated or even eliminated in the case of an air-breathing engine.) After entering the expansion section, the exhaust gases expand rapidly from compound sonic to compound supersonic speeds and are then vented to the environment. The basic invention admits of numerous variations. For example, the combustion tubes can be arranged around the central axin a symmetrical or asymmetrical pattern other than the one shown in the figure. For another example, the flow cross-sectional area(s) of one or more of the intake ports in the interface section, of the common plenum, the throat, and/or the expansion section can be varied, either symmetrically or asymmetrically, to adjust dynamics of the exhaust stream or to direct the thrust vector away from the central axis.

Bratkovich, Thomas E.↗
DOE OSTI2024

A Mid-Century Net-Zero Scenario for the State of Wyoming and its Economic Impacts

Clean hydrogen has the potential to help achieve 10% economy-wide emissions reductions by 2050 relative to 2005, promote energy security and resilience, and develop a new economy in the United States. In 2030, the hydrogen economy could create about 100,000 new jobs to build new capital projects and clean hydrogen infrastructure. The Wyoming Energy Authority recently announced the state’s energy strategy, which establishes a goal of net-zero emissions by 2050. Under all likely scenarios, achieving a mid-century net-zero target will pose challenges and create opportunities for Wyoming’s energy sector. If executed properly, the transition could favorably affect the state’s economy overall in the long term. This research program examines the economic impact of fossil energy production in Wyoming and provides various predictions for future energy mixes to achieve net-zero emissions. Preliminary work suggests that Wyoming-based hydrogen production could have significant economic benefits and job creation implications for Wyoming. This study further assesses Wyoming’s opportunities to create hydrogen-based industries, assess economic impacts, identify knowledge gaps and research needs, and create a Hydrogen Center of Excellence to accelerate commercialization and deployment. This project helped to understand Wyoming's areas of focus for research and development and identified its areas of strength and potential challenges in creating a hydrogen ecosystem. As a result of this study, we estimate that for blue hydrogen produced from coal and gas resources, the overall cost reduction will be driven mainly by the carbon-sequestration tax credit and the improvement in carbon capture. Mature technologies, like SMR and PSA, will make limited contributions. They have no or limited reductions from an additional capacity deployment in future costs. We also understand the importance of continued support from public and private sectors for Carbon Capture and Storage (CCS)-related research, development, and demonstration programs at federal and state levels. The successful and efficient production of blue hydrogen requires a unique blend of energy resources, geology, regulation, law, and infrastructure. Wyoming has the distinction of meeting all these demands. The team also estimates that the availability and command of water resources accessible for hydrogen production are crucial for developing new projects. Water treatment, use, and disposal after treatment will also make projects possible. Primarily, this is relevant for hydrogen made using renewable energy. Wyoming has one of the best wind resource capacity in the nation. Harnessing this resource is challenging due to limited transmission line availability. Hydrogen could become one of the solutions to the stranded resource problem, primarily if the water availability challenge is addressed. Using produced oil & gas water could help to solve the problem. A commonly cited barrier to the expansion of hydrogen markets is the cost associated with constructing new pipelines, which typically require large amounts of capital to develop. Wyoming already possesses much of the export infrastructure needed to connect Wyoming’s hydrogen production with major markets across the West Coast, Pacific Northwest, Midwest, and Front Range regions of the United States, where a large portion of Wyoming’s natural gas is already transported. In addition to transportation by pipeline, rail transportation of hydrogen has also proven feasible. Wyoming uses its extensive railway system to transport large amounts of coal to its export partners across the United States. By using cryogenic or compressed-gas cars, Wyoming has the potential to add hydrogen to its existing network of railroad energy exports. The same technology may also be applied to hydrogen transport via trucks traveling interstate highways. Wyoming’s workforce is ready to meet the demands of clean hydrogen development. Many of the skills and training needed for hydrogen production are the same skills already possessed by Wyoming’s oil & gas and coal workforce. Many government and industry leaders expect clean hydrogen and other low-carbon energy projects to generate significant job growth and to recruit many already-trained oil & gas and coal workers whose jobs may be displaced. As energy companies seek to penetrate the markets for Wyoming hydrogen production, there is a natural mutual benefit to Wyoming’s workers and companies seeking to launch projects with the assistance of a trained workforce. Wyoming’s university and community college system have adopted several programs to ensure that highly qualified engineers and other technically skilled employees continue to graduate with skills to support the development of hydrogen and other innovative energy projects moving forward. Throughout the project, stakeholder outreach and education took many forms, including meetings with several major companies in the industry, collaborating with local government organizations, educational organizations, and national laboratories, tribal outreach and engagement, the sponsoring of several hydrogen-focused projects in many departments throughout the University of Wyoming, and developing a collaboration with international universities. The products of these collaborations consist of working relationships with several companies in the industry, educational institutions, national labs, and local government, as well as strong connections with individuals who will play an essential role in the success of the Hydrogen Energy Research Center.

08 HYDROGEN↗
NASA NTRS2017

Computational Modeling to Limit the Impact Displays and Indicator Lights Have on Habitable Volume Operational Lighting Constraints

NASA has demonstrated an interest in improving astronaut health and performance through the installment of a new lighting countermeasure on the International Space Station. The Solid State Lighting Assembly (SSLA) system is designed to positively influence astronaut health by providing a daily change to light spectrum to improve circadian entrainment. Unfortunately, existing NASA standards and requirements define ambient light level requirements for crew sleep and other tasks, yet the number of light-emitting diode (LED) indicators and displays within a habitable volume is currently uncontrolled. Because each of these light sources has its own unique spectral properties, the additive lighting environment ends up becoming something different from what was planned or researched. Restricting the use of displays and indicators is not a solution because these systems provide beneficial feedback to the crew. The research team for this grant used computer-based computational modeling and real-world lighting mockups to document the impact that light sources other than the ambient lighting system contribute to the ambient spectral lighting environment. In particular, the team was focused on understanding the impacts of long-term tasks located in front of avionics or computer displays. The team also wanted to understand options for mitigating the changes to the ambient light spectrum in the interest of maintaining the performance of a lighting countermeasure. The project utilized a variety of physical and computer-based simulations to determine direct relationships between system implementation and light spectrum. Using real-world data, computer models were built in the commercially available optics analysis software Zemax Optics Studio(c). The team also built a mockup test facility that had the same volume and configuration as one of the Zemax models. The team collected over 1200 spectral irradiance measurements, each representing a different configuration of the mockup. Analysis of the data showed a measurable impact on ambient light spectrum. This data showed that obvious design techniques exist that can be used to bind the ambient light spectrum closer to the planned spectral operating environment for the observer's eye point. The following observations should be considered when designing an operational environment that is dominated by computer displays. When more light is directed into the field of view of the observer, the greater the impact it will make on various human factors issues that depend on spectral shape and intensity. Because viewing angle has a large part to play in the amount of light flux on the crewmember's retina, beam shape, combined with light source location is an important factor for determining percent probable incident flux on the observer from any combination of light sources. Computer graphics design and display lumen output are major factors influencing the amount of spectrally intense light projected into the environment and in the viewer's direction. Use of adjustable white point display software was useful only if the predominant background color was white and if it matched the ambient light system's color. Display graphics that used a predominantly black background had the least influence on unplanned spectral energy projected into the environment. Percent reflectance makes a difference in total energy reflected back into an environment, and within certain architectural geometries, reflectance can be used to control the amount of a light spectrum that is allowed to perpetuate in the environment. Data showed that room volume and distance from significant light sources influence the total spectrum in a room. Smaller environments had a homogenizing effect on total light spectrum, whereas light from multiple sources in larger environments was less mixed. The findings indicated above should be considered when making recommendations for practice or standards for architectural systems. The ambient lighting system, surface reflectance, and display and indicator implementation all factor into the users' spectral environment. A variety of low-cost solutions exist to mitigate the impact of light from non-architectural lighting systems, and much potential for system automation and integration of display systems with the ambient environment. This team believes that proper planning can be used to avoid integration problems and also believes that human-in-the-loop evaluations, real-world test and measurement, and computer modeling can be used to determine how changes to a process, display graphics, and architecture will help maintain the planned spectral operating lighting environment.

Clark, T. A.↗
DOE OSTI2022

Plasma Wall Interaction with 3-D Plasma Boundaries

The interaction of the edge plasma and the material surfaces is one of the most critical challenges on the path to harness fusion power as new, fundamental energy source. This challenge typically combines the thrust to reach high density, low temperature (detached) plasmas in front of the divertor target plates as well as understanding the plasma material interaction (PMI) in particular in this regime. The combination of both research thrusts represents an extraordinarily challenging subject encompassing spatial scales spanning nanometers to meters in all states of matter and across a broad energy range. Modeling capabilities, which help to interpret data from nowadays experiments and enable extrapolation to future devices are urgently required. This is in particular true for toroidal magnetic confinement devices with three-dimensional (3D) plasma boundaries. Such plasma boundary geometries occur in tokamaks, when small amplitude magnetic perturbations are used to stabilize the unruly edge plasma or in stellarators, that are inherently 3D plasma confinement devices. In this project, the impact of 3D plasma boundaries on the plasma material interaction (PMI) was assessed. This work focused on plasma boundary conditions, in which high-density conditions at the material surfaces yield mitigation of the otherwise immense heat and particle loads that these materials would see. These so-called high recycling and eventually detached plasma regimes are of great interest for future reactor operation. In the project, key features that are unique to 3D boundaries were explored in comparison to canonically assumed axisymmetric plasma edge situations in tokamaks. In particular, the relevance of the 3D boundary situation in the extrapolation to the plasma boundary solution at ITER, the next step fusion energy experiment under construction as a multi-national, world-wide large-science experiment in southern France, has been explored. The EMC3-EIRENE plasma edge fluid and kinetic neutral transport code has been advanced to cope with the challenging and unprecedented conditions in the ITER boundary plasma including 3D fields that are planned to be used to suppress harmful edge instabilities, the so-called edge localized modes. This is a vital integration challenge for ITER and the results from this grant have provide a leading capability for this assessment. It was shown that the detachment process in a 3D edge solution for ITER follows the recycling regimes that are known from axisymmetric solutions, but that multiple plasma exhaust channels connected to the material surfaces are established which feature individual recycling characteristics. Because these channels touch the material surfaces in the divertor in a 3D geometry, the compatibility with the plasma material interaction (PMI), including erosion and impurity generation has been found to be an important part of the integration challenge. To address this, the fully 3D plasma material interaction code ERO2 has been adapted to these ITER specific geometries and a homogeneous mixing model was implemented, that allows to consider the mixing of Be and Was used at ITER in the PMI modeling. This model enhancement has been used to study non-local migration of Be in the JET ITER like wall configuration and it has been shown that with this model such complex migration processes in ITER relevant plasma shapes and with ITER relevant plasma boundary conditions can be addressed. The combined modeling approach using EMC3-EIRENE as a plasma boundary transport code and the ERO2 specialized PMI model will be an asset for the continued preparations of ITER operation as well as for Fusion Pilot Plant efforts that have emerged in the U.S. during the evolution of this grant. The predictive capability of this numerical tool has been validated at the DIII-D US national fusion facility. Here, dedicated plasma edge diagnostics were implemented to measure the impurity household around a 3D edge plasma during ELM suppression by 3D fields. Dedicated experiments with local material probes using these diagnostics and the state-of-the-art suite of boundary measurements at DIII-D have shown that the 3D perturbation of the plasma edge that is excreted by such 3D control fields yield a perturbation of the plasma boundary flux structure and hence also of the resulting PMI. The 3D boundary plasma is composed out of helical magnetic flux channels that intersect the divertor targets at an angle relative to the main guiding field, i.e., the toroidal magnetic field component of the tokamak. A similar effect has been measured as well on limiter surfaces during the startup campaign at the new stellarator experiment Wendelstein 7-X. These experiments ad initial analysis with the ERO plasma material interaction model, suggested that the place of erosion for a given particle from the surface and its re-deposition can be different in such 3D field geometries yielding potentially a significant level of net-erosion. This is not the case for axisymmetric solutions, where it was shown in the past that the eroded particles are effectively re-deposited into gaps produced by erosion at the same position and hence the net-erosion levels are small. For ITER, the quest to suppress the ELMs and at the same time maintain the integrity of the divertor is an issue, which these fundamental findings will help to resolve. The coupling of this work to the extrapolation in the ITER program has been addressed by both the PI and the lead numerical scientist being ITER Science Fellows in the duration of the contract and forward. A second focus in the exploration of 3D boundary effects on tokamaks and stellarators has been set on the measurement of helium exhaust features with such 3D fields. This is important because He represents the ash of the fusion process and needs to be exhausted. It was shown that 3D field application compatible with suppression of ELMs yields an increase of the helium exhaust performance. The ratio of the effective helium confinement time over the energy confinement time was reduced by almost 50% which demonstrated that the impact of helium accumulation in the plasma core with respect to the confinement of energy to sustain the fusion reaction is significantly improved with such 3D control fields. It was shown that this is the case for tokamaks as well as stellarators. At the Large helical Device in Japan, a similar enhancement of the helium exhaust features when small amplitude additional 3D fields were applied was measured. This is an important additional function of 3D field application and its impact on ITER is presently being studied in combination with investigations of helium exhaust in 3D field geometries of stellarator devices.

3D plasma edge transport↗
DOE OSTI2023

CCN data from TAMU TRACER campaign in the Houston TX region from July to September 2022

During TRACER, the Texas A&M Rapid Onsite Atmospheric Measurements Van (ROAM-V) was deployed to capture airmasses behind (maritime) and ahead (continental) of the passage of the sea-breeze front through Houston. On select sampling days, ROAM-V sampled in the morning/mid-day on the coast and then transited to a second inland site for the afternoon/evening. The suite of instruments deployed on ROAM-V included a Condensation Particle Counter (CPC; GRIMM Model 5.403 CPC), Scanning Mobility Particle Sizer (SMPS; TSI 3750 detector, TSI 3082 classifier, TSI 3088 neutralizer, TSI 3081A Differential Mobility Analyzer), Cloud Condensation Nuclei counter (Droplet Measurement Technologies CCN Counter), micro pulse lidar (Droplet Measurement Technologies Micro Pulse LiDAR (miniMPL)), and a Davis Rotating Uniform size-cut Monitor (DRUM; DRUMAir 4-DRUM). Before sampling at each location, the latitude and longitude were recorded using the GPS on the phone application “My Altitude”. Onboard the ROAM-V, aerosol samples are drawn through a shared isokinetic inlet at a flow rate ranging from 3.5 to 7.0 LPM. A portion of this flow is directed through a cyclone impactor (Brechtel, Inc. Model SCC 0.732) and 0.5 LPM is directed to the CCN. To calculate particle losses, we used a two-step method. First, the measured SMPS size distributions were used to calculate particle loss through the inlet during sampling. Second, the corrected SMPS data was used to calculate the average of the total losses per scan down the CCN line. Then, the correction was applied to the CCN data. This calculation was done separately for each deployment location due to changes in the measured size distributions between locations. Particle loss from diffusion (based on Kesten, 1991 and Gormley, 1949), inertial impaction in 90-degree bends (based on Aerosol Measurement, 2011 and Crane, 1977), and cyclone impactor efficiency (based on Dirgo, 1985) were included in the loss calculation. When the SMPS was not sampling at a location (in the case of an instrument malfunction or operator error), the reported CPC data was corrected with an average of the total losses for the entire campaign at the specified deployment location (e.g., if we needed to correct Galveston data, then the average of all calculated losses at Galveston was taken). These flatline corrections were used for all data on 22/07/13, 22/07/20, 22/07/22, and the data from Galveston on 22/08/09. The supersaturation uncertainty is estimated conservatively at +/- 0.03%, where variation in the inlet temperature, pressure, and calibration technique prevents a more accurate measurement. Confidence in the reported supersaturation measurements is based on a pre-campaign calibration (following the methods from our previous work and Deng, 2014 based on Rose, 2008) in addition to inter-comparisons with the DOE for two days (22/08/18 and 22/09/01) where TAMU was co-located with AMF1. The inter-comparisons show good agreement between our instrument and the DOEs instrument on both days at all supersaturations. After the last inter-comparison on 22/09/01, there was no indication of a malfunction by our instrument through the rest of the campaign. Unfortunately, the instrument was dropped during demobilization. A post-campaign calibration was conducted, which showed a substantial departure from the pre-campaign calibration. The drop may have damaged the instrument’s ability to produce the desired supersaturations. Therefore, we do consider the data after 22/09/01 to be correct, but it should be used with caution. The CCN counter sampled for 3 minutes at each supersaturation setpoint (0.2, 0.4, 0.6, 0.8, 1.0, and 1.2%). At the end of a cycle, the instrument was set to 0.01% supersaturation for 5 minutes. The data is comprised of the last 60 seconds of each supersaturation set point (0.2, 0.4, 0.6, 0.8, 1.0, and 1.2%) to ensure the instrument stabilized and was able to reach thermal equilibrium. We removed the data during the periods where there were operational difficulties, setup, or maintenance. This data was collected for ARM Field Campaign AFC07055 and supported by DOE ASR grant DE-SC0021047. For any further questions, please feel free to contact the instrument PI, Sarah D. Brooks, sbrooks@tamu.edu. Rose et. al. Calibration and Measurement Uncertainties of a Continuous-Flow Cloud Condensation Nuclei Counter (DMT-CCNC): CCN Activation of Ammonium Sulfate and Sodium Chloride Aerosol Particles in Theory and Experiment. Atmos. Chem. Phys., 8, 1153-1179, 2008. Deng et. al. Using Raman Microspectroscopy to Determine Chemical Composition and Mixing State of Airborne Marine Aerosols over the Pacific Ocean. Aerosol Science and Technology, Vol 48, Issue 2, 2014. Kesten et. al. Calibration of a TSI Model 3025 Ultrafine Condensation Particle Counter. Aerosol Science and Technology, 15:2, 107-111, 1991. Gormley et. al. Diffusion from a Stream Flowing through a Cylindrical Tube. Proceedings of the Royal Irish Academy, Vol 52, 163-169, 1948. Aerosol Measurement: Principles, Techniques, and Applications, Third Edition. John Wiley & Sons, Inc, 2011. Crane et. al. Inertial Deposition of Particles in a Bent Pipe. Journal of Aerosol Science, Vol 8, 161-170, 1977. Dirgo et. al. Cyclone Collection Efficiency: Comparison of Experimental Results with Theoretical Predictions. Aerosol Science and Technology, 4:4, 401-415, 1985.

54 ENVIRONMENTAL SCIENCES↗