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 901 records · Page 50

All sources NASA NTRS DOE OSTI

NASA NTRS1993

Nonlinear dynamic simulation of single- and multi-spool core engines

In this paper a new computational method for accurate simulation of the nonlinear dynamic behavior of single- and multi-spool core engines, turbofan engines, and power generation gas turbine engines is presented. In order to perform the simulation, a modularly structured computer code has been developed which includes individual mathematical modules representing various engine components. The generic structure of the code enables the dynamic simulation of arbitrary engine configurations ranging from single-spool thrust generation to multi-spool thrust/power generation engines under adverse dynamic operating conditions. For precise simulation of turbine and compressor components, row-by-row calculation procedures were implemented that account for the specific turbine and compressor cascade and blade geometry and characteristics. The dynamic behavior of the subject engine is calculated by solving a number of systems of partial differential equations, which describe the unsteady behavior of the individual components. In order to ensure the capability, accuracy, robustness, and reliability of the code, comprehensive critical performance assessment and validation tests were performed. As representatives, three different transient cases with single- and multi-spool thrust and power generation engines were simulated. The transient cases range from operating with a prescribed fuel schedule, to extreme load changes, to generator and turbine shut down.

Schobeiri, T.↗

NASA NTRS1950

Linearized compressible-flow theory for sonic flight speeds

The partial differential equation for the perturbation velocity potential is examined for free-stream Mach numbers close to and equal to one. It is found that, under the assumptions of linearized theory, solutions can be found consistent with the theory for lifting-surface problems both in stationary three-dimensional flow and in unsteady two-dimensional flow. Several examples are solved including a three dimensional swept-back wing and two dimensional harmonically-oscillating wing, both for a free stream Mach number equal to one. Momentum relations for the evaluation of wave and vortex drag are also discussed. (author)

Heaslet, Max A↗

NASA NTRS1950

Use of source distributions for evaluating theoretical aerodynamics of thin finite wings at supersonic speeds

A series of publications on the source-distribution methods for evaluating the aerodynamics of thin wings at supersonic speeds is summarized, extended, and unified. Included in the first part are the deviations of: (a) the linearized partial-differential equation for unsteady flow at a substantially constant Mach number. b) The source-distribution solution for the perturbation-velocity potential that satisfies the boundary conditions of tangential flow at the surface and in the plane of the wing; and (c) the integral equation for determining the strength and the location of sources to describe the interaction effects (as represented by upwash) of the bottom and top wing surfaces through the region between the finite wing boundary and the foremost Mach wave. The second part deals with steady-state thin-wing problems. The third part of the report approximates the integral equation for unsteady upwash and includes a solution of approximate equation. Expressions are then derived to evaluate the load distributions for time-dependent finite-wing motions.

Evvard, John C↗

NASA NTRS1950

The application of Green's theorem to the solution of boundary-value problems in linearized supersonic wing theory

Following the introduction of the linearized partial differential equation for nonsteady three-dimensional compressible flow, general methods of solution are given for the two and three-dimensional steady-state and two-dimensional unsteady-state equations. It is also pointed out that, in the absence of thickness effects, linear theory yields solutions consistent with the assumptions made when applied to lifting-surface problems for swept-back plan forms at sonic speeds. The solutions of the particular equations are determined in all cases by means of Green's theorem, and thus depend on the use of Green's equivalent layer of sources, sinks, and doublets. Improper integrals in the supersonic theory are treated by means of Hadamard's "finite part" technique.

Heaslet, Max A↗

NASA NTRS1953

A Vector Study of Linearized Supersonic Flow Applications to Nonplanar Problems

