Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Maxwell”

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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

Effect of Coulomb interaction on time of flight of cold antiprotons launched from an ion trap

Time-of-flight spectra for Maxwell-Boltzman (MB) distributions of antiprotons initially held in an ion trap and detected after being launched through a 50-cm-long shielding drift tube have been calculated. The distributions used are of temperature 0.4-40 K, cubic length 0.003-3.0 cm, and number 10-100 particles. The mutual Coulomb repulsion of the particles causes a reduction in the number of late arrival particles expected from the MB velocity distribution. The Coulomb energy is not equally divided among the particles during the expansion. The energy is transferred preferentially to the outer particles so that the reduction in the number of slow particles is not necessarily large. The reduction factor is found to be greater than unity when the potential energy of the trapped ions is greater than about 5 percent of the ions' kinetic energy and is about 2 for the launch parameters of the Los Alamos antiproton gravity experiment.

Camp, J. B.↗
NASA NTRS1993

Calculations and experiments concerning lifting force and power in TEMPUS

A critical comparison is reported between the theoretically predicted and experimentally measured values for the electromagnetic lifting force and the heating rates which may be achieved, under simulated microgravity conditions, using the TEMPUS electromagnetic levitation device. The experiments involved the suspending of a metallic sample from one arm of a recording balance, such that it was carefully positioned between the heating and the positioning coils of the levitation device. The net force exerted by the sample was measured as a function of position, the coil currents, and the nature of the sample. Some calculations are also reported regarding the power absorption by the sample. The theoretical predictions, based on the numerical solution of Maxwell's equations using the volume integral technique, were found to be in excellent agreement with the measurements. For the idealized case of a spherical sample, analytical solutions describing the lifting force were also found to agree very well with the computed results.

Zong, J. H.↗
NASA NTRS1993

Superposition of nonlinear plasma waves

We report results showing that spatially periodic Bernstein-Greene-Kruskal (BGK) waves, which are exact nonlinear traveling wave solutions of the Vlasov-Maxwell equations for collisionless plasmas, satisfy a nonlinear principle of superposition in the small amplitude limit. The analysis explicates the notion of superimposed BGK waves which, as recent numerical calculations suggest, is crucial in the proper description of the time-asymptotic state of a plasma when a large amplitude electrostatic wave undergoes nonlinear Landau damping.

Buchanan, Mark↗
NASA NTRS1992

Parallel decomposition methods for the solution of electromagnetic scattering problems

This paper contains a overview of the methods used in decomposing solutions to scattering problems onto coarse-grained parallel processors. Initially, a short summary of relevant computer architecture is presented as background to the subsequent discussion. After the introduction of a programming model for problem decomposition, specific decompositions of finite difference time domain, finite element, and integral equation solutions to Maxwell's equations are presented. The paper concludes with an outline of possible software-assisted decomposition methods and a summary.

Cwik, Tom↗
NASA NTRS1993

Flow resolution and domain of influence in rarefied hypersonic blunt-body flows

The study assesses the effects of the upstream domain size and grid resolution upon flow properties and body aerodynamics computed for rarefied flows over cold blunt bodies with a direct simulation Monte Carlo (DSMC) particle method. Empirical correlations are suggested for aerodynamic coefficients for two-dimensional flows past a perpendicular flat plate. Free-stream parameters which were varied in the study include the Mach number, Knudsen number, surface temperature, and intermolecular potential. Insufficient grid resolution leads to overprediction of aerodynamic heating and forces in the DSMC method. Solution accuracy correlates well with the Reynolds number defined at the wall temperature and the stagnation mean free path relative to the cell dimension. Insufficient upstream domain size in the DSMC method leads to overprediction of heating and drag. Errors in aerodynamic coefficients correlate well with the distance ahead of the body where flow temperature reaches half of its peak value. Simulation of a hard-sphere gas is more sensitive to grid resolution, while simulation of a Maxwell gas is more sensitive to upstream domain size.

Haas, Brian L.↗
NASA NTRS1993

Shape optimization for aerodynamic efficiency and low observability

