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,027 records · Page 57

Development & Validation of Low-Cost, Highly-Durable, Spinel-Based Materials for SOFC Cathode-Side Contact (Final Report)

A cathode-side contact layer is required to provide and maintain stable electrical conduction paths between the interconnect and cathode in a solid oxide fuel cell (SOFC) stack assembly and thus minimize the ohmic resistance and stack power loss. Current cathode-interconnect contact materials are based on noble metals, electrically-conductive perovskites, their composite materials, etc. These materials are either too expensive or do not possess the overall balanced performance required for the cathode-side contact application. To achieve the DOE SOFC system cost and performance stability goals, a new generation of low-cost, high-performance contact materials needs to be developed. In this project, spinel-based materials thermally converted from the Fe-Ni and Co-Mn based alloy precursors were developed and validated for the cathode-side contact application. The precursor alloy compositions were optimized via a combination of composition screening in the (Ni,Fe) 3 O 4 and (Mn,Co) 3 O 4 spinel system, alloy design using physical metallurgy principles, and cost considerations. The alloy powders with the desired composition and particle size were manufactured via gas atomization. The optimal process parameters for thermal conversion of these alloy precursor layers to a spinel-based layer were identified, i.e., 900°C x 2h in air, which is close to the initial stack firing condition. The area-specific resistances (ASRs) of the interconnect/contact/cathode test assemblies with the developed contact layer were determined for various durations (up to 5000 h) under simulated cathodic operation conditions. Some of the alloy-derived spinel contacts exhibited the lowest ASR and ASR degradation rate. The in-stack performance of the most promising alloy-derived contact layer is currently being evaluated via stack testing. To reduce the stack cost, the Co-Mn based alloy powders were utilized as the precursor for synthesis of dense spinel-based interconnect coating. By optimizing both the initial powder size/distribution and the alloy powder composition, a dense (Mn,Co) 3 O 4 -based spinel coating was achieved. Furthermore, co-sintering of the coating/contact dual-layer structure under the initial stack firing condition was realized by utilizing the tailored Co-Mn alloy precursors. Cost analysis of the developed technology indicated a total stack cost reduction of around 10.6% with the implementation of co-sintering of the interconnect coating and the contact layer during initial stack firing. Since low-cost processes such as screen printing is utilized in the precursor application and no reduction heat treatment is needed for the coating formation, the developed technology can be readily implemented at the industrial partner’s manufacturing facilities with no additional capital investment needed.

08 HYDROGEN↗

How Much Attenuation Extinguishes mm-Wave Vertically Pointing Radar Return Signals?

Vertically pointing radars (VPRs) operating at millimeter wavelengths measure the power return from raindrops enabling precipitation retrievals as a function of height. However, as the rain rate increases, there are combinations of rain rate and rain path length that produce sufficient attenuation to prevent the radar from detecting raindrops all the way through rain shafts. This study explores the question: Which rain rate and path length combinations completely extinguish radar return signals for VPRs operating between 3 and 200 GHz? An important step in these simulations is converting attenuated radar reflectivity factor into radar received signal-to-noise ratio (SNR) in order to determine the range where the SNR drops below the receiver detection threshold. Configuring the simulations to mimic a U.S. Department of Energy Atmospheric Radiation Mission (ARM) W-band (95 GHz) radar deployed in Brazil, the simulation results indicate that a W-band radar could observe raindrops above 3.5 km only when the rain rate was less than approximately 4 mm h -1 . The deployed W-band radar measurements confirm the simulation results with maximum observed heights ranging between 3 and 4.5 km when a surface disdrometer measured 4 mm h -1 rain rate (based on 25-to-75 percentiles from over 25,000 W-band radar profiles). In summary, this study contributes to our understanding of how rain and atmospheric gas attenuation impacts the performance of millimeter-wave VPRs and will help with the design and configuration of multi-frequency VPRs deployed in future field campaigns.

54 ENVIRONMENTAL SCIENCES↗

Ambulation Increases Decompression Sickness in Altitude Exposure