A vector study of the partial-differential equation of steady linearized supersonic flow is presented. General expressions which relate the velocity potential in the stream to the conditions on the disturbing surfaces, are derived. In connection with these general expressions the concept of the finite part of an integral is discussed. A discussion of problems dealing with planar bodies is given and the conditions for the solution to be unique are investigated. Problems concerning nonplanar systems are investigated, and methods are derived for the solution of some simple nonplanar bodies. The surface pressure distribution and the damping in roll are found for rolling tails consisting of four, six, and eight rectangular fins for the Mach number range where the region of interference between adjacent fins does not affect the fin tips.

Martin, John C↗

NASA NTRS1993

The Chebyshev-Legendre method: Implementing Legendre methods on Chebyshev points

We present a new collocation method for the numerical solution of partial differential equations. This method uses the Chebyshev collocation points, but because of the way the boundary conditions are implemented, it has all the advantages of the Legendre methods. In particular, L2 estimates can be obtained easily for hyperbolic and parabolic problems.

Don, Wai Sun↗

NASA NTRS1993

Sloshing dynamics on rotating helium dewar tank

The generalized mathematical formulation of sloshing dynamics for partially filled liquid of cryogenic superfluid helium II in dewar containers driven by both the gravity gradient and jitter accelerations applicable to scientific spacecraft which is eligible to carry out spinning motion and/or slew motion for the purpose to perform scientific observation during the normal spacecraft operation are investigated. An example is given with Gravity Probe-B (GP-B) spacecraft which is responsible for the sloshing dynamics. The jitter accelerations include slew motion, spinning motion, atmospheric drag on the spacecraft, spacecraft attitude motions arising from machinery vibrations, thruster firing, pointing control of spacecraft, crew motion, etc. Explicit mathematical expressions to cover these forces acting on the spacecraft fluid systems are derived. The numerical computation of sloshing dynamics were based on the non-inertia frame spacecraft bound coordinate, and solve time dependent, three-dimensional formulations of partial differential equations subject to initial and boundary conditions. The explicit mathematical expressions of boundary conditions to cover capillary force effect on the liquid vapor interface in microgravity environments are also derived. The formulations of fluid moment and angular moment fluctuations in fluid profiles induced by the sloshing dynamics, together with fluid stress and moment fluctuations exerted on the spacecraft dewar containers were derived. Results were widely published in the open journals.

Hung, R. J.↗

NASA NTRS1992

Comparison of level-1 and level-2 buckling and postbuckling solutions

The accuracy of the solutions obtained with the Level-1 computational modules of DISDECO (Delft Interactive Shell Design Code) are investigated by comparing them with solutions obtained with ANILISA, a Level-2 computational module. Whereas the Level-1 solutions are based on the use of trigonometric functions to represent the out-of-plane displacement W and involve Galerkin type approximations to reduce the governing partial differential equations to simple algebraic expressions, the Level-2 solutions employ a truncated Fourier expansion in the circumferential direction and solve the resulting ordinary differential equations forming a nonlinear eigenvalue problem by a Stodola like technique. This involves the numerical integration of the set of ordinary differential equations, which in turn makes it possible to satisfy the specified boundary conditions exactly. It is shown that in many applications the Level-1 predictions concerning the critical buckling load of the perfect structure and its imperfection sensitivity are quite accurate and reliable. For those cases where this is not true, it is shown how the pre and postprocessing facilities of DISDECO can be used efficiently to arrive at a reliable prediction. This involves always the use of the Level-2 computational module ANILISA with its rigorous prebuckling solution and its feature which makes it possible to satisfy the specified boundary conditions exactly.

Arbocz, J.↗

NASA NTRS1993

Inviscid instability of streamwise corner flow

Linear stability of the incompressible flow along a stream wise corner is studied by solving the two-dimensional eigenvalue problem governed by partial differential equations. It is found that this fully three-dimensional flow is subject to inviscid instability due to the inflectional nature of the stream wise velocity profile. The higher growth rates for the inviscid instability mode, which is symmetric about the corner bisector, as compared to the viscous Tollmien-Schlichting instability operative away from the corner is consistent with the experimental findings that the corner flow transitions to turbulence earlier than the two-dimensional Blasius flow away from the corner.

