Search NASASearch

SEARCH · Search NASA

Results for “VELOCITY”

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 91 records · Page 5

Thermal iron ions in high speed solar wind streams. II - Temperatures and bulk velocities

Mitchel and Roelof (1980) reported the detection of iron in high speed solar wind flows using the small, but finite sensitivity of solid state detectors to Fe ions in the low energy (50-200 keV protons) L1 channel of the NOAA/JHU energetic particle experiment (EPE). In the current investigation, the EPE response is modeled to a convected Maxwellian to obtain the thermal velocity, flow angle, and bulk velocity of the iron distribution. It is assumed that the iron bulk flow velocity can be represented as a vector sum of the hydrogen bulk velocity and an interplanetary magnetic field (IMF) aligned velocity increment. It is found that the velocity increment is smaller than the local Alfven speed in magnitude, and that the iron thermal velocity is comparable with or greater than the proton thermal velocity, with the 'thermal' velocity defined as the square root of 2kT/m.

Mitchell, D. G.

Impact and explosion crater ejecta, fragment size, and velocity

A model was developed for the mass distribution of fragments that are ejected at a given velocity for impact and explosion craters. The model is semi-empirical in nature and is derived from (1) numerical calculations of cratering and the resultant mass versus ejection velocity, (2) observed ejecta blanket particle size distributions, (3) an empirical relationship between maximum ejecta fragment size and crater diameter and an assumption on the functional form for the distribution of fragements ejected at a given velocity. This model implies that for planetary impacts into competent rock, the distribution of fragments ejected at a given velocity are nearly monodisperse, e.g., 20% of the mass of the ejecta at a given velocity contain fragments having a mass less than 0.1 times a mass of the largest fragment moving at that velocity. Using this model, the largest fragment that can be ejected from asteroids, the moon, Mars, and Earth is calculated as a function of crater diameter. In addition, the internal energy of ejecta versus ejecta velocity is found. The internal energy of fragments having velocities exceeding the escape velocity of the moon will exceed the energy required for incipient melting for solid silicates and thus, constrains the maximum ejected solid fragment size.

Okeefe, J. D.

The manifestation of Alfven's hypothesis of critical ionization velocity in the performance of MPD thrusters

The role of Alfven's critical ionization velocity in the performance of the self-field MPD thruster has been investigated. The existence of a well defined characteristic velocity can be attributed to an ionization process involving the production of a population of suprathermal electrons by an electrostatic instability. It is shown that for the MPD thruster plasma, suprathermalization of electrons via this electrostatic instability can only happen if ions are initially accelerated to velocities larger than the Alfven critical ionization velocity. When this occurs the mechanism will be initiated and the ions decelerated to velocities near the critical velocity. This mechanism ceases to be limiting when all neutrals are ionized. A model of MPD thruster terminal behavior, incorporating Alfven's hypothesis, is presented. Experiments with three different propellants reveal that operation at values of the current squared to total mass flow ratio corresponding to the Alfven critical velocity is marked by a transition wherein low frequency voltage oscillations and a notable change in the voltage-current dependence occurs. One major result of this study is the demonstration that the Alfven critical velocity is not a fundamental limitation on MPD exhaust velocity.

Choueiri, E. Y.

Computing induced velocity perturbations due to a helicopter fuselage in a free stream

The velocity field of a representative helicopter fuselage in a free stream is computed. Perturbation velocities due to the fuselage are computed in a plan above the location of the helicopter rotor (rotor removed). The velocity perturbations computed by a source-panel model of the fuselage are compared with experimental measurements taken with a laser velocimeter. Three paneled fuselage models are studied: fuselage shape, fuselage shape with hub shape, and a body of revolution. The velocity perturbations computed for both fuselage shape models agree well with the measured velocity field except in the close vicinity of the rotor hub. In the hub region, without knowing the extent of separation, modeling of the effective source shape is difficult. The effects of the fuselage perturbations are not well-predicted with a simplified ellipsoid fuselage. The velocity perturbations due to the fuselage at the plane of the measurements have magnitudes of less than 8 percent of free-stream velocity. The velocity perturbations computed by the panel method are tabulated for the same locations at which previously reported rotor-inflow velocity measurements were made.

Berry, John D.

Sea surface velocities from visible and infrared multispectral atmospheric mapping sensor imagery

