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 451 records · Page 25

Quantifying Leak Rates and Vacuum Impact on Granular Soil Column Pneumatic Seal

Introduction: All In-Situ Resource Utilization (ISRU) regolith reactor systems that aim to extract gases or vapors must maintain a seal along the flow path of the regolith from the inlet hopper to where extracted gases or vapors are generated. Using a column of regolith at the inlet has been proposed as an alternative to systems using series of mechanical valves to reduce mass and complexity. Preliminary investigations into the use of regolith itself as a barrier to escaping fluids have been conducted, which this research seeks to expand upon. Different vertical regolith simulant plug diameters and heights were tested under vacuum with dry nitrogen to observe leak rates and approximate steady state pressure gradient equilibria. The effect of absolute ambient pressure was also studied by repeating tests at incremental pressures ranging from 10e-5 Torr up to 1 atm to corroborate ambient bench tests from previous work. The impact of regolith dynamics was addressed to emulate continuous-feed regolith reaction processes by measuring leak rates through continuously draining simulant columns. OB-1 lunar highlands regolith simulant was used. Experimental Setup: Regolith simulant columns were constructed using a motorized butterfly valve separating upper acrylic tubes of various diameters from a larger, lower aluminum pipe plenum (See Figure 1). The butterfly valve was cracked open such that gas could freely pass through whereas regolith would bridge and stay suspended as a cohesive plug above in the acrylic portion. These assemblies were loaded into a vacuum chamber, and a compressed dry nitrogen line was routed into the lower plenum. Pressure transducers measured the absolute pressure of the lower plenum, which was compared to absolute chamber pressure to monitor the pressure gradient across the regolith simulant. A precision mass flow controller enabled fine control of plenum pressure. The plenum pressure was quasi-statically increased to set points informed by previous work [1] before shutting off the gas supply to allow the now pressurized gas volume to dissipate through the simulant into the chamber in what was referred to as a ‘leak test.’ In a separate set of tests, the pressure was allowed to quasi-statically build until visual disturbances in the soil were observed via remote camera feeds. Initial leak rate tests were repeated on the disturbed soil to observe changes in leak rates between initial as-poured and post-bubbling soil conditions. This test schedule was applied to simulant columns of three different diameters (approximately 2, 4, and 6 inches) as well as two different heights (around 8 and 16 inches) and was repeated at least three times for each configuration. Plenums were also pressurized while the butterfly valves were opened slightly further to facilitate slow, controlled, constant flowrate draining of the simulant to observe the impact that continuous-feed systems might have on leak rate through the regolith column. Finally, the regolith columns were reset and refilled and static leak rate tests were repeated approximately every decade in ambient pressure from 10e-5 Torr up to ambient conditions to observe the impact of transitioning out of the molecular flow regime on regolith column seal performance. Analysis: Trends with respect to each factor examined were constructed and reported alongside commentary on the most likely responsible physical phenomena. Confidence intervals were applied to the predictive findings of the data by quantifying instrumentation error and accounting for probabilistic impacts of curve fitting the generic exponential function onto pressure decay data. The leak rate is assumed directly proportional to the exponential term coefficient, while the steady-state static pressure gradient is assumed to be the asymptotic limit of each curve fit. Finally, a generic model developed by Ogino et al. [2] and later tuned for ambient testing was re-tuned using vacuum data to estimate leak rates on the lunar surface. Results: Pressure decay rates decreased up to 27-fold under vacuum when compared to ambient pressure, as shown in Figure 2. Simulant column diameter was found to exponentially correlate to both regolith column leak rate and steady-state equilibrium pressure delta. Height was found to be linearly related to leak rate –doubling the regolith column height consistently yielded around a 40% reduction in leak rate. Dynamic regolith draining was found to slightly increase leak rate, although to a lower degree than expected. Conclusion: Because of the transition into molecular gas flow through regolith, vacuum testing is required to properly estimate performance. Regolith columns should be designed as narrow as possible to avoid bridging and minimize leak rate, and as tall as is feasible to both further reduce leak rate and to help prevent bubbling, which has been shown to degrade the quality of the regolith seal. Disrupted soil was found to exhibit higher leak rates than pristine, as-poured soil, indicating the need to consider the impact of recovery from depressurization events. Regolith columns appear equally feasible for both dosed and continuous-feed reactors.

ISRU↗

Quasi-static magnetospheric MHD processes and the 'ground state' of the magnetosphere

