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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 829 records · Page 46

All sources NASA NTRS DOE OSTI

NASA NTRS1991

Computational results for 2-D and 3-D ramp flows with an upwind Navier-Stokes solver

An implicit, finite-difference, upwind, full Navier-Stokes solver was applied to supersonic/hypersonic flows over two-dimensional ramps and three-dimensional obstacle. Some of the computed results are presented. The numerical scheme used in the study is an implicit, spatially second order accurate, upwind, LU-ADI scheme based on Roe's approximate Reimann solver with MUSCL differencing of Van Leer. An algebraic grid generation scheme based on generalized interpolation scheme was used in generating the grids for the various 2-D and 3-D problems.

Venkatapathy, Ethiraj↗

NASA NTRS1991

Predicting the oxidative lifetime of beta NiAl-Zr alloys

Nickel aluminides containing 40 to 50 at. pct Al and 0.1 at. pct Zr were studied following cyclic oxidation at 1400 C. The selective oxidation of Al resulted in the formation of protective Al2O3 scales on each alloy composition. However, repeated cycling eventually resulted in the gradual formation of less-protective NiAl2O4, first appearing on the 40Al alloys followed at longer times on the 45Al alloys. The appearance of the NiAl2O4, signaling the end of the protective scale-forming capability of the alloy, was related to the presence of gamma-prime (Ni3Al) which formed as a result of the loss of Al from the sample. A diffusion model, based on finite-difference techniques, was developed to predict the protective life of beta Ni-Al alloys. This model predicts Ni and Al concentration profiles after various oxidation exposures. The model can predict the oxidative lifetime due to Al depletion when the Al concentration decreases to a critical concentration. Measured Al concentration profiles on two alloys after various oxidation exposures are compared to those predicted by the diffusion model. The time to the appearance of the NiAl2O4 and that predicted by the diffusion model are compared and discussed.

Nesbitt, J. A.↗

NASA NTRS1991

Simulations of curved turbulent boundary layers

The objective of this work is to develop a space-time accurate numerical method for the solution of incompressible Navier-Stokes equations in generalized coordinates. The resulting code is to be used for direct and large-eddy simulation of turbulence in complex geometries. In a previous paper, the system of Navier-Stokes equations in general curvilinear coordinates was solved by a second-order accurate finite-difference scheme. Satisfactory results were obtained for several flows in two and three dimensions. The system of Navier-Stokes for the fluxes are given in Orlandi (1989). The main deficiency of the numerical scheme was the large CPU time required for the solution of the Poisson equation for the 'pressure' field. The point SOR relaxation, in conjunction with a multigrid scheme, was used for the Poisson equation. In some cases, particularly with very fine grids, it was impossible to obtain a divergent-free flow. A preliminary attempt is made to compute the spatially evolving flow of Swearingen & Blackwelder. To reduce the streamwise distance, the inflow was at a distance x = 60 cm from the leading edge.

Orlandi, Paolo↗

NASA NTRS1947

Determination of elastic stresses in gas-turbine disks

A method is presented for the calculation of elastic stresses in symmetrical disks typical of those of a high-temperature gas turbine. The method is essentially a finite-difference solution of the equilibrium and compatibility equations for elastic stresses in a symmetrical disk. Account can be taken of point-to-point variations in disk thickness, in temperature, in elastic modulus, in coefficient of thermal expansion, in material density, and in Poisson's ratio. No numerical integration or trial-and-error procedures are involved and the computations can be performed in rapid and routine fashion by nontechnical computers with little engineering supervision. Checks on problems for which exact mathematical solutions are known indicate that the method yields results of high accuracy. Illustrative examples are presented to show the manner of treating solid disks, disks with central holes, and disks constructed either of a single material or two or more welded materials. The effect of shrink fitting is taken into account by a very simple device.

Manson, S S↗

NASA NTRS1950

Direct method of design and stress analysis of rotating disks with temperature gradient

