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 289 records · Page 16

Equilibrium stellar systems with spindle singularities

Equilibrium sequences of axisymmetric Newtonian clusters that tend toward singular states are constructed. The distribution functions are chosen to be of the form f = f(E, Jz). The numerical method then determines the density and gravitational potential self-consistently to satisfy Poisson's equation. For the prolate models, spindle singularities arise from the depletion of angular momentum near the symmetry axis. While the resulting density enhancement is confined to the region near the axis, the influence of the spindle extends much further out through its tidal gravitational field. Centrally condensed prolate clusters may contain strong-field regions even though the spindle mass is small and the mean cluster eccentricity is not extreme. While the calculations performed here are entirely Newtonian, the issue of singularities is an important topic in general relativity. Equilibrium solutions for relativistic star clusters can provide a testing ground for exploring this issue. The methods used in this paper for building nonspherical clusters can be extended to relativistic systems.

Shapiro, Stuart L.↗

Drag Minimization for Wings and Bodies in Supersonic Flow

The minimization of inviscid fluid drag is studied for aerodynamic shapes satisfying the conditions of linearized theory, and subject to imposed constraints on lift, pitching moment, base area, or volume. The problem is transformed to one of determining two-dimensional potential flows satisfying either Laplace's or Poisson's equations with boundary values fixed by the imposed conditions. A general method for determining integral relations between perturbation velocity components is developed. This analysis is not restricted in application to optimum cases; it may be used for any supersonic wing problem.

Heaslet, Max A↗

Algorithms for parallel and vector computations

This is a final report on work performed under NASA grant NAG-1-1112-FOP during the period March, 1990 through February 1995. Four major topics are covered: (1) solution of nonlinear poisson-type equations; (2) parallel reduced system conjugate gradient method; (3) orderings for conjugate gradient preconditioners, and (4) SOR as a preconditioner.

Ortega, James M.↗

A comparison of two algorithms for simulating collisionless systems

Two completely different simulation algorithms are compared by applying them to the same stellar dynamical problems: one is a self-consistent field (SCF) method for solving Poisson's equation and the other is a phase-space method for integrating the collisionless Boltzmann equation. We consider simulations of spherical stellar systems which are initially far from equilibrium and relax to their final states by gravitational collapse. The initial conditions consist of either uniform-density spheres or nonequilibrium models having Plummer density profiles, in which velocity dispersions are assigned according to given virial ratios. If a few tens of radial expansion terms with hundreds of thousands of particles are used in the SCF code, excellent agreement is found between the results it generates and those obtained with the phase-space solver, provided that a sufficiently large number of grid cells are employed with the latter. These findings imply that for simulating collisionless systmes over many dynamical times, the SCF approach based on sampling phase space is competitive with the approach treating phase space as a continuous fluid. The results of our tests make it possible to estimate the number of particles and basis functions required in situations like those modeled. Limitations of the SCF method and the choice of an optimal set of basis functions are also discussed.

Hozumi, Shunsuke↗

Multiple Coarse Grid Multigrid Methods for Solving Elliptic Problems

In this paper we describe some classes of multigrid methods for solving large linear systems arising in the solution by finite difference methods of certain boundary value problems involving Poisson's equation on rectangular regions. If parallel computing systems are used, then with standard multigrid methods many of the processors will be idle when one is working at the coarsest grid levels. We describe the use of Multiple Coarse Grid MultiGrid (MCGMG) methods. Here one first constructs a periodic set of equations corresponding to the given system. One then constructs a set of coarse grids such that for each grid corresponding to the grid size h there are four grids corresponding to the grid size 2*h. Multigrid operations such as restriction of residuals and interpolation of corrections are done in parallel at each grid level. For suitable choices of the multigrid operators the MCGMG method is equivalent to the Parallel Superconvergent MultiGrid (PSMG) method of Frederickson and McBryan. The convergence properties of MCGMG methods can be accurately analyzed using spectral methods.

Xiao, Shengyou↗

Locations of Halide Ions in Tetragonal Lysozyme Crystals

