Engineering Papers⌕ Search

Engineering topics

Shebalin, John V.

Publications and source records attributed to Shebalin, John V..

At least 37 records · Page 2

Global Aspects of Charged Particle Motion in Axially Symmetric Multipole Magnetic Fields

The motion of a single charged particle in the space outside of a compact region of steady currents is investigated. The charged particle is assumed to produce negligible electromagnetic radiation, so that its energy is conserved. The source of the magnetic field is represented as a point multipole. After a general description, attention is focused on magnetic fields with axial symmetry. Lagrangian dynamical theory is utilized to identify constants of the motion as well as the equations of motion themselves. The qualitative method of Stonner is used to examine charged particle motion in axisymmetric multipole fields of all orders. Although the equations of motion generally have no analytical solutions and must be integrated numerically to produce a specific orbit, a topological examination of dynamics is possible, and can be used, d la Stonner, to completely describe the global aspects of the motion of a single charged particle in a space with an axisymmetric multipole magnetic field.

Shebalin, John V.↗

The Statistical Mechanics of Ideal Homogeneous Turbulence

Plasmas, such as those found in the space environment or in plasma confinement devices, are often modeled as electrically conducting fluids. When fluids and plasmas are energetically stirred, regions of highly nonlinear, chaotic behavior known as turbulence arise. Understanding the fundamental nature of turbulence is a long-standing theoretical challenge. The present work describes a statistical theory concerning a certain class of nonlinear, finite dimensional, dynamical models of turbulence. These models arise when the partial differential equations describing incompressible, ideal (i.e., nondissipative) homogeneous fluid and magnetofluid (i.e., plasma) turbulence are Fourier transformed into a very large set of ordinary differential equations. These equations define a divergenceless flow in a high-dimensional phase space, which allows for the existence of a Liouville theorem, guaranteeing a distribution function based on constants of the motion (integral invariants). The novelty of these particular dynamical systems is that there are integral invariants other than the energy, and that some of these invariants behave like pseudoscalars under two of the discrete symmetry transformations of physics, parity, and charge conjugation. In this work the 'rugged invariants' of ideal homogeneous turbulence are shown to be the only significant scalar and pseudoscalar invariants. The discovery that pseudoscalar invariants cause symmetries of the original equations to be dynamically broken and induce a nonergodic structure on the associated phase space is the primary result presented here. Applicability of this result to dissipative turbulence is also discussed.

Shebalin, John V.↗

A Spectral Algorithm for Solving the Relativistic Vlasov-Maxwell Equations

A spectral method algorithm is developed for the numerical solution of the full six-dimensional Vlasov-Maxwell system of equations. Here, the focus is on the electron distribution function, with positive ions providing a constant background. The algorithm consists of a Jacobi polynomial-spherical harmonic formulation in velocity space and a trigonometric formulation in position space. A transform procedure is used to evaluate nonlinear terms. The algorithm is suitable for performing moderate resolution simulations on currently available supercomputers for both scientific and engineering applications.

Shebalin, John V.↗

Absolute Equilibrium Entropy

The entropy associated with absolute equilibrium ensemble theories of ideal, homogeneous, fluid and magneto-fluid turbulence is discussed and the three-dimensional fluid case is examined in detail. A sigma-function is defined, whose minimum value with respect to global parameters is the entropy. A comparison is made between the use of global functions sigma and phase functions H (associated with the development of various H-theorems of ideal turbulence). It is shown that the two approaches are complimentary though conceptually different: H-theorems show that an isolated system tends to equilibrium while sigma-functions allow the demonstration that entropy never decreases when two previously isolated systems are combined. This provides a more complete picture of entropy in the statistical mechanics of ideal fluids.

Shebalin, John V.↗

Kolmogorov Flow in Three Dimensions