The possible existence of a theoretically well-defined ground state of the earth's magnetosphere is discussed in the context of magnetic equilibrium and convection theories. The quasi-static MHD theory is reviewed and a two-dimensional model is used to show that the convecting magnetosphere can reach a steady state. It is suggested that, under the influence of convection, magnetic substorms occur periodically in the magnetosphere and are an integral part of the entire convection cycle. The concept of the ground state of the magnetosphere is examined. The term 'average magnetosphere', is proposed for defining a baseline configuration that corresponds to average solar wind conditions.

Voigt, G.-H.↗

Venusian hydrology: Steady state reconsidered

In 1987, Grinspoon proposed that the data on hydrogen abundance, isotopic composition, and escape rate were consistent with the hypothesis that water on Venus might be in steady state rather than monotonic decline since the dawn of time. This conclusion was partially based on a derived water lifetime against nonthermal escape of approximately 10(exp 8) yr. De Bergh et al., preferring the earlier Pioneer Venus value of 200 ppm water to the significantly lower value detected by Bezard et al., found H2O lifetimes of greater than 10(exp 9) yr. Donahue and Hodges derived H2O lifetimes of 0.4-5 x 10 (exp 9) yr. Both these analyses used estimates of H escape flux between 0.4 x 10(exp 7) and 1 x 10(exp 7) cm(exp -2)s(exp -1) from Rodriguez et al. Yet in more recent Monte Carlo modeling, Hodges and Tinsley found an escape flux due to charge exchange with hot H(+) of 2.8 x 10(exp 7) cm(exp -2)s(exp -1). McElroy et al. estimated an escape flux of 8 x 10(exp 6) cm(exp -2)s(exp -1) from collisions with hot O produced by dissociative recombination of O2(+). Brace et al. estimated an escape flux of 5 x 10(exp 6) cm(exp -2)s(exp -1) from ion escape from the ionotail of Venus. The combined estimated escape flux from all these processes is approximately 4 x 10(exp 7) cm(exp -2)s(exp -1). The most sophisticated analysis to date of near-IR radiation from Venus' nightside reveals a water mixing ratio of approximately 30 ppm, suggesting a lifetime against escape for water of less than 10(exp 8) yr. Large uncertainties remain in these quantities, yet the data point toward a steady state. Further evaluation of these uncertainties, and new evolutionary modeling incorporating estimates of the outgassing rate from post-Magellan estimates of the volcanic resurfacing rate are presented.

Grinspoon, David H.↗

Markov reward processes

Numerous applications in the area of computer system analysis can be effectively studied with Markov reward models. These models describe the behavior of the system with a continuous-time Markov chain, where a reward rate is associated with each state. In a reliability/availability model, upstates may have reward rate 1 and down states may have reward rate zero associated with them. In a queueing model, the number of jobs of certain type in a given state may be the reward rate attached to that state. In a combined model of performance and reliability, the reward rate of a state may be the computational capacity, or a related performance measure. Expected steady-state reward rate and expected instantaneous reward rate are clearly useful measures of the Markov reward model. More generally, the distribution of accumulated reward or time-averaged reward over a finite time interval may be determined from the solution of the Markov reward model. This information is of great practical significance in situations where the workload can be well characterized (deterministically, or by continuous functions e.g., distributions). The design process in the development of a computer system is an expensive and long term endeavor. For aerospace applications the reliability of the computer system is essential, as is the ability to complete critical workloads in a well defined real time interval. Consequently, effective modeling of such systems must take into account both performance and reliability. This fact motivates our use of Markov reward models to aid in the development and evaluation of fault tolerant computer systems.

Smith, R. M.↗

The thermal structure and thermal evolution of the continental lithosphere

The thermal structure and evolution of the continental lithosphere are examined. Surface heat flow data and the factors which modify them are addressed, and the diversity of thermal phenomena in the lithosphere is discussed in the framework of plate interactions. The lithosphere is divided into three sections for the purposes of discussion. In the upper, near-surface zone, temperatures can be strongly affected by near-surface processes, which must be taken into account in the measurement and evaluation of surface heat flow. The thermal structure of the middle, internal zone of the lithosphere responds to the heat balance and thermal properties of the lithosphere, which define its steady state thermal structure. Internal deformation and magmatic intrusion within this zone, and interaction between the lithosphere and the asthenosphere in the lower boundary zone of the lithosphere cause transient thermal disturbances in the lithosphere. The criteria for defining the base of the thermal lithosphere are briefly discussed.