High resolution (100 m), sequential Multispectral Atmospheric Mapping Sensor (MAMS) images were used in a study to calculate advective surface velocities using the Maximum Cross Correlation (MCC) technique. Radiance and brightness temperature gradient magnitude images were formed from visible (0.48 microns) and infrared (11.12 microns) image pairs, respectively, of Chandeleur Sound, which is a shallow body of water northeast of the Mississippi delta, at 145546 GMT and 170701 GMT on 30 Mar. 1989. The gradient magnitude images enhanced the surface water feature boundaries, and a lower cutoff on the gradient magnitudes calculated allowed the undesirable sunglare and backscatter gradients in the visible images, and the water vapor absorption gradients in the infrared images, to be reduced in strength. Requiring high (greater than 0.4) maximum cross correlation coefficients and spatial coherence of the vector field aided in the selection of an optimal template size of 10 x 10 pixels (first image) and search limit of 20 pixels (second image) to use in the MCC technique. Use of these optimum input parameters to the MCC algorithm, and high correlation and spatial coherence filtering of the resulting velocity field from the MCC calculation yielded a clustered velocity distribution over the visible and infrared gradient images. The velocity field calculated from the visible gradient image pair agreed well with a subjective analysis of the motion, but the velocity field from the infrared gradient image pair did not. This was attributed to the changing shapes of the gradient features, their nonuniqueness, and large displacements relative to the mean distance between them. These problems implied a lower repeat time for the imagery was needed in order to improve the velocity field derived from gradient imagery. Suggestions are given for optimizing the repeat time of sequential imagery when using the MCC method for motion studies. Applying the MCC method to the infrared brightness temperature imagery yielded a velocity field which did agree with the subjective analysis of the motion and that derived from the visible gradient imagery. Differences between the visible and infrared derived velocities were 14.9 cm/s in speed and 56.7 degrees in direction. Both of these velocity fields also agreed well with the motion expected from considerations of the ocean bottom topography and wind and tidal forcing in the study area during the 2.175 hour time interval.

Pope, P. A.

Burning Velocity Measurements in Aluminum-Air Suspensions using Bunsen Type Dust Flames

Laminar burning velocity (sometimes also referred in literature as fundamental or normal flame propagation speed) is probably the most important combustion characteristic of the premixed combustible mixture. The majority of experimental data on burning velocities in gaseous mixtures was obtained with the help of the Bunsen conical flame. The Bunsen cone method was found to be sufficiently accurate for gaseous mixtures with burning velocities higher than 10-15 cm/s at normal pressure. Hans Cassel was the first to demonstrate that suspensions of micron-size solid fuel particles in a gaseous oxidizer can also form self-sustained Bunsen flames. He was able to stabilize Bunsen flames in a number of suspensions of different nonvolatile solid fuels (aluminum, carbon, and boron). Using the Bunsen cone method he estimated burning velocities in the premixed aluminum-air mixtures (particle size less than 10 microns) to be in the range of 30-40 cm/s. Cassel also found, that the burning velocity in dust clouds is a function of the burner diameter. In our recent work, we have used the Bunsen cone method to investigate dependence of burning velocity on dust concentration in fuel-rich aluminum dust clouds. Burning velocities in stoichiometric and fuel-rich aluminum dust suspensions with average particle sizes of about 5 microns were found to be in the range of 20-25 cm/s and largely independent on dust concentration. These results raise the question to what degree burning velocities derived from Bunsen flame specifically and other dust flame configurations in general, are indeed fundamental characteristics of the mixture and to what degree are they apparatus dependent. Dust flames in comparison to gas combustion, are thicker, may be influenced by radiation heat transfer in the flame front, respond differently to heat losses, and are fundamentally influenced by the particular flow configuration due to the particles inertia. Since characteristic spatial scales of dust flames are larger, one can expect that they will also be more sensitive than homogeneous combustion to a particular experimental geometric configuration of the flame and the flow. With such sensitivity the introduction of the very concept of the fundamental flame speed may be problematic for dust combustion. With this in mind, the objective of the present work is to further investigate Bunsen dust flames and evaluate to what degree burning velocities derived from Bunsen cone depend on experimental conditions (i.e. flow rate and nozzle diameter).

Lee, John

Concentration and Velocity Gradients in Fluidized Beds

