Search NASA⌕ Search

SEARCH · Search NASA

Results for “Numerical 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 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Numerical simulation of transition in boundary layers

Different implications of employment of the temporal or spatial model in direct numerical simulations of transition based on the complete Navier-Stokes equations are outlined, and it is noted that in situations where approximations are not justified, the temporal model cannot be employed. The discussion is then concentrated on developing numerical methods for solving the Navier-Stokes equations based on the spatial model that are applicable for realistic situations for which the temporal model would not be valid. It is concluded that in the future numerical simulations based on the spatial approach will replace, to a large extent, the temporal simulations.

Fasel, H.↗
NASA NTRS1992

Aerodynamics of loaded cascades in subsonic flows subject to unsteady three-dimensional vortical disturbances

A highly efficient numerical method is developed for three-dimensional, periodic, vortical flows around a cascade of loaded airfoils. The method linearizes Euler's equations about the mean flow of the cascade and thus fully accounts for the effects of distortion of the vortical disturbances as they propogate and interact with the cascade mean flow. The numerical scheme is based on splitting the unsteady velocity into vortical and potential parts. The latter is governed by a non-constant coefficient inhomogeneous convective wave equation. A new and computationally suitable out-flow conditions are derived and avoid the difficulties associated with the singular velocity downstream. Solutions were obtained in the frequency domain by using a body-fitted coordinate system. Results are presented to demonstrate the effects of the out-flow boundary conditions, cascade spacing, mean blade loading and gust upstream conditions on the aerodynamics response and unsteady pressure field of a cascade.

Fang, J.↗
NASA NTRS1993

Non-oscillatory and non-diffusive solution of convection problems by the iteratively reweighted least-squares finite element method

A comparative description is presented for the least-squares FEM (LSFEM) for 2D steady-state pure convection problems. In addition to exhibiting better control of the streamline derivative than the streamline upwinding Petrov-Galerkin method, numerical convergence rates are obtained which show the LSFEM to be virtually optimal. The LSFEM is used as a framework for an iteratively reweighted LSFEM yielding nonoscillatory and nondiffusive solutions for problems with contact discontinuities; this method is shown to convect contact discontinuities without error when using triangular and bilinear elements.

Jiang, Bo-Nan↗
NASA NTRS1999

A Parallel, Finite-Volume Algorithm for Large-Eddy Simulation of Turbulent Flows

A parallel, finite-volume algorithm has been developed for large-eddy simulation (LES) of compressible turbulent flows. This algorithm includes piecewise linear least-square reconstruction, trilinear finite-element interpolation, Roe flux-difference splitting, and second-order MacCormack time marching. Parallel implementation is done using the message-passing programming model. In this paper, the numerical algorithm is described. To validate the numerical method for turbulence simulation, LES of fully developed turbulent flow in a square duct is performed for a Reynolds number of 320 based on the average friction velocity and the hydraulic diameter of the duct. Direct numerical simulation (DNS) results are available for this test case, and the accuracy of this algorithm for turbulence simulations can be ascertained by comparing the LES solutions with the DNS results. The effects of grid resolution, upwind numerical dissipation, and subgrid-scale dissipation on the accuracy of the LES are examined. Comparison with DNS results shows that the standard Roe flux-difference splitting dissipation adversely affects the accuracy of the turbulence simulation. For accurate turbulence simulations, only 3-5 percent of the standard Roe flux-difference splitting dissipation is needed.

Bui, Trong T.↗
NASA NTRS1983

Symposium on Numerical and Physical Aspects of Aerodynamic Flows, 2nd, California State University, Long Beach, CA, January 17-20, 1983, Proceedings

The present conference covers topics concerning the measurement and calculation of interactive flows, together with problems posed by subsonic and transonic wings, missiles, and ships. Discussions are presented on the time-dependent finite difference simulation of unsteady interactive flows, Navier-Stokes equation methods, numerical solutions for spatially periodic boundary layers, the application of unsteady laminar tripple deck theory to viscous-inviscid interaction, the coupling of boundary layer and Euler equation solutions, and viscous-inviscid flow interactions. Also discussed are leading and trailing edge flows, three-dimensional wing flows, small disturbance calculations including entropy corrections, an inviscid computational method for tactical missiles, and boundary layer and flow separation characteristics of bodies of revolution at incidence.

Source record↗
NASA NTRS1973

Cracked finite elements proposed for NASTRAN

The recent introduction of special crack-tip singularity elements, usually referred to as cracked elements, has brought the power and flexibility of the finite-element method to bear much more effectively on fracture mechanics problems. This paper recalls the development of two cracked elements and presents the results of some applications proving their accuracy and economy. Judging from the available literature on numerical methods in fracture mechanics, it seems clear that the elements described have been used more extensively than any others in practical fracture mechanics applications.

Aberson, J. A.↗
NASA NTRS1977

Theory of viscous transonic flow over airfoils at high Reynolds number

