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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

A SINDA thermal model using CAD/CAE technologies

The approach to thermal analysis described by this paper is a technique that incorporates Computer Aided Design (CAD) and Computer Aided Engineering (CAE) to develop a thermal model that has the advantages of Finite Element Methods (FEM) without abandoning the unique advantages of Finite Difference Methods (FDM) in the analysis of thermal systems. The incorporation of existing CAD geometry, the powerful use of a pre and post processor and the ability to do interdisciplinary analysis, will be described.

Rodriguez, Jose A.↗
NASA NTRS1987

Development of computational methods for unsteady aerodynamics at the NASA Langley Research Center

The current scope, recent progress, and plans for research and development of computational methods for unsteady aerodynamics at the NASA Langley Research Center are reviewed. Both integral-equations and finite-difference method for inviscid and viscous flows are discussed. Although the great bulk of the effort has focused on finite-difference solution of the transonic small-perturbation equation, the integral-equation program is given primary emphasis here because it is less well known.

Yates, E. Carson, Jr.↗
NASA NTRS1988

Accuracy analysis of the semi-analytical method for shape sensitivity calculation

The semianalytical method, widely used for calculating derivatives of static response with respect to design variables for structures modeled by finite elements, is studied in this paper. The paper shows that the method can have serious accuracy problems for shape design variables in structures modeled by beam, plate, truss, frame, and solid elements. The errors are shown to be associated with the structural model. An error index is developed to test the accuracy of the semianalytical method. It characterizes the difference in errors between a general finite-difference method and the semianalytical method. Moreover, a method for improving the accuracy of the semianalytical method (when possible) is provided. Examples are presented to demonstrate the use of the error index.

Barthelemy, Bruno↗
NASA NTRS1990

Accuracy analysis of the semi-analytical method for shape sensitivity calculation

The semianalytical method of design sensitivity analysis that is widely used for calculating derivatives of static response with respect to design variables for structures modeled by finite elements is studied. It is shown that the method can have serious accuracy problems for shape design variables in structures modeled by beam, plate, truss, frame, and solid elements. Errors are shown to be associated with an incompatibility of the sensitivity field with the structure. An error index is developed to test the accuracy of the semianalytical method. It characterizes the difference in errors between a general finite-difference method and the semianalytical method. A method for improving the accuracy of the semianalytical method (when possible) is provided. Examples are presented to demonstrate the use of the error index.

Barthelemy, Bruno↗
NASA NTRS1992

The solution of coupled Schroedinger equations using an extrapolation method

In this paper, extrapolation to the limit in a finite-difference method is applied to solve a system of coupled Schroedinger equations. This combination results in a method that only requires knowledge of the potential energy functions for the system. This numerical procedure has several distinct advantages over the more conventional methods. Namely, initial guesses for the term values are not needed; assumptions need be made about the behavior of the wavefunctions, such as the slope or magnitude in the nonclassical region; and the algorithm is easy to implement, has a firm mathematical foundation, and provides error estimates. Moreover, the method is less sensitive to round-off error than other methods since a small number of mesh points is used and it can be implemented on small computers. A comparison of the method with another numerical method shows results agreeing within 1 part in 10 exp 4.

Goorvitch, D.↗
NASA NTRS2003

High Order Filter Methods for Shock/Turbulence MHD Flows

Low-dissipative high order filter finite difference methods for shock/turbulence/combustion compressible viscous MHD flows has been constructed. Several variants of the filter approach that cater to different flow types are proposed. These filters provide a natural and efficient way for the minimization of the divergence of the magnetic field (del (raided dot) B) numerical error in the sense that no standard divergence cleaning is required. For certain 2-D MHD test problems, divergence free preservation of the magnetic fields of these filter schemes has been achieved.

Yee, H. C.↗
NASA NTRS1987

Geometrical acoustics and transonic helicopter sound

