Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “variable coefficients”

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 37 records · Page 2

Surface properties and effective resolution from ground-based observations of Mercury

A detailed study to evaluate ground-based photographs of Mercury has been carried out. Models of the surface scattering properties have been assumed and smeared with a Gaussian function for direct comparison with center-to-limb scans along Mercury's intensity equator. Data from a range of phase angles from 31 to 92 deg have been compared with smeared models assuming a Lambert surface, a surface which obeys the Lommel-Seeliger law and one which is Minnaertian, having a variable coefficient. Within the limits of the observations a lunar Minnaert surface yields the most consistent interpretation. An objective evaluation of the resolution of the photographs is obtained in terms of Gaussian half-widths.

Beebe, R. F.↗

Power series solutions of ordinary differential equations in MACSYMA

A program is described which extends the differential equation solving capability of MACSYMA to power series solutions and is available via the SHARE library. The program is directed toward those classes of equations with variable coefficients (in particular, those with singularities) and uses the method of Frobenius. Probably the most important distinction between this package and others currently available or being developed is that, wherever possible, this program will attempt to provide a complete solution to the equation rather than an approximation, i.e., a finite number of terms. This solution will take the form of a sum of infinite series.

Lafferty, E. L.↗

On improving the iterative convergence properties of an implicit approximate-factorization finite difference algorithm

The iterative convergence properties of an approximate-factorization implicit finite-difference algorithm are analyzed both theoretically and numerically. Modifications to the base algorithm were made to remove the inconsistency in the original implementation of artificial dissipation. In this way, the steady-state solution became independent of the time-step, and much larger time-steps can be used stably. To accelerate the iterative convergence, large time-steps and a cyclic sequence of time-steps were used. For a model transonic flow problem governed by the Euler equations, convergence was achieved with 10 times fewer time-steps using the modified differencing scheme. A particular form of instability due to variable coefficients is also analyzed.

Desideri, J. A.↗

General heating distributions on the tumbling infinite cylinder

The Shuttle program plans to recover all flight hardware except the External Tank (ET). The ET reentry task requires the quaranteed safe disposal in a remote section and specific footprint of the Indian Ocean or the Antarctic Pacific. The ET is separated from the Orbiter at near orbital speeds and must reenter tumbling, to neutralize lift in an essentially ballistic trajectory. The aerothermodynamic heating becomes an intriguing, complex, and sophisticated series of flow field, boundary layer, and gas dynamic problems that is an important new technology in the reentering disposal task. The singular aspect given in this paper is the general heating distributions for all tumbling attitudes. The objective is to maintain simplicity in form but with variable coefficients to be analytically scaled to account for free molecule, convective laminar and turbulent flow, real gas dissociation or any Reynolds number or Mach number effects.

Carroll, H. R.↗

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

A transformation of the boundary layer equations for free convection past a vertical flat plate with arbitrary blowing and wall temperature variations

A transformation of the laminar boundary layer equations, similar to the transformation of Kao et al. (1977), is presented which allows arbitrary distributions of both wall temperature and blowing. The procedure yields constant boundary conditions, but variable coefficients appear in the differential equations. The proposed transformation is applicable to the problem of downward burning of vertical pieces of condensed-phase combustibles.

Vedhanayagam, M.↗

Development of an unsteady aerodynamic analysis for finite-deflection subsonic cascades

An unsteady potential flow analysis, which accounts for the effects of blade geometry and steady turning, was developed to predict aerodynamic forces and moments associated with free vibration or flutter phenomena in the fan, compressor, or turbine stages of modern jet engines. Based on the assumption of small amplitude blade motions, the unsteady flow is governed by linear equations with variable coefficients which depend on the underlying steady low. These equations were approximated using difference expressions determined from an implicit least squares development and applicable on arbitrary grids. The resulting linear system of algebraic equations is block tridiagonal, which permits an efficient, direct (i.e., noniterative) solution. The solution procedure was extended to treat blades with rounded or blunt edges at incidence relative to the inlet flow.

Verdon, J. M.↗

The stability of pseudospectral-Chebyshev methods

The pseudospectral-Chebyshev methods are shown to be convergent in variable coefficient problems and, in some cases, hyperbolic problems. The analysis demonstrates that the rate of convergence is greater for finite difference methods or the finite element method. For a single first-order hyperbolic equation, the method is seen as remaining stable even when the coefficient changes sign, although in this case it is specified that care must be taken to have adequate spatial resolution. It is noted that this fact, combined with the fact that collocation methods are easy to apply in the nonlinear case, shows that the pseudospectral method is in general preferable to the Galerkin or Tau methods.

Gottlieb, D.↗

Stability theorems for multidimensional linear systems with variable parameters

A Liapunov-type approach is used to derive two equivalent theorems which govern the stability of coupled linear systems with varying multiple parameters. The theorems generalize some of the existing theorems applicable to systems with constant parameters and the Sonin-Polya theorem applicable to a single-degree-of-freedom system with variable coefficients. As an illustration, the proposed theorems are applied to mechanical systems with varying inertia, stiffness, gyroscopic, and damping terms, and velocity and position-dependent forces.

Shrivastava, S. K.↗

Iterative spectral methods and spectral solutions to compressible flows

