SEARCH · Search NASA

Results for “numerical algorithm computing”

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 1,351 records · Page 75

All sources NASA NTRS DOE OSTI

NASA NTRS2001

Weather Prediction Improvement Using Advanced Satellite Technology

We discuss in this paper some of the problems that exist today in the fall utilization of satellite data to improve weather forecasts and we propose specific recommendations to solve them. This discussion can be viewed as an aspect of the general debate on how best to organize the transition from research to operational satellites and how to evaluate the impact of a research instrument on numerical weather predictions. A method for providing this transition is offered by the National Polar-Orbiting Operational Environmental Satellite System (NPOESS) Preparatory Project (NPP). This mission will bridge the time between the present NOAA and Department of Defense (DOD) polar orbiting missions and the initiation of the converged NPOESS series and will evaluate some of the Earth Observing System (EOS) instruments as appropriate for operational missions. Thus, this mission can be viewed as an effort to meet the operational requirements of NOAA and DOD and the research requirements of NASA. More generally, however, it can be said that the process of going from the conception of new, more advanced instruments to their operational implementation and full utilization by the weather forecast communities is not optimal. Instruments developed for research purposes may have insufficient funding to explore their potential operational capabilities. Furthermore, instrument development programs designed for operational satellites typically have insufficient funding for assimilation algorithms needed to transform the satellite observations into data that can be used by sophisticated global weather forecast models. As a result, years often go by before satellite data are efficiently used for operational forecasts. NASA and NOAA each have unique expertise in the design of satellite instruments, their use for basic and applied research and their utilization in weather and climate research. At a time of limited resources, the two agencies must combine their efforts to work toward common goals of full utilization of satellite data. This is a challenge that requires the assimilation of myriad new data into increasingly sophisticated numerical forecast models that run on increasingly sophisticated computer systems. In section II, we briefly outline the impact of satellite data on the quality of the National Centers for Environmental Prediction (NCEP) forecasts. In section III, we describe the present status of the utilization of satellite data in NCEP models and the challenges that lie ahead. In section IV, we propose solutions whose goals are summarized in section V.

Einaudi, Franco↗

NASA NTRS1980

CAS2D: FORTRAN program for nonrotating blade-to-blade, steady, potential transonic cascade flows

An exact, full-potential-equation (FPE) model for the steady, irrotational, homentropic and homoenergetic flow of a compressible, homocompositional, inviscid fluid through two dimensional planar cascades of airfoils was derived, together with its appropriate boundary conditions. A computer program, CAS2D, was developed that numerically solves an artificially time-dependent form of the actual FPE. The governing equation was discretized by using type-dependent, rotated finite differencing and the finite area technique. The flow field was discretized by providing a boundary-fitted, nonuniform computational mesh. The mesh was generated by using a sequence of conforming mapping, nonorthogonal coordinate stretching, and local, isoparametric, bilinear mapping functions. The discretized form of the FPE was solved iteratively by using successive line overrelaxation. The possible isentropic shocks were correctly captured by adding explicitly an artificial viscosity in a conservative form. In addition, a three-level consecutive, mesh refinement feature makes CAS2D a reliable and fast algorithm for the analysis of transonic, two dimensional cascade flows.

Dulikravich, D. S.↗

NASA NTRS1989

Numerical Simulation of Flow Through an Artificial Heart

A solution procedure was developed that solves the unsteady, incompressible Navier-Stokes equations, and was used to numerically simulate viscous incompressible flow through a model of the Pennsylvania State artificial heart. The solution algorithm is based on the artificial compressibility method, and uses flux-difference splitting to upwind the convective terms; a line-relaxation scheme is used to solve the equations. The time-accuracy of the method is obtained by iteratively solving the equations at each physical time step. The artificial heart geometry involves a piston-type action with a moving solid wall. A single H-grid is fit inside the heart chamber. The grid is continuously compressed and expanded with a constant number of grid points to accommodate the moving piston. The computational domain ends at the valve openings where nonreflective boundary conditions based on the method of characteristics are applied. Although a number of simplifing assumptions were made regarding the geometry, the computational results agreed reasonably well with an experimental picture. The computer time requirements for this flow simulation, however, are quite extensive. Computational study of this type of geometry would benefit greatly from improvements in computer hardware speed and algorithm efficiency enhancements.