A new method is presented for predicting the impulsive noise generated by a transonic rotor blade. The method is a combined approach involving computational fluid dynamics and geometrical acoustics. A full-potential finite-difference method is used to obtain the pressure field close to the blade. A Kirchhoff integral formulation is then used to extend these finite-difference results into the far field. This Kirchhoff formula is based on geometrical acoustics approximations. It requires initial data across a plane at the sonic radius in a blade-fixed coordinate system. This data is provided by the finite-difference solution. Acoustic pressure predictions show good agreement with hover experimental data for cases with hover tip Mach numbers of 0.88 through 0.96. The cases above 0.92 tip Mach number are dominated by non-linear transonic effects seen as strong shocks on and off the blade tip. This paper gives the first successful predictions of far-field acoustic pressures for high-speed impulsive noise over a range of Mach numbers after delocalization.

Isom, Morris↗
NASA NTRS1987

Motion and decay of vortex rings submerged in a rotational flow

The interaction between vortex rings of finite strength and an axisymmetric rotational background flow is studied by a singular perturbation method, because it is difficult to use a finite-difference method to analyze the viscous decay in the small core of a vortex ring. The analysis is carried out by combining a composite solution of a vortex ring and an unsteady Euler solution for the background rotational flow. Using the method of averaging, a numerical scheme is developed to obtain an Euler solution in which the grid and time-step sizes depend solely on the length and velocity scales of the background flow. Numerical results are presented to illustrate the interaction between the trajectories and decay rates of the vortex rings and the background rotational flow.

Ishii, K.↗
NASA NTRS1980

The Split Coefficient Matrix method for hyperbolic systems of gasdynamic equations

The Split Coefficient Matrix (SCM) finite difference method for solving hyperbolic systems of equations is presented. This new method is based on the mathematical theory of characteristics. The development of the method from characteristic theory is presented. Boundary point calculation procedures consistent with the SCM method used at interior points are explained. The split coefficient matrices that define the method for steady supersonic and unsteady inviscid flows are given for several examples. The SCM method is used to compute several flow fields to demonstrate its accuracy and versatility. The similarities and differences between the SCM method and the lambda-scheme are discussed.

Chakravarthy, S. R.↗
NASA NTRS1990

A numerical method for phase-change problems

A highly accurate and efficient finite-difference method for phase-change problems with multiple moving boundaries of irregular shape is developed by employing a coordinate transformation that immobilizes moving boundaries and preserves the conservative forms of the original governing equations. The numerical method is first presented for one-dimensional phase-change problems (involving large density variation between phases, heat generation, and multiple moving boundaries) and then extended to solve two-dimensional problems (without change of densities between phases). Numerical solutions are obtained non-iteratively using an explicit treatment of the interfacial mass and energy balances and an implicit treatment of the temperature field equations. The accuracy and flexibility of the present numerical method are verified by solving some phase-change problems and comparing the results with existing analytical, semi-analytical and numerical solutions. Results indicate that one- and two-dimensional phase-change problems can be handled easily with excellent accuracies.

Kim, Charn-Jung↗
NASA NTRS2004

Adaptive Low Dissipative High Order Filter Methods for Multiscale MHD Flows

Adaptive low-dissipative high order filter finite difference methods for long time wave propagation of shock/turbulence/combustion compressible viscous MHD flows has been constructed. Several variants of the filter approach that cater to different flow types are proposed. These filters provide a natural and efficient way for the minimization of the divergence of the magnetic field [divergence of B] numerical error in the sense that no standard divergence cleaning is required. For certain 2-D MHD test problems, divergence free preservation of the magnetic fields of these filter schemes has been achieved.

Yee, H. C.↗
NASA NTRS1979

A vector-continuous loading concept for aerodynamic panel methods

An approach to the reduction of discretization errors in aerodynamic panel methods is presented. The approach is based on preventing the occurence of induced velocity singularities at panel slope discontinuities by maintaining continuity of the velocity jump vector across the panels. The approach was implemented in a two-dimensional incompressible panel method formulation and evaluated by application to several external and internal flow problems. The method is shown to exhibit a second order accuracy trend and to produce smaller errors with velocity component boundary conditions imposed on the real flow than with equipotential boundary conditions imposed on the imaginary flow behind the panels. For flows around airfoil sections with either sharp or blunt trailing edges, the method gives excellent agreement with results from a well developed finite difference method. The method is well behaved and is insensitive to irregularities in panel size distribution.

