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.

SourceResearch connectionFrom yearThrough year

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

At least 1,621 records · Page 90

All sources NASA NTRS DOE OSTI

NASA NTRS1992

Incompressible viscous flow computations for the pump components and the artificial heart

A finite difference, three dimensional incompressible Navier-Stokes formulation to calculate the flow through turbopump components is utilized. The solution method is based on the pseudo compressibility approach and uses an implicit upwind differencing scheme together with the Gauss-Seidel line relaxation method. Both steady and unsteady flow calculations can be performed using the current algorithm. Here, equations are solved in steadily rotating reference frames by using the steady state formulation in order to simulate the flow through a turbopump inducer. Eddy viscosity is computed by using an algebraic mixing-length turbulence model. Numerical results are compared with experimental measurements and a good agreement is found between the two.

Kiris, Cetin↗

NASA NTRS1990

Comparisons of thermochemical nonequilibrium viscous flowfield predictions for AFE vehicle

Numerical solutions for the Aeroassist Flight Experiment vehicle were obtained from three methods at a trajectory point corresponding to the maximum aerodynamic heating. The flow regime and vehicle's speed require a viscous model of laminar flow and finite-rate chemical and thermal modeling. The computational domain covers both forebody and base such that the shock layer and near wake flowfield are included. Because of differences in computational grids, methods of solution, and models of rate equations, the results are generally in poor agreement. Temperature and species concentrations are strongly affected by the physical model equations and associated parameters. The chemistry model based on 11 species is found to yield lower translational temperature profile near stagnation than those from a seven-species model. The vibrational temperature varies according to the modeling details. All solutions indicate that strong neutral and/or molecular dissociation and weak ionization take place at the forebody and vibrational freezing is present in the afterbody expansion region where the vibrational temperature is higher than the translational temperature. Some forms of shear layer emanating from the aerobrake skirt coalesce in the region of reversed flow behind the vehicle.

Tam, Luen T.↗

NASA NTRS1988

Transonic stability and control of aircraft using CFD methods

Implementation of a capability to calculate longitudinal short-period response in the CAP-TSD (Computational Aeroelasticity Program - Transonic Small Disturbance) finite-difference code is described. The code, developed recently at the NASA Langley Research Center, is capable of solving steady and unsteady flows about complete aircraft configurations and is used primarily for aeroelastic calculations in the critical transonic speed range. The longitudinal short-period equations of motion in state-space form have been coupled to the time-accurate lift and moment calculated by the program. Transient responses to an elevator pulse for free-flying aircraft demonstrate the new capability. A trim routine is also added to the code to obtain trim automatically during steady-state flow field convergence. Stability and control derivatives are estimated from the calculated transient response by a maximum likelihood estimation program. Results for a fighter configuration and a general aviation configuration are presented to assess the capability.

Vinh, Lam-Son↗

NASA NTRS1974

A quasi-vortex-lattice method in thin wing theory

A quasi-continuous method is developed for solving thin-wing problems. For the purpose of satisfying the wing boundary conditions, the spanwise vortex distribution is assumed to be stepwise-constant, while the chordwise vortex integral is reduced to a finite sum through a modified trapezoidal rule and the theory of Chebyshev polynomials. Wing-edge and Cauchy singularities are acounted for. The total aerodynamic characteristics are obtained by an appropriate quadrature integration. The two-dimensional results for airfoils without flap deflection reproduce the exact solutions in lift and pitching moment coefficients, the leading edge suction, and the pressure difference at a finite number of points. For a flapped airfoil, the present results are more accurate than those given by the vortex-lattice method. The three-dimensional results also show an improvement over the results of the vortex-lattice method. Extension to nonplanar applications is discussed.

Lan, C. E.↗

NASA NTRS1988

Method to predict external store carriage characteristics at transonic speeds

Development of a computational method for prediction of external store carriage characteristics at transonic speeds is described. The geometric flexibility required for treatment of pylon-mounted stores is achieved by computing finite difference solutions on a five-level embedded grid arrangement. A completely automated grid generation procedure facilitates applications. Store modeling capability consists of bodies of revolution with multiple fore and aft fins. A body-conforming grid improves the accuracy of the computed store body flow field. A nonlinear relaxation scheme developed specifically for modified transonic small disturbance flow equations enhances the method's numerical stability and accuracy. As a result, treatment of lower aspect ratio, more highly swept and tapered wings is possible. A limited supersonic freestream capability is also provided. Pressure, load distribution, and force/moment correlations show good agreement with experimental data for several test cases. A detailed computer program description for the Transonic Store Carriage Loads Prediction (TSCLP) Code is included.