Anions play an important role in the crystallization of lysozyme, and are known to bind to the crystalline protein. Previous studies employing X-ray crystallography had found one chloride ion binding site in the tetragonal crystal form of the protein and four nitrate ion binding sites in the monoclinic form. Studies using other approaches have reported more chloride ion binding sites, but their locations were not known. Knowing the precise location of these anions is also useful in determining the correct electrostatic fields surrounding the protein. In the first part of this study the anion positions in the tetragonal form were determined from the difference Fourier map obtained from the lysozyme crystals grown in bromide and chloride solutions under identical conditions. The anion locations were then obtained from standard crystallographic methods and five possible anion binding sites were found in this manner. The sole chloride ion binding site found in previous studies was confirmed. The remaining four sites were new ones for tetragonal lysozyme crystals. However, three of these new sites and the previously found one corresponded to the four unique binding sites found for nitrate ions in monoclinic crystals. This suggests that most of the anion binding sites in lysozyme remain unchanged, even when different anions and different crystal forms of lysozyme are employed. It is unlikely that there are many more anions in the tetragonal lysozyme crystal structure. Assuming osmotic equilibrium it can be shown that there are at most three more anions in the crystal channels. Some of the new anion binding sites found in this study were, as expected, in pockets containing basic residues. However, some of them were near neutral, but polar, residues. Thus, the study also showed the importance of uncharged, but polar groups, on the protein surface in determining its electrostatic field. This was important for the second part of this study where the electrostatic field surrounding the protein was accurately determined. This was achieved by solving the linearized version of the Poisson-Boltzmann equation for the protein in solution. The solution was computed employing the commercial code Delphi which uses a finite difference technique. This has recently become available as a module in the general protein visualization code Insight II. Partial charges were assigned to the polar groups of lysozyme for the calculations done here. The calculations showed the complexity of the electrostatic field surrounding the protein. Although most of the region near the protein surface had a positive field strength, the active site cleft was negatively charged and this was projected a considerable distance. This might explain the occurrence of "head-to-side" interactions in the formation of lysozyme aggregates in solution. Pockets of high positive field strength were also found in the vicinity of the anion locations obtained from the crystallographic part of this study, confirming the validity of these calculations. This study clearly shows not only the importance of determining the counterion locations in protein crystals and the electrostatic fields surrounding the protein, but also the advantage of performing them together.

Lim, Kap↗

Statistically Reliable 'Atomistic' Simulation of Sub 100 nm MOSFETs

A 3D 'atomistic' simulation technique to study random impurity induced threshold voltage lowering and fluctuations in sub 0. 1 micron MOSFETs is presented. It allows statistical analysis of random impurity effects down to the individual impurity level-Efficient algorithms based on a single solution of Poisson's equation, followed by the solution of a simplified current continuity equation are used in the simulations.

Asenov, Asen↗

Efficient 3D 'Atomistic' Simulation Technique for Studying of Random Dopant Induced Threshold Voltage Lowering and Fluctuations in Decanano MOSFETs

A 3D 'atomistic' simulation technique to study random dopant induced threshold voltage lowering and fluctuations in sub 0.1 micron MOSFETs is presented. It allows statistical analysis of random impurity effects down to the individual impurity level. Efficient algorithms based on a single solution of Poisson's equation, followed by the solution of a simplified current continuity equation are used in the simulations.

Asenov, Asen↗

Plasma Reactor Modeling and Validation Experiments

Plasma processing is a key processing stop in integrated circuit manufacturing. Low pressure, high density plum reactors are widely used for etching and deposition. Inductively coupled plasma (ICP) source has become popular recently in many processing applications. In order to accelerate equipment and process design, an understanding of the physics and chemistry, particularly, plasma power coupling, plasma and processing uniformity and mechanism is important. This understanding is facilitated by comprehensive modeling and simulation as well as plasma diagnostics to provide the necessary data for model validation which are addressed in this presentation. We have developed a complete code for simulating an ICP reactor and the model consists of transport of electrons, ions, and neutrals, Poisson's equation, and Maxwell's equation along with gas flow and energy equations. Results will be presented for chlorine and fluorocarbon plasmas and compared with data from Langmuir probe, mass spectrometry and FTIR.

