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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

A parallel iterative solution method for systems of nonlinear hyperbolic equations

An iterative algorithm suitable for the solution of a system of nonlinear hyperbolic partial differentiation equations in multiple dimensions is discussed. Current numerical methods for systems of nonlinear PDEs have limited parallelism due to strong coupling between the equations. This method decouples the PDEs by linearizing the convention coefficient for a space-time domain. This provides large grain parallelism. The linearization also allows the treatment of some terms in the equations as source terms, providing more freedom to choose from a wider variety of numerical methods. Smaller grain parallelism may be exploited within the solves for each equation. Thus, the method has potential for parallelism at several levels.

Scroggs, Jeffrey S.↗
NASA NTRS2001

B-spline Method in Fluid Dynamics

B-spline functions are bases for piecewise polynomials that possess attractive properties for complex flow simulations : they have compact support, provide a straightforward handling of boundary conditions and grid nonuniformities, and yield numerical schemes with high resolving power, where the order of accuracy is a mere input parameter. This paper reviews the progress made on the development and application of B-spline numerical methods to computational fluid dynamics problems. Basic B-spline approximation properties is investigated, and their relationship with conventional numerical methods is reviewed. Some fundamental developments towards efficient complex geometry spline methods are covered, such as local interpolation methods, fast solution algorithms on cartesian grid, non-conformal block-structured discretization, formulation of spline bases of higher continuity over triangulation, and treatment of pressure oscillations in Navier-Stokes equations. Application of some of these techniques to the computation of viscous incompressible flows is presented.

Botella, Olivier↗
NASA NTRS1984

Aerothermal environment of a blunted three-dimensional nonaxisymmetric body at Mach 6.8

Measurements of surface pressure and heat flux at Mach 6.8 for a nonaxisymmetric body, typical of the forward fuselage of an entry vehicle, have been obtained in the Langley 8-ft high-temperature tunnel. Comparisons of these data with predictions obtained from two numerical methods, which treat the test medium as an ideal gas with a laminar boundary layer, are presented. In general, the numerical results are consistent with each other and follow the data in regions of laminar flow; however, the data for angle of attack 14.9 deg indicate transitional flow on most of the windward surface. Also, the two numerical methods disagree slightly in regions of zero or nearly zero pressure gradients.

Martin, F. W., Jr.↗
NASA NTRS1990

Polyatomic molecular Dirac-Hartree-Fock calculations with Gaussian basis sets

Numerical methods have been used successfully in atomic Dirac-Hartree-Fock (DHF) calculations for many years. Some DHF calculations using numerical methods have been done on diatomic molecules, but while these serve a useful purpose for calibration, the computational effort in extending this approach to polyatomic molecules is prohibitive. An alternative more in line with traditional quantum chemistry is to use an analytical basis set expansion of the wave function. This approach fell into disrepute in the early 1980's due to problems with variational collapse and intruder states, but has recently been put on firm theoretical foundations. In particular, the problems of variational collapse are well understood, and prescriptions for avoiding the most serious failures have been developed. Consequently, it is now possible to develop reliable molecular programs using basis set methods. This paper describes such a program and reports results of test calculations to demonstrate the convergence and stability of the method.

Dyall, Kenneth G.↗
NASA NTRS1972

Analytic partial derivatives for estimating low-thrust parameters.

Analytic partial derivatives for estimating orbital low-thrust parameters via differential correction are developed and compared with two different numerical methods. The formulation is independent of the particular thrust model used and is applicable to all physically possible elliptic orbits. The starting point for the development is the set of variational equations of the elliptic orbital elements in the form due to Lagrange. The first time derivatives of the elements are transformed to derivatives with respect to the space variable, true anomaly, and integrated to first order in closed form in a straightforward general perturbations approach, with one exception: particular attention is given to the mean anomaly as influenced by thrust perturbations in the semimajor axis so that the complete first-order effect is included. The partials of the elements are then taken with respect to any given thrust parameter. Two comparisons are made with numerical methods for computing these thrust partials: numerical quotients and numerical integration of the variational equations for thrust.

Cunningham, G. W.↗
NASA NTRS1990

Optimal helicopter trajectory planning for terrain following flight

