Search NASA⌕ Search

SEARCH · Search NASA

Results for “numerical algorithm”

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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Parabolized Navier-Stokes methods for hypersonic flows

The development and numerical solution of the Parabolized Navier-Stokes (PNS) equations is discussed. Special emphasis is placed on numerical integration algorithms for PNS solvers, most notably recent applications of upwind schemes. Application of this technology to high speed, high Reynolds number flows was the focus, and modeling techniques for real gas effects and turbulence were discussed. Experiences in validating the UPS upwind PNS flow solver for relatively simple body shapes are presented. The test cases range in complexity from two dimensional perfect gas cases to three dimensional finite rate chemistry cases.

Lawrence, S. L.↗
NASA NTRS1993

Stability of mixing layers

The research program for the first year of this project (see the original research proposal) consists of developing an explicit marching scheme for solving the parabolized stability equations (PSE). Performing mathematical analysis of the computational algorithm including numerical stability analysis and the determination of the proper boundary conditions needed at the boundary of the computation domain are implicit in the task. Before one can solve the parabolized stability equations for high-speed mixing layers, the mean flow must first be found. In the past, instability analysis of high-speed mixing layer has mostly been performed on mean flow profiles calculated by the boundary layer equations. In carrying out this project, it is believed that the boundary layer equations might not give an accurate enough nonparallel, nonlinear mean flow needed for parabolized stability analysis. A more accurate mean flow can, however, be found by solving the parabolized Navier-Stokes equations. The advantage of the parabolized Navier-Stokes equations is that its accuracy is consistent with the PSE method. Furthermore, the method of solution is similar. Hence, the major part of the effort of the work of this year has been devoted to the development of an explicit numerical marching scheme for the solution of the Parabolized Navier-Stokes equation as applied to the high-seed mixing layer problem.

Tam, Christopher↗
NASA NTRS1993

Optimal design of composite hip implants using NASA technology

Using an adaptation of NASA software, we have investigated the use of numerical optimization techniques for the shape and material optimization of fiber composite hip implants. The original NASA inhouse codes, were originally developed for the optimization of aerospace structures. The adapted code, which was called OPORIM, couples numerical optimization algorithms with finite element analysis and composite laminate theory to perform design optimization using both shape and material design variables. The external and internal geometry of the implant and the surrounding bone is described with quintic spline curves. This geometric representation is then used to create an equivalent 2-D finite element model of the structure. Using laminate theory and the 3-D geometric information, equivalent stiffnesses are generated for each element of the 2-D finite element model, so that the 3-D stiffness of the structure can be approximated. The geometric information to construct the model of the femur was obtained from a CT scan. A variety of test cases were examined, incorporating several implant constructions and design variable sets. Typically the code was able to produce optimized shape and/or material parameters which substantially reduced stress concentrations in the bone adjacent of the implant. The results indicate that this technology can provide meaningful insight into the design of fiber composite hip implants.

Blake, T. A.↗
NASA NTRS1993

Radio jets in colliding galaxies: Testing the interaction-activity connection

The idea that galaxy interactions and merging are related to the generation of starburst and AGN activity in galactic nuclei has been the subject of intensive investigations over the past several years and is still a matter of lively debate. Peculiar morphologies, indicative of tidal interactions, have been detected in high-luminosity radio galaxies, in quasars, and in ultraluminous IRAS galaxies. In addition, low-luminosity radio and active galaxies show similar evidence for a recent merger or for nearby companions. In a recent CCD optical study of galaxies selected on the basis that they all contain well defined radio jets, it was found that almost half of the sample consists of pairs of elliptical galaxies. Many of these low-luminosity radio galaxies with companions show a well defined distorted radio jet structure at the VLA scale with an S- or C-shaped morphology. We have developed a general numerical simulation algorithm for ballistic radio jets with the intention of applying this model to the study of the bent jets seen in colliding pairs of galaxies and with the hope of testing the well documented interaction-activity connection. In our model the morphological evolution of the jets is determined by their response to the simple mechanical forces (i.e., gravity and ram pressure) imposed on them from both the host and the companion galaxies. Radiative losses, jet precession, magnetic effects, relativistic terms, and hydrodynamic instabilities have all been ignored. Starting with a previously derived collision model for the interacting pair of elliptical galaxies NGC 4782/4783, we have used our algorithm to simulate the specific two-sided jet morphology seen in the radio source 3C 278, associated with NGC 4782. This is the first time that such jet simulations have been produced for a galaxy pair whose relative orbit was determined independently from the jet modeling.

