Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “power gas simulation”

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 1,081 records · Page 60
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Cambium 2023 Scenario Descriptions and Documentation

The National Renewable Energy Laboratory's (NREL's) Cambium data sets are annually released sets of simulated hourly emission, cost, and operational data for a range of modeled futures of the U.S. electric sector with metrics designed to be useful for long-term decision- making. The 2023 Cambium data set is the fourth annual release. The data sets are a companion product to NREL's Standard Scenarios, which are likewise released annually and are a set of projections of how the U.S. electric sector could evolve across a suite of different potential futures, but covering more scenarios with less temporal granularity (Gagnon et al. 2024). Information about Cambium and related publications can be found at https://www.nrel.gov/analysis/cambium.html, and the Cambium data sets can be viewed and downloaded at https://scenarioviewer.nrel.gov/. In this documentation, we describe Cambium 2023's scenarios, define the metrics, and document the Cambium-specific methods for calculating those metrics.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS2004

Epoxy Adhesives for Stator Magnet Assembly in Stirling Radioisotope Generators (SRG)

As NASA seeks to fulfill its goals of exploration and understanding through missions planned to visit the moons of Saturn and beyond, a number of challenges arise from the idea of deep space flight. One of the first problems associated with deep space travel is electrical power production for systems on the spacecraft. Conventional methods such as solar power are not practical because efficiency decreases substantially as the craft moves away from the Sun. The criterion for power generation during deep space missions are very specific, the main points requiring high reliability, low mass, minimal vibration and a long lifespan. A Stirling generator, although fairly old in concept, is considered to be a potential solution for electrical power generation for deep space flight. A Stirling generator works on the same electromagnetic principles of a standard generator, using the linear motion of the alternator through the stationary stator which produces electric induction. The motion of the alternator, however, is produced by the heating and cooling dynamics of pressurized gases. Essentially heating one end and cooling another of a contained gas will cause a periodic expansion and compression of the gas from one side to the other, which a displacer translates into linear mechanical motion. NASA needs to confirm that the materials used in the generator will be able to withstand the rigors of space and the life expectancy of the mission. I am working on the verification of the epoxy adhesives used to bond magnets to the steel lamination stack to complete the stator; in terms of in-service performance and durability under various space environments. Understanding the proper curing conditions, high temperature properties, and degassing problems as well as production difficulties are crucial to the long term success of the generator. system and steel substrate used in the stator. To optimize the curing conditions of the epoxies, modulated differential scanning calorimetry analysis was done as a function of cure time and temperatures. Adhesion bond strength was tested at various temperatures with lap shear samples using Hiperco 50 substrate to ensure that the proper adhesive is being used. To try and solve the problem of bondline thickness, micro glass beads of 0.0017" in diameter were investigated to see if any other physical properties of the epoxy were affected. Efforts will be made to develop a standard, optimized, fabrication process/procedure of sub-scale magnet-stator assemblies for various adhesive performance evaluation studies under simulated generator conditions. Also, accelerated aging testing will be done in a pressurized canister with stator assembly samples for three years to verify if any degassing or thermal degradation of the epoxy occurs. The necessity of verifying the correct epoxy adhesive system for the stator magnet in the SRG is crucial because failure of the stator assembly would jeopardize the electrical system, and thereby the entire mission itself. My work involves specimen fabrications, testing, and data analyses of the epoxy adhesive system for the stator magnet in the SRG is crucial because failure of the stator assembly would jeopardize the electrical system, and thereby the entire mission itself.

Cater, George M.↗
DOE OSTI2025

Thermal-Fluid and Thermal-Structural Response of the T-Tube Modular Divertor to Spatiotemporally Varying Heat Loads