Rogers, Stuart E.↗

NASA NTRS1990

Multi-blocking strategies for the INS3D incompressible Navier-Stokes code

With the continuing development of bigger and faster supercomputers, computational fluid dynamics (CFD) has become a useful tool for real-world engineering design and analysis. However, the number of grid points necessary to resolve realistic flow fields numerically can easily exceed the memory capacity of available computers. In addition, geometric shapes of flow fields, such as those in the Space Shuttle Main Engine (SSME) power head, may be impossible to fill with continuous grids upon which to obtain numerical solutions to the equations of fluid motion. The solution to this dilemma is simply to decompose the computational domain into subblocks of manageable size. Computer codes that are single-block by construction can be modified to handle multiple blocks, but ad-hoc changes in the FORTRAN have to be made for each geometry treated. For engineering design and analysis, what is needed is generalization so that the blocking arrangement can be specified by the user. INS3D is a computer program for the solution of steady, incompressible flow problems. It is used frequently to solve engineering problems in the CFD Branch at Marshall Space Flight Center. INS3D uses an implicit solution algorithm and the concept of artificial compressibility to provide the necessary coupling between the pressure field and the velocity field. The development of generalized multi-block capability in INS3D is described.

Gatlin, Boyd↗

NASA NTRS1993

Three-dimensional solidification and melting using magnetic field control

A new two-fluid mathematical model for fully three dimensional steady solidification under the influence of an arbitrary acceleration vector and with or without an arbitrary externally applied steady magnetic field have been formulated and integrated numerically. The model includes Joule heating and allows for separate temperature dependent physical properties within the melt and the solid. Latent heat of phase change during melting/solidification was incorporated using an enthalpy method. Mushy region was automatically captured by varying viscosity orders of magnitude between liquidus and solidus temperature. Computational results were obtained for silicon melt solidification in a parallelepiped container cooled from above and from a side. The results confirm that the magnetic field has a profound influence on the solidifying melt flow field thus changing convective heat transfer through the boundaries and the amount and shape of the solid accrued. This suggests that development of a quick-response algorithm for active control of three dimensional solidification is feasible since it would require low strength magnetic fields.

Dulikravich, George S.↗

NASA NTRS2005

Best-Fit Conic Approximation of Spacecraft Trajectory

A computer program calculates a best conic fit of a given spacecraft trajectory. Spacecraft trajectories are often propagated as conics onboard. The conic-section parameters as a result of the best-conic-fit are uplinked to computers aboard the spacecraft for use in updating predictions of the spacecraft trajectory for operational purposes. In the initial application for which this program was written, there is a requirement to fit a single conic section (necessitated by onboard memory constraints) accurate within 200 microradians to a sequence of positions measured over a 4.7-hour interval. The present program supplants a prior one that could not cover the interval with fewer than four successive conic sections. The present program is based on formulating the best-fit conic problem as a parameter-optimization problem and solving the problem numerically, on the ground, by use of a modified steepest-descent algorithm. For the purpose of this algorithm, optimization is defined as minimization of the maximum directional propagation error across the fit interval. In the specific initial application, the program generates a single 4.7-hour conic, the directional propagation of which is accurate to within 34 microradians easily exceeding the mission constraints by a wide margin.

Singh, Gurkipal↗

NASA NTRSDate not supplied

Investigating the Feasibility and Stability for Modeling Broadband Acoustic Wave Scattering using a Time-Domain Boundary Integral Equation

Reducing aircraft noise is a major objective in the field of computational aeroacoustics. Acoustic liners are an effective tool for reducing aircraft noise and are characterized by a frequency-dependent impedance value. Converted into the time-domain using Fourier transforms, an impedance boundary condition can be used to simulate the acoustic wave scattering by geometric bodies treated with acoustic liners. A Broadband Impedance Model is discussed in which the liner impedance is specified along a wide range of frequencies. An impedance boundary condition is derived and coupled with a time-domain boundary integral equation to model acoustic scattering by a flat plate. It is assumed the flat plate has surfaces treated with acoustic liners. The stability of the numerical algorithm is assessed using eigenvalue analysis.

