Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Steady state process”

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,117 records · Page 62

High-resolution radiation detection using Ni/SiO 2 /n-4H-SiC vertical metal-oxide-semiconductor capacitor

In this article, we demonstrate the radiation detection performance of vertical metal-oxide-semiconductor (MOS) capacitors fabricated on 20 μm thick n-4H-SiC epitaxial layers with the highest energy resolution ever reported. The 100 nm SiO2 layer was achieved on the Si face of n-4H-SiC epilayers using dry oxidation in air. The Ni/SiO2/n-4H-SiC MOS detectors not only demonstrated an excellent energy resolution of 0.42% (ΔE/E 100) for 5.48 MeV alpha particles but also caused a lower enhancement in the electronic noise components of the spectrometer compared with that observed for the best high-resolution Schottky barrier detectors. The MOS detectors also exhibited a high charge collection efficiency (CCE) of 96% at the optimized operating bias despite the presence of the oxide layer. A drift-diffusion model applied to the CCE vs gate bias voltage data revealed a minority (hole) carrier diffusion length of 24 μm. Capacitance mode deep level transient spectroscopy (C-DLTS) scans in the temperature range 84–800 K were carried out to identify the resolution limiting electrically active defects. Interestingly, the C-DLTS spectra revealed both positive and negative peaks, indicating the simultaneous presence of electron (majority) and hole (minority) trap centers. It has been inferred that at the steady-state bias for the C-DLTS measurement, the MOS detector operates in the inversion mode at certain device temperatures, causing holes to populate the minority trap centers and, hence, manifests minority carrier peaks as well.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Particle-in-cell modeling of electron beam generated plasma

Plasmas generated using energetic electron beams are well known for their low electron temperature (T e ) and plasma potential, which makes them attractive for atomic-precision plasma processing applications such as atomic layer etch and deposition. A 2-dimensional particle-in-cell model for an electron beam-generated plasma in argon confined by a constant applied magnetic field is described here in this article. Plasma production primarily occurs in the path of the beam electrons in the center of the chamber. The resulting plasma spreads out in the chamber through non-ambipolar diffusion with a short-circuiting effect allowing unequal electron and ion fluxes to different regions of the bounding conductive chamber walls. The cross-field transport of the electrons (and thus the steady-state characteristics of the plasma) are strongly impacted by the magnetic field. T e is anisotropic in the electron beam region, but low and isotropic away from the plasma production zone. The plasma density increases and the plasma becomes more confined near the region of production when the magnetic field strengthens. The magnetic field reduces both electron physical and energy transport perpendicular to the magnetic field. T e is uniform along the magnetic field lines and slowly decreases perpendicular to it. Electrons are less energetic in the sheath regions where the sheath electric field repels and confines the low-energy electrons from the bulk plasma. Even though electron and ion densities are similar in the bulk plasma due to quasi-neutrality, electron and ion fluxes on the grounded chamber walls are unequal at most locations. Electron confinement by the magnetic field weakens with increasing pressure, and the plasma spread out farther from the electron beam region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Extension of the high-resolution thermal-hydraulics code ESCOT to hexagonal core geometries for multi-physics calculations

The extension of the capabilities of the pin-level nuclear reactor core thermal-hydraulics (T/H) code ESCOT to analyze hexagonal fueled cores and its performance are presented. ESCOT is an accurate yet fast core thermal-hydraulics solution aiming at high-fidelity and high-resolution multi-physics core analysis in the framework of massively parallel computing platforms. Its algorithm solution is based on the four-equation drift-flux model for two-phase calculations, these are numerically solved by applying the Finite Volume Method (FVM) and the Semi-Implicit Method for Pressure-Linked Equation (SIMPLE)-like algorithm in a staggered grid system. Constitutive models such as turbulent mixing, pressure drop, and vapor generation are employed to simulate key phenomena in subchannel-scale analysis. ESCOT is parallelized by a double (radial and axial) domain decomposition that enables its highly parallelized execution. The coupling of the code with the neutronics whole core solver for hexagonal geometries nTRACER is described. The newly implemented ESCOT features are validated by comparing single assembly and full core steady state nTRACER-ESCOT solutions with nTRACER standalone internal one-dimensional T/H solver results. The validation problems are based on the VVER 440 and VVER 1000 cores. ESCOT results show differences within an acceptable range with respect to the simple 1D nTRACER built-in solver. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Nonlinear displacement analysis of advanced propeller structures using NASTRAN

