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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 685 records · Page 38

All sources NASA NTRS DOE OSTI

NASA NTRS1982

Algebraic grid generation

The numerical solution of partial differential equations about irregular geometries and with varying characteristic scales has created the need for coordinate systems and associated transformations which reflect both geometric and physical requirements. The process of finding coordinate transformations in discrete representations is called 'grid generation'. The present investigation is concerned with three algebraic grid generation methods. The methods include transfinite interpolation, the multisurface method, and the two-boundary technique. Interpolation formulas, in terms of homotopic mappings and constraints in terms of point positions and/or derivatives are the essential elements of the techniques. The methods are relatively simple to understand, they are explicit and do not require extensive computational effort, and they have a high degree of generality.

Smith, R. E.↗

NASA NTRS1983

Calculation of the polarization potential for e-N2 collisions

A polarization potential V(pol) for e-N2 collisions is calculated by the generalization of the static part of the method of polarized orbitals to molecular targets. Partial differential equations (PDE) are derived for polarized orbitals, which are functions of the distance (r) from the molecular center and angle (theta) from the molecular axis. The equations are solved with the use of the noniterative PDE technique. From the polarized orbitals a polarization potential can be constructed whose r and theta dependence is found to be significantly different from the well-known phenomenological one. The two potentials are further compared by carrying out limited hybrid-theory scattering calculations. Only those scattering results based on the calculated V(pol) are in satisfactory accord with experiment.

Onda, K.↗

NASA NTRS1982

Flow structure in continuous flow electrophoresis chambers

There are at least two ways that hydrodynamic processes can limit continiuous flow electrophoresis. One arises from the sensitivity of the flow to small temerature gradients, especially at low flow rates and power levels. This sensitivity can be suppressed, at least in principle, by providing a carefully tailored, stabilizing temperature gradient in the cooling system that surrounds the flow channel. At higher power levels another limitation arises due to a restructuring of the main flow. This restructuring is caused by buoyancy, which is in turn affected by the electro-osmotic crossflow. Approximate solutions to appropriate partial differential equations have been computed by finite difference methods. One set of results is described here to illustrate the strong coupling between the structure of the main (axial) flow and the electro-osmotic flow.

Deiber, J. A.↗

NASA NTRS1983

Cubic spline approximation techniques for parameter estimation in distributed systems

Approximation schemes employing cubic splines in the context of a linear semigroup framework are developed for both parabolic and hyperbolic second-order partial differential equation parameter estimation problems. Convergence results are established for problems with linear and nonlinear systems, and a summary of numerical experiments with the techniques proposed is given.

Banks, H. T.↗

NASA NTRS1983

Finite-dimensional discrete-time control of linear distributed parameter systems

Systems described by partial differential equations have an infinite-dimensional state space. Feedback control of such distributed parameter systems must be accomplished by finite-dimensional controllers to be implemented by on-line digital computers. A further practical constraint is that the controller must operate in discrete (rather than continuous) time. This paper investigates the stability of such distributed parameter feedback controllers in closed-loop with the actual system.

Balas, M. J.↗

NASA NTRS1983

A lifting surface theory in rotational flow

The partial differential equation for small disturbance steady rotational flow in three dimensions is solved through an integral equation approach. The solution is obtained by using the method of weighted residuals. Specific applications are directed to wings in nonuniform subsonic parallel streams with velocity varying in vertical and spanwise directions and to airfoils in nonuniform freestream. Comparison with limited known results indicates that the present method is reasonably accurate. Numerical results for the lifting pressure of airfoil, lift, induced drag, and pitching moments of airfoil, lift, induced drag, and pitching moments of elliptic, rectangular, and delta wings in a jet, wake, or monotonic sheared stream are presented. It is shown that, in addition to the effect of local dynamic pressures, a positive velocity gradient tends to enhance the lift.

Shiau, M. J.↗

NASA NTRS1983

The theory of an active magnetic regenerative refrigerator

The adiabatic temperature change with field which is limited to about 2 K/Tesla for ferromagnets near their Curie temperatures by the change of magnetization with temperature and the lattice heat capacity is discussed. Practical magnetic refrigerators operate on a regenerative cycle such as the Brayton cycle. This cycle can be executed through the use of an active magnetic regenerator, i.e., a regenerator composed of magnetic material that is cycled in an out of a magnetic field with appropriate fluid flows. The theory of these devices is predicted by solving the partial differential equations that describe fluid and the magnetic solid. The active magnetic regenerator is described along with the method of calculation. Temperature profiles for a normal regenerator and a magnetic regenerative refrigerator are shown.