Borne, Kirk D.↗
NASA NTRS1992

Parallel-vector design sensitivity analysis in structural dynamics

This paper presents a parallel-vector algorithm for sensitivity calculations in linear structural dynamics. The proposed alternative formulation works efficiently with the reduced system of dynamic equations, since it eliminates the need for expensive and complicated based-vector derivatives, which are required in the conventional reduced system formulation. The relationship between the alternative formulation and the conventional reduced system formulation has been established, and it has been proven analytically that the two approaches are identical when all the mode shapes are included. This paper validates the proposed alternative algorithm through numerical experiments, where only a small number of mode shapes are used. In addition, a modified mode acceleration method is presented, thus not only the displacements but also the velocities and accelerations are shown to be improved.

Zhang, Y.↗
NASA NTRS1992

Accurate ab initio quartic force fields for borane and BeH2

The quartic force fields of BH3 and BeH2 have been computed ab initio using an augmented coupled cluster (CCSD(T)) method and basis sets of spdf and spdfg quality. For BH3, the computed spectroscopic constants are in very good agreement with recent experimental data, and definitively confirm misassignments in some older work, in agreement with recent ab initio studies. Using the computed spectroscopic constants, the rovibrational partition function for both molecules has been constructed using a modified direct numerical summation algorithm, and JANAF-style thermochemical tables are presented.

Martin, J. M. L.↗
NASA NTRS1993

Symmetric Galerkin boundary formulations employing curved elements

Accounts of the symmetric Galerkin approach to boundary element analysis (BEA) have recently been published. This paper attempts to add to the understanding of this method by addressing a series of fundamental issues associated with its potential computational efficiency. A new symmetric Galerkin theoretical formulation for both the (harmonic) heat conduction and the (biharmonic) elasticity problem that employs regularized singular and hypersingular boundary integral equations (BIEs) is presented. The novel use of regularized BIEs in the Galerkin context is shown to allow straightforward incorporation of curved, isoparametric elements. A symmetric reusable intrinsic sample point (RISP) numerical integration algorithm is shown to produce a Galerkin (i.e., double) integration strategy that is competitive with its counterpart (i.e., singular) integration procedure in the collocation BEA approach when the time saved in the symmetric equation solution phase is also taken into account. This new formulation is shown to be capable of employing hypersingular BIEs while obviating the requirement of C 1 continuity, a fact that allows the employment of the popular continuous element technology. The behavior of the symmetric Galerkin BEA method with regard to both direct and iterative equation solution operations is also addressed. A series of example problems are presented to quantify the performance of this symmetric approach, relative to the more conventional unsymmetric BEA, in terms of both accuracy and efficiency. It is concluded that appropriate implementations of the symmetric Galerkin approach to BEA indeed have the potential to be competitive with, if not superior to, collocation-based BEA, for large-scale problems.

Kane, J. H.↗
NASA NTRS1990

On domains of convergence in optimization problems

Numerical optimization algorithms require the knowledge of an initial set of design variables. Starting from an initial design x(sup 0), improved solutions are obtained by updating the design iteratively in a way prescribed by the particular algorithm used. If the algorithm is successful, convergence is achieved to a local optimal solution. Let A denote the iterative procedure that characterizes a typical optimization algorithm, applied to the problem: Find x belonging to R(sup n) that maximizes f(x) subject to x belonging to Omega contained in R(sup n). We are interested in problems with several local maxima (x(sub j))(sup *), j=1, ..., m, in the feasible design space Omega. In general, convergence of the algorithm A to a specific solution (x(sub j))(sup *) is determined by the choice of initial design x(sup 0). The domain of convergence D(sub j) of A associated with a local maximum (x(sub j))(sup *) is a subset of initial designs x(sup 0) in Omega such that the sequence (x(sup k)), k=0,1,2,... defined by x(sup k+1) = A(x(sup k)), k=0,1,... converges to (x(sub j))(sup *). The set D(sub j) is also called the basin of attraction of (x(sub j))(sup *). Cayley first proposed the problem of finding the basin of attraction for Newton's method in 1897. It has been shown that the basin of attraction for Newton's method exhibits chaotic behavior in problems with polynomial objective. This implies that there may be regions in the feasible design space where arbitrarily close starting points will converge to different local optimal solutions. Furthermore, the boundaries of the domains of convergence may have a very complex, even fractal structure. In this paper we show that even simple structural optimization problems solved using standard gradient based (first order) algorithms exhibit similar features.

