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

A vectorized Poisson solver over a spherical shell and its application to the quasi-geostrophic omega-equation

The quasi-geostrophic omega-equation in flux form is developed as an example of a Poisson problem over a spherical shell. Solutions of this equation are obtained by applying a two-parameter Chebyshev solver in vector layout for CDC 200 series computers. The performance of this vectorized algorithm greatly exceeds the performance of its scalar analog. The algorithm generates solutions of the omega-equation which are compared with the omega fields calculated with the aid of the mass continuity equation.

Mullenmeister, Paul↗

The equilibrium and stability of the gaseous component of the galaxy, 2

A time-independent, linear, plane and axially-symmetric stability analysis was performed on a self-gravitating, plane-parallel, isothermal layer of nonmagnetic, nonrotating gas. The gas layer was immersed in a plane-stratified field isothermal layer of stars which supply a self-consistent gravitational field. Only the gaseous component was perturbed. Expressions were derived for the perturbed gas potential and perturbed gas density that satisfied both the Poisson and hydrostatic equilibrium equations. The equation governing the size of the perturbations in the mid-plane was found to be analogous to the one-dimensional time-independent Schrodinger equation for a particle bound by a potential well, and with similar boundary conditions. The radius of the neutral state was computed numerically and compared with the Jeans' and Ledoux radius. The inclusion of a rigid stellar component increased the Ledoux radius, though only slightly. Isodensity contours of the neutrual or marginally unstable state were constructed.

Kellman, S. A.↗

Systems improved numerical differencing analyzer

Analyzer is software system that possesses capabilities which make it well suited for solving lumped-parameter representations of physical problems governed by diffusion-type equations, such as Fourier, Poisson, or La Place.

Fink, L. C.↗

Cell volume control at a surface for three-dimensional grid generation packages

An alternate method of calculating the cell size for orthogonality control in the solution of Poisson's 3D space equations is presented. The method provides the capability to enforce a better initial guess for the grid distribution required for boundary layer resolution. This grid point distribution is accomplished by enforcing grid spacing from a grid block boundary where orthogonality is required. The actual grid spacing or cell size for that boundary is determined by the two or four adjacent boundaries in the grid block definition, which are two dimensional grids. These two dimensional grids are in turn defined by the user using insight into the flow field and boundary layer characteristics. The adjoining boundaries are extended using a multifunctional blending scheme, with user control of the blending and interpolating functions to be used. This grid generation procedure results in an enhanced computational fluid dynamics calculation by allowing a quicker resolution of the configuration's boundary layer and flow field and by limiting the number of grid re-adaptations. The cell size specification calculation was applied to a variety of configurations ranging from axisymmetric to complex three-dimensional configurations. Representative grids are shown for the Space Shuttle and the Langley Lifting Body (HL-20).

Alter, Stephen J.↗

An Elongated Tetrakaidecahedron Model for Open-Celled Foams

A micro-mechanics model for non-isotropic, open-celled foams is developed using an elongated tetrakaidecahedron (Kelvin model) as the repeating unit cell. The micro-mechanics model employs an elongated Kelvin model geometry which is more general than that employed by previous authors. Assuming the cell edges possess axial and bending rigidity, the mechanics of deformation of the elongated tetrakaidecahedron lead to a set of equations for the Young's modulus, Poisson's ratio and strength of the foam in the principal material directions. These equations are written as a function of the cell edge lengths and cross-section properties, the inclination angle and the strength and stiffness of the solid material. The model is applied to predict the strength and stiffness of several polymeric foams. Good agreement is observed between the model results and the experimental measurements.

Sullivan, Roy M.↗

A stochastic-dynamic model for global atmospheric mass-field statistics

Global atmospheric mass field error correlations based on satellite observations and on numerical forecasts show strong and systematic latitude dependence. A model for the latitude dependent spatial correlation structure of mass field forecast errors is derived from dynamical considerations. Three methods of solution were tested. In the first method, the equation was solved by expansion in spherical harmonics, and the correlation function was computed analytically using the expansion coefficients. In the second method, the finite difference equivalent of the equation was solved using a fast poisson solver. The correlation function was computed using stratified sampling of the individual realizations. In the third method, a higher order equation was derived, and solved directly in finite differences by two successive applications of the fast poisson solver. The three methods were compared for accuracy and efficiency, and the third method was chosen as clearly superior.

Ghil, M.↗

Computation of solar perturbations with Poisson series

