SEARCH · Search NASA

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.

SourceResearch connectionFrom yearThrough year

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

At least 487 records · Page 27

All sources NASA NTRS DOE OSTI

NASA NTRS1985

The numerical simulation of steady transonic rotational flow using a dual potential formulation

A finite-difference method is presented that simulates steady transonic rotational flow of an inviscid fluid by representing the velocity field as the sum of scalar and vector potentials. This dual potential velocity decomposition extends the validity of the scalar (full) velocity potential to include vorticity. The inclusion of a vector potential also permits an alternate treatment of lift that does not require a circulation wake cut. This is accomplished by specifying the vector potential as a constant on the airfoil surface in order to satisfy a Kutta condition. The governing equations are solved as iteratively decoupled scalar equations using approximate factorization techniques, and the overall efficiency approaches that of the full potential equation. The governing equations are able to convect entropy and vorticity throughout the flow field and are equivalent to the Euler equations in continuous flow domains, however at shocks the Rankine-Hugoniot entropy jump must be supplied. An entropy correction method is presented and verified with transonic airfoil solutions of the Euler equations.

Chaderjian, N. M.↗

NASA NTRS1985

Full two-dimensional transient solutions of electrothermal aircraft blade deicing

Two finite difference methods are presented for the analysis of transient, two-dimensional responses of an electrothermal de-icer pad of an aircraft wing or blade with attached variable ice layer thickness. Both models employ a Crank-Nicholson iterative scheme, and use an enthalpy formulation to handle the phase change in the ice layer. The first technique makes use of a 'staircase' approach, fitting the irregular ice boundary with square computational cells. The second technique uses a body fitted coordinate transform, and maps the exact shape of the irregular boundary into a rectangular body, with uniformally square computational cells. The numerical solution takes place in the transformed plane. Initial results accounting for variable ice layer thickness are presented. Details of planned de-icing tests at NASA-Lewis, which will provide empirical verification for the above two methods, are also presented.

Masiulaniec, K. C.↗

NASA NTRS1985

A fully implicit scheme for the barotropic primitive equations

An efficient implicit finite-difference method is developed and tested for a global barotropic model. The scheme requires, at each time step, the solution of only one-dimensional block-tridiagonal linear systems. This additional computation is offset by the use of a time step chosen independently of the mesh spacing. The method is second-order accurate in time and fourth-order accurate in space. Present experience indicates that this implicit method is practical for numerical simulation on fine meshes.

Cohn, S. E.↗

NASA NTRS1985

Large deflections of circular isotropic membranes subjected to arbitrary axisymmetric loading

Circular membranes with fixed peripheral edges, subjected to arbitrary axisymmetric loading were analyzed. A single governing differential equation in terms of radial stress was used. This nonlinear governing equation was solved using the finite difference method in conjunction with Newton-Raphson method. Three loading cases, namely: (1) uniformly loaded membrane, (2) a membrane with uniform load over an inner portion, and (3) a membrane with ring load, were analyzed. Calculated central displacement and the central and edge radial stresses for uniformly loaded membrane, agreed extremely well with the classical solution.

Kelkar, A.↗

NASA NTRS1986

Multiple Grids in Finite-Difference Flow Analysis

Multiple solutions superimposed to resolve flows about complex configurations. Use of multiple, overset grids in computational fluid dynamics extends application of finite-difference methods to more complex configurations. Rather than trying to generate single mesh about all components of configuration, multiple, individual meshes used, then overset on major grid. Major grid used to resolve flow field or wrapped around main component.

Dougherty, F. C.↗

NASA NTRS1985

Algebraic grid generation using tensor product B-splines

Finite difference methods are more successful if the accompanying grid has lines which are smooth and nearly orthogonal. The development of an algorithm which produces such a grid when given the boundary description. Topological considerations in structuring the grid generation mapping are discussed. The concept of the degree of a mapping and how it can be used to determine what requirements are necessary if a mapping is to produce a suitable grid is examined. The grid generation algorithm uses a mapping composed of bicubic B-splines. Boundary coefficients are chosen so that the splines produce Schoenberg's variation diminishing spline approximation to the boundary. Interior coefficients are initially chosen to give a variation diminishing approximation to the transfinite bilinear interpolant of the function mapping the boundary of the unit square onto the boundary grid. The practicality of optimizing the grid by minimizing a functional involving the Jacobian of the grid generation mapping at each interior grid point and the dot product of vectors tangent to the grid lines is investigated. Grids generated by using the algorithm are presented.

