Search NASA⌕ Search

SEARCH · Search NASA

Results for “Finite difference scheme”

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 469 records · Page 26

A procedure for calculating laminar, transitional and/or turbulent three-dimensional boundary layers including surface curvature effects about arbitrary blunt-nosed bodies

The governing equations for laminar, transitional and/or turbulent, incompressible, three-dimensional boundary layers are solved numerically. The equations are developed in terms of an orthogonal surface coordinate system and include surface curvature effects. The coordinates are generated numerically, while the inviscid flow is obtained from a general potential flow procedure. The only restrictions on body geometry are that it possess a blunt nose and a plane of symmetry. The boundary-layer equations, after being transformed into similarity type variables, are solved using the implicit finite-difference Krause scheme. Various test cases are presented to establish the accuracy of the resulting computer program.

Dwoyer, D. L.↗

Simulated transonic flows for aircraft with nacelles, pylons, and winglets

A computational method which simulates transonic flow about wing-fuselage configurations has been extended to include the treatment of multiple body and non-planar wing surfaces. The finite difference relaxation scheme is characterized by a modified small disturbance flow equation and multiple embedded grid system. Wing-body combinations with as many as four nacelles/pods, four pylons, and wing-tip-mounted winglets can be analyzed. A scheme for modeling inlet spillage and engine exhaust interference effects has been included. Computed results are correlated with experimental data for three transport configurations.

Boppe, C. W.↗

Finite element methods for nonlinear elastostatic problems in rubber elasticity

A number of finite element methods for the analysis of nonlinear problems in rubber elasticity are outlined. Several different finite element schemes are discussed. These include the augmented Lagrangian method, continuation or incremental loading methods, and associated Riks-type methods which have the capability of incorporating limit point behavior and bifurcations. Algorithms for the analysis of limit point behavior and bifurcations are described and the results of several numerical experiments are presented. In addition, a brief survey of some recent work on modelling contact and friction in elasticity problems is given. These results pertain to the use of new nonlocal and nonlinear friction laws.

Oden, J. T.↗

Application of a viscous-inviscid interaction procedure to predict separated flows with heat transfer

A viscous-inviscid interaction procedure is described for predicting heat transfer in separated flows. The separating flow in a rearward-facing step/asymmetric channel expansion is considered. For viscous regions, the boundary-layer momentum and continuity equations are solved inversely in a coupled manner by a finite-difference numerical scheme. The streamwise convective term is altered to permit marching the solution through regions of reversed flow. The inviscid flow is computed by numerically solving the Laplace equation for stream-function in the region bounded by the displacement surfaces used in the inverse boundary-layer solution. The viscous and inviscid solutions are repeated iteratively until the edge velocities obtained from both solutions are in agreement. Predictions using this method compare favorably with experimental data and other predictions.

Hall, E. J.↗

A three-dimensional boundary-layer analysis including heat-transfer and blade-rotation effects

An analysis has been developed for calculating three-dimensional viscous flows over arbitrary bodies through the solution of complete three-dimensional boundary-layer equations. The governing equations include the effect of blade rotation and are written in nonorthogonal surface-fitted coordinates to allow maximum flexibility. A coupled finite-difference numerical scheme is used for solving the governing equations. The accuracy of this scheme has been assessed through comparisons with known analytical/experimental results for several test cases. The present analysis has then been applied to compute viscous flows over a turbine endwall and airfoil suction surface, and fairly good correlation has been obtained with available experimental data for surface heat transfer and viscous streamline skewing.

Vatsa, V. N.↗

Limitations of the boundary-layer equations for predicting laminar symmetric sudden expansion flows

A finite-difference solution scheme is used to study the limitations and capabilities of the boundary-layer equation model for flow through abrupt, symmetric expansions. Solutions of the boundary-layer equations are compared with previous numerical predictions and experimental measurements. Some flow parameters are not well predicted for Reynolds numbers below 200. Global iteration over the flow field to include upstream effects does not significantly influence the predictions. Axisymmetric and two-dimensional flows are investigated. The effect of initial conditions is discussed.

Lewis, J. P.↗

Aerodynamic analysis for aircraft with nacelles, pylons, and winglets at transonic speeds

A computational method has been developed to provide an analysis for complex realistic aircraft configurations at transonic speeds. Wing-fuselage configurations with various combinations of pods, pylons, nacelles, and winglets can be analyzed along with simpler shapes such as airfoils, isolated wings, and isolated bodies. The flexibility required for the treatment of such diverse geometries is obtained by using a multiple nested grid approach in the finite-difference relaxation scheme. Aircraft components (and their grid systems) can be added or removed as required. As a result, the computational method can be used in the same manner as a wind tunnel to study high-speed aerodynamic interference effects. The multiple grid approach also provides high boundary point density/cost ratio. High resolution pressure distributions can be obtained. Computed results are correlated with wind tunnel and flight data using four different transport configurations. Experimental/computational component interference effects are included for cases where data are available. The computer code used for these comparisons is described in the appendices.

Boppe, Charles W.↗

External store carriage loads prediction at transonic speeds