In this work we focus on the height dependence of particle concentration, average velocity components, fluctuations in these velocities and, with the flow turned off, the sedimentation velocity. The latter quantities are measured using Particle Imaging Velocimetry (PIV). The PIV technique uses a 1-megapixel camera to capture two time-displaced images of particles in the bed. The depth of field of the imaging system is approximately 0.5 cm. The camera images a region with characteristic length of 2.6 cm for the small particles and 4.7 cm. for the large particles. The local direction of particle flow is determined by calculating the correlation function for sub-regions of 32 x 32 pixels. The velocity vector map is created from this correlation function using the time between images (we use 15 to 30 ms). The software is sensitive variations of 1/64th of a pixel. We produce velocity maps at various heights, each consisting of 3844 velocities. We break this map into three vertical zones for increased height information. The concentration profile is measured using an expanded (1 cm diameter) linearly polarized HeNe Laser incident on the fluidized bed. A COHU camera (gamma=1, AGC off) with a lens and a polarizer images the transmitted linearly polarized light to minimize the effects of multiply scattered light. The intensity profile (640 X 480 pixels) is well described by a Gaussian fit and the height of the Gaussian is used to characterize the concentration. This value is compared to the heights found for known concentrations. The sedimentation velocity is estimated using by imaging a region near the bottom of the bed and using PIV to measure the velocity as a function of time. With a nearly uniform concentration profile, the time can be converted to height information. The stable fluidized beds are made from large pseudo-monodisperse particles (silica spheres with radii (250-300) microns and (425-500) microns) dispersed in a glycerin/water mix. The Peclet number is sufficiently large that Brownian motion of the particles can be ignored and the Reynolds number sufficiently small that particle inertia is negligible. A packed particle bed is used to randomize and disperse the flowing fluid introduced by a peristaltic pump. The bed itself is a rectangular glass cell 8 cm wide (x), 0.8 cm deep and a height of 30.5 cm (z). The depth of field of the camera is approximately 0.5 cm so depth information is averaged. Over flow fluid is returned to the reservoir making a closed loop system. In these experiments the particles form a sediment approximately 5.7 cm high with the pump off and expand to 22 cm with the pump on. For the smaller particles the pump velocity is .5 millimeters per second and 1.1 millimeters per second for the large particles. At this concentration the bed has a very well defined top where particle concentration rapidly drops to zero.

McClymer, James P.

Ultraviolet Molecular Rayleigh Scattering Used to Measure Velocity in High-Speed Flow

Molecular Rayleigh scattering offers a means to measure gas flow parameters including density, temperature, and velocity. No seeding of the flow is necessary. The Rayleigh scattered power is proportional to the gas density, the spectral width is related to the gas temperature, and the shift in the frequency of the spectral peak is proportional to one component of the fluid velocity. Velocity measurements based on Rayleigh scattering are more suitable for high-speed flow, where the bulk fluid velocity is on the order of, or larger than, the molecular thermal velocities. Use of ultraviolet wavelengths for Rayleigh scattering diagnostics is attractive for two reasons. First, the Rayleigh scattering cross section is proportional to the inverse 4th power of the wavelength. And second, the reflectivity of metallic surfaces is generally less than it is at longer wavelengths. This is of particular interest in confined flow situations, such as in small wind tunnels and aircraft engine components, where the stray laser light scattered from the windows and internal surfaces in the test facility limits the application of Rayleigh scattering diagnostics. In this work at the NASA Lewis Research Center, molecular Rayleigh scattering of the 266-nm fourth harmonic of a pulsed, injection seeded Nd:YAG (neodymium:yttriumaluminum- garnet) laser was used to measure velocity in a supersonic free air jet with a 9.3- mm exit diameter. The frequency of the Rayleigh scattered light was analyzed with a planar mirror Fabry-Perot interferometer used in a static imaging mode, with the images recorded on a cooled, high-quantum-efficiency charge-coupled discharge (CCD) camera. In addition, some unshifted light from the same laser pulse was imaged through the interferometer to generate a reference. Data were obtained with single laser pulses at velocities up to Mach 1.3. The measured velocities were in good agreement with velocities calculated from isentropic flow relations. Our conclusion from this study was that ultraviolet Rayleigh scattering is preferable in confined flow situations because of the increase in the ratio of Rayleigh scattering signal to stray laser light. On the other hand, in open flows, such as free jets and larger wind tunnels where stray laser light can be controlled, visible Rayleigh scattering is preferable.

