Engineering PapersSearch

SEARCH · Engineering Papers

Results for “Poisson equation”

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 145 records · Page 8

A numerical solution algorithm for prediction of turbulent aerodynamic corner flows

A numerical solution algorithm is established for prediction of subsonic turbulent three-dimensional flows in aerodynamic configuration juncture regions. In concert with a complete three-dimensional exterior potential flow solution, the developed parabolic algorithm yields prediction of the details of the corner region flowfield. Turbulence closure is established using the complete Reynolds stress. Pressure coupling is accomplished using the concepts of complementary and particular solutions to a Poisson equation. Numerical results for three-dimensional turbulent flow in the juncture of two intersecting parabolic arc airfoils are presented.

Baker, A. J.

Protostellar formation in rotating interstellar clouds. I - Numerical methods and tests

Attention is given to numerical methods and tests of a series of gravitational hydrodynamics computer codes constructed in order to numerically follow the dynamic collapse of spherically symmetric (1-D), axisymmetric (2-D), and non-axisymmetric (3-D) isothermal interstellar clouds. A spherical harmonic expansion is used to solve the Poisson equation for the gravitational potential. It is shown that the use of explicit donor-cell hydrodynamics on a moving spherical coordinate grid ensures mass and momentum conservation and allows the grid to follow the collapse of the fluid. Finally, the performance of the codes is examined.

Boss, A. P.

Prediction and measurement of turbulent aerodynamic trailing edge flows

A viscous-inviscid interaction algorithm is developed for prediction of two-dimensional mean and fluctuating velocity distributions in the wake immediately downstream of an airfoil trailing edge. A composite pressure field is defined, and a Poisson equation solved for transverse pressure variations. A parabolized form of the time-averaged steady Navier-Stokes equations are solved in conjunction with a viscous-augmented two-dimensional inviscid potential flow analysis. A tensor constitutive equation is employed to predict Reynolds stress distributions from solutions of a turbulence kinetic energy two equation closure model. Numerical predictions compared favorably with detailed experimental data for mean and fluctuating velocities, and Reynolds shear stress distributions, in the trailing edge region of a NACA 63-012 airfoil.

Baker, A. J.

Generation of orthogonal boundary-fitted coordinate systems

A method is presented for computing orthogonal boundary fitted coordinate systems for geometries with coordinate distributions specified on all boundaries. The system which has found most extensive use in generating boundary fitted grids is made up of Poisson equations, of which the functions P and Q provide a means for controlling the spacing and density of grid lines in the coordinate system. While questions remain concerning the existence and uniqueness of orthogonal systems, the generating method presented adds to the available, useful techniques for constructing these systems.

Coleman, R. M.

Vortex methods for two- and three-dimensional flow simulations

The point vortex and vortex blob methods for two dimensional flows are presented. Several results are discussed concerning the numerical analysis of the latter scheme, e.g., the preservation of globally conserved quantities and the analysis of the spatial discretization error resulting from the convection of fixed blobs of vorticity. An application to the two dimensional mixing layer is briefly described. The contour dynamics method is also discussed. The simulation of three dimensional flows with vortex methods is discussed. A natural way to represent the vorticity is in the form of closed tubes of filaments of vorticity, although other schemes are examined. Applications to aircraft trailing vortices and to a turbulent spot in a laminar boundary layer are presented. Hybrid schemes that use an Eulerian mesh to solve the Poisson equation for the velocity field are discussed. The goal of these schemes is to avoid the high cost of the Biot-Savart integration if many vortex elements are used while enjoying most of the advantages of pure Lagrangian schemes.

Leonard, A.

Trends and Techniques for Space Base Electronics

Simulations of various phosphorus and boron diffusions in SOS were completed and a sputtering system, furnaces, and photolithography related equipment were set up. Double layer metal experiments initially utilized wet chemistry techniques. By incorporating ultrasonic etching of the vias, premetal cleaning a modified buffered HF, phosphorus doped vapox, and extended sintering, yields of 98% were obtained using the standard test pattern. A two dimensional modeling program was written for simulating short channel MOSFETs with nonuniform substrate doping. A key simplifying assumption used is that the majority carriers can be represented by a sheet charge at the silicon dioxide silicon interface. Although the program is incomplete, the two dimensional Poisson equation for the potential distribution was achieved. The status of other Z-D MOSFET simulation programs is summarized.