Meyyappan, M.↗

Electrical double layers at the oil/water interface

This review presents the historical development and current status of the theory of the electrical double layer at a liquid/liquid interface. It gives rigorous thermodynamic definitions of all basic concepts related to liquid interfaces and to the electrical double layer. The difference between the surface of a solid electrode and the interface of two immiscible electrolyte solutions (ITIES) is analyzed in connection to their electrical properties. The most important classical relationships for the electrical double layer are presented and critically discussed. The generalized adsorption isotherm is derived. After a short review of the classical Gouy-Chapman and Verwey-Niessen models, more recent developments of the double layer theory are presented. These include effects of variable dielectric permittivity, nonlocal electrostatics, hydration forces, the modified Poisson-Boltzmann equation and the ion-dipole plasma. The relative merits of different theories are estimated by comparing them with computer simulation of the ITIES and electrical double layer. Special attention is given to the structure of ITIES and its variation due to adsorption of ions and amphiphilic molecules.

Non-NASA Center↗

A fast Poisson solver for unsteady incompressible Navier-Stokes equations on the half-staggered grid

In this paper, a fast Poisson solver for unsteady, incompressible Navier-Stokes equations with finite difference methods on the non-uniform, half-staggered grid is presented. To achieve this, new algorithms for diagonalizing a semi-definite pair are developed. Our fast solver can also be extended to the three dimensional case. The motivation and related issues in using this second kind of staggered grid are also discussed. Numerical testing has indicated the effectiveness of this algorithm.

Golub, G. H.↗

Characteristics of 0.8- and 0.2-microns gate length In(x)Ga(1-x) As/In(0.52)Al(0.48)As/InP (0.53 less than or equal to x less than or equal to 0.70) modulation-doped field-effect transistors at cryogenic temperatures