Seasholtz, Richard G.

Stochastic modal velocity field in rough-wall turbulence

Stochastically generated instantaneous velocity profiles are used to reproduce the outer region of rough-wall turbulent boundary layers in a range of Reynolds numbers extending from the wind tunnel to field conditions. Each profile consists in a sequence of steps, defined by the modal velocities and representing uniform momentum zones (UMZs), separated by velocity jumps representing the internal shear layers. Height-dependent UMZ is described by a minimal set of attributes: thickness, mid-height elevation, and streamwise (modal) and vertical velocities. These are informed by experimental observations and reproducing the statistical behaviour of rough-wall turbulence and attached eddy scaling, consistent with the corresponding experimental datasets. Sets of independently generated profiles are reorganized in the streamwise direction to form a spatially consistent modal velocity field, starting from any randomly selected profile. The operation allows one to stretch or compress the velocity field in space, increases the size of the domain and adjusts the size of the largest emerging structures to the Reynolds number of the simulated flow. By imposing the autocorrelation function of the modal velocity field to be anchored on the experimental measurements, we obtain a physically based spatial resolution, which is employed in the computation of the velocity spectrum, and second-order structure functions. The results reproduce the Kolmogorov inertial range extending from the UMZ and their attached-eddy vertical organization to the very-large-scale motions (VLSMs) introduced with the reordering process. The dynamic role of VLSM is confirmed in the –u'w' co-spectra and in their vertical derivative, representing a scale-dependent pressure gradient contribution.

42 ENGINEERING

A three-point velocity estimation method for two-dimensional coarse-grained imaging data

Time delay and velocity estimation methods have been widely studied subjects in the context of signal processing, with applications in many different fields of physics. The velocity of waves or coherent fluctuation structures is commonly estimated as the distance between two measurement points divided by the time lag that maximizes the cross correlation function between the measured signals, but this is demonstrated to result in erroneous estimates for two spatial dimensions. We present an improved method to accurately estimate both components of the velocity vector, relying on three non-aligned measurement points. We introduce a stochastic process describing the fluctuations as a superposition of uncorrelated pulses moving in two dimensions. Using this model, we show that the three-point velocity estimation method, using time delays calculated through cross correlations, yields the exact velocity components when all pulses have the same velocity. The two- and three-point methods are tested on synthetic data generated from realizations of such processes for which the underlying velocity components are known. The results reveal the superiority of the three-point technique. Finally, we demonstrate the applicability of the velocity estimation on gas puff imaging data of strongly intermittent plasma fluctuations due to the radial motion of coherent, blob-like structures at the boundary of the Alcator C-Mod tokamak.

Materials Science

Graph-learning approach to combine multiresolution seismic velocity models

SUMMARY The resolution of velocity models obtained by tomography varies due to multiple factors and variables, such as the inversion approach, ray coverage, data quality, etc. Combining velocity models with different resolutions can enable more accurate ground motion simulations. Toward this goal, we present a novel methodology to fuse multiresolution seismic velocity maps with probabilistic graphical models (PGMs). The PGMs provide segmentation results, corresponding to various velocity intervals, in seismic velocity models with different resolutions. Further, by considering physical information (such as ray path density), we introduce physics-informed probabilistic graphical models (PIPGMs). These models provide data-driven relations between subdomains with low (LR) and high (HR) resolutions. Transferring (segmented) distribution information from the HR regions enhances the details in the LR regions by solving a maximum likelihood problem with prior knowledge from HR models. When updating areas bordering HR and LR regions, a patch-scanning policy is adopted to consider local patterns and avoid sharp boundaries. To evaluate the efficacy of the proposed PGM fusion method, we tested the fusion approach on both a synthetic checkerboard model and a fault zone structure imaged from the 2019 Ridgecrest, CA, earthquake sequence. The Ridgecrest fault zone image consists of a shallow (top 1 km) high-resolution shear-wave velocity model obtained from ambient noise tomography, which is embedded into the coarser Statewide California Earthquake Center Community Velocity Model version S4.26-M01. The model efficacy is underscored by the deviation between observed and calculated traveltimes along the boundaries between HR and LR regions, 38 per cent less than obtained by conventional Gaussian interpolation. The proposed PGM fusion method can merge any gridded multiresolution velocity model, a valuable tool for computational seismology and ground motion estimation.

