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.

SourceResearch connectionFrom yearThrough year

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

At least 289 records · Page 16

All sources NASA NTRS DOE OSTI

NASA NTRS1982

A computer program for the calculation of the flow field including boundary layer effects for mixed-compression inlets at angle of attack

A computer program was developed which is capable of calculating the flow field in the supersonic portion of a mixed compression aircraft inlet operating at angle of attack. The supersonic core flow is computed using a second-order three dimensional method-of-characteristics algorithm. The bow shock and the internal shock train are treated discretely using a three dimensional shock fitting procedure. The boundary layer flows are computed using a second-order implicit finite difference method. The shock wave-boundary layer interaction is computed using an integral formulation. The general structure of the computer program is discussed, and a brief description of each subroutine is given. All program input parameters are defined, and a brief discussion on interpretation of the output is provided. A number of sample cases, complete with data listings, are provided.

Vadyak, J.↗

NASA NTRS1975

Application of boundary integral equations to elastoplastic problems

The application of the boundary integral equation method (BIE) to the elastoplastic torsion problem is considered. It is found that the BIE is very suitable for the elastoplastic analysis of the torsion of prismatic bars. A comparison of the BIE with the finite difference method shows savings for the BIE concerning the number of unknowns which have to be determined and also a much faster convergence rate. Attention is given to the problem of an edge-notched beam in pure bending, taking into account a biharmonic formulation and a displacement formulation.

Mendelson, A.↗

NASA NTRS1978

CAVE3: A general transient heat transfer computer code utilizing eigenvectors and eigenvalues

The method of solution is a hybrid analytical numerical technique which utilizes eigenvalues and eigenvectors. The method is inherently stable, permitting large time steps even with the best of conductors with the finest of mesh sizes which can provide a factor of five reduction in machine time compared to conventional explicit finite difference methods when structures with small time constants are analyzed over long time periods. This code will find utility in analyzing hypersonic missile and aircraft structures which fall naturally into this class. The code is a completely general one in that problems involving any geometry, boundary conditions and materials can be analyzed. This is made possible by requiring the user to establish the thermal network conductances between nodes. Dynamic storage allocation is used to minimize core storage requirements. This report is primarily a user's manual for CAVE3 code. Input and output formats are presented and explained. Sample problems are included which illustrate the usage of the code as well as establish the validity and accuracy of the method.

Palmieri, J. V.↗

NASA NTRS1979

Finite element methodology for transient conduction/forced-convection thermal analysis

Finite element methodology for steady state thermal analysis of convectively cooled structures has been extended for transient analysis. The finite elements are based on representing the fluid passages by fluid bulk-temperature nodes and fluid-solid interface nodes. The formulation of the finite element equations for a typical flow passage is based on the weighted residual method with upwind weighting functions. Computer implementation of the convective finite element methodology using explicit and implicit time integration algorithms is described. Accuracy and efficiency of the methodology is evaluated by comparisons with analytical solutions and finite-difference lumped-parameter analyses. The comparative analyses demonstrate that finite element conduction/conduction methodology may be used to predict transient temperatures with an accuracy equal or superior to the lumped-parameter finite-difference method.

Thornton, E. A.↗

NASA NTRS1978

A study of trends and techniques for space base electronics

A sputtering system was developed to deposit aluminum and aluminum alloys by the dc sputtering technique. This system is designed for a high level of cleanliness and for monitoring the deposition parameters during film preparation. This system is now ready for studying the deposition and annealing parameters upon double-level metal preparation. A technique recently applied for semiconductor analysis, the finite element method, was studied for use in the computer modeling of two dimensional MOS transistor structures. It was concluded that the method has not been sufficiently well developed for confident use at this time. An algorithm was developed for confident use at this time. An algorithm was developed for implementing a computer study which is based upon the finite difference method. The program which was developed was modified and used to calculate redistribution data for boron and phosphorous which had been predeposited by ion implantation with range and straggle conditions. Data were generated for 111 oriented SOS films with redistribution in N2, dry O2 and steam ambients.

Trotter, J. D.↗

NASA NTRS1981

Efficient computation of the stability of three-dimensional compressible boundary layers

