Search NASA⌕ Search

SEARCH · Search NASA

Results for “PARABOLIC 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

Towards developing robust algorithms for solving partial differential equations on MIMD machines

Methods for efficient computation of numerical algorithms on a wide variety of MIMD machines are proposed. These techniques reorganize the data dependency patterns to improve the processor utilization. The model problem finds the time-accurate solution to a parabolic partial differential equation discretized in space and implicitly marched forward in time. The algorithms are extensions of Jacobi and SOR. The extensions consist of iterating over a window of several timesteps, allowing efficient overlap of computation with communication. The methods increase the degree to which work can be performed while data are communicated between processors. The effect of the window size and of domain partitioning on the system performance is examined both by implementing the algorithm on a simulated multiprocessor system.

Saltz, J. H.↗

Towards developing robust algorithms for solving partial differential equations on MIMD machines

Methods for efficient computation of numerical algorithms on a wide variety of MIMD machines are proposed. These techniques reorganize the data dependency patterns to improve the processor utilization. The model problem finds the time-accurate solution to a parabolic partial differential equation discretized in space and implicitly marched forward in time. The algorithms are extensions of Jacobi and SOR. The extensions consist of iterating over a window of several timesteps, allowing efficient overlap of computation with communication. The methods increase the degree to which work can be performed while data are communicated between processors. The effect of the window size and of domain partitioning on the system performance is examined both by implementing the algorithm on a simulated multiprocessor system.

Saltz, Joel H.↗

Parameter estimation problems for distributed systems using a multigrid method

The problem of estimating spatially varying coefficients of partial differential equations is considered from observation of the solution and of the right hand side of the equation. It is assumed that the observations are distributed in the domain and that enough observations are given. A method of discretization and an efficient multigrid method for solving the resulting discrete systems are described. Numerical results are presented for estimation of coefficients in an elliptic and a parabolic partial differential equation.

Taasan, S.↗

Parameter estimation problems for distributed systems using a multigrid method

The problem of estimating spatially varying coefficients of partial differential equations is considered from observation of the solution and of the right hand side of the equation. It is assumed that the observations are distributed in the domain and that enough observations are given. A method of discretization and an efficient multigrid method for solving the resulting discrete systems are described. Numerical results are presented for estimation of coefficients in an elliptic and a parabolic partial differential equation.

Ta'asan, Shlomo↗

Interface procedures for overlapping grids

The solution of partial differential equations by means of either implicit or conservative explicit finite-difference methods, on the one hand, or multicomponent overlapping grid schemes on the other, is presently subjected to interpolation at grid boundaries in an attempt to eliminate the time lag. The accuracy of the transient solution of a hyperbolic or parabolic partial differential equation is dependent on the procedures used to transfer information between grids in a composite grid system. The present techniques reduce numerical solution errors.

Mastin, C. Wayne↗

CAN-DO, CFD-based Aerodynamic Nozzle Design and Optimization program for supersonic/hypersonic wind tunnels

A design program is developed which incorporates a modern approach to the design of supersonic/hypersonic wind-tunnel nozzles. The approach is obtained by the coupling of computational fluid dynamics (CFD) with design optimization. The program can be used to design a 2D or axisymmetric, supersonic or hypersonic, wind-tunnel nozzles that can be modeled with a calorically perfect gas. The nozzle design is obtained by solving a nonlinear least-squares optimization problem (LSOP). The LSOP is solved using an iterative procedure which requires intermediate flowfield solutions. The nozzle flowfield is simulated by solving the Navier-Stokes equations for the subsonic and transonic flow regions and the parabolized Navier-Stokes equations for the supersonic flow regions. The advantages of this method are that the design is based on the solution of the viscous equations eliminating the need to make separate corrections to a design contour, and the flexibility of applying the procedure to different types of nozzle design problems.

Korte, John J.↗

The development of a three-dimensional partially elliptic flow computer program for combustor research

A three dimensional, partially elliptic, computer program was developed. Without requiring three dimensional computer storage locations for all flow variables, the partially elliptic program is capable of predicting three dimensional combustor flow fields with large downstream effects. The program requires only slight increase of computer storage over the parabolic flow program from which it was developed. A finite difference formulation for a three dimensional, fully elliptic, turbulent, reacting, flow field was derived. Because of the negligible diffusion effects in the main flow direction in a supersonic combustor, the set of finite-difference equations can be reduced to a partially elliptic form. Only the pressure field was governed by an elliptic equation and requires three dimensional storage; all other dependent variables are governed by parabolic equations. A numerical procedure which combines a marching integration scheme with an iterative scheme for solving the elliptic pressure was adopted.

Pan, Y. S.↗

Explicit-implicit staggered procedure for multibody dynamics analysis