A numerical study of the long-time evolution of incompressible Navier-Stokes turbulence forced at a single long-wavelength Fourier mode, i.e., a Kolmogorov flow, has been completed. The boundary conditions are periodic in three dimensions and the forcing is effected by imposing a steady, two-dimensional, sinusoidal shear velocity which is directed along the x-direction and varies along the z-direction. A comparison with experimental data shows agreement with measured cross-correlations of the turbulent velocity components which lie in the mean-flow plane. A statistical analysis reveals that the shear-driven turbulence studied here has significant spectral anisotropy which increases with wave number.

Shebalin, John V.↗

Numerical simulation of three-dimensional self-gravitating flow

The three-dimensional flow of a self-gravitating fluid is numerically simulated using a Fourier pseudospectral method with a logarithmic variable formulation. Two cases with zero total angular momentum are studied in detail, a 323 simulation (Run B). Other than the grid size, the primary difference between the two cases are that Run A modeled atomic hydrogen and had considerably more compressible motion initially than Run B, which modeled molecular hydrogen. The numerical results indicate that gravitational collapse can proceed in a variety of ways. In the Run A, collapse led to an elongated tube-like structure, while in the Run B, collapse led to a flatter, disklike structure.

Shebalin, John V.↗

Broken symmetry in ideal magnetohydrodynamic turbulence

A numerical study of the long-time evolution of a number of cases of inviscid, isotropic, incompressible, three-dimensional fluid, and magneto-fluid turbulence has been completed. The results confirm that ideal magnetohydrodynamic turbulence is non-ergodic if there is no external magnetic field present. This is due essentially to a canonical symmetry being broken in an arbitrary dynamical representation. The broken symmetry manifests itself as a coherent structure, i.e., a non-zero time-averaged part of the turbulent magnetic field. The coherent structure is observed, in one case, to contain about eighteen percent of the total energy.

Shebalin, John V.↗

Pseudospectral simulation of compressible turbulence using logarithmic variables

The direct numerical simulation of dissipative, highly compressible turbulent flow is performed using a pseudospectral Fourier technique. The governing equations are cast in a form where the important physical variables are the fluid velocity and the natural logarithms of the fluid density and temperature. Bulk viscosity is utilized to model polyatomic gases more accurately and to ensure numerical stability in the presence of strong shocks. Numerical examples include three-dimensional supersonic homogeneous turbulence and two-dimensional shock-turbulence interactions.

Shebalin, John V.↗

Electromagnetic potential vectors and the Lagrangian of a charged particle

Maxwell's equations can be shown to imply the existence of two independent three-dimensional potential vectors. A comparison between the potential vectors and the electric and magnetic field vectors, using a spatial Fourier transformation, reveals six independent potential components but only four independent electromagnetic field components for each mode. Although the electromagnetic fields determined by Maxwell's equations give a complete description of all possible classical electromagnetic phenomena, potential vectors contains more information and allow for a description of such quantum mechanical phenomena as the Aharonov-Bohm effect. A new result is that a charged particle Lagrangian written in terms of potential vectors automatically contains a 'spontaneous symmetry breaking' potential.

Shebalin, John V.↗

Homogeneous quantum electrodynamic turbulence

The electromagnetic field equations and Dirac equations for oppositely charged wave functions are numerically time-integrated using a spatial Fourier method. The numerical approach used, a spectral transform technique, is based on a continuum representation of physical space. The coupled classical field equations contain a dimensionless parameter which sets the strength of the nonlinear interaction (as the parameter increases, interaction volume decreases). For a parameter value of unity, highly nonlinear behavior in the time-evolution of an individual wave function, analogous to ideal fluid turbulence, is observed. In the truncated Fourier representation which is numerically implemented here, the quantum turbulence is homogeneous but anisotropic and manifests itself in the nonlinear evolution of equilibrium modal spatial spectra for the probability density of each particle and also for the electromagnetic energy density. The results show that nonlinearly interacting fermionic wave functions quickly approach a multi-mode, dynamic equilibrium state, and that this state can be determined by numerical means.

Shebalin, John V.↗

A Fourier collocation time domain method for numerically solving Maxwell's equations

