SEARCH · Engineering Papers

Results for “numerical method”

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 1,675 records · Page 93

All sources NASA NTRS DOE OSTI

NASA NTRS1985

Shape design sensitivity analysis using domain information

A numerical method for obtaining accurate shape design sensitivity information for built-up structures is developed and demonstrated through analysis of examples. The basic character of the finite element method, which gives more accurate domain information than boundary information, is utilized for shape design sensitivity improvement. A domain approach for shape design sensitivity analysis of built-up structures is derived using the material derivative idea of structural mechanics and the adjoint variable method of design sensitivity analysis. Velocity elements and B-spline curves are introduced to alleviate difficulties in generating domain velocity fields. The regularity requirements of the design velocity field are studied.

Seong, Hwal-Gyeong↗

NASA NTRS1985

The design and analysis of simple low speed flap systems with the aid of linearized theory computer programs

The purpose here is to show how two linearized theory computer programs in combination may be used for the design of low speed wing flap systems capable of high levels of aerodynamic efficiency. A fundamental premise of the study is that high levels of aerodynamic performance for flap systems can be achieved only if the flow about the wing remains predominantly attached. Based on this premise, a wing design program is used to provide idealized attached flow camber surfaces from which candidate flap systems may be derived, and, in a following step, a wing evaluation program is used to provide estimates of the aerodynamic performance of the candidate systems. Design strategies and techniques that may be employed are illustrated through a series of examples. Applicability of the numerical methods to the analysis of a representative flap system (although not a system designed by the process described here) is demonstrated in a comparison with experimental data.

Carlson, Harry W.↗

NASA NTRS1985

High-latitude filtering in a global grid-point model using model normal modes

A normal modes expansion technique is applied to perform high latitude filtering in the GLAS fourth order global shallow water model with orography. The maximum permissible time step in the solution code is controlled by the frequency of the fastest propagating mode, which can be a gravity wave. Numerical methods are defined for filtering the data to identify the number of gravity modes to be included in the computations in order to obtain the appropriate zonal wavenumbers. The performances of the model with and without the filter, and with a time tendency and a prognostic field filter are tested with simulations of the Northern Hemisphere winter. The normal modes expansion technique is shown to leave the Rossby modes intact and permit 3-5 day predictions, a range not possible with the other high-latitude filters.

Takacs, L. L.↗

NASA NTRS1985

Numerical solution of the Euler equation for compressible inviscid fluids

An effort is made to develop a satisfactory numerical method for the calculation of steady solutions of the Euler equations for inviscid compressible gas flows. The intended application is the prediction of the aerodynamic properties of aircraft flying at transonic speeds. Particular consideration is given to the modification of the equations to improve convergence to a steady state; finite volume formulation; adaptive dissipation; schemes designed to improve the resolution of shock waves; hybrid multistage time stepping schemes; residual averaging; a multigrid scheme; and trials with Burgers' equation. Results clearly demonstrate that the convergence of a time-dependent hyperbolic system to a steady state can be substantially accelerated by the introduction of multiple grids.

Jameson, A.↗

NASA NTRS1986

Connections between conventional and singular-value-based multi-variable flight control system design techniques

It is proposed that frequency-domain multivariable robustness techniques, when combined with classical multivariable procedures, can offer an additional means of evaluating FCS designs. A lateral-directional FCS for an advanced fighter is used as an example. Robustness to unstructured aircraft-input uncertainties is assessed using purely numerical singular-value procedures. Literal approximations for the singular values of the open-loop plant and controller and for the inverse return difference are shown to provide a means of decomposing and diagnosing robustness problems that are insoluble via purely numerical methods.

Mcruer, D. T.↗

NASA NTRS1987

Numerical simulation of transition, compressible turbulence, and reacting flows

Some of the recent work at NASA Langley on transition, turbulence, and reacting flows is summarized. Much of this effort has been motivated by outstanding technological problems in high-speed flow. A class of numerical algorithms for these problems has been developed and a variety of physical problems have been simulated. Descriptions are provided of the basic mathematical models, the nature of the numerical methods, and some of the recent simulations.

Zang, T. A.↗

NASA NTRS1987

An L-U implicit multigrid algorithm for the three-dimensional Euler equations