Geochemistry & Geophysics

Symmetry Determining Equations of the Rankine-Hugoniot Equations for Variable Velocity Shock Waves

The “constant velocity piston” problem (Fig. 1), also known as the “piston problem,” is a standard model for a one dimensional, in our case linear, symmetric shock wave moving through an inviscid, perfect gas. The model can be divided into two regions - a perturbed section on the left and an unperturbed section on the right - by a moving shock wave moving left to right. Both the perturbed and unperturbed sections, i.e. the shocked and unshocked regions, respectively, obey the Eulerian conservation equations; however, at the exact location of the shock, there is a mathematical discontinuity not satisfied by the Euler equations. To ensure continuity and conservation of certain quantities when crossing between the unshocked and shocked regions, we evoke a series of equations derived from the Eulerian conservation equations, called the Rankine-Hugoniot equations, or “jump” equations as it is often referred to in the literature on the topic. The classical constant-velocity piston problem assumes the piston features a constant driving velocity (among many other willing suspensions of belief required in the pursuit of a first principles equation model); consequent to this assumption is a constant-velocity shock and a constant-velocity shocked flow state. However, using Lie Group Theory (LGT), also known as symmetry analysis, we can attempt to reinterpret the model with a shock wave of variable velocity in time and space. An extension of the model in this way opens up the possibility for obtaining new analytical solutions to the piston problem for certain shock velocity models. In this report, we use LGT to derive the symmetry determining equations (SDEs), whose solutions are Lie groups, which permit analytical solutions. In the future, we can then use the SDEs to define constraint equations on the shock velocity model and what the successive solutions to the Euler equations might be based off such constraints. This report is structured as follows: Section 2 provides a brief derivation of the Rankine-Hugoniot (“jump”) equations; Section 3 gives an overview of Lie group theory; Section 4 derives the SDEs of the jump equations; Section 5 derives the Euler conservation equations for fluids; and Section 6 presents concluding remarks and opportunities for future studies.

42 ENGINEERING

Derivation and evaluation of an approximate analysis for three-dimensional viscous subsonic flow with large secondary velocities

An approximate analysis is presented for calculating three-dimensional, low Mach number, laminar viscous flows in curved passages with large secondary flows and corner boundary layers. The analysis is based on the decomposition of the overall velocity field into inviscid and viscous components with the overall velocity being determined from superposition. An incompressible vorticity transport equation is used to estimate inviscid secondary flow velocities to be used as corrections to the potential flow velocity field. A parabolized streamwise momentum equation coupled to an adiabatic energy equation and global continuity equation is used to obtain an approximate viscous correction to the pressure and longitudinal velocity fields. A collateral flow assumption is invoked to estimate the viscous correction to the transverse velocity fields. The approximate analysis is solved numerically using an implicit ADI solution for the viscous pressure and velocity fields. An iterative ADI procedure is used to solve for the inviscid secondary vorticity and velocity fields. This method was applied to computing the flow within a turbine vane passage with inlet flow conditions of M = 0.1 and M = 0.25, Re = 1000 and adiabatic walls, and for a constant radius curved rectangular duct with R/D = 12 and 14 and with inlet flow conditions of M = 0.1, Re = 1000, and adiabatic walls.

Anderson, O. L.

The size distributions of fragments ejected at a given velocity from impact craters

The mass distribution of fragments that are ejected at a given velocity for impact craters is modeled to allow extrapolation of laboratory, field, and numerical results to large scale planetary events. The model is semi-empirical in nature and is derived from: (1) numerical calculations of cratering and the resultant mass versus ejection velocity, (2) observed ejecta blanket particle size distributions, (3) an empirical relationship between maximum ejecta fragment size and crater diameter, (4) measurements and theory of maximum ejecta size versus ejecta velocity, and (5) an assumption on the functional form for the distribution of fragments ejected at a given velocity. This model implies that or planetary impacts into competent rock, the distribution of fragments ejected at a given velocity is broad, e.g., 68% of the mass of the ejecta at a given velocity contains fragments having a mass less than 0.1 times a mass of the largest fragment moving at that velocity. The broad distribution suggests that in impact processes, additional comminution of ejecta occurs after the upward initial shock has passed in the process of the ejecta velocity vector rotating from an initially downward orientation. This additional comminution produces the broader size distribution in impact ejecta as compared to that obtained in simple brittle failure experiments.

