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.

SourceResearch connectionFrom yearThrough year

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

At least 1,657 records · Page 92

All sources NASA NTRS DOE OSTI

NASA NTRS1976

Calculation of wind-driven surface currents in the North Atlantic Ocean

Calculations to simulate the wind driven near surface currents of the North Atlantic Ocean are described. The primitive equations were integrated on a finite difference grid with a horizontal resolution of 2.5 deg in longitude and latitude. The model ocean was homogeneous with a uniform depth of 100 m and with five levels in the vertical direction. A form of the rigid-lid approximation was applied. Generally, the computed surface current patterns agreed with observed currents. The development of a subsurface equatorial countercurrent was observed.

Rees, T. H.↗

NASA NTRS1976

Impact ejecta on the moon

The partitioning of energy and the distribution of the resultant ejecta on the moon is numerically modeled using a Eulerian finite difference grid. The impact of an iron meteoroid at 15 km/sec on a gabbroic anorthosite lunar crust is examined. The high speed impact induced flow is described over the entire hydrodynamic regime from a time where the peak pressures are 6 Mbar until the stresses everywhere in the flow are linearly elastic, and less than 5 kbar. Shock-induced polymorphic phase changes, (plagioclase and pyroxene to hollandite and perovskite), and the subsequent reversion to low pressure phases are demonstrated to enhance shock wave attenuation. A rate-dependent equation of state is used for describing the hysteretic effect of the phase change. Ballistic equations for a spherical planet are then applied to material with net velocity away from the moon.

Okeefe, J. D.↗

NASA NTRS1976

Historical evolution of vortex-lattice methods

A review of the beginning and some orientation of the vortex-lattice method were given. The historical course of this method was followed in conjunction with its field of computational fluid dynamics, spanning the period from L.F. Richardson's paper in 1910 to 1975. The following landmarks were pointed out: numerical analysis of partial differential equations, lifting-line theory, finite-difference method, 1/4-3/4 rule, block relaxation technique, application of electronic computers, and advanced panel methods.

Deyoung, J.↗

NASA NTRS1976

Higher-order numerical methods derived from three-point polynomial interpolation

Higher-order collocation procedures resulting in tridiagonal matrix systems are derived from polynomial spline interpolation and Hermitian finite-difference discretization. The equations generally apply for both uniform and variable meshes. Hybrid schemes resulting from different polynomial approximations for first and second derivatives lead to the nonuniform mesh extension of the so-called compact or Pade difference techniques. A variety of fourth-order methods are described and this concept is extended to sixth-order. Solutions with these procedures are presented for the similar and non-similar boundary layer equations with and without mass transfer, the Burgers equation, and the incompressible viscous flow in a driven cavity. Finally, the interpolation procedure is used to derive higher-order temporal integration schemes and results are shown for the diffusion equation.

Rubin, S. G.↗

NASA NTRS1976

A steady and oscillatory kernel function method for interfering surfaces in subsonic, transonic and supersonic flow

The theory, results and user instructions for an aerodynamic computer program are presented. The theory is based on linear lifting surface theory, and the method is the kernel function. The program is applicable to multiple interfering surfaces which may be coplanar or noncoplanar. Local linearization was used to treat nonuniform flow problems without shocks. For cases with imbedded shocks, the appropriate boundary conditions were added to account for the flow discontinuities. The data describing nonuniform flow fields must be input from some other source such as an experiment or a finite difference solution. The results are in the form of small linear perturbations about nonlinear flow fields. The method was applied to a wide variety of problems for which it is demonstrated to be significantly superior to the uniform flow method. Program user instructions are given for easy access.

Cunningham, A. M., Jr.↗

NASA NTRS1975

Internal flow calculations for axisymmetric supersonic inlets at angle of attack

Several solutions for the internal flow in axisymmetric supersonic inlets at angle of attack were obtained. A shock-capturing finite-difference technique was used to obtain the solutions. Good overall agreement between the shock-capturing solutions and experimental data was obtained, except in regions of strong viscous effects or boundary-layer removal. The required centerbody translation with changes in freestream Mach number and angle of attack was obtained and agrees qualitatively with experimental data.