The performance characteristics of InP-based In(x)Ga(1-x)As/In(0.52)Al(0.48)As (0.53 is less than or equal to x is less than or equal to 0.70) pseudomorphic modulation-doped field-effect transistors (MODFET's) as a function of strain in the channel, gate, length, and temperature were investigated analytically and experimentally. The strain in the channel was varied by varying the In composition x. The temperature was varied in the range of 40-300 K and the devices have gate lengths L(sub g) of 0.8 and 0.2 microns. Analysis of the device was done using a one-dimensional self consistent solution of the Poisson and Schroedinger equations in the channel, a two-dimensional Poisson solver to obtain the channel electric field, and a Monte Carlo simulation to estimate the carrier transit times in the channel. An increase in the value of the cutoff frequency is predicted for an increase in In composition, a decrease in temperature, and a decrease in gate length. The improvements seen with decreasing temperature, decreasing gate length, and increased In composition were smaller than those predicted by analysis. The experimental results on pseudomorphic InGaAs/InAlAs MODFET's showed that there is a 15-30 percent improvement in cutoff frequency in both the 0.8- and 0.2-micron gate length devices when the temperature is lowered from 300 to 40 K.

Lai, Richard↗

A numerical method for solving the Vlasov equation

A numerical procedure is derived for the solution of the Vlasov-Poisson system of equations in two phase-space variables. Derivatives with respect to the phase-space variables are approximated by a weighted sum of the values of the distribution function at property chosen neighboring points. The resulting set of ordinary differential equations is then solved by using an appropriate time intergration scheme. The accuracy of the proposed method is tested with some simple model problems. The results for the free streaming case, linear Landau damping, and nonlinear Landau damping are investigated and compared with those of the splitting scheme. The proposed method is found to be very accurate and efficient.

Satofuka, N.↗

Consistent boundary conditions for Reduced Navier-Stokes (RNS) scheme applied to 3-dimensional internal viscous flows

A consistent and efficient set of boundary conditions is developed for the multi-sweep space-marching pressure-elliptic Reduced Navier-Stokes (RNS) scheme as applied for 3-D internal viscous flow problems. No-slip boundary conditions are directly imposed on the solid walls. There is no iteration procedure required in the cross plane to ensure mass conservation across each marching plane. The finite difference equations forming the coefficient matrix are ordered such that the surface normal velocity is specified on all the solid walls; unlike external flows, a pressure boundary condition in the cross plane is not required. Since continuity is directly satisfied at all points in the flow domain, the first order momentum equations can be solved directly for the pressure without the need for a Poisson pressure correction equation. The procedure developed herein can also be applied with periodic boundary conditions. The analysis is given for general compressible flows. Incompressible flow solutions are obtained, for straight and curved ducts of square cross section, to validate the procedure. These solutions are used to demonstrate the applicability of the RNS scheme, with the improved boundary conditions for internal flows with strong interaction, as would be encountered in ducts and turbomachinery geometries.

Reddy, D. R.↗

Train small, model big: Scalable physics simulators via reduced order modeling and domain decomposition

Numerous cutting-edge scientific technologies originate at the laboratory scale, but transitioning them to practical industry applications is a formidable challenge. Traditional pilot projects at intermediate scales are costly and time-consuming. An alternative, the pilot-scale model, relies on high-fidelity numerical simulations, but even these simulations can be computationally prohibitive at larger scales. To overcome these limitations, we propose a scalable, physics-constrained reduced order model (ROM) method. The ROM identifies critical physics modes from small-scale unit components, projecting governing equations onto these modes to create a reduced model that retains essential physics details. We also employ Discontinuous Galerkin Domain Decomposition (DG-DD) to apply ROM to unit components and interfaces, enabling the construction of large-scale global systems without data at such large scales. Here this method is demonstrated on the Poisson and Stokes flow equations, showing that it can solve equations about 15–40 times faster with only ~1% relative error. Furthermore, ROM takes one order of magnitude less memory than the full order model, enabling larger scale predictions at a given memory limitation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

SINDA '85/FLUINT, Version 2.2

SINDA, Systems Improved Numerical Differencing Analyzer, software system developed for solving physical problems governed by diffusion-type equations modeled by lumped-parameter representations. Used as general thermal analyzer with resistor-and-capacitor-network representations, adapted to wide range of problems represented by such differential equations as Fourier, Poisson, or LaPlace. Software system updated. FLUINT, FLUid INTegrator, is advanced, one-dimensional fluid-analysis program solving equations of arbitrary fluid-flow networks. Written in FORTRAN and ASSEMBLER.

Cullimore, Brent↗

Incompressible spectral-element method: Derivation of equations

A fractional-step splitting scheme breaks the full Navier-Stokes equations into explicit and implicit portions amenable to the calculus of variations. Beginning with the functional forms of the Poisson and Helmholtz equations, we substitute finite expansion series for the dependent variables and derive the matrix equations for the unknown expansion coefficients. This method employs a new splitting scheme which differs from conventional three-step (nonlinear, pressure, viscous) schemes. The nonlinear step appears in the conventional, explicit manner, the difference occurs in the pressure step. Instead of solving for the pressure gradient using the nonlinear velocity, we add the viscous portion of the Navier-Stokes equation from the previous time step to the velocity before solving for the pressure gradient. By combining this 'predicted' pressure gradient with the nonlinear velocity in an explicit term, and the Crank-Nicholson method for the viscous terms, we develop a Helmholtz equation for the final velocity.

Deanna, Russell G.↗

A tensor train-based isogeometric solver for large-scale 3D poisson problems

We introduce a three-dimensional (3D), fully tensor train (TT) assembled isogeometric analysis (IGA) framework, TT-IGA, for solving partial differential equations (PDEs). Our method reformulates IGA discrete operators into TT format, enabling efficient compression and computation. Geometry evaluations use the original NURBS description at sampling points and TT approximation is applied to geometry-derived coefficient fields and discrete operators. We demonstrate the effectiveness of the proposed TT-IGA framework on the three-dimensional Poisson equation, achieving substantial reductions in memory and computational cost without compromising solution quality.

97 MATHEMATICS AND COMPUTING↗