Balachandar, S.↗

NASA NTRS1993

Linear and nonlinear PSE for compressible boundary layers

Compressible stability of growing boundary layers is studied by numerically solving the partial differential equations under a parabolizing approximation. The resulting parabolized stability equations (PSE) account for nonparallel as well as nonlinear effects. Evolution of disturbances in compressible flat-plate boundary layers are studied for freestream Mach numbers ranging from 0 to 4.5. Results indicate that the effect of boundary-layer growth is important for linear disturbances. Nonlinear calculations are performed for various Mach numbers. Two-dimensional nonlinear results using the PSE approach agree well with those from direct numerical simulations using the full Navier-Stokes equations while the required computational time is less by an order of magnitude. Spatial simulation using PSE were carried out for both the fundamental and subharmonic type breakdown for a Mach 1.6 boundary layer. The promising results obtained show that the PSE method is a powerful tool for studying boundary-layer instabilities and for predicting transition over a wide range of Mach numbers.

Chang, Chau-Lyan↗

NASA NTRS1993

A multi-level solution algorithm for steady-state Markov chains

A new iterative algorithm, the multi-level algorithm, for the numerical solution of steady state Markov chains is presented. The method utilizes a set of recursively coarsened representations of the original system to achieve accelerated convergence. It is motivated by multigrid methods, which are widely used for fast solution of partial differential equations. Initial results of numerical experiments are reported, showing significant reductions in computation time, often an order of magnitude or more, relative to the Gauss-Seidel and optimal SOR algorithms for a variety of test problems. The multi-level method is compared and contrasted with the iterative aggregation-disaggregation algorithm of Takahashi.

Horton, Graham↗

NASA NTRS1993

MGGHAT: Elliptic PDE software with adaptive refinement, multigrid and high order finite elements

MGGHAT (MultiGrid Galerkin Hierarchical Adaptive Triangles) is a program for the solution of linear second order elliptic partial differential equations in two dimensional polygonal domains. This program is now available for public use. It is a finite element method with linear, quadratic or cubic elements over triangles. The adaptive refinement via newest vertex bisection and the multigrid iteration are both based on a hierarchical basis formulation. Visualization is available at run time through an X Window display, and a posteriori through output files that can be used as GNUPLOT input. In this paper, we describe the methods used by MGGHAT, define the problem domain for which it is appropriate, illustrate use of the program, show numerical and graphical examples, and explain how to obtain the software.

Mitchell, William F.↗

NASA NTRS1993

An analytic study of a two-phase laminar airfoil in simulated heavy rain

A mathematical model for a two-phase flow laminar airfoil in simulated heavy rain has been established. The set of non-linear partial differential equations has been converted into a set of finite difference equations; appropriate initial and boundary conditions are provided. The numerical results are compared with the experimental measurements. They show good agreement in quality.

Hsu, Yu-Kao↗

NASA NTRS1993

A dual reciprocal boundary element formulation for viscous flows

The advantages inherent in the boundary element method (BEM) for potential flows are exploited to solve viscous flow problems. The trick is the introduction of a so-called dual reciprocal technique in which the convective terms are represented by a global function whose unknown coefficients are determined by collocation. The approach, which is necessarily iterative, converts the governing partial differential equations into integral equations via the distribution of fictitious sources or dipoles of unknown strength on the boundary. These integral equations consist of two parts. The first is a boundary integral term, whose kernel is the unknown strength of the fictitious sources and the fundamental solution of a convection-free flow problem. The second part is a domain integral term whose kernel is the convective portion of the governing PDEs. The domain integration can be transformed to the boundary by using the dual reciprocal (DR) concept. The resulting formulation is a pure boundary integral computational process.

Lafe, Olu↗

NASA NTRS1993

A theorem regarding roots of the zero-order Bessel function of the first kind