Diaz, Alejandro R.↗
NASA NTRS1994

Design and implementation of a parallel unstructured Euler solver using software primitives

This paper is concerned with the implementation of a three-dimensional unstructured-grid Euler solver on massively parallel distributed-memory computer architectures. The goal is to minimize solution time by achieving high computational rates with a numerically efficient algorithm. An unstructured multigrid algorithm with an edge-based data structure has been adopted, and a number of optimizations have been devised and implemented to accelerate the parallel computational rates. The implementation is carried out by creating a set of software tools, which provide an interface between the parallelization issues and the sequential code, while providing a basis for future automatic run-time compilation support. Large practical unstructured grid problems are solved on the Intel iPSC/860 hypercube and Intel Touchstone Delta machine. The quantitative effects of the various optimizations are demonstrated, and we show that the combined effect of these optimizations leads to roughly a factor of 3 performance improvement. The overall solution efficiency is compared with that obtained on the Cray Y-MP vector supercomputer.

Das, R.↗
NASA NTRS1996

Aerodynamic shape optimization using control theory

Aerodynamic shape design has long persisted as a difficult scientific challenge due its highly nonlinear flow physics and daunting geometric complexity. However, with the emergence of Computational Fluid Dynamics (CFD) it has become possible to make accurate predictions of flows which are not dominated by viscous effects. It is thus worthwhile to explore the extension of CFD methods for flow analysis to the treatment of aerodynamic shape design. Two new aerodynamic shape design methods are developed which combine existing CFD technology, optimal control theory, and numerical optimization techniques. Flow analysis methods for the potential flow equation and the Euler equations form the basis of the two respective design methods. In each case, optimal control theory is used to derive the adjoint differential equations, the solution of which provides the necessary gradient information to a numerical optimization method much more efficiently then by conventional finite differencing. Each technique uses a quasi-Newton numerical optimization algorithm to drive an aerodynamic objective function toward a minimum. An analytic grid perturbation method is developed to modify body fitted meshes to accommodate shape changes during the design process. Both Hicks-Henne perturbation functions and B-spline control points are explored as suitable design variables. The new methods prove to be computationally efficient and robust, and can be used for practical airfoil design including geometric and aerodynamic constraints. Objective functions are chosen to allow both inverse design to a target pressure distribution and wave drag minimization. Several design cases are presented for each method illustrating its practicality and efficiency. These include non-lifting and lifting airfoils operating at both subsonic and transonic conditions.

Reuther, James↗
NASA NTRS1996

Development of a Vorticity-Velocity Navier-Stokes Formulation for the Study of Compressibility Effects on Dynamic Stall

The first major area of this study was to develop a vorticity-velocity formulation and numerical solution algorithms suitable for the analyses of incompressible as well as low-to- moderate-speed compressible flows. Research performed towards contributing to the determination of the appropriate vorticity and dilation creation boundary conditions suggested to temporarily set aside this approach and use a primitive-variable approach other than the pseudo-compressibility approach used. The second major area of study was initiated to comprehensively examine the INS-2D and INS-3D programs from the point of view of boundary conditions. The research carried out was documented in the form of two technical papers which are included in Appendices A and B; the boundary-condition related issues for INS-3D are briefly mentioned.

Ghia, K. N.↗
NASA NTRS1999

The Krigifier: A Procedure for Generating Pseudorandom Nonlinear Objective Functions for Computational Experimentation