Barclay, J. A.↗

NASA NTRS1984

Mathematical algorithms to maximize performance in numerical weather prediction

Numerical weather prediction models, which involve the solution of non-linear partial differential equations at points on an extensive three dimensional grid, are ideally suited for processing on vector machines. It was logical therefore that the new global forecast model to be implemented at the Meteorological Office should be written in vector code for the CYBER 205. In order to achieve full efficiency and to reduce storage requirements the model used 32-bit arithmetic which was found to provide high enough precision. Unfortunately, however, the trigonometrical and logarithmic functions provided by CDC could only handle 64-bit vectors and, although written in efficient scalar code, did not take advantage of the special facilities of a vector processor. It was therefore necessary to rewrite the functions in vector code to handle both 32 and 64-bit vectors. There was also no half-precision compiler available for the Cyber 205 at that time and so the functions, like the model, had to make extensive use of the special call syntax. This made the code more difficult to write but it allowed much greater flexibility in that it became possible to access the exponent of a floating-point number independently of its coefficient. A description is given of the technique and the results which were achieved are summarized.

Foreman, A.↗

NASA NTRS1984

Analytical results for postbuckling behavior of plates in compression and shear

The postbuckling behavior of long rectangular isotropic and orthotropic plates is determined. By assuming trigonometric functions in one direction, the nonlinear partial differential equations of von Karman large deflection plate theory are converted into nonlinear ordinary differential equations. The ordinary differential equations are solved numerically using an available boundary value problem solver which makes use of Newton's method. Results for longitudinal compression show different postbuckling behavior between isotropic and orthotropic plates. Results for shear show that change in inplane edge constraints can cause large change in postbuckling stiffness.

Stein, M.↗

NASA NTRS1984

Space structure vibration modes: How many exist? Which ones are important?

This report attempts to shed some light on the two issues raised in the title, namely, how many vibration modes does a real structure have, and which of these modes are important? The surprise-free answers to these two questions are, respectively, an infinite number and the first several modes. The author argues that the absurd subspace (all but the first billion modes) is not a strength of continuum modeling, but, in fact, a weakness. Partial differential equations are not real structures, only mathematical models. This note also explains (1) that the PDE model and the finite element model are, in fact, the same model, the latter being a numerical method for dealing with the former, (2) that modes may be selected on dynamical grounds other than frequency alone, and (3) that long slender rods are useful as primitive cases but dangerous to extrapolate from.

Hughes, P. C.↗

NASA NTRS1984

Analysis of lateral and torsional vibration characteristics of beams and shafts with end located rotational masses

Partial differential equations are derived for free lateral and torsional vibration of a uniform free-free beam with a rotational mass attached to each extremity. For appropriate boundary conditions, nonlinear algebraic equations are obtained using a symbolic manipulation computer program, the solutions of which enable the computation of the neutral frequencies and mode-shapes. The mode-shapes are linear combinations of trigonometric and hyperbolic sine and cosine functions. A computer program is written for the numerical solution of the algebraic equations mentioned above, which can compute the natural frequencies, mode-shapes, and node points for any given set of parameters, for any given number of modes.

Robertson, D. K.↗

NASA NTRS1979

Solution of 3-dimensional time-dependent viscous flows. Part 1: Investigation of candidate algoriths

There is considerable interest in developing a numerical scheme for solving the time dependent viscous compressible three dimensional flow to aid in the design of helicopter rotors. Numerical algorithms are examined to determine their overall suitability for the efficient and routine solution of an appropriate system of partial differential equations. It is concluded that a consistently split time linearized block implicit scheme using either quintic B spline collocation or the generalized operator compact implicit approach to generate a fourth order accurate algorithm is particularly well suited for use on the present problem. High cell Reynolds number behavior leads to favoring the generalized operator compact implicit approach over the quintic B spline collocation method.

Weinberg, B. C.↗

NASA NTRS1981

PAN AIR: A computer program for predicting subsonic or supersonic linear potential flows about arbitrary configurations using a higher order panel method. Volume 1: Theory document (version 1.1)