mathmatics↗

NASA NTRS1983

Numerical computation of transonic flow governed by the full-potential equation

Numerical solution techniques for solving transonic flow fields governed by the full potential equation are discussed. In a general sense relaxation schemes suitable for the numerical solution of elliptic partial differential equations are presented and discussed with emphasis on transonic flow applications. The presentation can be divided into two general categories: An introductory treatment of the basic concepts associated with the numerical solution of elliptic partial differential equations and a more advanced treatment of current procedures used to solve the full potential equation for transonic flow fields. The introductory material is presented for completeness and includes a brief introduction (Chapter 1), governing equations (Chapter 2), classical relaxation schemes (Chapter 3), and early concepts regarding transonic full potential equation algorithms (Chapter 4).

Holst, T. L.↗

NASA NTRS1992

Multigrid calculation of internal flows in complex geometries

The development, validation, and application of a general purpose multigrid solution algorithm and computer program for the computation of elliptic flows in complex geometries is presented. This computer program combines several desirable features including a curvilinear coordinate system, collocated arrangement of the variables, and Full Multi-Grid/Full Approximation Scheme (FMG/FAS). Provisions are made for the inclusion of embedded obstacles and baffles inside the flow domain. The momentum and continuity equations are solved in a decoupled manner and a pressure corrective equation is used to update the pressures such that the fluxes at the cell faces satisfy local mass continuity. Despite the computational overhead required in the restriction and prolongation phases of the multigrid cycling, the superior convergence results in reduced overall CPU time. The numerical scheme and selected results of several validation flows are presented. Finally, the procedure is applied to study the flowfield in a side-inlet dump combustor and twin jet impingement from a simulated aircraft fuselage.

Smith, K. M.↗

NASA NTRS1995

Analysis of composite ablators using massively parallel computation

In this work, the feasibility of using massively parallel computation to study the response of ablative materials is investigated. Explicit and implicit finite difference methods are used on a massively parallel computer, the Thinking Machines CM-5. The governing equations are a set of nonlinear partial differential equations. The governing equations are developed for three sample problems: (1) transpiration cooling, (2) ablative composite plate, and (3) restrained thermal growth testing. The transpiration cooling problem is solved using a solution scheme based solely on the explicit finite difference method. The results are compared with available analytical steady-state through-thickness temperature and pressure distributions and good agreement between the numerical and analytical solutions is found. It is also found that a solution scheme based on the explicit finite difference method has the following advantages: incorporates complex physics easily, results in a simple algorithm, and is easily parallelizable. However, a solution scheme of this kind needs very small time steps to maintain stability. A solution scheme based on the implicit finite difference method has the advantage that it does not require very small times steps to maintain stability. However, this kind of solution scheme has the disadvantages that complex physics cannot be easily incorporated into the algorithm and that the solution scheme is difficult to parallelize. A hybrid solution scheme is then developed to combine the strengths of the explicit and implicit finite difference methods and minimize their weaknesses. This is achieved by identifying the critical time scale associated with the governing equations and applying the appropriate finite difference method according to this critical time scale. The hybrid solution scheme is then applied to the ablative composite plate and restrained thermal growth problems. The gas storage term is included in the explicit pressure calculation of both problems. Results from ablative composite plate problems are compared with previous numerical results which did not include the gas storage term. It is found that the through-thickness temperature distribution is not affected much by the gas storage term. However, the through-thickness pressure and stress distributions, and the extent of chemical reactions are different from the previous numerical results. Two types of chemical reaction models are used in the restrained thermal growth testing problem: (1) pressure-independent Arrhenius type rate equations and (2) pressure-dependent Arrhenius type rate equations. The numerical results are compared to experimental results and the pressure-dependent model is able to capture the trend better than the pressure-independent one. Finally, a performance study is done on the hybrid algorithm using the ablative composite plate problem. It is found that there is a good speedup of performance on the CM-5. For 32 CPU's, the speedup of performance is 20. The efficiency of the algorithm is found to be a function of the size and execution time of a given problem and the effective parallelization of the algorithm. It also seems that there is an optimum number of CPU's to use for a given problem.