Helicopters operating in high threat areas have to fly close to the earth surface to minimize the risk of being detected by the adversaries. Techniques are presented for low altitude helicopter trajectory planning. These methods are based on optimal control theory and appear to be implementable onboard in realtime. Second order necessary conditions are obtained to provide a criterion for finding the optimal trajectory when more than one extremal passes through a given point. A second trajectory planning method incorporating a quadratic performance index is also discussed. Trajectory planning problem is formulated as a differential game. The objective is to synthesize optimal trajectories in the presence of an actively maneuvering adversary. Numerical methods for obtaining solutions to these problems are outlined. As an alternative to numerical method, feedback linearizing transformations are combined with the linear quadratic game results to synthesize explicit nonlinear feedback strategies for helicopter pursuit-evasion. Some of the trajectories generated from this research are evaluated on a six-degree-of-freedom helicopter simulation incorporating an advanced autopilot. The optimal trajectory planning methods presented are also useful for autonomous land vehicle guidance.

Menon, P. K. A.↗
NASA NTRS2013

Advanced Magnetic Materials Methods and Numerical Models for Fluidization in Microgravity and Hypogravity

To support long-duration manned missions in space such as a permanent lunar base, Mars transit, or Mars Surface Mission, improved methods for the treatment of solid wastes, particularly methods that recover valuable resources, are needed. The ability to operate under microgravity and hypogravity conditions is essential to meet this objective. The utilization of magnetic forces to manipulate granular magnetic media has provided the means to treat solid wastes under variable gravity conditions by filtration using a consolidated magnetic media bed followed by thermal processing of the solid wastes in a fluidized bed reactor. Non-uniform magnetic fields will produce a magnetic field gradient in a bed of magnetically susceptible media toward the distributor plate of a fluidized bed reactor. A correctly oriented magnetic field gradient will generate a downward direct force on magnetic media that can substitute for gravitational force in microgravity, or which may augment low levels of gravity, such as on the Moon or Mars. This approach is termed Gradient Magnetically Assisted Fluidization (G-MAFB), in which the magnitude of the force on the fluidized media depends upon the intensity of the magnetic field (H), the intensity of the field gradient (dH/dz), and the magnetic susceptibility of the media. Fluidized beds based on the G-MAFB process can operate in any gravitational environment by tuning the magnetic field appropriately. Magnetic materials and methods have been developed that enable G-MAFB operation under variable gravity conditions.

Atwater, James↗
NASA NTRS2013

Extended Analytic Device Optimization Employing Asymptotic Expansion

Analytic optimization of a thermoelectric junction often introduces several simplifying assumptionsincluding constant material properties, fixed known hot and cold shoe temperatures, and thermallyinsulated leg sides. In fact all of these simplifications will have an effect on device performance,ranging from negligible to significant depending on conditions. Numerical methods, such as FiniteElement Analysis or iterative techniques, are often used to perform more detailed analysis andaccount for these simplifications. While numerical methods may stand as a suitable solution scheme,they are weak in gaining physical understanding and only serve to optimize through iterativesearching techniques. Analytic and asymptotic expansion techniques can be used to solve thegoverning system of thermoelectric differential equations with fewer or less severe assumptionsthan the classic case. Analytic methods can provide meaningful closed form solutions and generatebetter physical understanding of the conditions for when simplifying assumptions may be valid.In obtaining the analytic solutions a set of dimensionless parameters, which characterize allthermoelectric couples, is formulated and provide the limiting cases for validating assumptions.Presentation includes optimization of both classic rectangular couples as well as practically andtheoretically interesting cylindrical couples using optimization parameters physically meaningful toa cylindrical couple. Solutions incorporate the physical behavior for i) thermal resistance of hot andcold shoes, ii) variable material properties with temperature, and iii) lateral heat transfer through legsides.

Thermoelectric Generators↗
NASA NTRS1987

An Exponential Finite Difference Technique for Solving Partial Differential Equations