Comprehensive computational experiments to assess the performance of algorithms for numerical optimization require (among other things) a practical procedure for generating pseudorandom nonlinear objective functions. We propose a procedure that is based on the convenient fiction that objective functions are realizations of stochastic processes. This report details the calculations necessary to implement our procedure for the case of certain stationary Gaussian processes and presents a specific implementation in the statistical programming language S-PLUS.

Trosset, Michael W.↗
NASA NTRS1999

Implementation and Validation of the Chien k-epsilon Turbulence Model in the Wind Navier-Stokes Code

The two equation k-epsilon turbulence model of Chien has been implemented in the WIND Navier-Stokes flow solver. Details of the numerical solution algorithm, initialization procedure, and stability enhancements are described. Results obtained with this version of the model are compared with those from the Chien k-epsilon model in the NPARC Navier-Stokes code and from the WIND SST model for three validation cases: the incompressible flow over a smooth flat plate, the incompressible flow over a backward facing step, and the shock-induced flow separation inside a transonic diffuser. The k-epsilon model results indicate that the WIND model functions very similarly to that in NPARC, though the WIND code appears to he slightly more accurate in the treatment of the near-wall region. Comparisons of the k-epsilon model results with those from the SST model were less definitive, as each model exhibited strengths and weaknesses for each particular case.

Yoder, Dennis A.↗
NASA NTRS1997

Airfoil Design Using a Coupled Euler and Integral Boundary Layer Method with Adjoint Based Sensitivities

The objective of this paper is to present a control theory approach for the design of airfoils in the presence of viscous compressible flows. A coupled system of the integral boundary layer and the Euler equations is solved to provide rapid flow simulations. An adjunct approach consistent with the complete coupled state equations is employed to obtain the sensitivities needed to drive a numerical optimization algorithm. Design to target pressure distribution is demonstrated on an RAE 2822 airfoil at transonic speed.

Edwards, S.↗
NASA NTRS1993

Ballistic Models for Radio Jets in Colliding Galaxies: 3C 278 (NGC 4782/4783)

We present here a general numerical simulation algorithm for ballistic radio jets. The evolution of the jets is determined by their response to mechanical forces (i.e., gravity and ram pressure). Using a previously derived collision model for the interacting pair of elliptical galaxies NGC 4782/4783, we have used this algorithm to simulate the specific two-sided jet morphology seen in the radio source associated with NGC 4782 (3C 278). This is the first time such simulations have been attempted for a galaxy pair whose relative orbit was determined independently from the jet modeling. The general morphological features of the jets in 3C 278 (e.g., position angles, lengths, curvature, and deflection angles) are well matched by the simulations, indicating that our model of the mechanical forces acting on the jets can indeed reproduce most of the basic details of the jet morphology and thereby provide an opportunity to model a combination of detailed optical, radio, and X-ray data for a wide variety of interacting radio-jet galaxies. Our models constrain the initial jet parameters (i.e., the ejection speed, direction, and starting time), the properties of the gaseous medium into which the jets are ejected, and the relative importance of gravitational deflection versus ram pressure bending in influencing the jet morphology. For 3C 278, we find that the effects of ram pressure dominate the structural evolution of the jets. In our best-fit simulation the jet ejection speed is of order 10(exp 4) km/s, to within 50%, suggesting that similar low-luminosity radio jets are also significantly subluminal. We also find that the jets are ejected within approximately 5 deg of the line of sight, the ISM in the nonjet galaxy has a much larger effect on the jet deflection than does the host galaxy, and the jet activity began just over 70 million years ago, roughly 50 million years before the pericenter passage of the two galaxies. Therefore, the onset of nuclear activity appears in 3C 278 to be related to the kinematically observed tidal shock that has been induced in the core of NGC 4782 as a result of its deeply penetrating collision with NGC 4783.

Borne, Kirk D.↗
NASA NTRS2003

Self-Consistent Model of Magnetospheric Ring Current and Electromagnetic Ion Cyclotron Waves: The May 2-7, 1998, Storm