Shia, David↗

NASA NTRS1996

Parallel Newton-Krylov-Schwarz algorithms for the transonic full potential equation

We study parallel two-level overlapping Schwarz algorithms for solving nonlinear finite element problems, in particular, for the full potential equation of aerodynamics discretized in two dimensions with bilinear elements. The overall algorithm, Newton-Krylov-Schwarz (NKS), employs an inexact finite-difference Newton method and a Krylov space iterative method, with a two-level overlapping Schwarz method as a preconditioner. We demonstrate that NKS, combined with a density upwinding continuation strategy for problems with weak shocks, is robust and, economical for this class of mixed elliptic-hyperbolic nonlinear partial differential equations, with proper specification of several parameters. We study upwinding parameters, inner convergence tolerance, coarse grid density, subdomain overlap, and the level of fill-in in the incomplete factorization, and report their effect on numerical convergence rate, overall execution time, and parallel efficiency on a distributed-memory parallel computer.

Cai, Xiao-Chuan↗

NASA NTRS1994

Parallel Computational Fluid Dynamics: Current Status and Future Requirements

One or the key objectives of the Applied Research Branch in the Numerical Aerodynamic Simulation (NAS) Systems Division at NASA Allies Research Center is the accelerated introduction of highly parallel machines into a full operational environment. In this report we discuss the performance results obtained from the implementation of some computational fluid dynamics (CFD) applications on the Connection Machine CM-2 and the Intel iPSC/860. We summarize some of the experiences made so far with the parallel testbed machines at the NAS Applied Research Branch. Then we discuss the long term computational requirements for accomplishing some of the grand challenge problems in computational aerosciences. We argue that only massively parallel machines will be able to meet these grand challenge requirements, and we outline the computer science and algorithm research challenges ahead.

Simon, Horst D.↗

NASA NTRS1989

Modeling of high speed chemically reacting flow-fields

The SPARK3D and SPARK3D-PNS computer programs were developed to model 3-D supersonic, chemically reacting flow-fields. The SPARK3D code is a full Navier-Stokes solver, and is suitable for use in scramjet combustors and other regions where recirculation may be present. The SPARK3D-PNS is a parabolized Navier-Stokes solver and provides an efficient means of calculating steady-state combustor far-fields and nozzles. Each code has a generalized chemistry package, making modeling of any chemically reacting flow possible. Research activities by the Langley group range from addressing fundamental theoretical issues to simulating problems of practical importance. Algorithmic development includes work on higher order and upwind spatial difference schemes. Direct numerical simulations employ these algorithms to address the fundamental issues of flow stability and transition, and the chemical reaction of supersonic mixing layers and jets. It is believed that this work will lend greater insight into phenomenological model development for simulating supersonic chemically reacting flows in practical combustors. Currently, the SPARK3D and SPARK3D-PNS codes are used to study problems of engineering interest, including various injector designs and 3-D combustor-nozzle configurations. Examples, which demonstrate the capabilities of each code are presented.

Drummond, J. P.↗

NASA NTRS1997

Nonlinear Behavior of a Typical Airfoil Section with Control Surface Freeplay: A Numerical and Experimental Study

A three degree-of-freedom aeroelastic typical section with control surface freeplay is modeled theoretically as a system of piecewise linear state-space models. The system response is determined by time marching of the governing equations using a standard Runge-Kutta algorithm in conjunction with Henon's method for integrating a system of equations to a prescribed surface of phase space section. Henon's method is used to locate the "switching points" accurately and efficiently as the system moves from one linear region into another. An experimental model which closely approximates the three degree-of-freedom, typical section in two-dimensional, incompressible flow has been created to validate the theoretical model. Consideration is given to modeling realistically the structural damping present in the experimental system. The effect of the freeplay on the system response is examined numerically and experimentally. The development of the state-space model offers a low-order, computationally efficient means of modeling fully the freeplay nonlinearity and may offer advantages in future research which will investigate the effects of freeplay on the control of flutter in the typical section.

