Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “partial differential 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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
NASA NTRS1987

Algebraic multigrid

The state of the art in algebraic multgrid (AMG) methods is discussed. The interaction between the relaxation process and the coarse grid correction necessary for proper behavior of the solution probes is discussed in detail. Sufficient conditions on relaxation and interpolation for the convergence of the V-cycle are given. The relaxation used in AMG, what smoothing means in an algebraic setting, and how it relates to the existing theory are considered. Some properties of the coarse grid operator are discussed, and results on the convergence of two-level and multilevel convergence are given. Details of an algorithm particularly studied for problems obtained by discretizing a single elliptic, second order partial differential equation are given. Results of experiments with such problems using both finite difference and finite element discretizations are presented.

Ruge, J. W.↗
NASA NTRS1988

A comparative study of theoretical methods on Goertler instability

Goertler vortices arise in boundary layers along concave surfaces due to centrifugal effects and these vortices in combination with Tollmien-Schlichting (TS) waves and crossflow vortices may play an important role in triggering early transition. There have been two distinct theoretical approaches to solve the problem. In the classical, normal-mode approach (NMA) the linear development of Goertler vortices is reduced to an eigenvalue problem. In the other method, initiated by Hall, the governing partial differential equations for the perturbations which are parabolic in the streamwise direction, are solved as an initial value problem. This method predicts multiple neutral curves depending on how and where the mean flow is perturbed, unlike the NMA where a unique but different neutral curve is predicted by each investigator. The present paper attempts to compare the two techniques and shows that mutually compatible results can be obtained when physically consistent assumptions are made.

Kalburgi, Vijay↗
NASA NTRS1988

The computation of hypersonic ionized flows in chemical and thermal nonequlibrium

A numerical method to compute a two-dimensional hypersonic flowfield that is ionized and in thermochemical nonequilibrium has been developed. Such a flowfield is described by coupled time-dependent partial differential equations for the conservation of species mass, mass-average momentum, vibrational energy of each diatomic species, electron energy, and total mass-averaged energy. The steady-state solution to these fully coupled equations is obtained using an implicit Gauss-Seidel line relaxation technique. The computed electron densities in the flowfield compare well with experimental results.

Maccormack, Robert W.↗
NASA NTRS1987

Aeroelastic modeling of swept tip rotor blades using finite elements

A special finite element for the modeling of a swept tip rotor blade is derived. The swept tip blade undergoes moderate deflections in flap, lag and torsion. The nonlinear, partial differential equations of motion are discretized using a Galerkin finite element method. Tip sweep introduces flap-torsion and lag-axial couplings, and may lead to aeroelastic instabilities associated with frequency coalescence. When frequency coalescence does not occur, sweep is usually stabilizing. A comparison of the approximate model of a swept tip blade, using a blade with straight elastic axis and offsets of aerodynamic centers of gravity, with the exact swept tip model developed in this paper indicates that the approximate model can produce inaccurate results for the case of hingeless rotor blades.

Celi, R.↗
NASA NTRS1987

Theories of multiple equilibria and weather regimes : A critical reexamination. II - Baroclinic two-layer models

Previous results based on low- and intermediate-order truncations of the two-layer model suggest the existence of multiple equilibria and/or multiple weather regimes for the extratropical large-scale flow. The importance of the transient waves in the synoptic scales in organizing the large-scale flow and in the maintenance of weather regimes was emphasized. The result shows that multiple equilibria/weather regimes that are present in lower-order models examined disappear when a sufficient number of modes are kept in the spectral expansion of the solution to the governing partial differential equations. Much of the chaotic behavior of the large-scale flow that is present in intermediate-order models is now found to be spurious. Physical reasons for the drastic modification are offered. A peculiarity in the formulation of most existing two-layer models is noted that also tends to exaggerate the importance of baroclinic processes and increase the degree of unpredictability of the large-scale flow.

Cehelsky, Priscilla↗
NASA NTRS1988

Far field boundary conditions for computation over long time

A new class of computational far field boundary conditions for hyperbolic partial differential equations is developed. These boundary conditions combine properties of absorbing boundary conditions for transient solutions and properties of far field boundary conditions for steady-state problems. The conditions can be used to limit the computational domain when both traveling waves and evanescent waves are present. Boundary conditions for scalar wave equations are derived and analyzed. Extensions to systems of equations are discussed and results from numerical experiments are presented.

Engquist, Bjorn↗
NASA NTRS1987

