Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “finite difference”

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 NTRS1984

Improved finite difference schemes for transonic potential calculations

Engquist and Osher (1980) have introduced a finite difference scheme for solving the transonic small disturbance equation, taking into account cases in which only compression shocks are admitted. Osher et al. (1983) studied a class of schemes for the full potential equation. It is proved that these schemes satisfy a new discrete 'entropy inequality' which rules out expansion shocks. However, the conducted analysis is restricted to steady two-dimensional flows. The present investigation is concerned with the adoption of a heuristic approach. The full potential equation in conservation form is solved with the aid of a modified artificial density method, based on flux biasing. It is shown that, with the current scheme, expansion shocks are not possible.

Hafez, M.↗
NASA NTRS1985

Propagation and stability of wavelike solutions of finite difference equations with variable coefficients

The propagation and dissipation of wavelike solutions to finite difference equations is analyzed on the basis of an asymptotic approach in which a wave solution is expressed as a product of a complex amplitude and an oscillatory phase function whose frequency and wavenumber may also be complex. An asymptotic expansion leads to a local dispersion relation for wavenumber and frequency; the first-order terms produce an equation for the amplitude in which the local group velocity appears as the convection velocity of the amplitude. Equations for the motion of wavepackets and their interaction at boundaries are derived, and a global stability analysis is carried out.

Giles, M. B.↗
NASA NTRS1987

A two-dimensional finite difference program for thermal analysis of rocket thrust chambers

A two-dimensional finite difference computer model for thermal analysis of rocket thrust chambers has been developed. The model uses an iterative scheme for calculating the temperature distribution within the chamber wall and implements a successive overrelaxation formula for a quick convergence. The inputs of the models are the dimensions of the thrust chamber wall, types of materials used, heat transfer coefficients and temperatures of the hot gas and the coolant. The resulting output of the program consists of the nodal temperature distribution, heat transfer to the coolant and heat transfer from the hot gas.

Naraghi, Mohammad H.↗
NASA NTRS1991

Direct Finite-Difference Simulations Of Turbulent Flow

Report discusses use of upwind-biased finite-difference numerical-integration scheme to simulate evolution of small disturbances and fully developed turbulence in three-dimensional flow of viscous, incompressible fluid in channel. Involves use of computational grid sufficiently fine to resolve motion of fluid at all relevant length scales.

Rai, Man Mohan↗
NASA NTRS1989

Nonlinear triggered lightning models for use in finite difference calculations

Two nonlinear triggered lightning models have been developed for use in finite difference calculations. Both are based on three species of air chemistry physics and couple nonlinearly calculated air conductivity to Maxwell's equations. The first model is suitable for use in three-dimensional modeling and has been applied to the analysis of triggered lightning on the NASA F106B Thunderstorm Research Aircraft. The model calculates number densities of positive ions, negative ions, and electrons as a function of time and space through continuity equations, including convective derivative terms. The set of equations is closed by using experimentally determined mobilities, and the mobilities are also used to determine the air conductivity. Results from the model's application to the F106B are shown. The second model is two-dimensional and incorporates an enhanced air chemistry formulation. Momentum conservation equations replace the mobility assumption of the first model. Energy conservation equations for neutrals, heavy ions, and electrons are also used. Energy transfer into molecular vibrational modes is accounted for. The purpose for the enhanced model is to include the effects of temperature into the air breakdown, a necessary step if the model is to simulate more than the very earliest stages of breakdown. Therefore, the model also incorporates a temperature-dependent electron avalanche rate. Results from the model's application to breakdown around a conducting ellipsoid placed in an electric field are shown.

Rudolph, Terence↗
NASA NTRS2012

An Analysis Technique/Automated Tool for Comparing and Tracking Analysis Modes of Different Finite Element Models

