Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Finite-differences”

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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
NASA NTRS1987

Computation and analysis of a cylinder wake flow

The Karman vortex wake of a circular cylinder at low Reynolds number was computed by a time-accurate, two-dimensional compressible Navier-Stokes equation solver which uses the MacCormack predictor-corrector finite-difference scheme and a nonreflecting boundary condition on the outer flow boundary. The results from a large number of time steps were analyzed using Fast Fourier Transform techniques to identify the important frequency components for comparison with published experimental data. A strong low-frequency component was found below the vortex shedding frequency and not harmonically related to it. The experimentally discovered low-frequency fluctuations in the cylinder wake are considered possibly to be precursors to transition from laminar to turbulent flow conditions. The present finding of similar frequencies in a computed wake tends to confirm their existence as a real wake phenomenon. This computational work provides a complementary means to experimental investigations of wake phenomena.

Townsend, J. C.↗
NASA NTRS1987

Steady and transient least square solvers for thermal problems

This paper develops a hierarchical least square solution algorithm for highly nonlinear heat transfer problems. The methodology's capability is such that both steady and transient implicit formulations can be handled. This includes problems arising from highly nonlinear heat transfer systems modeled by either finite-element or finite-difference schemes. The overall procedure developed enables localized updating, iteration, and convergence checking as well as constraint application. The localized updating can be performed at a variety of hierarchical levels, i.e., degree of freedom, substructural, material-nonlinear groups, and/or boundary groups. The choice of such partitions can be made via energy partitioning or nonlinearity levels as well as by user selection. Overall, this leads to extremely robust computational characteristics. To demonstrate the methodology, problems are drawn from nonlinear heat conduction. These are used to quantify the robust capabilities of the hierarchical least square scheme.

Padovan, Joe↗
NASA NTRS1987

Implicit finite difference methods on composite grids

Techniques for eliminating time lags in the implicit finite-difference solution of partial differential equations are investigated analytically, with a focus on transient fluid dynamics problems on overlapping multicomponent grids. The fundamental principles of the approach are explained, and the method is shown to be applicable to both rectangular and curvilinear grids. Numerical results for sample problems are compared with exact solutions in graphs, and good agreement is demonstrated.

Mastin, C. Wayne↗
NASA NTRS1987

Impact on the earth, ocean and atmosphere

On the basis of finite-difference techniques, cratering flow calculations are used to obtain the spatial attenuation of shock pressure with radius along the impact axis for the impact of silicate rock and iron impactors on a silicate half-space at speeds of 5 to 45 km/sec. Upon impact of a 10 to 30 km diameter silicate or water object onto a 5 km deep ocean overlying a silicate half-space planet at 30 km/sec, it is found that from 12 to 15 percent of the incident energy is coupled into the water. The mass of atmosphere lost due to impacts of 1 to 5 km radius projectiles is calculated.

Ahrens, Thomas J.↗
NASA NTRS1988

The calculation of flow over iced airfoils

Progress toward the development of a method for predicting the flowfield of an iced airfoil is described and shown to offer the prospect of a priori calculations of the effects of ice accretion and roughness on airfoil performance. The approach is based on interaction of inviscid flow solutions obtained by a panel method and improved upon by a finite-difference boundary-layer method which, operating in an inverse mode, incorporates viscous effects including those associated with separated flows. Results are presented for smooth, rough and iced airfoils as a function of angle of attack. Those for smooth and rough airfoils confirm the accuracy of the method and its applicability to surfaces with roughness similar to that associated with insect deposition and some forms of ice. Two procedures have been developed to deal with large ice accretion and their performance is examined and shown to be appropriate to the engineering requirements.

Cebeci, Tuncer↗
NASA NTRS1988

Multi-airfoil Navier-Stokes simulations of turbine rotor-stator interaction

An accurate numerical analysis of the flows associated with rotor-stator configurations in turbomachinery can be extremely helpful in optimizing the performance of turbomachinery. In this study the unsteady, thin-layer, Navier-Stokes equations in two spatial dimensions are solved on a system of patched and overlaid grids for a rotor-stator configuration from an axial turbine. The governing equations are solved using a finite-difference, upwind algorithm that is set in an iterative, implicit framework. Results in the form of pressure contours, time-averaged pressures, unsteady pressure amplitudes and phase are presented. The numerical results are compared with experimental data and the agreement is found to be good. The results are also compared with those of an earlier study which used only one rotor and one stator. The current study uses multiple rotors and stators and a pitch ratio that is much closer to the experimental ratio. Consequently the results of this study are found to be closer to the experimental data.