Methods for the computer analysis of the stability of three-dimensional compressible boundary layers are discussed and the user-oriented Compressible Stability Analysis (COSAL) computer code is described. The COSAL code uses a matrix finite-difference method for local eigenvalue solution when a good guess for the eigenvalue is available and is significantly more computationally efficient than the commonly used initial-value approach. The local eigenvalue search procedure also results in eigenfunctions and, at little extra work, group velocities. A globally convergent eigenvalue procedure is also developed which may be used when no guess for the eigenvalue is available. The global problem is formulated in such a way that no unstable spurious modes appear so that the method is suitable for use in a black-box stability code. Sample stability calculations are presented for the boundary layer profiles of an LFC swept wing.

Malik, M. R.↗

NASA NTRS1982

Nonlinear Green's function method for unsteady transonic flows

Advantages to employing Green's function in describing unsteady three-dimensional transonic flows are explored. The development of the function for application to linear subsonic and supersonic unsteady aerodynamics is reviewed. It is shown that unique solutions are possible for external flows, with all functional expressions being defined in Prandtl-Glauert space. The development of methods of using the Green's function for transonic flows is traced, noting the necessity of including the effects of significant nonlinear terms. The steady-state problem is considered to demonstrate the shock-capturing ability of the method and the usefulness of the function in the incompressible, subsonic, transonic, and supersonic areas of potential unsteady three-dimensional flows around complex configurations. Computational time is asserted to be an order of magnitude less than with finite difference methods.

Tseng, K.↗

NASA NTRS1984

Application of the Green's function method for 2- and 3-dimensional steady transonic flows

A Time-Domain Green's function method for the nonlinear time-dependent three-dimensional aerodynamic potential equation is presented. The Green's theorem is being used to transform the partial differential equation into an integro-differential-delay equation. Finite-element and finite-difference methods are employed for the spatial and time discretizations to approximate the integral equation by a system of differential-delay equations. Solution may be obtained by solving for this nonlinear simultaneous system of equations in time. This paper discusses the application of the method to the Transonic Small Disturbance Equation and numerical results for lifting and nonlifting airfoils and wings in steady flows are presented.

Tseng, K.↗

NASA NTRS2001

Computation of Nonlinear Backscattering Using a High-Order Numerical Method

The nonlinear Schrodinger equation (NLS) is the standard model for propagation of intense laser beams in Kerr media. The NLS is derived from the nonlinear Helmholtz equation (NLH) by employing the paraxial approximation and neglecting the backscattered waves. In this study we use a fourth-order finite-difference method supplemented by special two-way artificial boundary conditions (ABCs) to solve the NLH as a boundary value problem. Our numerical methodology allows for a direct comparison of the NLH and NLS models and for an accurate quantitative assessment of the backscattered signal.

Fibich, G.↗

NASA NTRS1976

A comparison of matrix methods for calculating eigenvalues in acoustically lined ducts

Three approximate methods - finite differences, weighted residuals, and finite elements - were used to solve the eigenvalue problem which arises in finding the acoustic modes and propagation constants in an absorptively lined two-dimensional duct without airflow. The matrix equations derived for each of these methods were solved for the eigenvalues corresponding to various values of wall impedance. Two matrix orders, 20 x 20 and 40 x 40, were used. The cases considered included values of wall admittance for which exact eigenvalues were known and for which several nearly equal roots were present. Ten of the lower order eigenvalues obtained from the three approximate methods were compared with solutions calculated from the exact characteristic equation in order to make an assessment of the relative accuracy and reliability of the three methods. The best results were given by the finite element method using a cubic polynomial. Excellent accuracy was consistently obtained, even for nearly equal eigenvalues, by using a 20 x 20 order matrix.

Watson, W.↗

NASA NTRS2019

Predictions of Slat Noise from the 30P30N at High Angles of Attack Using Zonal Hybrid RANS-LES

Aeroacoustic predictions of slat noise from the 30P30N three-element high-lift system at high angles of attack are presented using a zonal hybrid RANS-LES method. The simulations are part of the 5th AIAA Benchmark problems for Airframe Noise Computations (BANC-V) Workshop. An economical approach utilizing structured overset grids with spatially varying span-wise grid resolution and a high-order accurate finite difference method is described. The method is utilized for near-field predictions at three angles of attack: α = 5.5, 9.5, and 14.0 degrees. Far-field noise is obtained by propagating the near-field solution using a permeable surface Ffowcs Williams-Hawkings (FWH) method. Good agreement is obtained with both near-field and far-field Power Spectral Density (PSD) data from an experimental study of the 30P30N in the 2m x 2m Kevlar-wall wind tunnel at the Japan Aerospace Exploration Agency (JAXA). Specifically, the reduction in narrow band peaks and overall broadband noise levels with increasing angle of attack is captured well using the zonal hybrid RANS-LES method.

