Results for “Model initialization” — Page 33

NASA NTRSDate not supplied

Comparison of the NAIRAS Trajectory Dose Model With ISS Measurements: Effects of Trapped Particles and Solar Energetic Particle Events.

The NAIRAS(Nowcast of Aerospace Ionizing RAdiation system) model was initially developed for fast computation of GCR and SEP events in the Earth's atmosphere. It was recently improved for computing the effects of trapped particles and its domain of validity was extended to the near space environments, which notably includes the international space station (ISS). To validate the new capabilities of the model, we compared its outputs with some measurements made in the ISS. Historical data contained several SEP events while newer measurements made with the ARMAS instrument flying in the experience bay of KIBO show the influence of trapped particles and GCR over the total dose received.

Guillaume Gronoff↗

NASA NTRS2014

Expansion of the Real-time Sport-land Information System for NOAA/National Weather Service Situational Awareness and Local Modeling Applications

The NASA Short-term Prediction Research and Transition (SPoRT) Center has been running a real-time version of the Land Information System (LIS) since summer 2010 (hereafter, SPoRTLIS). The real-time SPoRT-LIS runs the Noah land surface model (LSM) in an offline capacity apart from a numerical weather prediction model, using input atmospheric and precipitation analyses (i.e., "forcings") to drive the Noah LSM integration at 3-km resolution. Its objectives are to (1) produce local-scale information about the soil state for NOAA/National Weather Service (NWS) situational awareness applications such as drought monitoring and assessing flood potential, and (2) provide land surface initialization fields for local modeling initiatives. The current domain extent has been limited by the input atmospheric analyses that drive the Noah LSM integration within SPoRT-LIS, specifically the National Centers for Environmental Prediction (NCEP) Stage IV precipitation analyses. Due to the nature of the geographical edges of the Stage IV precipitation grid and its limitations in the western U.S., the SPoRT-LIS was originally confined to a domain fully nested within the Stage IV grid, over the southeastern half of the Conterminous United States (CONUS). In order to expand the real-time SPoRT-LIS to a full CONUS domain, alternative precipitation forcing datasets were explored in year-long, offline comparison runs of the Noah LSM. Based on results of these comparison simulations, we chose to implement the radar/gauge-based precipitation analyses from the National Severe Storms Laboratory as a replacement to the Stage IV product. The Multi-Radar Multi-Sensor (MRMS; formerly known as the National Mosaic and multi-sensor Quantitative precipitation estimate) product has full CONUS coverage at higher-resolution, thereby providing better coverage and greater detail than that of the Stage IV product. This paper will describe the expanded/upgraded SPoRT-LIS, present comparisons between the original and upgraded SPoRT-LIS, and discuss the path forward for future collaboration opportunities with SPoRT partners in the NWS.

Case, Jonathan L.↗

NASA NTRSDate not supplied

Harnessing the Digital Transformation for Development of Electrified Aircraft Propulsion Control Systems

Hybrid electric aircraft propulsion is an emerging technology that presents a variety of potential benefits along with technical integration challenges. Developing these new propulsion architectures with their complex control systems, and ultimately proving their benefit, is a multistep process that goes from concept, to analysis, to dynamic simulation, to hardware in the loop testing, to full scale testing and beyond. This effort is being revolutionized and indeed enabled by new digital tools that support the increasing technology readiness level throughout the maturation process. NASA has developed and made available a suite of digital tools that ease the path from concept to implementation. The three packages are the Toolbox for the Modeling and Analysis of Thermodynamic Systems (T-MATS), the Electrical Modeling and Thermal Analysis Toolbox (EMTAT), and the Thermal Systems Analysis Toolbox (TSAT). These tools are interactive, complementary, and compatible with each other. T-MATS is a modular thermodynamic modeling framework designed for creating custom component level models of jet engines. EMTAT is a modeling framework used to simulate a variety of power electronic devices, using both physics-based and power flow calculations. TSAT is a framework for modeling and analysis of dynamic heat transfer. These packages all consist of graphical, drag-and-drop, parameterizable building blocks representing various components of the system to be modeled, e.g., compressors, turbines, motors, energy storage devices, etc. They are designed to enable the user to model and simulate the end-to-end dynamic operation of a hybrid powertrain at the timescale of the turbomachinery, capturing mechanical, electrical, and thermal interactions. This paper will describe through multiple examples how these tools have been used successfully in several of the early stages of hybrid electric propulsion system development, from the initial system modeling to real-time interactive pilot-in-the-loop simulation to physical hardware-in-the-loop testing, each step bringing the technology closer to fruition.