An LU implicit multigrid scheme is developed for the calculation of three-dimensional transonic flow through rotating cascades. This numerical method solves the unsteady Euler equations of gas dynamics in a finite-volume form. The implicit scheme makes it possible to take a much larger time step than is normally permitted in most explicit schemes, while the multigrid method is incorporated to accelerate the convergence rate for steady state calculations. Using this method, computational storage requirements are comparable to those of explicit schemes, while operation counts are considerably less than those found in the more widely-used ADI schemes.

Yokota, Jeffrey W.↗

NASA NTRS1987

Multigrid acceleration of the isenthalpic form of the compressible flow equations

A numerical method for solving the isenthalpic form of the governing equations for compressible inviscid flows was developed. The method is based on the concept of flux vector splitting in its implicit form and was tested on several demanding configurations. Time marching to steady state was accelerated by the implementation of the multigrid procedure which very effectively increased the rate of convergence. High quality steady-state results were obtained for various test cases and required only short computational times due to the relative efficiency of the basic method.

Melson, N. Duane↗

NASA NTRS1987

Free-wake analysis of a rotor in hover

A numerical method based on the axisymmetric, incompressible Navier-Stokes equations is combined with a lifting surface code to predict the vortex wake of hovering rotors. The lifting surface code, AMI Hover, is used to obtain the circulation distribution on the blade. This circulation distribution is fed into the Navier-Stokes code to compute the vortex wake under this specified circulation distribution. An iteration approach is used between these two codes to converge the circulation distribution and the shape of the vortex wake. A relaxation scheme is developed to resolve the instability encountered among the tip vortices. A reconcentration scheme is used to solve the diffusion problem due to the strong artificial viscosity. The results from the present method are compared with experimental data obtained by smoke-flow visualization and hot-wire measurements for several rotor blade configurations. The comparisons show that the present method is able to predict the complex wake system shed by a hovering rotor.

Chen, C. S.↗

NASA NTRS1987

Numerical solution of the radiation transport equation in disk geometry

An efficient numerical method for solving the problem of radiation transport in a dusty medium with two dimensional (2-D) disk geometry is described. It is a generalization of the one-dimensional quasi-diffusion method in which the transport equation is cast in diffusion form and then solved as a boundary value problem. The method should be applicable to a variety of astronomical sources, the dynamics of which are angular-momentum dominated and hence not accurately treated by spherical geometry, e.g., protoplanetary nebulae, circumstellar disks, interstellar molecular clouds, accretion disks, and disk galaxies. The computational procedure and practical considerations for implementing the method are described in detail. To illustrate the effects of 2-D radiation transport, some model results (dust temperature distributions and IR flux spectra) for externally heated, interstellar dust clouds with spherically symmetric and disk geometry are compared.

Spagna, George F., Jr.↗

NASA NTRS1987

Nonlinear steady-state response of a rotor-support system

A numerical method is presented to determine the steady-state nonlinear response of a rotor-support system due to deadband and rubbing using discrete Fourier transformation and inverse discrete Fourier transformation. Damaging subharmonic and superharmonic responses are found to occur in presence of a side force. The calculated results agree with the general trends which have been observed experimentally by other investigators. The effects of selected nondimensionalized parameters on rotor response are studied.

Noah, Sherif T.↗

NASA NTRS1987

A numerical study of transient, thermally-conductive solar wind

A numerical analysis of transient solar wind starting at the solar surface and arriving at 1 AU is performed by an implicit numerical method. The model hydrodynamic equations include thermal conduction terms for both steady and unsteady simulations. Simulation results show significant influence of thermal conduction on both steady and time-dependent solar wind. Higher thermal conduction results in higher solar wind speed, higher temperature, but lower plasma density at 1 AU. Higher base temperature at the solar surface gives lower plasma speed, lower temperature, but higher density at 1 AU. Higher base density, on the other hand, gives lower velocity, lower temperature, but higher density at 1 AU.

Han, S. M.↗

NASA NTRS1986

Numerical simulation of a controlled boundary layer