An explicit-implicit staggered time-integration procedure is presented for the solution of multibody dynamical equations involving large rotations and constraints. The algorithms adopts a two-stage modification of the central difference algorithm for integrating the translational coordinates and the angular velocity vector, and the midpoint implicit algorithm to solve the kinematical relation in terms of the Euler parameters for updating the angular orientations. The Lagrange multipliers to enforce the system constraints are obtained by implicitly integrating a parabolically regularized differential equation for the multipliers. The performance of the present procedure has been evaluated to applying the procedure to solve several sample problems. The results indicate that the procedure is robust in dealing with a variety of constraints and spatial kinematic motions, hence it is recommended for applications to general multibody dynamics analyses.

Park, K. C.↗

The parabolic approximation for sound propagation in a stratified moving medium

Propagation of sound in a stratified moving medium is discussed through an extension of the parabolic approximation to the acoustic equations of motion for short wavelengths. The parabolic approximation is related to the theory of geometric acoustics, and it is shown that it yields an improvement in accuracy over geometric theory. Also, the approximation corrects cumulative failures of geometric theory which occur when sound propagates many wavelengths from its source. The theory is illustrated by application to simple examples of quasi-plane wave propagation.

Myers, M. K.↗

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.↗

Two examples in the time optimal control theory of distributed parameter systems.

The behavior of hyperbolic and parabolic partial differential equations is contrasted by studying the point-to-point time-optimal control problem for the equation of heat conduction and the equation of motion of a vibrating string. A maximal principle is obtained for the time-optimal control of the one-dimensional heat equation, and it is proven that time optimal controls are weakly bang-bang. The bang-bang principle is proven to be invalid for hyperbolic equations because of the finite speed of wave propagation. In the case of boundary value control of the vibrating spring, the latter is demonstrated by deriving an explicit formula for the time optimal control.

Quinn, J. P.↗

Propagation of pulse trains through a random medium

Applying the parabolic approximation, the equations for two-frequency symmetric and antisymmetric mutual intensity functions for waves propagating through a random medium are derived, including the multiple scattering effects. These functions are applied to derive the general formulas for the covariance functions of narrow-band pulses. They are used to compute the signal intensities for pulse trains passing through an ionospheric irregularity slab.

Liu, C. H.↗

Steady-state solutions of a diffusive energy-balance climate model and their stability

A diffusive energy-balance climate model, governed by a nonlinear parabolic partial differential equation, was studied. Three positive steady-state solutions of this equation are found; they correspond to three possible climates of our planet: an interglacial (nearly identical to the present climate), a glacial, and a completely ice-covered earth. Models similar to the main one are considered, and the number of their steady states was determined. All the models have albedo continuously varying with latitude and temperature, and entirely diffusive horizontal heat transfer. The stability under small perturbations of the main model's climates was investigated. A stability criterion is derived, and its application shows that the present climate and the deep freeze are stable, whereas the model's glacial is unstable. The dependence was examined of the number of steady states and of their stability on the average solar radiation.

Ghil, M.↗

A formulation consideration for orthogonal collocation procedures

Numerical instabilities often arise in the use of high-ordered collocation approximations for numerically solving parabolic partial differential equations. These problems may be reduced by formulations involving evaluation of collocation polynomials rather than combination of the polynomials into a power series. As an illustration, two formulations using shifted Legendre polynomials of order 26 and less are compared.

Lashmet, P. K.↗

Climate stability for a Sellers-type model

We study a diffusive energy-balance climate model governed by a nonlinear parabolic partial differential equation. Three positive steady-state solutions of this equation are found; they correspond to three possible climates of our planet: an interglacial (nearly identical to the present climate), a glacial, and a completely ice-covered earth. We consider also models similar to the main one studied, and determine the number of their steady states. All the models have albedo continuously varying with latitude and temperature, and entirely diffusive horizontal heat transfer. The diffusion is taken to be nonlinear as well as linear. We investigate the stability under small perturbations of the main model's climates. A stability criterion is derived, and its application shows that the 'present climate' and the 'deep freeze' are stable, whereas the model's glacial is unstable. A variational principle is introduced to confirm the results of this stability analysis. For a sufficient decrease in solar radiation (about 2%) the glacial and interglacial solutions disappear, leaving the ice-covered earth as the only possible climate.

Ghil, M.↗

An extension of A-stability to alternating direction implicit methods

An alternating direction implicit (ADI) scheme was constructed by the method of approximate factorization. An A-stable linear multistep method (LMM) was used to integrate a model two-dimensional hyperbolic-parabolic partial differential equation. Sufficient conditions for the A-stability of the LMM were determined by applying the theory of positive real functions to reduce the stability analysis of the partial differential equations to a simple algebraic test. A linear test equation for partial differential equations is defined and then used to analyze the stability of approximate factorization schemes. An ADI method for the three-dimensional heat equation is also presented.

Warming, R. F.↗

Initial boundary value problems for the method of lines

This paper treats the stability of the initial boundary value problem for the method of lines applied to hyperbolic and parabolic partial differential equations in one space dimension. The theory treats the case of variable coefficients and allows for very general boundary conditions. Several examples are given which illustrate the theory. The theory is analogous to that developed by Gustafsson, Kreiss, and Sundstrom for finite-difference methods.

Strikwerda, J. C.↗