The steady state displacements of a rotating advanced turboprop are computed using the geometrically nonlinear capabilities of COSMIC NASTRAN Rigid Format 4 and MSC NASTRAN Solution 64. A description of the modified Newton-Raphson algorithm used by Solution 64 and the iterative scheme used by Rigid Format 4 is provided. A representative advanced turboprop, SR3, was used for the study. Displacements for SR3 are computed for rotational speeds up to 10,000 rpm. The results show Solution 64 to be superior for computating displacements of flexible rotating structures. This is attributed to its ability to update the displacement dependent centrifugal force during the solution process.

Lawrence, C.↗

Convection in the Physical Vapor Transport Process-I: Thermal

The effects of convection on diffusive-convective physical vapor transport process are examined computationally. We analyze conditions ranging from typical laboratory conditions to conditions achievable only in a low gravity environment. This corresponds to thermal Rayleigh numbers Ra, ranging from 1.80 x 10 to 1.92 x 10(exp 6). Our results indicate that the effect of the sublimation and condensation fluxes at the boundaries is to increase the threshold of instability. For typical ground based conditions, time dependent oscillatory convection can occur. This results in unsteady transport, and non- uniform temperature and concentration gradients at the crystal interface. Spectral analysis of the flow field shows parametric regions exhibiting both an oscillatory approach to steady state and a chaotic transient to a periodic state. Low gravity conditions stabilize the flow field. Convective effects are effectively reduced, thus resulting in uniform temperature and concentration gradients at the interface, a desirable condition for crystal growth.

Duval, Walter M. B.↗

Protoplanetary Nebula Evolution using the Beta Viscosity Model

The evolutionary dynamics of a protoplanetary disk is an important component of the planet formation process. In particular, the dynamic and thermodynamic field plays a critical role in chemical evolution, the migration of dust particles in the nebula, and the radial transport of meteoritic components. The dynamic evolution is investigated using analytical solutions of the surface density transport equations using a turbulence model based on hydrodynamic generation of turbulence. It captures the major properties of the disk including region of separation between radial inflow and-outflow and the evolution of the central plane temperature. The analytical formulas are compared with available numerical solutions based on the alpha viscosity model. The beta viscosity model, heretofore used for steady-state disks, is shown to be a useful approximation for unsteady problems.

Davis, Sanford S.↗

Flight Test of a 30-Foot Nominal-Diameter Disk-Gap-Band Parachute Deployed at Mach 1.56 and Dynamic Pressure of 11.4 Pounds per Square Foot

