Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 1,243 records · Page 69
All sourcesNASA NTRSDOE OSTI
NASA NTRS1980

Viscous compressible flow about blunt bodies using a numerically generated orthogonal coordinate system

A numerical solution to the Navier-Stokes equations was obtained for blunt axisymmetric entry bodies of arbitrary shape in supersonic flow. These equations are solved on a finite difference mesh obtained from a simple numerical technique which generates orthogonal coordinates between arbitrary boundaries. The governing equations are solved in time dependent form using Stetter's improved stability three step predictor corrector method. For the present application, the metric coefficients were obtained numerically using fourth order accurate, finite difference relations and proved to be totally reliable for the highly stretched mesh used to resolve the thin viscous boundary layer. Solutions are obtained for a range of blunt body nose shapes including concavities.

Graves, R. A., Jr.↗
NASA NTRS1978

Stress analysis study in cooled radial inflow turbine

With increased turbine inlet temperatures, numerical methods of thermal and stress analysis are becoming more valuable in the design of air-cooled turbines. This paper presents a study of the stresses associated with different cooling patterns in a radial inflow turbine rotor. The finite element method is used in the stress calculations taking into consideration centrifugal, thermal and aerodynamic loading. The effects of temperature distribution and the presence of internal cooling passages are discussed.

Hamed, A.↗
NASA NTRS1994

Radiation boundary condition and anisotropy correction for finite difference solutions of the Helmholtz equation

In this paper finite-difference solutions of the Helmholtz equation in an open domain are considered. By using a second-order central difference scheme and the Bayliss-Turkel radiation boundary condition, reasonably accurate solutions can be obtained when the number of grid points per acoustic wavelength used is large. However, when a smaller number of grid points per wavelength is used excessive reflections occur which tend to overwhelm the computed solutions. Excessive reflections are due to the incompability between the governing finite difference equation and the Bayliss-Turkel radiation boundary condition. The Bayliss-Turkel radiation boundary condition was developed from the asymptotic solution of the partial differential equation. To obtain compatibility, the radiation boundary condition should be constructed from the asymptotic solution of the finite difference equation instead. Examples are provided using the improved radiation boundary condition based on the asymptotic solution of the governing finite difference equation. The computed results are free of reflections even when only five grid points per wavelength are used. The improved radiation boundary condition has also been tested for problems with complex acoustic sources and sources embedded in a uniform mean flow. The present method of developing a radiation boundary condition is also applicable to higher order finite difference schemes. In all these cases no reflected waves could be detected. The use of finite difference approximation inevita bly introduces anisotropy into the governing field equation. The effect of anisotropy is to distort the directional distribution of the amplitude and phase of the computed solution. It can be quite large when the number of grid points per wavelength used in the computation is small. A way to correct this effect is proposed. The correction factor developed from the asymptotic solutions is source independent and, hence, can be determined once and for all. The effectiveness of the correction factor in providing improvements to the computed solution is demonstrated in this paper.

Tam, Christopher K. W.↗
NASA NTRS1990

Application of implicit numerical techniques to the solution of the three-dimensional diffusion equation

Implicit techniques for calculating three-dimensional, time-dependent heat diffusion in a cube are tested with emphasis on storage efficiency, accuracy, and speed of calculation. For this purpose, a tensor product technique with both Chebyshev collocation and finite differences and a generalized conjugate gradient technique with finite differences are used in conjunction with Crank-Nicolson discretization. An Euler explicit finite difference calculation is performed for use as a benchmark. The implicit techniques are found to be competitive with the Euler explicit method in terms of storage efficiency and speed of calculation and offer advantages both in accuracy and stability. Mesh stretching in the finite difference calculations is shown to markedly improve the accuracy of the solution.

Peltier, Leonard Joel↗
NASA NTRS1987

A mapped finite difference study of noise propagation in nonuniform ducts with mean flow