Adaptive grid generation

The fundamental principles of adaptive grid generation for the numerical analysis of physical phenomena described by systems of partial differential equations are examined in an analytical review. Topics addressed include weight functions, equidistribution in one dimension, the specification of coefficients in the linear weight, the attraction to a given grid on a curve, evolutionary forces, and metric notation. Consideration is given to curve-by-curve methods, finite-volume methods, variational methods, and temporal aspects.

Eiseman, Peter R.↗
NASA NTRS1988

A new linearized theory of laminar film condensation of two phase annular flow in a capillary pumped loop

Future large space based facilities, such as Space Station, will require energy management systems capable of transporting tens of kilowatts of heat over a hundred meters or more. This represents better than an order of magnitude improvement over current technology. Two-phase thermal systems are currently being developed to meet this challenge. Condensation heat transfer plays a very important role in this system. The present study attempts an analytic solution to the set of linearized partial differential equations. The axial velocity and temperature functions were found to be Bessel functions which have oscillatory behavior. This result agrees qualitatively with the experimental evidence from tests at both NASA Goddard Space Flight Center and elsewhere.

Hsu, Y. K.↗
NASA NTRS1988

The identification of a distributed parameter model for a flexible structure

A computational method is developed for the estimation of parameters in a distributed model for a flexible structure. The structure we consider (part of the RPL experiment) consists of a cantilevered beam with a thruster and linear accelerometer at the free end. The thruster is fed by a pressurized hose whose horizontal motion effects the transverse vibration of the beam. The Euler-Bernoulli theory is used to model the vibration of the beam and treat the hose thruster assembly as a lumped or point mass dashpot spring system at the tip. Measurements of linear acceleration at the tip are used to estimate the hose parameters (mass, stiffness, damping) and a Voigt-Kelvin viscoelastic structural damping parameter for the beam using a least squares fit to the data. Spline based approximations are considered to the hybrid (coupled ordinary and partial differential equations) systems; theoretical convergence results and numerical studies with both simulation and actual experimental data obtained from the structure are presented and discussed.

Banks, H. T.↗
NASA NTRS1988

Goertler vortices in supersonic boundary layers

The problem of Goertler vortices in compressible boundary layers over concave walls is studied by solving parabolic partial differential equations. The numerical scheme used is based upon a fourth-order accurate compact difference scheme. The results show that both the velocity and temperature fluctuations may lead to a Goertler vortex. The vortex growth rates determined from the present method are found to differ somewhat from those given by normal mode solution. The effects of cooling, heating and pressure gradient are also studied. At a Mach number of 3.5, cooling has a small destabilizing effect. The effect of adverse pressure gradient is also found to be destabilizing on a Mach 8 boundary layer.

Spall, R. E.↗
NASA NTRS1988

Problem size, parallel architecture, and optimal speedup

The communication and synchronization overhead inherent in parallel processing can lead to situations where adding processors to the solution method actually increases execution time. Problem type, problem size, and architecture type all affect the optimal number of processors to employ. The numerical solution of an elliptic partial differential equation is examined in order to study the relationship between problem size and architecture. The equation's domain is discretized into n sup 2 grid points which are divided into partitions and mapped onto the individual processor memories. The relationships between grid size, stencil type, partitioning strategy, processor execution time, and communication network type are analytically quantified. In so doing, the optimal number of processors was determined to assign to the solution, and identified (1) the smallest grid size which fully benefits from using all available processors, (2) the leverage on performance given by increasing processor speed or communication network speed, and (3) the suitability of various architectures for large numerical problems.

Nicol, David M.↗
NASA NTRS1989

Exponential Finite-Difference Technique

Report discusses use of explicit exponential finite-difference technique to solve various diffusion-type partial differential equations. Study extends technique to transient-heat-transfer problems in one dimensional cylindrical coordinates and two and three dimensional Cartesian coordinates and to some nonlinear problems in one or two Cartesian coordinates.

Handschuh, Robert F.↗
NASA NTRS1988

Domain decomposition methods for the parallel computation of reacting flows