A method is presented for the determination of the contour of disks, typified by those of aircraft gas turbines, to incorporate arbitrary elastic-stress distributions resulting from either centrifugal or combined centrifugal and thermal effects. The specified stress may be radial, tangential, or any combination of the two. Use is made of the finite-difference approach in solving the stress equations, the amount of computation necessary in the evolution of a design being greatly reduced by the judicious selection of point stations by the aid of a design chart. Use of the charts and of a preselected schedule of point stations is also applied to the direct problem of finding the elastic and plastic stress distribution in disks of a given design, thereby effecting a great reduction in the amount of calculation. Illustrative examples are presented to show computational procedures in the determination of a new design and in analyzing an existing design for elastic stress and for stresses resulting from plastic flow.

Manson, S S↗

NASA NTRS1990

Constrained optimization using design of experiment surfaces

An algorithm for solving constrained optimization problems is presented. First, design of experiment techniques are used to survey the design space. After evaluating the objective and constraint functions, as specified by Taguchi orthogonal arrays, analytical models of these functions are generated using a least-squares regression analysis. Next, a nonlinear programming package is used to optimize the analytical model. Based on the optimization information, the design space is reduced so as to close in around the minimum, and the entire procedure is repeated until convergence. An important feature of the algorithm is that function gradients are not required; therefore, for problems in which gradients would have to be estimated using finite-differences the number of function evaluations required for the optimization is significantly reduced, when compared with traditional nonlinear programming techniques. In addition, there is no requirement that the gradients must be smooth and continuous.

Bolt, Marvin Vance↗

NASA NTRS1992

Numerical simulation of entry flow over blunt swept-wing planes

Viscous, heat-conducting flow with chemical and vibrational relaxation processes of the constituent gases surrounding winged spacecraft is considered in the continuum regime. The Navier-Stokes equations are appended by additional vibrational energy and species rate equations and supplemented by the equations of state and the phenomenological laws based on mixture rules or collisional cross sections. Numerical convective flux can be obtained from several forms of one-dimensional Riemann solver, with or without entropy correction. High-order accuracy is obtained from two types of reconstructive interpolation. A number of explicit and implicit numerical schemes have been implemented as a means to yield converged solutions. Both shock-fitting, finite-difference and shock-capturing, finite-volume techniques have been tested for configurations such as a sphere, double ellipsoid, blunt-edge delta wing, a European Hermes vehicle, and the U.S. Shuttle Orbiter. The shock-fitting code provides excellent results only for simple configurations, whereas the shock-capturing code leads to overall satisfying solutions for complex geometries.

Li, C. P.↗

NASA NTRS1985

An Algorithm for Efficient Maximum Likelihood Estimation and Confidence Interval Determination in Nonlinear Estimation Problems

An algorithm for maximum likelihood (ML) estimation is developed with an efficient method for approximating the sensitivities. The algorithm was developed for airplane parameter estimation problems but is well suited for most nonlinear, multivariable, dynamic systems. The ML algorithm relies on a new optimization method referred to as a modified Newton-Raphson with estimated sensitivities (MNRES). MNRES determines sensitivities by using slope information from local surface approximations of each output variable in parameter space. The fitted surface allows sensitivity information to be updated at each iteration with a significant reduction in computational effort. MNRES determines the sensitivities with less computational effort than using either a finite-difference method or integrating the analytically determined sensitivity equations. MNRES eliminates the need to derive sensitivity equations for each new model, thus eliminating algorithm reformulation with each new model and providing flexibility to use model equations in any format that is convenient. A random search technique for determining the confidence limits of ML parameter estimates is applied to nonlinear estimation problems for airplanes. The confidence intervals obtained by the search are compared with Cramer-Rao (CR) bounds at the same confidence level. It is observed that the degree of nonlinearity in the estimation problem is an important factor in the relationship between CR bounds and the error bounds determined by the search technique. The CR bounds were found to be close to the bounds determined by the search when the degree of nonlinearity was small. Beale's measure of nonlinearity is developed in this study for airplane identification problems; it is used to empirically correct confidence levels for the parameter confidence limits. The primary utility of the measure, however, was found to be in predicting the degree of agreement between Cramer-Rao bounds and search estimates.