This paper considers viscous flows with unseparated turbulent boundary layers over two-dimensional airfoils at transonic speeds. Conventional theoretical methods are based on boundary layer formulations which do not account for the effect of the curved wake and static pressure variations across the boundary layer in the trailing edge region. In this investigation an extended viscous theory is developed that accounts for both effects. The theory is based on a rational analysis of the strong turbulent interaction at airfoil trailing edges. The method of matched asymptotic expansions is employed to develop formal series solutions of the full Reynolds equations in the limit of Reynolds numbers tending to infinity. Procedures are developed for combining the local trailing edge solution with numerical methods for solving the full potential flow and boundary layer equations. Theoretical results indicate that conventional boundary layer methods account for only about 50% of the viscous effect on lift, the remaining contribution arising from wake curvature and normal pressure gradient effects.

Melnik, R. E.↗
NASA NTRS1978

Inviscid, nonadiabatic flow fields over blunt, sonic corner bodies for outer planet entry conditions by a method of integral relations

An investigation has been made into the ability of a method of integral relations to calculate inviscid zero degree angle of attack, radiative heating distributions over blunt, sonic corner bodies for some representative outer planet entry conditions is investigated. Comparisons have been made with a more detailed numerical method, a time asymptotic technique, using the same equilibrium chemistry and radiation transport subroutines. An effort to produce a second order approximation (two-strip) method of integral relations code to aid in this investigation is also described and a modified two-strip routine is presented. Results indicate that the one-strip method of integral relations cannot be used to obtain accurate estimates of the radiative heating distribution because of its inability to resolve thermal gradients near the wall. The two-strip method can sometimes be used to improve these estimates; however, the two-strip method has only a small range of conditions over which it will yield significant improvement over the one-strip method.

Gnoffo, P. A.↗
NASA NTRS1980

Geometric methods in computational fluid dynamics

General methods for the construction of geometric computational fluid dynamic algorithms are presented which simulate a variety of flow fields in various nontrivial regions. Included are: basic developments with tensors; various forms for the equations of motion; generalized numerical methods and boundary conditions; and methods for mesh generation to meet the strong geometric constraints of turbomachines. Coordinate generation is shown generally to yield mesh descriptions from one or more transformations that are smoothly joined together to form a composite mesh.

Eiseman, P. R.↗
NASA NTRS1981

Scale effects on turbulent boundary layer development and flow separation around V/STOL inlets at high incidence

Numerical methods for calculating laminar and turbulent boundary layers development around V/STOL engine inlets at high incidence angles, along with the procedures for predicting flow separation, are presented. Results of scale-effects, which are obtained by a numerical scaling procedure on the boundary layer characteristics and incidence angles at onset of separation are discussed. The interesting 'cross-over' phenomena, where the full-scale models actually exhibit earlier separation than the scaled-models, is illustrated for a typical V/STOL inlet at a certain operating condition. Some of the numerical results are compared with the existing wind-tunnel test data for a 1/6 scale inlet model to demonstrate the validity of the numerical approach.

Chou, D. C.↗
NASA NTRS1996

Mixing and Transition Control Studied

Considerable progress in understanding nonlinear phenomena in both unbounded and wallbounded shear flow transition has been made through the use of a combination of high- Reynolds-number asymptotic and numerical methods. The objective of this continuing work is to fully understand the nonlinear dynamics so that ultimately (1) an effective means of mixing and transition control can be developed and (2) the source terms in the aeroacoustic noise problem can be modeled more accurately. Two important aspects of the work are that (1) the disturbances evolve from strictly linear instability waves on weakly nonparallel mean flows so that the proper upstream conditions are applied in the nonlinear or wave-interaction streamwise region and (2) the asymptotic formulations lead to parabolic problems so that the question of proper out-flow boundary conditions--still a research issue for direct numerical simulations of convectively unstable shear flows--does not arise. Composite expansion techniques are used to obtain solutions that account for both mean-flow-evolution and nonlinear effects. A previously derived theory for the amplitude evolution of a two-dimensional instability wave in an incompressible mixing layer (which is in quantitative agreement with available experimental data for the first nonlinear saturation stage for a plane-jet shear layer, a circular-jet shear layer, and a mixing layer behind a splitter plate) have been extended to include a wave-interaction stage with a three-dimensional subharmonic. The ultimate wave interaction effects can either give rise to explosive growth or an equilibrium solution, both of which are intimately associated with the nonlinear self-interaction of the three dimensional component. The extended theory is being evaluated numerically. In contrast to the mixing-layer situation, earlier comparisons of theoretical predictions based on asymptotic methods and experiments in wall-bounded shear-flow transition have been somewhat lacking in one aspect or another. The current work strongly suggests that the main weakness is the underlying asymptotic representation of the linear "part" of the problem and not the explicit modeling of the nonlinear/wave-interaction effects. Consequently, the long-wave-length/high-Reynolds-number asymptotic limit for the Blasius boundary-layer stability problem was reexamined, and a new dispersion relationship for the instability waves that is uniformly valid for both the upper- and lower branch regions to the required order of approximation was obtained. A comparison with numerical results, obtained by solving the Orr-Sommerfeld stability problem, shows that the asymptotic formula provides surprisingly good results, even for values of the frequency parameter usually encountered in experimental investigations. This is particularly evident in the dynamically important upper-branch region, where much of the nonlinear interactions in transition experiments are believed to take place. The result is important in that it can be used to greatly improve the accuracy of weakly nonlinear critical-layer-based theories, and a consistent nonlinear theory is currently under evaluation.

