SEARCH · Engineering Papers

Results for “numerical algorithms”

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 505 records · Page 28

All sources NASA NTRS DOE OSTI

NASA NTRS1991

Aerodynamic shape optimization of arbitrary hypersonic vehicles

A new method was developed to optimize, in terms of aerodynamic wave drag minimization, arbitrary (nonaxisymmetric) hypersonic vehicles in modified Newtonian flow, while maintaining the initial volume and length of the vehicle. This new method uses either a surface fitted Fourier series to represent the vehicle's geometry or an independent point motion algorithm. In either case, the coefficients of the Fourier series or the spatial locations of the points defining each cross section were varied and a numerical optimization algorithm based on a quasi-Newton gradient search concept was used to determine the new optimal configuration. Results indicate a significant decrease in aerodynamic wave drag for simple and complex geometries at relatively low CPU costs. In the case of a cone, the results agreed well with known analytical optimum ogive shapes. The procedure is capable of accepting more complex flow field analysis codes.

Dulikravich, George S.↗

NASA NTRS1995

Wave propagation and scattering in computational aeroacoustics

This paper describes predictions of model problems in computational aeroacoustics. Two problem classes are considered: The first (Workshop Category 3, Problem 1 and 2) considers two-dimensional wave propagation and non-reflecting boundary conditions in the presence of a mean flow. The second (Workshop Category 4, Problem 1) examines wall boundary conditions. For the last problem we introduce the Impedance Mismatch Method (IMM) to treat the solid wall boundaries. In this method the solid wall is simulated using a wall region in which the characteristic impedance is set to a different value to that in the fluid region. This method has advantages over traditional solid wall boundary conditions including simplicity of coding, speed of computations, and the ability to treat curved boundaries efficiently. Several numerical examples are given in addition to the Workshop Problems. The discretization of the Euler equations is performed in all cases with a Dispersion-Relation-Preserving (DRP) algorithm. The numerical results are compared with either analytical solutions or solutions obtained using traditional solid wall boundary conditions.

Chung, Cathy↗

NASA NTRS2017

Fast Optimization for Aircraft Descent and Approach Trajectory

We address problem of on-line scheduling of the aircraft descent and approach trajectory. We formulate a general multiphase optimal control problem for optimization of the descent trajectory and review available methods of its solution. We develop a fast algorithm for solution of this problem using two key components: (i) fast inference of the dynamical and control variables of the descending trajectory from the low dimensional flight profile data and (ii) efficient local search for the resulting reduced dimensionality non-linear optimization problem. We compare the performance of the proposed algorithm with numerical solution obtained using optimal control toolbox General Pseudospectral Optimal Control Software. We present results of the solution of the scheduling problem for aircraft descent using novel fast algorithm and discuss its future applications.

Luchinsky, Dmitry G.↗

NASA NTRS1985

Numerical synthesis of tri-variate velocity realizations of turbulence

An approach for synthesizing trivariate turbulence velocity field spatial realizations is presented. Some of the spatial frequency characteristics of the random velocity field are described by the von Karman spectrum. The simulation algorithm is based on an efficient autoregressive-moving average (ARMA) scheme involving coefficient square matrices of order three. The determination of the efficient low order ARMA algorithm is preceded by the determination of a suitable high order autoregressive (AR) simulation algorithm. The numerical results are presented in a dimensionless form. Thus, they are applicable for any scale of turbulence.

Spanos, P.-T. D.↗

NASA NTRS1979

A fast, conservative algorithm for solving the transonic full-potential equation

A fast, fully implicit approximate factorization (AF) algorithm designed to solve the conservative transonic full-potential equation in either two or three dimensions is described. The algorithm uses an upwind bias of the density coefficient for stability in supersonic regions. This provides an effective upwind difference of the streamwise terms for any orientation of the velocity vector (i.e., 'rotated differencing'), and thereby greatly enhances the reliability of the present algorithm. A numerical transformation is used to establish an arbitrary body-fitted finite-difference mesh. Computed results for both airfoils and simplified wings demonstrate substantial improvement in convergence speed for the new algorithm relative to standard successive-line overrelaxation algorithms.