The problem of interest is the boundary layer over a flat plate. The three standard laminar flow control (LFC) techniques are pressure gradient, suction, and heating. The parameters used to describe the amount of control in the context of the boundary layer equations are introduced. The numerical method required to find the mean flow, the linear eigenvalues of the Orr-Sommerfeld equation, and the full, nonlinear, 3-D solution of the Navier-Stokes equations are outlined. A secondary instability exists for the parallel boundary subject to uniform pressure gradient, suction, or heating. Selective control of the spanwise mode reduces the secondary instability in the parallel boundary layer at low Reynolds number.

Zang, Thomas A.↗

NASA NTRS1987

Navier-Stokes calculations and turbulence modeling in the trailing edge region of a circulation control airfoil

The accurate prediction of turbulent flows over curved surfaces in general and over the trailing edge region of circulation control airfoils in particular requires the coupled efforts of turbulence modelers, numerical analysts and experimentalists. The purpose of the research program in this area is described. Then, the influence on turbulence modeling of the flow characteristics over a typical circulation control wing is discussed. Next, the scope of this effort to study turbulence in the trailing edge region of a circulation control airfoil is presented. This is followed by a brief overview of the computation scheme, including the grid, governing equations, numerical method, boundary conditions and turbulence models applied to date. Then, examples of applications of two algebraic eddy viscosity models to the trailing edge region of a circulation control airfoil is presented. The results from the calculations is summarized, and conclusions drawn based on examples. Finally, the future directions of the program is outlined.

Viegas, John R.↗

NASA NTRS1988

A constitutive material model for nonlinear finite element structural analysis using an iterative matrix approach

A unified numerical method for the integration of stiff time dependent constitutive equations is presented. The solution process is directly applied to a constitutive model proposed by Bodner. The theory confronts time dependent inelastic behavior coupled with both isotropic hardening and directional hardening behaviors. Predicted stress-strain responses from this model are compared to experimental data from cyclic tests on uniaxial specimens. An algorithm is developed for the efficient integration of the Bodner flow equation. A comparison is made with the Euler integration method. An analysis of computational time is presented for the three algorithms.

Koenig, Herbert A.↗

NASA NTRS1988

A decade of computer simulations for space shuttle aerodynamics

Ten years ago computer simulations of the flow field around the Space Shuttle Orbiter were limited to inviscid calculations for the windward side of the forebody and viscous calculations for selected 2-D problems. Advances in computer hardware and numerical methods during the past ten years made it possible to calculate viscous flow over the complete orbiter configuration at angle of attack. The equations solved are the Reynolds-averaged, Navier-Stokes equation, simplified by either the thin-layer or parabolized approximation. An algebraic eddy viscosity model is used for turbulent flow. The free stream is assumed to be a perfect gas for wind tunnel conditions and a real gas in thermodynamic equilibrium for flight conditions. Four examples of recent computer simulations are presented. Flow field results include oil flow patterns on the surface and Mach number contours, isobars, and cross-flow velocity vectors in the shock layer.

Inouye, Mamoru↗

NASA NTRS1988

Numerical solutions of the complete Navier-Stokes equations

The physical phenomena within supersonic flows that sustain chemical reactions are investigated. An earlier study to develop accurate physical models for supersonic reacting flowfields focused on 2-D laminar shear layers. The objective is to examine the mixing and subsequent combustion within turbulent reacting shear layers. To conduct this study, a computer program has been written to solve the axisymmetric Reynolds averaged Navier-Stokes equations. The numerical method uses a cell-centered finite volume approach and a Runge Kutta time stepping scheme. The Reynolds averaged equations are closed using the eddy viscosity concept. Several zero-equation models have been tested by making calculations for an H2-air nonreacting coaxial jet flow. Comparisons made with experimental data show that Cohen's eddy viscosity model provides best agreement. The finite rate chemistry model used in the study of 2-D laminar shear layers is incorporated into the computer program and data is compared from a recent experiment performed at NASA Langley.

Hassan, H. A.↗

NASA NTRS1988

Three-dimensional adaptive grid generation for body-fitted coordinate system

This report describes a numerical method for generating 3-D grids for general configurations. The basic method involves the solution of a set of quasi-linear elliptic partial differential equations via pointwise relaxation with a local relaxation factor. It allows specification of the grid spacing off the boundary surfaces and the grid orthogonality at the boundary surfaces. It includes adaptive mechanisms to improve smoothness, orthogonality, and flow resolution in the grid interior.

Chen, S. C.↗