An analysis technique was developed to compare and track mode shapes for different Finite Element Models. The technique may be applied to a variety of structural dynamics analyses, including model reduction validation (comparing unreduced and reduced models), mode tracking for various parametric analyses (e.g., launch vehicle model dispersion analysis to identify sensitivities to modal gain for Guidance, Navigation, and Control), comparing models of different mesh fidelity (e.g., a coarse model for a preliminary analysis compared to a higher-fidelity model for a detailed analysis) and mode tracking for a structure with properties that change over time (e.g., a launch vehicle from liftoff through end-of-burn, with propellant being expended during the flight). Mode shapes for different models are compared and tracked using several numerical indicators, including traditional Cross-Orthogonality and Modal Assurance Criteria approaches, as well as numerical indicators obtained by comparing modal strain energy and kinetic energy distributions. This analysis technique has been used to reliably identify correlated mode shapes for complex Finite Element Models that would otherwise be difficult to compare using traditional techniques. This improved approach also utilizes an adaptive mode tracking algorithm that allows for automated tracking when working with complex models and/or comparing a large group of models.

Towner, Robert L.↗
NASA NTRS1974

Finite-difference solution for turbulent swirling compressible flow in axisymmetric ducts with struts

A finite-difference procedure for computing the turbulent, swirling, compressible flow in axisymmetric ducts is described. Arbitrary distributions of heat and mass transfer at the boundaries can be treated, and the effects of struts, inlet guide vanes, and flow straightening vanes can be calculated. The calculation procedure is programmed in FORTRAN 4 and has operated successfully on the UNIVAC 1108, IBM 360, and CDC 6600 computers. The analysis which forms the basis of the procedure, a detailed description of the computer program, and the input/output formats are presented. The results of sample calculations performed with the computer program are compared with experimental data.

Anderson, O. L.↗
NASA NTRS1974

A finite difference method for the solution of the transonic flow around harmonically oscillating wings

A finite difference method for the solution of the transonic flow about a harmonically oscillating wing is presented. The partial differential equation for the unsteady transonic flow was linearized by dividing the flow into separate steady and unsteady perturbation velocity potentials and by assuming small amplitudes of harmonic oscillation. The resulting linear differential equation is of mixed type, being elliptic or hyperbolic whereever the steady flow equation is elliptic or hyperbolic. Central differences were used for all derivatives except at supersonic points where backward differencing was used for the streamwise direction. Detailed formulas and procedures are described in sufficient detail for programming on high speed computers. To test the method, the problem of the oscillating flap on a NACA 64A006 airfoil was programmed. The numerical procedure was found to be stable and convergent even in regions of local supersonic flow with shocks.

Ehlers, E. F.↗
NASA NTRS1973

A study of three finite difference schemes and their role in asynoptic meteorological data assimilation

An investigation of three finite difference methods and their responses to the insertion of simulated satellite data is presented. A simple-one-level barotropic model is used as the forecast model, while the Mintz-Arakawa two-layer model is used to furnish the initial field, the verification fields, and the simulated satellite data. The schemes tested are the Shuman, the Matsuno-TASU, and an implicit scheme devised by McPherson. Results indicate that the schemes react to inserted data as they would react to unfiltered initial fields. Schemes which contain significant implicit viscosity are capable of damping the high frequency oscillations which occur after insertions, but such schemes may cause a loss of information. Schemes which contain less damping capability produce shock waves which damage the forecasts. It is also found that insertion of winds along with temperature data improves the forecast considerably.

Halberstan, I.↗
NASA NTRS1975

A finite-difference analysis of the nozzle starting process in an expansion tunnel

A suitable finite-difference method for computing the quasi-one-dimensional unsteady flow in an expansion tunnel nozzle was identified. The difference equations are presented along with the appropriate stability limits. A parametric study of the starting process in an expansion tunnel nozzle is made, and acceptable operating conditions were determined.

Weilmuenster, K. J.↗
NASA NTRS1977

Application of a finite difference method to the analysis of transonic flow over oscillating airfoils and wings

