Search NASA⌕ Search

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.

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

At least 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS1978

Computation of three-dimensional nozzle-exhaust flow fields with the GIM code

A methodology is introduced for constructing numerical analogs of the partial differential equations of continuum mechanics. A general formulation is provided which permits classical finite element and many of the finite difference methods to be derived directly. The approach, termed the General Interpolants Method (GIM), can combined the best features of finite element and finite difference methods. A quasi-variational procedure is used to formulate the element equations, to introduce boundary conditions into the method and to provide a natural assembly sequence. A derivation is given in terms of general interpolation functions from this procedure. Example computations for transonic and supersonic flows in two and three dimensions are given to illustrate the utility of GIM. A three-dimensional nozzle-exhaust flow field is solved including interaction with the freestream and a coupled treatment of the shear layer. Potential applications of the GIM code to a variety of computational fluid dynamics problems is then discussed in terms of existing capability or by extension of the methodology.

Spradley, L. W.↗
NASA NTRS1985

Symposium on Numerical and Physical Aspects of Aerodynamic Flows, 3rd, California State University, Long Beach, CA, January 21-24, 1985, Proceedings

Among the topics discussed are hydrodynamic stability and turbulent transition, separated flow simulations employing finite difference methods, a quasi-simultaneous finite difference approach for strongly interacting flow, the Newton solution of coupled Euler and boundary layer equations, multielement airfoil aerodynamics, vorticity dynamics, the significance of the thin layer Navier-Stokes approximation, a comparison of interactive and Navier-Stokes calculations of separating boundary layer flows, and a Navier-Stokes analysis for the mechanism of shock wave/boundary layer interactions. Also discussed are transonic interactions of unsteady vortical flows, boundary layers on oscillating airfoils, a laser-interferometer skin friction meter, flow over surface-mounted semibluff bodies, computations for spinning and spin-stabilized supersonic projectiles, transition calculations in three-dimensional flows, the effect of projectile boattailing in transonic flow, flow characteristics of a body of revolution at incidence, the influence of the boundary layer on the leading edge flow of a swept wing, and an interactive scheme for three-dimensional transonic flows.

Source record↗
NASA NTRS1986

Shock-wave boundary layer interactions

Presented is a comprehensive, up-to-date review of the shock-wave boundary-layer interaction problem. A detailed physical description of the phenomena for transonic and supersonic speed regimes is given based on experimental observations, correlations, and theoretical concepts. Approaches for solving the problem are then reviewed in depth. Specifically, these include: global methods developed to predict sudden changes in boundary-layer properties; integral or finite-difference methods developed to predict the continuous evolution of a boundary-layer encountering a pressure field induced by a shock wave; coupling methods to predict entire flow fields; analytical methods such as multi-deck techniques; and finite-difference methods for solving the time-dependent Reynolds-averaged Navier-Stokes equations used to predict the development of entire flow fields. Examples are presented to illustrate the status of the various methods and some discussion is devoted to delineating their advantages and shortcomings. Reference citations for the wide variety of subject material are provided for readers interested in further study.

Delery, J.↗
NASA NTRS1999

Unified CFD Methods Via Flowfield-Dependent Variation Theory

This paper addresses the flowfield-dependent variation (FDV) methods in which complex physical phenomena are taken into account in the final form of partial differential equations to be solved so that finite difference methods (FDM) or finite element methods (FEM) themselves will not dictate the physics, but rather are no more than simply the options how to discretize between adjacent nodal points or within an element. The variation parameters introduced in the formulation are calculated from the current flowfield based on changes of Mach numbers, Reynolds numbers, Peclet numbers, and Damkohler numbers between adjacent nodal points, which play many significant roles such as adjusting the governing equations (hyperbolic, parabolic, and/or e!liptic), resolving various physical phenomena, and controlling the accuracy and stability of the numerical solution. The theory is verified by a number of example problems addressing the physical implications of the variation parameters which resemble the flowfield itself, shock capturing mechanism, transitions and interactions between inviscid/viscous, compressibility/incompressibility, and laminar/turbulent flows.

Chung, T. J.↗
NASA NTRS1973

Convergence of finite difference transient response computations for thin shells.