Presley, L. L.↗

NASA NTRS1975

Numerical solutions of the unsteady Navier-Stokes equations for arbitrary bodies using boundary-fitted curvilinear coordinates

A method of automatic body-fitted curvilinear coordinate generation is described and used to construct a finite-difference solution of the full incompressible time-dependent Navier-Stokes equations for the unsteady laminar viscous flow arbitrary two-dimensional airfoils or any other two-dimensional body. A method of controlling the spacing of the coordinate lines encircling the body is developed in order to treat high Reynolds number flows, since the coordinate lines must concentrate near the surface to a greater degree as the Reynolds number increases. Multiple airfoils and submerged hydrofoils are treated as illustrative examples. The solution shows good agreement with the Blasius boundary layer solution for the flow past a semi-infinite flat plate.

Thompson, J. F.↗

NASA NTRS1975

Some examples of unsteady transonic flows over airfoils

A finite difference flutter analysis is presented for the NACA 64A-410 airfoil at M equals 0.72, where the incidence is abruptly changed from 2 to 4 degrees. The effect of gust loads is studied, and the unsteady flow adjusting process is displayed. The semi-implicit procedure of Ballhaus and Lomax (1974) is used to solve the small disturbance transonic potential equation. The physical aspects of the results, rather than the numerical details, are emphasized.

Ballhaus, W. F.↗

NASA NTRS1975

Diagonalization and simultaneous symmetrization of the gas-dynamic matrices

The hyperbolic nature of the unsteady, inviscid, gas-dynamic equations implies the existence of a similarity transformation for diagonalizing an arbitrary linear combination of coefficient matrices. It is shown that the individual matrices are simultaneously symmetrized by the similarity transformation. The transformations and their norms can be applied to the well-posedness of the Cauchy problem, linear stability theory for finite-difference approximations, and simplification of block-tridiagonal systems that arise in implicit time-split algorithms.

Warming, R. F.↗

NASA NTRS1975

A two-level trajectory decomposition algorithm featuring optimal intermediate target selection

A decomposition algorithm is presented that optimizes complex missions by partitioning the trajectory into natural segments such as ascent or entry. Each segment defines a full-rank targeting subproblem. These are solved sequentially using the Newton-Raphson algorithm. The master problem, representing the complete mission, is to determine subproblem targets and master-problem controls that optimize the mission objective subject to intersegment constraints. The gradient projection algorithm solves this problem using derivatives obtained analytically from finite-difference subproblem sensitivities. Thus, the mission is optimized by coordinating the solution of tractible subproblems. Computational results for a synchronous equatorial mission are included.

Petersen, F. M.↗

NASA NTRS1975

Multiple slot skin friction reduction

The paper investigates analytically the effect of multiple slot injection on skin friction for a representative fuselage shape (ogive-cylinder body) and evaluates the potential of slot injection as a drag reduction system in subsonic flow. Typical CTOL cruise flight conditions (Mach number equals 0.82 at altitudes of 11 km) were adopted for a fuselage 67.06 m in length and with maximum diameter of 7.32 m. The numerical method of Price and Harris (1972) was used to calculate the boundary-layer characteristics up to the first slot, while the finite-difference method of Beckwith and Bushnell (1971) was used to calculate the velocity profile downstream of one, three, five, or ten slots. An integral expression is proposed for characterizing skin friction reduction effectiveness, and it is seen that large reductions in viscous drag (50%) are available through slot injection. Skin friction reduction is improved by increasing the number of injection slots but at a diminishing rate.

Howard, F. G.↗

NASA NTRS1976

On the computation of the transonic perturbation flow field around two- and three-dimensional oscillating 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 linearizing the resulting unsteady differential equation for small disturbances. Solutions are obtained using relaxation procedures. It is determined that there is a limit on reduced frequency, which is a function of Mach number and size of mesh region, above which the relaxation procedures will not converge. It is found that row line relaxation is more efficient than column relaxation and results are presented for a rectangular wing in harmonic pitch.

Weatherill, W. H.↗

NASA NTRS1976

Direct measurements at Mach 6 of turbulent skin-friction reduction by injection from single and multiple flush slots

