One-dimensional space-charge theory. II - Relativistic Child law.
Space-charge theory for steady one-dimensional current flow across evaluated space between plane electrodes
SEARCH · Search NASA
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.
Space-charge theory for steady one-dimensional current flow across evaluated space between plane electrodes
Ion currents may be several times Child law limit if initial kinetic energies of one-dimensional bipolar space-charge current are of order of potential energy
Recent studies have applied variational calculus, conformal mapping, and point transformations to generalize the one-dimensional (1D) space-charge limited current density (SCLCD) and electron emission mechanisms to nonplanar geometries; however, these assessments have focused on extending the Child–Langmuir law (CLL) for SCLCD in vacuum. Since the charge in the diode is independent of the coordinate system (i.e., covariant), we apply bijective point transformations to extend the Mott–Gurney law (MGL) for the SCLCD in a collisional or semiconductor gap to nonplanar 1D geometries. This yields a modified MGL that replaces the Cartesian gap distance with a canonical gap distance that may be written generally in terms of geometric scale factors that are known for multiple geometries. We tabulate results for common geometries. Such an approach may be applied to any current density, including non-space-charge limited gaps and SCLCD that may fall between the CLL and MGL.
We identify the time-dependent physics responsible for the critical reduction of current losses in magnetically insulated transmission lines (MITLs) due to uninsulated space charge-limited currents of electrons emitted by field stress. A drive current of sufficiently short pulse length introduces a strong enough time dependence that steady-state results alone become inadequate for the complete understanding of current losses. The time-dependent physics can be described as a nonlocal, diamagnetic electromagnetic response of space charge limited currents. As the pulse length is increased or equivalently, the MITL length reduced, these time-dependent effects diminish and current losses converge to those predicted by the well-known Child–Langmuir law in the external (vacuum) fields. We present a simple one-dimensional (1D) model that encapsulates the essence of this physics. We find excellent agreement with 2D particle-in-cell simulations for two MITL geometries, Cartesian parallel plate and azimuthally symmetric straight coaxial. Based on the 1D model, we explore various scaling dependencies of MITL losses with relevant parameters, e.g., peak current, pulse length, geometrical dimensions, etc. We propose an improved physics model of magnetic insulation in the form of a Hull curve, which could also help improve predictions of current losses by common circuit element codes, such as BERTHA. Finally, we describe how to calculate the temperature rise due to electron impact within the 1D model.
The space-charge-limited current (SCLC) in a vacuum diode is given by the Child-Langmuir law (CLL), whose electric potential ϕ(x) ∝ (x/D) 4/3 , where x is the spatial coordinate across the gap and D is the gap separation distance. For a collisional diode, SCLC is given by the Mott-Gurney law (MGL) and ϕ(x) ∝ (x/D) 3/2 . Here, we apply a capacitance argument for SCLC and use the transit time from a recent exact solution for collisional SCLC to show that ϕ(x) ∝ (x/D)ξ for a general collisional gap, where 4/3 ≤ ξ ≤ 3/2 . Furthermore, ξ is strictly a function of νT, where ν is the collision frequency and T is the electron transit time. Using this definition of ξ, we estimate the spatial dependence of the electron velocity and use the gap capacitance to derive an analytic equation for collisional SCLC that agrees within ~4.5% of the exact solution that requires solving parametrically through T. This analytic equation for general ξ asymptotically recovers the CLL as ν → 0 and the MGL as ν → ∞. As a result, matching these limits shows that ξ ≈ 1.40 and V ∝ D 2 ν 2 at the transition from a vacuum to a collisional diode for any device condition.
Oppositely directed flows of positive and negative charges constitute an ambipolar current. A general case is analyzed in which the charges traverse the evacuated space between plane-parallel boundaries and all particles of a given charge species are monoenergetic. It is assumed that the charges may possess non-vanishing initial kinetic energy. This represents an extension of analyses by Langmuir and Muller-Lubeck for vanishing initial kinetic energy. Sample calculations of dimensionless current densities, electric and potential fields, and charge-density distributions are exhibited for cases where the electric field is assumed to vanish at one boundary but the species initial velocities may not. It is shown that resulting ion currents may be several times the Child's law limit if initial kinetic energies are of the order of the potential energy.
Space charge theory of one-dimensional, monoenergetic, unipolar and bipolar currents traversing on otherwise evacuated space between infinite plane parallel electrodes
Full scale reactor model based on fluid equations is widely used to analyze high density plasma reactors. It is well known that the submillimeter scale sheath in front of a biased electrode supporting the wafer is difficult to resolve in numerical simulations, and the common practice is to use results for electric field from some form of analytical sheath model as boundary conditions for full scale reactor simulation. There are several sheath models in the literature ranging from Child's law to a recent unified sheath model [P. A. Miller and M. E. Riley, J. Appl. Phys. 82, 3689 (1997)l. In the present work, the cold ion fluid equations in the radio frequency sheath are solved numerically to show that the spatiotemporal variation of ion flux inside the sheath, commonly ignored in analytical models, is important in determining the electric field and ion energy at the electrode. Consequently, a semianalytical model that includes the spatiotemporal variation of ion flux is developed for use as boundary condition in reactor simulations. This semianalytical model is shown to yield results for sheath properties in close agreement with numerical solutions.
A model describing the doubly charged ion densities in argon and xenon ion thrusters is presented. Doubly charged ions are shown to be produced in significant numbers from both the neutral and singly ionized states. Doubly-to-singly charged ion density ratios calculated using this model are compared to experimental values of this ratio measured in the 15 cm multipole thruster using both propellants. Agreement between theory and experiment is shown to be good. Using this doubly charged ion model, and a similar one obtained previously for mercury, together with the constant neutral loss rate theory, Child's current density law and accelerator grid fabrication limitations, a set of scaling relationships are developed. These relationships show the propellant utilizations, thrust densities and discharge power levels that can be expected as thruster diameter and/or specific impulse are increased and doubly charged ion densities are held at acceptably low values.
This memorandum summarizes a study undertaken to apply the program MAFIA to the modeling of an electron gun in a traveling wave tube (TWT). The basic problem is to emit particles from the cathode in the proper manner. The electrons are emitted with the classical Maxwell-Boltzmann (M-B) energy distribution; and for a small patch of emitting surface; the distribution with angle obeys Lambert's law. This states that the current density drops off as the cosine of the angle from the normal. The motivation for the work is to extend the analysis beyond that which has been done using older codes. Some existing programs use the Child-Langmuir, or 3/2 power law, for the description of the gun. This means the current varies as the 3/2 power of the anode voltage. The proportionality constant is termed the perveance of the gun. This is limited, however, since the 3/2 variation is only an approximation. Also, if the cathode is near saturation, the 3/2 law definitely will not hold. In most of the older codes, the electron beam is decomposed into current tubes, which imply laminar flow in the beam; even though experiments show the flow to be turbulent. Also, the proper inclusion of noise in the beam is not possible. These older methods of calculation do, however, give reasonable values for parameters of the electron beam and the overall gun, and these values will be used as the starting point for a more precise particle-in-cell (PIC) calculation. To minimize the time needed for a given computer run, all beams will use the same number of particles in a simulation. This is accomplished by varying the mass and charge of the emitted particles (macroparticles) in a certain manner, to be consistent with the desired beam current.
Green's theorem and Green's formula for the diffusive cosmic-ray transport equation in relativistic flows are derived. Green's formula gives the solution of the transport equation in terms of the Green's function of the adjoint transport equation, and in terms of distributed sources throughout the region R of interest, plus terms involving the particle intensity and streaming on the boundary. The adjoint transport equation describes the time-reversed particle transport. An Euler-Lagrange variational principle is then obtained for both the mean scattering frame distribution function f, and its adjoint f(dagger). Variations of the variational functional with respect to f(dagger) yield the transport equation, whereas variations of f yield the adjoint transport equation. The variational principle, when combined with Noether's theorem, yields the conservation law associated with Green's theorem. An investigation of the transport equation for steady, azimuthal, rotating flows suggests the introduction of a new independent variable H to replace the comoving frame momentum variable p'. For the case of rigid rotating flows, H is conserved and is shown to be analogous to the Hamiltonian for a bead on a rigidly rotating wire. The variable H corresponds to a balance between the centrifugal force and the particle inertia in the rotating frame. The physical interpretation of H includes a discussion of nonrelativistic and special relativistic rotating flows as well as the cases of aziuthal, differentially rotating flows about Schwarzs-child and Kerr black holes. Green's formula is then applied to the problem of the acceleration of ultra-high-energy cosmic rays by galactic rotation. The model for galactic rotation assumes an angular velocity law Omega = Omega(sub 0)(omega(sub 0)/omega), where omega denotes radial distance from the axis of rotation. Green's functions for the galactic rotation problem are used to investigate the spectrum of accelerated particles arising from monoenergetic and truncated power-law sources. We conclude that it is possible to accelerate particles beyond the knee by galactic rotation, but not in sufficient number to adequately explain the observed spectrum.
The law of the wake and the law of the wall in incompressible turbulent boundary layers formulated by Coles (1956) and their use by Mathews et al. (1970) in the development of a wall-wake representation of the velocity profile in a form applicable for isoenergetic compressible boundary layers are extended to a modified wall-wake velocity profile for turbulent compressible boundary layers. The modified wall-wake profile is shown to provide good representations of experimental velocity distributions.
BACKGROUND: New methods of R-R interval variability based on fractal scaling and nonlinear dynamics ("chaos theory") may give new insights into heart rate dynamics. The aims of this study were to (1) systematically characterize and quantify the effects of aging from early childhood to advanced age on 24-hour heart rate dynamics in healthy subjects; (2) compare age-related changes in conventional time- and frequency-domain measures with changes in newly derived measures based on fractal scaling and complexity (chaos) theory; and (3) further test the hypothesis that there is loss of complexity and altered fractal scaling of heart rate dynamics with advanced age. METHODS AND RESULTS: The relationship between age and cardiac interbeat (R-R) interval dynamics from childhood to senescence was studied in 114 healthy subjects (age range, 1 to 82 years) by measurement of the slope, beta, of the power-law regression line (log power-log frequency) of R-R interval variability (10(-4) to 10(-2) Hz), approximate entropy (ApEn), short-term (alpha(1)) and intermediate-term (alpha(2)) fractal scaling exponents obtained by detrended fluctuation analysis, and traditional time- and frequency-domain measures from 24-hour ECG recordings. Compared with young adults (<40 years old, n=29), children (<15 years old, n=27) showed similar complexity (ApEn) and fractal correlation properties (alpha(1), alpha(2), beta) of R-R interval dynamics despite lower spectral and time-domain measures. Progressive loss of complexity (decreased ApEn, r=-0.69, P<0.001) and alterations of long-term fractal-like heart rate behavior (increased alpha(2), r=0.63, decreased beta, r=-0.60, P<0.001 for both) were observed thereafter from middle age (40 to 60 years, n=29) to old age (>60 years, n=29). CONCLUSIONS: Cardiac interbeat interval dynamics change markedly from childhood to old age in healthy subjects. Children show complexity and fractal correlation properties of R-R interval time series comparable to those of young adults, despite lower overall heart rate variability. Healthy aging is associated with R-R interval dynamics showing higher regularity and altered fractal scaling consistent with a loss of complex variability.