Saunders, B. V.↗

NASA NTRS1985

Numerical solutions of Navier-Stokes equations for a Butler wing

The flow field is simulated on the surface of a given delta wing (Butler wing) at zero incident in a uniform stream. The simulation is done by integrating a set of flow field equations. This set of equations governs the unsteady, viscous, compressible, heat conducting flow of an ideal gas. The equations are written in curvilinear coordinates so that the wing surface is represented accurately. These equations are solved by the finite difference method, and results obtained for high-speed freestream conditions are compared with theoretical and experimental results. In this study, the Navier-Stokes equations are solved numerically. These equations are unsteady, compressible, viscous, and three-dimensional without neglecting any terms. The time dependency of the governing equations allows the solution to progress naturally for an arbitrary initial initial guess to an asymptotic steady state, if one exists. The equations are transformed from physical coordinates to the computational coordinates, allowing the solution of the governing equations in a rectangular parallel-piped domain. The equations are solved by the MacCormack time-split technique which is vectorized and programmed to run on the CDC VPS 32 computer.

Abolhassani, J. S.↗

NASA NTRS1985

Interface procedures for overlapping grids

Interpolation at grid boundaries is studied for the purpose of solving partial differential equations using either implicit or conservative explicit finite-difference methods on multi-component overlapping grid systems.

Mastin, C. W.↗

NASA NTRS1984

User's guide for NASCRIN: A vectorized code for calculating two-dimensional supersonic internal flow fields

A computer program NASCRIN has been developed for analyzing two-dimensional flow fields in high-speed inlets. It solves the two-dimensional Euler or Navier-Stokes equations in conservation form by an explicit, two-step finite-difference method. An explicit-implicit method can also be used at the user's discretion for viscous flow calculations. For turbulent flow, an algebraic, two-layer eddy-viscosity model is used. The code is operational on the CDC CYBER 203 computer system and is highly vectorized to take full advantage of the vector-processing capability of the system. It is highly user oriented and is structured in such a way that for most supersonic flow problems, the user has to make only a few changes. Although the code is primarily written for supersonic internal flow, it can be used with suitable changes in the boundary conditions for a variety of other problems.

Kumar, A.↗

NASA NTRS1986

Navier-Stokes computations for circulation controlled airfoils

Navier-Stokes computations of subsonic to transonic flow past airfoils with augmented lift due to rearward jet blowing over a curved trailing edge are presented. The approach uses a spiral grid topology. Solutions are obtained using a Navier-Stokes code which employs an implicit finite difference method, an algebraic turbulence model, and developments which improve stability, convergence, and accuracy. Results are compared against experiments for no jet blowing and moderate jet pressures and demonstrate the capability to compute these complicated flows.

Pulliam, T. H.↗

NASA NTRS1986

Wall reflection of a viscous vortex ring

The behavior of a viscous axisymmetric vortex ring being reflected from a wall is investigated. The incompressible Navier-Stokes equations formulated in terms of the vorticity function and vector potential are numerically integrated by implicit finite difference methods. To specify the vector potential at a far boundary from the wall, the existing integral method used so far only for an unbounded domain is modified by a kind of image method. The trajectory of the vortex ring calcualted as a result closely resembles that observable from the experiment.

Sa, J. Y.↗

NASA NTRS1985

Principles of magnetohydrodynamic simulation in space plasmas

Attention is given to the philosophical as well as physical principles that are essential to the establishment of MHD simulation studies for solar plasma research, assuming the capabilities of state-of-the-art computers and emphasizing the importance of 'local' MHD simulation. Solar-terrestrial plasma space is divided into several elementary regions where a macroscopic elementary energy conversion process could conceivably occur; the local MHD simulation is defined as self-contained in each of the regions. The importance of, and the difficulties associated with, the boundary condition are discussed in detail. The roles of diagnostics and of the finite difference method are noted.