Hybrid electric propulsion↗

NASA NTRSDate not supplied

Harnessing the Digital Transformation for Development of Hybrid Electric Aircraft Propulsion Control Systems

Hybrid electric aircraft propulsion is an emerging technology that presents a variety of potential benefits along with technical integration challenges. Developing these new propulsion architectures with their complex control systems, and ultimately proving their benefit, is a multistep process that goes from concept, to analysis, to dynamic simulation, to hardware in the loop testing, to full scale testing and beyond. This effort is being revolutionized and indeed enabled by new digital tools that support the increasing technology readiness level throughout the maturation process. NASA has developed and made available a suite of digital tools that ease the path from concept to implementation. The three packages are the Toolbox for the Modeling and Analysis of Thermodynamic Systems (T-MATS), the Electrical Modeling and Thermal Analysis Toolbox (EMTAT), and the Thermal Systems Analysis Toolbox (TSAT). These tools are interactive, complementary, and compatible with each other. T-MATS is a modular thermodynamic modeling framework designed for creating custom component level models of jet engines. EMTAT is a modeling framework used to simulate a variety of power electronic devices, using both physics-based and power flow calculations. TSAT is a framework for modeling and analysis of dynamic heat transfer. These packages all consist of graphical, drag-and-drop, parameterizable building blocks representing various components of the system to be modeled, e.g., compressors, turbines, motors, energy storage devices, etc. They are designed to enable the user to model and simulate the end-to-end dynamic operation of a hybrid powertrain at the timescale of the turbomachinery, capturing mechanical, electrical, and thermal interactions. This paper will describe through multiple examples how these tools have been used successfully in several of the early stages of hybrid electric propulsion system development, from the initial system modeling to real-time interactive pilot-in-the-loop simulation to physical hardware-in-the-loop testing, each step bringing the technology closer to fruition.

Electrified Propulsion↗

NASA NTRS2003

Interface Finite Elements for the Analysis of Fracture Initiation and Progression

Progressive failure analyses (PFA) are important for the prediction of residual strength and damage tolerance of vehicle structures, and to predict the energy absorbing capability of vehicle structures under crash-type loads. Typically continuum damage mechanics (CDM) and fracture mechanics (FM) are the methods used for PFA. The method of interfacial damage mechanics (IDM) is used for PFA in this research. IDM has capabilities intermediate between CDM and FM, and is used to numerically model the initiation, growth, and arrest of cracks. IDM smooths the stress singularity at the crack tip, and is easily adaptable with other nonlinearties such as plasticity and material damage. IDM is implemented by user-defined interface elements in the ABAQUS/ Standard structural analysis software package. The structural components selected to demonstrate the effectiveness PFA using interface elements are, for the most part, those with published test data. These structural components were subjected to quasi-static loading in the tests. Thus, the ABAQUS analyses are used to predict geometrically and materially nonlinear equilibrium states. Impact loading, dynamic fracture, reflected stress waves, inertia, and time dependent material behavior are not considered.

Ambur, Damodar R.↗

NASA NTRS1990

Experimental observations and finite element analysis of the initiation of fiber microbuckling in notched composite laminates