Okeefe, J. D.

A new instrument to measure charged and neutral cometary dust particles at low and high impact velocities

A new class of dust particle detector, the PVDF dust detector, was designed for space missions such as the Halley Comet missions where the particle impact velocity is very high. It is demonstrated that this same PVDF detector (operating in a different mode) also has the capability of detecting dust particles having low velocity (approx. 100 m/s). This low velocity detection capability is extremely important in terms of planned missions requiring measurement of low velocity dust particles such as comet rendezvous missions. An additional detecting element (charge induction cylinder) was also developed which, when combined with a PVDF detector, yields a system which will measure the charge (magnitude and sign) carried by a cometary particle as well as the particle velocity and mass for impact velocities in the range 100 to 500 m/s. Since the cylinder-PVDF detector system has a relatively small geometry factors, an array of PVDF detectors was included having a total sensing area of 0.1 sq m for measurements in regions of space where the dust flux is expected to be low. The characteristics of the detectors in this array have been chosen to provide optimum mass sensitivity for both low-velocity cometary dust as well as high-velocity asteroid associated and interplanetary dust.

Economon, T.

Spatially growing disturbances in a high velocity ratio two-stream, coplanar jet

The influence of cold and heated secondary flow on the instability of a two-stream, coplanar jet having a 0.7 Mach number heated primary jet for a nominal fan to primary velocity ratio of 0.68 was investigated by means of inviscid linearized stability theory. The instability properties of spatially growing axisymmetric and first order azimuthal disturbances were studied. The instability characteristics of the two-stream jet with a velocity ratio of 0.68 are very different from those of a single stream jet, and a two-stream, coplanar jet having a 0.9 Mach number heated primary jet and a cold secondary jet for a fan to primary velocity ratio of 0.30. For X/D = 1 and in comparison to the case where the velocity ratio was 0.3, the presence of the fan stream with a velocity ratio of 0.68 enhanced the instability of the jet and increased the unstable frequency range. However, the axisymmetric mode (m = 0) and the first order azimuthal mode (m = 1) have similar spatial growth rates where the velocity ratio is 0.68 while for a velocity ratio of 0.3 the growth rate of the first order azimuthal mode (m = 1) is greater. Comparing the cold and hot secondary flow results showed that for a velocity ratio of 0.68 the growth rate is greater for cold.

Miles, J. H.

Spatially growing disturbances in a high velocity ratio two-stream, coplanar jet

The influence of cold and heated secondary flow on the instability of a two-stream, coplanar jet having a 0.7 Mach number heated primary jet for a nominal fan to primary velocity ratio of 0.68 was investigated by means of inviscid linearized stability theory. The instability properties of spatially growing axisymmetric and first order azimuthal disturbances were studied. The instability characteristics of the two-stream jet with a velocity ratio of 0.68 are very different from those of a single stream jet, and a two-stream, coplanar jet having a 0.9 Mach number heated primary jet and a cold secondary jet for a fan to primary velocity ratio of 0.30. For X/D = 1 and in comparison to the case where the velocity ratio was 0.3, the presence of the fan stream with a velocity ratio of 0.68 enhanced the instability of the jet and increased the unstable frequency range. However, the axisymmetric mode (m = 0) and the first order azimuthal mode (m = 1) have similar spatial growth rates where the velocity ratio is 0.68 while for a velocity ratio of 0.3 the growth rate of the first order azimuthal mode (m = 1) is greater. Comparing the cold and hot secondary flow results showed that for a velocity ratio of 0.68 the growth rate is greater for cold.

Miles, Jeffrey H.

On the turbulent heating and the threshold condition in the critical ionization velocity interaction

On the basis of the nonlinear treatment of the ion beam instability and the modified two-stream instability, threshold conditions for the critical ionization velocity interaction are derived. There are three different regimes of interaction: (1) additional ionization for relative velocities smaller than the critical velocity, (2) the self-sustained discharge for velocities greater than the critical velocity which indeed turns out as a sharp threshold, and (3) an explosive growth regime for velocities exceeding 1.5 times the critical velocity. Additional charge exchange collisions of ions and energy loss of electrons due to excitation do not change the basic threshold behavior but modify the value of the critical velocity.

Moebius, E.