Search NASA⌕ Search

SEARCH · Search NASA

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Research on Orbital Plasma-Electrodynamics (ROPE)

Since the development of probe theory by Langmuir and Blodgett, the problem of current collection by a charged spherically or cylindrically symmetric body has been investigated by a number of authors. This paper overviews the development of a fully three-dimensional particle simulation code which can be used to understand the physics of current collection in three dimensions and can be used to analyze data resulting from the future tethered satellite system (TSS). According to the TSS configurations, two types of particle simulation models were constructed: a simple particle simulation (SIPS) and a super particle simulation (SUPS). The models study the electron transient response and its asymptotic behavior around a three dimensional, highly biased satellite. The potential distribution surrounding the satellite is determined by solving Laplace's equation in the SIPS model and by solving Poisson's equation in the SUPS model. Thus, the potential distribution in space is independent of the density distribution of the particles in the SUPS model but it does depend on the density distribution of the particles in the SUPS model. The evolution of the potential distribution in the SUPS model is described. When the spherical satellite is charged to a highly positive potential and immersed in a plasma with a uniform magnetic field, the formation of an electron torus in the equatorial plane (the plane in perpendicular to the magnetic field) and elongation of the torus along the magnetic field are found in both the SIPS and the SUPS models but the shape of the torus is different. The areas of high potential that exist in the polar regions in the SUPS model exaggerate the elongation of the electron torus along the magnetic field. The current collected by the satellite for different magentic field strengths is investigated in both models. Due to the nonlinear effects present in SUPS, the oscillating phenomenon of the current collection curve during the first 10 plasma periods can be seen (this does not appear in SIPS). From the parametric studies, it appears that the oscillating phenomenon of the current collection curve occurs only when the magnetic field strength is less than 0.2 gauss for the present model.

Wu, S. T.↗
NASA NTRS1975

Advances in Application of Fast Semidirect Computational Methods in Transonic Flow

This paper is intended as a review and summary of the advances made in a recently developed approach for rapid numerical solution of the equations of inviscid transonic aerodynamics. The investigation has been limited to two-dimensional, steady, inviscid flow over airfoils in a subsonic free stream, with emphasis on development of a rapid computational technique, rather than on generality of application. The approach uses finite-difference algorithms called "fast direct elliptic solvers" within an iteration scheme. "Direct" means that the entire computation field is solved at once, rather than in successive traverses over the field as in a point- or line-relaxation method. Such an iterative method is referred to as "semidirect." The iterative convergence can be faster than in other relaxation methods because changes are felt simultaneously at all points in each succeeding iteration. Direct elliptic solvers and semidirect methods have restrictions, but these are gradually being removed. Direct solvers were first developed for solving Poisson's equation on a rectangle without interior boundaries. A method to treat first-order systems, a direct Cauchy-Riemann solver has also been developed. Numerical treatment of part of a system of nonlinear equations by a Poisson solver has been reported. Also Poisson solvers in semidirect methods were used for nonseparable elliptic equations. The semidirect method was extended to the solution of a problem of mixed type, where the improved Murman-Cole transonic small-disturbance difference equations were solved. A slightly supercritical flow over a biconvex airfoil was treated successfully, but the iterations did not converge for more strongly supercritical conditions In another work the addition of terms ot both sides of the difference equations stabilized the iteration for supercritical conditions with large supersonic zones. For this, the Cauchy-Riemann solver was revised to incl,ude the needed terms. Most recently, the evaluation of parameters for rapid convergence and comparisons, with Murman's line-relaxation method was described. The method was extended to full second order accuracy in a fully conservative formulation in another work.

Martin, E. Dale↗
NASA NTRS1979

Analytical study of the time dependent spacecraft-plasma interaction

A study of the time dependent interaction of an initially uniform equilibrium plasma with a plane conducting surface was made in order to achieve a more complete understanding of the dynamics of the charging process and of the approach to the floating potential on the surface. Numerical solutions of the cold ion equations of motion in conjunction with equilibrium electrons and Poisson's equation show the formation of an ion-rich sheath near the surface and the coupling of the non-neutral region to the undisturbed plasma through a quasineutral rarefaction. Analytical treatment of the quasineutral region showed excellent agreement (within 1 percent) with the numerical results.

Cipolla, J. W., Jr.↗
NASA NTRS1977

An efficient numerical technique for calculating thermal spreading resistance

An efficient numerical technique for solving the equations resulting from finite difference analyses of fields governed by Poisson's equation is presented. The method is direct (noniterative)and the computer work required varies with the square of the order of the coefficient matrix. The computational work required varies with the cube of this order for standard inversion techniques, e.g., Gaussian elimination, Jordan, Doolittle, etc.

Gale, E. H., Jr.↗
NASA NTRS1996

Multigrid, Fractional-Step Computation Of Flow