Trotter, J. D.

Some features of inverted-V events as seen from simulated double layers

Results from a numerical simulation of a double layer show some features similar to those of inverted-V events. The strong heating of thermal and precipitating electrons is observed along with extremely low frequency fluctuations found during inverted-V events. It is suggested that after the acceleration of auroral electrons by the double layer, the precipitating free electrons are heated by the nonlinear effects of the electron beam plasma instability. Fluctuations and pulsations of auroral electron fluxes during auroral events are caused by a relaxation type of oscillation. The finite extent of one dimensional plasma is simulated by solving the Vlasov and Poisson equations as an initial and boundary value problem.

Singh, N.

A new global operator for two-particle delta functions

A new type of global operator to be used in evaluating matrix elements of two-particle delta functions is introduced. It is based, like the Trivedi one-particle operator, on the Poisson equation and is easier to apply than the method of Hiller, Sucher and Feinberg. After a test in the helium isoelectronic sequence, the new method is applied successfully to the interesting problem of hyperfine structure in muonic helium.

Drachman, R. J.

Double layers on auroral field lines

Time-stationary solutions to the Vlasov-Poisson equation for ion holes and double layers were examined along with particle simulations which pertain to recent observations of small amplitude (e phi)/t sub e approx. 1 electric field structures on auroral field lines. Both the time-stationary analysis and the simulations suggest that double layers evolve from holes in ion phase space when their amplitude reaches (e phi)/t sub e approx. 1. Multiple small amplitude double layers which are seen in long simulation systems and are seen to propagate past spacecraft may account for the acceleration of plasma sheet electrons to produce the discrete aurora.

Hudson, M. K.

Effects of auroral-particle anisotropies and mirror forces on high-latitude electric fields

It is noted that, for most of the mechanisms for the strong electric fields that characterize the narrow regions in which there is acceleration and precipitation of ring current and/or plasma-sheet plasma, certain effects must be taken into account in simulations of auroral electric fields. The effects are those of auroral particle anisotropy, of mirror forces due to the inhomogeneous geomagnetic field, of auroral electron backscatter by the atmosphere, and of electron trapping by the combination of magnetic mirroring and electrostatic forces. What is more, the effects of the very strong perpendicular electric field must also be taken into account in a kinetic description of the Poisson equation in order to achieve a unified theory of the auroral electrostatic structure. Progress in these areas during the past few years is reviewed. It is shown that particle anisotropies and mirror forces can account for some basic electrostatic features of the quiet arc, while additional effects may be occurring in strong events in which the parallel potential drop is more than about 10 kV.

Chiu, Y. T.

Dynamical features of moving double layers

Numerical simulations of the dynamics of double layers that form in response to an applied potential drop across a plasma are carried out for plasmas of lengths much greater than previously considered. The evolution of finite-size plasmas of lengths 100, 200 and 400 Debye lengths at the low-potential end of the simulation region was followed by the solution of the Vlasov and Poisson equations. Results show that the double layer moves towards the high-potential side of the layer, at a speed directly dependent on the length of the plasma. The moving double layer is accompanied by a moving density front which is similar to an ion acoustic shock created by the expansion of a high-density plasma into a low-density plasma, as well as plasma heating and evacuation, electron and ion current interruptions and recovery. Double layer formation, motion, and accompanying phenomena are found to repeat periodically. Results are in good agreement with laboratory experiments, and may be used to explain certain phenomena observed in auroral plasmas.

Singh, N.

Tenth NASTRAN User's Colloquium

The development of the NASTRAN computer program, a general purpose finite element computer code for structural analysis, was discussed. The application and development of NASTRAN is presented in the following topics: improvements and enhancements; developments of pre and postprocessors; interactive review system; the use of harmonic expansions in magnetic field problems; improving a dynamic model with test data using Linwood; solution of axisymmetric fluid structure interaction problems; large displacements and stability analysis of nonlinear propeller structures; prediction of bead area contact load at the tire wheel interface; elastic plastic analysis of an overloaded breech ring; finite element solution of torsion and other 2-D Poisson equations; new capability for elastic aircraft airloads; usage of substructuring analysis in the get away special program; solving symmetric structures with nonsymmetric loads; evaluation and reduction of errors induced by Guyan transformation.

Source record

Test results of modified electrical charged particle generator for application to fog dispersal