INTRODUCTION - Exercise accelerates inert gas elimination during oxygen breathing prior to decompression (prebreathe), but may also promote bubble formation and increase the risk of decompression sickness (DCS). The timing, pattern and intensity of exercise are likely critical to the net effect. The NASA Prebreathe Reduction Program (PRP) combined oxygen prebreathe and exercise preceding a 4.3 psi exposure in non-ambulatory subjects (a microgravity analog) to produce two protocols now used by astronauts preparing for extravehicular activity (CEVIS and ISLE). Additional work is required to investigate whether exercise normal to 1 G environments increases the risk of DCS over microgravity simulation. METHODS - The CEVIS protocol was replicated with one exception. Our subjects completed controlled ambulation (walking in place with fixed cadence and step height) during both preflight and at 4.3 psi instead of remaining non-ambulatory throughout. Decompression stress was graded with aural Doppler (Spencer 0-IV scale). Two-dimensional echocardiographic imaging was used to look for left heart gas emboli (the presence of which prompted test termination). Venous blood was collected at three points to correlate Doppler measures of decompression stress with microparticle (cell fragment) accumulation. Fisher Exact Tests compared test and control groups. Trial suspension would occur when DCS risk >15% or grade IV venous gas emboli (VGE) risk >20% (at 70% confidence). RESULTS - Eleven person-trials were completed (9 male, 2 female) when DCS prompted suspension. DCS was greater than in CEVIS trials (3/11 [27%] vs. 0/45 [0%], respectively, p=0.03). Statistical significance was not reached for peak grade IV VGE (2/11 [18%] vs. 3/45 [7%], p=0.149) or cumulative grade IV VGE observations per subject-trial (8/128 [6%] vs. 26/630 [4%], p=0.151). Microparticle data were collected for 5/11 trials (3 with DCS outcomes), with widely varying patterns that could not be resolved statistically. CONCLUSION - We did find that that ambulation increases decompression stress. Additional trials would improve the statistical power to assess differences in VGE and to evaluate the relationship between decompression stress and microparticles.

Conkin, Johnny↗

Phase II Field Demonstration at Lansing Smith Generating Plant, Southport, Florida (Final Report)

The Final Technical Report: Field Test Design and Pressure Management Strategies for Phase II Field Demonstration of Optimal Pressure Control, Plume Management, and Produced Water Strategies presents the culmination of multi-year efforts under the U.S. Department of Energy’s Brine Extraction and Storage Test (BEST) program, led by EPRI in partnership with Gulf Power at the Plant Smith site near Panama City, Florida. The project was designed to evaluate and demonstrate the technical feasibility of managing subsurface pressures and fluid movement associated with large-scale CO₂ injection, using low-salinity water as a proxy fluid. Through a combination of field injection testing, reservoir modeling, and optimization studies, the research team developed and refined active and passive brine extraction strategies aimed at controlling injection-induced pressure buildup, mitigating risks of fault activation, and managing plume migration. The field demonstration incorporated a new injection well (TIW-2), a new monitoring/extraction well (TEMW-A), and an existing well (TIW-1) repurposed for passive pressure relief. Complementary geophysical monitoring was designed to track plume development and evaluate the effectiveness of pressure management strategies. The report details the integrated workflow encompassing wellfield development, regulatory permitting, model calibration, and survey design. It includes a comprehensive electromagnetic (EM) modeling and inversion study used to develop a cost-effective, time-lapse geophysical monitoring plan capable of imaging the evolving low-salinity plume within the highly saline Lower Tuscaloosa formation. Reservoir simulation results guided the design of the 17-month injection program and the timing of active extraction to maintain formation pressures below a hypothetical fault reactivation threshold. Supporting analyses evaluated potential injectivity risks related to clay fines migration and geochemical incompatibility, identifying practical mitigation measures such as salinity and pH control. The final design recommends a series of crosswell EM surveys—conducted before, during, and after injection—as the most effective and economical approach for plume imaging, supplemented by continuous downhole pressure and flow monitoring. Collectively, the study provides a field-ready framework for cost-effective pressure management and monitoring in support of future CO₂ storage projects.

01 COAL, LIGNITE, AND PEAT↗

Passive Energy-Saving Solutions for Clothes Dryers: A Modeling and Experimental Study for Improved Efficiency and Affordability