Speed of computer code solving three-dimensional Navier-Stokes equations of flow of incompressible fluid by fractional-step method increased significantly by use of multigrid procedures. In method, equations solved on general nonorthogonal curvilinear coordinate grid, using volume fluxes. At each time step, computations performed in two fractional steps. In first step, equations of conservation of momentum solved by use of gradient of pressure from previous time step via explicit approximate-factorization method, yielding approximate flow field that does not satisfy equation of conservation of mass. In second step, discrete Poisson-like equation with Neumann-type boundary conditions, formed by combining equations of conservation of momentum and mass, solved iteratively.

Kwak, Dochan↗
NASA NTRS1988

A solution method for the unsteady incompressible Navier-Stokes equations in generalized coordinate systems

A solution method based on a fractional step approach is developed for obtaining time-dependent solutions of the three-dimensional, incompressible Navier-Stokes equations in generalized coordinate systems. The governing equations are discretized conservatively by finite volumes using a staggered mesh system. The primitive variable formulation uses the volume fluxes across the faces of each computational cell as dependent variables. This procedure, combined with accurate and consistent approximations of geometric parameters, is done to satisfy the discretized mass conservation equation to machine accuracy as well as to gain favorable convergence properties of the Poisson solver. The discretized equations are second-order-accurate in time and space and no smoothing terms are added. An approximate-factorization scheme is implemented in solving the momentum equations. A novel ZEBRA scheme with four-color ordering is devised for the efficient solution of the Poisson equation. Several two and three-dimensional solutions are compared with other numerical and experimental results to validate the present method.

Rosenfeld, Moshe↗
NASA NTRS1970

High frequency Hall current instability

The Hall current electrojet instability was studied to determine if it generates VLF hiss in the auroral zone of magnetic substorms. Using the kinetic equations with particle conserving Krook collision operator for electrons and ions, and Poisson's equation, it is shown that in the low frequency region, the range of frequency of unstable modes increases as the electron density increases, and that the high frequency component must be taken into account in calculations of turbulent saturation of the Hall current instability driven beyond threshold. It is concluded that Hall currents can generate waves in the VLF hiss frequency band which are related to local auroral electrojet activity.

Lee, K.↗
NASA NTRS1995

Two-dimensional mesh embedding for Galerkin B-spline methods

A number of advantages result from using B-splines as basis functions in a Galerkin method for solving partial differential equations. Among them are arbitrary order of accuracy and high resolution similar to that of compact schemes but without the aliasing error. This work develops another property, namely, the ability to treat semi-structured embedded or zonal meshes for two-dimensional geometries. This can drastically reduce the number of grid points in many applications. Both integer and non-integer refinement ratios are allowed. The report begins by developing an algorithm for choosing basis functions that yield the desired mesh resolution. These functions are suitable products of one-dimensional B-splines. Finally, test cases for linear scalar equations such as the Poisson and advection equation are presented. The scheme is conservative and has uniformly high order of accuracy throughout the domain.

Shariff, Karim↗
NASA NTRS1988

Dynamics of localized ion-acoustic waves in a magnetized plasma

The evolution of negative potential pulses in a magnetized plasma is studied. A three-dimensional nonlinear ion-acoustic wave equation, including nonstationary effects of reflected electrons, has been derived from the Poisson-Vlasov equations with uniform magnetic field. The studies show that a negative potential pulse can be enhanced by drifting electrons. The growing pulse develops asymmetrically with an oscillatory precursor and a local potential jump resembling the early phase of weak double layer formation.

Qian, S.↗
NASA NTRS1975

A rigorous solution of the Navier-Stokes equations for unsteady viscous flow at high Reynolds numbers around oscillating airfoils

A method based on the Navier-Stokes equations was developed for analyzing the unsteady incompressible viscous flow around oscillating airfoils at high Reynolds numbers. The Navier-Stokes equations have been integrated in their classical Helmholtz vorticity transport equation form, and the instantaneous velocity field at each time step was determined by the solution of Poisson's equation. A refined finite element was utilized to allow for a conformable solution of the stream function and its first space derivatives at the element interfaces. A corresponding set of accurate boundary conditions was applied; thus obtaining a rigorous solution for the velocity field. The details of the computational procedure and examples of computed results describing the unsteady flow characteristics around the airfoil are presented.

Bratanow, T.↗
NASA NTRS1986

Fourier analysis of the SOR iteration

The SOR iteration for solving linear systems of equations depends upon an overrelaxation factor omega. It is shown that for the standard model problem of Poisson's equation on a rectangle, the optimal omega and corresponding convergence rate can be rigorously obtained by Fourier analysis. The trick is to tilt the space-time grid so that the SOR stencil becomes symmetrical. The tilted grid also gives insight into the relation between convergence rates of several variants.

Leveque, R. J.↗
NASA NTRS1991

Particle flows to shape and voltage surface discontinuities in the electron sheath surrounding a high voltage solar array in LEO