Holst, T. L.↗

NASA NTRS1980

Fast, Conservative Algorithm for Solving the Transonic Full-Potential Equation

A fast, fully implicit approximate factorization algorithm designed to solve the conservative, transonic, full-potential equation in either two or three dimensions is described. The algorithm uses an upwind bias of the density coefficient for stability in supersonic regions. This provides an effective upwind difference of the streamwise terms for any orientation of the velocity vector (i.e., rotated differencing), thereby greatly enhancing the reliability of the present algorithm. A numerical transformation is used to establish an arbitrary body-fitted, finite-difference mesh. Computed results for both airfoils and simplified wings demonstrate substantial improvement in convergence speed for the new algorithm relative to standard successive-line over-relaxation algorithms.

Holst, Terry L.↗

NASA NTRS2020

The AERONET Version 3 Aerosol Retrieval Algorithm, Associated Uncertainties and Comparisons to Version 2

The Aerosol Robotic Network (AERONET) Version 3 (V3) aerosol retrieval algorithm is described, which is based on the Version 2 (V2) algorithm with numerous updates. Comparisons of V3 aerosol retrievals to those of V2 are presented, along with a new approach to estimate uncertainties in many of the retrieved aerosol parameters. Changes in the V3 aerosol retrieval algorithm include (1) a new polarized radiative transfer code (RTC), which replaced the scalar RTC of V2, (2) detailed characterization of gas absorption by adding NO2 and H2O to specify total gas absorption in the atmospheric column, specification of vertical profiles of all the atmospheric species, (3) new bidirectional reflectance distribution function (BRDF) parameters for land sites adopted from the MODIS BRDF/Albedo product, (4) a new version of the extraterrestrial solar flux spectrum, and (5) a new temperature correction procedure of both direct Sun and sky radiance measurements. The potential effect of each change in V3 on single scattering albedo (SSA) retrievals was analyzed. The operational almucantar retrievals of V2 versus V3 were compared for four AERONET sites: GSFC, Mezaira, Mongu, and Kanpur. Analysis showed very good agreement in retrieved parameters of the size distributions. Comparisons of SSA retrievals for dust aerosols (Mezaira) showed a good agreement in 440 nm SSA, while for longer wavelengths V3 SSAs are systematically higher than those of V2, with the largest mean difference at 675 nm due to cumulative effects of both extraterrestrial solar flux and BRDF changes. For non-dust aerosols, the largest SSA deviation is at 675 nm due to differences in extraterrestrial solar flux spectrums used in each version. Further, the SSA 675 nm mean differences are very different for weakly (GSFC) and strongly (Mongu) absorbing aerosols, which is explained by the lower sensitivity to a bias in aerosol scattering optical depth by less absorbing aerosols. A new hybrid (HYB) sky radiance measurement scan is introduced and discussed. The HYB combines features of scans in two different planes to maximize the range of scattering angles and achieve scan symmetry, thereby allowing for cloud screening and spatial averaging, which is an advantage over the principal plane scan that lacks robust symmetry. We show that due to an extended range of scattering angles, HYB SSA retrievals for dust aerosols exhibit smaller variability with solar zenith angles (SZAs) than those of almucantar (ALM), which allows extension of HYB SSA retrievals to SZAs less than 50∘ to as small as 25∘. The comparison of SSA retrievals from closely time-matched HYB and ALM scans in the 50 to 75∘ SZA range showed good agreement with the differences below ∼0.005. We also present an approach to estimate retrieval uncertainties which utilizes the variability in retrieved parameters generated by perturbing both measurements and auxiliary input parameters as a proxy for retrieval uncertainty. The perturbations in measurements and auxiliary inputs are assumed as estimated biases in aerosol optical depth (AOD), radiometric calibration of sky radiances combined with solar spectral irradiance, and surface reflectance. For each set of Level 2 Sun/sky radiometer observations, 27 inputs corresponding to 27 combinations of biases were produced and separately inverted to generate the following statistics of the inversion results: average, standard deviation, minimum and maximum values. From these statistics, standard deviation (labeled U27) is used as a proxy for estimated uncertainty, and a lookup table (LUT) approach was implemented to reduce the computational time. The U27 climatological LUT was generated from the entire AERONET almucantar (1993–2018) and hybrid (2014–2018) scan databases by binning U27s in AOD (440 nm), Angström exponent (AE, 440–870 nm), and SSA (440, 675, 870, 1020 nm). Using this LUT approach, the uncertainty estimates U27 for each individual V3 Level 2 retrieval can be obtained by interpolation using the corresponding measured and inverted combination of AOD, AE, and SSA.