Source record↗
NASA NTRS1972

Analysis of supersonic combustion flow fields with embedded subsonic regions

The viscous characteristic analysis for supersonic chemically reacting flows was extended to include provisions for analyzing embedded subsonic regions. The numerical method developed to analyze this mixed subsonic-supersonic flow fields is described. The boundary conditions are discussed related to the supersonic-subsonic and subsonic-supersonic transition, as well as a heuristic description of several other numerical schemes for analyzing this problem. An analysis of shock waves generated either by pressure mismatch between the injected fluid and surrounding flow or by chemical heat release is also described.

Dash, S.↗
NASA NTRS1977

A piecewise linear approximation scheme for hereditary optimal control problems

An approximation scheme based on 'piecewise linear' approximations of L2 spaces is employed to formulate a numerical method for solving quadratic optimal control problems governed by linear retarded functional differential equations. This piecewise linear method is an extension of the so called averaging technique. It is shown that the Riccati equation for the linear approximation is solved by simple transformation of the averaging solution. Thus, the computational requirements are essentially the same. Numerical results are given.

Cliff, E. M.↗
NASA NTRS1979

Simulation of three-dimensional compressible viscous flow on the Illiac IV computer

Complicated three-dimensional viscous transonic flows about bodies at high angles of attack are solved on the Illiac IV computer. It is shown that certain approximate forms of the compressible Reynolds-averaged Navier-Stokes equations can be computed about realistic three-dimensional geometries with relative ease on the Illiac IV. The ease and efficiency with which this can be done depend on the approximations made in the basic equations, the choice of the numerical algorithm used for the solution, and the data-base system that controls the data management and identifies and manipulates the vectors. A pencil data-base system is found to be particularly suitable for the approximations and numerical method chosen to produce the results presented. In addition, some comparisons are made of computer predictions with experimental results for various lows about hemisphere-cylinders in both subsonic and supersonic free streams. The same viscous model and numerical model are used, showing good qualitative agreement in the location of separation lines and pressure distributions.

Pulliam, T. H.↗
NASA NTRS1994

A Parallel Incompressible Navier-Stokes Solver with Multigrid Iterations

We developed a parallel, numerically accurate and stable, and computationally efficient finate-difference incompressible Navier-Stokes (N-S) fluid flow solver. The solver runs on both sequential and massively parallel computers. The numerical method used here is a second-order projection method on a staggered grid. The code is highly modular and it can be used either as a stand-alone flow solver and or a template code which can be adapted or expanded to a specific application. Numerical results and parallel performances of our code on Intel Delta and Paragon are reported.

Navier-Stokes solver↗
NASA NTRS2011

Obtaining the Grobner Initialization for the Ground Flash Fraction Retrieval Algorithm

At optical wavelengths and from the vantage point of space, the multiple scattering cloud medium obscures one's view and prevents one from easily determining what flashes strike the ground. However, recent investigations have made some progress examining the (easier, but still difficult) problem of estimating the ground flash fraction in a set of N flashes observed from space In the study by Koshak, a Bayesian inversion method was introduced for retrieving the fraction of ground flashes in a set of flashes observed from a (low earth orbiting or geostationary) satellite lightning imager. The method employed a constrained mixed exponential distribution model to describe the lightning optical measurements. To obtain the optimum model parameters, a scalar function of three variables (one of which is the ground flash fraction) was minimized by a numerical method. This method has formed the basis of a Ground Flash Fraction Retrieval Algorithm (GoFFRA) that is being tested as part of GOES-R GLM risk reduction.

Solakiewicz, R.↗
DOE OSTI2024

Anti-symmetric and positivity preserving formulation of a spectral method for Vlasov-Poisson equations

We analyze the anti-symmetric properties of a spectral discretization for the one-dimensional Vlasov-Poisson equations. The discretization is based on a spectral expansion in velocity with the symmetrically weighted Hermite basis functions, central finite differencing in space, and an implicit Runge Kutta integrator in time. The proposed discretization preserves the anti-symmetric structure of the advection operator in the Vlasov equation, resulting in a stable numerical method. We apply such discretization to two formulations: the canonical Vlasov-Poisson equations and their continuously transformed square-root representation. The latter preserves the positivity of the particle distribution function. We derive analytically the conservation properties of both formulations, including particle number, momentum, and energy, which are verified numerically on the following benchmark problems: manufactured solution, linear and nonlinear Landau damping, two-stream instability, bump-on-tail instability, and ion-acoustic wave.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2002

High Order Finite Difference Methods for Multiscale Complex Compressible Flows

The classical way of analyzing finite difference schemes for hyperbolic problems is to investigate as many as possible of the following points: (1) Linear stability for constant coefficients; (2) Linear stability for variable coefficients; (3) Non-linear stability; and (4) Stability at discontinuities. We will build a new numerical method, which satisfies all types of stability, by dealing with each of the points above step by step.

Sjoegreen, Bjoern↗