This paper discusses the numerical modeling of electron flows from the sheath surrounding high positively biased objects in LEO (Low Earth Orbit) to regions of voltage or shape discontinuity on the biased surfaces. The sheath equations are derived from the Two-fluid, Warm Plasma Model. An equipotential corner and a plane containing strips of alternating voltage bias are treated in two dimensions. A self-consistent field solution of the sheath equations is outlined and is pursued through one cycle. The electron density field is determined by numerical solution of Poisson's equation for the electrostatic potential in the sheath using the NASCAP-LEO relation between electrostatic potential and charge density. Electron flows are calculated numerically from the electron continuity equation. Magnetic field effects are not treated.

Metz, Roger N.↗
NASA NTRS1980

Electrostatic model of a quiet auroral arc

The physics and equations used in the proposed electrostatic model of a quiet auroral arc are discussed. The procedure for obtaining approximate solutions of Poisson's equation coupled to ionosphere current conservation is given. A comparison with a specific discrete auroral arc as seen from the auroral satellite S3-3 is presented.

Chiu, Y. T.↗
NASA NTRS1988

I-BIEM, an iterative boundary integral equation method for computer solutions of current distribution problems with complex boundaries: A new algorithm. I - Theoretical

A new algorithm for an iterative computation of solutions of Laplace's or Poisson's equations in two dimensions, using Green's second identity, is presented. This algorithm converges strongly and geometrically and can be applied to curved, irregular, or moving boundaries with nonlinear and/or discontinuous boundary conditions. It has been implemented in Pascal on a number of micro- and minicomputers and applied to several geometries. Cases with known analytic solutions have been tested. Convergence to within 0.1 percent to 0.01 percent of the theoretical values are obtained in a few minutes on a microcomputer.

Cahan, B. D.↗
NASA NTRS1973

Finite element analysis of unsteady incompressible flow around an oscillating obstacle of arbitrary shape.

An analytical procedure based on Navier-Stokes equations was developed for representing unsteady flow patterns around oscillating obstacles. A variational formulation of the Helmholtz vorticity equation was discretized in finite element form and integrated numerically. At each step of the numerical integration the velocity field around the obstacle was determined from the finite element solution of Poisson's equation. The time-dependent boundary conditions around the oscillating obstacle were introduced as external constraints at each time step of the numerical integration. The obtained results for a cylinder and an airfoil were illustrated in the form of streamlines and vorticity and pressure distributions.

Bratanow, T.↗
NASA NTRS1981

Three-dimensional structures and turbulence closure of the wake developing in a wall shear layer

The turbulent wake interacting with the rotating wall shear layer is investigated analytically and numerically. The turbulent wakes of the rotating blades in a compressor which are interacting with the rotating hub-wall boundary layer are analyzed. A modified version of the closure model of the pressure-strain correlation term in the Reynolds stress transport equation is developed to predict the effect of rotation, which is appreciable for the present flow because the thick hub-wall boundary layer is interacting with the rotor wake. It is noted that the Poisson type equation for the pressure-strain correlation has an extra rotation term when the entire flow field is rotating. This extra rotation term is modeled to accommodate the effect of rotation. In addition, the standard correction for the wall effect is incorporated for the utilized Reynolds stress closure model. The rotation-modified Reynolds stress closure model is used to predict the present flow, and the predictions are compared with the experimental data. The experimental data reveal that the characteristics of the three-dimensional turbulent wake interacting with the wall shear layer are considerably altered by the effects of the wall and the rotation. These features are predicted with good accuracy by the turbulence closure model developed.

Hah, C.↗
NASA NTRS1974

Numerical calculations of velocity and pressure distribution around oscillating airfoils

An analytical procedure based on the Navier-Stokes equations was developed for analyzing and representing properties of unsteady viscous flow around oscillating obstacles. A variational formulation of the vorticity transport equation was discretized in finite element form and integrated numerically. At each time step of the numerical integration, the velocity field around the obstacle was determined for the instantaneous vorticity distribution from the finite element solution of Poisson's equation. The time-dependent boundary conditions around the oscillating obstacle were introduced as external constraints, using the Lagrangian Multiplier Technique, at each time step of the numerical integration. The procedure was then applied for determining pressures around obstacles oscillating in unsteady flow. The obtained results for a cylinder and an airfoil were illustrated in the form of streamlines and vorticity and pressure distributions.

Bratanow, T.↗
NASA NTRS1986

Numerical simulation of transonic separated flows over low-aspect ratio wings

Transonic flow fields about a low-aspect-ratio advanced technology wing have been computed using a viscous/inviscid zonal approach. The flow field near the wing where viscous effects are important was solved using the 'Reynolds-Averaged Navier-Stokes Equations' in 'thin-layer' form. The Euler equations were used to determine the flow field in regions away from the wing where viscous effects are insignificant. A zonal grid using an H-H topology was generated around the wing by first solving a set of Poisson's equations for the global grid. This grid was then subdivided into separate zones of viscous or inviscid flow as suggested by the flow physics. A series of flow cases were computed and compared with corresponding sets of experimental data. All cases showed good agreement with experiment in terms of the pressure field. Also, a good correlation between computed separated surface flow and experimental oil flow was obtained.

Kaynak, U.↗