A 30-foot (9.1 meter) nominal-diameter disk-gap-band parachute (reference area 707 sq ft (65.7 m(exp 2)) was flight tested with a 200-pound (90.7 kg) instrumented payload as part of the NASA Planetary Entry Parachute Program. A deployment mortar ejected the test parachute when the payload was at a Mach number of 1.56 and a dynamic pressure of 11.4 lb/sq ft (546 newtons per m 2 ) at an altitude of 127,500 feet (38.86 km). The parachute reached suspension line stretch in 0.37 second resulting in a snatch force loading of 1270 pounds (5650 N). Canopy inflation began 0.10 second after line stretch. A delay in the opening process occurred and was apparently due to a momentary interference of the glass-fiber shroud used in packing the parachute bag in the mortar. Continuous canopy inflation began 0.73 second after initiation of deployment and 0.21 second later full inflation was attained for a total elapsed time from mortar fire of 0.94 second. The maximum opening load of 3915 pounds (17,400 newtons) occurred at the time the canopy was first fully opened. The parachute exhibited an average drag coefficient of 0.52 during the deceleration period and pitch-yaw oscillations of the canopy were less than 5 degrees. During the steady-state descent portion of the test period, the average effective drag coefficient was about 0.47 (based on vertical descent velocity and total system weight).

EXPERIMENTS↗

Mixing Cell Model: A One-Dimensional Numerical Model for Assessment of Water Flow and Contaminant Transport in the Unsaturated Zone

This report describes the Mixing Cell Model code, a one dimensional model for water flow and solute transport in the unsaturated zone under steady state or transient flow conditions. The model is based on the principles and assumptions underlying mixing-cell model formulations. The unsaturated zone is discretized into a series of independent mixing cells. Each cell may have unique hydrologic, lithologic, and sorptive properties. Ordinary differential equations describe the material (water and solute) balance within each cell. Water-flow equations are derived from the continuity equation, assuming that unit gradient conditions exist at all times in each cell. Pressure gradients are considered implicitly through model discretization. Unsaturated hydraulic conductivity and moisture contents are determined by the material specific moisture-characteristic curves. Solute-transport processes include explicit treatment of advective processes, first order chain decay, and linear sorption reactions. Dispersion is addressed through implicit and explicit dispersion. Implicit dispersion is an inherent feature of all mixing-cell models and originates from the formulation of the problem in terms of mass balance around fully mixed volume elements. Expressions are provided that relate implicit dispersion to the physical dispersion of the system. Two FORTRAN codes were developed to solve the water flow and solute-transport equations: (1) the Mixing Cell Model for Flow (MCMF) solves transient water-flow problems and (2) the Mixing Cell Model for Transport (MCMT) solves the solute-transport problem. The transient water-flow problem is typically solved first by estimating the water flux through each cell in the model domain as a function of time using the MCMF code. These data are stored in either ASCII or binary files that are later read by the solute transport code (MCMT). Code output includes solute pore water concentrations, water and solute inventories in each cell and at each specified output time, and water and solute fluxes through each cell and specified output time. Computer run times for coupled transient water flow and solute transport were typically several seconds on a 2 GHz Intel Pentium IV desktop computer. The model was benchmarked against analytical solutions and finite element approximations to the partial differential equations (PDE) describing unsaturated flow and transport. Differences between the maximum solute flux estimated by the mixing cell model and the PDE models were typically less than two percent. This revision includes an option for a fixed concentration lower boundary condition for diffusive fluxes for versions 020321 and later.

54 ENVIRONMENTAL SCIENCES↗

Convection in the Physical Vapor Transport Process: Thermosolutal - Part 2

We consider the effect of an inert gas on the diffusive-convective physical vapor transport process. We investigate the case when the temperature gradient is stabilizing and the concentration gradient is destabilizing for a wide parametric range. When an inert gas is present, the thermal and solutal convection oppose each other. The solutal field is destabilizing while the thermal field and the advective-diffusive flux stabilize the flow field. When the pressure of the inert component is increased, the stabilizing effect of the advective-diffusive flux is decreased. The intensity of convection as well as the oscillatory transient time increases. Below, the critical Rayleigh number, the nonlinear dynamics of the flow field show an oscillatory approach to steady state. For parametric values in the neighborhood of the critical Rayleigh number, the flow field undergoes a chaotic transient which settles to a periodic state. The asymptotic state of the flow field shows that growth and amalgamation of cells yields an overturning motion which results in an asymmetric cellular structure. The low gravity environment yields the stabilizing advective-diffusive flow which results in uniform temperature and concentration gradients near the crystal interface.

Duval, Walter M. B.↗

Axisymmetric Two-Dimensional Computation of Magnetic Field Dragging in Accretion Disks

In this paper we model a geometrically thin accretion disk interacting with an externally imposed, uniform, vertical magnetic field. The accretion flow in the disk drags and distorts field lines, amplifying the magnetic field in the process. Inside the disk the radial component of the field is sheared into a toroidal component. The aim of this work is to establish the character of the resultant magnetic field and its dependence on the disk's parameters. We concentrate on alpha-disks driven by turbulent viscosity. Axisymmetric, two-dimensional solutions are obtained without taking into account the back-reaction of the magnetic field on the structure of the disk. The character of the magnetic field depends strongly on the magnitude of the magnetic Prandtl number, P . We present two illustrative examples of viscous disks: a so-called 'standard' steady state model of a disk around a compact star (e.g., cataclysmic variable), and a steady state model of a proto-planetary disk. In both cases, P = 1, P = 10(sup -1), and P = 10(sup -2) scenarios are calculated. Significant bending and magnification of the magnetic field is possible only for disks characterized by P of the order of 10(sup -2). In such a case, the field lines are bent sufficiently to allow the development of a centrifugally driven wind. Inside the disk the field is dominated by its toroidal component. We also investigate the dragging of the magnetic field by a nonviscous protoplanetary disk described by a phenomenological model. This scenario leads to large distortion and magnification of the magnetic field.

Reyes-Ruiz, Mauricio↗

Development of k{sub eff} in the unfolding process of the LOTUS reactor design

An innovative small thermal reactor concept named LOTUS, which is a pool-type reactor aimed at heat supply. The passive concept applied in the LOTUS and it has a separable core. It means that its core can passively unfold into four parts along its central axis only with the buoyancy and gravity when some accidents happen, especially the wall rupture of the pool caused by earthquake, aircraft impact and so on. This idea can effectively enhance the ability of pool-type reactor for resisting some severe accidents. We have executed some pre-computations in order to check a non-supercritical system for the transient simulation. In this paper, the development of k{sub eff} have been calculated by both stepwise steady-state and transient modes using Serpent, as well as including some sensitivity analyses about the k{sub eff} to unfolding velocity, inserting velocity of control rod (CR) and time step. The result showed that k{sub eff} value decreased definitely as the reactor unfolded and the consistent change exposed between the stepwise steady state and transient result when the unfolding angle varied from 0 to 1 degree. However, the stepwise steady state k{sub eff} negatively declined after that angle, which demonstrated that for such a dynamic unfolding change the transient mode was more suitable, peculiarly for considering delayed neutrons and neutron moderation. For sensitivity calculation, the k{sub eff} decreased more rapidly as the unfolding velocity and inserting velocity of CR increased, and a smaller time step corresponded to a more stable change. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

A Concept of a Convection–Cloud Chamber to Study Aerosol–Cloud–Drizzle Interactions

Understanding and quantifying the full chain of processes from aerosol activation to drizzle formation, and the associated feedbacks to the aerosol chemical and physical properties, all within a turbulent cloud are some of the toughest challenges in atmospheric chemistry and physics and are keys to the cloud–precipitation puzzle. This paper describes a concept for a new type of research facility consisting of a cloud chamber plus associated instrumentation and computational models, to explore aerosol–cloud interactions and processing, cloud optical properties, entrainment–cloud interactions, and quantitative assessment of drizzle onset. The envisioned design is for a 3 m × 3 m × 9 m chamber, such that the height is sufficient to achieve long lifetimes for aerosol processing and for significant drizzle growth by collision and coalescence. A suite of computational tools for simulating microphysical properties in the chamber provides a digital twin for designing the chamber and a range of example experiments. Theory and test results from novel remote sensing systems for exploring chemical and physical interactions and evolution of aerosols, cloud droplets, and drizzle within turbulent clouds are described. Testing of technology needed for the operation of a large-volume chamber, including aerosol generation methods and novel materials for water vapor boundary conditions, is described. Simulations suggest that spatially uniform turbulence and microphysical properties can be sustained in a steady state, with reasonable aerosol and water vapor fluxes, and that substantial drizzle can be produced through collision and coalescence of cloud droplets. Remaining challenges for more detailed engineering design and a discussion of possible first-light experiments are described.

54 ENVIRONMENTAL SCIENCES↗

Steady-State ALPS for Real-Valued Problems

The two objectives of this paper are to describe a steady-state version of the Age-Layered Population Structure (ALPS) Evolutionary Algorithm (EA) and to compare it against other GAs on real-valued problems. Motivation for this work comes from our previous success in demonstrating that a generational version of ALPS greatly improves search performance on a Genetic Programming problem. In making steady-state ALPS some modifications were made to the method for calculating age and the method for moving individuals up layers. To demonstrate that ALPS works well on real-valued problems we compare it against CMA-ES and Differential Evolution (DE) on five challenging, real-valued functions and on one real-world problem. While CMA-ES and DE outperform ALPS on the two unimodal test functions, ALPS is much better on the three multimodal test problems and on the real-world problem. Further examination shows that, unlike the other GAs, ALPS maintains a genotypically diverse population throughout the entire search process. These findings strongly suggest that the ALPS paradigm is better able to avoid premature convergence then the other GAs.

Hornby, Gregory S.↗

Vitrification of Hanford Tank Wastes for Condensate Recycle and Feed Composition Changeover Testing (Rev. 1)

During the vitrification of Hanford Site nuclear waste at the Waste Treatment and Immobilization Plant (WTP), which is a part of the safe and efficient retrieval, treatment, and disposal mission of the U.S. Department of Energy Office of River Protection, the offgas condensate generated from the waste-to-glass conversion is currently planned to be concentrated by evaporation in the Effluent Management Facility (EMF). This concentrated condensate can then be recycled back to the incoming waste and vitrified. To test the recycle process, an apparatus was designed and built to mimic the EMF evaporator and was then used to concentrate a volume of condensate produced during the vitrification of a sample of Hanford tank 241-AP-107 (referred to herein as AP-107) waste in a continuous laboratory-scale melter (CLSM). The concentrated condensate was added to an additional sample of AP-107 waste, to mimic one round of the recycle process, and the combined solution was vitrified, producing a second round of recycle condensate. In the current study, the EMF test apparatus was used to concentrate the second-round recycle condensate under evaporation conditions (at 45 °C and 1.4 psia) designed to emulate EMF operation. The condensate was successfully concentrated by a factor of ~10 while retaining over 95 % of the technetium-99 (99Tc), Cs, and I inventories in the concentrate. Another portion of AP-107 waste was retrieved by Washington River Protection Solutions, LLC (WRPS) and transferred to Pacific Northwest National Laboratory (PNNL), where it was pretreated and then combined with the second-round recycle AP-107 condensate concentrate and glass-forming chemicals (GFCs) to form the two-time recycle AP-107 melter feed, approximating a second round to the recycling action to be performed at the WTP. A portion of AP-105 waste was also retrieved by WRPS and provided to PNNL for pretreatment and combining with GFCs to form AP-105 melter feed. The two-time recycle AP-107 and AP-105 melter feeds were processed consecutively in the CLSM. The CLSM run proceeded for 13.63 hours, producing 9.70 kg of glass for an average glass production rate of 1464 kg m 2 d -1 during the two-time recycle AP-107 feed charging and 1568 kg m 2 d -1 during the AP-105 feed charging. The rate during AP-107 charging was essentially equivalent to the rate when processing no-recycle AP-107 feed and lower than that achieved when processing one-time recycle AP-107 feed. However, all rates were within the potential range of variability when processing melter feeds with similar composition in the CLSM. Likewise, the rate during AP-105 charging was higher than the previous rate processing AP-105, but within the potential CLSM range. The cold-cap characteristics changed from the typically thin AP-107 cold cap to a foamy-edged cold cap as previously seen with AP-105 shortly after transitioning to the AP-105 melter feed. The glass produced during the CLSM run was within 10 % of its target composition for the primary glass components. The CaO and Li 2 O targets varied by more than 1 wt% between the two-time recycle AP-107 and AP-105 glass targets and it took about 2 turnovers of the CLSM glass inventory to reach a relative chemical steady state in the glass for CaO and Li 2 O after the melter feed inputs were switched.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Convection in the Physical Vapor Transport Process: Thermal - 1

The effects of convection on diffusive-convective physical vapor transport process are examined computationally. We analyze conditions ranging from typical laboratory conditions to conditions achievable only in a low gravity environment. This corresponds to thermal Rayleigh numbers Ra(sub tau) ranging from 1.80 x 10 to 1.92 x 10(exp 6). Our results indicate that the effect of the sublimation and condensation fluxes at the boundaries is to increase the threshold of instability. For typical ground based conditions, time dependent oscillatory convection can occur. This results in unsteady transport, and non-uniform temperature and concentration gradients at the crystal interface. Spectral analysis of the flow field shows parametric regions exhibiting both an oscillatory approach to steady state and a chaotic transient to a periodic state. Low gravity conditions stabilize the flow field. Convective effects are effectively reduced, thus resulting in uniform temperature and concentration gradients at the interface, a desirable condition for crystal growth.

Duval, Walter M. B.↗

Convergence Acceleration and Documentation of CFD Codes for Turbomachinery Applications

The development and analysis of turbomachinery components for industrial and aerospace applications has been greatly enhanced in recent years through the advent of computational fluid dynamics (CFD) codes and techniques. Although the use of this technology has greatly reduced the time required to perform analysis and design, there still remains much room for improvement in the process. In particular, there is a steep learning curve associated with most turbomachinery CFD codes, and the computation times need to be reduced in order to facilitate their integration into standard work processes. Two turbomachinery codes have recently been developed by Dr. Daniel Dorney (MSFC) and Dr. Douglas Sondak (Boston University). These codes are entitled Aardvark (for 2-D and quasi 3-D simulations) and Phantom (for 3-D simulations). The codes utilize the General Equation Set (GES), structured grid methodology, and overset O- and H-grids. The codes have been used with success by Drs. Dorney and Sondak, as well as others within the turbomachinery community, to analyze engine components and other geometries. One of the primary objectives of this study was to establish a set of parametric input values which will enhance convergence rates for steady state simulations, as well as reduce the runtime required for unsteady cases. The goal is to reduce the turnaround time for CFD simulations, thus permitting more design parametrics to be run within a given time period. In addition, other code enhancements to reduce runtimes were investigated and implemented. The other primary goal of the study was to develop enhanced users manuals for Aardvark and Phantom. These manuals are intended to answer most questions for new users, as well as provide valuable detailed information for the experienced user. The existence of detailed user s manuals will enable new users to become proficient with the codes, as well as reducing the dependency of new users on the code authors. In order to achieve the objectives listed, the following tasks were accomplished: 1) Parametric Study Of Preconditioning Parameters And Other Code Inputs; 2) Code Modifications To Reduce Runtimes; 3) Investigation Of Compiler Options To Reduce Code Runtime; and 4) Development/Enhancement of Users Manuals for Aardvark and Phantom