Complete description of a self-consistent model for magnetospheric ring current interacting with electromagnetic ion cyclotron waves is presented. The model is based on the system of two kinetic equations; one equation describes the ring current ion dynamics, and another equation describes the wave evolution. The effects on ring current ions interacting with electromagnetic ion cyclotron waves, and back on waves, are considered self-consistently by solving both equations on a global magnetospheric scale under non steady-state conditions. In the paper by Khazanov et al. [2002] this self-consistent model has only been shortly outlined, and discussions of many the model related details have been omitted. For example, in present study for the first time a new algorithm for numerical finding of the resonant numbers for quasilinear wave-particle interaction is described, or it is demonstrated that in order to describe quasilinear interaction in a multi-ion thermal plasma correctly, both e and He(+) modes of electromagnetic ion cyclotron waves should be employed. The developed model is used to simulate the entire May 2-7, 1998 storm period. Trapped number fluxes of the ring current protons are calculated and presented along with their comparison with the data measured by the 3D hot plasma instrument Polar/HYDRA. Examining of the wave (MLT, L shell) distributions produced during the storm progress reveals an essential intensification of the wave emissions in about two days after main phase of storm. This result is well consistent with the earlier ground-based observations. Also the theoretical shapes and the occurrence rates for power spectral densities of electromagnetic ion cyclotron waves are studied. It is found that in about 2 days after the storm main phase on May 4, mainly non Gaussian shapes of power spectral densities are produced.

Khazanov, G. V.↗
NASA NTRS1999

The NOAA-9 Earth Radiation Budget Experiment Wide Field-of-View Data Set

The Earth Radiation Budget Experiment (ERBE) consisted of wide field-of-view (WFOV) radiometers and scanning radiometers for measuring outgoing longwave radiation and solar radiation reflected from the Earth. These instruments were carried by the dedicated Earth Radiation Budget Satellite (ERBS) and by the NOAA-9 and -10 operational spacecraft. The WFOV radiometers provided data from which instantaneous fluxes at the top of the atmosphere (TOA) are computed by use of a numerical filter algorithm. Monthly mean fluxes over a 5-degree equal angle grid are computed from the instantaneous TOA fluxes. The WFOV radiometers aboard the NOAA-9 spacecraft operated from February 1985 through December 1992, at which time a failure of the shortwave radiometer ended the usable data after nearly 8 years. This paper examines the monthly mean products from that data set.

Bush, Kathryn A.↗
NASA NTRS2004

Fourth Computational Aeroacoustics (CAA) Workshop on Benchmark Problems

This publication contains the proceedings of the Fourth Computational Aeroacoustics (CAA) Workshop on Benchmark Problems. In this workshop, as in previous workshops, the problems were devised to gauge the technological advancement of computational techniques to calculate all aspects of sound generation and propagation in air directly from the fundamental governing equations. A variety of benchmark problems have been previously solved ranging from simple geometries with idealized acoustic conditions to test the accuracy and effectiveness of computational algorithms and numerical boundary conditions; to sound radiation from a duct; to gust interaction with a cascade of airfoils; to the sound generated by a separating, turbulent viscous flow. By solving these and similar problems, workshop participants have shown the technical progress from the basic challenges to accurate CAA calculations to the solution of CAA problems of increasing complexity and difficulty. The fourth CAA workshop emphasized the application of CAA methods to the solution of realistic problems. The workshop was held at the Ohio Aerospace Institute in Cleveland, Ohio, on October 20 to 22, 2003. At that time, workshop participants presented their solutions to problems in one or more of five categories. Their solutions are presented in this proceedings along with the comparisons of their solutions to the benchmark solutions or experimental data. The five categories for the benchmark problems were as follows: Category 1:Basic Methods. The numerical computation of sound is affected by, among other issues, the choice of grid used and by the boundary conditions. Category 2:Complex Geometry. The ability to compute the sound in the presence of complex geometric surfaces is important in practical applications of CAA. Category 3:Sound Generation by Interacting With a Gust. The practical application of CAA for computing noise generated by turbomachinery involves the modeling of the noise source mechanism as a vortical gust interacting with an airfoil. Category 4:Sound Transmission and Radiation. Category 5:Sound Generation in Viscous Problems. Sound is generated under certain conditions by a viscous flow as the flow passes an object or a cavity.

Dahl, Milo D.↗