Murphy, Patrick Charles↗

NASA NTRS1980

Fast, Conservative Algorithm for Solving the Transonic Full-Potential Equation

A fast, fully implicit approximate factorization algorithm designed to solve the conservative, transonic, full-potential equation in either two or three dimensions is described. The algorithm uses an upwind bias of the density coefficient for stability in supersonic regions. This provides an effective upwind difference of the streamwise terms for any orientation of the velocity vector (i.e., rotated differencing), thereby greatly enhancing the reliability of the present algorithm. A numerical transformation is used to establish an arbitrary body-fitted, finite-difference mesh. Computed results for both airfoils and simplified wings demonstrate substantial improvement in convergence speed for the new algorithm relative to standard successive-line over-relaxation algorithms.

Holst, Terry L.↗

NASA NTRS1981

Two-Dimensional Inlet Simulation Using a Diagonal Implicit Algorithm

A modification of an implicit approximate-factorization finite-difference algorithm applied to the two-dimensional Euler and Navier-Stokes equations in general curvilinear coordinates is presented for supersonic freestream flow about and through inlets. The modification transforms the coupled system of equations Into an uncoupled diagonal form which requires less computation work. For steady-state applications the resulting diagonal algorithm retains the stability and accuracy characteristics of the original algorithm. Solutions are given for inviscid and laminar flow about a two-dimensional wedge inlet configuration. Comparisons are made between computed results and exact theory.

Chaussee, D.S.↗

NASA NTRS1961

Cross-Sectional Deformations of Monocoque Beams and Their Effects on the Natural Vibration Frequencies

The variational principle, differential equations, and boundary conditions governing the cross-sectional distortions due to inertia loading of a two-dimensional model of a thin monocoque wing are shown. A theoretical analysis of this simplified model is made in order to determine the nature of the coupling between the cross-sectional modes and the spanwise deformation modes. General solutions are obtained in finite-difference form for arbitrary cross sections and an exact solution is presented for a parabolic-arc cross section of constant cover thickness. The application of these results in evaluating the coupled frequencies of the actual structure is discussed. Frequencies evaluated for a parabolic-arc monocoque beam show good agreement with experimental values.

Thomson, Robert G.↗

NASA NTRS1981

Calculation of Unsteady Transonic Flow over an Airfoil

An implicit finite-difference solver for either the Euler equations or the "thin-layer" Navier-Stokes equations was used to calculate a transonic flow over the NACA 64A010 airfoil pitching about its one-quarter chord. An unsteady automatic grid-generation procedure that will improve significantly the computational efficiency of various unsteady flow problems is described. The calculated results for both inviscid and viscous flows at Much number 0.8 over the airfoil oscillating with reduced frequency referenced to one-half chord, 0.2, are compared with experimental data measured in the Ames 11 x 11-ft Transonic Wind Tunnel. Nonlinear, unsteady effects of the flow on the surface pressure variations, shock-wave excursions, and overall airloads are examined. Good agreements between the results of computations and experiments were obtained. In the shock-wave region, however, the results of the viscous-flow computations showed closer agreement with the experimental data.

Chyu, W. J.↗

NASA NTRS1981

Calculation of Unsteady Transonic Flow Over an Airfoil

An implicit finite-difference solver for either the Euler equations or the "thin-layer" Navier-Stokes equations was used to calculate a transonic flow over the NACA 64A010 airfoil pitching about its one-quarter chord. An unsteady automatic grid-generation procedure that will improve significantly the computational efficiency of various unsteady flow problems is described. The calculated results for both inviscid and viscous flows at Mach number 0.8 over the airfoil oscillating with reduced frequency referenced to one-half chord, 0.2, are compared with experimental data measured in the Ames 11 x 11 ft Transonic Wind Tunnel. Nonlinear, unsteady effects of the flow on the surface pressure variations, shock-wave excursions, and overall airloads are examined. Good agreements between the results of computations and experiments were obtained. In the shock-wave region, however, the results of the viscous-flow computations showed closer agreement with the experimental data.