Field methods based on the finite-difference approximations of the time-domain Maxwell's equations and the potential-flow equation have been developed to solve the multidisciplinary problem of airfoil shaping for aerodynamic efficiency and low radar cross section (RCS). A parametric study and an optimization study employing the two analysis methods are presented to illustrate their combined capabilities. The parametric study shows that for frontal radar illumination, the RCS of an airfoil is independent of the chordwise location of maximum thickness but depends strongly on the maximum thickness, leading-edge radius, and leadingedge shape. In addition, this study shows that the RCS of an airfoil can be reduced without significant effects on its transonic aerodynamic efficiency by reducing the leading-edge radius and/or modifying the shape of the leading edge. The optimization study involves the minimization of wave drag for a non-lifting, symmetrical airfoil with constraints on the airfoil maximum thickness and monostatic RCS. This optimization study shows that the two analysis methods can be used effectively to design aerodynamically efficient airfoils with certain desired RCS characteristics.

Vinh, Hoang↗
NASA NTRS1993

Analysis of microwave-heated rocket engines for space propulsion

Microwave-heated plasmas in convergent nozzles are analyzed using a coupled Maxwell and Navier-Stokes solver to examine relevant issues associated with microwave thermal propulsion. Parametric studies are conducted to understand the effect of power, pressure, and plasma location with respect to the nozzle throat. For nozzles in the 0.5 to 3 N range with helium flow, results show that specific impulses up to 550-650 sec are possible, with further increases being limited by severe wall-heating. Coupling efficiencies of over 90 percent are consistently obtained, with overall efficiencies ranging from 40 percent to 80 percent. Size scale-up studies (done by scaling the frequency from 2.45 GHz to 0.91 GHz) indicate that plasma migration toward the walls occurs more frequently for the lower frequency. Increasing the cavity aspect ratio and detuning the cavity are found to be effective ways of keeping the plasma on-axis.

Schwer, D. A.↗
NASA NTRS1993

Interaction with the lower ionosphere of electromagnetic pulses from lightning - Heating, attachment, and ionization

A Boltzmann formulation of the electron distribution function and Maxwell's equations for the EM fields are used to simulate the interaction of lightning radiated EM pulses with the lower ionosphere. Ionization and dissociative attachment induced by the heated electrons cause significant changes in the local electron density, N(e). Due to 'slow' field changes of typical lightning EM pulses over time scales of tens of microsec, the distribution function follows the quasi-equilibrium solution of the Boltzmann equation in the altitude range of interest (70 to 100 km). The EM pulse is simulated as a planar 100 microsec long single period oscillation of a 10 kHz wave injected at 70 km. Under nighttime conditions, individual pulses of intensity 10-20 V/m (normalized to 100 km horizontal distance) produce changes in N(e) of 1-30 percent while a sequence of pulses leads to strong modification of N(e) at altitudes less than 95 km. The N(e) changes produce a 'sharpening' of the lower ionospheric boundary by causing a reduction in electron density at 75-85 km (due to attachment) and a substantial increase at 85-95 km (due to ionization) (e.g., the scale height decreases by a factor of about 2 at about 85 km for a single 20 V/m EM pulse). No substantial N(e) changes occur during daytime.

Taranenko, Y. N.↗
NASA NTRS1993

Galileo/NIMS near-infrared thermal imagery of the surface of Venus

Numerous highland and lowland features on the surface of Venus are observed in multispectral imagery acquired at approximately 50 km spatial resolution by the Near-Infrared Mapping Spectrometer (NIMS) on board the Galileo spacecraft in Feb. 1990. Specifically, such features are observed at 1.18 microns, a wavelength particularly sensitive to thermal emission from the hot, lower atmosphere (less than 10 km) and surface, and show up particularly well when this image is 'de-clouded' using a simultaneously-acquired 2.3-microns image of the upper, cloudy atmosphere. Due to the steep atmospheric temperature gradient (approximately 8 degrees per kilometer), hot lowland areas appear relatively bright, while cooler, highland areas appear dark (due to the steep atmospheric temperature gradient - approximately 8 degrees per kilometer - surface temperatures span approximately 100 K over the 13 kilometer range of surface altitudes observed in this image). Prominent highland features include Maxwell Montes (approximately 12 km altitude), Alpha Regio (2.5 km), Eistla Regio (approximately 2.0 km), Bell Regio (2-3 km), and the western edge of Aphrodite Terra (2-2.5 km). Low-lying regions include Sedna Planitia (-1.0 km), Tinatin Planitia (-0.5 km), and the Bereghinya Planitia (0 km). From correlations with radar altimetry maps, such imagery may place useful constraints on surface emissivity and temperature variations, as well as on the nature of continuum opacity of CO2 in the 1-micron region.