Morgan, P.↗

Quasi-static MHD processes in earth's magnetosphere

An attempt is made to use the MHD equilibrium theory to describe the global magnetic field configuration of earth's magnetosphere and its time evolution under the influence of magnetospheric convection. To circumvent the difficulties inherent in today's MHD codes, use is made of a restriction to slowly time-dependent convection processes with convective velocities well below the typical Alfven speed. This restriction leads to a quasi-static MHD theory. The two-dimensional theory is outlined, and it is shown how sequences of two-dimensional equilibria evolve into a steady state configuration that is likely to become tearing mode unstable. It is then concluded that magnetospheric substorms occur periodically in earth's magnetosphere, thus being an integral part of the entire convection cycle.

Voigt, Gerd-Hannes↗

Navier-Stokes simulation of the supersonic combustion flowfield in a ram accelerator

A computational study of the ram accelerator, a ramjet-in-tube device for accelerating projectiles to ultrahigh velocities, is presented. The analysis is performed using a fully implicit TVD scheme that efficiently solves the Reynolds-averaged Navier-Stokes equations and the species continuity equations associated with a finite rate combustion model. Previous analyses of this concept were based on inviscid assumptions. The present results indicate that viscous effects are of primary importance; in all the cases studied, shock-induced combustion always started in the boundary layer. The effects of Mach number, mixture composition, pressure, and turbulence are investigated for various configurations. Two types of combustion processes, one stable and the other unstable, were observed depending on the inflow conditions. In the unstable case, a detonation wave is formed, which propagates upstream and unstarts the ram accelerator. In the stable case, a solution that converges to steady-state is obtained, in which the combustion wave remains stationary with respect to the ram accelerator projectile. The possibility of stabilizing the detonation wave by means of a backward facing step is also investigated. In addition to these studies, two numerical techniques were tested. These two techniques are vector extrapolation to accelerate convergence, and a diagonal formulation that eliminates the expense of inverting large block matrices that arise in chemically reacting flows.

Yungster, Shaye↗

Characterization of Urban Air Mobility Vehicle Operational Noise and Community Noise Impact

This presentation focuses on the process for predicting urban air mobility (UAM) vehicle noise at the source and how that can be used to estimate and mitigate the impact on the community. Following conceptual design, in which the vehicle is appropriately sized for its intended mission, a comprehensive analysis must be performed for a range of operating conditions spanning the flight envelope to determine the corresponding configurations of the vehicle, that is, the trimmed states. Many UAM vehicles have redundant controls, so the trimmed state for any particular operating condition may not be unique, with some states producing more noise than others. For each trimmed state, the noise produced by each source, for example, steady and unsteady rotor noise, inclusive of propulsion airframe aeroacoustic effects, may be computed and so-called source noise (hemi)spheres generated. The data may be used for analyses supporting noise certification or serve as input to auralizations (turning the numerical data into an audible sound) that can subsequently be used as part of a perception-influenced acoustic design process. Land use planning tools, including those using simulation and integrated modeling approaches, also use these (or derived) data to generate noise exposure maps on the ground. Finally, the data may be used to support development of low noise flight operations using an acoustic flight simulator or acoustically aware flight control.

urban air mobility↗

Internal strain analysis of ceramics using scanning laser acoustic microscopy