Aliaksandr Siniuk↗

NASA NTRS2019

Pterodactyl: An Uncoupled Range Control Approach to Fully Numerical Predictor-Corrector Entry Guidance

Entry, descent, and landing (EDL) has been identified as a core area of investment in NASA's Strategic Technology Investment Plan (NASA STIP). STIP lists the space technologies needed to help achieve NASA's science, technology, and exploration goals across the agency. Within the EDL core area, deployable hypersonic decelerators, also known as deployable entry vehicles (DEVs), have been identified as an area of investment, due to its potential to revolutionize payload delivery methods to Earth and other planets. These vehicles, which can deploy their heat shields or alter their shape before entry, exploit an increased and more effective drag ratio by using less mass than traditional blunt body vehicles with rigid aeroshells. DEVs like Adaptive Deployable Entry and Placement Technology (ADEPT) and Hypersonic Inflatable Aerodynamic Decelerator (HIAD) have demonstrated the capability of transporting the equivalent science payloads of blunt body rigid aeroshells, while using a significantly smaller diameter when stowed within a launch vehicle. While DEVs' increased energy dissipation for less mass is an attractive feature, their ability to contract and expand would require advancements in the current state-of-the-art guidance and control (G&C) architectures used by traditional rigid vehicles. Pterodactyl, a project funded by NASA's Space Technology Mission Directorate (STMD), aims to provide feasible integrated G&C solutions for DEVs, complete with optimized vehicle designs and packaging analyses. Structural and aerodynamic analyses for the explored control systems suggested a need for a bank angle guidance algorithm, a heritage guidance approach that has been used in many entry precision targeting vehicles, as well as an additional need for the development of a non-bank angle guidance. For this reason, Pterodactyl will consider four different G&C configurations during its design phase: i) a reaction control system for bank (sigma) control, ii) a mass movement system for angle of attack (alpha) sideslip (beta) control, iii) flaps for alpha - beta control, and iv) flaps for sigma control. To increase the applicability of each proposed integrated G&C architecture, an 11 km/s lunar return demonstration mission is selected to stress the developed technology capability. The Lifting Nano-ADEPT (LNA) vehicle is chosen as the DEV to demonstrate the integrated solutions. This paper will detail the trajectory design for a lunar return mission, using the validated bank control guidance algorithm Fully Numerical Predictor-Corrector Entry Guidance (FNPEG) and a newly developed guidance algorithm: FNPEG Uncoupled Range Control (URC). FNPEG-URC diverges from traditional bank angle guidances by producing alpha and beta commands to thereby decouple downrange and crossrange control. This presentation will discuss the development and overall performance of FNPEG and FNPEG-URC for each of the four G&C configurations. Successful G&C configurations are defined as those that can deliver payloads to the intended descent and landing site while abiding by trajectory constraints in the face of dispersions.

Johnson, Breanna↗

NASA NTRS1976

A numerical comparison of discrete Kalman filtering algorithms - An orbit determination case study

