Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Maxwell equations”

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 307 records · Page 17

An Implicit Characteristic Based Method for Electromagnetics

An implicit characteristic-based approach for numerical solution of Maxwell's time-dependent curl equations in flux conservative form is introduced. This method combines a characteristic based finite difference spatial approximation with an implicit lower-upper approximate factorization (LU/AF) time integration scheme. This approach is advantageous for three-dimensional applications because the characteristic differencing enables a two-factor approximate factorization that retains its unconditional stability in three space dimensions, and it does not require solution of tridiagonal systems. Results are given both for a Fourier analysis of stability, damping and dispersion properties, and for one-dimensional model problems involving propagation and scattering for free space and dielectric materials using both uniform and nonuniform grids. The explicit Finite Difference Time Domain Method (FDTD) algorithm is used as a convenient reference algorithm for comparison. The one-dimensional results indicate that for low frequency problems on a highly resolved uniform or nonuniform grid, this LU/AF algorithm can produce accurate solutions at Courant numbers significantly greater than one, with a corresponding improvement in efficiency for simulating a given period of time. This approach appears promising for development of dispersion optimized LU/AF schemes for three dimensional applications.

Beggs, John H.↗

Expendable Launch Vehicle Studies

Analytical support studies of expendable launch vehicles concentrates on the stability of the dynamics during launch especially during or near the region of maximum dynamic pressure. The in-plane linearized dynamic equations of a generic launch vehicle with multiple flexible bending and fuel sloshing modes are developed. The design of a robust LQR controller based on the reduced order system is accomplished using the parameter perturbation technique. The ELV modeling and analysis team has spent the past three years working on the theoretical development and application of sensitivity analysis to solve eigenvalue problems associated with structural dynamics. Specific application areas include stochastic vibrations, viscously damped vibrations and nonlinear dynamics. In stochastic linear vibrations, sensitivity analysis methods have been developed which determined the expectation and variance of the response for given probability density functions of various stochastic parameters. Examples worked include beam vibrations with up to six stochastic parameters; and the eigensensitivity results differ little from those obtained, with much effect, from Monte Carlo techniques. For viscously damped vibrations, eigensensitivity analysis has given excellent results for homogeneous beams, modeled by either quadratic or quartic eigenvalue equations. Specific applications include Kevin and Maxwell-type viscoelastic beams.

Bainum, P. M.↗

A case study of the Thunderstorm Research International Project storm of July 11, 1978. I - Analysis of the data base

The data obtained at the time of the Thunderstorm Research International Project storm at the Kennedy Space Center on July 11, 1978 are analyzed in a model-independent manner. The data base included data from three Doppler radars, a lightning detection and ranging system and a network of 25 electric field mills, and rain gages. Electric field measurements were used to analyze the charge moments transferred by lightning flashes, and the data were fitted to Weibull distributions; these were used to estimate statistical parameters of the lightning for both intracloud and cloud-to-ground flashes and to estimate the fraction of the flashes which were below the observation threshold. The displacement and the conduction current densities were calculated from electric field measurements between flashes. These values were used to derive the magnitudes and the locations of dipole and monopole generators by least squares fitting the measured Maxwell current densities to the displacement-dominated equations.

Nisbet, John S.↗

Temperature-dependent creep buckling of plates

Time-dependent lateral deflection of flat rectangular plates is predicted by the Norton-Bailey (Norton 1929, Bailey 1935) power law for material creep. The plates have a through-thickness steady-state temperature distribution, and the effects are considered by using Maxwell's law to modify the power creep law. Equations are derived for creep exponents of 3 and 5, using the sandwich plate element to predict creep buckling of plates. Predictions of creep buckling with a temperature variation between the inner and outer plate surfaces are found to be somewhat dependent on the creep buckling relationship assumed. When significant scatter justifies a variation in the creep constants up to an order of magnitude, discrepancies in predictions using the two exponents are reasonable, and for one engineering material, the predictions have the same degree of agreement with experimental data as have the respective creep laws.

Ross, D. A.↗

Testing continuum descriptions of low-Mach-number shock structures