Rai, Man Mohan↗
NASA NTRS1987

Prediction of SA 349/2 GV blade loads in high speed flight using several rotor analyses

The influence of blade dynamics, dynamic stall, and transonic aerodynamics on the predictions of rotor loads in high-speed flight are presented. Data were obtained from an Aerospatiale Gazelle SA 349/2 helicopter with three Grande Vitesse blades. Several analyses are used for this investigation. First, blade dynamics effects on the correlation are studied using three rotor analyses which differ mainly in the method of calculating the blade elastic response. Next, an ONERA dynamic stall model is used to predict retreating blade stall. Finally, advancing blade aerodynamic loads are calculated using a NASA-developed rotorcraft analysis coupled with two transonic finite-difference analyses.

Gaubert, Michel↗
NASA NTRS1987

Aerodynamic interaction between vortical wakes and lifting two-dimensional bodies

Unsteady rotor wake interactions with the empenage, tail boom, and other aerodynamic surfaces of a helicopter have a significant influence on its aerodynamic performance, the ride quality, and amount of vibration. A numerical method for computing the aerodynamic interaction between an interacting vortex wake and the viscous flow about arbitrary two-dimensional bodies has been developed to address this helicopter problem. The method solves for the flow field velocities on a body-fitted computational mesh using finite-difference techniques. The interaction of a rotor wake with the flow about a 4:1 elliptic cylinder at 45-deg incidence was calculated for a Reynolds number of 3000.

Stremel, Paul M.↗
NASA NTRS1987

Interdiffusion in Ni-rich, Ni-Cr-Al alloys at 1100 and 1200 C. I - Diffusion paths and microstructures. II - Diffusion coefficients and predicted concentration profiles

Interdiffusion in Ni-rich Ni-Cr-Al alloys is investigated experimentally after annealing at 1100 and 1200 C using gamma/gamma, gamma/gamma+beta, gamma/gamma+gamma prime, and gamma/gamma+alpha diffusion couples. The amount and location of Kirkendall porosity suggests that Al diffuses more rapidly than Cr which diffuses more rapidly than Ni in the gamma phase of Ni-Cr-Al alloys. The location and extent of maxima and minima in the concentration profiles of the diffusion couples indicate that both cross-term diffusion coefficients are positive. Measurements are also presented of the ternary interdiffusion coefficients of the gamma phase in the Ni-Cr-Al system. It is shown that the interdiffusion coefficients can be accurately predicted by using a ternary finite-difference interdiffusion model.

Nesbitt, J. A.↗
NASA NTRS1986

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↗
NASA NTRS1988

Implications of truncating semi-infinite physical domains on the accuracy of the solutions to the N-S equations

Numerical techniques for solving the Navier-Stokes equations of laminar incompressible flow over a forward-facing step on a semiinfinite flat plate are evaluated analytically and by means of sample computations. A primitive-variable method with pseudocompressibility and full mapping of the semiinfinite domain is compared with a vorticity/stream-function method with a truncated domain, applying a flux-corrected explicit finite-difference procedure to solve the discretized transport equations. The results are presented in graphs and compared with published experimental data. The solutions obtained with the truncated-domain methods are shown to be inaccurate, even in cases when they are smooth.

Dekruif, Jeffery S.↗
NASA NTRS1987

A numerical method for the solution of three dimensional, incompressible, viscous flows over slender bodies

A marching iterative method for the solution of the three dimensional, incompressibhle, steady and parabolized Navier-Stokes equations is described. The equations are written in primitive variables and discretized in general axisymmetric orthogonal coordinate systems. The coupled set of finite-difference equations are solved without any splitting or factorization errors. Moreover, the continuity equation and the two crossflow momentum equations are exactly satisfied at every step of the iterative process. The solution scheme is equivalent to the solution of one Poisson equation by the Successive Plane Over Relaxation method and has good convergence properties. Other existing solution methods resemble a Jacobi-type iterative scheme and therefore are less efficient. Numerical experiments include the laminar, incompressible flow over prolate spheroids at incidence.

Rosenfeld, Moshe↗
NASA NTRS1987

Wave envelope technique for multimode wave guide problems