A finite difference method for solving the unsteady flow about harmonically oscillating wings is investigated. The procedure is based on separating the velocity potential into steady and unsteady parts and linearing the resulting unsteady differential equation for small disturbances. Solutions are obtained using relaxation procedures. The means for improving the solution stability characteristics of the relaxation process are explored. A direct procedure is formulated which permits obtaining solutions for combinations of Mach number and reduced frequency for which the relaxation process has proved unstable. The pressure distribution for an aspect ratio 5 rectangular wing oscillating in pitch is presented.

Weatherill, W. H.↗
NASA NTRS1982

Analysis of a finite difference grid

Some means of assessing the suitability of a mesh network for a finite difference calculation are investigated in this study. This has been done by a study of the nonlinear truncation errors of the scheme. It turns out that the mesh can not be properly assessed a priori. The effect of the mesh on the numerical solution depends on several factors including the mesh itself, the numerical algorithm, and the solution. Several recommendations are made with regard to generating the mesh and to assessing its suitability for a particular numerical calculation.

Klopfer, G. H.↗
NASA NTRS1983

Explicit finite difference predictor and convex corrector with applications to hyperbolic partial differential equations

An explicit finite difference scheme consisting of a predictor and a corrector has been developed and applied to solve some hyperbolic partial differential equations (PDEs). The corrector is a convex-type function which is applied at each time level and at each mesh point. It consists of a parameter which may be estimated such that for larger time steps the algorithm should remain stable and generate a fast speed of convergence to the steady-state solution. Some examples have been given.

Dey, C.↗
NASA NTRS1982

An exploratory study of finite difference grids for transonic unsteady aerodynamics

Unsteady aerodynamic forces are calculated by the XTRAN2L finite difference program which solves the complete two dimensional unsteady transonic small perturbation equation. The unsteady forces are obtained using a pulse transfer function technique which assumes the flow field behaves in a locally linear fashion about a mean condition. Forces are calculated for a linear flat plate using the default grids from the LTRAN2-NLR, LTRAN2-HI, and XTRAN3S programs. The forces are compared to the exact theoretical values for flat plate, and grid generated boundary and internal numerical reflections are observed to cause significant errors in the unsteady airloads. Grids are presented that alleviate the reflections while reducing computational time up to fifty-three percent and program size up to twenty-eight percent. Forces are presented for a six percent thick parabolic arc airfoil which demonstrate that the transform technique may be successfully applied to nonlinear transonic flows.

Seidel, D. A.↗
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 NTRS1990

Finite-difference solutions of three-dimensional rotor blade-vortex interactions

A numerical finite-difference procedure has been developed for the modeling of rotor self-generated blade-vortex interactions (BVI). The interaction velocity field is obtained through a nonlinear superposition of the rotor and the vortex wake flow fields. Vortex effects are simulated using the velocity 'transpiration' approach. Potential blade-vortex encounters are identified and tracked in time at equal increments of rotor azimuth using the lifting-line helicopter/rotor trim code CAMRAD. Accuracy of the numerical results were found to rely heavily on the CAMRAD-predicted vortex orientation, vortex strength, and user specified vortex core radius. Results for the subcritical self-generated BVI for the model OLS rotor are presented. They show good agreement with the experimental wind tunnel data.

Hassan, A. A.↗
NASA NTRS1998

Nonlinear Comparison of High-Order and Optimized Finite-Difference Schemes

The effect of reducing the formal order of accuracy of a finite-difference scheme in order to optimize its high-frequency performance is investigated using the I-D nonlinear unsteady inviscid Burgers'equation. It is found that the benefits of optimization do carry over into nonlinear applications. Both explicit and compact schemes are compared to Tam and Webb's explicit 7-point Dispersion Relation Preserving scheme as well as a Spectral-like compact scheme derived following Lele's work. Results are given for the absolute and L2 errors as a function of time.

Hixon, R.↗