The primary objective of this work is to study noise propagation in acoustically lined variable area ducts with mean fluid flow. The method of study is numerical in nature and involves a body-fitted grid mapping procedure in conjunction with a factored-implicit finite difference technique. The mean fluid flow model used is two-dimensional, inviscid, irrotational, incompressible, and nonheat conducting. Fully-coupled solutions of the linearized gasdynamic equations are obtained for both positive and negative Mach numbers as well as for hard and soft wall conditions. The factored-implicit finite difference technique used did give rise to short wavelength perturbations, but these were dampened by the introduction of higher order artificial dissipation terms into the scheme. Results compared favorably with available numerical and experimental data.

Raad, Peter E.↗
DOE OSTI2021

Ab initio molecular dynamics on quantum computers

Ab initio molecular dynamics (AIMD) is a valuable technique for studying molecules and materials at finite temperatures where the nuclei evolve on potential energy surfaces obtained from accurate electronic structure calculations. In this work, we present an approach to running AIMD simulations on noisy intermediate-scale quantum (NISQ)-era quantum computers. The electronic energies are calculated on a quantum computer using the variational quantum eigensolver (VQE) method. Algorithms for computation of analytical gradients entirely on a quantum computer require quantum fault-tolerant hardware, which is beyond NISQ-era. Therefore, we compute the energy gradients numerically using finite differences, the Hellmann-Feynman theorem, and a correlated sampling technique. This method only requires additional classical calculations of electron integrals for each degree of freedom without any additional computations on a quantum computer beyond the initial VQE run. As a proof of concept, AIMD simulations are demonstrated for the H-2 molecule on IBM quantum devices. In addition, we demonstrate the validity of the method for larger molecules using full configuration interaction wave functions. As quantum hardware and noise mitigation techniques continue to improve, the method can be utilized for studying larger molecular systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1979

A computational fluid dynamic technique valid at the centerline for non-axisymmetric problems in cylindrical coordinates

A technique is described for the numerical solution of non-axisymmetric flow problems posed in cylindrical coordinates when the z-axis is included in the flowfield. The highlight of the technique is the manner in which the singularities at the centerline are handled. Specifically, the governing flowfield equations at r = 0 are put in a special form by applying L'Hospital's Rule. The required radial derivatives are evaluated using a one-sided, second-order accurate, first-difference. This leads to a smooth, convergent calculation of the flowfield at the centerline. This appears to be the first generally applicable numerical method for avoiding coordinate system singularities in the context of a finite-difference scheme, and could have application to many nonaxisymmetric flows. The technique is illustrated by specific results for the time-dependent flowfield inside an internal combustion engine.

Griffin, M. D.↗
NASA NTRS1973

Extravehicular mobility unit thermal simulator

The analytical methods, thermal model, and user's instructions for the SIM bay extravehicular mobility unit (EMU) routine are presented. This digital computer program was developed for detailed thermal performance predictions of the crewman performing a command module extravehicular activity during transearth coast. It accounts for conductive, convective, and radiative heat transfer as well as fluid flow and associated flow control components. The program is a derivative of the Apollo lunar surface EMU digital simulator. It has the operational flexibility to accept card or magnetic tape for both the input data and program logic. Output can be tabular and/or plotted and the mission simulation can be stopped and restarted at the discretion of the user. The program was developed for the NASA-JSC Univac 1108 computer system and several of the capabilities represent utilization of unique features of that system. Analytical methods used in the computer routine are based on finite difference approximations to differential heat and mass balance equations which account for temperature or time dependent thermo-physical properties.

Hixon, C. W.↗
NASA NTRS1976

Temperature distribution study in a cooled radial inflow turbine rotor

A numerical study to determine the temperature distribution in the rotor of a radial inflow turbine is presented. The study is based on the use of the finite element method in the three dimensional heat conduction problem. Different cooling techniques with various coolant to primary mass flow ratios are investigated. The resulting temperature distribution in the rotor are presented for comparison.

Hamed, A.↗
NASA NTRS1977

Varieties of operator manipulation

