SEARCH · Engineering Papers

Results for “Numerical partial differential equations”

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 127 records · Page 7

All sources NASA NTRS DOE OSTI

NASA NTRS1987

Numerical observations of dynamic behavior in two-dimensional compressible convection

Using numerical simulations of the anelastic partial differential equations, convective flow in a two-dimensional polytropic atmosphere is investigated. Fixing the aspect ratio at unity and choosing the viscosity as the order parameter, the numerical solutions have indicated a variety of dynamic behavior including (1) oscillatory motion with a period of order the eddy turnover time, (2) quasi-periodic motion with a long period much longer than an eddy turnover time, and (3) at least two branches of solutions in certain regimes.

Ginet, G. P.↗

NASA NTRS1980

An alternative approach to the numerical simulation of steady inviscid flow

A numerical procedure for the efficient simulation of steady inviscid flow is described and its utility demonstrated. Application of the surrogate equation technique allows the formulation of stable, fully conservative, type dependent finite difference equations for use in obtaining numerical solutions to systems of first order partial differential equations, such as the steady state Euler equations or their various approximations. Computational results are presented for the full Euler equations and for the transonic disturbance equations. For the latter case, a computational efficiency greater than that obtained by means of the standard perturbation potential approach is indicated.

Johnson, G. M.↗

NASA NTRS1981

Surrogate-equation technique for simulation of steady inviscid flow

A numerical procedure for the iterative solution of inviscid flow problems is described, and its utility for the calculation of steady subsonic and transonic flow fields is demonstrated. Application of the surrogate equation technique defined herein allows the formulation of stable, fully conservative, type dependent finite difference equations for use in obtaining numerical solutions to systems of first order partial differential equations, such as the steady state Euler equations. Steady, two dimensional solutions to the Euler equations for both subsonic, rotational flow and supersonic flow and to the small disturbance equations for transonic flow are presented.

Johnson, G. M.↗

NASA NTRS1996

Multigrid Methods for Fully Implicit Oil Reservoir Simulation