A better understanding of the factors that affect the semi-circular edge-notched compressive strength is developed, and the associated failure mode(s) of thermoplastic composite laminates with multidirectional stacking sequences are identified. The primary variables in this investigation are the resin nonlinear shear constitutive behavior, stacking sequence (orientation of plies adjacent to the 0 degree plies), resin-rich regions between the 0 degree plies and the off-axis supporting plies, fiber/matrix interfacial bond strength, and initial fiber waviness. Two thermoplastic composite material systems are used in this investigation. The materials are the commercial APC-2 (AS4/PEEK) and a poor interface experimental material, AU4U/PEEK, designed for this investigation. Notched compression specimens are studied at 21, 77, and 132 C. Geometric and material nonlinear two-dimensional finite element analysis is used to model the initiation of fiber microbuckling of both the ideal straight fiber and the more realistic initially wavy fiber. The effects of free surface, fiber constitutive properties, matrix constitutive behavior, initial fiber curvature, and fiber/matrix interfacial bond strength on fiber microbuckling initiation strain levels are considered.

Guynn, E. Gail↗

NASA NTRS1977

The multiple nesting of mesoscale, submesoscale, and nonhydrostatic microscale numerical models for a case of mesocyclogenesis and severe storm development

A fundamental approach involving the integration of two different numerical models over several different spatial mesh lengths is considered. Model I is a hydrostatic, incompressible, and adiabatic primitive (Newtonian) equation set with no boundary layer forcing, i.e., surface heat or moisture flux or external frictional stresses. Typically this model is initialized from conventional radiosonde data and run for a 12-15 hour forecast period over a 40 km horizontal mesh. If phase I of the tornado developmental processes develops, the data is 'frozen' in space and time, interpolated to a 20 km mesh, and integrated for a shorter time period (3 hours) with a shorter time step. Model II is a nonhydrostatic compressible 18-level complement of model I. If model II indicates the development of the intensification of phase V dynamical processes, the data is frozen and interpolated to finer mesh lengths until the tornadic circulation system is approximated

Kaplan, M. L.↗

NASA NTRS1986

Vertical normal modes of a mesoscale model using a scaled height coordinate

Vertical modes were derived for a version of the Colorado State Regional Atmospheric Mesoscale Modeling System. The impacts of three options for dealing with the upper boundary of the model were studied. The standard model formulation holds pressure constant at a fixed altitude near the model top, and produces a fastest mode with a speed of about 90 m/sec. An alternative formulation, which allows for an external mode, could require recomputation of vertical modes for every surface elevation on the horizontal grid unless the modes are derived in a particular way. These results have bearing on the feasibility of applying vertical mode initialization to models with scaled height coordinates.

Lipton, A. E.↗

NASA NTRS1980

Spherical harmonic analysis of a synoptic climatology generated with a global general circulation model

Spherical harmonic analysis was used to analyze the observed climatological (C) fields of temperature at 850 mb, geopotential height at 500 mb, and sea level pressure. The spherical harmonic method was also applied to the corresponding "model climatological" fields (M) generated by a general circulation model, the "GISS climate model." The climate model was initialized with observed data for the first of December 1976 at 00. GMT and allowed to generate five years of meteorological history. Monthly means of the above fields for the five years were computed and subjected to spherical harmonic analysis. It was found from the comparison of the spectral components of both sets, M and C, that the climate model generated reasonable 500 mb geopotential heights. The model temperature field at 850 mb exhibited a generally correct structure. However, the meridional temperature gradient was overestimated and overheating of the continents was observed in summer.

Christidis, Z. D.↗

NASA NTRS1988

Modelling a C-type flare observed in microwaves and hard X-rays

Using the very large array (VLA) at 6 and 20 cm wavelength and the hard X-ray burst spectrometer on the Solar Maximum Mission, a two-ribbon flare was observed from the onset phase through the maximum and decline on November 14, 1981. Because of the extensive size of the microwave source and the gradual variations in hard X-rays whose spectrum becomes progressively flatter with time, the flare is classified as a C-type flare. Considering the hardening of the X-ray spectrum and its non-impulsive nature, a coronal trap model was invoked for the energetic electrons. The microwave emission is easily accounted for by gyrosynchronous radiation from mildly relativistic electrons. It was found that the source must be optically thick at 20 cm during the maximum phase, but as the source evolved toward an optically thin regime, the intensity decreased while the degree of circular polarization increased. In an initial homogeneous model, we found that the computed microwave spectrum was too narrow to match the patrol spectrum from 606 to 15400 MHz. In the model, the magnetic field consists of a dipolar arcade bridging the H alpha ribbons, and extending to heights of order 40,000 to 50,000 km. The variation of the magnetic field strength from footpoints to apex causes the gyrosynchrotron spectrum to be broader. Preliminary conclusions regarding the electron distributions producing the hard X-rays and the microwaves, and the suitability of this model for C-type flares is presented.