Marquart, Jed E.↗

Integrated protonic ceramic electrochemical cell for sustainable energy economy using water-energy nexus framework

Reliance on fossil fuels will continue for the next decades even though there are global pushes away from it to mitigate the overarching climate challenge, most especially by its highest consumers and availability. While there is a hastening global shift away from fossil fuel, integrating its assets into this technology helps limit the risk and future losses of stranded assets and reduce the cost of investment in the new technologies. Moreover, the generation of electricity from intermittent renewable sources like solar and wind has witnessed a significant surge in recent years, leading to a pressing demand for practical energy storage systems. Electrical energy storage is anticipated to play a pivotal role in the future global energy system, facilitating load-leveling operations to support the greater integration of renewable and distributed generation. Reversible electrochemical cells (RECs) offer a promising option for addressing the fossil fuel assets integration and energy storage challenges through the interconversion between electrical and chemical energy and concurrent utilizing carbon emission. In their electrolysis mode, the RECs convert electricity into durable, storable, and portable valuable chemical fuels such as syngas and methane. Conversely, the produced chemical fuels can be used as reactants in the fuel cell mode to generate electricity on demand with minimal (hydrocarbons) or zero when H2 or NH3 is used emissions. However, a challenging goal for this type of technology remains to achieve optimal operation and high roundtrip efficiencies, which has hindered the deployment of previous electrochemical cells. This dissertation demonstrates how reversible protonic ceramic electrochemical cells (RePCECs) can be integrated with fossil fuel power plants and renewable energy sources as a potential energy storage system. In this work, integrated RePCEC systems are designed and examined using computational modeling at scales to determine appropriate system configurations and operating conditions that achieve high roundtrip efficiencies. Cell level design of the PCEC is the first approach, several cells are assembled for the stack level model that is integrated into combined cycle powerplant and solar photovoltaic for the system level model. After critical literature review, this answered the operational and integration research questions proposed to address these challenges. The designed systems perform two functions, utilizing captured CO2 and storing renewable energy through co-electrolysis of steam and CO2. The co-electrolysis reaction involves endothermic water electrolysis and exothermic methanation reaction. To enhance high roundtrip efficiency, there is a need for thermal balance and management in the electrolysis mode. This involves operating the RePCEC stack under conditions that favor methane production to balance out heat needed by water electrolysis, it crucial for the RePCEC system operation. Methanation is enhanced by low temperatures. Leveraging on fabricated BCZYYb-electrolyte RePCEC, the cell model designed revealed that the optimum temperature for methane production is 450℃ at atmospheric pressure. Thus, to achieve optimum system performance, operating in the temperature range 450-525℃ is recommended at the given configuration, combining between the optimum temperature for methane production and temperature for the optimum stack roundtrip efficiency. Configuration with carbon capture system and purge stream is the optimum configuration from the seven conceptualized and evaluated. The modeling outcomes include a thermodynamic examination of integrated RePCEC systems, calibration of cell and stack level models, and steady-state simulation and integration into a 600MW combined cycle power plant retrofitted with two two-stage membrane-based carbon capture system and a wastewater treatment and recovery unit. At 100% powerplant loading, the stack and system roundtrip efficiencies are 72% and 51.37% respectively. Adding a purge stream for produced hydrogen at the system downstream improves the efficiencies to 74 and 55.48% respectively. At atmospheric pressure and 525℃, the system model suggests that a stack roundtrip of 82% is achievable, and overall system efficiency increases by reducing the energy consumption by the balance of plant components for steam generation and storage. Economic analysis of the process gives levelized cost of methane as $2.24/MMBtu lower than the conventional production route that range between $3.46/MMBtu and $9.85/MMBtu. The lifecycle analysis shows that the global warming potential for the production of methane and hydrogen from the RePCEC system is 3.83 kg CO2 eq which is lower than 9.35 kg CO2 eq emission during steam methane reforming for hydrogen production. This answered both the environmental and economic concerns in the raised research question. The proposed RePCEC configuration and analysis carried out in this dissertation to address the surge in renewable energy and challenges with PCEC technology hold significant potential in achieving large-scale energy storage while simultaneously reducing carbon emissions. These advancements, coupled with suitable governmental policies and incentive programs, have the potential to economically disrupt the natural gas industries by using RePCEC systems for methane production, thereby making them more favorable for eventual implementation and commercialization.