The surface skin-friction reduction downstream of one to four successive flush slots injecting at an angle of 10 deg into a turbulent Mach 6 boundary layer has been determined from direct measurements of surface shear. Increasing the number of succeeding slots increases the skin-friction reduction, but this trend reverses at high injection rates. The incremental improvement in skin-friction reduction decreases with increasing number of slots. Comparison with previously reported step slot data indicates that step slots with tangential injection are more effective in reducing skin friction than the present flush slot configuration. Finite-difference predictions are in reasonable agreement with skin-friction and boundary-layer profile data.

Srokowski, A. J.↗

NASA NTRS1976

Advances in the application of fast semidirect computational methods in transonic flow

The paper uses finite-difference algorithms called 'fast direct elliptic solvers' within an iteration scheme for the rapid solution of the equations of inviscid transonic aerodynamics. The methods are called 'direct' (or semidirect) because the entire computational field is solved at once rather than in successive traverses over the field. These semidirect iterative methods have been limited here to the investigation of two-dimensional steady inviscid flow over airfoils in subsonic free stream.

Martin, E. D.↗

NASA NTRS1976

Comparison of a two-dimensional shock impingement computation with experiment

Results of computations of two-dimensional viscous blunt-body flowfields with an impinging shock wave, with a time-dependent finite-difference method employed to solve the complete set of Navier-Stokes equations, are compared with experimental results. The experimental results were obtained in a 20-inch hypersonic tunnel with a planar shock impinging on the cylindrical leading edge of a fin, hence with the shock parallel to the centerline of the leading edge, so that type III and type IV interference patterns were generated. Close agreement is found. The overall effects of smoothing and grid size on the calculations are determined. A 31 x 51 mesh is adequate for wall pressure values (except in peaked regions).

Tannehill, J. C.↗

NASA NTRS1976

Computation of the inviscid supersonic flow over an external axial corner

A second-order finite-difference procedure is used to evaluate the inviscid supersonic flowfield surrounding an external axial corner composed of swept planar compression surfaces and representing the inlets on existing high-speed aircraft. The governing partial differential equations in conservation-law form are hyperbolic with respect to the axial coordinate and are solved iteratively by means of MacCormack's algorithm. The procedure treats both the peripheral shock wave and vortical singularities as discontinuities. Numerical results are presented for two parametric studies regarding the effects on the flowfield of varying the free-stream Mach number and the leading edge sweep of the horizontal wedge. Results of parametric Mach number study agree with the Mach number independence principle in that as the Mach number increases, such characteristics as shock shape, cross-flow sonic line location, and vortical singularity position approach an asymptote.

Kutler, P.↗

NASA NTRS1976

Diffraction of a shock wave by a compression corner. I - Regular reflection

The unsteady, two-dimensional flowfield resulting from the interaction of a moving planar shock wave with a compression corner is determined using a second-order, discontinuity-fitting, finite-difference approach. The time-dependent Euler equations are transformed to normalize the distance between the body and peripheral shock and to include the existing self-similar property of the flow. The resulting set of partial differential equations in conservation-law form is then solved in a time-dependent fashion using MacCormack's scheme. The vortical singularity, which lies on the body surface, and the single reflected shock are both treated as discontinuities in the numerical procedure. The results of the numerical simulation compare quite favorably with existing experimental interferograms and yield better flowfield resolution than previous first-order, shock-capturing, numerical solutions.

Kutler, P.↗

NASA NTRS1976

A numerical study of viscous flow around an airfoil

An integrodifferential method, previously formulated in terms of velocity and vorticity vectors, is reformulated in terms of stream function and vorticity for two-dimensional incompressible viscous flows. The reformulated integrodifferential method is shown to retain the distinguishing feature of the previous formulation in permitting the confinement of the solution field to the viscous region of the flow and consequently offers great computational advantages. The application of this procedure in a study of an incompressible flow around an impulsively started 9% thick symmetric Joukowski airfoil at an angle of attack of 15 deg and a Reynolds number of 1000 is discussed. Numerical results are presented and compared with available finite-difference results.

Wu, J. C.↗