Quantitative studies of material behavior characteristics are essential for predicting the functionality of a material under its operating conditions. A nonintrusive methodology for measuring the in situ strain of small dimeter (to 11 microns) ceramic fibers under uniaxial tensile loading and the local internal strains of ceramics and ceramic composites under flexural loading is introduced. The strain measurements and experimentally observed mechanical behavior are analyzed in terms of the microstructural development and fracture behavior of each test specimen evaluated. Measurement and analysis of Nicalon silicon carbide (SiC) fiber (15 microns diameter) indicate that the mean elastic modulus of the individual fiber is 185.3 GPa. Deviations observed in the experimentally determined elastic modulus values between specimens were attributed to microstructural variations which occur during processing. Corresponding variations in the fracture surface morphology were also observed. The observed local mechanical behavior of a lithium alumino-silicate (LAS) glass ceramic, a LAS/SiC monofilament composite, and a calcium alumino-silicate (CAS)/SiC fully reinforced composite exhibits nonlinearities and apparent hysteresis due to the subcritical mechanical loading. Local hysteresis in the LAS matrices coincided with the occurrence of multiple fracture initiation sites, localized microcracking, and secondary cracking. The observed microcracking phenomenon was attributed to stress relaxation of residual stresses developed during processing, and local interaction of the crack front with the microstructure. The relaxation strain and stress predicted on apparent mechanical hysteresis effects were defined and correlated with the magnitude of the measured fracture stress for each specimen studied. This quantitative correlation indicated a repeatable measure of the stress at which matrix microcracking occurred for stress relief of each material system. Stress relaxation occurred prior to the onset of steady state cracking conditions. The relaxation stress occurred at 18.5 percent of the fracture stress in LAS and 11.0 percent of the yield stress in CAS/SiC. The relaxation stress ratio was dependent upon the dominant fracture mode of the LAS/SiC specimens. Relaxation stress ratios greater than 0.30 were observed for specimens which fractured due to shear at the fiber matrix interface; specimens which fracture due to tensile cracking had relaxation stress ratios less than 0.30. The stress relaxation ratio appeared to be a specific characteristic of the glass ceramic material. The measured stress relaxation for LAS indicated a measure of the inherent residual stresses in the material due to processing and suggested localized toughening mechanisms for brittle material structures.

Kent, Renee M.↗

A numerical simulation of the NFAC (National Full-scale Aerodynamics Complex) open-return wind tunnel inlet flow

The flow into an open return wind tunnel inlet was simulated using Euler equations. An explicit predictor-corrector method was employed to solve the system. The calculation is time-accurate and was performed to achieve a steady-state solution. The predictions are in reasonable agreement with the experimental data. Wall pressures are accurately predicted except in a region of recirculating flow. Flow-field surveys agree qualitatively with laser velocimeter measurements. The method can be used in the design process for open return wind tunnels.

Kaul, U. K.↗

A numerical simulation of the NFAC (National Full-scale Aerodynamics Complex) open-return wind tunnel inlet flow

The flow into an open return wind tunnel inlet was simulated using Euler equations. An explicit predictor-corrector method was employed to solve the system. The calculation is time-accurate and was performed to achieve a steady-state solution. The predictions are in reasonable agreement with the experimental data. Wall pressures are accurately predicted except in a region of recirculating flow. Flow-field surveys agree qualitatively with laser velocimeter measurements. The method can be used in the design process for open-return wind tunnels.

Kaul, U. K.↗

Fundamental processes in pair plasmas

The various processes that produce and destroy electron-positron pairs are reviewed, and the timescales of these processes are compared to thermalization, accretion, and cooling timescales. The various radiation spectra produced by relativistic, thermal plasmas are considered. Recent results for the equilibria available to finite, thermal relativistic plasmas with and without embedded magnetic fields are reviewed. Such plasmas, in steady state, have maximum temperatures, luminosities, and field strengths, useful diagnostics for interpreting quasars and active galaxies.

Lightman, A. P.↗

An experimental study using flow visualization on the effect of an acoustic field on heat transfer from spheres

The physical mechanisms responsible for the heat transfer process in a thermal-acoustic field were investigated using the technique of holographic interferometry for flow visualization. Experimental results were obtained with sound pressure levels in the range of 120 to 150 decibels, relative to a pressure of 0.0002 dynes/sq cm. Steady state laminar flow was observed when the vibrational Reynolds number was below 400; separated flow was observed when it was above 400. In the presence of a horizontal sound field, the data indicate that the relation between the vibrational Nusselt number, Nu(v) and the vibrational Reynolds number, Re(v) is given by Nu(v) = Re(v) exp 0.22. In the presence of a vertical sound field, the corresponding relation is Nu(v) = Re(v) exp 0.15.

Leung, W. W.↗

Rocket engine failure detection using system identification techiques

The theoretical foundation and application of two univariate failure detection algorithms to Space Shuttle Main Engine (SSME) test firing data is presented. Both algorithms were applied to data collected during steady state operation of the engine. One algorithm, the time series algorithm, is based on time series techniques and involves the computation of autoregressive models. Times series techniques have been previously applied to SSME data. The second algorithm is based on standard signal processing techniques. It consists of tracking the variations in the average signal power with time. The average signal power algorithm is a newly proposed SSME failure detection algorithm. Seven nominal test firings were used to develop failure indication thresholds for each algorithm. These thresholds were tested using four anomalous firings and one additional nominal firing. Both algorithms provided significantly earlier failure indication times than did the current redline limit system. Neither algorithm gave false failure indications for the nominal firing. The strengths and weaknesses of the two algorithms are discussed and compared. The average signal algorithm was found to have several advantages over the time series algorithm.

