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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Spectral methods for the Navier-Stokes equations with one infinite and two periodic directions

The time-dependent, 3D incompressible Navier-Stokes equations in (1) boundary layers, the semiinfinite domain, and (2) mixing layers or wakes, the fully infinite domain, are respectively solved by two numerical methods which employ rapidly decaying spectral basis functions to approximate the vertical dependence of the solutions. These functions are then combined with one, for method (1), and two, for method (2), slowly decaying 'extra functions' for each wave vector. Each extra function can exactly represent the solution's irrotational component at large distances. The two methods have been applied to extensive direct-numerical simulation of transition and turbulence.

Spalart, Philippe R.↗
NASA NTRS2013

Numerical Modeling of Ablation Heat Transfer

A unique numerical method has been developed for solving one-dimensional ablation heat transfer problems. This paper provides a comprehensive description of the method, along with detailed derivations of the governing equations. This methodology supports solutions for traditional ablation modeling including such effects as heat transfer, material decomposition, pyrolysis gas permeation and heat exchange, and thermochemical surface erosion. The numerical scheme utilizes a control-volume approach with a variable grid to account for surface movement. This method directly supports implementation of nontraditional models such as material swelling and mechanical erosion, extending capabilities for modeling complex ablation phenomena. Verifications of the numerical implementation are provided using analytical solutions, code comparisons, and the method of manufactured solutions. These verifications are used to demonstrate solution accuracy and proper error convergence rates. A simple demonstration of a mechanical erosion (spallation) model is also provided to illustrate the unique capabilities of the method.

Ewing, Mark E.↗
NASA NTRS1985

Large-scale computations in fluid mechanics; Proceedings of the Fifteenth Summer Seminar on Applied Mathematics, University of California, La Jolla, CA, June 27-July 8, 1983. Parts 1 & 2

Papers are presented on such topics as the use of semi-Lagrangian advective schemes in meteorological modeling; computation with high-resolution upwind schemes for hyperbolic equations; dynamics of flame propagation in a turbulent field; a modified finite element method for solving the incompressible Navier-Stokes equations; computational fusion magnetohydrodynamics; and a nonoscillatory shock capturing scheme using flux-limited dissipation. Consideration is also given to the use of spectral techniques in numerical weather prediction; numerical methods for the incorporation of mountains in atmospheric models; techniques for the numerical simulation of large-scale eddies in geophysical fluid dynamics; high-resolution TVD schemes using flux limiters; upwind-difference methods for aerodynamic problems governed by the Euler equations; and an MHD model of the earth's magnetosphere.

Engquist, B. E.↗
NASA NTRS1977

Displacement incrementation in non-linear structural analysis by the self-correcting method

The prediction of nonlinear structural behavior by the finite element method wherein buckling does not occur has received considerable attention and, with it, reasonable success has been achieved. However, the post-buckling problem has been less actively pursued probably because of the inherent numerical difficulties encountered. This note reviews very briefly the numerical methods currently being used for pre- and post-buckling analysis and presents a self-correcting approach based on load and displacement incrementation which is shown to be efficient, reliable, and easy to program. Numerical solutions are presented which demonstrate the effectiveness of the method.

Haisler, W. E.↗
NASA NTRS1983

Formulation and solution of rotary-wing aeroelastic stability and response problems

The state of the art in the formulation and solution of rotary-wing aeroelastic stability and response problems is reviewed in detail. The approximations used in the structural, inertia and aerodynamic operators are discussed. The important role of geometric nonlinearities, due to moderate deflections, and aerodynamic stall in the aeroelastic stability and response problem are identified. It is also shown that geometric nonlinearities are of primary importance in aeroelastic stability calculations, and have a more limited, though important, role in response calculations. Next, formulation of coupled rotor/fuselage problems is described, for both air and ground resonance type problems. Both topics, the isolated blade problem and the coupled rotor/fuselage problem, are treated for both hover and forward flight. Solution of aeroelastic stability and response problems proceeds in two stages. First, the spatial dependence is eliminated by using Galerkin's method, or by using the finite element method. Next the nonlinear, or linear, ordinary differential equation with periodic coefficients have to be solved for stability or response. Efficient numerical methods for accomplishing these objectives are presented in a comprehensive manner. The paper contains a number of illustrative numerical results which are intended to clarify various aspects of the modeling process and serve as representative results for both aeroelastic stability and response calculations for a variety of blade and rotor configurations.