Carlson, R. W.↗
NASA NTRS1993

Measurement of birefringence for optical recording disk substrates

The birefringence of bare and coated substrates for magneto-optical recording is experimentally investigated using ellipsometry at the wavelengths of 632.8 nm and 780 nm. The rotation and ellipticity of the polarization state of the reflected or transmitted light is measured for different incident angles and different orientations of the incident linear polarization. The measured data is then fitted by a computer program which solves the Maxwell equations for the plane-wave propagation in a multilayer structure and minimizes the error between the measured and calculated data by adjusting the unknown parameters of the multilayer. This approach enables us to determine orientations of the three principal axes in the substrate and the corresponding refractive indices. A special feature of our ellipsometers is that a glass hemisphere is placed in contact with the substrate, which eliminates the refraction of the incident light and enables a maximum propagation angle of 70 degrees (with respect to the normal) in the substrate. This increases the sensitivity of the measurement. Certain anomalies were observed, which we believe are associated with the presence of grooves on these substrates.

Fu, Hong↗
NASA NTRS1993

Aspects of modelling the tectonics of large volcanoes on the terrestrial planets

Analytic solutions for the responses of planetary lithospheres to volcanic loads have been used to model faulting and infer elastic plate thicknesses. Predictions of the distribution of faulting around volcanic loads, based on the application of Anderson's criteria for faulting to the results of the models, do not agree well with observations. Such models do not give the stress state in the load itself, but only suggest a state of horizontal compressive stress there. Further, these models have considered only the effect of an instantaneously emplaced load. They do not address the time evolution of stresses, nor do they consider the effect of a load which grows. A finite element approach allows us to assign elements to the load itself, and thus permits calculation of the stress state and stress history within the edifice. The effects of episodic load growth can also be treated. When these effects are included, models give much better agreement with observations. We use the finite element code TECTON to construct axisymmetric models of volcanoes resting on an elastic lithospheric plate overlying a viscoelastic asthenosphere. We have implemented time-dependent material properties in order to simulate incremental volcano growth. The viscoelastic layer was taken to extend to a sufficient depth so that a rigid lower boundary has no significant influence on the results. The code first calculates elastic deformations and stresses and then determines the time-dependent viscous deformations and stresses. Time in the model scales as the Maxwell time tau(m) in the asthenosphere. We consider a volcano 25 km in height and 200 km in radius on an elastic lithosphere 40 km thick (parameters approximately appropriate to Ascraeus Mons). The volcano consists of three load increments applied at intervals of 1000 tau(m). Contours of maximum deviatoric stress in the fully-grown edifice at the conclusion of flexure (t = 3000 tau(m)) are shown.

Mcgovern, Patrick J.↗
NASA NTRS1993

Relief of some small landforms on Venus

Three sets of radar images have been acquired under different viewing conditions by the Magellan synthetic aperture radar: (1) left-looking with varied incidence angles (cycle 1); (2) right-looking with nearly constant incidence angles (cycle 2); and (3) left-looking with varied incidence angles, most of which were smaller than those in (1) except for those acquired on passes across Maxwell Montes with incidence angles larger than those in (1) (cycle 3). Image displacements in the radar images that are caused by the relief of landforms provide several methods of estimating this relief: (1) monoscopic measurements of foreshortening of landforms that are symmetrical in the plane of the look-direction of the radar (includes radial symmetry); (2) stereoscopic measurements of parallax in same-side image pairs (cycles 1-2 and 3); and (3) measurements of parallax in opposite-side image pairs (cycles 1-2 and/or 2-3). Success in methods 2 and 3 (especially 3) depends on identifying conjugate image points in the two images. Here, we report our preliminary results for five impact craters, seven small volcanic edifices, and two lava flows. The three methods mentioned above lead to the interesting result that Venusian impact craters have depth-diameter ratios like those on Mars rather than those on Earth, but some appear partly filled. Our results for de Lalande and Melba also suggest filling, but there may be other causes for their relatively small depth-diameter ratios. A host of small volcanic edifices have relief that can be crudely estimated using the above methods. Relief/diameter ratios for our cratered cones are about the same as those of Icelandic lava shields; some Venusian cones resemble the Martian shields of Mareotis-Tempe and Ceraunius Fossae, but the Venusian relief diameter ratios are larger. The smallest cratered dome is similar in size and profile to a Martian dome north of Uranius Patera; the smallest cratered cone resembles one in Chryse Planitia. Lava flows on Venus that are thick enough to measure are rare, but we have applied methods 1 and 3 to the huge flow of Ovda Regio and flows of an unusual volcano, Mahuea Tholus.