Panel aerodynamics (PAN AIR) is a system of computer programs designed to analyze subsonic and supersonic inviscid flows about arbitrary configurations. A panel method is a program which solves a linear partial differential equation by approximating the configuration surface by a set of panels. An overview of the theory of potential flow in general and PAN AIR in particular is given along with detailed mathematical formulations. Fluid dynamics, the Navier-Stokes equation, and the theory of panel methods were also discussed.

Magnus, A. E.↗

NASA NTRS1984

Three-dimensional transonic potential flow about complex 3-dimensional configurations

An analysis has been developed and a computer code written to predict three-dimensional subsonic or transonic potential flow fields about lifting or nonlifting configurations. Possible condfigurations include inlets, nacelles, nacelles with ground planes, S-ducts, turboprop nacelles, wings, and wing-pylon-nacelle combinations. The solution of the full partial differential equation for compressible potential flow written in terms of a velocity potential is obtained using finite differences, line relaxation, and multigrid. The analysis uses either a cylindrical or Cartesian coordinate system. The computational mesh is not body fitted. The analysis has been programmed in FORTRAN for both the CDC CYBER 203 and the CRAY-1 computers. Comparisons of computed results with experimental measurement are presented. Descriptions of the program input and output formats are included.

Reyhner, T. A.↗

NASA NTRS1984

Two-dimensional computer simulation of EMVJ and grating solar cells under AMO illumination

A computer program, SCAP2D (Solar Cell Analysis Program in 2-Dimensions), is used to evaluate the Etched Multiple Vertical Junction (EMVJ) and grating solar cells. The aim is to demonstrate how SCAP2D can be used to evaluate cell designs. The cell designs studied are by no means optimal designs. The SCAP2D program solves the three coupled, nonlinear partial differential equations, Poisson's Equation and the hole and electron continuity equations, simultaneously in two-dimensions using finite differences to discretize the equations and Newton's Method to linearize them. The variables solved for are the electrostatic potential and the hole and electron concentrations. Each linear system of equations is solved directly by Gaussian Elimination. Convergence of the Newton Iteration is assumed when the largest correction to the electrostatic potential or hole or electron quasi-potential is less than some predetermined error. A typical problem involves 2000 nodes with a Jacobi matrix of order 6000 and a bandwidth of 243.

Gray, J. L.↗

NASA NTRS1984

A numerical scheme for the identification of hybrid systems describing the vibration of flexible beams with tip bodies

A cubic spline based Galerkin-like method is developed for the identification of a class of hybrid systems which describe the transverse vibration to flexible beams with attached tip bodies. The identification problem is formulated as a least squares fit to data subject to the system dynamics given by a coupled system of ordnary and partial differential equations recast as an abstract evolution equation (AEE) in an appropriate infinite dimensional Hilbert space. Projecting the AEE into spline-based subspaces leads naturally to a sequence of approximating finite dimensional identification problems. The solutions to these problems are shown to exist, are relatively easily computed, and are shown to, in some sense, converge to solutions to the original identification problem. Numerical results for a variety of examples are discussed.

Rosen, I. G.↗

NASA NTRS1984

The influence of a wind tunnel on helicopter rotational noise: Formulation of analysis

An analytical model is discussed that can be used to examine the effects of wind tunnel walls on helicopter rotational noise. A complete physical model of an acoustic source in a wind tunnel is described and a simplified version is then developed. This simplified model retains the important physical processes involved, yet it is more amenable to analysis. The simplified physical model is then modeled as a mathematical problem. An inhomogeneous partial differential equation with mixed boundary conditions is set up and then transformed into an integral equation. Details of generating a suitable Green's function and integral equation are included and the equation is discussed and also given for a two-dimensional case.

Mosher, M.↗

NASA NTRS1983

Implicit, nonswitching, vector-oriented algorithm for steady transonic flow

A rapid computation of a sequence of transonic flow solutions has to be performed in many areas of aerodynamic technology. The employment of low-cost vector array processors makes the conduction of such calculations economically feasible. However, for a full utilization of the new hardware, the developed algorithms must take advantage of the special characteristics of the vector array processor. The present investigation has the objective to develop an efficient algorithm for solving transonic flow problems governed by mixed partial differential equations on an array processor.

Lottati, I.↗