Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite-Difference Methods”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
NASA NTRS1976

A two-dimensional finite difference solution for the transient thermal behavior of tubular solar collector

A numerical procedure was established using the finite-difference technique in the determination of the time-varying temperature distribution of a tubular solar collector under changing solar radiancy and ambient temperature. Three types of spatial discretization processes were considered and compared for their accuracy of computations and for selection of the shortest computer time and cost. The stability criteria of this technique was analyzed in detail to give the critical time increment to ensure stable computations. The results of the numerical analysis were in good agreement with the analytical solution previously reported. The numerical method proved to be a powerful tool in the investigation of the collector sensitivity to two different flow patterns and several flow control mechanisms.

Lansing, F. L.↗
NASA NTRS1980

A two-dimensional finite-difference solution for the transient thermal behavior of a tubular solar collector

A numerical procedure was established using the finite-difference technique in the determination of the time-varying temperature distribution of a tubular solar collector under changing solar radiancy and ambient temperature. Three types of spatial discretization processes were considered and compared for their accuracy of computations and for selection of the shortest computer time and cost. The stability criteria of this technique were analyzed in detail to give the critical time increment to ensure stable computations. The results of the numerical analysis were in good agreement with the analytical solution previously reported. The numerical method proved to be a powerful tool in the investigation of the collector sensitivity to two different flow patterns and several flow control mechanisms.

Lansing, F. L.↗
NASA NTRS1981

An alternating-direction implicit algorithm for unsteady potential flow

An implicit finite-difference scheme is presented for the efficient computation of unsteady potential flow about airfoils. The formulation uses density and the velocity potential as dependent variables and is cast in conservation form to assure the theoretically correct determination of shockwave location and speed. To enable boundary conditions to be imposed directly on the airfoil surface, a time-varying sheared-rectilinear coordinate transformation is employed. Calculated time-history solutions on a pulsating airfoil are compared with the results of other unsteady transonic codes, including a previous method of the authors. The present method is demonstrated to be unconditionally stable and to give accurate solutions with sharply resolved shocks.

Chipman, R.↗
NASA NTRS1977

A perturbation solution for interlaminar stresses in bidirectional laminates

The paper considers a thin elastic bidirectional symmetric composite laminate subjected to a uniform axial strain. All interlaminar stresses are determined by a zeroth-order perturbation analysis of the governing equations. The method of solution considers the laminate to be composed of an interior region (removed from the free edge) and a boundary layer region. The interior region is analyzed by dropping terms related to powers of the thickness-to-width ratio in the dimensionless elasticity equations. The boundary layer region is studied by introducing a stretching transformation to the governing equations. A uniform solution is then formed by satisfying Prandtl's matching principle of perturbation theory. Results for a four-layer bidirectional graphite-epoxy laminate are compared with existing finite-difference results. It is shown that the perturbation analysis provides more insight into the stress distribution near the free edge.

Hsu, P. W.↗
NASA NTRS1982

A conservative finite difference algorithm for the unsteady transonic potential equation in generalized coordinates

An implicit, approximate-factorization, finite-difference algorithm has been developed for the computation of unsteady, inviscid transonic flows in two and three dimensions. The computer program solves the full-potential equation in generalized coordinates in conservation-law form in order to properly capture shock-wave position and speed. A body-fitted coordinate system is employed for the simple and accurate treatment of boundary conditions on the body surface. The time-accurate algorithm is modified to a conventional ADI relaxation scheme for steady-state computations. Results from two- and three-dimensional steady and two-dimensional unsteady calculations are compared with existing methods.

Bridgeman, J. O.↗
NASA NTRS1974

The numerical calculation of laminar boundary-layer separation

Iterative finite-difference techniques are developed for integrating the boundary-layer equations, without approximation, through a region of reversed flow. The numerical procedures are used to calculate incompressible laminar separated flows and to investigate the conditions for regular behavior at the point of separation. Regular flows are shown to be characterized by an integrable saddle-type singularity that makes it difficult to obtain numerical solutions which pass continuously into the separated region. The singularity is removed and continuous solutions ensured by specifying the wall shear distribution and computing the pressure gradient as part of the solution. Calculated results are presented for several separated flows and the accuracy of the method is verified. A computer program listing and complete solution case are included.