Numerical experiments have been performed on normal shock waves with Monte Carlo Direct Simulations (MCDS's) to investigate the validity of continuum theories at very low Mach numbers. Results from the Navier-Stokes and the Burnett equations are compared to MCDS's for both hard-sphere and Maxwell gases. It is found that the maximum-slope shock thicknesses are described equally well (within the MCDS computational scatter) by either of the continuum formulations for Mach numbers smaller than about 1.2. For Mach numbers greater that 1.2, the Burnett predictions are more accurate than the Navier-Stokes results. Temperature-density profile separations are best described by the Burnett equations for Mach numbers greater than about 1.3. At lower Mach numbers the MCDS scatter is too great to differentiate between the two continuum theories. For all Mach numbers above one, the shock shapes are more accurately described by the Burnett equations.

Pham-Van-diep, Gerald C.↗

The effects of molecular weight on the single lap shear creep and constant strain rate behavior of thermoplastic polyimidesulfone adhesive

The bonded shear creep and constant strain rate behaviors of zero, one, and three percent endcapped thermoplastic polyimidesulfone adhesive were examined at room and elevated temperatures. Endcapping was accomplished by the addition of phthalic anhydrides. The primary objective was to determine the effects of molecular weight on the mechanical properties of the adhesive. Viscoelastic and nonlinear elastic constitutive equations were utilized to model the adhesive. Ludwik's and Crochet's relations were used to describe the experimental failure data. The effects of molecular weight changes on the above mentioned mechanical behavior were assessed. The viscoelastic Chase-Goldsmith and elastic nonlinear relations gave a good fit to the experimental stress strain behavior. Crochet's relations based on Maxwell and Chase-Goldsmith models were fit to delayed failure data. Ludwik's equations revealed negligible rate dependence. Ultimate stress levels and the safe levels for creep stresses were found to decrease as molecular weight was reduced.

Dembosky, Stanley K.↗

Gas Flux and Density Surrounding a Cylindrical Aperture in the Free Molecular Flow Regime

The equations for rigorously calculating the particle flux and density surrounding a cylindrical aperture in the free molecular flow regime are developed and presented. The fundamental equations for particle flux and density from a reservoir and a diffusely reflecting surface will initially be developed. Assumptions will include a Maxwell-Boltzmann speed distribution, equal particle and wall temperatures, and a linear flux distribution along the cylindrical aperture walls. With this information, the equations for axial flux and density surrounding a cylindrical aperture will be developed. The cylindrical aperture will be divided into multiple volumes and regions to rigorously determine the surrounding axial flux and density, and appropriate limits of integration will be determined. The results of these equations will then be evaluated. The linear wall flux distribution assumption will be assessed. The axial flux and density surrounding a cylindrical aperture with a thickness-to-radius ratio of 1.25 will be presented. Finally, the equations determined in this study will be verified using multiple methods.

Soulas, George C.↗

Formation of the 0.511 MeV line in solar flares

The slowing down and annihilation of positrons and the formation of positronium in a solar flare plasma are investigated to determine how the width of the 0.511-MeV line and its strength relative to the three-photon continuum from positronium decay depend on the temperature and density of the medium in which the positron comes to rest. The calculations are limited to the cases of annihilation in a completely ionized plasma, in a partially ionized plasma with an electron/neutron density ratio of 1.0 or 0.1, and in an atomic gas with a very small ion density. Thermally averaged rate coefficients are obtained for the free annihilation of positrons and for positronium formation through radiative recombination in the fully ionized plasma. Positronium formation rates and the resultant energy distributions of the positronium atoms in the partially ionized medium are determined by numerically solving the Fokker-Planck equation in a medium where the ambient free electrons have a Maxwell-Boltzmann distribution of finite temperature but the density of the medium is sufficiently low that positronium atoms decay without further collisions following their formation. A Monte Carlo calculation is performed for the positron energy loss, positronium formation through charge exchange, and positronium breakup in the weakly ionized medium. The energy distributions of decaying positronium atoms and the relative number of triplet to singlet positronium decays are evaluated for ion concentrations and densities characteristic of the solar photosphere.

Crannell, C. J.↗

Convection zone origins of solar atmospheric heating

Spicules are examined as a means for supplying the corona with mass, energy, and magnetic field. It is suggested that spicules form from the supersonic upward expansion of material on nearly evacuated network flux tubes embedded within the sun's convection zone. This allows supersonic but subescape velocities to be attained by the material as it flows outward through the photosphere. Although supersonic, the kinetic energy (subescape) of the spicule material, as observed, is insufficient for coronal heating. It is suggested that, through buoyancy changes on evacuated flux tubes, the magnetic field first 'wicks' material flow into the solar atmosphere. Subsequently, the magnetic field energizes the gaseous material to form the conventional hot, dynamically expanding, solar corona. This occurs through momentum and energy transport by Alfven waves and associated Maxwell stresses concurrently flowing upward on these 'geysers' (spicules). The vertical momentum equation governing fluid flow is examined, and a particular equipartition solution is presented for the flow velocity along a simple field geometry.

Schatten, Kenneth H.↗

Recovery of the Navier-Stokes equations using a lattice-gas Boltzmann method

A lattice Boltzmann model is presented which gives the complete Navier-Stokes equation and may provide an efficient parallel numerical method for solving various fluid problems. The model uses the single-time relaxation approximation and a particular Maxwell-type distribution. The model eliminates exactly (1) the non-Galilean invariance caused by a density-dependent coefficient in the convection term and (2) a velocity-dependent equation of state.

Chen, Hudong↗

Investigation of mass transfer between two parallel walls at different temperatures by a moment method

One-dimensional flow between two fixed parallel walls composed of the same substance but at different temperatures and spaced a distance 1 apart is considered. The hot plate is the evaporating surface (source) and the cold plate is the condensing surface (sink). The vapor between the two plates is assumed to be a monatomic gas consisting of Maxwell molecules. Lee's moment method is used to obtain a set of six nonlinear equations. Both the nonlinear equations and a linearized approximation to them are solved.

Sloat, T. N.↗

A Transport Model for Non-Local Heating of Electrons in ICP Reactors

A new model has been developed for non-local heating of electrons in ICP reactors, based on a hydrodynamic approach. The model has been derived using the electron momentum conservation in azimuthal direction with electromagnetic and frictional forces respectively as driving force and damper of harmonic oscillatory motion of electrons. The resulting transport equations include the convection of azimuthal electron momentum in radial and axial directions, thereby accounting for the non-local effects. The azimuthal velocity of electrons and the resulting electrical current are coupled to the Maxwell's relations, thus forming a self-consistent model for non-local heating. This model is being implemented along with a set of Navier-Stokes equations for plasma dynamics and gas flow to simulate low-pressure (few mTorr's) ICP discharges. Characteristics of nitrogen plasma in a TCP 300mm etch reactor is being studied. The results will be compared against the available Langmuir probe measurements.