Dryers are integral appliances in modern households, yet their significant power consumption remains a critical challenge for reducing energy bills and upgrades. This work presents the outcomes of a comprehensive investigation aimed at lowering the operating energy costs of the dryer through modeling and experimental approaches. Specifically, the objectives of this research are twofold: (1) to reduce energy consumption without negatively impacting drying performance or time; and (2) to ensure affordability by developing retrofittable solutions characterized by low cost and a quick payback period. These strategies promise universal applicability to all dryer categories, encompassing both gas and electric models, by tackling core issues such as unnecessary heat loss and excess energy supply, both prevalent across dryer designs. Our methodology combines robust modeling frameworks and experimental validation to target energy efficiency improvements through three primary pathways: (1) effective heat loss management using ultralow-cost insulation materials tailored for dryer systems; (2) heat management across the drying cycle, facilitated by passive heat transfer mechanisms, and (3) optimization of heat supply to the load to ensure precise energy delivery, minimizing waste. These innovations are designed to seamlessly integrate with existing dryer configurations, providing a scalable and retrofittable solution that ensures affordability without requiring substantial redesigns or expensive components. Using these methods, the study demonstrates improvement opportunities in dryer energy consumption while maintaining the desired performance. The modeling component utilizes computational simulations to evaluate the thermal and energy performance of these innovations under varied operational conditions, providing foundational insights for experimental design. By addressing critical areas such as heat loss, energy recovery, and supply optimization, this work proposes impactful solutions for lowering the energy costs in domestic clothes dryers. The ultralow-cost, retrofittable nature of the proposed strategies ensures widespread adoption potential across diverse dryer categories, making significant strides toward affordable household energy practices.