Numerical studies pertaining to the limits of applicability of the finite difference method in the solution of linear transient shell response problems are performed, and a computational procedure for the use of the method is recommended. It is found that the only inherent limitation of the finite difference method is its inability to reproduce accurately response discontinuities. This is not a serious limitation in view of natural constraints imposed by the extension of Saint Venant's principle to transient response problems. It is also found that the short wavelength limitations of thin shell (Bernoulli-Euler) theory create significant convergence difficulties in computed response to certain types of transverse excitations. These difficulties may be overcome, however, through proper selection of finite difference mesh dimensions and temporal smoothing of the excitation.

Sobel, L. H.↗
NASA NTRS1987

Comparison of three explicit multigrid methods for the Euler and Navier-Stokes equations

Three explicit multigrid methods, Ni's method, Jameson's finite-volume method, and a finite-difference method based on Brandt's work, are described and compared for two model problems. All three methods use an explicit multistage Runge-Kutta scheme on the fine grid, and this scheme is also described. Convergence histories for inviscid flow over a bump in a channel for the fine-grid scheme alone show that convergence rate is proportional to Courant number and that implicit residual smoothing can significantly accelerate the scheme. Ni's method was slightly slower than the implicitly-smoothed scheme alone. Brandt's and Jameson's methods are shown to be equivalent in form but differ in their node versus cell-centered implementations. They are about 8.5 times faster than Ni's method in terms of CPU time. Results for an oblique shock/boundary layer interaction problem verify the accuracy of the finite-difference code. All methods slowed considerably on the stretched viscous grid but Brandt's method was still 2.1 times faster than Ni's method.

Chima, Rodrick V.↗
NASA NTRS1987

Comparison of three explicit multigrid methods for the Euler and Navier-Stokes equations

Three explicit multigrid methods, Ni's method, Jameson's finite-volume method, and a finite-difference method based on Brandt's work, are described and compared for two model problems. All three methods use an explicit multistage Runge-Kutta scheme on the fine grid, and this scheme is also described. Convergence histories for inviscid flow over a bump in a channel for the fine-grid scheme alone show that convergence rate is proportional to Courant number and that implicit residual smoothing can significantly accelerate the scheme. Ni's method was slightly slower than the implicitly-smoothed scheme alone. Brandt's and Jameson's methods are shown to be equivalent in form but differ in their node versus cell-centered implementations. They are about 8.5 times faster than Ni's method in terms of CPU time. Results for an oblique shock/boundary layer interaction problem verify the accuracy of the finite-difference code. All methods slowed considerably on the stretched viscous grid but Brandt's method was still 2.1 times faster than Ni's method.

Chima, Rodrick V.↗
NASA NTRS1981

The stability of pseudospectral-Chebyshev methods

The pseudospectral-Chebyshev methods are shown to be convergent in variable coefficient problems and, in some cases, hyperbolic problems. The analysis demonstrates that the rate of convergence is greater for finite difference methods or the finite element method. For a single first-order hyperbolic equation, the method is seen as remaining stable even when the coefficient changes sign, although in this case it is specified that care must be taken to have adequate spatial resolution. It is noted that this fact, combined with the fact that collocation methods are easy to apply in the nonlinear case, shows that the pseudospectral method is in general preferable to the Galerkin or Tau methods.

Gottlieb, D.↗
NASA NTRS1988

A review of various flow-solid interaction analysis methods with emphasis on recent advances in turbulence models and flow analysis methods

A review on various experimental, analytical, and numerical methods for the flow-solid interaction of a nest of cylinders subjected to cross flow is presented in this paper. In order to assess the possibility of initiating a development of major numerical analysis method for flow-solid interaction of a nest of cylinders in cross flows, a review on various turbulence models, numerical methods such as the finite element method and the finite difference method, and the available numerical methods for flow-solid interactions have been included in this paper. A multiple-time-scale turbulence model for complex turbulent flows and a finite element method for high Reynolds number flows have been introduced herein for the same reason.

Kim, S.-W.↗
NASA NTRS1982

Shock-fitted Euler solutions to shock vortex interactions

The interaction of a planar shock wave with one or more vortexes is computed using a pseudospectral method and a finite difference method. The development of the spectral method is emphasized. In both methods the shock wave is fitted as a boundary of the computational domain. The results show good agreement between both computational methods. The spectral method is, however, restricted to smaller time steps and requires use of filtering techniques.

Salas, M. D.↗
NASA NTRS1982

Shock-fitted Euler solutions to shock-vortex interactions