A change in MACSYMA syntax is proposed to accommodate the operator manipulators necessary to implement direct and indirect methods for the solution of differential equations, calculus of finite differences, and the fractional calculus, as well as their modern counterparts. To illustrate the benefits and convenience of this syntax extension, an example is given to show how MACSYMA's pattern-matching capability can be used to implement a particular set of operator identities which can then be used to obtain exact solutions to nonlinear differential equations.

Doohovskoy, A.↗
NASA NTRS1978

Method for calculating convective heat-transfer coefficients over turbine vane surfaces

A method for calculating laminar, transitional, and turbulent convective heat-transfer coefficients for turbine vane surfaces is described. An approximate integral solution method produced results in good agreement with a finite-difference solution. Comparisons between the two are presented. The integral solution results agreed well with the finite-difference solution results in the laminar and turbulent regions. Differences in calculating the start of transition produced a later starting point for the approximate integral solution's transitional flow regime.

Gauntner, D. J.↗
NASA NTRS1982

Applications of perturbation techniques

Two perturbation techniques were applied to two singular perturbation problems in heat transfer to obtain uniformly valid solutions which can serve as benchmarks for finite difference and finite element techniques. In the first problem, the method of strained parameters coupled with the application of a solvability condition is used to obtain a uniform solution for the problem of unsteady heat conduction in a long nearly circular cylinder. In the second problem, the method of matched asymptotic expansion coupled with Van Dyke's matching principle is used to obtain a uniform solution for the problem of one dimensional conduction-convection heat transfer of a uniform fluid flow.

Kandil, O. A.↗
NASA NTRS1984

Static and dynamic structural-sensitivity derivative calculations in the finite-element-based Engineering Analysis Language (EAL) system

The implementation of static and dynamic structural-sensitivity derivative calculations in a general purpose, finite-element computer program denoted the Engineering Analysis Language (EAL) System is described. Derivatives are calculated with respect to structural parameters, specifically, member sectional properties including thicknesses, cross-sectional areas, and moments of inertia. Derivatives are obtained for displacements, stresses, vibration frequencies and mode shapes, and buckling loads and mode shapes. Three methods for calculating derivatives are implemented (analytical, semianalytical, and finite differences), and comparisons of computer time and accuracy are made. Results are presented for four examples: a swept wing, a box beam, a stiffened cylinder with a cutout, and a space radiometer-antenna truss.

Camarda, C. J.↗
NASA NTRS1985

Spectral methods for the Euler equations. II - Chebyshev methods and shock fitting

The Chebyshev spectral collocation method for the Euler gasdynamic equations is described. It is used with shock fitting to compute several two-dimensional gasdynamic flows. Examples include a shock/acoustic wave interaction, a shock/vortex interaction, and the classical blunt-body problem. With shock fitting, the spectral method has a clear advantage over second-order finite differences in that equivalent accuracy can be obtained with far fewer grid points.

Hussaini, M. Y.↗
NASA NTRS1985

Computation of transonic vortex flows past delta wings Integral equation approach

The steady full-potential equation is written in the form of Poisson's equation, and the solution of the velocity field is expressed in terms of an integral equation. The solution consists of a surface integral of vorticity distribution on the wing and its free-vortex sheets and a volume integral of source distribution within a volume around the wing and its free-vortex sheets. The solution is obtained through successive iteration cycles. The source distribution is computed by using a mixed finite-difference scheme of the Murman-Cole type. The method is applied to delta wings. Numerical examples show that a conical shock is captured on the suction side of the wing. It is attached to the lower surface of the leading-edge vortex but does not necessarily reach to the wing surface.

Kandil, O. A.↗
NASA NTRS1985

Boundary layers on oscillating airfoils

A two-point finite-difference unsteady laminar and turbulent boundary-layer computational method has been used to investigate the properties of the flow around an airfoil (NACA 0012) oscillating through angles of attack up to 18 degrees, and for reduced frequencies of 0.01 and 0.20. The unsteady potential flow was determined using the method of Geissler. The influence of transition location on stall behavior has been investigated, using both experimentally determined transition information, and transition located at the pressure peak; the results show the need for viscous-inviscid interaction in future computations of such flows.

Carr, L. W.↗