Friedmann, P. P.↗
NASA NTRS1991

On the wall-normal velocity of the compressible boundary-layer equations

Numerical methods for the compressible boundary-layer equations are facilitated by transformation from the physical (x,y) plane to a computational (xi,eta) plane in which the evolution of the flow is 'slow' in the time-like xi direction. The commonly used Levy-Lees transformation results in a computationally well-behaved problem for a wide class of non-similar boundary-layer flows, but it complicates interpretation of the solution in physical space. Specifically, the transformation is inherently nonlinear, and the physical wall-normal velocity is transformed out of the problem and is not readily recovered. In light of recent research which shows mean-flow non-parallelism to significantly influence the stability of high-speed compressible flows, the contribution of the wall-normal velocity in the analysis of stability should not be routinely neglected. Conventional methods extract the wall-normal velocity in physical space from the continuity equation, using finite-difference techniques and interpolation procedures. The present spectrally-accurate method extracts the wall-normal velocity directly from the transformation itself, without interpolation, leaving the continuity equation free as a check on the quality of the solution. The present method for recovering wall-normal velocity, when used in conjunction with a highly-accurate spectral collocation method for solving the compressible boundary-layer equations, results in a discrete solution which is extraordinarily smooth and accurate, and which satisfies the continuity equation nearly to machine precision. These qualities make the method well suited to the computation of the non-parallel mean flows needed by spatial direct numerical simulations (DNS) and parabolized stability equation (PSE) approaches to the analysis of stability.

Pruett, C. David↗
NASA NTRS1993

Computation of unsteady shock-induced combustion using logarithmic species conservation equations

Numerical simulations are used to investigate periodic combustion instabilities observed in ballistic-range experiments of blunt bodies flying at supersonic speeds through hydrogen-air mixtures. The computations are validated by comparing experimental shadowgraphs with shadowgraphs created from the computed flowfields and by comparing the experimentally measured instability frequencies with computed frequencies. The numerical simulations use a logarithmic transformation of the species conservation equations as a way to reduce the grid requirements for computing shock-induced combustion. The transformation is applied to the Euler equations coupled to a detailed hydrogen-air chemical reaction mechanism with 13 species and 33 reactions. The resulting differential equations are solved using a finite volume formulation and a two-step predictor-corrector scheme to advance the solution in time. Results are presented and compared for both a flux-vector splitting scheme and an upwind TVD scheme. The computations add insight to the physical processes observed in the experiments and the numerical methods needed to simulate them. The usefulness of the ballistic-range experiments for the validation of numerical techniques and chemical kinetic models is also demonstrated.

Wilson, Gregory J.↗
NASA NTRS1999

Real Time Optima Tracking Using Harvesting Models of the Genetic Algorithm

Tracking optima in real time propulsion control, particularly for non-stationary optimization problems is a challenging task. Several approaches have been put forward for such a study including the numerical method called the genetic algorithm. In brief, this approach is built upon Darwinian-style competition between numerical alternatives displayed in the form of binary strings, or by analogy to 'pseudogenes'. Breeding of improved solution is an often cited parallel to natural selection in.evolutionary or soft computing. In this report we present our results of applying a novel model of a genetic algorithm for tracking optima in propulsion engineering and in real time control. We specialize the algorithm to mission profiling and planning optimizations, both to select reduced propulsion needs through trajectory planning and to explore time or fuel conservation strategies.

Baskaran, Subbiah↗
NASA NTRS2003