The interaction of a planar shock wave with one or more vortexes is computed using a pseudospectral method and a finite difference method. The development of the spectral method is emphasized. In both methods the shock wave is fitted as a boundary of the computational domain. The results show good agreement between both computational methods. The spectral method is, however, restricted to smaller time steps and requires use of filtering techniques. Previously announced in STAR as N82-28061

Salas, M. D.↗
NASA NTRS1991

Eigenvalue sensitivity analysis of planar frames with variable joint and support locations

Two sensitivity equations are derived in this study based upon the continuum approach for eigenvalue sensitivity analysis of planar frame structures with variable joint and support locations. A variational form of an eigenvalue equation is first derived in which all of the quantities are expressed in the local coordinate system attached to each member. Material derivative of this variational equation is then sought to account for changes in member's length and orientation resulting form the perturbation of joint and support locations. Finally, eigenvalue sensitivity equations are formulated in either domain quantities (by the domain method) or boundary quantities (by the boundary method). It is concluded that the sensitivity equation derived by the boundary method is more efficient in computation but less accurate than that of the domain method. Nevertheless, both of them in terms of computational efficiency are superior to the conventional direct differentiation method and the finite difference method.

Chuang, Ching H.↗
NASA NTRS1992

Divergence preserving discrete surface integral methods for Maxwell's curl equations using non-orthogonal unstructured grids

Several new discrete surface integral (DSI) methods for solving Maxwell's equations in the time-domain are presented. These methods, which allow the use of general nonorthogonal mixed-polyhedral unstructured grids, are direct generalizations of the canonical staggered-grid finite difference method. These methods are conservative in that they locally preserve divergence or charge. Employing mixed polyhedral cells, (hexahedral, tetrahedral, etc.) these methods allow more accurate modeling of non-rectangular structures and objects because the traditional stair-stepped boundary approximations associated with the orthogonal grid based finite difference methods can be avoided. Numerical results demonstrating the accuracy of these new methods are presented.

Madsen, Niel K.↗
NASA NTRS1991

Computational aspects of sensitivity calculations in linear transient structural analysis

The calculation of sensitivities of displacements, velocities, accelerations, and stresses in linear, structural, and transient response problems is studied. Several existing sensitivity calculation methods and two new methods are compared for three example problems. Approximation vectors such as vibration mode shapes are used to reduce the dimensionality of the finite model. To be computationally practical in large-order problems, the overall finite difference methods must use the approximation vectors from the original design in the analyses of the perturbed models. This was found to result in poor convergence of stress sensitivities in several cases. Two semianalytical techniques are developed to overcome this poor convergence. Both new methods result in very good convergence of the stress sensitivities; the computational cost is much less than would result if the vibration modes were recalculated and then used in an overall finite difference method.

Greene, W. H.↗
NASA NTRS1989

Computational Aspects of Sensitivity Calculations in Linear Transient Structural Analysis

A study has been performed focusing on the calculation of sensitivities of displacements, velocities, accelerations, and stresses in linear, structural, transient response problems. One significant goal was to develop and evaluate sensitivity calculation techniques suitable for large-order finite element analyses. Accordingly, approximation vectors such as vibration mode shapes are used to reduce the dimensionality of the finite element model. Much of the research focused on the accuracy of both response quantities and sensitivities as a function of number of vectors used. Two types of sensitivity calculation techniques were developed and evaluated. The first type of technique is an overall finite difference method where the analysis is repeated for perturbed designs. The second type of technique is termed semianalytical because it involves direct, analytical differentiation of the equations of motion with finite difference approximation of the coefficient matrices. To be computationally practical in large-order problems, the overall finite difference methods must use the approximation vectors from the original design in the analyses of the perturbed models.

Greene, William H.↗
NASA NTRS1948

Determination of stresses in gas-turbine disks subjected to plastic flow and creep

A finite-difference method previously presented for computing elastic stresses in rotating disks is extended to include the computation of the disk stresses when plastic flow and creep are considered. A finite-difference method is employed to eliminate numerical integration and to permit nontechnical personnel to make the calculations with a minimum of engineering supervision. Illustrative examples are included to facilitate explanation of the procedure by carrying out the computations on a typical gas-turbine disk through a complete running cycle. The results of the numerical examples presented indicate that plastic flow markedly alters the elastic-stress distribution.

Millenson, M B↗