Kemp, W. B., Jr.↗
NASA NTRS1978

Finite elements and finite differences for transonic flow calculations

The paper reviews the chief finite difference and finite element techniques used for numerical solution of nonlinear mixed elliptic-hyperbolic equations governing transonic flow. The forms of the governing equations for unsteady two-dimensional transonic flow considered are the Euler equation, the full potential equation in both conservative and nonconservative form, the transonic small-disturbance equation in both conservative and nonconservative form, and the hodograph equations for the small-disturbance case and the full-potential case. Finite difference methods considered include time-dependent methods, relaxation methods, semidirect methods, and hybrid methods. Finite element methods include finite element Lax-Wendroff schemes, implicit Galerkin method, mixed variational principles, dual iterative procedures, optimal control methods and least squares.

Hafez, M. M.↗
NASA NTRS1983

A baroclinic quasigeostrophic open ocean model

A baroclinic quasigeostrophic open ocean model is presented, calibrated by a series of test problems, and demonstrated to be feasible and efficient for application to realistic mid-oceanic mesoscale eddy flow regimes. Two methods of treating the depth dependence of the flow, a finite difference method and a collocation method, are tested and intercompared. Sample Rossby wave calculations with and without advection are performed with constant stratification and two levels of nonlinearity, one weaker than and one typical of real ocean flows. Using exact analytical solutions for comparison, the accuracy and efficiency of the model is tabulated as a function of the computational parameters and stability limits set; typically, errors were controlled between 1 percent and 10 percent RMS after two wave periods. Further Rossby wave tests with realistic stratification and wave parameters chosen to mimic real ocean conditions were performed to determine computational parameters for use with real and simulated data. Finally, a prototype calculation with quasiturbulent simulated data was performed successfully, which demonstrates the practicality of the model for scientific use.

Miller, R. N.↗
NASA NTRS1981

Viscous-inviscid interaction method including wake effects for three-dimensional wing-body configurations

A viscous-inviscid interaction method has been developed by using a three-dimensional integral boundary-layer method which produces results in good agreement with a finite-difference method in a fraction of the computer time. The integral method is stable and robust and incorporates a model for computation in a small region of streamwise separation. A locally two-dimensional wake model, accounting for thickness and curvature effects, is also included in the interaction procedure. Computation time spent in converging an interacted result is, many times, only slightly greater than that required to converge an inviscid calculation. Results are shown from the interaction method, run at experimental angle of attack, Reynolds number, and Mach number, on a wing-body test case for which viscous effects are large. Agreement with experiment is good; in particular, the present wake model improves prediction of the spanwise lift distribution and lower surface cove pressure.

Streett, C. L.↗
NASA NTRS1993

Double-diffusive fingering convection in a porous medium

The characteristics of nonlinear two-dimensional horizontally periodic double-diffusive fingering convection in a saturated porous medium is investigated, using the Darcy equation including Brinkman and Forchheimer terms for the momentum equation. To solve the equations and the corresponding initial and boundary conditions, a Galerkin method is applied in the horizontal direction, and a hybrid finite difference method is used in the vertical direction. The developed computer code was used to compute the thermal convection case, and the results were found to be in good agreement with existing results.

Chen, Falin↗
NASA NTRS1987

Simulation of charged particle transport on the MPP

Computations of cosmic ray transport based upon finite difference methods are afflicted by instabilities, inaccuracies, and artifacts. To avoid these problems, a Monte Carlo formulation was developed which is closely related not only to the finite difference formulation, but also to the underlying physics of transport phenomena. Implementations of this approach are currently running on the Massively Parallel Processor at Goddard, whose enormous computing power overwhelms the poor statistical accuracy that usually limits the use of stochastic methods. In a Monte Carlo simulation of rectilinear transport, the coherent and diffusive effects that appeared are in good quantitative agreement with both finite difference and analytical calculations.

Earl, James A.↗