Tungsten (W) is the leading candidate for divertor target plates because of its high melting point (>3000°C), thermal conductivity, and ultimate tensile stress. While W and its alloys are the only solid materials that can survive the high heat fluxes incident on the divertor, W’s low-ductility high ductile-to-brittle transition temperature of ~600°C and relatively low recrystallization temperature (RT) of ~1300°C pose structural (among other) challenges. The objective of this work is to estimate the thermal-fluid and thermal-structural performance of the helium (He)-cooled T-tube divertor, which was originally developed by the Advanced Reactor Innovation and Evaluation Study (ARIES) using numerical simulations. Here, predictions of temperature distributions across the plasma-facing structural component and surface pressures from computational fluid dynamics simulations are used to determine stress distributions using commercial structural finite element modeling software over a range of fusion-relevant conditions. The maximum allowable incident heat fluxes are determined based on the temperature limits imposed by the ITER elastic Structural Design Criteria for In-vessel Components (SDC-IC) and the maximum RT over a range of He mass flow rates and presented in the form of performance design charts. Our recent work found that thermal- structural criteria accounting for the low ductility of W in a finger-type modular divertor constrain the maximum incident heat fluxes to values well below the ITER specifications, and those based on considering only the RT demonstrate that integrated thermal-fluid and elastic structural performance evaluation are required for accurate assessment of divertor performance. This novel analysis of the T-tube considers how nonuniform and transient incident heat fluxes affect its thermal-fluid and thermal-structural performance, as well as the effect of volumetric heating, which can be as great as 27% of the power incident on the divertor surface. The W tile of the T-tube, with its relatively large plasma-facing area of ~15 cm 2 , will likely experience significant spatial variations in incident heat flux. This work therefore assesses whether steady-state incident heat flux profiles with a peak of 10 MW/m 2 and maximum heat flux gradients of 200 MW/m 2 per m exceed the structural limits imposed by the ITER elastic SDC-IC and the maximum RT over a range of fusion-relevant conditions. The effect of transient heat fluxes typical of plasma detachment and reattachment from the target plate due, for example, to gas injection are also evaluated

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Computational Evaluation of Carriers in Facilitated Transport Membranes for Postcombustion Carbon Capture

Incorporating molecular amines as mobile carriers in facilitated transport membranes (FTMs) has been demonstrated to significantly enhance the CO 2 permeance and CO 2 /N 2 selectivity of the membrane for CO 2 capture from flue gas. In this study, by employing computational techniques including density functional theory calculations and molecular simulations, the role of mobile carriers has been systematically studied at a molecular level from the perspectives of the amine–CO 2 reaction chemistry, diffusivities of carriers and gases, and N 2 solubility. The latter two properties were also investigated as a function of water uptake. The water uptake values of FTMs were experimentally quantified too. The introduction of mobile carriers was shown to substantially enhance the diffusivities of CO 2 reaction products compared to FTMs without mobile carriers. The choice of mobile carriers was also demonstrated to influence the separation performance. Computationally, 2-(1-piperazinyl)ethylamine sarcosinate (PZEA-Sar) exhibited a faster reaction kinetics and slightly higher CO 2 absorption capacity as compared to piperazine glycinate (PZ-Gly). Experimentally, the FTM incorporating PZEA-Sar mobile carriers also showed a higher CO 2 permeance. The good agreement validated the computational models employed and insights generated in this study. In conclusion, the outcomes of this work shed light on the future design and selection of carrier structures, and the adopted computational approaches can be employed to discover promising mobile carrier candidates.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2023

Scale-up Testing of Advanced Polaris Membrane in CO2 Capture Technology