25 ENERGY STORAGE↗

Detonation to Deflagration-Mode Transitions in Pulsed Plasma Accelerators

PULSED plasma accelerators typically operate by storing energy in a capacitor bank and then discharging this energy through a gas. The current in such discharges will ionize the gas and produce a strong magnetic field, which interacts with the flowing current to accelerate the plasma through the Lorentz body force. For the present work, two plasma accelerator types employing this general scheme are of interest: the gas-fed pulsed plasma thruster (PPT) 1 and the quasi-steady magnetoplasmadynamic (MPD) accelerator. 2 The gas-fed pulsed plasma accelerator is generally understood as a completely transient device discharging in ∼1-10 μs. When the capacitor bank is discharged through the gas, a current sheet forms at the breech of the thruster and propagates forward under a j × B body force, entraining and accelerating propellant it encounters. This process is sometimes referred to in literature as 'snowplowing' the propellant or accelerating the gas in a detonation-mode because the current sheet representation approximates that of a strong detonation shockwave propagating through the gas. For these devices, acceleration of the initial current sheet ceases when either the current sheet reaches the end of the device and is ejected or when the current in the circuit reverses, striking a new 'crowbar' discharge at the breech and depriving the initial sheet of additional acceleration. In general, PPTs typically claim thrust efficiencies (ratio of jet kinetic energy to input electrical energy) registering in the teens or lower. 3 In the quasi-steady MPD accelerator, the pulse is lengthened to ∼1 ms or longer and maintained at an approximately constant level during discharge through the use of a pulse-forming network (PFN) of capacitors. After an initial transient discharge, which is typically short relative to the overall discharge period, the plasma assumes a relatively steady-state configuration, known as 'quasi-steady' MPD operation. 2 In this state, ionized gas flows through a stationary current channel in a manner that is sometimes referred to as deflagration-mode operation owing to the similarities to deflagration waves in gases. The plasma experiences electromagnetic acceleration as the plasma flows through the current channel towards the exit of the device. Quasi-steady MPD thrusters claim efficiencies up to 50% for certain propellants.4 There has been significant and sustained research over several decades on both gas-fed PPTs and quasi-steady MPD thrusters, however there have been pulsed thrusters that do not appear to exactly fit either classification, instead possessing a mixture of operational qualities characteristic of both thruster variants. The Coaxial High ENerGy (CHENG) thruster by Cheng, et al.5 operated on the short 10 μs timescales characteristic of PPTs, but claimed the high thrust densities, high efficiencies, and low electrode erosion rates that are more consistent with the MPD/deflagration mode of plasma acceleration. Gas-fed PPT research by Ziemer, et al. 3, 6 identified two separate regimes of performance in those thrusters. The regime at higher mass bits (termed Mode I in that work) possessed relatively constant thrust efficiency as a function of mass bit, while the second regime at very low mass bits (termed Mode II) exhibited an increase in efficiency with decreasing mass bit. Work by Poehlmann et al.7 and by Sitaraman and Raja8 sought to understand the performance of the CHENG thruster and the Mode I/Mode II performance in PPTs by modeling the acceleration using the Hugoniot Relation, with the detonation and deflagration modes of plasma acceleration representing two distinct sets of solutions to the relevant conservation laws. In these works, it was proposed that the values of the various controllable parameters determined whether the accelerator would operate in detonation or deflagration mode. Our hypothesized view of the acceleration process in the CHENG thruster and in PPTs experiencing a transition from Mode I to Mode II is inspired by observations of the transition from the PPT mode of operation to the quasi-steady MPD mode. Specifically, the quasi-steady MPD was discovered by driving a PPT to extended pulse lengths. Above a certain pulse length threshold the transient plasma current sheet transitions into a stable plasma acceleration mode that 'replicates in every observable detail steady flow self-field magnetoplasmadynamic acceleration.' 9 In the present work, instead of treating the accelerator as if it were only operating in a single mode during a pulse, we consider the initial stage of the discharge in all cases as a current sheet forming at the breach of the accelerator and moving towards the exit as a detonation wave. If the current sheet reaches the exit of the accelerator before the discharge is completed, the view of the acceleration mode transitions to the deflagration mode-type found in quasi-steady MPD thrusters. In previous work10 we presented a modeling framework that first captured the time-evolution of the current sheet (detonation) mode of the thruster and then transitioned into the quasi-steady MPD (deflagration) mode of plasma acceleration. In the present work, variations of the controllable parameters - specifically the pulsed circuit properties, the amount of mass injected into the thruster, and the relative timing between the initial gas injection and the initiation of the plasma current sheet - will be used to explore the thruster performance. A range of parameters are explored to demonstrate that standard gas-fed pulsed plasma accelerators, the CHENG thruster, and the quasi-steady MPD accelerator are variations of the same device, with the overall acceleration of the plasma depending upon the behavior of the plasma discharge during initial transient phase and the relative lengths of the detonation and deflagration modes of operation.

Polzin, Kurt A.↗