Schmahl, E. J.↗

NASA NTRS2017

Daily Kilometer-Scale MODIS Satellite Maps of PM2.5 Describe Wintertime Episodes

The San Joaquin Valley (SJV) suffers from severe health-endangering episodes of PM2.5 aerosol loadings in wintertime; episodes last approximately 5 days and differ in geographical distribution and composition. PM2.5 stations are scattered; consequently the use of remote sensing to map variable regional patterns of these varying respirable aerosol concentrations is desirable. High-precision AOT retrievals can capture column particulate loading. However,PM2.5 mapping is challenging due to several reasons: particularly thin mixed layers (ML) and thus relatively low aerosol optical thickness (AOT) close to current measurement limits, variable and a typical composition of the aerosols, and complex surface bidirectional reflectance. However, the West does present some advantages in analysis. Air basins are isolated from long-distance transport, and experience predominant strong meteorological subsidence. Thus these Western basin regions have fewer problematic cases of overriding aerosol layers detached from the surface. To counter such local overriding, Chu et al. have described an approach for the Eastern US, and He et al have described a synoptic classification approach useful in Shanghai. The Bay Area Air Quality Management District (BAAQMD) expands our experience with the use of AOT, with lower PM2.5 and several isolated sub-basins. We have prepared daily maps of episodes in each region. We present also a sequence of increasingly detailed statistical models, AOT initially appears to contribute little information; however, inclusion of weather information reveals its utility. Lyapustin and Wang's MultiAngle Implementation of Atmospheric Correction (MAIAC) retrieval for AOT provided the most useful operational remote sensing information for these regions. It provides high (1-km) spatial resolution maps and a high percentage of availability. Empirical regression methods have found that random effects regression models (aka mixed effects models, ME) employing AOT provide good estimates of ground PM2.5 concentrations.Here, we attempt to extend these methods and evaluate the usefulness of AOT with greater physical analysis, based on DISCOVER-AQ4 experience.

Chatfield, Robert B.↗

NASA NTRS1994

Bias and high-order galaxy correlation functions in the APM galaxy survey

On large scales, the higher order moments of the mass distribution, S(sub J) = bar-zeta(sub J)/bar-zeta(sup J-1)(sub 2), e.g., the skewness S(sub 3) and kurtosis S(sub 4), can be predicted using nonlinear perturbation theory. Comparison of these predictions with moments of the observed galaxy distribution probes the bias between galaxies and mass. Applying this method to models with initially Gaussian fluctuations and power spectra P(k) similar to that of galaxies in the Automatic Plate Measuring (APM) survey, we find that the predicted higher order moments S(sub J)(R) are in good agreement with those directly inferred from the APM survey in the absence of bias. We use this result to place limits on the linear and nonlinear bias parameters. Models in which the extra power observed on large scales (with respect to the standard cold dark matter (CDM) model) is produced by scale-dependent bias match the APM higher order amplitudes only if nonlinear bias (rather than nonlinear gravity) generates the observed higher order moments. When normalized to Cosmic Background Explorer Differential Microwave Radiometer (COBE DMR), these models are siginificantly ruled out by the S(sub 3) observations. The cold plus hot dark matter model normalized to COBE can reproduce the APM higher order correlations if one introduces nonlinear bias terms, while the low-density CDM model with a cosmological constant does not require any bias to fit the large-scale amplitudes.

Gaztanaga, Enrique↗

NASA NTRS2005

Uncertainty Evaluation of Computational Model Used to Support the Integrated Powerhead Demonstration Project