Description of a project for computing first-order perturbations of natural or artificial satellites by integrating the equations of motion on a computer with automatic Poisson series expansions. A basic feature of the method of solution is that the classical variation-of-parameters formulation is used rather than rectangular coordinates. However, the variation-of-parameters formulation uses the three rectangular components of the disturbing force rather than the classical disturbing function, so that there is no problem in expanding the disturbing function in series. Another characteristic of the variation-of-parameters formulation employed is that six rather unusual variables are used in order to avoid singularities at the zero eccentricity and zero (or 90 deg) inclination. The integration process starts by assuming that all the orbit elements present on the right-hand sides of the equations of motion are constants. These right-hand sides are then simple Poisson series which can be obtained with the use of the Bessel expansions of the two-body problem in conjunction with certain interation methods. These Poisson series can then be integrated term by term, and a first-order solution is obtained.

Broucke, R.↗

Self-consistent analysis of lattice-matched and pseudomorphic quantum-well emission transistors

A self-consistent analysis of the quantum-well emission transistor (QWET) is presented allowing an exact calculation of the device quantum properties. Poisson's and Schroedinger's equation are solved numerically using a finite-difference method on a self-consistent basis. Pseudomorphic AlGaAs/InGaAs designs with 15-20 percent excess In are suggested for improving the device performance. Design with doping in various parts of the QWET are also studied. This analysis reveals that the device performance is less optimistic than previously predicted by analytic approaches. By introducing the pseudomorphic channel principle, while maintaining a reasonably low Al content for the gate and collector layers, it is, however, possible to obtain satisfactory performance. Optimum pseudomorphic designs showed high current driving capability (200,000 A/sq cm), high transconductance (3S/mm), and small intrinsic delay time (2 ps).

Hong, Kyushik↗

An improved plate theory of order (1,2) for thick composite laminates

A new (1,2)-order theory is proposed for the linear elasto-static analysis of laminated composite plates. The basic assumptions are those concerning the distribution through the laminate thickness of the displacements, transverse shear strains and the transverse normal stress, with these quantities regarded as some weighted averages of their exact elasticity theory representations. The displacement expansions are linear for the inplane components and quadratic for the transverse component, whereas the transverse shear strains and transverse normal stress are respectively quadratic and cubic through the thickness. The main distinguishing feature of the theory is that all strain and stress components are expressed in terms of the assumed displacements prior to the application of a variational principle. This is accomplished by an a priori least-square compatibility requirement for the transverse strains and by requiring exact stress boundary conditions at the top and bottom plate surfaces. Equations of equilibrium and associated Poisson boundary conditions are derived from the virtual work principle. It is shown that the theory is particularly suited for finite element discretization as it requires simple C(sup 0)- and C(sup -1)-continuous displacement interpolation fields. Analytic solutions for the problem of cylindrical bending are derived and compared with the exact elasticity solutions and those of our earlier (1,2)-order theory based on the assumed displacements and transverse strains.

Tessler, A.↗

Parallel Unsteady Turbopump Flow Simulations for Reusable Launch Vehicles

An efficient solution procedure for time-accurate solutions of Incompressible Navier-Stokes equation is obtained. Artificial compressibility method requires a fast convergence scheme. Pressure projection method is efficient when small time-step is required. The number of sub-iteration is reduced significantly when Poisson solver employed with the continuity equation. Both computing time and memory usage are reduced (at least 3 times). Other work includes Multi Level Parallelism (MLP) of INS3D, overset connectivity for the validation case, experimental measurements, and computational model for boost pump.

Kiris, Cetin↗

Computations of Wall Distances Based on Differential Equations

The use of differential equations such as Eikonal, Hamilton-Jacobi and Poisson for the economical calculation of the nearest wall distance d, which is needed by some turbulence models, is explored. Modifications that could palliate some turbulence-modeling anomalies are also discussed. Economy is of especial value for deforming/adaptive grid problems. For these, ideally, d is repeatedly computed. It is shown that the Eikonal and Hamilton-Jacobi equations can be easy to implement when written in implicit (or iterated) advection and advection-diffusion equation analogous forms, respectively. These, like the Poisson Laplacian term, are commonly occurring in CFD solvers, allowing the re-use of efficient algorithms and code components. The use of the NASA CFL3D CFD program to solve the implicit Eikonal and Hamilton-Jacobi equations is explored. The re-formulated d equations are easy to implement, and are found to have robust convergence. For accurate Eikonal solutions, upwind metric differences are required. The Poisson approach is also found effective, and easiest to implement. Modified distances are not found to affect global outputs such as lift and drag significantly, at least in common situations such as airfoil flows.