Housman, Jeffrey A.↗

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 NTRS2004

Comparison of Numerical Schemes for a Realistic Computational Aeroacoustics Benchmark Problem

In this work, a nonlinear structured-multiblock CAA solver, the NASA GRC BASS code, will be tested on a realistic CAA benchmark problem. The purpose of this test is to ascertain what effect the high-accuracy solution methods used in CAA have on a realistic test problem, where both the mean flow and the unsteady waves are simultaneously computed on a fully curvilinear grid from a commercial grid generator. The proposed test will compare the solutions obtained using several finite-difference methods on identical grids to determine whether high-accuracy schemes have advantages for this benchmark problem.

Hixon, R.↗

NASA NTRS2016

Thermodynamics of Gas Turbine Cycles with Analytic Derivatives in OpenMDAO

A new equilibrium thermodynamics analysis tool was built based on the CEA method using the OpenMDAO framework. The new tool provides forward and adjoint analytic derivatives for use with gradient based optimization algorithms. The new tool was validated against the original CEA code to ensure an accurate analysis and the analytic derivatives were validated against finite-difference approximations. Performance comparisons between analytic and finite difference methods showed a significant speed advantage for the analytic methods. To further test the new analysis tool, a sample optimization was performed to find the optimal air-fuel equivalence ratio, , maximizing combustion temperature for a range of different pressures. Collectively, the results demonstrate the viability of the new tool to serve as the thermodynamic backbone for future work on a full propulsion modeling tool.

Propulsion↗

NASA NTRS1984

On the use of parametric differentiation to predict weak shock waves

The method of parametric differentiation is shown to be an effective means to transform both the time-linearized and the fully nonlinear transonic small disturbance potential flow equations to a linear system of equations. The globally linearized equations are then solved in phase space by finite difference methods and integral equation methods. The predicted birth, growth, and motion of weak shock waves in two-dimensional transonic flows are compared with existing theory and experiment.

Harris, W. L.↗

NASA NTRS1977

Separating and turbulent boundary layer calculations using polynomial interpretation

Higher order numerical methods derived from polynomial spline interpolation or Hermitian differencing are applied to a separating laminar boundary layer, i.e., the Howarth problem, and the turbulent flat plate boundary layer flow. Preliminary results are presented. It is found that accuracy equal to that of conventional second order accurate finite difference methods is achieved with many fewer mesh points and with reduced computer storage and time requirements.

Rubin, S. G.↗

NASA NTRS1975

Numerical computation of three-dimensional blunt body flow fields with an impinging shock

A time-marching finite-difference method was used to solve the compressible Navier-Stokes equations for the three-dimensional wing-leading-edge shock impingement problem. The bow shock was treated as a discontinuity across which the exact shock jump conditions were applied. All interior shock layer detail such as shear layers, shock waves, jets, and the wall boundary layer were automatically captured in the solution. The impinging shock was introduced by discontinuously changing the freestream conditions across the intersection line at the bow shock. A special storage-saving procedure for sweeping through the finite-difference mesh was developed which reduces the required amount of computer storage by at least a factor of two without sacrificing the execution time. Numerical results are presented for infinite cylinder blunt body cases as well as the three-dimensional shock impingement case. The numerical results are compared with existing experimental and theoretical results.

Holst, T. L.↗

NASA NTRS1985

Simulation of rotating stall by the vortex method

The vortex method, coupled to a boundary-layer solver, is applied to the numerical simulation of high Reynolds number incompressible flow in two-dimensional cascades. Periodic conditions are imposed along the plane of the cascade, with several blades per period. Good agreement is found with two finite-difference methods for a single-blade case. When a staggered cascade is treated with five independent blades, the simulation predicts rotating stall, for a range of angles of attack and stagger, and the essential features of the flow are correct. The stall cell steadily propagates along the cascade. The sensitivity of this phenomenon to two parameters is studied, and the stall boundary is found. Quantitative results and visualizations are presented.

Spalart, P. R.↗