Modifications to a charged particle generator for use in fog dispersal applications were made and additional testing carried out. The modified nozzle, however, did not work as planned, and reported results are the unmodified nozzle. The addition of a positive displacement pump to supply the liquid water was highly successful. Measurements of the generator output current were made with a cylindrical collector system as well as with the needle probe used in previous studies. Measurements with the cylindrical collector and the needle probe showed identical agreement within the variability of the experiment. A high-voltage prove was purchased, and measurements of the corona voltage as well as the voltage variation in the charged particle jet were made. Electric fields in the vertical direction on the order of 1,000,000 v/m were measured. The voltage distribution along the centerline of the jet was compared with the numerical solutions of the Poisson equation and showed very good agreement. Velocity measurements using a pitot tube were made. The resulting measurements were compared with theoretical and other reported experimental results. The measured data showed the appropriate trends and agreed well with reported results. Based on the measured current-to-mass ratio from the charged particle generator, a calculation of the average droplet size was made. Droplet sizes were estimated to range between 0.8 and 0.4 microns. Using measured data, an analysis of the height to which the droplet can be dispersed by the charged particle generator was made. Although the mathematical model is highly simplified, the results indicated that particles would achieve heights on the order of 80 m.

Frost, W.

Solitary waves and double layers on auroral field lines

Time stationary solutions to the Vlasov-Poisson equations for ion holes and double layers are examined along with particle simulations that pertain to recent observations of small amplitude electric field structures on auroral field lines. Both the time stationary analysis and the simulations suggest that the observed double layers evolve from holes in ion phase space. Multiple small amplitude double layers, as seen in long simulation systems, are observed to propagate past the spacecraft and may account for the acceleration of plasma sheet electrons to produce inverted-V precipitation.

Hudson, M. K.

Vectorized multigrid Poisson solver for the CDC CYBER 205

The full multigrid (FMG) method is applied to the two dimensional Poisson equation with Dirichlet boundary conditions. This has been chosen as a relatively simple test case for examining the efficiency of fully vectorizing of the multigrid method. Data structure and programming considerations and techniques are discussed, accompanied by performance details.

Barkai, D.

An alternative to reduction of surface pressure to sea level

The pitfalls of the present method of reducing surface pressure to sea level are reviewed, and an alternative, adjusted pressure, P, is proposed. P is obtained from solution of a Poisson equation over a continental region, using the simplest boundary condition along the perimeter or coastline where P equals the sea level pressure. The use of P would avoid the empiricisms and disadvantages of pressure reduction to sea level, and would produce surface pressure charts which depict the true geostrophic wind at the surface.

Deardorff, J. W.

Spacecraft-environment interaction: The environmental plasma aspect

The effects involved in the interaction between an obstacle and a space plasma can be divided into: (1) effects on the obstacle itselt (i.e., its charging); and (2) effects on the environmental plasma due to the motion of the obstacle (i.e., the creation of shocks ahead of the obstacle and complicated wakes behind the obstacle). In the wake (or antisolar direction), plasma oscillations are excited and instabilities, wave-particle interactions, turbulence, etc., are believed to take place. The effects on the obstacle and on the environmental space plasma are coupled. Hence, simultaneous solutions to the Vlasov (or Boltzmann) and Poisson equations are sought. To obtain realistic solutions of practical use, three-dimensional and time-dependent models of the interaction are needed. Achieving the latter is indeed not simple.

Samir, U.

Numerical simulation of spacecraft charging by impact-induced plasmas during a cometary flyby

A numerical model is developed for the interaction of a cometary probe, such as Giotto, with its environment, i.e., dust particles and gas. The spacecraft was set on a course to pass the comet at a velocity of 69 km/sec, so a chance existed that a potential field would form around the spacecraft and block lower energy particles from reaching the spacecraft instruments. The motion of electrons and ions is traced as a function of time to examine the evolution of the electric field, electric potential and the total space charge distributions on the surface of the spacecraft and its environment. Account is taken of the density of the particles and gas molecules at various distances from the comet, the collision energies involved, and the Giotto geometry. A solution is defined for the Poisson equation to describe the evolution of the plasma around Giotto, including the effects of ion collisions with the Al bumper protecting the spacecraft. The simulation predicts formation of an ion wake behind Giotto and the evolution of a positive potential on the order of 10 V around the spacecraft, i.e., sufficient for a positive potential barrier near the surface of the spacecraft.

Thiemann, H.