Meyer, Claudia M.↗

Rocket engine failure detection using system identification techniques

The theoretical foundation and application of two univariate failure detection algorithms to Space Shuttle Main Engine (SSME) test firing data is presented. Both algorithms were applied to data collected during steady state operation of the engine. One algorithm, the time series algorithm, is based on time series techniques and involves the computation of autoregressive models. Time series techniques have been previously applied to SSME data. The second algorithm is based on standard signal processing techniques. It consists of tracking the variations in the average signal power with time. The average signal power algorithm is a newly proposed SSME failure detection algorithm. Seven nominal test firings were used to develop failure indication thresholds for each algorithm. These thresholds were tested using four anomalous firings and one additional nominal firing. Both algorithms provided significantly earlier failure indication times than did the current redline limit system. Neither algorithm gave false failure indications for the nominal firing. The strengths and weaknesses of the two algorithms are discussed and compared. The average signal algorithm was found to have several advantages over the time series algorithm.

Meyer, Claudia M.↗

Task scheduling in dataflow computer architectures

Dataflow computers provide a platform for the solution of a large class of computational problems, which includes digital signal processing and image processing. Many typical applications are represented by a set of tasks which can be repetitively executed in parallel as specified by an associated dataflow graph. Research in this area aims to model these architectures, develop scheduling procedures, and predict the transient and steady state performance. Researchers at NASA have created a model and developed associated software tools which are capable of analyzing a dataflow graph and predicting its runtime performance under various resource and timing constraints. These models and tools were extended and used in this work. Experiments using these tools revealed certain properties of such graphs that require further study. Specifically, the transient behavior at the beginning of the execution of a graph can have a significant effect on the steady state performance. Transformation and retiming of the application algorithm and its initial conditions can produce a different transient behavior and consequently different steady state performance. The effect of such transformations on the resource requirements or under resource constraints requires extensive study. Task scheduling to obtain maximum performance (based on user-defined criteria), or to satisfy a set of resource constraints, can also be significantly affected by a transformation of the application algorithm. Since task scheduling is performed by heuristic algorithms, further research is needed to determine if new scheduling heuristics can be developed that can exploit such transformations. This work has provided the initial development for further long-term research efforts. A simulation tool was completed to provide insight into the transient and steady state execution of a dataflow graph. A set of scheduling algorithms was completed which can operate in conjunction with the modeling and performance tools previously developed. Initial studies on the performance of these algorithms were done to examine the effects of application algorithm transformations as measured by such quantities as number of processors, time between outputs, time between input and output, communication time, and memory size.

Katsinis, Constantine↗

Numerical Studies of Boundary-Layer Receptivity

Direct numerical simulations (DNS) of the acoustic receptivity process on a semi-infinite flat plate with a modified-super-elliptic (MSE) leading edge are performed. The incompressible Navier-Stokes equations are solved in stream-function/vorticity form in a general curvilinear coordinate system. The steady basic-state solution is found by solving the governing equations using an alternating direction implicit (ADI) procedure which takes advantage of the parallelism present in line-splitting techniques. Time-harmonic oscillations of the farfield velocity are applied as unsteady boundary conditions to the unsteady disturbance equations. An efficient time-harmonic scheme is used to produce the disturbance solutions. Buffer-zone techniques have been applied to eliminate wave reflection from the outflow boundary. The spatial evolution of Tollmien-Schlichting (T-S) waves is analyzed and compared with experiment and theory. The effects of nose-radius, frequency, Reynolds number, angle of attack, and amplitude of the acoustic wave are investigated. This work is being performed in conjunction with the experiments at the Arizona State University Unsteady Wind Tunnel under the direction of Professor William Saric. The simulations are of the same configuration and parameters used in the wind-tunnel experiments.

Reed, Helen L.↗

"Coulombic Viscosity" In Granular Materials: Planetary and Astrophysical Implications