This final technical report describes work conducted by Membrane Technology and Research, Inc. (MTR) for the U.S. Department of Energy (DOE), National Energy Technology Lab (NETL) on the scale-up and testing of advanced Polaris™ membrane CO2 capture technology at the Technology Centre Mongstad (TCM) under award number DE-FE0031591. The work was performed from August 1, 2018 through January 31, 2023. The overall goal of this project was to design, build and operate an advanced Polaris membrane CO2 capture system at TCM. MTR was assisted in this project by Trimeric Corporation (Trimeric), an engineering design services company, the Carbon Capture Simulation for Industry Impact (CCSI2), a partnership among national laboratories, industry, and academic institutions, and the Technology Centre Mongstad (TCM), who provided the host site for the slipstream field test. This report details the work conducted to scale-up MTR’s second-generation (Gen-2) Polaris membrane and advanced planar membrane modules to a final form factor optimized for commercial use; validate their performance in an engineering-scale field test at TCM; and to show the potential of the MTR process to meet DOE CO2 capture targets from large source point emitters. Work for this project included membrane optimization and scale-up, advanced planar module design and fabrication, design and fabrication of an engineering-scale field test membrane skid, operation of the field test skid processing Residue Fluid Catalytic Cracker (RFCC) industrial flue gas at TCM, and a detailed techno-economic analysis (TEA) of the MTR membrane post-combustion process for CO2 capture. This project validated recent membrane technology advancements at the engineering-scale, moves the MTR advanced post-combustion capture technology to TRL-6, and mitigates risk in future Large Pilot or Demonstration scale-up activities.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2025

Transformational Nano-confined Ionic Liquid Membrane for Greater than or Equal to 97 Percent Carbon Dioxide Capture from Natural Gas Combined Cycle Flue Gas

A transformational process based on nano-confined ionic liquid (NCIL) membranes was developed for capturing ≥97% CO 2 from natural gas combined cycle (NCCC) flue gas. The NCIL membranes were prepared by loading amino acid ionic liquid into a framework composed of single-walled carbon nanotube mesh filled with graphene oxide quantum dots. The membranes exhibited CO 2 permeance as high as 2,000 GPU with a CO 2 /N2 selectivity of 2,300 for a typical NGCC flue gas composition. When H 2 O vapor sweep was applied in the permeate side, 96.6% CO2 dry-basis purity and 97.6% CO 2 capture rate were achieved for a simulated NGCC flue gas with single stage. In the process design, a highly H 2 O-selective membrane would be needed to recover majority of the H 2 O vapor, and the recovered H 2 O vapor could be recycled to the permeate side of the NCIL membrane. Sulfonated poly(ether ether ketone) membranes were successfully developed for this purpose. These membranes exhibited H 2 O permeance great than 11,000 GPU and H 2 O/CO 2 selectivity greater than 1,000 at 70ºC for a feed mixture consisting of 14.5 vol% H2O and balanced CO 2 . A standalone membrane model using MATLAB platform was developed for process simulation. The model was validated with experimental data. Techno-economic analysis based on the testing data collected during the current program suggests this transformational membrane process can achieve 97% CO 2 capture efficiency with a cost of $47.8/tonne of CO 2 , which is a 21% reduction versus DOE’s reference case B31B.97.

03 NATURAL GAS↗
DOE OSTI2019

Experimental and Computational Analysis of NEAMS Pebble Bed Reactors

The Advanced High-Temperature Pebble Bed Reactor concept leverages a particle-based fuel format consisting of discrete spherical graphite pebbles arrayed in a packed bed architecture. Thermal regulation achieved via flow of gas (e.g., helium) or liquid (e.g., molten salt) coolants through the void spaces between pebbles in the bed (characteristic pebble diameters: ~ 6 cm, gas cooled; ~ 3 cm, liquid cooled). This design is of interest owing to inherent advantages including passive safety features and highly efficient heat transfer characteristics that enable high power densities to be achieved. Fully leveraging these advantages requires a fundamental understanding of the complex coolant flow structure within the pebble-shaped fuel bed. Unfortunately, complexities associated with non-invasively probing the micro-scale flow phenomena within the complex network of void spaces between pebbles have hindered efforts to characterize the underlying transport phenomena. Without this fundamental foundation, it becomes challenging to rationally construct macroscopic transport models based upon volume averages of micro-scale parameters that accurately capture the observed flow behavior. Consequently, this lack of fundamental knowledge has limited previous work to coarse comparisons between theory and extremely limited experimental data. Moreover, these comparisons have been made only in terms of integral or macroscopic results. Therefore, multipoint measurements with a high level of spatial and temporal resolution are critically needed to fully map the complex flow patterns and to provide data at high spatial density to permit accurate volume averaging in the pebble bed. For this project, we performed a coordinated experimental and computational effort to quantitatively map the full-field velocity and temperature fields in the interstitial spaces within the pebble bed. The major accomplishments of this project included the important knowledge acquired for the flow mixing and heat transfer in the interstitial pore scales, the effects of flow mixing in the near-wall and far-wall regions, and the by-pass flows of pebble bed systems. The produced high-resolution experimental data (not broadly available in existing literature) were then used for the numerical studies enabling the refinement and validation of standard simulation tools (Nek5000).