An improved Kalman filter algorithm based on a modified Givens matrix triangularization technique is proposed for solving a nonstationary discrete-time linear filtering problem. The proposed U-D covariance factorization filter uses orthogonal transformation technique; measurement and time updating of the U-D factors involve separate application of Gentleman's fast square-root-free Givens rotations. Numerical stability and accuracy of the algorithm are compared with those of the conventional and stabilized Kalman filters and the Potter-Schmidt square-root filter, by applying these techniques to a realistic planetary navigation problem (orbit determination for the Saturn approach phase of the Mariner Jupiter-Saturn Mission, 1977). The new algorithm is shown to combine the numerical precision of square root filtering with the efficiency of the original Kalman algorithm.

Thornton, C. L.↗

NASA NTRS1991

An unconditionally stable staggered algorithm for transient finite element analysis of coupled thermoelastic problems

An unconditionally stable second order accurate implicit-implicit staggered procedure for the finite element solution of fully coupled thermoelasticity transient problems is proposed. The procedure is stabilized with a semi-algebraic augmentation technique. A comparative cost analysis reveals the superiority of the proposed computational strategy to other conventional staggered procedures. Numerical examples of one and two-dimensional thermomechanical coupled problems demonstrate the accuracy of the proposed numerical solution algorithm.

Farhat, C.↗

NASA NTRS1991

An unconditionally stable staggered algorithm for transient finite element analysis of coupled thermoelastic problems

An unconditionally stable second order accurate implicit-implicit staggered procedure for the finite element solution of fully coupled thermoelasticity transient problems is proposed. The procedure is stabilized with a semi-algebraic augmentation technique. A comparative cost analysis reveals the superiority of the proposed computational strategy to other conventional staggered procedures. Numerical examples of one- and two-dimensional thermomechanical coupled problems demonstrate the accuracy of the proposed numerical solution algorithm.

Farhat, Charbel↗

NASA NTRS2012

Adaptive Aft Signature Shaping of a Low-Boom Supersonic Aircraft Using Off-Body Pressures

The design and optimization of a low-boom supersonic aircraft using the state-of-the- art o -body aerodynamics and sonic boom analysis has long been a challenging problem. The focus of this paper is to demonstrate an e ective geometry parameterization scheme and a numerical optimization approach for the aft shaping of a low-boom supersonic aircraft using o -body pressure calculations. A gradient-based numerical optimization algorithm that models the objective and constraints as response surface equations is used to drive the aft ground signature toward a ramp shape. The design objective is the minimization of the variation between the ground signature and the target signature subject to several geometric and signature constraints. The target signature is computed by using a least-squares regression of the aft portion of the ground signature. The parameterization and the deformation of the geometry is performed with a NASA in- house shaping tool. The optimization algorithm uses the shaping tool to drive the geometric deformation of a horizontal tail with a parameterization scheme that consists of seven camber design variables and an additional design variable that describes the spanwise location of the midspan section. The demonstration cases show that numerical optimization using the state-of-the-art o -body aerodynamic calculations is not only feasible and repeatable but also allows the exploration of complex design spaces for which a knowledge-based design method becomes less effective.

Ordaz, Irian↗

NASA NTRS1975

A general algorithm using finite element method for aerodynamic configurations at low speeds

A finite element algorithm for numerical simulation of two-dimensional, incompressible, viscous flows was developed. The Navier-Stokes equations are suitably modelled to facilitate direct solution for the essential flow parameters. A leap-frog time differencing and Galerkin minimization of these model equations yields the finite element algorithm. The finite elements are triangular with bicubic shape functions approximating the solution space. The finite element matrices are unsymmetrically banded to facilitate savings in storage. An unsymmetric L-U decomposition is performed on the finite element matrices to obtain the solution for the boundary value problem.

Balasubramanian, R.↗

NASA NTRS1991

Optimal design of solidification processes