In this paper we consider the simultaneous flow of oil and water in reservoir rock. This displacement process is modeled by two basic equations: the material balance or continuity equations and the equation of motion (Darcy's law). For the numerical solution of this system of nonlinear partial differential equations there are two approaches: the fully implicit or simultaneous solution method and the sequential solution method. In the sequential solution method the system of partial differential equations is manipulated to give an elliptic pressure equation and a hyperbolic (or parabolic) saturation equation. In the IMPES approach the pressure equation is first solved, using values for the saturation from the previous time level. Next the saturations are updated by some explicit time stepping method; this implies that the method is only conditionally stable. For the numerical solution of the linear, elliptic pressure equation multigrid methods have become an accepted technique. On the other hand, the fully implicit method is unconditionally stable, but it has the disadvantage that in every time step a large system of nonlinear algebraic equations has to be solved. The most time-consuming part of any fully implicit reservoir simulator is the solution of this large system of equations. Usually this is done by Newton's method. The resulting systems of linear equations are then either solved by a direct method or by some conjugate gradient type method. In this paper we consider the possibility of applying multigrid methods for the iterative solution of the systems of nonlinear equations. There are two ways of using multigrid for this job: either we use a nonlinear multigrid method or we use a linear multigrid method to deal with the linear systems that arise in Newton's method. So far only a few authors have reported on the use of multigrid methods for fully implicit simulations. Two-level FAS algorithm is presented for the black-oil equations, and linear multigrid for two-phase flow problems with strong heterogeneities and anisotropies is studied. Here we consider both possibilities. Moreover we present a novel way for constructing the coarse grid correction operator in linear multigrid algorithms. This approach has the advantage in that it preserves the sparsity pattern of the fine grid matrix and it can be extended to systems of equations in a straightforward manner. We compare the linear and nonlinear multigrid algorithms by means of a numerical experiment.

Molenaar, J.↗

NASA NTRS1980

An interaction algorithm for prediction of mean and fluctuating velocities in two-dimensional aerodynamic wake flows

A theoretical analysis is presented yielding sets of partial differential equations for determination of turbulent aerodynamic flowfields in the vicinity of an airfoil trailing edge. A four phase interaction algorithm is derived to complete the analysis. Following input, the first computational phase is an elementary viscous corrected two dimensional potential flow solution yielding an estimate of the inviscid-flow induced pressure distribution. Phase C involves solution of the turbulent two dimensional boundary layer equations over the trailing edge, with transition to a two dimensional parabolic Navier-Stokes equation system describing the near-wake merging of the upper and lower surface boundary layers. An iteration provides refinement of the potential flow induced pressure coupling to the viscous flow solutions. The final phase is a complete two dimensional Navier-Stokes analysis of the wake flow in the vicinity of a blunt-bases airfoil. A finite element numerical algorithm is presented which is applicable to solution of all partial differential equation sets of inviscid-viscous aerodynamic interaction algorithm. Numerical results are discussed.

Baker, A. J.↗

NASA NTRS1985

Generation of Surface Coordinates by Elliptic Partial Differential Equations

The problem of generating spatial coordinates by numerical methods through carefully selected mathematical models is of current interest both in mechanics and physics. The problem of generation of a desired system of coordinates in a given surface was considered, which essentially is an effort directed to the problem of grid generation in a two-dimensional non-Euclidean space. The mathematical model selected for this purpose is based on the formulae of Gauss for a surface. The proposed equations can be used to generate a new coordinate system from the data of an already given coordinate system in a surface. If the coefficients of the first and second fundamental forms have been given, then the proposed equations can be used to generate a surface satisfying the given data (surface fitting). The proposed equations can also be used to generate surfaces in the space between two arbitrary given surfaces, thus providing 3D grids in an Euclidean space.

Warsi, Z. U. A.↗

NASA NTRS1989

Nonlinear interactions in mixing layers and compressible heated round jets

The nonlinear interactions between a fundamental instability mode and both its harmonics and the changing mean flow are studied using the weakly nonlinear stability theory of Stuart and Watson, and numerical solutions of coupled nonlinear partial differential equations. The first part focuses on incompressible cold (or isothermal; constant temperature throughout) mixing layers, and for these, the first and second Landau constants are calculated as functions of wavenumber and Reynolds number. It is found that the dominant contribution to the Landau constants arises from the mean flow changes and not from the higher harmonics. In order to establish the range of validity of the weakly nonlinear theory, the weakly nonlinear and numerical solutions are compared and the limitation of each is discussed. At small amplitudes and at low-to-moderate Reynolds numbers, the two results compare well in describing the saturation of the fundamental, the distortion of the mean flow, and the initial stages of vorticity roll-up. At larger amplitudes, the interaction between the fundamental, second harmonic, and the mean flow is strongly nonlinear and the numerical solution predicts flow oscillations, whereas the weakly nonlinear theory yields saturation. In the second part, the weakly nonlinear theory is extended to heated (or nonisothermal; mean temperature distribution) subsonic round jets where quadratic and cubic nonlinear interactions are present, and the Landau constants also depend on jet temperature ratio, Mach number and azimuthal mode number. Under exponential growth and nonlinear saturation, it is found that heating and compressibility suppress the growth of instability waves, that the first azimuthal mode is the dominant instability mode, and that the weakly nonlinear solution describes the early stages of the roll-up of an axisymmetric shear layer. The receptivity of a typical jet flow to pulse type input disturbance is also studied by solving the initial value problem and then examining the behavior of the long-time solution.

Jarrah, Yousef Mohd↗

NASA NTRS1999

A Parallel Compact Multi-Dimensional Numerical Algorithm with Aeroacoustics Applications

In this study we propose a novel method to parallelize high-order compact numerical algorithms for the solution of three-dimensional PDEs (Partial Differential Equations) in a space-time domain. For this numerical integration most of the computer time is spent in computation of spatial derivatives at each stage of the Runge-Kutta temporal update. The most efficient direct method to compute spatial derivatives on a serial computer is a version of Gaussian elimination for narrow linear banded systems known as the Thomas algorithm. In a straightforward pipelined implementation of the Thomas algorithm processors are idle due to the forward and backward recurrences of the Thomas algorithm. To utilize processors during this time, we propose to use them for either non-local data independent computations, solving lines in the next spatial direction, or local data-dependent computations by the Runge-Kutta method. To achieve this goal, control of processor communication and computations by a static schedule is adopted. Thus, our parallel code is driven by a communication and computation schedule instead of the usual "creative, programming" approach. The obtained parallelization speed-up of the novel algorithm is about twice as much as that for the standard pipelined algorithm and close to that for the explicit DRP algorithm.

Povitsky, Alex↗

NASA NTRS1987

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 NTRS1988

Problem size, parallel architecture, and optimal speedup

Nicol, David M.↗

NASA NTRS1980

An alternative approach to the numerical simulation of steady inviscid flow

A numerical procedure for the efficient simulation of steady inviscid flow is described and its utility is demonstrated. The method is uniformly valid for application in the subsonic, transonic and supersonic flow regimes. It does not rely on the introduction of additional assumptions beyond those necessary to obtain the Euler equations from the Navier-Stokes equations, nor does it make use of a time-asymptotic solution of the unsteady equations of motion. Application of the herein-defined surrogate equation technique allows the formulation of stable, fully-conservative, type-dependent finite difference equations for use in obtaining numerical solutions to systems of first-order partial differential equations, such as the steady-state Euler equations or their various approximations. Computational results are presented for the full Euler equations used to simulate rotational subsonic flow and for the transonic small disturbance equations. For the latter case, a computational efficiency greater than that obtained by means of the standard perturbation potential approach is indicated.

Johnson, G. M.↗

NASA NTRS1986

Adaptive computational methods for SSME internal flow analysis

Adaptive finite element methods for the analysis of classes of problems in compressible and incompressible flow of interest in SSME (space shuttle main engine) analysis and design are described. The general objective of the adaptive methods is to improve and to quantify the quality of numerical solutions to the governing partial differential equations of fluid dynamics in two-dimensional cases. There are several different families of adaptive schemes that can be used to improve the quality of solutions in complex flow simulations. Among these are: (1) r-methods (node-redistribution or moving mesh methods) in which a fixed number of nodal points is allowed to migrate to points in the mesh where high error is detected; (2) h-methods, in which the mesh size h is automatically refined to reduce local error; and (3) p-methods, in which the local degree p of the finite element approximation is increased to reduce local error. Two of the three basic techniques have been studied in this project: an r-method for steady Euler equations in two dimensions and a p-method for transient, laminar, viscous incompressible flow. Numerical results are presented. A brief introduction to residual methods of a-posterior error estimation is also given and some pertinent conclusions of the study are listed.

Oden, J. T.↗

NASA NTRS1983

Finite element or Galerkin type semidiscrete schemes

A finite element of Galerkin type semidiscrete method is proposed for numerical solution of a linear hyperbolic partial differential equation. The question of stability is reduced to the stability of a system of ordinary differential equations for which Dahlquist theory applied. Results of separating the part of numerical solution which causes the spurious oscillation near shock-like response of semidiscrete scheme to a step function initial condition are presented. In general all methods produce such oscillatory overshoots on either side of shocks. This overshoot pathology, which displays a behavior similar to Gibb's phenomena of Fourier series, is explained on the basis of dispersion of separated Fourier components which relies on linearized theory to be satisfactory. Expository results represented.

Durgun, K.↗

NASA NTRS1979

Aperture amplitude and phase control of offset dual reflectors

An approximate synthesis method is presented for solving the problem of offset dual-shaped reflectors. It is shown that the equations obtained by this method are partial differential equations that can be integrated numerically, except that these equations are generally not total and therefore do not have a 'smooth' solution. However, it is demonstrated that these equations form a nearly total differential that enables one to generate a smooth subreflector by integrating the differential equations and then synthesizing a main reflector which gives an exact solution for the specified aperture phase distribution. Several important illustrative examples are treated, which validate the adequacy of the proposed method.

Galindo-Israel, V.↗

NASA NTRS1991

Modeling flow at the nozzle of a solid rocket motor

The mechanical behavior of a rocket motor internal flow field results in a system of nonlinear partial differential equations which can be solved numerically. The accuracy and the convergence of the solution of the system of equations depends largely on how precisely the sharp gradients can be resolved. An adaptive grid generation scheme is incorporated into the computer algorithm to enhance the capability of numerical modeling. With this scheme, the grid is refined as the solution evolves. This scheme significantly improves the methodology of solving flow problems in rocket nozzle by putting the refinement part of grid generation into the computer algorithm.

Chow, Alan S.↗

NASA NTRS1991

Stability of the flow around a cylinder: The spin-up problem

A concern is the flow around an infinite cylinder, which at a certain instant impulsively starts to spin. The growth of vortices in the resulting boundary layer occurring outside the cylinder is investigated. This layer is essentially a Rayleigh layer which grows with time, so the mechanism involved is similar to that studied in Hall (1983). Vortices with wavenumber comparable to the layer thickness are shown to be described by partial differential equations that govern the system numerically. It is assumed that the Rayleigh layer is thin, so particles are confined to move in a path with radius of curvature the same as the cylinder. The Goertler number is a function of time, so the time scale which produces an order, is considered one Goertler number. The right hand branch calculation is considered by letting the time tend to infinity, also inviscid Goertler modes are considered.

Otto, S. R.↗

NASA NTRS1990

On the nature of the interaction between structures and actuators in vibration suppression

This lecture presents a summary of work on the analysis of the interaction between a structure, an actuator used to suppress the vibration of the structure, and the control law implemented by the actuator. Two control applications are considered. First, a proof-mass actuator with experimentally verified dynamics, capable of being used in a space structure configuration is examined. This is connected to a cantilevered beam modeled as a single degree of freedom system. Secondly, an electric motor, again with experimentally verified dynamics, is used to slew a beam modeled by a partial differential equation. Both the experimental and numerical configuration predict the presence of potential instabilities in system performance if proper consideration is not given to interactions between the control law, the structure and the actuator. In addition, an understanding of the interaction is shown to greatly effect performance. In particular it is shown that judicious choice of actuator parameters greatly improves closed loop performance.

Inman, Daniel J.↗

NASA NTRS1995

Wavenumber selection for small-wavelength Goertler vortices in curved channel flows

The problem of wavenumber selection for fully nonlinear, small-wavelength Goertler vortices in a curved channel flow is considered. These types of Goertler vortices were first considered by Hall & Lakin (1988) for an external boundary layer flow. They proved particularly amenable to asymptotic description, it was possible to consider vortices large enough so that the mean flow correction driven by them is as large as the basic state, and this prompted the authors to consider them in a curved channel flow as an initial application of the phase-equation approach to Goertler vortices. This involves the assumption that the phase variable of these Goertler vortices varies on slow spanwise and time scales, then an analysis of both inside and outside the core region, to which vortex activity is restricted, leads to a system of partial differential equations which can be solved numerically for the wavenumber. The authors consider in particular the effect on the wavenumber of the outer channel wall varying on the same slow spanwise scale as the phase variable.

Dando, Andrew↗