A spectral multigrid scheme is described which can solve pseudospectral discretizations of self-adjoint elliptic problems in O(N log N) operations. An iterative technique for efficiently implementing semi-implicit time-stepping for pseudospectral discretizations of Navier-Stokes equations is discussed. This approach can handle variable coefficient terms in an effective manner. Pseudospectral solutions of compressible flow problems are presented. These include one dimensional problems and two dimensional Euler solutions. Results are given both for shock-capturing approaches and for shock-fitting ones.

Hussaini, M. Y.↗

A solution procedure for behavior of thick plates on a nonlinear foundation and postbuckling behavior of long plates

Approximate solutions for three nonlinear orthotropic plate problems are presented: (1) a thick plate attached to a pad having nonlinear material properties which, in turn, is attached to a substructure which is then deformed; (2) a long plate loaded in inplane longitudinal compression beyond its buckling load; and (3) a long plate loaded in inplane shear beyond its buckling load. For all three problems, the two dimensional plate equations are reduced to one dimensional equations in the y-direction by using a one dimensional trigonometric approximation in the x-direction. Each problem uses different trigonometric terms. Solutions are obtained using an existing algorithm for simultaneous, first order, nonlinear, ordinary differential equations subject to two point boundary conditions. Ordinary differential equations are derived to determine the variable coefficients of the trigonometric terms.

Stein, M.↗

An approximate buckling analysis for rectangular orthotropic plates with centrally located cutouts

An approximate analysis for predicting buckling of rectangular orthotropic composite plates with centrally located cutouts is presented. In this analysis, prebuckling and buckling problems are converted from a two-dimensional to a one-dimensional system of linear differential equations with variable coefficients. The conversion is accomplished by expressing the displacements as series with each element containing a trigonometric function of one coordinate and a coefficient that is an arbitrary function of the other coordinate. Ordinary differential equations are then obtained from a variational principle. Analytical results obtained from the approximate analysis are compared with finite element analyses for isotropic plates and for specially orthotropic plates with central circular cutouts of various sizes. Experimental results for the specially orthotropic plates are also presented. In nearly all cases, the approximate analysis predicts the buckling mode shapes correctly and predicts the buckling loads to within a few percent of the finite element and experimental results.

Nemeth, M. P.↗

Regularity estimates up to the boundary for elliptic systems of difference equations

Regularity estimates up to the boundary for solutions of elliptic systems of finite difference equations were proved. The regularity estimates, obtained for boundary fitted coordinate systems on domains with smooth boundary, involve discrete Sobolev norms and are proved using pseudo-difference operators to treat systems with variable coefficients. The elliptic systems of difference equations and the boundary conditions which are considered are very general in form. The regularity of a regular elliptic system of difference equations was proved equivalent to the nonexistence of eigensolutions. The regularity estimates obtained are analogous to those in the theory of elliptic systems of partial differential equations, and to the results of Gustafsson, Kreiss, and Sundstrom (1972) and others for hyperbolic difference equations.

Strikwerda, J. C.↗

A linearized Euler analysis of unsteady flows in turbomachinery

A method for calculating unsteady flows in cascades is presented. The model, which is based on the linearized unsteady Euler equations, accounts for blade loading shock motion, wake motion, and blade geometry. The mean flow through the cascade is determined by solving the full nonlinear Euler equations. Assuming the unsteadiness in the flow is small, then the Euler equations are linearized about the mean flow to obtain a set of linear variable coefficient equations which describe the small amplitude, harmonic motion of the flow. These equations are discretized on a computational grid via a finite volume operator and solved directly subject to an appropriate set of linearized boundary conditions. The steady flow, which is calculated prior to the unsteady flow, is found via a Newton iteration procedure. An important feature of the analysis is the use of shock fitting to model steady and unsteady shocks. Use of the Euler equations with the unsteady Rankine-Hugoniot shock jump conditions correctly models the generation of steady and unsteady entropy and vorticity at shocks. In particular, the low frequency shock displacement is correctly predicted. Results of this method are presented for a variety of test cases. Predicted unsteady transonic flows in channels are compared to full nonlinear Euler solutions obtained using time-accurate, time-marching methods. The agreement between the two methods is excellent for small to moderate levels of flow unsteadiness. The method is also used to predict unsteady flows in cascades due to blade motion (flutter problem) and incoming disturbances (gust response problem).

Hall, Kenneth C.↗

Improving the convergence rate to steady state of parabolic ADI methods

The present, residuals' L(2)-norms analysis of the rate of convergence to steady state for parabolic ADI solvers allows the prediction of the number of iterations required for convergence, as a function of the Courant number alpha. A modification of current ADI codes is presented which significantly improves the convergence rate and is insensitive to the Courant number over a large range of alpha. This corrected algorithm is tested for the cases of Dirichlet problems for uniform grids of many mesh sizes, mixed Dirichlet-Neumann problems, and problems defined on stretched grids and/or problems with variable coefficients.

Abarbanel, Saul S.↗

Stability with large step sizes for multistep discretizations of stiff ordinary differential equations

One-leg and multistep discretizations of variable-coefficient linear systems of ODEs having both slow and fast time scales are investigated analytically. The stability properties of these discretizations are obtained independent of ODE stiffness and compared. The results of numerical computations are presented in tables, and it is shown that for large step sizes the stability of one-leg methods is better than that of the corresponding linear multistep methods.

Majda, George↗