Klineberg, J. M.↗
NASA NTRS1974

On laminar boundary-layer separation

Iterative finite-difference techniques are developed for integrating the boundary-layer equations, without approximation, through a region of reversed flow. The numerical procedures are used to calculate incompressible laminar separated flows and to investigate the conditions for regular behavior at the point of separation. Regular flows are shown to be characterized by an integrable saddle-type singularity that makes it difficult to obtain numerical solutions which pass continuously into the separated region. The singularity is removed and continuous solutions ensured by specifying the wall shear distribution and computing the pressure gradient as part of the solution. Calculated results are presented for a number of separated flows and the accuracy of the method is verified.

Klineberg, J. M.↗
NASA NTRS1975

Time-dependent transonic flow solutions for axial turbomachinery

Three-dimensional unsteady transonic flow through an axial turbomachine stage is described in terms of a pair of two-dimensional formulations pertaining to orthogonal surfaces, namely, a blade-to-blade surface and a hub-to-casing surface. The resulting systems of nonlinear, inviscid, compressible equations of motion are solved by an explicit finite-difference technique. The blade-to-blade program includes the periodic interaction between rotor and stator blade rows. Treatment of the boundary conditions and of the blade slipstream motion by a characteristic type procedure is discussed in detail. Harmonic analysis of the acoustic far field produced by the blade row interaction, including an arbitrary initial transient, is outlined. Results from the blade-to-blade program are compared with experimental measurements of the rotating pressure field at the tip of a high-speed fan. The hub-to-casing program determines circumferentially averaged flow properties on a meridional plane. Blade row interactions are neglected in this formulation, but the force distributions over the entire blade surface for both the rotor and stator are obtained. Results from the hub-to-casing program are compared with a relaxation method solution for a subsonic rotor. Results are also presented for a quiet fan stage which includes transonic flow in both the rotor and stator and a normal shock in the stator.

Erdos, J.↗
NASA NTRS1972

Finite element analysis of large transient elastic-plastic deformations of simple structures, with application to the engine rotor fragment containment/deflection problem

Assumed-displacement versions of the finite-element method are developed to predict large-deformation elastic-plastic transient deformations of structures. Both the conventional and a new improved finite-element variational formulation are derived. These formulations are then developed in detail for straight-beam and curved-beam elements undergoing (1) Bernoulli-Euler-Kirchhoff or (2) Timoshenko deformation behavior, in one plane. For each of these categories, several types of assumed-displacement finite elements are developed, and transient response predictions are compared with available exact solutions for small-deflection, linear-elastic transient responses. The present finite-element predictions for large-deflection elastic-plastic transient responses are evaluated via several beam and ring examples for which experimental measurements of transient strains and large transient deformations and independent finite-difference predictions are available.

Wu, R. W.↗
NASA NTRS1974

Liquefaction dynamics of n-octadecane in cylindrical coordinates

The study presents one method of mathematically modeling the effects of gravity-induced free convection upon the melting phenomena of a finite paraffin slab. The theoretical model is based on the numerical computer solution of the solid phase and liquid phase energy equations, coupled with an imposed ideal-viscous flow approximation. An implicit alternating-direction technique is used as the finite-difference approximation for the basic energy equations in the numerical solution. The theoretical predictions are compared to experimental data from an annular test cell, in which n-octadecane, C18H38, was used as the test material. The results show that the model is a good first approximation to the solution of solid-liquid phase change involving gravity-induced free convection.

Bain, R. L.↗
NASA NTRS1975

A finite-difference program for stresses in anisotropic, layered plates in bending

The interlaminar stresses induced in a layered laminate that is bent into a cylindrical surface are studied. The laminate is modeled as a continuum, and the resulting elasticity equations are solved using the finite difference method. The report sets forth the mathematical framework, presents some preliminary results, and provides a listing and explanation of the computer program. Significant among the results are apparent symmetry relationships that will reduce the numerical size of certain problems and an interlaminar stress behavior having a sharp rise at the free edges.

Salamon, N. J.↗
NASA NTRS1975

An investigation of several numerical procedures for time-asymptotic compressible Navier-Stokes solutions