NASA and the U.S. Air Force are working on a joint project to develop a new hydrogen-fueled, full-flow, staged combustion rocket engine. The initial testing and modeling work for the Integrated Powerhead Demonstrator (IPD) project is being performed by NASA Marshall and Stennis Space Centers. A key factor in the testing of this engine is the ability to predict and measure the transient fluid flow during engine start and shutdown phases of operation. A model built by NASA Marshall in the ROCket Engine Transient Simulation (ROCETS) program is used to predict transient engine fluid flows. The model is initially calibrated to data from previous tests on the Stennis E1 test stand. The model is then used to predict the next run. Data from this run can then be used to recalibrate the model providing a tool to guide the test program in incremental steps to reduce the risk to the prototype engine. In this paper, they define this type of model as a calibrated model. This paper proposes a method to estimate the uncertainty of a model calibrated to a set of experimental test data. The method is similar to that used in the calibration of experiment instrumentation. For the IPD example used in this paper, the model uncertainty is determined for both LOX and LH flow rates using previous data. The successful use of this model is then demonstrated to predict another similar test run within the uncertainty bounds. The paper summarizes the uncertainty methodology when a model is continually recalibrated with new test data. The methodology is general and can be applied to other calibrated models.

Steele, W. G.↗

NASA NTRS2023

Harnessing the Digital Transformation for Development of Electrified Aircraft Propulsion Control Systems

Hybrid electric aircraft propulsion is an emerging technology that presents a variety of potential benefits along with technical integration challenges. Developing these new propulsion architectures with their complex control systems, and ultimately proving their benefit, is a multistep process. This process includes concept development and analysis, dynamic simulation, hardware-in-the-loop testing, full-scale testing, and so on. This effort is being revolutionized and indeed enabled by new digital tools that support increasing the technology readiness level throughout the maturation process. As part of this Digital Transformation, NASA has developed a suite of publicly available digital tools that facilitate the path from concept to implementation. This paper describes the NASA-developed tools and puts them in the context of control system development for hybrid electric aircraft propulsion. The three MATLAB®-based software packages are the Toolbox for the Modeling and Analysis of Thermodynamic Systems (T-MATS), the Electrical Modeling and Thermal Analysis Toolbox (EMTAT), and the Thermal Systems Analysis Toolbox (TSAT). These tools are interactive, complementary, and compatible with each other. T-MATS is a modular thermodynamic modeling framework designed for creating custom component level models of jet engines. EMTAT is a modeling framework used to simulate a variety of power electronic devices, using both physics-based and power flow calculations. TSAT is a framework for modeling and analysis of dynamic heat transfer. These packages all consist of graphical, drag-and-drop, parameterizable building blocks representing various components of the system to be modeled, e.g., compressors, turbines, motors, energy storage devices, etc. They are designed to enable the user to model and simulate the end-to-end dynamic operation of a hybrid electric gas turbine engine powertrain at the timescale of the turbomachinery, capturing mechanical, electrical, and thermal interactions. This paper demonstrates through multiple examples how these tools have been used successfully in a variety of applications, including several of the early stages of hybrid electric gas turbine engine propulsion system development, from the initial system modeling to real-time interactive pilot-in-the-loop simulation to physical hardware-in-the-loop testing, each step bringing the technology closer to fruition.

hybrid electric propulsion↗

NASA NTRS2017

Monte Carlo Computational Modeling of Atomic Oxygen Interactions