42 ENGINEERING↗
NASA NTRS2010

Ballistics Model for Particles on a Horizontal Plane in a Vacuum Propelled by a Vertically Impinging Gas Jet

A simple trajectory model has been developed and is presented. The particle trajectory path is estimated by computing the vertical position as a function of the horizontal position using a constant horizontal velocity and a vertical acceleration approximated as a power law. The vertical particle position is then found by solving the differential equation of motion using a double integral of vertical acceleration divided by the square of the horizontal velocity, integrated over the horizontal position. The input parameters are: x(sub 0) and y(sub 0), the initial particle starting point; the derivative of the trajectory at x(sub 0) and y(sub 0), s(sub 0) = s(x(sub 0))= dx(y)/dy conditional expectation y = y((sub 0); and b where bx(sub 0)/y(sub 0) is the final trajectory angle before gravity pulls the particle down. The final parameter v(sub 0) is an approximation to a constant horizontal velocity. This model is time independent, providing vertical position x as a function of horizontal distance y: x(y) = (x(sub 0) + s(sub 0) (y-y(sub 0))) + bx(sub 0) -(s(sub 0)y(sub 0) ((y - y(sub 0)/y(sub 0) - ln((y/y(sub 0)))-((g(y-y(sub 0)(exp 2))/ 2((v(sub 0)(exp 2). The first term on the right in the above equation is due to simple ballistics and a spherically expanding gas so that the trajectory is a straight line intersecting (0,0), which is the point at the center of the gas impingement on the surface. The second term on the right is due to vertical acceleration, which may be positive or negative. The last term on the right is the gravity term, which for a particle with velocities less than escape velocity will eventually bring the particle back to the ground. The parameters b, s(sub 0), and in some cases v(sub 0), are taken from an interpolation of similar parameters determined from a CFD simulation matrix, coupled with complete particle trajectory simulations.

Lane, J. E.↗
DOE OSTI2024

Development of an Improved RELAP5-3D Model for the High Temperature Test Facility

High-temperature gas-cooled reactors (HTGRs) are rapidly approaching deployment. Confidence in transient analysis of these systems for design, optimization, and licensing calculations requires modeling and simulation tools that have been validated against data relevant to HTGR conditions. The High Temperature Test Facility (HTTF) is an integral effects thermal hydraulics test facility for prismatic HTGRs. In spring and summer of 2019, HTTF was used for a series of experiments that now serve as the basis for the OECD/NEA Thermal Hydraulic Code Validation Benchmark for High Temperature Gas-Cooled Reactors using HTTF Data (HTGR T/H Benchmark). This benchmark contains problems for systems code, computational fluid dynamics (CFD), and coupled systems code/CFD modeling representing lower plenum mixing and both the depressurized and pressurized conduction cooldown (DCC and PCC respectively) transients. Benchmark problems include exercises for code-to-code and code-to-data comparisons as well as an exercise for error scaling between HTTF and the Modular High Temperature Gas-Cooled Reactor, which serves as the basis for the HTTF design. Previous analysis as part of the HTGR T/H benchmark used a RELAP5-3D model developed at Idaho National Laboratory (INL) and demonstrated an ability to reproduce trends in the measured data but difficulties reproducing experimental values within their uncertainty. These difficulties were largely attributed to assumptions made during the development of the initial RELAP5-3D model, which predated the HTTF experiments. A significant cause of difficulty reproducing the measured temperatures may be the radial nodalization of the previous RELAP5-3D model. The new model provides a finer nodalization to assess the impact of radial nodalization and allows for asymmetric heating within the core, which was a feature of multiple HTTF experiments. In this paper, we present the new RELAP5-3D model of HTTF. In addition to describing the new model, this paper compares the new and old models and provides results for a full-power steady state, a DCC, and a PCC in HTTF. These analyses are based on the code-to-code comparison exercises for the DCC and PCC problems of the HTGR T/H benchmark. We present the results of these exercises from the new model and compare them to the results of the old model.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Xenon Trapping in Metal-Supported Silica Nanocages

Xenon (Xe) is a valuable and scarce noble gas used in various applications, including lighting, electronics, and anesthetics, among many others. It is also a volatile byproduct of the nuclear fission of uranium. In this work, a novel material architecture consisting of silicate nanocages in contact with a metal surface and an approach for trapping single Xe atoms in these cages is presented. The trapping is done at low Xe pressures and temperatures between 400 and 600 K, and the process is monitored in situ using synchrotron-based ambient pressure X-ray photoelectron spectroscopy. Release of the Xe from the cages occurs only when heating to temperatures above 750 K. A model that explains the experimental trapping kinetics is proposed and tested using Monte Carlo methods. Density functional theory calculations show activation energies for Xe exiting the cages consistent with experiments. This work can have significant implications in various fields, including Xe production, nuclear power, nuclear waste remediation, and nonproliferation of nuclear weapons. The results are also expected to apply to argon, krypton, and radon, opening an even more comprehensive range of applications.

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

Temperature sensitivity of the equilibrium neutronics and accident analysis of the HTR-10

Pebble-Bed High-Temperature Gas-cooled Reactors (PB-HTGR) are moderated by the graphite in the fuel pebbles and the graphite reflector surrounding the pebble-bed. Because graphite is by far the most abundant material in PB-HTGRs and the primary moderator, accurate modeling of the graphite material, including density, impurities, and temperatures, is crucial for accurate computational modeling and simulation of these reactors. While main characteristics of the graphite components are often known, the local temperature is less well known and often averaged over all components. Here, this work studies the impact of considering accurate temperature profiles in the graphite material on the generation of a small PB-HTGR model at the state of equilibrium operation and on short-term accident progression. The fuel compositions for the PB-HTGR were determined using a jump-in equilibrium modeling method, the Axial Radial Zone Equilibrium Modeling (AR-ZEM) method. In contrast to previous work, the AR-ZEM method was used considering thermal-hydraulic feedback from the MELCOR code to determine temperatures of the fuel pebbles and the surrounding graphite reflector. The consideration of an axial and radial temperature profile in the core and reflector, as opposed to uniform material temperatures, had an impact of almost 1,300 pcm on the equilibrium core eigenvalue and caused significant differences in the discharged plutonium fuel inventory with up to 4.9% and 11.0% for Pu-239 and Pu-242, respectively. To assess the impact on short-term accident progression, two Anticipated Transient Without SCRAM (ATWS) events, a Pressurized Loss of Forced Coolant (PLOFC) and a Control Rod Withdrawal (CRW) with loss of flow, were simulated with MELCOR. The use of temperature profiles in the equilibrium core models did not reveal a significant impact on the temperature, power, or reactivity responses during the transients. In conclusion, a need for consideration of accurate temperature profiles, in particular for the graphite reflector, was found for the generation of equilibrium PB-HTGRs core models using jump-in methods, but detailed temperature profiles may not be necessary when performing conservative transient analysis.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
NASA NTRS2001

Optical Diagnosis of Gas Turbine Combustors Being Conducted

Researchers at the NASA Glenn Research Center, in collaboration with industry, are reducing gas turbine engine emissions by studying visually the air-fuel interactions and combustion processes in combustors. This is especially critical for next generation engines that, in order to be more fuel-efficient, operate at higher temperatures and pressures than the current fleet engines. Optically based experiments were conducted in support of the Ultra-Efficient Engine Technology program in Glenn's unique, world-class, advanced subsonic combustion rig (ASCR) facility. The ASCR can supply air and jet fuel at the flow rates, temperatures, and pressures that simulate the conditions expected in the combustors of high-performance, civilian aircraft engines. In addition, this facility is large enough to support true sectors ("pie" slices of a full annular combustor). Sectors enable one to test true shapes rather than rectangular approximations of the actual hardware. Therefore, there is no compromise to actual engine geometry. A schematic drawing of the sector test stand is shown. The test hardware is mounted just upstream of the instrumentation section. The test stand can accommodate hardware up to 0.76-m diameter by 1.2-m long; thus sectors or small full annular combustors can be examined in this facility. Planar (two-dimensional) imaging using laser-induced fluorescence and Mie scattering, chemiluminescence, and video imagery were obtained for a variety of engine cycle conditions. The hardware tested was a double annular sector (two adjacent fuel injectors aligned radially) representing approximately 15 of a full annular combustor. An example of the two-dimensional data obtained for this configuration is also shown. The fluorescence data show the location of fuel and hydroxyl radical (OH) along the centerline of the fuel injectors. The chemiluminescence data show C2 within the total observable volume. The top row of this figure shows images obtained at an engine low-power condition, and the bottom row shows data from a higher power operating point. The data show distinctly the differences in flame structure between low-power and high-power engine conditions, in both location and amount of species produced (OH, C2) or consumed (fuel). The unique capability of the facility coupled with its optical accessibility helps to eliminate the need for high-pressure performance extrapolations. Tests such as described here have been used successfully to assess the performance of fuel-injection concepts and to modify those designs, if needed.

Hicks, Yolanda R.↗
DOE OSTI2020

Development of the INL Thermal Energy Distribution System (TEDS) in the Modelica Eco-System for Validation and Verification

This report provides an overview on the development of the Thermal Energy Distribution System (TEDS) in the Modelica process model eco-system at Idaho National Laboratory as part of the integrated energy system initiative. Model development has led to the creation of a dynamic process model of the experimental TEDS facility housed within the Engineering Science Laboratory (ESL) at INL. The model was then used during the pre-construction phase of the experimental effort to inform experimental design (insulation requirements, bypass line placement, expected performance of components) and to test innovative control schemes prior to the initial operation. The TEDS model developed in Modelica includes the primary components of the TEDS experimental unit including: a 200kW Chromalox heater, a single-tank packed bed thermal energy storage system filled with 0.125 inch Alumina (Al 2 O 3 ), an ethylene-glycol to Therminol-66 heat exchanger, system piping, six control valves, and all associated temperature, pressure, and mass flow sensors. The model does not include nitrogen fill gas tanks, or associated overfill tanks, as these are not part of standard system control. Two simulation sets were run. The first is a five-hour test that operates as a shakedown test for the facility. This simulation puts the facility through all five potential operating modes and showcases the ability of valving, control sensors, and component controllers to meet the system demands. The second case imposes a typical summer day demand on the system from a region with mixed commercial and residential electrical needs where the generator alone cannot meet peak demand but instead requires the thermal storage unit to act as a peaking unit. Over the course of the fifteen-hour discharge cycle heat losses account for a heat loss of approximately 15°C in outlet temperature. Such heat losses would lead to a decrease in overall system efficiency and may lead to additional design changes for such designs in the future that may include a topping heater. Through commencement of this work, a systems level model of TEDS with associated control systems, sensors, piping diameters, and component capabilities has been created. This model has been utilized in the pre-experimental phase to inform system design, insulation thicknesses, and potential control schemes to operate the system effectively and safely. Once TEDS begins operation, this model will be refined, tuned, and used for validation and verification purposes.

24 POWER TRANSMISSION AND DISTRIBUTION↗
NASA NTRS2016

Generalized Fluid System Simulation Program, Version 6.0

The Generalized Fluid System Simulation Program (GFSSP) is a general purpose computer program for analyzing steady state and time-dependent flow rates, pressures, temperatures, and concentrations in a complex flow network. The program is capable of modeling real fluids with phase changes, compressibility, mixture thermodynamics, conjugate heat transfer between solid and fluid, fluid transients, pumps, compressors, and external body forces such as gravity and centrifugal. The thermofluid system to be analyzed is discretized into nodes, branches, and conductors. The scalar properties such as pressure, temperature, and concentrations are calculated at nodes. Mass flow rates and heat transfer rates are computed in branches and conductors. The graphical user interface allows users to build their models using the 'point, drag, and click' method; the users can also run their models and post-process the results in the same environment. Two thermodynamic property programs (GASP/WASP and GASPAK) provide required thermodynamic and thermophysical properties for 36 fluids: helium, methane, neon, nitrogen, carbon monoxide, oxygen, argon, carbon dioxide, fluorine, hydrogen, parahydrogen, water, kerosene (RP-1), isobutene, butane, deuterium, ethane, ethylene, hydrogen sulfide, krypton, propane, xenon, R-11, R-12, R-22, R-32, R-123, R-124, R-125, R-134A, R-152A, nitrogen trifluoride, ammonia, hydrogen peroxide, and air. The program also provides the options of using any incompressible fluid with constant density and viscosity or ideal gas. The users can also supply property tables for fluids that are not in the library. Twenty-four different resistance/source options are provided for modeling momentum sources or sinks in the branches. These options include pipe flow, flow through a restriction, noncircular duct, pipe flow with entrance and/or exit losses, thin sharp orifice, thick orifice, square edge reduction, square edge expansion, rotating annular duct, rotating radial duct, labyrinth seal, parallel plates, common fittings and valves, pump characteristics, pump power, valve with a given loss coefficient, Joule-Thompson device, control valve, heat exchanger core, parallel tube, and compressible orifice. The program has the provision of including additional resistance options through User Subroutines. GFSSP employs a finite volume formulation of mass, momentum, and energy conservation equations in conjunction with the thermodynamic equations of state for real fluids as well as energy conservation equations for the solid. The system of equations describing the fluid network is solved by a hybrid numerical method that is a combination of the Newton-Raphson and successive substitution methods. The application and verification of the code has been demonstrated through 30 example problems.

Majumdar, A. K.↗
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↗
NASA NTRSDate not supplied

Large Vehicle Lunar Landing Surface Interaction and In-Situ Resource Based Risk Mitigation: Landing & Launch Pads

A key capability required for the exploration of planetary bodies is the ability to land on the surface. Previous work performed by NASA and other institutions has primarily focused on landing small spacecraft on planetary surfaces and the associated small-to-medium thrusters required for the soft landing. In the case of human exploration—particularly the establishment of long duration exploration and habitation outposts—the ability to land large landers, such as the SpaceX Starship, is necessary. These larger landing systems require the use of more powerful engines, with higher engine exhaust temperatures and higher landing loads. Understanding the excavation of material by the engines, as well as the potential for the landing legs to sink into the subsurface, is key in ensuring reliable and safe landings. A further improvement in landing reliability can be achieved by constructing landing / launch pads, especially with in-situ resources. Some material excavation by the plume is inevitable, leaving at least a portion of the surface scoured and uneven under the lander and ejecting regolith particles and rocks at very high velocities. One possible solution would be to robotically build landing / launch pads (ideally autonomously) at the destination using in-situ materials. In this case, the first one or few landers will need to land on unimproved surfaces at higher risk; however, they would bring the required equipment to build the landing pads with mostly local resources, thus increasing the reliability of safe landing for subsequent larger landers. A number of methods to build in-situ landing and launch pads have already been developed. These methods include no, or some, addition of required binder additives to the local regolith material, different processing approaches and result in varying landing pad strengths. A sub-scale rocket engine plume, was used to simulate some of the conditions of a landing on the Moon to assess the effectiveness of various materials for an in-situ built landing pad, The GO2/GCH4rocket engine fired on a 1m2area coupons of representative pad materials. The results will allow continued development towards materials that satisfy the landing pad properties required for the effective risk reduction and increased reliability for landing people and equipment on the lunar surface. This work contained two parts: (1) computer modeling of a large rocket engine plume interacting with regolith on the Moon, using the Granular Gas Flow Solver (GGFS) provided by CFD Research Corporation as well as other computational fluid dynamics codes (CFD) such as Loci/CHEM. (2) Developing landing/launch pad materials that could be used for in-situ construction on the lunar surface in the future, to mitigate the calculated effects of a large vehicle rocket engine landing and launching on the Moon.

Lunar↗
DOE OSTI2026

A Sensitivity-driven Wide Area Protection (SWAP) Coordination Tool for High Penetration of Inverter-based Resources (IBR)

Traditionally, power system generation sources have been composed of synchronous generators, of which the fault current behavior is understood with minimal differences between generation size and types due to the physics of their construction. Present protection schemes and modeling methods are based upon these understood characteristics. Most renewable generation is composed of inverter-based resources (IBR), in which fault current is determined by switching control software and hardware limitations, each of which can vary between manufacturers and even between models of the same manufacturer. The resulting fault current is low in magnitude, low in negative-sequence current, unpredictable phase angles, and is a challenge to model. These characteristics also result in a challenge to traditional protection schemes and fault simulation software. To address several of these concerns, the project has the following goals: 1. Improve IBR models: Improve IBR models used in short circuit (SC) programs to accurately capture the response of IBRs at the bulk power system (BPS) level for fault and protection studies. 2. Develop automation tool: Develop an automation tool that allows engineers to identify protection coordination and sensitivity issues by performing SC and protection coordination studies in a high IBR-penetrated grid by applying variations to the IBR models, faults, contingencies, etc. 3. Develop schemes: Develop new protection mitigation solution schemes that complement the existing protection systems to ensure safe operation of the BPS with higher IBR penetration levels. The project team did not achieve this final goal, as the Department of Energy (DOE) stopped the project early due to changes in DOE funding priorities. The termination notice came at the beginning of the final project phase, while the team was identifying and beginning to investigate protection issues. It should be noted that the team discussed a 100% penetration scenario. However, this scenario would require the use of grid-forming IBR models that are not presently available. Since developing these models requires additional effort, the 100% penetration scenario was not pursued during this project. In the future, developing the methodology and models for the 100% scenario could benefit the industry.

14 SOLAR ENERGY↗
DOE OSTI2023

Analysis of thermal radiation of a gas-particulate-cloud in a pressurized oxy-coal combustor

Pressurized oxy-combustion (POC) aims to reduce carbon emissions from coal-fired power plants. In the POC process, pulverized coal is burned at elevated pressure in an O2-CO2 environment. The elevated pressure and high CO2 and H2O concentrations strongly impact thermal radiation, thereby significantly distinguishing POC system from conventional atmospheric-pressure air-fired and oxy-fuel combustion. Consequently, a thorough understanding of thermal radiation in POC is required to promote the development of novel combustors. The present computational work analyzes thermal radiation of the gaseous and particle phases in the pilot-scale POC combustor in Washington University in St. Louis. Particularly, the radiation characteristics of the gaseous and particulate cloud were estimated by employing the statistical narrow band model and the large-particle model. It is found that thermal radiation of a particulate cloud dominates in the combustor under the conditions of a furnace temperature of 1500 K and no substantial particle loss. Another important result is that the gas-and-particulate cloud can be approximately treated as a graybody under the same conditions. To be more specific, the results of spectral radiation intensity of a gas comprised of 40% (vol) H2O and 60% CO2 show that when the pressure is increased to 15 bar, and the radiation pathlength to 100 cm, the spectral radiation profile of the gas phase approaches that of a blackbody at the respective temperature. In addition, the emissivity of the particulate cloud has been evaluated as a function of the particle concentration and diameter by means of the large particle model. It is shown that the emissivity grows with the particle concentration but decreases with the particle size for the same mass of a particulate cloud. The emissivity of a particulate cloud in oxy- combustion exceeds that in air-fired combustion. This work could be used to validate the assumption of the gray gas model adopted in numerical simulations.

Li, Lei↗