Tucker, Paul G.↗

Generic Karman-Rostovstev plate equations in an affine space

The Karman-Rostovstev plate equations are affinely transformed into equations that depend on the generic plate constants generalized rigidity, generalized extensional rigidity, and generalized Poisson ratio associated with the well-known linearized orthotropic plate/slab equations. Furthermore, since the generalized extensional rigidity is a function of generalized rigidity and generalized Poisson ratio only, with generalized Poisson ratio being a weak parameter, the only fundamental constant is generalized rigidity. This means that regardless of what analysis technique is used, the results can be cast in terms of only one strong material parameter (whose range is from zero to one). The affine stretching is performed on the dependent variables as well as the independent variables, so that the dependence of all quantities of interest on generalized rigidity and generalized Poisson ratio are easily calculated (e.g., even auxiliary quantities such as stress, strain, all weighted stress resultants, etc.). Some typical results (static side load response, collapse load limited by rib buckling) are presented by modifying an approximate technique due to Donnell (1976).

Brunelle, E. J.↗

Special perturbations employing osculating reference states

The concept of employing osculating reference position and velocity vectors in the numerical integration of the equations of motion of a satellite is examined. The choice of the reference point is shown to have a significant effect upon numerical efficiency and the class of trajectories described by the differential equations of motion. For example, when the position and velocity vectors on the osculating orbit at a fixed reference time are chosen, a universal formulation is yielded. For elliptical orbits, however, this formulation is unattractive for numerical integration purposes due to Poisson terms (mixed secular) appearing in the equations of motion. Other choices for the reference point eliminate this problem but usually at the expense of universality. A number of these formulations, including a universal one, are considered here. Comparisons of the numerical characteristics of these techniques with those of the Encke method are presented.

Born, G. H.↗

Viscoplastic Characterization of Ti-6-4: Experiments

As part of a continued effort to improve the understanding of material time-dependent response, a series of mechanical tests have been conducted on the titanium alloy, Ti-6Al-4V. Tensile, creep, and stress relaxation tests were performed over a wide range of temperatures and strain rates to engage various amounts of time-dependent behavior. Additional tests were conducted that involved loading steps, overloads, dwell periods, and block loading segments to characterize the interaction between plasticity and time-dependent behavior. These data will be used to characterize a recently developed, viscoelastoplastic constitutive model with a goal toward better estimates of aerospace component behavior, resulting in improved safety.

mechanical properties↗

Nonlinear modes of the tensor Dirac equation and CPT violation

Recently, it has been shown that Dirac's bispinor equation can be expressed, in an equivalent tensor form, as a constrained Yang-Mills equation in the limit of an infinitely large coupling constant. It was also shown that the free tensor Dirac equation is a completely integrable Hamiltonian system with Lie algebra type Poisson brackets, from which Fermi quantization can be derived directly without using bispinors. The Yang-Mills equation for a finite coupling constant is investigated. It is shown that the nonlinear Yang-Mills equation has exact plane wave solutions in one-to-one correspondence with the plane wave solutions of Dirac's bispinor equation. The theory of nonlinear dispersive waves is applied to establish the existence of wave packets. The CPT violation of these nonlinear wave packets, which could lead to new observable effects consistent with current experimental bounds, is investigated.

Reifler, Frank J.↗

Particle localization, spinor two-valuedness, and Fermi quantization of tensor systems

Recent studies of particle localization shows that square-integrable positive energy bispinor fields in a Minkowski space-time cannot be physically distinguished from constrained tensor fields. In this paper we generalize this result by characterizing all classical tensor systems, which admit Fermi quantization, as those having unitary Lie-Poisson brackets. Examples include Euler's tensor equation for a rigid body and Dirac's equation in tensor form.

Reifler, Frank↗

Calculation of sheath and wake structure about a pillbox-shaped spacecraft in a flowing plasma

A computer program was used for studies of the disturbed zones around bodies in flowing plasmas, particularly spacecraft and their associated sheaths and wakes. The program solved a coupled Poisson-Vlasov system of nonlinear partial differential integral equations to obtain distributions of electric potential and ion and electron density about a finite length cylinder in a plasma flow at arbitrary ion Mach numbers. The approach was applicable to a larger range of parameters than other available approaches. In sample calculations, bodies up to 100 Debye lengths in radius were treated, that is, larger than any previously treated realistically. Applications were made to in-situ satellite experiments.

Parker, L. W.↗