Numerical simulation of a methane-oxygen rotating detonation rocket engine
Not provided.
SEARCH · Engineering Papers
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.
Not provided.
Azimuthal magnetic field measurements obtained during the implosion phase of an oxygen gas-puff Z-pinch on a 500 kA peak current and 180 ns rise time linear transformer driver are presented. While a fraction of the driver current was measured within the imploding plasma, key initial conditions were found to significantly impact the delivery of current to the plasma load. The electrode geometry was modified to assist the initial dielectric breakdown and resulted in improved shot reproducibility. Optimization of the gas injection plenum pressure and timing resulted in an increase in the current coupling parameter, defined as the ratio of the measured value of $B$ $θ$ to the expected value, from 50% to 75%. Furthermore, the degree of radial expansion of the gas puff in the load region, which is suspected to lead to the observed current loss during the implosion, was reduced by shortening the valve opening duration. Additionally, a pre-embedded axial magnetic field of up to 0.2 T was found to have no significant impact on the plasma-current coupling of the oxygen implosions.
We present simulations of an oxygen gas puff Z-pinch on a University scale generator at the Weizmann Institute of Science. The work accounts for the detailed geometry of the nozzle, the initial neutral gas density distribution, and the subsequent implosion. The modeling results show significant improvement with data for the current at the time of stagnation in comparison with a previous effort [Rosenzweig et al., Phys. Plasmas 27, 022705 (2020)]. As a first step, we performed simulations of the flow of neutral diatomic oxygen from a plenum through a nozzle within a recessed cathode, across a gap, and into the anode with a recessed grounded honeycomb. These simulations show an agreement with the measured initial gas density profiles within the region not blocked by the recesses and accessible to visible measurements. The computed neutral gas flow profile serves as the initial condition for a radiation magnetohydrodynamic simulation of the implosion using the MACH2-TCRE code. By considering the specific details of the nozzle and chamber geometry, we find agreement with the measured current profile, including the inductive notch. The simulations predict that the plasma undergoes a strong pinch within the hidden anode recess. The simulations also predict the strongest radiation pulse occurs within the anode recess and at the time of the observed inductive notch.
This paper presents the structural design and supporting analysis for the Microreactor Applications Research Validation and Evaluation (MARVEL) primary coolant system (PCS) using the American Society of Mechanical Engineers (ASME) Boiler and Pressure Vessel Code Section III, Division 5, rules. MARVEL is a liquid metal–cooled microreactor intended to provide experimental capabilities for the rapid testing and development of microreactor technologies. The PCS utilizes high-temperature sodium-potassium liquid metal as the primary coolant and operates at a design temperature of 570°C, necessitating the consideration of creep-related failure mechanisms. The base metal for the PCS is 316H stainless steel, and the weldments are made with a 16-8-2 filler. The design approach incorporates the current base code rules along with ASME code cases N-924, N-861, and N-862 to address primary load, ratcheting, and creep-fatigue evaluations, respectively. The reactor’s operation involves complex thermal and mechanical interactions due to natural convective flow and differential thermal expansion between components. In conclusion, this paper discusses the structural engineering challenges encountered, such as managing thermal stresses in the distribution plenum and guard vessel, and outlines the strategies implemented to meet the code requirements, including design modifications and operational constraints.
The Transformational Challenge Reactor is a 3-MW(thermal) helium-cooled experimental nuclear reactor designed using an additive manufacturing–informed agile design process. This design process leverages rapid prototyping and advanced materials from emerging additive manufacturing technologies, key characteristics that enable rapid design maturation. The resulting core design incorporates a blend of advanced reactor technologies into an intermediate-spectrum microreactor, including conventionally manufactured tristructural isotropic (TRISO) fuel particles in an advanced manufactured SiC fuel element and a solid yttrium hydride moderator encapsulated in steel. Matured during the design effort, these technologies are incorporated with additively manufactured steel support and fluidic structures to form a 75-cm-outer-diameter cylindrical active core region. Below and above the active core region are axial SiC reflectors, which are housed inside the reactor pressure vessel. The reactor is controlled with an annular shroud actuated external to the pressure vessel in the gap between the pressure vessel and a steel radial reflector. A safety rod is at the center of the core to shut down the reactor when necessary. Helium pressurized at 5 MPa is forced into the pressure vessel below the core and around the core to the top plenum before it is forced down through the axial reflectors and the active core region. The primary pressurized helium loop is operated up to 500°C and includes the pressure vessel, the circulator, and the hot side of a helium-to-air heat exchanger. The secondary loop rejects all heat from the primary loop to ambient air through a heat exchanger. A vented temporary confinement building contains the entire primary loop, with penetrations for a stack, cooling, and the secondary ambient air loop. Finally, this is the first advanced nuclear microreactor designed using additive manufacturing technologies, demonstrating their applicability in an accelerated advanced design process.
High-burnup and extended enrichment fuels are of interest for extending the cycle lengths of pressurized water reactors from 18 to 24 months. Changes to the fuel design and core loading scheme have potentially significant safety implications due to power distribution effects, reduced thermal conductivity, and/or increased plenum pressures due to additional burnable poison loadings. Additionally, higher burnups result in increased material degradation and risk of fuel fragmentation, relocation, and dispersal (FFRD). A representative 24-month core design using gadolinia-doped UO 2 was analyzed for performance under large-break (LB) loss-of-coolant accident (LOCA) conditions using PARCS, RELAP5-3D, and BISON. Furthermore, all considered acceptance criteria were met, with no cases exceeding the 1477 K maximum cladding temperature or the post-quench ductility oxidation limit of 17% equivalent cladding reacted. Full-core FFRD susceptibility was estimated to be approximately 455 kg, though high uncertainties exist with current approaches for computing susceptibility. Undoped fuel rods are more likely to be limiting due to higher linear heat rates. Relatively high burnup and linear heat rate rods located in second batch assemblies are of greatest safety significance during LB LOCA for this high-burnup core design.
The performance of the refractory lining in glass melters used for nuclear waste vitrification is critical to the melter reliability for long-term continuous operation. Monofrax® K-3, a high Cr 2 O 3 fused cast refractory material, has been widely used to build the liners of nuclear waste glass melters in the United States. Corrosion behaviour of Monofrax® K-3 refractory has been evaluated based on crucible-scale testing, inspection of the refractory components following scaled melter testing, and inspections of the Defense Waste Processing Facility (DWPF) melter refractory after service. The literature generally consists of empirical models based on short-term testing to describe refractory corrosion dependence on glass composition. Corrosion data from tests with longer testing times, at various temperatures, in the presence of molten salts, and with different redox reactions in the plenum atmosphere exist, may be insufficient to provide accurate refractory service life estimates. Additionally, the corrosion data collected under actual and scaled melter operating conditions are limited. Recommendations to achieve more direct correlation between the laboratory refractory corrosion data predictions and the observed melter service life are discussed to allow for more accurate predictions of the useful life of melter refractory linings.
Future fusion reactors must be able to breed the tritium they will consume. Several breeding and tritium extraction technologies are under investigation internationally. PbLi is of particular interest as a breeder material. A new forced convection PbLi loop, the Tritium Extraction eXperiment (TEX), is being designed and constructed to investigate tritium extraction from PbLi. Specifically, TEX will serve to investigate the extraction efficiency of various vacuum permeator membranes and configurations. The major components of the loop include a moving magnet pump, reverse permeator, furnace test section, analysis chamber, supply tank, and plenum. As significant hazards are present in such experiments, safety is an integral focus of the experimental design.
Effective helium confinement time, τ p * ,He , and its ratio with energy confinement time, τ p * ,He / τ E , are key metrics quantifying the suitability of fusion plasmas for a continuous burn. Comparisons in the DIII-D tokamak of discharges with suppression of edge localized modes (ELMs) by resonant magnetic field perturbations (RMPs) to the corresponding unperturbed ELMy discharges found that both of these metrics were strongly reduced, by a factor of approximately 2, after application of the RMPs. This reduction in τ p * ,He during RMP ELM suppression was observed in the plasma core, edge, and pumping plenum, where higher neutral helium concentration during RMPs was also measured. In conclusion, these findings provide evidence that future devices employing RMP ELM suppression may meet or even exceed the helium exhaust provided by the ELMs themselves, reducing helium ash, and thus maintaining high fusion gain.
An experimental discharge from the DIII-D tokamak is modeled using the SOLPS-ITER code suite and compared against measurements in the pumped and relatively closed upper divertor. Uncertainties of boundary plasma simulations are identified by attempting to match code inputs to experimental conditions, including iteratively solving transport coefficients to match upstream experimental profiles using varying quantities of core particle flux, different pumping models, and various assumptions of ion thermal transport. Simulated boundary conditions for particle injection at the core interface are shown to be relevant to the plasma solution at the divertor targets, even if upstream transport is modified so that plasma profiles are comparatively similar, although seperatrix density is not held constant. When upstream plasma profiles are matched to experimental measurements by varying diffusive transport coefficients, using either poloidally symmetric or ballooning structure, the model finds a majority of injected energy being transported radially off the computational domain, in conflict with experimental radiated power measurements and heat flux measurements at the divertor target. Imposing a maximum thermal diffusivity or radially shifting the experimental separatrix location of the fitted profiles to increase power conducted to the targets by increasing the upstream electron temperature does not significantly modify this result. Including a thermalizing plenum volume in the simulation domain is shown to maintain the experimental volumetric pumping rate without knowing the neutral energy distribution incident on the pump duct a priori. By modifying transport parameters to match different assumptions for ion temperature, downstream neutral pressure changes by more than a factor of two, suggesting that attention to ion thermal transport may be a critical parameter for simulations to accurately resolve recycling and neutral transport, particularly in a closed divertor geometry. In addition to quantifying various modeling uncertainties, this work motivates both further experimental study and modeling improvements to improve predictive capabilities.
Tracer gas experiments were conducted in a 158 m 3 room with overhead supply diffusers to study dispersion of contaminants from simulated speaking in physically distanced meeting and classroom configurations. The room was contained within a 237 m 3 cell with open plenum return to the HVAC system. Heated manikins at desks and a researcher operating the tracer release apparatus presented 8-9 thermal plumes. Experiments were conducted under conditions of no forced air and neutral, cooled, or heated air supplied at 980-1100 cmh, and with/out 20% outdoor air. CO 2 was released at the head of one manikin in each experiment to simulate small (< 5 µm diameter) respiratory aerosols. The metric of exposure relative to perfectly mixed (ERM) is introduced to quantify impacts, based on measurements at manikin heads and at three heights in the center and corners of the room. Chilled or neutral supply air provided good mixing with ERMs close to one. Thermal stratification during heating produced higher ERMs at most manikins: 25% were ≥ 2.5 and the highest were > 5× perfectly mixed conditions. Operation of two within-zone air cleaners together moving ≥400 cmh vertically in the room provided enough mixing to mitigate elevated exposure variations.
Slurry feed charged into an electric melter creates a layer of reacting and melting material (termed cold cap) that floats on the surface of molten glass. The rheological behavior of heated melter feed affects the spreading of slurry at the top of the cold cap and the stability of the primary foam, affecting cold-cap coverage and melter plenum temperatures. The apparent viscosity of a high-alumina high-level waste melter feed was assessed by thermomechanical analysis, high-temperature viscometer, and the hot stage microscopy method, yielding viscosity estimates from ≈10 7.5 Pa s at 550°C to ≈10 2.5 Pa s at 1050°C. As the temperature of feed materials increased, their state changed from rigid solid to dilatant fluid, to pseudoplastic bubbly liquid with dissolving solids, to fully developed foam, and finally to Newtonian glass melt.
In the present work, a novel computational fluid dynamics (CFD) methodology was developed to simulate full-scale non-premixed rotating detonation engines (RDEs). A unique feature of the modeling approach was the incorporation of adaptive mesh refinement (AMR) to achieve a good trade-off between model accuracy and computational expense. Here, unsteady Reynolds-averaged Navier–Stokes (RANS) simulations were performed for an Air Force Research Laboratory (AFRL) non-premixed RDE configuration with hydrogen as fuel and air as the oxidizer. The finite-rate chemistry model, along with a ten-species detailed kinetic mechanism, was employed to describe the H 2 -Air combustion chemistry. Three distinct operating conditions were simulated, corresponding to the same global equivalence ratio of unity but different fuel/air mass flowrates. For all conditions, the capability of the model to capture essential detonation wave dynamics was assessed. An exhaustive verification and validation study was performed against experimental data in terms of a number of waves, wave frequency, wave height, reactant fill height, oblique shock angle, axial pressure distribution in the channel, and fuel/air plenum pressure. The CFD model was demonstrated to accurately predict the sensitivity of these wave characteristics to the operating conditions, both qualitatively and quantitatively. A comprehensive heat release analysis was also conducted to quantify detonative versus deflagrative burning for the three simulated cases. The present CFD model offers a potential capability to perform rapid design space exploration and/or performance optimization studies for realistic full-scale RDE configurations.
Downstream vortex generators that involve a pair of rectangular plates arranged in an open V-shape placed just downstream of each film-cooling hole were shown to create flow and vortical structures that entrain lifted film-cooling flow back to the surface and increase its lateral spreading on the surface (GT2020–14317). In this study, computations and measurements were performed to examine the flow mechanisms this vortex generator induces to improve film-cooling effectiveness of a flat plate with the cooling flow emanating from one row of inclined holes. Parameters studied include blowing ratio (BR = 0.75 and 1.0) and temperature ratio (TR = 1.07 and 1.9). Here, the computational study is based on steady Reynolds-averaged Navier–Stokes (RANS) closed by the shear-stress transport (SST) turbulence model with and without conjugate analysis. The experimental study was conducted by using a conjugate heat transfer test rig with a plenum, where cooling flow is introduced. Measurements made include velocity and temperature profiles upstream and downstream of the film-cooling holes as well as the temperature at several locations on the hot and cold sides of the film-cooled flat plate. The computational study was validated by comparing computed results with those from measurements at BR = 0.75 and 1.0 and TR = 1.9. Computational and experimental results are presented to show the effects of BR and TR on the flow structures and how those structures improve the effectiveness of film cooling with and without the downstream vortex generators and with and without conjugate heat transfer.
In the present work, a first-of-its-kind three-dimensional (3D) large-eddy simulation (LES) study is conducted to numerically investigate the combustion dynamics as well as aero-thermal phenomena in a full-scale nonpremixed hydrogen–air rotating detonation engine (RDE) (with a diverging-shaped lower-end wall), when integrated with nozzle guide vanes (NGV) acting as the turbine stator. The wall-modeled LES framework incorporates hydrogen–air detailed chemical kinetics and adaptive mesh refinement (AMR). A comparative analysis is carried out for two operating conditions with different fuel/air mass flow rates but global equivalence ratio of unity, and considering RDE configurations without and with stator. The LES model is validated against available experimental data for the low mass flux condition with respect to detonation wave speed/height, wave dynamics, and axial static pressure distribution. Numerical results indicate significant deflagrative combustion occurring in the fill region near the inner wall due to formation of recirculation zones in the injection near-field driven by the backward facing step. The leading detonation wave is found to be trailed by an azimuthal reflected-shock combustion (ARSC) wave, consistent with experimental observations, which consumes unburned vitiated reactants that leak through the main detonation wave. The main detonation wave characteristics, such as detonation wave speed/height and combustion efficiency, do not change appreciably with the presence of NGV. A novel combustion diagnostic technique based on chemical explosive mode analysis (CEMA) is employed to quantify the fraction of heat release occurring in the detonative mode versus deflagrative mode for the simulated conditions. The exit flow is found to be nearly fully subsonic and supersonic for the low and high mass flux conditions, respectively. Further analysis of the exit flow profiles shows that the presence of NGV renders the flow more axial and significantly impacts the exit Mach number and total pressure, while the total temperature shows negligible change. In addition, the low mass flux operating point, despite exhibiting more deflagrative losses within the combustor, yields overall lower pressure drop from plenum to exhaust, which is mainly attributed to lower pressure drop across the injectors. Lastly, the rotating detonation engine-nozzle guide vanes (RDE-NGV) configuration exhibits higher total pressure loss compared to rotating detonation engine (RDE) without stator across both the mass flux conditions. In conclusion, this study extends the state-of-the-art in numerical modeling of pressure gain combustion (PGC) systems by demonstrating high-fidelity 3D reacting LES of full-scale RDE-NGV systems relevant to RDE-turbine integration for stationary power generation.
A plethora of theoretical, numerical, and experimental investigations have relied on updraft forced convection of operating fuel assemblies for light water-cooled fuel assemblies. The inertial scales of turbulence in updraft flow through PWR fuel assemblies are primarily influenced by upstream mixing patterns within the lower plenum and consequent acceleration of flow through the lower core and support plates into the bottom nozzles. In secondary heat transfer applications such as steam generation, recuperation and economization, downdraft flow is utilized to retrieve maximal sensible heat from the primary coolant. Sub-channel volumes aligned vertically could also benefit from force of gravity in developing turbulent flows. A large reduction in pumping power, associated cost of operation and maintenance, and improved balance of plant is suggested to be possible through downdraft forced convection in PWR-type reactors. This article explores the impact of downdraft forced convection within the well-studied 3x3 sub-channel within a range of 20 hydraulic diameters across a representative intermediate flow mixing grid. A range of inflow Reynolds number from 10,000 to 60,000 is simulated over a range of equivalent heat loads of 20% to 100% spanning 20 hydraulic diameters. The fuel rods are represented in true geometric detail, and a constant-value heat profile is assumed for the investigation. At each inflow condition, the turbulence intensity, peak velocity, local heat transfer coefficient, peak wall heat flux, and peak wall temperature are collected across the sub-channel and compared with an updraft configuration evaluated at the same conditions.
Potential iodine gas release from failed fuel pins is a critical factor in the source term analysis of oxide fuel-loaded sodium fast reactors (SFRs). The accumulated iodine-containing gas mixtures inside pin plenums are expected to be ejected during pin failures and rise through sodium pool, with potential release of gaseous iodine to the cover gas region. Due to its potential radiological impacts, a proper assessment of iodine behavior is necessary for an accurate source term assessment. Throughout the bubble rise trajectory in the sodium pool, iodine gas is continuously removed or transformed at the bubble interface by diffusion, as the combining reaction between the iodine and sodium to form sodium iodide (NaI) is a chemically preferred process. As the final amount of iodine released from the facility is strongly influenced by the removal phenomenon inside the sodium pool, experiments were previously performed by PNC (Power Reactor and Nuclear Fuel Development Corporation) to provide insight into this phenomenon. To assess the accuracy of present methods for predicting iodine gas removal within sodium pools, several candidate approaches, available in source term analysis codes, have been summarized and evaluated in this study. Spherical cap bubbles and spherical bubbles are considered in accordance with the methods adopted in each approach, and different forms of correlations for major parameters have been implemented in accordance with the original adoptions. Based on the summarized results, important aspects to be considered have been derived.
Uranium–molybdenum (U–Mo) alloys in monolithic fuel foil are the primary candidate for the conversion of high-performance research reactors in the USA. Monolithic fuel is utilized in a plate-type design with a zirconium diffusion barrier and aluminum cladding. These fuel types are unique in that they contain no plenum for the release of fission gases, which, in conjunction with the aluminum cladding, can lead to large stress states within the fuel. The nature of how fundamental processes of radiation damage, including the evolution of point defects, under such stresses occur is unknown. In this work, we present molecular dynamics simulations of the formation energy of point defects under applied stress. We report this work will allow for the implementation of stress-dependent microstructural evolution models of nuclear fuels, including those for both fission gas bubble growth and creep, which are critical to ensure the stable and predictable behavior of research reactor fuels.