Moore, H. J.↗
NASA NTRS1993

Finite difference time domain analysis of chirped dielectric gratings

The finite difference time domain (FDTD) method for solving Maxwell's time-dependent curl equations is accurate, computationally efficient, and straight-forward to implement. Since both time and space derivatives are employed, the propagation of an electromagnetic wave can be treated as an initial-value problem. Second-order central-difference approximations are applied to the space and time derivatives of the electric and magnetic fields providing a discretization of the fields in a volume of space, for a period of time. The solution to this system of equations is stepped through time, thus, simulating the propagation of the incident wave. If the simulation is continued until a steady-state is reached, an appropriate far-field transformation can be applied to the time-domain scattered fields to obtain reflected and transmitted powers. From this information diffraction efficiencies can also be determined. In analyzing the chirped structure, a mesh is applied only to the area immediately around the grating. The size of the mesh is then proportional to the electric size of the grating. Doing this, however, imposes an artificial boundary around the area of interest. An absorbing boundary condition must be applied along the artificial boundary so that the outgoing waves are absorbed as if the boundary were absent. Many such boundary conditions have been developed that give near-perfect absorption. In this analysis, the Mur absorbing boundary conditions are employed. Several grating structures were analyzed using the FDTD method.

Hochmuth, Diane H.↗
NASA NTRS1993

Mathmatical modeling for diffractive optics

We consider a 'diffractive optic' to be a biperiodic surface separating two half-spaces, each having constant constitutive parameters; within a unit cell of the periodic surface and across the transition zone between the two half-spaces, the constitutive parameters can be a continuous, complex-valued function. Mathematical models for diffractive optics have been developed, and implemented as numerical codes, both for the 'direct' problem and for the 'inverse' problem. In problems of the 'direct' class, the diffractive optic is specified, and the full set of Maxwell's equations is cast in a variational form and solved numerically by a finite element approach. This approach is well-posed in the sense that existence and uniqueness of the solution can be proved and specific convergence conditions can be derived. An example of a metallic grating at a Wood anomaly is presented as a case where other approaches are known to have convergence problems. In problems of the 'inverse' class, some information about the diffracted field (e.g., the far-field intensity) is given, and the problem is to find the periodic structure in some optimal sense. Two approaches are described: phase reconstruction in the far-field approximation; and relaxed optimal design based on the Helmholtz equation. Practical examples are discussed for each approach to the inverse problem, including array generators in the far-field case and antireflective structures for the relaxed optimal design.

Dobson, David↗
NASA NTRS1993

Analysis of microstrip dipoles and slots transversely coupled to a microstrip line using the FDTD method

Printed dipole elements and their complement, linear slots, are elementary radiators that have found use in low-profile antenna arrays. Low-profile antenna arrays, in addition to their small size and low weight characteristics, offer the potential advantage of low-cost, high-volume production with easy integration with active integrated circuit components. The design of such arrays requires that the radiation and impedance characteristics of the radiating elements be known. The FDTD (Finite-Difference Time-Domain) method is a general, straight-forward implementation of Maxwell's equations and offers a relatively simple way of analyzing both printed dipole and slot elements. Investigated in this work is the application of the FDTD method to the analysis of printed dipole and slot elements transversely coupled to an infinite transmission line in a multilayered configuration. Such dipole and slot elements may be used in dipole and slot series-fed-type linear arrays, where element offsets and interelement line lengths are used to obtain the desired amplitude distribution and beam direction, respectively. The design of such arrays is achieved using transmission line theory with equivalent circuit models for the radiating elements. In an equivalent circuit model, the dipole represents a shunt impedance to the transmission line, where the impedance is a function of dipole offset, length, and width. Similarly, the slot represents a series impedance to the transmission line. The FDTD method is applied to single dipole and slot elements transversely coupled to an infinite microstrip line using a fixed rectangular grid with Mur's second order absorbing boundary conditions. Frequency-dependent circuit and scattering parameters are obtained by saving desired time-domain quantities and using the Fourier transform. A Gaussian pulse excitation is applied to the microstrip transmission line, where the resulting reflected signal due to the presence of the radiating element is used to determine the equivalent element impedance.