Rosen, Bruce S.↗

NASA NTRS1985

Hovering rotor airload prediction using a full potential flow analysis with realistic wake geometry

A three-dimensional, full potential flow analysis with realistic hover wake geometry is presented for the prediction of hovering rotor airloads. The method of analysis is based on the concept of matching inner and outer domain solutions in three dimensions. The inner domain nonlinear solution is obtained using a finite difference analysis and the outer domain solution is based on prescribed wake methodology. This formulation which includes three-dimensional wake influence, was initially validated using a fixed-wing analysis, and has been extended to hovering rotor flight. Detailed chordwise and spanwise loading results are compared with subsonic and transonic test results from two rotor configurations to illustrate the predictive capabilities of the analysis. The extension of the method to steady-level forward flight is also discussed.

Egolf, T. A.↗

DOE OSTI2025

Optimal Polynomial Smoothers and One‐Sided V‐Cycles for Poisson Problems

The solution to the Poisson equation arising from the spectral element discretization of the incompressible Navier‐Stokes equations needs robust preconditioning strategies. One such strategy is multigrid. To realize the potential of multigrid methods, effective smoothing strategies are needed. Chebyshev polynomial smoothers, in conjunction with pointwise Jacobi or additive Schwarz methods (ASMs), prove to be an effective smoother. Other polynomial smoothers, however, may provide superior convergence to the multigrid preconditioner. The authors compare the standard Chebyshev polynomial smoothers to both the novel fourth‐kind Chebyshev polynomial smoothers proposed by Lottes as well as smoothers based on the polynomial of best uniform approximation to as proposed by Kraus, Vassilevski, and Zikatanov. At the cost of symmetry, further improvements may be made. For example, a order polynomial smoother on both sides of the V‐cycle may be substituted with an order polynomial smoother on one side at no additional cost. The choice of omitting the postsmoother in favor of higher‐order polynomial presmoothing is advantageous in cases where the multigrid approximation property constant is large. The authors consider a 2D model problem based on finite differences to motivate the choice of polynomial smoother, order, and whether to apply postsmoothing for the target application of high‐order ‐geometric multigrid methods for GPU architectures. Results from both domains demonstrate the substantial improvement of these approaches over the standard Chebyshev polynomial smoother with a symmetric V‐cycle.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1978

A finite-step method for estimating the spanwise lift distribution of wings in symmetric, yawed, and rotary flight at low speeds

The finite-step method was programmed for computing the span loading and stability derivatives of trapezoidal shaped wings in symmetric, yawed, and rotary flight. Calculations were made for a series of different wing planforms and the results compared with several available methods for estimating these derivatives in the linear angle of attack range. The agreement shown was generally good except in a few cases. An attempt was made to estimate the nonlinear variation of lift with angle of attack in the high alpha range by introducing the measured airfoil section data into the finite-step method. The numerical procedure was found to be stable only at low angles of attack.

Krenkel, A. R.↗

DOE OSTI2024

A Characteristics Approach to the Finite Element Method

Herein, we present a new method for solving the linear Boltzmann transport equation. Two commonly used and well-understood methods for solving partial differential equations are the method of characteristics (MOC) and the finite element method (FEM). We propose a new method that combines the fundamental concept of the FEM with the analytic solution from the MOC to obtain coefficients for the FEM basis function expansion. Traditionally, coefficients for the FEM basis function expansion are obtained via matrix inversion. Instead, we solve for the coefficients with the MOC and represent the underlying fields with the basis function expansion using these coefficients. We provide a convergence study for our method with results from two sets of FEM basis functions: Gauss-Legendre and Gauss-Lobatto sets. We also compare two different variations of our method categorized as short characteristics and intermediate characteristics.

42 ENGINEERING↗

DOE OSTI2021

DPC Disposal Thermal Scoping Analysis

This is a progress report on thermal modeling for dual-purpose canister (DPCs) direct disposal that covers several available calculation methods and addresses creep and temperature-dependent properties in a salt repository. Three modeling approaches are demonstrated: A semi-analytical calculation method that uses linear solutions with superposition and imaging, to represent a central waste package in a larger array; A finite difference model of coupled thermal creep, implemented in FLAC2D; and An integrated finite difference thermal-hydrologic modeling approach for repositories in different generic host media, implemented in PFLOTRAN. These approaches are at different levels of maturity, and future work is expected to add refinements and establish the best applications for each.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