Computational modeling of the erosion of polymers caused by atomic oxygen in low Earth orbit (LEO) is useful for determining areas of concern for spacecraft environment durability. Successful modeling requires that the characteristics of the environment such as atomic oxygen energy distribution, flux, and angular distribution be properly represented in the model. Thus whether the atomic oxygen is arriving normal to or inclined to a surface and whether it arrives in a consistent direction or is sweeping across the surface such as in the case of polymeric solar array blankets is important to determine durability. When atomic oxygen impacts a polymer surface it can react removing a certain volume per incident atom (called the erosion yield), recombine, or be ejected as an active oxygen atom to potentially either react with other polymer atoms or exit into space. Scattered atoms can also have a lower energy as a result of partial or total thermal accommodation. Many solutions to polymer durability in LEO involve protective thin films of metal oxides such as SiO2 to prevent atomic oxygen erosion. Such protective films also have their own interaction characteristics. A Monte Carlo computational model has been developed which takes into account the various types of atomic oxygen arrival and how it reacts with a representative polymer (polyimide Kapton H) and how it reacts at defect sites in an oxide protective coating, such as SiO2 on that polymer. Although this model was initially intended to determine atomic oxygen erosion behavior at defect sites for the International Space Station solar arrays, it has been used to predict atomic oxygen erosion or oxidation behavior on many other spacecraft components including erosion of polymeric joints, durability of solar array blanket box covers, and scattering of atomic oxygen into telescopes and microwave cavities where oxidation of critical component surfaces can take place. The computational model is a two dimensional model which has the capability to tune the interactions of how the atomic oxygen reacts, scatters, or recombines on polymer or nonreactive surfaces. In addition to the specification of atomic oxygen arrival details, a total of 15 atomic oxygen interaction parameters have been identified as necessary to properly simulate observed interactions and resulting polymer erosion that have been observed in LEO. The tuning of the Monte Carlo model has been accomplished by adjusting interaction parameters so the erosion patterns produced by the model match those from several actual LEO space experiments. Surface texturing in LEO can also be predicted by the model. Such comparison of space tests with ground laboratory experiments have enabled confidence in ground laboratory lifetime prediction of protected polymers. Results of Monte Carlo tuning, examples of surface texturing and undercutting erosion prediction, and several examples of how the model can be used to predict other LEO and Mars orbital space results are presented.

Monte Carlo Modeling↗

NASA NTRS2017

Monte Carlo Computational Modeling of Atomic Oxygen Interactions

Computational modeling of the erosion of polymers caused by atomic oxygen in low Earth orbit (LEO) is useful for determining areas of concern for spacecraft environment durability. Successful modeling requires that the characteristics of the environment such as atomic oxygen energy distribution, flux, and angular distribution be properly represented in the model. Thus whether the atomic oxygen is arriving normal to or inclined to a surface and whether it arrives in a consistent direction or is sweeping across the surface such as in the case of polymeric solar array blankets is important to determine durability. When atomic oxygen impacts a polymer surface it can react removing a certain volume per incident atom (called the erosion yield), recombine, or be ejected as an active oxygen atom to potentially either react with other polymer atoms or exit into space. Scattered atoms can also have a lower energy as a result of partial or total thermal accommodation. Many solutions to polymer durability in LEO involve protective thin films of metal oxides such as SiO2 to prevent atomic oxygen erosion. Such protective films also have their own interaction characteristics. A Monte Carlo computational model has been developed which takes into account the various types of atomic oxygen arrival and how it reacts with a representative polymer (polyimide Kapton H) and how it reacts at defect sites in an oxide protective coating, such as SiO2 on that polymer. Although this model was initially intended to determine atomic oxygen erosion behavior at defect sites for the International Space Station solar arrays, it has been used to predict atomic oxygen erosion or oxidation behavior on many other spacecraft components including erosion of polymeric joints, durability of solar array blanket box covers, and scattering of atomic oxygen into telescopes and microwave cavities where oxidation of critical component surfaces can take place. The computational model is a two dimensional model which has the capability to tune the interactions of how the atomic oxygen reacts, scatters, or recombines on polymer or nonreactive surfaces. In addition to the specification of atomic oxygen arrival details, a total of 15 atomic oxygen interaction parameters have been identified as necessary to properly simulate observed interactions and resulting polymer erosion that have been observed in LEO. The tuning of the Monte Carlo model has been accomplished by adjusting interaction parameters so the erosion patterns produced by the model match those from several actual LEO space experiments. Surface texturing in LEO can also be predicted by the model. Such comparison of space tests with ground laboratory experiments have enabled confidence in ground laboratory lifetime prediction of protected polymers. Results of Monte Carlo tuning, examples of surface texturing and undercutting erosion prediction, and several examples of how the model can be used to predict other LEO and Mars orbital space results are presented.