Tulintseff, A. N.↗
NASA NTRS1993

Analytical evaluation of a surface integral expressing the coupling between interior and exterior volumes in a FE-IE approach

Recently an approach which combines the finite element technique and an integral equation to determine the fields scattered by inhomogeneous bodies of complicated shape has been proposed. Basically, a mathematical surface which encloses the scatterers is introduced, thus dividing the space into an interior and an exterior volume, in which the finite element technique and an integral equation for EM scattering, respectively, are applied. The integral equation is set up for the tangential components of the fields at the surface, while the interior volume the unknowns are the total fields. Continuity of the tangential fields at the boundary, as required by Maxwell's equations, is imposed, thus coupling the two methods to obtain a consistent solution. The coupling term is expressed by a surface integral formed by the dot product of a FE basis function and an IE testing function, or viceversa. By choosing the boundary to be a surface of revolution and by making a convenient selection of IE basis (testing) functions, it is possible to evaluate the integrals analytically on surfaces such as curved triangles, curved quadrilaterals and curved pentagons. We will illustrate the salient steps involved in setting up and carrying out these integrals and discuss what class of basis (testing) functions and analytic surfaces of revolution they are applicable to. Analytic calculations offer the advantage of better accuracy than purely numerical ones, and, when combined with them, often shed light on issues of numerical convergence and limiting values. Furthermore, they may reduce computation time and storage requirements.

Zuffada, C.↗
NASA NTRS1993

Time-domain analysis of planar microstrip devices using a generalized Yee-algorithm based on unstructured grids

The generalized Yee-algorithm is presented for the temporal full-wave analysis of planar microstrip devices. This algorithm has the significant advantage over the traditional Yee-algorithm in that it is based on unstructured and irregular grids. The robustness of the generalized Yee-algorithm is that structures that contain curved conductors or complex three-dimensional geometries can be more accurately, and much more conveniently modeled using standard automatic grid generation techniques. This generalized Yee-algorithm is based on the the time-marching solution of the discrete form of Maxwell's equations in their integral form. To this end, the electric and magnetic fields are discretized over a dual, irregular, and unstructured grid. The primary grid is assumed to be composed of general fitted polyhedra distributed throughout the volume. The secondary grid (or dual grid) is built up of the closed polyhedra whose edges connect the centroid's of adjacent primary cells, penetrating shared faces. Faraday's law and Ampere's law are used to update the fields normal to the primary and secondary grid faces, respectively. Subsequently, a correction scheme is introduced to project the normal fields onto the grid edges. It is shown that this scheme is stable, maintains second-order accuracy, and preserves the divergenceless nature of the flux densities. Finally, for computational efficiency the algorithm is structured as a series of sparse matrix-vector multiplications. Based on this scheme, the generalized Yee-algorithm has been implemented on vector and parallel high performance computers in a highly efficient manner.

Gedney, Stephen D.↗
NASA NTRS1993

Extinction efficiencies from DDA calculations solved for finite circular cylinders and disks

One of the most commonly noted uncertainties with respect to the modeling of cirrus clouds and their effect upon the planetary radiation balance is the disputed validity of the use of Mie scattering results as an approximation to the scattering results of the hexagonal plates and columns found in cirrus clouds. This approximation has historically been a kind of default, a result of the lack of an appropriate analytical solution of Maxwell's equations to particles other than infinite cylinders and spheroids. Recently, however, the use of such approximate techniques as the Discrete Dipole Approximation has made scattering solutions on such particles a computationally intensive but feasible possibility. In this study, the Discrete Dipole Approximation (DDA) developed by Flatau (1992) is used to find such solutions for homogeneous, circular cylinders and disks. This can serve to not only assess the validity of the current radiative transfer schemes which are available for the study of cirrus but also to extend the current approximation of equivalent spheres to an approximation of second order, homogeneous finite circular cylinders and disks. The results will be presented in the form of a single variable, the extinction efficiency.

Withrow, J. R.↗