SEARCH · Engineering Papers
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.
Approximate field values computed rapidly. MODEL computer program developed to calculate electromagnetic field values of large loop antenna at all distances to observation point. Antenna assumed to be in x-y plane with center at origin of coordinate system. Calculates field values in both rectangular and spherical components. Also solves for wave impedance. Written in MicroSoft FORTRAN 77.
A method is presented for calculating the electromagnetic wave field on the earth's surface associated with the operation of an electrodynamic tethered satellite system of constant or slowly varying current in the upper ionosphere. The wave field at the ionospheric boundary and on the earth's surface is obtained by numerical integration. The results suggest that the ionospheric waves do not propagate into the atmosphere and that the image of the Alfven wings from a steady-current tether should be greatly broadened on the earth's surface.
Numerical evaluation of Kirchhoff-Kottler formulation of electromagnetic field backscatter and application to antenna surfaces
Some results obtained with a digital computer program written at Goddard Space Flight Center to obtain electromagnetic fields scattered by perfectly reflecting surfaces are presented. For purposes of illustration a paraboloidal reflector was illuminated at radio frequencies in the simulation for both receiving and transmitting modes of operation. Fields were computed in the Fresnel and Fraunhofer regions. A dual-reflector system (Cassegrain) was also simulated for the transmitting case, and fields were computed in the Fraunhofer region. Appended results include derivations which show that the vector Kirchhoff-Kottler formulation has an equivalent form requiring only incident magnetic fields as a driving function. Satisfaction of the radiation conditions at infinity by the equivalent form is demonstrated by a conversion from Cartesian to spherical vector operators. A subsequent development presents the formulation by which Fresnel or Fraunhofer patterns are obtainable for dual-reflector systems. A discussion of the time-average Poynting vector is also appended.
A fully 3-D numerical model is developed to represent magnetic damping of complex fluid flow, heat transfer and electromagnetic field distributions in a melt cavity. The model is developed based on our in-house finite element code for the fluid flow, heat transfer and electromagnetic field calculations. The computer code has been tested against benchmark test problems that are solved by other commercial codes as well as analytical solutions whenever available. The numerical model is tested against numerical and experimental results for water reported in literature. With the model so tested, various numerical simulations are carried out for the Sn-35.5% Pb melt convection and temperature distribution in a cylindrical cavity with and without the presence of a transverse magnetic field. Numerical results show that magnetic damping can be effectively applied to reduce turbulence and flow levels in the melt undergoing solidification and over a certain threshold value a higher magnetic field resulted in a higher velocity reduction. It is found also that for a fully 3-D representation of the magnetic damping effects, the electric field induced in the melt by the applied DC magnetic field does not vanish, as some researchers suggested, and must be included even for molten metal and semiconductors. Also, for the study of the melt flow instability, a long enough time has to be applied to ensure the final fluid flow recirculation pattern. Moreover, our numerical results suggested that there seems to exist a threshold value of applied magnetic field, above which magnetic damping becomes possible and below which the convection in the melt is actually enhanced. Because of the limited financial resource allocated for the project, we are unable to carry out extensive study on this effect, which should warrant further theoretical and experimental study. In that endeavor, the developed numerical model should be very useful; and the model should serve as a useful tool for exploring necessary design parameters for planning magnetic damping experiments and interpreting the experimental results.
A fully 3-D, transient finite element model is developed to represent the magnetic damping effects on complex fluid flow, heat transfer and electromagnetic field distributions in a Sn- 35.5%Pb melt undergoing unidirectional solidification. The model is developed based on our in- house finite element code for the fluid flow, heat transfer and electromagnetic field calculations. The numerical model is tested against numerical and experimental results for water as reported in literature. Various numerical simulations are carried out for the melt convection and temperature distribution with and without the presence of a transverse magnetic field. Numerical results show that magnetic damping can be effectively applied to stabilize melt flow, reduce turbulence and flow levels in the melt and over a certain threshold value a higher magnetic field resulted in a greater reduction in velocity. Also, for the study of melt flow instability, a long enough running time is needed to ensure the final fluid flow recirculation pattern. Moreover, numerical results suggest that there seems to exist a threshold value of applied magnetic field, above which magnetic damping becomes possible and below which the 0 convection in the melt is actually enhanced.
We suggest a merging procedure for the Particle-Mesh Complex Particle Kinetic (PMCPK) method in case of inter-penetrating flow (multiple plasma beams). We examine the standard particle-in-cell (PIC) and the PMCPK methods in the case of particle acceleration by shock surfing for a wide range of the control numerical parameters. The plasma dynamics is described by a hybrid (particle-ion-fluid-electron) model. Note that one may need a mesh if modeling with the computation of an electromagnetic field. Our calculations use specified, time-independent electromagnetic fields for the shock, rather than self-consistently generated fields. While a particle-mesh method is a well-verified approach, the CPK method seems to be a good approach for multiscale modeling that includes multiple regions with various particle/fluid plasma behavior. However, the CPK method is still in need of a verification for studying the basic plasma phenomena: particle heating and acceleration by collisionless shocks, magnetic field reconnection, beam dynamics, etc.
The technical feasibility is studied of growing single-crystal silicon ribbon in the space environment. Procedures are described for calculating the electromagnetic fields produced in a silicon ribbon by an rf shaping coil. The forces on the ribbon and the degree of shaping to be expected are determined. The expected steady-state temperature distribution in the ribbon is calculated in the one-dimensional approximation. Calculations on simplified models indicate, that lack of flatness of the shaped ribbon and excessive heating of the melt by the eddy currents induced by the shaping fields may pose problems. An analysis of the relative effects of various kinds of forces other than electromagnetic showed that in the space environment capillarity forces would dominate, and that the shape of the melt is thus principally determined by the shape of any solids with which it comes in contact. This suggests that ribbon may be produced simply by drawing between parallel wires. A concept is developed for a process of off-angle growth, in which the ribbon is pulled at an angle to the solidification front. Such a process promises to offer increased growth rate, better homogeneity, and thinner ribbon.
Supraluminous longitudinal plasma wave scattering by localized inhomogeneity in thermal plasma density, calculating far field electromagnetic radiation production
Computer program for calculating scattered electromagnetic field patterns from asymmetrical reflectors
System analyses of cylindrical plasma centrifuges are presented, for which the velocity field and electromagnetic fields are calculated. The effects of different electrode geometrics, induced magnetic fields, Hall-effect, and secondary flows are discussed. It is shown that speeds of 10000 m/sec can be achieved in plasma centrifuges, and that an efficient separation of U238 and U235 in uranium plasmas is feasible. The external boundary-value problem for the deposition of sputtering products is reduced to a Fredholm integral equation, which is solved analytically by means of the method of successive approximations.
It is often desirable to calculate the electromagnetic fields inside and about a complicated system of scattering bodies, as well as in their far-field region. The finite element method (FE) is well suited to solving the interior problem, but the domain has to be limited to a manageable size. At the truncation of the FE mesh one can either impose approximate (absorbing) boundary conditions or set up an integral equation (IE) for the fields scattered from the bodies. The latter approach is preferable since it results in higher accuracy. Hence, the two techniques can be successfully combined by introducing a surface that encloses the scatterers, applying a FE model to the inner volume and setting up an IE for the tangential fields components on the surface. Here the continuity of the tangential fields is used bo obtain a consistent solution. A few coupled FE-IE methods have recently appeared in the literature. The approach presented here has the advantage of using edge-based finite elements, a type of finite elements with degrees of freedom associated with edges of the mesh. Because of their properties, they are better suited than the conventional node based elements to represent electromagnetic fields, particularly when inhomogeneous regions are modeled, since the node based elements impose an unnatural continuity of all field components across boundaries of mesh elements. Additionally, our approach is well suited to handle large size problems and lends itself to code parallelization. We will discuss the salient features that make our approach very efficient from the standpoint of numerical computation, and the fields and RCS of a few objects are illustrated as examples.
The development of improved and more efficient plasma reactors is a costly process for the semiconductor industry. Until five years ago, the Industry made most of its advancements through a trial and error approach. More recently, the role of computational modeling in the design process has increased. Both conventional computational fluid dynamics (CFD) techniques like Navier-Stokes solvers as well as particle simulation methods are used to model plasma reactor flowfields. However, since high-density plasma reactors generally operate at low gas pressures on the order of 1 to 10 mTorr, a particle simulation may be necessary because of the failure of CFD techniques to model rarefaction effects. The direct simulation Monte Carlo method is the most widely accepted and employed particle simulation tool and has previously been used to investigate plasma reactor flowfields. A plasma DSMC code is currently under development at NASA Ames Research Center with its foundation as the object-oriented parallel Cornell DSMC code, MONACO. The present investigation is a follow up of a neutral flow investigation of the effects of process parameters as well as reactor design on etch rate and etch rate uniformity. The previous work concentrated on silicon etch of a chlorine flow in a configuration typical of electron cyclotron resonance (ECR) or helical resonator type reactors. The effects of the plasma on the dissociation chemistry were modeled by making assumptions about the electron temperature and number density. The electrons or ions themselves were not simulated.The present work extends these results by simulating the charged species.The electromagnetic fields are calculated such that power deposition is modeled self-consistently. Electron impact reactions are modeled along with mechanisms for charge exchange. An bipolar diffusion assumption is made whereby electrons remain tied to the ions. However, the velocities of tile electrons are allowed to be modified during collisions and are not confined to a Maxwellian distribution. The interaction between the neutral flow and plasma is examined, and results for etch rate uniformity from the previous research and the present plasma simulations are compared.
We describe a new parallel, non-orthogonal grid, three-dimensional electromagnetic particle-in-cell (EMPIC) code based on a finite-volume formulation. This code uses a logically Cartesian grid of deformable hexahedral cells, a discrete surface integral (DSI) algorithm to calculate the electromagnetic field, and a hybrid logical-physical space algorithm to push particles.
Experimental transient electromagnetic field measurements were performed on a Transall C160 aircraft during in-flight lightning strikes. The data allow a test of the predictive capabilities of a three dimensional time domain finite difference code (ALICE) developed at ONERA in order to investigate lightning-aircraft interactions. Using a transfer function technique in the 3D code, it is shown that a bi-leader attached to an aircraft can be simulated by a linear model, and so the electromagnetic fields can be calculated anywhere on the vehicle. Comparison of experimental and numerical results were made for several lightning strikes. Skin current density and electromagnetic field distributions are discussed in detail.
A method to calculate the electromagnetic scattered fields from a dielectric strap wrapped around convex, conducting structure is presented. A moment method technique is used to find the current excited within the strap by the incident plane wave. Then, Uniform Geometrical Theory of Diffraction (UTD) is used to compute the fields scattered by the strap. Reasonable agreement was obtained between the computed and the measured results. The results found in this study are useful in evaluating straps as a target support structure for scattering measurements.
An improved computational technique for calculating the electromagnetic force field, the power absorption and the deformation of an electromagnetically levitated metal sample is described. The technique is based on the volume integral method, but represents a substantial refinement; the coordinate transformation employed allows the efficient treatment of a broad class of rotationally symmetrical bodies. Computed results are presented to represent the behavior of levitation melted metal samples in a multi-coil, multi-frequency levitation unit to be used in microgravity experiments. The theoretical predictions are compared with both analytical solutions and with the results of previous computational efforts for the spherical samples and the agreement has been very good. The treatment of problems involving deformed surfaces and actually predicting the deformed shape of the specimens breaks new ground and should be the major usefulness of the proposed method.
A mathematical formulation is given and computed results are presented describing the behavior of electromagnetically-levitated metal droplets under the conditions of microgravity. In the formulation the electromagnetic force field is calculated using a modification of the volume integral method and these results are then combined with the FIDAP code to calculate the steady state melt velocities. The specific computational results are presented for the conditions corresponding to the planned IML-2 Space Shuttle experiment, using the TEMPUS device, which has separate 'heating' and 'positioning' coils. While the computed results are necessarily specific to the input conditions, some general conclusions may be drawn from this work. These include the fact that for the planned TEMPUS experiments to positioning coils will produce only a weak melt circulation, while the heating coils are like to produce a mildly turbulent recirculating flow pattern within the samples. The computed results also allow us to assess the effect of sample size, material properties and the applied current on these phenomena.