A computational method for prediction of external store carriage loads at transonic speeds is described. The geometric flexibility required for treatment of isolated and underwing, pylon-mounted stores is achieved by computing 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 finite difference relaxation scheme, developed specifically for modified transonic small disturbance flow equations, enhances numerical stability and accuracy. As a result, more accurate treatment of low aspect ratio, highly swept and tapered wing planforms is possible. A limited supersonic freestream capability is also provided. Pressure, load distribution, force and moment correlations show good agreement for several test cases.

Rosen, Bruce S.↗

Method to predict external store carriage loads at transonic speeds

A computational method for prediction of external store carriage loads at transonic speeds is described. The geometric flexibility required for treatment of isolated and underwing, pylon mounted stores is achieved by computing 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 finite difference relaxation scheme, developed specifically for modified transonic small disturbance flow equations, enhances numerical stability and accuracy. As a result, more accurate treatment of low aspect ratio, highly swept and tapered wing planforms is possible. A limited supersonic freestream capability is also provided. Pressure, load distribution, force and moment correlation show good agreement for several test cases.

Rosen, Bruce S.↗

A finite-difference, frequency-domain numerical scheme for the solution of the linearized unsteady Euler equations

A numerical method is developed for solving periodic, three-dimensional, vortical flows around lifting airfoils in subsonic flow. The first-order method, that is presented, fully accounts for the distortion effects of the nonuniform mean flow on the convected upstream vortical disturbances. The unsteady velocity is split into a vortical component which is a known function of the upstream flow conditions and the Lagrangian coordinates of the mean flow, and an irrotational field whose potential satisfies a nonconstant-coefficient, inhomogeneous, convective wave equation. Using an elliptic coordinate transformation, the unsteady boundary value problem is solved in the frequency domain on grids which are determined as a function of the Mach number and reduced frequency. Extensive comparisons are made with known solutions to unsteady vortical flow problems, and it is seen that the agreement is generally very good for reduced frequencies ranging from 0 up to 4.

Scott, James R.↗

Navier-Stokes Flowfield Simulation of Boeing 747-200 as Platform for SOFIA

Steady and unsteady viscous, three-dimensional flowfields are calculated using a thin layer approximation of Navier-Stokes equations in conjunction with Chimera overset grids. The finite-difference numerical scheme uses structured grids and a pentadiagonal flow solver called "OVERFLOW". The configuration of Boeing 747-200 has been chosen as one of configurations to be used as a platform for the SOFIA (Stratospheric Observatory For Infrared Astronomy). Initially, the steady flowfield of the full aircraft is calculated for the clean configuration (without a cavity to house telescope). This solution is then used to start the unsteady flowfield of a configuration containing cavity housing the observation telescope and its peripheral units. Analysis of unsteady flowfield in the cavity and its influence on the tail empennage, as well as the noise due to turbulence and optical quality of the flow are the main focus of this study. For the configuration considered here, the telescope housing cavity is located slightly downstream of the portwing. The entire flow-field is carefully constructed using 45 overset grids and consists of nearly 4 million grid points. All the computations axe done at one freestream flow condition of M(sub infinity) = 0.85, alpha = 2.5deg, and a Reynolds of Re = 1.85x10deg

Srinivasan, G.R.↗

A Robust Absorbing Boundary Condition for Compressible Flows

An absorbing non-reflecting boundary condition (NRBC) for practical computations in fluid dynamics and aeroacoustics is presented with theoretical proof. This paper is a continuation and improvement of a previous paper by the author. The absorbing NRBC technique is based on a first principle of non reflecting, which contains the essential physics that a plane wave solution of the Euler equations remains intact across the boundary. The technique is theoretically shown to work for a large class of finite volume approaches. When combined with the hyperbolic conservation laws, the NRBC is simple, robust and truly multi-dimensional; no additional implementation is needed except the prescribed physical boundary conditions. Several numerical examples in multi-dimensional spaces using two different finite volume schemes are illustrated to demonstrate its robustness in practical computations. Limitations and remedies of the technique are also discussed.

Loh, Ching Y.↗

Combined Bluntness and Roughness Effects on Cones at Hypersonic Speeds

This computational study investigates the effects of discrete roughness elements on a blunt cone at zero degrees angle of attack in a Mach 6 flow. Motivation was provided by experiments conducted in the Air Force Research Laboratory Mach 6 High Reynolds Number facility on a 7-degree half-angle cone with a roughness array located at 45 degrees from the apex on two nose tips of different blutness but equivalent roughness Reynolds number. Transition was only affected on the blunter cone, indicating that the transition onset is associated with the combined effects of bluntness and roughness. The present study investigates the 15.24 mm nose radius, 420 azimuthal wavenumber case via Navier-Stokes computations of the laminar base flow and instability analysis. Plane-marching parabolized stability equations (PSE) and inflow-resolvent analysis based on the three-dimensional, harmonic linearized Navier-Stokes equations (HLNSE) are used to calculate the amplification of disturbances along the roughness wake as well as over the roughness nearfield. Results show that the roughness shape can have a great impact on the characteristics of the most amplified wake instabilities. For the experimental configuration with cubic roughness elements of 15 μ m height, the flow is marginally unstable. For prismatic elements of 20 μ m height, the PSE predicts a logarithmic disturbance amplification ratio of N = 5.6 along its wake, but this ratio increases to N = 9.6 when the amplification over the roughness and separation regions is included in the inflow-resolvent analysis.