An exponential finite difference algorithm, as first presented by Bhattacharya for one-dimensianal steady-state, heat conduction in Cartesian coordinates, has been extended. The finite difference algorithm developed was used to solve the diffusion equation in one-dimensional cylindrical coordinates and applied to two- and three-dimensional problems in Cartesian coordinates. The method was also used to solve nonlinear partial differential equations in one (Burger's equation) and two (Boundary Layer equations) dimensional Cartesian coordinates. Predicted results were compared to exact solutions where available, or to results obtained by other numerical methods. It was found that the exponential finite difference method produced results that were more accurate than those obtained by other numerical methods, especially during the initial transient portion of the solution. Other applications made using the exponential finite difference technique included unsteady one-dimensional heat transfer with temperature varying thermal conductivity and the development of the temperature field in a laminar Couette flow.

Handschuh, Robert F.↗
NASA NTRS1985

An implicit time-marching method for studying unsteady flow with massive separation

A fully implicit time-marching method is developed such that all spatial derivatives are approximated using central differences, but no use is made of any artificial dissipation. The numerical method solves the discretized equations using Alternating Direction Implicit-Block Gaussian Elimination technique. The method is implemented in the unsteady analysis, which solves the incompressible Navier-Stokes equations in terms of vorticity and stream function in generalized orthogonal coordinates. A clustered conformal C-grid is employed, and every effort is made to resolve the various length scales in the flow problem. The metric discontinuity at the branch-cut is treated appropriately using analytic continuation. Introduction of the BGE reordering permits implicit treatment of the branch cut in the numerical method. The vorticity singularity at the cusped trailing edge is also appropriately treated. This accurate and efficient implicit method is used to study flow at Re = 1000, past a 12-percent thick symmetric Joukowski airfoil at high angle of attack 30 and 53 deg.

Osswald, G. A.↗
NASA NTRS1973

Propagation of plane waves in flow through a variable area duct

Duct geometry is an important factor that influences the transmission of sound in a duct. Nonconstant area produces variations in steady flow quantities through the duct, which can cause reflections of acoustic disturbances and the creation of standing wave patterns, and which can attenuate and disperse propagating waves. To study the problem, a second-order-accurate numerical method has been developed. The method of characteristics is used to solve the acoustic equations with the assumption of quasi-one-dimensional flow in the duct. Numerical results are presented over a broad range of frequencies. Low-frequency results are compared with approximate solutions for compact nozzles and high-frequency calculations are compared with results from the short-wave theory. The numerical results are also compared with an analytical solution of the acoustic equations for exponential horn sections with zero mean flow to further establish the validity of the numerical method.

King, L. S.↗
NASA NTRS1997

The Accuracy of Shock Capturing in Two Spatial Dimensions

An assessment of the accuracy of shock capturing schemes is made for two-dimensional steady flow around a cylindrical projectile. Both a linear fourth-order method and a nonlinear third-order method are used in this study. It is shown, contrary to conventional wisdom, that captured two-dimensional shocks are asymptotically first-order, regardless of the design accuracy of the numerical method. The practical implications of this finding are discussed in the context of the efficacy of high-order numerical methods for discontinuous flows.

Carpenter, Mark H.↗
NASA NTRS2008

Neutron Transport Models and Methods for HZETRN and Coupling to Low Energy Light Ion Transport

Exposure estimates inside space vehicles, surface habitats, and high altitude aircraft exposed to space radiation are highly influenced by secondary neutron production. The deterministic transport code HZETRN has been identified as a reliable and efficient tool for such studies, but improvements to the underlying transport models and numerical methods are still necessary. In this paper, the forward-backward (FB) and directionally coupled forward-backward (DC) neutron transport models are derived, numerical methods for the FB model are reviewed, and a computationally efficient numerical solution is presented for the DC model. Both models are compared to the Monte Carlo codes HETCHEDS and FLUKA, and the DC model is shown to agree closely with the Monte Carlo results. Finally, it is found in the development of either model that the decoupling of low energy neutrons from the light ion (A<4) transport procedure adversely affects low energy light ion fluence spectra and exposure quantities. A first order correction is presented to resolve the problem, and it is shown to be both accurate and efficient.

Blattnig, S.R.↗
NASA NTRS1982

Development of a locally mass flux conservative computer code for calculating 3-D viscous flow in turbomachines

The VANS successive approximation numerical method was extended to the computation of three dimensional, viscous, transonic flows in turbomachines. A cross-sectional computer code, which conserves mass flux at each point of the cross-sectional surface of computation was developed. In the VANS numerical method, the cross-sectional computation follows a blade-to-blade calculation. Numerical calculations were made for an axial annular turbine cascade and a transonic, centrifugal impeller with splitter vanes. The subsonic turbine cascade computation was generated in blade-to-blade surface to evaluate the accuracy of the blade-to-blade mode of marching. Calculated blade pressures at the hub, mid, and tip radii of the cascade agreed with corresponding measurements. The transonic impeller computation was conducted to test the newly developed locally mass flux conservative cross-sectional computer code. Both blade-to-blade and cross sectional modes of calculation were implemented for this problem. A triplet point shock structure was computed in the inducer region of the impeller. In addition, time-averaged shroud static pressures generally agreed with measured shroud pressures. It is concluded that the blade-to-blade computation produces a useful engineering flow field in regions of subsonic relative flow; and cross-sectional computation, with a locally mass flux conservative continuity equation, is required to compute the shock waves in regions of supersonic relative flow.

Walitt, L.↗
NASA NTRS2003

The Langley Stability and Transition Analysis Code (LASTRAC) : LST, Linear and Nonlinear PSE for 2-D, Axisymmetric, and Infinite Swept Wing Boundary Layers

During the past two decades, our understanding of laminar-turbulent transition flow physics has advanced significantly owing to, in a large part, the NASA program support such as the National Aerospace Plane (NASP), High-speed Civil Transport (HSCT), and Advanced Subsonic Technology (AST). Experimental, theoretical, as well as computational efforts on various issues such as receptivity and linear and nonlinear evolution of instability waves take part in broadening our knowledge base for this intricate flow phenomenon. Despite all these advances, transition prediction remains a nontrivial task for engineers due to the lack of a widely available, robust, and efficient prediction tool. The design and development of the LASTRAC code is aimed at providing one such engineering tool that is easy to use and yet capable of dealing with a broad range of transition related issues. LASTRAC was written from scratch based on the state-of-the-art numerical methods for stability analysis and modem software technologies. At low fidelity, it allows users to perform linear stability analysis and N-factor transition correlation for a broad range of flow regimes and configurations by using either the linear stability theory (LST) or linear parabolized stability equations (LPSE) method. At high fidelity, users may use nonlinear PSE to track finite-amplitude disturbances until the skin friction rise. Coupled with the built-in receptivity model that is currently under development, the nonlinear PSE method offers a synergistic approach to predict transition onset for a given disturbance environment based on first principles. This paper describes the governing equations, numerical methods, code development, and case studies for the current release of LASTRAC. Practical applications of LASTRAC are demonstrated for linear stability calculations, N-factor transition correlation, non-linear breakdown simulations, and controls of stationary crossflow instability in supersonic swept wing boundary layers.

Chang, Chau-Lyan↗
NASA NTRS1983

Calculation of unsteady small disturbance transonic flow at arbitrary reduced frequency using an extended integral equation method.

A novel analytical-numerical method for calculating unsteady small disturbance transonic flow over airfoils has been developed. The method uses an extended integral equation technique, based on both the velocity potential and the acceleration potential, to predict unsteady aerodynamic loading on airfoils oscillating in subcritical transonic free stream conditions. The formulation is an extension of the work of Sivaneri and Harris (1980) for steady, non-lifting flows and utilizes the linear theory of Landahl (1961) for decoupling of steady and unsteady components. The analytical-numerical procedure involves several intnegrating schemes and applies to general frequencies of oscillations. The technique is illustrated by computing the transonic flow about parabolic arc airfoils. Specific unsteady results for reduced frequencies based on semi-chord of 0.01, 0.1, 0.3, 0.4 and 0.6 are given. Comparison of results with those obtained by an ADI finite difference scheme is made.

Lamah, C. A.↗
NASA NTRS2004

Langley Stability and Transition Analysis Code (LASTRAC) Version 1.2 User Manual

LASTRAC is a general-purposed, physics-based transition prediction code released by NASA for Laminar Flow Control studies and transition research. The design and development of the LASTRAC code is aimed at providing an engineering tool that is easy to use and yet capable of dealing with a broad range of transition related issues. It was written from scratch based on the state-of-the-art numerical methods for stability analysis and modern software technologies. At low fidelity, it allows users to perform linear stability analysis and N-factor transition correlation for a broad range of flow regimes and configurations by using either the linear stability theory or linear parabolized stability equations method. At high fidelity, users may use nonlinear PSE to track finite-amplitude disturbances until the skin friction rise. This document describes the governing equations, numerical methods, code development, detailed description of input/output parameters, and case studies for the current release of LASTRAC.

Chang, Chau-Lyan↗
NASA NTRS1994

Numerical solution of the two-dimensional time-dependent incompressible Euler equations

A numerical method is presented for solving the artificial compressibility form of the 2D time-dependent incompressible Euler equations. The approach is based on using an approximate Riemann solver for the cell face numerical flux of a finite volume discretization. Characteristic variable boundary conditions are developed and presented for all boundaries and in-flow out-flow situations. The system of algebraic equations is solved using the discretized Newton-relaxation (DNR) implicit method. Numerical results are presented for both steady and unsteady flow.

Whitfield, David L.↗