A new method for solving Maxwell's equations in the time domain for arbitrary values of permittivity, conductivity, and permeability is presented. Spatial derivatives are found by a Fourier transform method and time integration is performed using a second order, semi-implicit procedure. Electric and magnetic fields are collocated on the same grid points, rather than on interleaved points, as in the Finite Difference Time Domain (FDTD) method. Numerical results are presented for the propagation of a 2-D Transverse Electromagnetic (TEM) mode out of a parallel plate waveguide and into a dielectric and conducting medium.

Shebalin, John V.↗

Thermoelectric effects in decaying homogeneous magneto-gas turbulence

In the formulation of compressible MHD (i.e., magneto-gas dynamics), a 'generalized Ohm's law is required. In particular, an electron pressure term and a 'Hall effect' term may appear as non-negligible additions to the Ohm's law that is conventionally used for incompressible MHD. In 'high-beta' (i.e., relatively low magnetic energy) situations, the Hall term may be neglected (at least initially) but, as it turns out, the electron pressure term cannot be neglected. Here, three-dimensional, high-beta, homogeneous, decaying, magneto-gas turbulence is examined with regard to this additional term. Through numerical simulation, it is found that 'thermoelectric effects' are produced that significantly alter the evolution of the magnetic field and electric current strengths.

Shebalin, John V.↗

Thermoelectric effects and magnetic field amplification in magnetogasdynamic turbulence

It will be shown that thermoelectric effects amplify magnetic fields in compressible magnetogasdynamic turbulence (though not nearly as much as occurs across a curved reently bowshock). The importance of this result lies in the recognition that thermoelectric effects (in addition to kinetic effects) provide a real mechanism for the amplification of magnetic field strength (and total energy dissipation through ohmic losses) in a compressible, turbulent plasma.

Shebalin, John V.↗

Aerobrake plasmadynamics - Macroscopic effects

The flow around an aerobraking spacecraft (such as the Aeroassist Flight Experiment reentry vehicle) will contain a region of partially ionized gas, that is, a plasma. It is shown here by numerical simulation that macroscopic plasmadynamic effects (which are not included in standard aerothermodynamic simulations) will have a noticeable effect on the reentry flow field. In particular, there are thermoelectric phenomena which can have a major influence in flow dynamics at the front of an ionizing bowshock. These thermoelectric phenomena arise because of the presence of large density and temperature gradients at the front of a reentry bowshock, and they include strong local magnetic fields, electric currents, and ohmic heating. One important result is the dramatic increase in temperature (over the case where plasma effects are neglected) at a reentry shock front; the implication is that macroscopic plasmadynamic effects can no longer be neglected in simulations of hypersonic reentry flow fields.

Shebalin, John V.↗

Magnetic dynamo activity in mechanically driven compressible magnetohydrodynamic turbulence

Magnetic dynamo activity in a homogeneous, dissipative, polytropic, two-dimensional, turbulent magneto-fluid is simulated numerically. The magneto-fluid is simulated numerically. The magneto-fluid is, in a number of cases, mechanically forced so that energy input balances dissipation, thereby maintaining constant energy. In the presence of a mean magnetic field, a magneto-fluid whose initial turbulent magnetic energy is zero quickly arrives at a state of non-zero turbulent magnetic energy. If the mean magnetic field energy density is small, the turbulent magnetic field can achieve a local energy density more than four hundred times larger; if the mean magnetic field energy density is large, then equipartition between the turbulent magnetic and kinetic energy is achieved. Compared to the presence of a mean magnetic field, compressibility appears to have only a marginal effect in mediating the transfer of turbulent kinetic energy into magnetic energy.

Shebalin, John V.↗

An exact solution to the relativistic equation of motion of a charged particle driven by a linearly polarized electromagnetic wave

An exact analytic solution is found for a basic electromagnetic wave-charged particle interaction by solving the nonlinear equations of motion. The particle position, velocity, and corresponding time are found to be explicit functions of the total phase of the wave. Particle position and velocity are thus implicit functions of time. Applications include describing the motion of a free electron driven by an intense laser beam..

Shebalin, John V.↗