Chang, C. H.↗

Thundercloud current determination from measurements at the earth's surface

For the analysis of the actual dominant mechanism or mechanisms of electrification in a given storm cell, it is of great importance to make observations which can be used to determine the configuration and the magnitudes of the currents and impedances of the generator systems in active thunderclouds. Krider and Musser (1982) have measured Maxwell currents below Florida thunderclouds and shown that their patterns are consistent with simple generator geometries. It is, therefore, of great interest to study the extent to which a set of Maxwell current measurements can be inverted to obtain the generator currents and locations. Nisbet (1983) has developed an electrical model for active thunderclouds which solves equations for current densities and electric fields. The present investigation is concerned with the extent to which such a model can be used to derive data about the cloud generator magnitudes and locations, taking into account also the conductivity profile from ground-based measurements or estimates of the Maxwell current density.

Nisbet, J. S.↗

Finite grid radius and thickness effects on retarding potential analyzer measured suprathermal electron density and temperature

The effect of finite grid radius and thickness on the electron current measured by planar retarding potential analyzers (RPAs) is analyzed numerically. Depending on the plasma environment, the current is significantly reduced below that which is calculated using a theoretical equation derived for an idealized RPA having grids with infinite radius and vanishingly small thickness. A correction factor to the idealized theoretical equation is derived for the Pioneer Venus (PV) orbiter RPA (ORPA) for electron gasses consisting of one or more components obeying Maxwell statistics. The error in density and temperature of Maxwellian electron distributions previously derived from ORPA data using the theoretical expression for the idealized ORPA is evaluated by comparing the densities and temperatures derived from a sample of PV ORPA data using the theoretical expression with and without the correction factor.

Knudsen, William C.↗