Atomic Oxygen↗

NASA NTRS2004

Sensor Web Dynamic Measurement Techniques and Adaptive Observing Strategies

Sensor Web observing systems may have the potential to significantly improve our ability to monitor, understand, and predict the evolution of rapidly evolving, transient, or variable environmental features and events. This improvement will come about by integrating novel data collection techniques, new or improved instruments, emerging communications technologies and protocols, sensor mark-up languages, and interoperable planning and scheduling systems. In contrast to today's observing systems, "event-driven" sensor webs will synthesize real- or near-real time measurements and information from other platforms and then react by reconfiguring the platforms and instruments to invoke new measurement modes and adaptive observation strategies. Similarly, "model-driven" sensor webs will utilize environmental prediction models to initiate targeted sensor measurements or to use a new observing strategy. The sensor web concept contrasts with today's data collection techniques and observing system operations concepts where independent measurements are made by remote sensing and in situ platforms that do not share, and therefore cannot act upon, potentially useful complementary sensor measurement data and platform state information. This presentation describes NASA's view of event-driven and model-driven Sensor Webs and highlights several research and development activities at the Goddard Space Flight Center.

Talabac, Stephen J.↗

NASA NTRS2018

Acceleration Modes and 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, ionizing and accelerating it through the Lorentz body force. Two plasma accelerator types employing this general scheme have typically been studied: the gas-fed pulsed plasma thruster and the quasi-steady magnetoplasmadynamic (MPD) accelerator. The gas-fed pulsed plasma accelerator is generally represented as a completely transient device discharging in approximately 1-10 microseconds. 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 (current density) by B (magnetic field) body force, entraining propellant it encounters. This process is sometimes referred to as detonation-mode acceleration because the current sheet representation approximates that of a strong shock propagating through the gas. 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 current sheet at the breech and depriving the initial sheet of additional acceleration. In the quasi-steady MPD accelerator, the pulse is lengthened to approximately 1 millisecond or longer and maintained at an approximately constant level during discharge. The time over which the transient phenomena experienced during startup typically occur is short relative to the overall discharge time, which is now long enough for the plasma to assume a relatively steady-state configuration. The ionized gas flows through a stationary current channel in a manner that is sometimes referred to as the deflagration-mode of operation. The plasma experiences electromagnetic acceleration as it flows through the current channel towards the exit of the device. A device that had a short pulse length but appeared to operate in a plasma acceleration regime different from the gas-fed pulsed plasma accelerators was developed by Cheng, et al. The Coaxial High ENerGy (CHENG) thruster operated on the 10-microseconds timescales of pulsed plasma thrusters, but claimed high thrust density, high efficiency and low electrode erosion rates, which are more consistent with the deflagration mode of acceleration. Separate work on gas-fed pulsed plasma thrusters (PPTs) by Ziemer, et al. identified two separate regimes of performance. The regime at higher mass bits (termed Mode I in that work) possessed relatively constant thrust efficiency (ratio of jet kinetic energy to input electrical energy) as a function of mass bit. In the second regime at very low mass bits (termed Mode II), the efficiency increased with decreasing mass bit. Work by Poehlmann et al. and by Sitaraman and Raja 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 representing two distinct sets of solutions to the relevant conservation laws. These works studied the proposal that, depending upon the values of the various controllable parameters, the accelerator would operate in either the detonation or deflagration mode. In the present work, we propose a variation on the explanation for the differences in performance between the various pulsed plasma accelerators. Instead of treating the accelerator as if it were only operating in one mode or the other during a pulse, we model the initial stage of the discharge in all cases as an accelerating current sheet (detonation mode). If the current sheet reaches the exit of the accelerator before the discharge is completed, the acceleration mode transitions to the deflagration mode type found in the quasi-steady MPD thrusters. This modeling method is used 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.↗