Conner, M. D.↗

NASA NTRS1998

Algorithms for Efficient Computation of Transfer Functions for Large Order Flexible Systems

An efficient and robust computational scheme is given for the calculation of the frequency response function of a large order, flexible system implemented with a linear, time invariant control system. Advantage is taken of the highly structured sparsity of the system matrix of the plant based on a model of the structure using normal mode coordinates. The computational time per frequency point of the new computational scheme is a linear function of system size, a significant improvement over traditional, still-matrix techniques whose computational times per frequency point range from quadratic to cubic functions of system size. This permits the practical frequency domain analysis of systems of much larger order than by traditional, full-matrix techniques. Formulations are given for both open- and closed-loop systems. Numerical examples are presented showing the advantages of the present formulation over traditional approaches, both in speed and in accuracy. Using a model with 703 structural modes, the present method was up to two orders of magnitude faster than a traditional method. The present method generally showed good to excellent accuracy throughout the range of test frequencies, while traditional methods gave adequate accuracy for lower frequencies, but generally deteriorated in performance at higher frequencies with worst case errors being many orders of magnitude times the correct values.

Maghami, Peiman G.↗

NASA NTRS1994

Numerical Wave Propagation 211 Based on Wave Primitives

Compact higher order finite difference equations are applied to a sequence of problems in wave propagation and aeroacoustics. Systems of PDE's are reduced to a sequence of simple wave primitives using a local eigenvector decomposition. The wave primitives are first order PDE's in two independent variable and allow natural boundary conditions to be imposed for both single- and multidimensional problems. The method uses a "discrete dispersion relation" approach to obtain high order approximations to the wave primitives on a 3 spatial point / 2 time level computational molecule. The scheme is fourth order accurate for the class of system with constant coefficients, e.g., those that support exponential solutions. Weakly non-linear PDE's are solved in a similar manner using a variant of the "method of frozen coefficients." Experience with the new algorithm for linear, non-linear, and multi-dimensional test problems will be described.

Davis, Sanford S.↗

NASA NTRS1973

Free vibration analysis of spinning structural systems.

This article presents an efficient digital computer procedure, along with the complete listing of the associated computer program, which may be conveniently utilized for the accurate solution of a wide range of practical eigenvalue problems. Important applications of the present work are envisaged in the natural frequency analysis of spinning structures discretized by the finite element technique, and in the determination of transfer functions associated with the dynamic blocks of control systems of spacecraft utilizing gas jets or reaction wheels for attitude control, as well as of spin-stabilized and dual-spin-stabilized satellites. The validity of the Sturm sequence property is first established for the related matrix formulation involving Hermitian and real symmetric, positive-definite matrices, both being usually of highly banded configuration. A numerically stable algorithm based on the Sturm sequence method is then developed which fully exploits the banded form of the associated matrices.

Gupta, K. K.↗

NASA NTRS1991

Time-domain finite elements in optimal control with application to launch-vehicle guidance

A time-domain finite element method is developed for optimal control problems. The theory derived is general enough to handle a large class of problems including optimal control problems that are continuous in the states and controls, problems with discontinuities in the states and/or system equations, problems with control inequality constraints, problems with state inequality constraints, or problems involving any combination of the above. The theory is developed in such a way that no numerical quadrature is necessary regardless of the degree of nonlinearity in the equations. Also, the same shape functions may be employed for every problem because all strong boundary conditions are transformed into natural or weak boundary conditions. In addition, the resulting nonlinear algebraic equations are very sparse. Use of sparse matrix solvers allows for the rapid and accurate solution of very difficult optimization problems. The formulation is applied to launch-vehicle trajectory optimization problems, and results show that real-time optimal guidance is realizable with this method. Finally, a general problem solving environment is created for solving a large class of optimal control problems. The algorithm uses both FORTRAN and a symbolic computation program to solve problems with a minimum of user interaction. The use of symbolic computation eliminates the need for user-written subroutines which greatly reduces the setup time for solving problems.

Bless, Robert R.↗