An optimal design algorithm is presented for the analysis of general solidification processes, and is demonstrated for the growth of GaAs crystals in a Bridgman furnace. The system is optimal in the sense that the prespecified temperature distribution in the solidifying materials is obtained to maximize product quality. The optimization uses traditional numerical programming techniques which require the evaluation of cost and constraint functions and their sensitivities. The finite element method is incorporated to analyze the crystal solidification problem, evaluate the cost and constraint functions, and compute the sensitivities. These techniques are demonstrated in the crystal growth application by determining an optimal furnace wall temperature distribution to obtain the desired temperature profile in the crystal, and hence to maximize the crystal's quality. Several numerical optimization algorithms are studied to determine the proper convergence criteria, effective 1-D search strategies, appropriate forms of the cost and constraint functions, etc. In particular, we incorporate the conjugate gradient and quasi-Newton methods for unconstrained problems. The efficiency and effectiveness of each algorithm is presented in the example problem.

Dantzig, Jonathan A.↗

NASA NTRS1989

Navier-Stokes analysis of solid propellant rocket motor internal flows

A multidimensional implicit Navier-Stokes analysis that uses numerical solution of the ensemble-averaged Navier-Stokes equations in a nonorthogonal, body-fitted, cylindrical coordinate system has been applied to the simulation of the steady mean flow in solid propellant rocket motor chambers. The calculation procedure incorporates a two-equation (k-epsilon) turbulence model and utilizes a consistently split, linearized block-implicit algorithm for numerical solution of the governing equations. The code was validated by comparing computed results with the experimental data obtained in cylindrical-port cold-flow tests. The agreement between the computed and experimentally measured mean axial velocities is excellent. The axial location of transition to turbulent flow predicted by the two-equation (k-epsilon) turbulence model used in the computations also agrees well with the experimental data. Computations performed to simulate the axisymmetric flowfield in the vicinity of the aft field joint in the Space Shuttle solid rocket motor using 14,725 grid points show the presence of a region of reversed axial flow near the downstream edge of the slot.

Sabnis, J. S.↗

NASA NTRS1987

Navier-Stokes analysis of two- and three-dimensional flow field in solid rocket motors with segment joints

A multidimensional implicit Navier-Stokes analysis which uses numerical solution of ensemble-averaged Navier-Stokes equations in a nonorthogonal bodyfitted cylindrical-polar coordinate system has been applied to simulation of the internal flow field in solid-propellant rocket motor chambers with segment joints. The calculation procedure incorporates a two-equation (k-epsilon) turbulence model and utilizes a consistently split, linearized block-implicit algorithm for numerical solution of the governing equations. Computations performed to simulate the axisymmetric flow field in the vicinity of the aft field joint in the Space Shuttle SRB using 14,725 grid points show the presence of a region of reversed axial flow near the downstream edge of the slot. Calculations were also performed for two cases involving asymmetric three-dimensional flow in the vicinity of the aft field joint in the SRB using 721,525 grid points to estimate circumferential velocities and pressure gradients at the joint.

Sabnis, J. S.↗

NASA NTRS1985

Expert systems and their use in augmenting design optimization

The challenging requirements that are evolving for future aircraft demand that each design be optimally integrated, for the penalties imposed by nonoptimal performance are significant. Classic numerical optimization algorithms have been and will continue to be important tools for aircraft designers. These methods are, however, limited to certain categories of aircraft design variables, leaving the remainder to be determined by the user. A method that makes use of knowledge-based expert systems offers the potential for aiding the conceptual design process in a way that is similar to that of numerical optimization, except that it would address discrete, discontinuous, abstract, or any other unoptimized aspect of vehicle design and integration. Other unique capabilities such as automatic discovery and learning in design may also be achievable in the near term. This paper discusses current practice in conceptual aircraft design and knowledge-based systems, and how knowledge-based systems can be used in conceptual design.

Kidwell, G. H.↗

NASA NTRS1981

Numerical simulation of steady supersonic flow

A noniterative, implicit, space-marching, finite-difference algorithm was developed for the steady thin-layer Navier-Stokes equations in conservation-law form. The numerical algorithm is applicable to steady supersonic viscous flow over bodies of arbitrary shape. In addition, the same code can be used to compute supersonic inviscid flow or three-dimensional boundary layers. Computed results from two-dimensional and three-dimensional versions of the numerical algorithm are in good agreement with those obtained from more costly time-marching techniques.

Schiff, L. B.↗