A numerical solution technique for two-dimensional EM waveguide problems is developed analytically and demonstrated. The frequency-domain waveguide problem is transformed to a time-dependent problem and treated following the limiting-amplitude principle; the wave-envelope technique is applied; and the resulting problem is solved numerically using a finite-difference scheme based on a variation of the finite-area method (Baumeister, 1977 and 1984). Results for several sample problems, obtained in 10 min or less of execution time on PCs, are presented in tables.

Hariharan, S. I.↗
NASA NTRS1988

Mach number effects on transonic aeroelastic forces and flutter characteristics

Transonic aeroelastic stability analysis and flutter calculations are presented for a generic transport-type wing based on the use of the CAP-TSD (Computational Aeroelasticity Program - Transonic Small Disturbance) finite-difference code. The CAP-TSD code was recently developed for transonic unsteady aerodynamic and aeroelastic analysis of complete aircraft configurations. A binary aeroelastic system consisting of simple bending and torsion modes was used to study aeroelastic behavior at transonic speeds. Generalized aerodynamic forces are presented for a wide range of Mach number and reduced frequency. Aeroelastic characteristics are presented for variations in freestream Mach number, mass ratio, and bending-torsion frequency ratio. Flutter boundaries are presented which have two transonic dips in flutter speed. The first dip is the usual transonic dip involving a bending-dominated flutter mode. The second dip is characterized by a single degree-of-freedom torsion oscillation. These aeroelastic results are physically interpreted and shown to be related to the steady state shock location and changes in generalized aerodynamic forces due to freestream Mach number.

Mohr, Ross W.↗
NASA NTRS1988

An improved flux-split algorithm applied to hypersonic flows in chemical equilibrium

An explicit, finite-difference, shock-capturing numerical algorithm is presented and applied to hypersonic flows assumed to be in thermochemical equilibrium. Real-gas chemistry is either loosely coupled to the gasdynamics by way of a Gibbs free energy minimization package or fully coupled using species mass conservation equations with finite-rate chemical reactions. A scheme is developed that maintains stability in the explicit, finite-rate formulation while allowing relatively high time steps. The codes use flux vector splitting to difference the inviscid fluxes and employ real-gas corrections to viscosity and thermal conductivity. Numerical results are compared against existing ballistic range and flight data. Flows about complex geometries are also computed.

Palmer, Grant↗
NASA NTRS1988

Properties of finite difference models of non-linear conservative oscillators

Finite-difference (FD) approaches to the numerical solution of the differential equations describing the motion of a nonlinear conservative oscillator are investigated analytically. A generalized formulation of the Duffing and modified Duffing equations is derived and analyzed using several FD techniques, and it is concluded that, although it is always possible to contstruct FD models of conservative oscillators which are themselves conservative, caution is required to avoid numerical solutions which do not accurately reflect the properties of the original equation.

Mickens, R. E.↗
NASA NTRS1988

Transonic stability and control of aircraft using CFD methods

Implementation of a capability to calculate longitudinal short-period response in the CAP-TSD (Computational Aeroelasticity Program - Transonic Small Disturbance) finite-difference code is described. The code, developed recently at the NASA Langley Research Center, is capable of solving steady and unsteady flows about complete aircraft configurations and is used primarily for aeroelastic calculations in the critical transonic speed range. The longitudinal short-period equations of motion in state-space form have been coupled to the time-accurate lift and moment calculated by the program. Transient responses to an elevator pulse for free-flying aircraft demonstrate the new capability. A trim routine is also added to the code to obtain trim automatically during steady-state flow field convergence. Stability and control derivatives are estimated from the calculated transient response by a maximum likelihood estimation program. Results for a fighter configuration and a general aviation configuration are presented to assess the capability.

Vinh, Lam-Son↗
NASA NTRS1988

The effects of turbulence and stator/rotor interactions on turbine heat transfer. I - Design operating conditions

A combined experimental and analytical program was conducted to examine the effects of inlet turbulence, stator-rotor axial spacing, and relative circumferential spacing of first and second stators on turbine airfoil heat transfer. The experimental portion of the study was conducted in a large-scale (approximately 5X engine), ambient temperature, stage-and-a-half rotating turbine model. The data indicate that while turbine inlet turbulence can have a very strong impact on the first stator heat transfer, its impact in downstream rows is minimal. The effects on heat transfer produced by relatively large changes in stator/rotor spacing or by changing the relative row-to-row circumferential positions of stators were very small. Analytical results consist of airfoil heat transfer distributions computed with a finite-difference boundary layer code.

Blair, M. F.↗