This paper investigates a problem on the steady-state, conduction-convection heat transfer process in cylindrical porous heat exchangers. The governing partial differential equations for the system are obtained using the energy conservation law. Solution of these equations and the concept of enthalpy lead to a new approach to prove a theorem that the sum of inverse squares of all the positive roots of the zero order Bessel function of the first kind equals to one-forth. As a corollary, it is shown that the sum of one over pth power (p greater than or equal to 2) of the roots converges to some constant.

Lin, X.-A.↗

NASA NTRS1993

An object-oriented approach for parallel self adaptive mesh refinement on block structured grids

Self-adaptive mesh refinement dynamically matches the computational demands of a solver for partial differential equations to the activity in the application's domain. In this paper we present two C++ class libraries, P++ and AMR++, which significantly simplify the development of sophisticated adaptive mesh refinement codes on (massively) parallel distributed memory architectures. The development is based on our previous research in this area. The C++ class libraries provide abstractions to separate the issues of developing parallel adaptive mesh refinement applications into those of parallelism, abstracted by P++, and adaptive mesh refinement, abstracted by AMR++. P++ is a parallel array class library to permit efficient development of architecture independent codes for structured grid applications, and AMR++ provides support for self-adaptive mesh refinement on block-structured grids of rectangular non-overlapping blocks. Using these libraries, the application programmers' work is greatly simplified to primarily specifying the serial single grid application and obtaining the parallel and self-adaptive mesh refinement code with minimal effort. Initial results for simple singular perturbation problems solved by self-adaptive multilevel techniques (FAC, AFAC), being implemented on the basis of prototypes of the P++/AMR++ environment, are presented. Singular perturbation problems frequently arise in large applications, e.g. in the area of computational fluid dynamics. They usually have solutions with layers which require adaptive mesh refinement and fast basic solvers in order to be resolved efficiently.

Lemke, Max↗

NASA NTRS1993

The combustion program at CTR

Understanding and modeling of turbulent combustion are key problems in the computation of numerous practical systems. Because of the lack of analytical theories in this field and of the difficulty of performing precise experiments, direct numerical simulation (DNS) appears to be one of the most attractive tools to use in addressing this problem. The general objective of DNS of reacting flows is to improve our knowledge of turbulent combustion but also to use this information for turbulent combustion models. For the foreseeable future, numerical simulation of the full three-dimensional governing partial differential equations with variable density and transport properties as well as complex chemistry will remain intractable; thus, various levels of simplification will remain necessary. On one hand, the requirement to simplify is not necessarily a handicap: numerical simulations allow the researcher a degree of control in isolating specific physical phenomena that is inaccessible in experiments. CTR has pursued an intensive research program in the field of DNS for turbulent reacting flows since 1987. DNS of reacting flows is quite different from DNS of non-reacting flows: without reaction, the equations to solve are clearly the five conservation equations of the Navier Stokes system for compressible situations (four for incompressible cases), and the limitation of the approach is the Reynolds number (or in other words the number of points in the computation). For reacting flows, the choice of the equations, the species (each species will require one additional conservation equation), the chemical scheme, and the configuration itself is more complex.

Poinsot, Thierry J.↗

NASA NTRS1993

Projection techniques for iterative solution of A(bar x) equals (bar b) with successive right-hand sides

Two projection techniques for computing approximate solutions to linear systems of the form A(bar x)(sup n) = (bar b)(sup n), for a sequence n = 1, 2, ..., e.g., such as arises from time discretization of a partial differential equation, are presented. The inexpensive approximate solutions can be used as initial guesses for iterative solution of the system, resulting in significantly reduced computational expense. Examples of two- and three-dimensional incompressible Navier-Stokes calculations are presented in which x represents the pressure, and A is a discrete Poisson operator. In flows containing significant dynamic activity, these projection techniques lead to as much as a two-fold reduction in solution time.

Fischer, Paul F.↗