Domain decomposition is a natural route to parallel computing for partial differential equation solvers. Subdomains of which the original domain of definition is comprised are assigned to independent processors at the price of periodic coordination between processors to compute global parameters and maintain the requisite degree of continuity of the solution at the subdomain interfaces. In the domain-decomposed solution of steady multidimensional systems of PDEs by finite difference methods using a pseudo-transient version of Newton iteration, the only portion of the computation which generally stands in the way of efficient parallelization is the solution of the large, sparse linear systems arising at each Newton step. For some Jacobian matrices drawn from an actual two-dimensional reacting flow problem, comparisons are made between relaxation-based linear solvers and also preconditioned iterative methods of Conjugate Gradient and Chebyshev type, focusing attention on both iteration count and global inner product count. The generalized minimum residual method with block-ILU preconditioning is judged the best serial method among those considered, and parallel numerical experiments on the Encore Multimax demonstrate for it approximately 10-fold speedup on 16 processors.

Keyes, David E.↗
NASA NTRS1988

Nonconforming mortar element methods: Application to spectral discretizations

Spectral element methods are p-type weighted residual techniques for partial differential equations that combine the generality of finite element methods with the accuracy of spectral methods. Presented here is a new nonconforming discretization which greatly improves the flexibility of the spectral element approach as regards automatic mesh generation and non-propagating local mesh refinement. The method is based on the introduction of an auxiliary mortar trace space, and constitutes a new approach to discretization-driven domain decomposition characterized by a clean decoupling of the local, structure-preserving residual evaluations and the transmission of boundary and continuity conditions. The flexibility of the mortar method is illustrated by several nonconforming adaptive Navier-Stokes calculations in complex geometry.

Maday, Yvon↗
NASA NTRS1989

Role of computational fluid dynamics in unsteady aerodynamics for aeroelasticity

In the last two decades there have been extensive developments in computational unsteady transonic aerodynamics. Such developments are essential since the transonic regime plays an important role in the design of modern aircraft. Therefore, there has been a large effort to develop computational tools with which to accurately perform flutter analysis at transonic speeds. In the area of Computational Fluid Dynamics (CFD), unsteady transonic aerodynamics are characterized by the feature of modeling the motion of shock waves over aerodynamic bodies, such as wings. This modeling requires the solution of nonlinear partial differential equations. Most advanced codes such as XTRAN3S use the transonic small perturbation equation. Currently, XTRAN3S is being used for generic research in unsteady aerodynamics and aeroelasticity of almost full aircraft configurations. Use of Euler/Navier Stokes equations for simple typical sections has just begun. A brief history of the development of CFD for aeroelastic applications is summarized. The development of unsteady transonic aerodynamics and aeroelasticity are also summarized.

Guruswamy, Guru P.↗
NASA NTRS1989

Run-time scheduling and execution of loops on message passing machines

Sparse system solvers and general purpose codes for solving partial differential equations are examples of the many types of problems whose irregularity can result in poor performance on distributed memory machines. Often, the data structures used in these problems are very flexible. Crucial details concerning loop dependences are encoded in these structures rather than being explicitly represented in the program. Good methods for parallelizing and partitioning these types of problems require assignment of computations in rather arbitrary ways. Naive implementations of programs on distributed memory machines requiring general loop partitions can be extremely inefficient. Instead, the scheduling mechanism needs to capture the data reference patterns of the loops in order to partition the problem. First, the indices assigned to each processor must be locally numbered. Next, it is necessary to precompute what information is needed by each processor at various points in the computation. The precomputed information is then used to generate an execution template designed to carry out the computation, communication, and partitioning of data, in an optimized manner. The design is presented for a general preprocessor and schedule executer, the structures of which do not vary, even though the details of the computation and of the type of information are problem dependent.

Crowley, Kay↗
NASA NTRS1988

The physics of parallel machines

The idea is considered that architectures for massively parallel computers must be designed to go beyond supporting a particular class of algorithms to supporting the underlying physical processes being modelled. Physical processes modelled by partial differential equations (PDEs) are discussed. Also discussed is the idea that an efficient architecture must go beyond nearest neighbor mesh interconnections and support global and hierarchical communications.

Chan, Tony F.↗
NASA NTRS1989

On damping mechanisms in beams

A partial differential equation model of a cantilevered beam with a tip mass at its free end is used to study damping in a composite. Four separate damping mechanisms consisting of air damping, strain rate damping, spatial hysteresis and time hysteresis are considered experimentally. Dynamic tests were performed to produce time histories. The time history data is then used along with an approximate model to form a sequence of least squares problems. The solution of the least squares problem yields the estimated damping coefficients. The resulting experimentally determined analytical model is compared with the time histories via numerical simulation of the dynamic response. The procedure suggested here is compared with a standard modal damping ratio model commonly used in experimental modal analysis.

Banks, H. T.↗