Sato, T.↗

NASA NTRS1986

Finite difference simulations of seismic scattering Implications for the propagation of short-period seismic waves in the crust and models of crustal heterogeneity

Synthetic seismographs that were obtained by the finite difference method are presently applied to the study of elastic and acoustic wave scattering in two-dimensional media with random spatial variations in seismic velocity. The seismograms are analyzed to determine the variation in travel times and waveforms across arrays of receivers. The random media with Gaussian and exponential correlation functions considered differ in the spectral falloff of their velocity fluctuations at wavelengths smaller than 2pi times the correlation distance. It is found that alternative models of crustal heterogeneity can be tested by improved measurements of the frequency dependence of the crustal Q at frequencies greater than about 1 Hz, assuming that scattering is responsible for most of the attenuation at such frequencies.

Frankel, A.↗

NASA NTRS1987

Analysis of Scramjet Inlets

NASCRIN program analyzes two-dimensional flow fields in supersoniccombustion ramjet (scramjet) inlets. Solves two-dimensional Euler or Navier-Stokes equations in conservative form by unsplit, explicit, two-step finite-difference method. More recent explicit/implicit, two-step scheme incorporated by analysis of viscous flow. Algebraic, two-layer eddy-viscosity model used for turbulent-flow calculations. Vectorized version, written for CDC CYBER 205, whereas scalar version, can be run on CRAY or other scalar computers.

Kumar, A.↗

NASA NTRS1987

Calculating rotordynamic coefficients of seals by finite-difference techniques

For modelling the turbulent flow in a seal the Navier-Stokes equations in connection with a turbulence (kappa-epsilon) model are solved by a finite-difference method. A motion of the shaft round the centered position is assumed. After calculating the corresponding flow field and the pressure distribution, the rotor-dynamic coefficients of the seal can be determined. These coefficients are compared with results obtained by using the bulk flow theory of Childs and with experimental results.

Dietzen, F. J.↗

NASA NTRS1986

Finite-element reentry heat-transfer analysis of space shuttle Orbiter

A structural performance and resizing (SPAR) finite-element thermal analysis computer program was used in the heat-transfer analysis of the space shuttle orbiter subjected to reentry aerodynamic heating. Three wing cross sections and one midfuselage cross section were selected for the thermal analysis. The predicted thermal protection system temperatures were found to agree well with flight-measured temperatures. The calculated aluminum structural temperatures also agreed reasonably well with the flight data from reentry to touchdown. The effects of internal radiation and of internal convection were found to be significant. The SPAR finite-element solutions agreed reasonably well with those obtained from the conventional finite-difference method.

Ko, William L.↗

NASA NTRS1986

Dislocation dynamics during the growth of silicon ribbon

The thermal viscoplastic stresses and the dislocation densities in silicon ribbon are computed for an axially changing thermal profile by using an iterative finite difference method. A material constitutive equation (Haasen-Sumino model) which involves an internal variable (mobile dislocation density) is used. The results are interpreted as showing that there is a maximum width of silicon ribbon that can be grown when viscoplasticity and dislocations are considered. This maximum width limitation does not exist if the material behavior is elastic.

Dillon, O. W., Jr.↗

NASA NTRS1987

A new PNS code for chemical nonequilibrium flows

A new parabolized Navier-Stokes (PNS) code has been developed to compute the hypersonic laminar flow of a multicomponent, chemically reacting mixture of thermally perfect gases over two-dimensional and axisymmetric bodies. The new PNS code solves the gas dynamic and species conservation equations in a coupled manner using a noniterative, implicit, space-marching finite-difference method. The conditions for well-posedness of the space-marching method have been derived from an eigenvalue analysis of the governing equations. The code has been used to compute hypersonic laminar flow of chemically reacting air over wedges and cones. The results of these computations are in good agreement with the results of reacting boundary-layer calculations.

Prabhu, D. K.↗