SEARCH · Search NASA

Results for “Methods: numerical”

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

All sources NASA NTRS DOE OSTI

NASA NTRS1981

Trends and pacing items in computational aerodynamics

A perspective is presented of trends in computational aerodynamics, and of important technology development items that pace future advanced applications. From a survey of AIAA Journal papers published during the past two decades, the growth trends and the progressively increasing emphasis on code development for viscous, compressible, turbulent flow are illustrated. These trends are reflected in the chronology of introduction by the aerospace industry of new computational methods in aircraft design. Key pacing items outlined are: automatic grid generation for nonlinear inviscid computations; advanced computers, improved efficiency of numerical methods, and improved turbulence models for Reynolds-averaged Navier-Stokes computations; advanced computers, time-dependent three-dimensional law-of-the-wall, code development, improved efficiency of numerical methods, and improved subgrid-scale turbulence modeling for large eddy simulations.

Chapman, D. R.↗

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.↗

DOE OSTI2023

KiT-RT: An Extendable Framework for Radiative Transfer and Therapy

Here, in this article, we present Kinetic Transport Solver for Radiation Therapy (KiT-RT), an open-source C++-based framework for solving kinetic equations in therapy applications available at https://github.com/CSMMLab/KiT-RT . This software framework aims to provide a collection of classical deterministic solvers for unstructured meshes that allow for easy extendability. Therefore, KiT-RT is a convenient base to test new numerical methods in various applications and compare them against conventional solvers. The implementation includes spherical harmonics, minimal entropy, neural minimal entropy, and discrete ordinates methods. Solution characteristics and efficiency are presented through several test cases ranging from radiation transport to electron radiation therapy. Due to the variety of included numerical methods and easy extendability, the presented open-source code is attractive for both developers, who want a basis to build their numerical solvers, and users or application engineers, who want to gain experimental insights without directly interfering with the codebase.

97 MATHEMATICS AND COMPUTING↗

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.↗

DOE OSTI2022

Numerical schemes for a multi-species BGK model with velocity-dependent collision frequency

Here, we consider a kinetic description of multi-species gas mixture modeled with Bhatnagar-Gross-Krook (BGK) collision operators, in which the collision frequency varies not only in time and space but also with the microscopic velocity. In this model, the Maxwellians typically used in standard BGK operators are replaced by a generalization of such target functions, which are defined by a variational procedure. In this paper we present a numerical method for simulating this model, which uses an Implicit-Explicit (IMEX) scheme to minimize a certain potential function, mimicking the Lagrange functional that appears in the theoretical derivation. We show that theoretical properties such as conservation of mass, total momentum and total energy as well as positivity of the distribution functions are preserved by the numerical method, and illustrate its usefulness and effectiveness with numerical examples.

97 MATHEMATICS AND COMPUTING↗

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.↗

DOE OSTI2023

Simulation of Molecular Hydrogen Jets and Visualization of Associated Hydrodynamic Features [Slides]

A numerical method is presented for solving the Navier-Stokes equations to represent a steady-state hydrogen jet expanding into a vacuum chamber. This numerical method has been evaluated using three test problems with promising results. Simulations of the hydrogen jet indicate the presence of a low temperature region (less than 15 K) within the jet for the high flow rate case. The high flow rate case also gives rise to a large recirculation region within the jet. The conditions at the minimum temperature point within the jet approach the hydrogen sublimation line when the mass flow rate exceeds 2 g/s, which may allow for the formation of hydrogen ice crystals.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