Cheekatamarla, Praveen [ORNL] (ORCID:0000000248827↗

Testing the pressure-confined Ly alpha cloud model

The Ly-alpha absorption line forest, seen in quasar spectra, is generally interpreted as being due to cosmologically distributed 'clouds' of primordial gas. Analyses of the observations reveal that the number distribution can be described by power laws: dN/dz is proportional to (1 + z)(sup gamma) and dN/dN(sub HI) is proportional to N(sub HI)(sup -beyda), where N(sub HI) is the HI column density. The typical values for power law indices range between 2 is approximately less than gamma is approximately less than 2.6 and 1.7 is approximately less than gamma is approximately less than 1.9. One model postulates that the Ly-alpha clouds are optically thin entities, photoionized by the background UV flux, J(sub nu) is proportional to (1 + z)(sup j), and confined by an adiabatically evolving intercloud medium (ICM): P(z) is proportional to (1 + z)(sup 5). Analytic studies of this model suggest that the ensuing Ly-alpha line statistics can account for the observations (in particular, the dN/dz and the dN/dN(sub HI) distributions) if the cloud mass spectrum is a power law dN/dN is proportional to M(sup -delta), delta is approximately 1.9, and j is approximately 4. One of the simplifying assumptions incorporated into these studies is the existence of a large mass range for the clouds at all epochs, the validity of which is questionable. The pressure-confined model is investigated using a 1-D spherically symmetric hydrodynamical code to simulate cloud evolution over the epoch 1.8 less than z less than 6. This enables us to relax many of the assumptions incorporated in the analytic studies.

Babul, Arif↗

Flow Characteristics of UHV Nitrogen Arcjet

Flow characteristics in the nozzle of a nitrogen arcjet have been simulated by a model. Electrons, ions, and neutral atoms and molecules are represented as separate species. Dissociation, ionization, and recombination are treated as separate reactions. Thermal non-equilibrium is represented by a two-temperature model. Energy input to plasma from the arc is modeled as a source determined by local current and electrical conductivity. Momentum and energy losses and recombinations at the nozzle wall are included as source/sink determined by fluxes to the wall. The results show essentially frozen gas-phase reactions and thermal nonequilibrium due to Joule heating of electrons and in the expansion part of the nozzle. Wall interactions have strong effects on the results, indicating that they play important roles in the flow due to very high surface area relative to the volume of the plasma. For example, plasma velocity decreases from 6000 m per s without wall interactions to 3000 m per s with modest amount of wall interactions. These results also show reasonable agreement with optical emission measurements, which confirms that the arcjet plasma is far from LTE. The spectra suggest nitrogen dissociation levels of 0.3% - 9%, depending on nitrogen flow rate and arc plasma current, which also determine the relative amounts of excited atomic and molecular nitrogen. Langmuir probe studies of the source show that electron and ion fluxes increase with increasing power, and that the ion energy distribution shifts to lower energies. Typical ion fluxes were on the order of 4E-9 A/cm**2 with a maximum ion kinetic energy of 3.5eV. The median electron energy was 1eV, with a maximum of less than 4eV.

Chang, C. H.↗

J-2X Abort System Development

The J-2X is an expendable liquid hydrogen (LH2)/liquid oxygen (LOX) gas generator cycle rocket engine that is currently being designed as the primary upper stage propulsion element for the new NASA Ares vehicle family. The J-2X engine will contain abort logic that functions as an integral component of the Ares vehicle abort system. This system is responsible for detecting and responding to conditions indicative of impending Loss of Mission (LOM), Loss of Vehicle (LOV), and/or catastrophic Loss of Crew (LOC) failure events. As an earth orbit ascent phase engine, the J-2X is a high power density propulsion element with non-negligible risk of fast propagation rate failures that can quickly lead to LOM, LOV, and/or LOC events. Aggressive reliability requirements for manned Ares missions and the risk of fast propagating J-2X failures dictate the need for on-engine abort condition monitoring and autonomous response capability as well as traditional abort agents such as the vehicle computer, flight crew, and ground control not located on the engine. This paper describes the baseline J-2X abort subsystem concept of operations, as well as the development process for this subsystem. A strategy that leverages heritage system experience and responds to an evolving engine design as well as J-2X specific test data to support abort system development is described. The utilization of performance and failure simulation models to support abort system sensor selection, failure detectability and discrimination studies, decision threshold definition, and abort system performance verification and validation is outlined. The basis for abort false positive and false negative performance constraints is described. Development challenges associated with information shortfalls in the design cycle, abort condition coverage and response assessment, engine-vehicle interface definition, and abort system performance verification and validation are also discussed.

Santi, Louis M.↗

Striving to translate shale physics across ten orders of magnitude: What have we learned?

Shales will play an important role in the successful transition of energy from fossil-based resources to renewables in the coming decades. Aside from being a significant source of low-carbon intensity fuels, like natural gas, they also serve as geologic seals of subsurface formations that may be used to isolate nuclear waste, sequester CO 2 , or store intermittent energy (e.g., solar hydrogen). Despite their importance, shales pose significant engineering and environmental challenges due to their nanoporous structure and extreme heterogeneity that spans at least ~10 orders of magnitude in spatial scale. Two challenges inhibit a system-level understanding: (1) the physics of fluid flow and phase behavior in shales are poorly understood due to the dominant molecular interactions between minerals and fluids under confinement, and (2) the apparent lack of scale separation that prevents a reliable (closed) description of the physics at any single scale of observation. In this review, we focus on the latter issue and discuss scale translation, which in its broadest sense is transforming data or simulations from one spatiotemporal scale to another. While effective scale translation is not exclusive to shales, but all geologic porous media, the need for it is especially acute in shales given their high degree of heterogeneity. Classical theories like homogenization, while indispensable, fail when scales are not separated. Other methods, like numerical upscaling, scale-translate in only one direction: small to large, but not the reverse, called downscaling. However, the confluence of advances in three areas are bringing challenging problems such as shales within reach: increased computational power and scalable algorithms; high-resolution imaging and multi-modal data acquisition; and machine learning to process massive amounts of data. While these advances equip geoscientists with a wide array of experimental and computational tools, no individual tool can probe the entire gamut of heterogeneity in shales. Their effective use, therefore, requires an ability to bridge between various data types obtained at different scales. The aim of this review is to present a coherent account of computational and experimental methods that may be used to achieve just that, i.e., to perform scale translation. We provide a broader definition of scale translation, one that transcends classical homogenization and upscaling methods, but is consistent with them and accommodates notions like downscaling and data translation. After a brief introduction to homogenization, we review hybrid methods, numerical upscaling and its recent extensions, multiscale computing, high-resolution imaging, and machine learning. We place particular emphasis on multiscale computing and propose an algorithmic framework to bridge between the pore (micro) and Darcy (macro) scales. Throughout the paper, we draw comparisons between the various methods and highlight their (often hidden) similarities, differences, benefits, and pitfalls. We finally conclude with two case studies on shales that exemplify some of the methods presented.

58 GEOSCIENCES↗

NASA Tech Briefs, October 2012

Topics discussed include: Detection of Chemical Precursors of Explosives; Detecting Methane From Leaking Pipelines and as Greenhouse Gas in the Atmosphere; Onboard Sensor Data Qualification in Human-Rated Launch Vehicles; Rugged, Portable, Real-Time Optical Gaseous Analyzer for Hydrogen Fluoride; A Probabilistic Mass Estimation Algorithm for a Novel 7-Channel Capacitive Sample Verification Sensor; Low-Power Architecture for an Optical Life Gas Analyzer; Online Cable Tester and Rerouter; A Three-Frequency Feed for Millimeter-Wave Radiometry; Capacitance Probe Resonator for Multichannel Electrometer; Inverted Three-Junction Tandem Thermophotovoltaic Modules; Fabrication of Single Crystal MgO Capsules; Inflatable Hangar for Assembly of Large Structures in Space; Mars Aqueous Processing System; Hybrid Filter Membrane; Design for the Structure and the Mechanics of Moballs; Pressure Dome for High-Pressure Electrolyzer; Cascading Tesla Oscillating Flow Diode for Stirling Engine Gas Bearings; Compact, Low-Force, Low-Noise Linear Actuator; Ultra-Compact Motor Controller; Extreme Ionizing-Radiation-Resistant Bacterium; Wideband Single-Crystal Transducer for Bone Characterization; Fluorescence-Activated Cell Sorting of Live Versus Dead Bacterial Cells and Spores; Nonhazardous Urine Pretreatment Method; Laser-Ranging Transponders for Science Investigations of the Moon and Mars; Ka-Band Waveguide Three-Way Serial Combiner for MMIC Amplifiers; Structural Health Monitoring with Fiber Bragg Grating and Piezo Arrays; Low-Gain Circularly Polarized Antenna with Torus-Shaped Pattern; Stereo and IMU- Assisted Visual Odometry for Small Robots; Global Swath and Gridded Data Tiling; GOES-R: Satellite Insight; Aquarius iPhone Application; Monitoring of International Space Station Telemetry Using Shewhart Control Charts; Theory of a Traveling Wave Feed for a Planar Slot Array Antenna; Time Manager Software for a Flight Processor; Simulation of Oxygen Disintegration and Mixing With Hydrogen or Helium at Supercritical Pressure; A Superfluid Pulse Tube Refrigerator Without Moving Parts for Sub-Kelvin Cooling; Sapphire Viewports for a Venus Probe; The Mobile Chamber; Electric Propulsion Induced Secondary Mass Spectroscopy; and Radiation-Tolerant DC-DC Converters.

Source record↗

Online randomized interpolative decomposition with a posteriori error estimator for temporal PDE data reduction

Traditional low-rank approximation is a powerful tool for compressing large data matrices that arise in simulations of partial differential equations (PDEs), but suffers from high computational cost and requires several passes over the PDE data. The compressed data may also lack interpretability thus making it difficult to identify feature patterns from the original data. Here, to address these issues, we present an online randomized algorithm to compute the interpolative decomposition (ID) of large-scale data matrices in situ. Compared to previous randomized IDs that used the QR decomposition to determine the column basis, we adopt a streaming ridge leverage score-based column subset selection algorithm that dynamically selects proper basis columns from the data and thus avoids an extra pass over the data to compute the coefficient matrix of the ID. In particular, we adopt a single-pass error estimator based on the non-adaptive Hutch++ algorithm to provide real-time error approximation for determining the best coefficients. As a result, our approach only needs a single pass over the original data and thus is suitable for large and high-dimensional matrices stored outside of core memory or generated in PDE simulations. A strategy to improve the accuracy of the reconstructed data gradient, when desired, within the ID framework is also presented. We provide numerical experiments on turbulent channel flow and ignition simulations, and on the NSTX Gas Puff Image dataset, comparing our algorithm with the offline ID algorithm to demonstrate its utility in real-world applications.

Column subset selection↗

Improving Efficiency of Off-Road Vehicles by Novel Integration of Electric Machines and Advanced Combustion Engines

Modern off-road equipment will increasingly rely on electrified implements that will deliver precision control with a smaller footprint than hydraulics. The primary energy converter, however, will be an onboard reciprocating internal combustion engine because of the power density of hydrocarbon fuels in comparison to electrical energy storage and the remote locations where much of this equipment is deployed. Adding energy storage and electric machines creates opportunities to improve efficiency while reducing emissions. This program investigated approaches to take advantage of the extra flexibility that an enhanced electrical system to enable high-efficiency, low-emissions combustion technologies. The current program evaluated hybridization of both the torque application and air-handling systems to maximize efficiency while minimizing cost. This program designed, analyzed and tested a hybrid off-road vehicle consisting of a series electric powertrain with energy storage, an electrified air system, and a 33% downsized diesel engine. Detailed comparisons were made between the base powertrain and the developed powertrain using powertrain simulations, engine testing, and vehicle testing. The key results of the study are: • The resulting vehicle reduces fuel consumption by 5% to 15% at equal productivity. The primary improvements were due to recovery of regenerative braking losses for cycles where large transients are encountered. • The electrified air system was an enabler to allow engine downsizing. This was most important for duty-cycles where the engine primarily operated at moderate loads and had few low speed high torque conditions. • Engine level improvements showed re-optimization of the powertrain is possible when electrified air handing is available. Increased exhaust gas recirculation and advanced injection timing allowed up to a 15% reduction in brake specific fuel consumption at equal NOx and transient response. Compared to the larger engine, the downsized engine achieves up to an 18% reduction in brake specific fuel consumption over the non-road transient cycle. • The life cycle analysis and total cost of ownership study showed that the hybrid powertrain has the potential to reduce 5-year CO2 and total cost of ownership by ~6% over the baseline vehicle. The results also indicated that a pure battery electric vehicle is not feasible in this application and is likely to increase the CO2 emissions due to the CO2 from battery production. A hybrid powertrain with low carbon fuels shows the potential to substantially reduce CO2 and total cost of ownership.

02 PETROLEUM↗

Divertor Plasma Detachment Control Neural Network

DivControlNN is a state-of-the-art software tool that leverages advanced machine learning techniques to predict and control divertor plasma behavior in fusion reactors. Plasma, a highly energetic and electrically charged gas, requires meticulous management to protect reactor components and maintain optimal energy production. Conventional simulation methods, although extremely detailed, typically demand extensive computational time-making them unsuitable for real-time control scenarios. DivControlNN addresses this challenge by learning from tens of thousands of high-fidelity simulations, thereby creating a rapid surrogate model that can deliver near-instantaneous predictions. At the core of its functionality is a sophisticated technique known as latent space mapping, which condenses complex, high-dimensional plasma data into a compact, lower-dimensional representation. This streamlined representation enables the system to quickly forecast essential plasma properties and determine the precise conditions required for effective detachment. Detachment is a crucial process in which the plasma is cooled before reaching the divertor plates, thereby reducing heat loads and mitigating material erosion. In recent experiments conducted on the KSTAR tokamak in South Korea, DivControlNN successfully guided the detachment process without any fine-tuning-even when applied to a new tungsten divertor configuration. By achieving a computational speed-up of over one hundred million times compared to traditional simulation methods while maintaining low prediction errors, DivControlNN stands to significantly enhance real-time control and diagnostic capabilities in future fusion reactors. This breakthrough paves the way for safer, more reliable reactor operation and represents a major advancement toward realizing fusion energy as a practical, sustainable, and clean power source.

Xu, Xueqiao [Lawrence Livermore National Laborator↗

Investigating instabilities in magnetized low-pressure capacitively coupled RF plasma using particle-in-cell (PIC) simulations

The effect of a uniform magnetic field on particle transport in low-pressure radio frequency (RF) capacitively coupled plasma (CCP) has been studied using a particle-in-cell model. Three distinct regimes of plasma behavior can be identified as a function of the magnetic field. In the first regime at low magnetic fields, asymmetric plasma profiles are observed within the CCP chamber due to the effect of $\overrightarrow{E}$ × $\overrightarrow{B}$ drift. As the magnetic field increases, instabilities develop and form self-organized spoke-shaped structures that are distinctly seen within the bulk plasma closer to the sheath. In this second regime, the spoke-shaped coherent structures rotate inside the plasma chamber in the -$\overrightarrow{E}$ × $\overrightarrow{B}$ direction, where $\overrightarrow{E}$ × $\overrightarrow{B}$ are the DC electric and magnetic field vectors, respectively, and the DC electric field exists in the sheath and pre-sheath regions. The spoke rotation frequency is in the megahertz range. As the magnetic field strength increases further, the rotating coherent spokes continue to exist near the sheath. The coherent structures are, however, accompanied by new small-scale incoherent structures originating and moving within the bulk plasma region away from the sheath. This is the third regime of plasma behavior. The threshold values of the magnetic field between these regimes were found not to vary with changing plasma reactor geometry (e.g., area ratio between ground and powered electrodes) or the use of an external capacitor between the RF-powered electrode and the RF source. The threshold values of the magnetic field between these regimes shift toward higher values with increasing gas pressure. Analysis of the results indicates that the rotating structures are due to the lower hybrid instability driven by density gradients and electron-neutral collisions. This paper provides guidance on the upper limit of the magnetic field for instability-free operation in low-pressure CCP-based semiconductor deposition and etch systems that use the external magnetic field for plasma uniformity control.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A Simple Condition for Sustained Super-Eddington Black Hole Growth

One of the most pressing questions in cosmology is how the black holes (BHs) powering quasars at high redshift grow to supermassive scales within a billion years of the Big Bang. Here we show that sustained super-Eddington accretion can be achieved for BHs with Eddington fractions f Edd ≳ 2/ϵ, where ϵ is the efficiency with which radiation is generated in the accretion process. In this regime, the radiation carries too little momentum to halt the accretion flow and the infalling gas traps the radiation. The BH growth then proceeds unimpeded until the gas supply is exhausted, in contrast to accretion at lower rates, which is limited by the radiation generated in the accretion process. The large gas supply available in massive high-redshift quasar host galaxies may be readily accreted onto seed BHs via this supply-limited mode of accretion, providing an explanation for how such supermassive BHs are assembled in the early universe. This sustained super-Eddington growth may also explain the short lifetimes inferred for the H ii regions surrounding high-redshift quasars, if the bulk of the BH growth occurs without the associated radiation escaping to ionize the intergalactic medium. It furthermore implies that a population of obscured rapidly growing BHs may be difficult to detect, perhaps explaining why so few quasars with Eddington fractions higher than a few have been observed. Finally, this simple condition for sustained super-Eddington growth can easily be implemented in cosmological simulations that can be used to assess in which environments it occurs.

79 ASTRONOMY AND ASTROPHYSICS↗

Plasma decay in hydrocarbons and hydrocarbon- and H 2 O-containing mixtures excited by high-voltage nanosecond discharge at elevated gas temperatures

Plasma decay after a high-voltage nanosecond discharge was experimentally and numerically studied in pure hydrocarbons (C 2 H 6 and C 3 H 8 ), and H 2 O:N 2 and C 3 H 8 :O 2 mixtures for pressures in the range 2–4 Torr and gas temperatures from 300 to 600 K. In a stoichiometric C 3 H 8 :O 2 mixture, plasma decay was also studied in a repetitively pulsed discharge for varing numbers of discharge pulses (varying degrees of fuel oxidation). The rate of plasma decay was determined from the temporal evolution of electron density measured using the microwave interferometer. It was observed that gas heating to 600 K leads to a decrease in the rate of plasma decay in all cases. The effect of heating on plasma decay was most profound in the H 2 O:N 2 mixture (after a single discharge pulse) and in the C 3 H 8 :O 2 mixture for high degrees of fuel oxidation. A kinetic scheme was developed to numerically simulate the plasma decay in hydrocarbons and combustible mixtures. Numerical analysis showed that, under the conditions studied, plasma decay was controlled by dissociative electron recombination with simple molecular and cluster ions. Gas heating led to a decrease in the rate of the electron-ion recombination and the rate of conversion of molecular ions to cluster ions. As a result, the gas temperature increase caused a decrease in the fraction of cluster ions for which the recombination coefficients are an order of magnitude higher than the recombination coefficients for molecular ions. The influence of gas heating on the decrease of the amount of cluster ions was more important when the ion composition was dominated by hydrated H 3 O + (H 2 O) k ions. The rates of the formation of these ions are extremely sensitive to any variations in gas temperature. A result, in agreement with our observations, gas heating led to an anomalous decrease in the rates of plasma decay in the H 2 O:N 2 mixture, as well as in the C 3 H 8 :O 2 mixture when H 2 O molecules were produced due to fuel oxidation.

20 FOSSIL-FUELED POWER PLANTS↗

Accuracy Enhancement of Nuclear Power Plant Simulators Utilizing High Accuracy Simulation Predictions

More recently, reactor core simulators for core designs associated with commercial nuclear power plants that utilize what is believed to be higher fidelity models have been developed. Features such as neutronics models that utilize transport equation solvers with fine spatial meshes and many energy-groups, thermal-hydraulic models that utilize sub-channel solvers with fine spatial mesh and capable of treating a wide range of fluid conditions, and fuel-coolant chemistry interaction models capable of treating CRUD deposition are to be found in these higher fidelity core simulators. These reactor core simulators require access to higher performance computers, characterized by many processors, cores and large memory. So associated with utilization of these simulators is access to high performance computers and ability to accommodate in one’s workflow longer execution times. By contrast, currently used core simulators by the nuclear industry can execute on engineering workstations and have execution times of seconds to minutes. The desirability for having short execution times is not only desired for support of time critical tasks but supports the mental process of decision making by engineers. The goal of the work reported upon here has the objective of retaining the fidelity of higher fidelity models while retaining the ability to utilize engineering workstations. Beyond the core simulator goal, additional goals of this work include incorporating the just described core simulator capability into a Nuclear Steam Supply System (NSSS) simulator, and to incorporate the resulting capability into an environment supportive of design and operational decision making associated with nuclear power stations. The model selected for the core neutronics model is the NESTLE code, for the core thermal-hydraulic model is the CTF code utilizing coarse mesh, and for the NSSS model is the RELAP5-3D code. WSC’s proprietary 3KEYMASTERTM platform is being used to provide software coupling, user interface, visualization, and reporting. The NESTLE core neutronics simulator was first integrated with the CTF core thermal-hydraulic simulator using CTF developed communication commands which are also used for CTF to communicate with RELAP5-3D under WSC’s proprietary 3KEYMASTERTM platform. To assure NESTLE prediction consistency with higher fidelity core neutronic simulators, buffer codes have been created to automatically generate from output files written by the VERA core simulator the NESTLE nodal neutronic parameter’ library, geometry, and pin-power reconstruction input files, thereby avoiding a number of challenges associated with utilizing lattice physics codes and providing consistency with VERA predictions. To treat absorber rod effects a multi-set library is utilized, where a set refers to a specific absorber rod fully inserted pattern. A coarse spatial mesh CTF model was developed with features added that support using CTF as envisioned in the engineering quality simulator. A hybrid meshing approach was implemented to allow for automated construction of models with mixed levels of refinement. Specifically, a core model could resolve some assemblies at a nodal level (4 subchannels per assembly) and others at a pin-resolution (one subchannel per coolant subchannel in the assembly). The intention is that this will allow for better resolution of limiting conditions such as DNBR and PCT, which are based on local rod and subchannel conditions. Further development was done of features that enhance the capabilities for the envisioned engineering quality simulator that has been developed, but now for RELAP-3D. The RELAP5-3D code development includes ability to model more than 999 components and the addition of the cross-channels turbulence mixing model and the void drift model that are implemented in CTF, aiming to achieve closer prediction agreement of the two codes for transient simulations, specifically, more accurate matches of the overall mass, momentum, and energy exchanges of both the liquid and gas phases between the neighboring core assemblies. Graphics were also developed for the Instructor Station for this project under WSC’s proprietary 3KEYMASTERTM platform to facilitate design and operational decision making.

42 ENGINEERING↗

Magnetohydrodynamic MACH Code Used to Simulate Magnetoplasmadynamic Thrusters

The On-Board Propulsion program at the NASA Glenn Research Center is utilizing a state of-the-art numerical simulation to model the performance of high-power electromagnetic plasma thrusters. Such thrusters are envisioned for use in lunar and Mars cargo transport, piloted interplanetary expeditions, and deep-space robotic exploration of the solar system. The experimental portion of this program is described in reference 1. This article describes the numerical modeling program used to guide the experimental research. The synergistic use of numerical simulations and experimental research has spurred the rapid advancement of high-power thruster technologies for a variety of bold new NASA missions. From its inception as a U.S. Department of Defense code in the mid-1980's, the Multiblock Arbitrary Coordinate Hydromagnetic (MACH) simulation tool has been used by the plasma physics community to model a diverse range of plasma problems--including plasma opening switches, inertial confinement fusion concepts, compact toroid formation and acceleration, z-pinch implosion physics, laser-target interactions, and a variety of plasma thrusters. The MACH2 code used at Glenn is a time-dependent, two-dimensional, axisymmetric, multimaterial code with a multiblock structure. MACH3, a more recent three-dimensional version of the code, is currently undergoing beta tests. The MACH computational mesh moves in an arbitrary Lagrangian-Eulerian (ALE) fashion that allows the simulation of diffusive-dominated and dispersive-dominated problems, and the mesh can be refined via a variety of adaptive schemes to capture regions of varying characteristic scale. The mass continuity and momentum equations model a compressible viscous fluid, and three energy equations are used to simulate nonthermal equilibrium between electrons, ions, and the radiation field. Magnetic fields are modeled by an induction equation that includes resistive diffusion, the Hall effect, and a thermal source for magnetic fields. Various models of plasma resistivity are included, along with ablation models and multiport circuit solvers. The set of equations is closed using either an ideal gas or real equation of state.

Mikellides, Pavlos G.↗