NASA NTRS1967

Non-stationary finite amplitude convection.

Method to describe finite amplitude convection time dependence

CONVECTIVE FLOW↗

NASA NTRS1982

Solution of the two-dimensional Euler equations with generalized coordinate transformation using flux vector splitting

An implicit finite difference code using flux vector splitting has been developed for solving the two-dimensional inviscid gas dynamics equations. The method is spatially second-order acurate, fully conservative, and uses body-conforming generalized coordinates for treating complex geometries. Numerical results have been obtained for transonic flow over a circular cylinder and airfoils. Steady results for a half cylinder (top and bottom symmetry-imposed) range from critical flow to a strong shock case with rotationally induced flow separation. Full cylinder solutions at freestream Mach number values of 0.5, however, show unsteady oscillation. A perturbation form of the method has also been developed and used to compute both fore and aft inviscid flow separation about a cylinder for a nonuniform incoming stream.

Buning, P. G.↗

NASA NTRS1983

A rotationally biased upwind difference scheme for the Euler equations

The upwind difference schemes of Godunov, Osher, Roe and van Leer are able to resolve one dimensional steady shocks for the Euler equations within one or two mesh intervals. Unfortunately, this resolution is lost in two dimensions when the shock crosses the computing grid at an oblique angle. To correct this problem, a numerical scheme was developed which automatically locates the angle at which a shock might be expected to cross the computing grid and then constructs separate finite difference formulas for the flux components normal and tangential to this direction. Numerical results which illustrate the ability of this method to resolve steady oblique shocks are presented.

Davis, S. F.↗

NASA NTRS1984

A rotationally biased upwind difference scheme for the Euler equations

Davis, S. F.↗

NASA NTRS1987

Numerical simulation of the hypersonic flow around lifting vehicles

A method for solving the viscous hypersonic flow field around realistic configurations is presented. The numerical procedure for generating the required finite difference grid and the two-factored implicit flow solver are described. Results are presented for the shuttle orbiter and a generic wing-body configuration at hypersonic Mach numbers.

Rizk, Yehia↗

NASA NTRS1987

Numerical simulation of the hypersonic flow around lifting vehicles

A method of solving the viscous hypersonic flow field around realistic configurations is presented. The numerical procedure for generating the required finite difference grid and the two-factored implicit flow solver are described. Results are presented for the shuttle orbiter and a generic wing-body configuration at hypersonic Mach numbers.

Rizk, Yehia↗

NASA NTRS1989

On the role of artificial viscosity in Navier-Stokes solvers

A method is proposed to determine directly the amount of artificial viscosity needed for stability using an eigenvalue analysis for a finite difference representation of the Navier-Stokes equations. The stability and growth of small perturbations about a steady flow over the airfoils are analyzed for various amounts of artificial viscosity. The eigenvalues were determined for a small perturbation about a steady inviscid flow over a NACA 0012 airfoil at a Mach number of 0.8 and angle of attack of 0 degrees. The movement of the eigenvalue constellation with respect to the amount of artificial viscosity is studied. The stability boundries as a function of the amount of artificial viscosity from both the eigenvalue analysis and the time marching scheme are also presented. This procedure not only allows for determining the effect of varying amounts of artificial viscosity, but also for the effects of different forms of terms for artificial viscosity.

Mahajan, Aparajit J.↗

NASA NTRS1989

A review of high-speed, convective, heat-transfer computation methods

The objective is to provide useful engineering formulations and to instill a modest degree of physical understanding of the phenomena governing convective aerodynamic heating at high flight speeds. Some physical insight is not only essential to the application of the information presented here, but also to the effective use of computer codes which may be available to the reader. Given first is a discussion of cold-wall, laminar boundary layer heating. A brief presentation of the complex boundary layer transition phenomenon follows. Next, cold-wall turbulent boundary layer heating is discussed. This topic is followed by a brief coverage of separated flow-region and shock-interaction heating. A review of heat protection methods follows, including the influence of mass addition on laminar and turbulent boundary layers. Next is a discussion of finite-difference computer codes and a comparison of some results from these codes. An extensive list of references is also provided from sources such as the various AIAA journals and NASA reports which are available in the open literature.

Tauber, Michael E.↗