Boundary-layer transition↗

Finite-Difference Algorithms For Computing Sound Waves

Governing equations considered as matrix system. Method variant of method described in "Scheme for Finite-Difference Computations of Waves" (ARC-12970). Present method begins with matrix-vector formulation of fundamental equations, involving first-order partial derivatives of primitive variables with respect to space and time. Particular matrix formulation places time and spatial coordinates on equal footing, so governing equations considered as matrix system and treated as unit. Spatial and temporal discretizations not treated separately as in other finite-difference methods, instead treated together by linking spatial-grid interval and time step via common scale factor related to speed of sound.

Davis, Sanford↗

The incremental unknowns-a multilevel scheme for the simulation of turbulent channel flows

In numerical simulation of complex flows, it is important to identify different length scales of the flow and treat them differently. In this report, we introduce a new multilevel scheme for simulating turbulent channel flows. Two different versions of the scheme, namely the spectral and finite difference versions, are presented. The spectral version of the scheme is based on a spectral-Galerkin formulation which provides a natural decomposition of the flow into small and large wavelength parts, and which leads to linear systems that can be solved with quasi-optimal computational complexity. In the finite difference version, the Incremental Unknown (IU) is used to separate the length scales. Preliminary numerical results indicate that the scheme is well suited for turbulence computations and provides results which are comparable to that by Direct Numerical Simulation (DNS) but with significantly less CPU time.

Chen, M.↗

Magnetic Field Solver

The Magnetic Field Solver computer program calculates the magnetic field generated by a group of collinear, cylindrical axisymmetric electromagnet coils. Given the current flowing in, and the number of turns, axial position, and axial and radial dimensions of each coil, the program calculates matrix coefficients for a finite-difference system of equations that approximates a two-dimensional partial differential equation for the magnetic potential contributed by the coil. The program iteratively solves these finite-difference equations by use of the modified incomplete Cholesky preconditioned-conjugate-gradient method. The total magnetic potential as a function of axial (z) and radial (r) position is then calculated as a sum of the magnetic potentials of the individual coils, using a high-accuracy interpolation scheme. Then the r and z components of the magnetic field as functions of r and z are calculated from the total magnetic potential by use of a high-accuracy finite-difference scheme. Notably, for the finite-difference calculations, the program generates nonuniform two-dimensional computational meshes from nonuniform one-dimensional meshes. Each mesh is generated in such a way as to minimize the numerical error for a benchmark one-dimensional magnetostatic problem.

Ilin, Andrew V.↗

On the dynamics of approximating schemes for dissipative nonlinear equations

Since one can rarely write down the analytical solutions to nonlinear dissipative partial differential equations (PDE's), it is important to understand whether, and in what sense, the behavior of approximating schemes to these equations reflects the true dynamics of the original equations. Further, because standard error estimates between approximations of the true solutions coming from spectral methods - finite difference or finite element schemes, for example - and the exact solutions grow exponentially in time, this analysis provides little value in understanding the infinite time behavior of a given approximating scheme. The notion of the global attractor has been useful in quantifying the infinite time behavior of dissipative PDEs, such as the Navier-Stokes equations. Loosely speaking, the global attractor is all that remains of a sufficiently large bounded set in phase space mapped infinitely forward in time under the evolution of the PDE. Though the attractor has been shown to have some nice properties - it is compact, connected, and finite dimensional, for example - it is in general quite complicated. Nevertheless, the global attractor gives a way to understand how the infinite time behavior of approximating schemes such as the ones coming from a finite difference, finite element, or spectral method relates to that of the original PDE. Indeed, one can often show that such approximations also have a global attractor. We therefore only need to understand how the structure of the attractor for the PDE behaves under approximation. This is by no means a trivial task. Several interesting results have been obtained in this direction. However, we will not go into the details. We mention here that approximations generally lose information about the system no matter how accurate they are. There are examples that show certain parts of the attractor may be lost by arbitrary small perturbations of the original equations.

Jones, Donald A.↗

Design of a Variational Multiscale Method for Turbulent Compressible Flows

A spectral-element framework is presented for the simulation of subsonic compressible high-Reynolds-number flows. The focus of the work is maximizing the efficiency of the computational schemes to enable unsteady simulations with a large number of spatial and temporal degrees of freedom. A collocation scheme is combined with optimized computational kernels to provide a residual evaluation with computational cost independent of order of accuracy up to 16th order. The optimized residual routines are used to develop a low-memory implicit scheme based on a matrix-free Newton-Krylov method. A preconditioner based on the finite-difference diagonalized ADI scheme is developed which maintains the low memory of the matrix-free implicit solver, while providing improved convergence properties. Emphasis on low memory usage throughout the solver development is leveraged to implement a coupled space-time DG solver which may offer further efficiency gains through adaptivity in both space and time.

Design↗