The term "Coulombic viscosity" is introduced here to define an empirically observed phenomenon from experiments conducted in both microgravity, and in ground-based 1-g conditions. In the latter case, a sand attrition device was employed to test the longevity of aeolian materials by creating two intersecting grain-circulation paths or cells that would lead to most of the grain energy being expended on grain-to-grain collisions (simulating dune systems). In the areas in the device where gravitationally-driven grain-slurries recycled the sand, the slurries moved with a boundary-layer impeded motion down the chamber walls. Excessive electrostatic charging of the grains during these experiments was prevented by the use of an a.c. corona (created by a Tesla coil) through which the grains passed on every cycle. This created both positive and negative ions which neutralized the triboelectrically-generated grain charges. When the corona was switched on, the velocity of the wall-attached slurries increased by a factor of two as approximately determined by direct observation. What appeared to be a freely-flowing slurry of grains impeded only by intergranular mechanical friction, had obviously been significantly retarded in its motion by electrostatic forces between the grains; with the charging reduced, the grains were able to move past one another without a flow "viscosity" imposed by the Coulombic intergranular forces. A similar phenomenon was observed during microgravity experiments aboard Space Shuttle in USML-1 & USML-2 spacelabs where freely-suspended clouds of sand were being investigated for their potential to for-m aggregates. In this environment, the grains were also charged electrostatically (by natural processes prior to flight), but were free from the intervention of gravity in their interactions. The grains were dispersed into dense clouds by bursts of air turbulence and allowed to form aggregates as the ballistic and turbulent motions damped out. During this very brief (30-60 sec) damping period, motion of the grains was observed to be retarded by the electrostatic interactions. The fact that the grains almost instantly formed aggregates was evidence that their ballistic motions had been constrained and redirected by the dipole-dipole interactions that led to filamentary aggregate development. Undoubtedly, the "Coulombic viscosity" of the cloud assisted in damping grain motion so rapidly. The electrostatically-induced grain-cloud viscosity or drag exerted on grain motion, is a complex function of three major parameters: charge magnitude, charge sign, and mean intergranular distance. The above experiments illustrate one particular type of granular behavior. The discussion here will therefore be restricted to drag relationships: (a) between grains that are naturally charged triboelectrically and thus exhibit dipole-dipole attractions between one another even if there are slight net charges present (which can be overwhelmed by dipole coupling at short distances), and (b) between grains that are densely spaced where the intergranular distance varies between zero and some value (usually tens or hundreds of grain diameters) that permits each grain to detect the dipole moment of another grain -- the distance is not so great that other grains appears as neutral electrical "singularities. I. Aeolian transport: During motion of grains in a saltation cloud (on Earth, Mars, or Venus), triboelectric charging must occur as a result of multiple grain contacts, and by friction with the entraining air. A situation might develop that is similar to the one described above in the attrition device: grain motion becoming significantly retarded (reduced flux) as grains find it increasingly difficult to either separate from the surface, or to pass one another without Coulombic retarding forces. A "Coulombic drag" will exist at flux initiation and increase with time to work in direct opposition to the aerodynamic drag that drives the grain motion. It is predicted that this will lead to an increase with time of both the aerodynamic and bed-dilatancy thresholds (3). Because of Paschen discharge effects in the martian atmosphere, the electrostatic charging in a saltation cloud may be partially abated, but this will lead to greater grain mobility, more charging, and thus to a charge-discharge steady state mediated by mechanical interactions. II. Dry colluvial systems: Sand avalanches on dunes, dry debris flows, talus flows, avalanches, and pyroclastic surges are examples of gravity-driven, dense granular flows where rock/grain fragmentation and grain-to-grain interactions cause triboelectrification (sometimes augmented by other electrical charging processes), and where the grain densities of the systems are such that strong dipole-dipole interactions between grains might be expected to be present. Because it is expected that the Coulombic forces between grains will cause a sluggishness or enhanced granular-flow viscosity, the motion of a grain mass will be retarded or damped so that this will assist, ultimately, in terminating the flow. The greatest Coulombic viscosity will be created in the most highly charged systems, which will also be the most energetic. Thus, grain flows have some tendency to be self-limiting by internal energy partitioning; gravitational potential is converted to Coulombic potential, which manifests itself as a drag force between the grains. III. Volcanic eruption plumes and impact ejecta curtains: The violence of these systems leads to powerful electrical charging of particulates. Lightning storms emanating from volcanic plumes are a testimony to the levels of charging. As pyroclastic grains interact forcefully and frequently within eruption plumes, it is reasonable to predict that the internal turbulent motions of the plume will be significantly damped by the Coulombic viscosity exerted by grain charges. Additional information is contained in the original.

Marshall, J. R.↗