The status of an investigation of four numerical techniques for the time-dependent compressible Navier-Stokes equations is presented. Results for free shear layer calculations in the Reynolds number range from 1000 to 81000 indicate that a sequential alternating-direction implicit (ADI) finite-difference procedure requires longer computing times to reach steady state than a low-storage hopscotch finite-difference procedure. A finite-element method with cubic approximating functions was found to require excessive computer storage and computation times. A fourth method, an alternating-direction cubic spline technique which is still being tested, is also described.

Rudy, D. H.↗
NASA NTRS1978

Numerical solution of the Navier-Stokes equations for super-sonic flows with strong shocks

The numerical solution of the full Navier-Stokes Equations for viscous flows with high Mach numbers and a strong detached bow shock was obtained. Two dimensional flows around a circular cylinder, and a circular cylinder with an aft-body in the form of a fairing, were considered. The solution of the compressible N.S. equations was accomplished by the method of finite differences. An implicit scheme of solution, the S.O.R., was used with the optimum acceleration parameters determined by trial and error. The tensor notation was used in writing the N-S Equations transformed into general curvilinear coordinates. The equations for the generation of the coordinate system were solved, followed by the solution of the N.S. equations, at the end of a set of given number of time steps. "Wiggles", constituted the one major problem that needed to be overcome. These oscillations give rise to quantities such as negative temperatures, which ultimately caused the computational program to break down. Certain dissipative finite-difference schemes damped these oscillations.

Devarayalu, K.↗
NASA NTRS1978

Numerical method for boundary layers with blowing - The exponential box scheme

The paper describes a new numerical scheme based on exponential difference operator concepts combined with Keller's (1968) box scheme approach to produce a stable second-order accurate finite-difference scheme for convection-diffusion problems arising in boundary layer flows in the presence of massive injection through a porous surface. The technique is demonstrated by application to the self-similar boundary layer equations with massive blowing at the surface.

El-Mistikawy, T. M.↗
NASA NTRS1979

Application of numerical methods to heat transfer and thermal stress analysis of aerospace vehicles

The paper describes a thermal-structural design analysis study of a fuel-injection strut for a hydrogen-cooled scramjet engine for a supersonic transport, utilizing finite-element methodology. Applications of finite-element and finite-difference codes to the thermal-structural design-analysis of space transports and structures are discussed. The interaction between the thermal and structural analyses has led to development of finite-element thermal methodology to improve the integration between these two disciplines. The integrated thermal-structural analysis capability developed within the framework of a computer code is outlined.

Wieting, A. R.↗
NASA NTRS1982

Finite-difference solution of the compressible stability eigenvalue problem

A compressible stability analysis computer code is developed. The code uses a matrix finite difference method for local eigenvalue solution when a good guess for the eigenvalue is available and is significantly more computationally efficient than the commonly used initial value approach. The local eigenvalue search procedure also results in eigenfunctions and, at little extra work, group velocities. A globally convergent eigenvalue procedure is also developed which may be used when no guess for the eigenvalue is available. The global problem is formulated in such a way that no unstable spurious modes appear so that the method is suitable for use in a black box stability code. Sample stability calculations are presented for the boundary layer profiles of a Laminar Flow Control (LFC) swept wing.

Malik, M. R.↗
NASA NTRS1982

Incompressible Navier-Stokes and parabolized Navier-Stokes solution procedures and computational techniques

Recent developments with finite-difference techniques are emphasized. The quotation marks reflect the fact that any finite discretization procedure can be included in this category. Many so-called finite element collocation and galerkin methods can be reproduced by appropriate forms of the differential equations and discretization formulas. Many of the difficulties encountered in early Navier-Stokes calculations were inherent not only in the choice of the different equations (accuracy), but also in the method of solution or choice of algorithm (convergence and stability, in the manner in which the dependent variables or discretized equations are related (coupling), in the manner that boundary conditions are applied, in the manner that the coordinate mesh is specified (grid generation), and finally, in recognizing that for many high Reynolds number flows not all contributions to the Navier-Stokes equations are necessarily of equal importance (parabolization, preferred direction, pressure interaction, asymptotic and mathematical character). It is these elements that are reviewed. Several Navier-Stokes and parabolized Navier-Stokes formulations are also presented.

Rubin, S. G.↗