Simulations of 6-DOF Motion with a Cartesian Method

Coupled 6-DOF/CFD trajectory predictions using an automated Cartesian method are demonstrated by simulating a GBU-32/JDAM store separating from an F-18C aircraft. Numerical simulations are performed at two Mach numbers near the sonic speed, and compared with flight-test telemetry and photographic-derived data. Simulation results obtained with a sequential-static series of flow solutions are contrasted with results using a time-dependent flow solver. Both numerical methods show good agreement with the flight-test data through the first half of the simulations. The sequential-static and time-dependent methods diverge over the last half of the trajectory prediction. after the store produces peak angular rates. A cost comparison for the Cartesian method is included, in terms of absolute cost and relative to computing uncoupled 6-DOF trajectories. A detailed description of the 6-DOF method, as well as a verification of its accuracy, is provided in an appendix.

Murman, Scott M.↗
NASA NTRS1985

A comparison of finite-difference and finite-element methods for calculating free edge stresses in composites

It is pointed out that edge delamination is a serious failure mechanism for laminated composite materials. Various numerical methods have been utilized in attempts to calculate the interlaminar stress components which precede delamination in a laminate. There are, however, discrepancies regarding the results provided by different methods, taking into account a finite-difference procedure, a perturbation procedure, and finite element approaches. The present investigation has the objective to assess the capacity of a finite difference method to predict the character and magnitude of the interlaminar stress distributions near an interface corner. A second purpose of the investigation is to determine if predictions by finite element method in-plane, interlaminar stress components near an interface corner represent actual laminate behavior.

Bauld, N. R., Jr.↗
NASA NTRS1996

Multigrid and Krylov Subspace Methods for the Discrete Stokes Equations

Discretization of the Stokes equations produces a symmetric indefinite system of linear equations. For stable discretizations, a variety of numerical methods have been proposed that have rates of convergence independent of the mesh size used in the discretization. In this paper, we compare the performance of four such methods: variants of the Uzawa, preconditioned conjugate gradient, preconditioned conjugate residual, and multigrid methods, for solving several two-dimensional model problems. The results indicate that where it is applicable, multigrid with smoothing based on incomplete factorization is more efficient than the other methods, but typically by no more than a factor of two. The conjugate residual method has the advantage of being both independent of iteration parameters and widely applicable.

Elman, Howard C.↗
NASA NTRS1973

Transient response of rotor-bearing systems.

The equations of motion necessary to calculate the transient response of a multimass flexible rotor supported by nonlinear, damped bearings are derived from energy principles. Rotor excitation may be the result of imbalance, internal friction, rotor acceleration, nonlinear forces due to any number of bearing or seal stations, and gyroscopic couples developed from skewed disk effects. The method of solution for transient response simulation is discussed in detail and is based on extensive evaluation of numerical methods available for transient analysis. Examples of the application of transient response for the analysis of rotor bearing systems are presented and compared to actual machine performance. Recommendations for the use and extension of the present system simulation model are discussed.

Kirk, R. G.↗
NASA NTRS1988

Parameter estimation for boundary value problems by integral equations of the second kind

This paper is concerned with the parameter estimation for boundary integral equations of the second kind. The parameter estimation technique through use of the spline collocation method is proposed. Based on the compactness assumption imposed on the parameter space, the convergence analysis for the numerical method of parameter estimation is discussed. The results obtained here are applied to a boundary parameter estimation for 2-D elliptic systems.

Kojima, Fumio↗
NASA NTRS1991

Application of computational aeroacoustic methodologies to advanced propeller configurations - A review

Current research in the area of advanced propeller configurations for performance and acoustics are briefly reviewed. Particular attention is given to the techniques of Lock and Theodorsen modified for use in the design of counterrotating propeller configurations; a numerical method known as SSTAGE, which is a Euler solver for the unducted fan concept; the NASPROP-E numerical analysis also based on a Euler solver and used to study the near acoustic fields for the SR series propfan configurations; and a counterrotating propeller test rig designed to obtain an experimental performance/acoustic data base for various propeller configurations.