Chyu, W. J.↗

NASA NTRS1988

Correlation of SA349/2 Helicopter Flight Test Data with a Comprehensive Rotorcraft Model

A comprehensive rotorcraft analysis model, CAMRAD, was used to predict blade aerodynamic and structural loads and performance for comparison with flight test data. The data were obtained from an SA349/2 helicopter with an advanced geometry rotor. Sensitivity of the correlation to wake geometry, blade dynamics, and blade aerodynamic effects was investigated. Blade chordwise pressure coefficients were predicted for the transonic regimes near the blade tip using the model coupled with two finite-difference codes.

Yamauchi, Gloria K.↗

NASA NTRS1981

Nonreflecting Far-Field Boundary Conditions for Unsteady Transonic Flow Computation

The approximate nonreflecting far-field boundary condition, as proposed by Engquisi and Majda, is implemented In the computer code LTRAN2. This code solves the Implicit finite-difference representation of the small-disturbance equations for unsteady transonic flows about airfoils. The nonreflecting boundary condition and the description of the algorithm for Implementing these conditions In LTRAN2 tire discussed. Various cases re computed and compared with results from the older, more conventional procedures. One concludes that the nonreflecting far-field boundary approximation allows the far-field boundary to be located closer to the airfoil; this permits a decrease in the computer lime required to obtain the solution through the use of fewer mesh points.

Kwak, D.↗

NASA NTRS1991

Thermal Modeling and Analysis of a Cryogenic Tank Design Exposed to Extreme Heating Profiles

A cryogenic test article, the Generic Research Cryogenic Tank, was designed to qualitatively simulate the thermal response of transatmospheric vehicle fuel tanks exposed to the environment of hypersonic flight. One-dimensional and two-dimensional finite-difference thermal models were developed to simulate the thermal response and assist in the design of the Generic Research Cryogenic Tank. The one-dimensional thermal analysis determined the required insulation thickness to meet the thermal design criteria and located the purge jacket to eliminate the liquefaction of air. The two-dimensional thermal analysis predicted the temperature gradients developed within the pressure-vessel wall, estimated the cryogen boiloff, and showed the effects the ullage condition has on pressure-vessel temperatures. The degree of ullage mixing, location of the applied high-temperature profile, and the purge gas influence on insulation thermal conductivity had significant effects on the thermal behavior of the Generic Research Cryogenic Tank. In addition to analysis results, a description of the Generic Research Cryogenic Tank and the role it will play in future thermal structures and transatmospheric vehicle research at the NASA Dryden Flight Research Facility is presented.

Stephens, Craig A.↗

NASA NTRS1986

The Critical Role of Computational Fluid Dynamics in Rotary-Wing Aerodynamics

Helicopter rotor aerodynamics is basically the study of unsteady aerodynamic flows in a rotating and translating coordinate system. Current trends in this field are briefly reviewed by examining recent advances in lifting-surface theory, wake modeling, panel methods, and finite-difference models' Examples are used to illustrate selected current methods and some indications of promising future directions are highlighted.

Davis, Sanford S.↗

NASA NTRS1989

A Computational Fluid Dynamics Algorithm on a Massively Parallel Computer

The discipline of computational fluid dynamics is demanding ever-increasing computational power to deal with complex fluid flow problems. We investigate the performance of a finite-difference computational fluid dynamics algorithm on a massively parallel computer, the Connection Machine. Of special interest is an implicit time-stepping algorithm; to obtain maximum performance from the Connection Machine, it is necessary to use a nonstandard algorithm to solve the linear systems that arise in the implicit algorithm. We find that the Connection Machine ran achieve very high computation rates on both explicit and implicit algorithms. The performance of the Connection Machine puts it in the same class as today's most powerful conventional supercomputers.

Jespersen, Dennis C.↗