Korkan, Kenneth D.↗
NASA NTRS2000

Perspectives on the Future of CFD

This viewgraph presentation gives an overview of the future of computational fluid dynamics (CFD), which in the past has pioneered the field of flow simulation. Over time CFD has progressed as computing power. Numerical methods have been advanced as CPU and memory capacity increases. Complex configurations are routinely computed now and direct numerical simulations (DNS) and large eddy simulations (LES) are used to study turbulence. As the computing resources changed to parallel and distributed platforms, computer science aspects such as scalability (algorithmic and implementation) and portability and transparent codings have advanced. Examples of potential future (or current) challenges include risk assessment, limitations of the heuristic model, and the development of CFD and information technology (IT) tools.

Kwak, Dochan↗
NASA NTRS2006

Adjoint Sensitivity Computations for an Embedded-Boundary Cartesian Mesh Method and CAD Geometry

Cartesian-mesh methods are perhaps the most promising approach for addressing the issues of flow solution automation for aerodynamic design problems. In these methods, the discretization of the wetted surface is decoupled from that of the volume mesh. This not only enables fast and robust mesh generation for geometry of arbitrary complexity, but also facilitates access to geometry modeling and manipulation using parametric Computer-Aided Design (CAD) tools. Our goal is to combine the automation capabilities of Cartesian methods with an eficient computation of design sensitivities. We address this issue using the adjoint method, where the computational cost of the design sensitivities, or objective function gradients, is esseutially indepeudent of the number of design variables. In previous work, we presented an accurate and efficient algorithm for the solution of the adjoint Euler equations discretized on Cartesian meshes with embedded, cut-cell boundaries. Novel aspects of the algorithm included the computation of surface shape sensitivities for triangulations based on parametric-CAD models and the linearization of the coupling between the surface triangulation and the cut-cells. The objective of the present work is to extend our adjoint formulation to problems involving general shape changes. Central to this development is the computation of volume-mesh sensitivities to obtain a reliable approximation of the objective finction gradient. Motivated by the success of mesh-perturbation schemes commonly used in body-fitted unstructured formulations, we propose an approach based on a local linearization of a mesh-perturbation scheme similar to the spring analogy. This approach circumvents most of the difficulties that arise due to non-smooth changes in the cut-cell layer as the boundary shape evolves and provides a consistent approximation tot he exact gradient of the discretized abjective function. A detailed gradient accurace study is presented to verify our approach. Thereafter, we focus on a shape optimization problem for an Apollo-like reentry capsule. The optimization seeks to enhance the lift-to-drag ratio of the capsule by modifyjing the shape of its heat-shield in conjunction with a center-of-gravity (c.g.) offset. This multipoint and multi-objective optimization problem is used to demonstrate the overall effectiveness of the Cartesian adjoint method for addressing the issues of complex aerodynamic design. This abstract presents only a brief outline of the numerical method and results; full details will be given in the final paper.

Nemec, Marian↗
NASA NTRS2014

Methods for Prediction of High-Speed Reacting Flows in Aerospace Propulsion

Research to develop high-speed airbreathing aerospace propulsion systems was underway in the late 1950s. A major part of the effort involved the supersonic combustion ramjet, or scramjet, engine. Work had also begun to develop computational techniques for solving the equations governing the flow through a scramjet engine. However, scramjet technology and the computational methods to assist in its evolution would remain apart for another decade. The principal barrier was that the computational methods needed for engine evolution lacked the computer technology required for solving the discrete equations resulting from the numerical methods. Even today, computer resources remain a major pacing item in overcoming this barrier. Significant advances have been made over the past 35 years, however, in modeling the supersonic chemically reacting flow in a scramjet combustor. To see how scramjet development and the required computational tools finally merged, we briefly trace the evolution of the technology in both areas.

Drummond, J. Philip↗