Search NASA⌕ Search

SEARCH · Search NASA

Results for “Computer Algebra”

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 793 records · Page 44
All sourcesNASA NTRSDOE OSTI
NASA NTRS1984

The cross section for double Compton scattering

Employing elementary methods in nonrelativistic quantum electrodynamics, the cross section for gamma sub 0 + e yields e + gamma + gamma is computed for arbitrary energy in the spectrum of the outgoing photons. The final result is given, differential in the energy of one of these photons, for the case where the incident photon is unpolarized and has energy E sub 0 much less than mc-squared, a polarization sum and angular integration being performed for the final-state photons. The cross section has a simple algebraic form resulting from contributions from the sum of squared direct and exchange amplitudes; interference terms from these amplitudes do not contribute to the angular-integrated cross section.

Gould, R. J.↗
NASA NTRS1986

Boundary-layer computational model for predicting the flow and heat transfer in sudden expansions

Fully developed turbulent and laminar flows through symmetric planar and axisymmetric expansions with heat transfer were modeled using a finite-difference discretization of the boundary-layer equations. By using the boundary-layer equations to model separated flow in place of the Navier-Stokes equations, computational effort was reduced permitting turbulence modelling studies to be economically carried out. For laminar flow, the reattachment length was well predicted for Reynolds numbers as low as 20 and the details of the trapped eddy were well predicted for Reynolds numbers above 200. For turbulent flows, the Boussinesq assumption was used to express the Reynolds stresses in terms of a turbulent viscosity. Near-wall algebraic turbulence models based on Prandtl's-mixing-length model and the maximum Reynolds shear stress were compared.

Lewis, J. P.↗
NASA NTRS1990

Preconditioned conjugate gradient methods for the compressible Navier-Stokes equations

The compressible Navier-Stokes equations are solved for a variety of two-dimensional inviscid and viscous problems by preconditioned conjugate gradient-like algorithms. Roe's flux difference splitting technique is used to discretize the inviscid fluxes. The viscous terms are discretized by using central differences. An algebraic turbulence model is also incorporated. The system of linear equations which arises out of the linearization of a fully implicit scheme is solved iteratively by the well known methods of GMRES (Generalized Minimum Residual technique) and Chebyschev iteration. Incomplete LU factorization and block diagonal factorization are used as preconditioners. The resulting algorithm is competitive with the best current schemes, but has wide applications in parallel computing and unstructured mesh computations.

Venkatakrishnan, V.↗
NASA NTRS1978

Numerical solution for supersonic flow near the trailing edge of a flat plate

In the present study, problems of laminar and turbulent two-dimensional flow of a viscous compressible fluid near the trailing edge of a thin flat plate are considered. The complete set of Navier-Stokes equations is solved by the finite-difference method of MacCormack (MacCormack and Baldwin, 1975). It is an explicit, predictor-corrector, time-splitting method of second order acuracy. The computational mesh employed has sufficient resolution for all the characteristic lengths suggested by theory. In the laminar case, the present results are compared with the triple deck solution of Daniels (1974). This comparison indicates that the asymptotic triple deck theory for supersonic trailing edge flow is accurate within five percent for Reynolds numbers greater than 1000. In the turbulent case, the Prandtl-Van Driest-Clauser algebraic eddy viscosity model is used. The numerical results show that the region of upstream influence is approximately of the order of the boundary layer thickness. The solutions for skin-friction, pressure and wake center-line velocity are presented.

Degani, D.↗
NASA NTRS1984

Algebraic grid generation for wing-fuselage bodies

An algebraic procedure for the generation of boundary-fitted grids about wing-fuselage configurations is presented. A wing-fuselage configuration is specified by cross sections and mathematically represented by Coons' patches. A configuration is divided into sections so that several grid blocks that either adjoin each other or partially overlap each other can be generated, and each grid has six surfaces that map into a computational cube. Grids are first determined on the six boundary surfaces and then in the interior. Grid curves that are on the surface of the configuration are derived using plane-patch intersections, and single-valued functions relating approximate arc lengths along the curves to computational coordinates define the distribution of grid points. The two-boundary technique and transfinite interpolation are used to determine the boundary surface grids that are not on the configuration, and transfinite interpolation with linear blending functions is used to determine the interior grids.

Smith, R. E.↗
NASA NTRS1986

Algebraic grid generation about wing-fuselage bodies

An algebraic procedure for the generation of boundary-fitted grids about wing-fuselage configurations is presented. A wing-fuselage configuration is specified by cross sections and mathematically represented by Coons' patches. A configuration is divided into sections so that several grid blocks that either adjoin each other or partially overlap each other can be generated. Each grid has six exterior surfaces that map into a computational cube. Grids are first determined on the six boundary surfaces and then in the interior. Grid curves that are on the surface of the configuration are derived from the intersection of planes with the Coons' patch definition. Single-valued functions relating approximate arc lengths along the grid curves to a computational coordinate define the distribution of grid points. The two-boundary technique and transfinite interpolation are used to determine the boundary surface grids that are not on the configuration, and transfinite interpolation with linear blending functions is used to determine the interior grid.

Smith, R. E.↗
NASA NTRS1994

Development of advanced Navier-Stokes solver

The objective of research was to develop and validate new computational algorithms for solving the steady and unsteady Euler and Navier-Stokes equations. The end-products are new three-dimensional Euler and Navier-Stokes codes that are faster, more reliable, more accurate, and easier to use. The three-dimensional Euler and full/thin-layer Reynolds-averaged Navier-Stokes equations for compressible/incompressible flows are solved on structured hexahedral grids. The Baldwin-Lomax algebraic turbulence model is used for closure. The space discretization is based on a cell-centered finite-volume method augmented by a variety of numerical dissipation models with optional total variation diminishing limiters. The governing equations are integrated in time by an implicit method based on lower-upper factorization and symmetric Gauss-Seidel relaxation. The algorithm is vectorized on diagonal planes of sweep using two-dimensional indices in three dimensions. Convergence rates and the robustness of the codes are enhanced by the use of an implicit full approximation storage multigrid method.

Yoon, Seokkwan↗
NASA NTRS1993

Convection equation modeling: A non-iterative direct matrix solution algorithm for use with SINDA

The determination of the boundary conditions for a component-level analysis, applying discrete finite element and finite difference modeling techniques often requires an analysis of complex coupled phenomenon that cannot be described algebraically. For example, an analysis of the temperature field of a coldplate surface with an integral fluid loop requires a solution to the parabolic heat equation and also requires the boundary conditions that describe the local fluid temperature. However, the local fluid temperature is described by a convection equation that can only be solved with the knowledge of the locally-coupled coldplate temperatures. Generally speaking, it is not computationally efficient, and sometimes, not even possible to perform a direct, coupled phenomenon analysis of the component-level and boundary condition models within a single analysis code. An alternative is to perform a disjoint analysis, but transmit the necessary information between models during the simulation to provide an indirect coupling. For this approach to be effective, the component-level model retains full detail while the boundary condition model is simplified to provide a fast, first-order prediction of the phenomenon in question. Specifically for the present study, the coldplate structure is analyzed with a discrete, numerical model (SINDA) while the fluid loop convection equation is analyzed with a discrete, analytical model (direct matrix solution). This indirect coupling allows a satisfactory prediction of the boundary condition, while not subjugating the overall computational efficiency of the component-level analysis. In the present study a discussion of the complete analysis of the derivation and direct matrix solution algorithm of the convection equation is presented. Discretization is analyzed and discussed to extend of solution accuracy, stability and computation speed. Case studies considering a pulsed and harmonic inlet disturbance to the fluid loop are analyzed to assist in the discussion of numerical dissipation and accuracy. In addition, the issues of code melding or integration with standard class solvers such as SINDA are discussed to advise the user of the potential problems to be encountered.

Schrage, Dean S.↗
NASA NTRSDate not supplied

Mars 2020 Reconstructed Aerothermal Environments and Design Margins

The Mars 2020 aeroshell's thermal protection system was nearly identical to the Mars Science Laboratory system that successfully completed its mission for the Curiosity rover's landing in 2012. It was predicted that, like Mars Science Laboratory, the Mars 2020 heatshield would experience boundary layer transition and that the thermal protection system would provide sufficient material thickness margins against the aerothermal environments, even after adding radiative heating that was not included for Mars Science Laboratory. The Mars 2020 flight instrumentation suite included heatshield sensors similar to Mars Science Laboratory and new backshell measurements. A full set of surface pressure and in-depth temperature data were collected during atmospheric entry. This paper provides an initial analysis of the flight data and focuses on the reconstructed total surface heat flux inferred from the measured temperatures. Turbulent boundary layer conditions again were observed on the Mars 2020 heatshield, and the temperatures at all eleven heatshield and six backshell thermocouple in-depth sensor locations were well within system capabilities due to favorable entry conditions and conservative uncertainties. New computational fluid dynamics results on the reconstructed entry trajectory are compared to the measured surface pressures and reconstructed total heat fluxes. The predicted heatshield surface pressures at six locations match the data qualitatively and quantitatively well, as expected. Smooth-wall laminar heating predictions prior to boundary layer transition fall above the reconstructed heating on the heatshield. After the observed boundary layer transition time, total heat flux based on algebraic turbulence model calculations generally match the heat flux trends relative to one another. The convective heat flux was predicted to constitute the majority of the total heating. On the backshell, smooth-wall laminar total heat flux predictions generally exceed the reconstructed total heating at six locations. At each of these locations, it is estimated that radiative heating provided almost all of the total heating. The Mars 2020 as-flown aerothermal environments were well below the design levels for all measurement locations due to conservative design assumptions and a stressing design trajectory.

Mars 2020↗
NASA NTRSDate not supplied

Mars 2020 Reconstructed Aerothermal Environments and Design Margins

The Mars 2020 aeroshell's thermal protection system was nearly identical to the Mars Science Laboratory system that successfully completed its mission for the Curiosity rover's landing in 2012. It was predicted that, like Mars Science Laboratory, the Mars 2020 heatshield would experience boundary layer transition and that the thermal protection system would provide sufficient material thickness margins against the aerothermal environments, even after adding radiative heating that was not included for Mars Science Laboratory. The Mars 2020 flight instrumentation suite included heatshield sensors similar to Mars Science Laboratory and new backshell measurements. A full set of surface pressure and in-depth temperature data were collected during atmospheric entry. This paper provides an initial analysis of the flight data and focuses on the reconstructed total surface heat flux inferred from the measured temperatures. Turbulent boundary layer conditions again were observed on the Mars 2020 heatshield, and the temperatures at all eleven heatshield and six backshell thermocouple in-depth sensor locations were well within system capabilities due to favorable entry conditions and conservative uncertainties. New computational fluid dynamics results on the reconstructed entry trajectory are compared to the measured surface pressures and reconstructed total heat fluxes. The predicted heatshield surface pressures at six locations match the data qualitatively and quantitatively well, as expected. Smooth-wall laminar heating predictions prior to boundary layer transition fall above the reconstructed heating on the heatshield. After the observed boundary layer transition time, total heat flux based on algebraic turbulence model calculations generally match the heat flux trends relative to one another. The convective heat flux was predicted to constitute the majority of the total heating. On the backshell, smooth-wall laminar total heat flux predictions generally exceed the reconstructed total heating at six locations. At each of these locations, it is estimated that radiative heating provided almost all of the total heating. The Mars 2020 as-flown aerothermal environments were well below the design levels for all measurement locations due to conservative design assumptions and a stressing design trajectory.

Mars 2020↗
NASA NTRS1976

New convergence criteria for the vortex-lattice models of the leading-edge separation

The convergence criterion for the vortex-lattice technique which deals with delta wings exhibiting significant leading-edge separation was studied. It was shown that one can predict pressure distributions without irregularities which agree fairly well with experimental data (which show some irregularities of their own) by replacing the system of discrete vortex lines with a single concentrated core. This core has a circulation equal to the algebraic sum of the circulations around the discrete lines and is located at the centroid of these lines. Moreover, there is a requirement that the position and strength of the core must converge as the number of elements increases. Because the calculation of the position and strength of the core is much less involved than the calculation of the loads, this approach has the additional desirable feature of requiring less computational time.

Kandil, O. A.↗
NASA NTRS2018

Development and Documentation of kL-Based Linear, Nonlinear, and Full Reynolds Stress Turbulence Models

The development and implementation of kL-based Reynolds Average Navier-Stokes (RANS) turbulence models are reported herein. The kL is based on Abdol-Hamid's closure and Menter's modi cation to Rotta's two-equation model. Rotta shows that a reliable transport equation can be formed from the turbulent length scale L, and the turbulent kinetic energy k. Rotta's kL equation is well suited for term-by-term modeling and displays useful features compared to other scale formulation. One of the important di erences is the inclusion of higher order velocity derivatives in the source terms of the scale equation. This can enhance the ability of RANS solvers to simulate unsteady ows in URANS mode. The present report documents the formulation of three model levels of turbulence models as implemented in the CFD code FUN3D. Methodology and calibration examples are shown in detail. The levels are the linear k-kL and the two-equation algebraic Reynolds stress model (ARSM) as well as the full Reynolds Stress Model (RSM). Attached, separated and corner ow cases are documented and compared with experimental, theoretical and other turbulence model data. The results show generally very good comparisons with canonical and experimental data. The results from this formulation are similar or better than results using the SST two- equation turbulence model. ARSM shows great promise with similar level of computational resources as general two equation turbulence models.

Abdol-Hamid, Khaled S.↗
NASA NTRS1978

An evaluation of several compressible turbulent boundary-layer models Effects of pressure gradient and Reynolds number

Computations, employing several turbulence models, are compared with a series of attached supersonic turbulent boundary-layer experiments over an extensive range of Reynolds numbers (11.7 x 10 to the 6th to 314 x 10 to the 6th). These experiments included measurements of surface pressure and skin friction for adverse pressure gradients ranging in magnitude from those of previous investigations to an order of magnitude greater. The turbulence models evaluated include algebraic and two-equation eddy-viscosity models and two full Reynolds stress models. In general all the models tested performed well independently of the magnitude of the pressure gradient or Reynolds number and could predict the measured skin friction for most cases with sufficient accuracy for engineering purposes.

Horstman, C. C.↗
NASA NTRS2001

Development of an Unstructured Mesh Code for Flows About Complete Vehicles

This report describes the research work undertaken at the Massachusetts Institute of Technology, under NASA Research Grant NAG4-157. The aim of this research is to identify effective algorithms and methodologies for the efficient and routine solution of flow simulations about complete vehicle configurations. For over ten years we have received support from NASA to develop unstructured mesh methods for Computational Fluid Dynamics. As a result of this effort a methodology based on the use of unstructured adapted meshes of tetrahedra and finite volume flow solvers has been developed. A number of gridding algorithms, flow solvers, and adaptive strategies have been proposed. The most successful algorithms developed from the basis of the unstructured mesh system FELISA. The FELISA system has been extensively for the analysis of transonic and hypersonic flows about complete vehicle configurations. The system is highly automatic and allows for the routine aerodynamic analysis of complex configurations starting from CAD data. The code has been parallelized and utilizes efficient solution algorithms. For hypersonic flows, a version of the code which incorporates real gas effects, has been produced. The FELISA system is also a component of the STARS aeroservoelastic system developed at NASA Dryden. One of the latest developments before the start of this grant was to extend the system to include viscous effects. This required the development of viscous generators, capable of generating the anisotropic grids required to represent boundary layers, and viscous flow solvers. We show some sample hypersonic viscous computations using the developed viscous generators and solvers. Although this initial results were encouraging it became apparent that in order to develop a fully functional capability for viscous flows, several advances in solution accuracy, robustness and efficiency were required. In this grant we set out to investigate some novel methodologies that could lead to the required improvements. In particular we focused on two fronts: (1) finite element methods and (2) iterative algebraic multigrid solution techniques.

Peraire, Jaime↗
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.↗
NASA NTRS1996

Application of Hybrid Optimization-Expert System for Optimal Power Management on Board Space Power Station

The space power system has two sources of energy: photo-voltaic blankets and batteries. The optimal power management problem on-board has two broad operations: off-line power scheduling to determine the load allocation schedule of the next several hours based on the forecast of load and solar power availability. The nature of this study puts less emphasis on speed requirement for computation and more importance on the optimality of the solution. The second category problem, on-line power rescheduling, is needed in the event of occurrence of a contingency to optimally reschedule the loads to minimize the 'unused' or 'wasted' energy while keeping the priority on certain type of load and minimum disturbance of the original optimal schedule determined in the first-stage off-line study. The computational performance of the on-line 'rescheduler' is an important criterion and plays a critical role in the selection of the appropriate tool. The Howard University Center for Energy Systems and Control has developed a hybrid optimization-expert systems based power management program. The pre-scheduler has been developed using a non-linear multi-objective optimization technique called the Outer Approximation method and implemented using the General Algebraic Modeling System (GAMS). The optimization model has the capability of dealing with multiple conflicting objectives viz. maximizing energy utilization, minimizing the variation of load over a day, etc. and incorporates several complex interaction between the loads in a space system. The rescheduling is performed using an expert system developed in PROLOG which utilizes a rule-base for reallocation of the loads in an emergency condition viz. shortage of power due to solar array failure, increase of base load, addition of new activity, repetition of old activity etc. Both the modules handle decision making on battery charging and discharging and allocation of loads over a time-horizon of a day divided into intervals of 10 minutes. The models have been extensively tested using a case study for the Space Station Freedom and the results for the case study will be presented. Several future enhancements of the pre-scheduler and the 'rescheduler' have been outlined which include graphic analyzer for the on-line module, incorporating probabilistic considerations, including spatial location of the loads and the connectivity using a direct current (DC) load flow model.

Momoh, James↗
NASA NTRS1995

Global Dynamic Modeling of Space-Geodetic Data

The proposal had outlined a year for program conversion, a year for testing and debugging, and two years for numerical experiments. We kept to that schedule. In first (partial) year, author designed a finite element for isostatic thin-shell deformation on a sphere, derived all of its algebraic and stiffness properties, and embedded it in a new finite element code which derives its basic solution strategy (and some critical subroutines) from earlier flat-Earth codes. Also designed and programmed a new fault element to represent faults along plate boundaries. Wrote a preliminary version of a spherical graphics program for the display of output. Tested this new code for accuracy on individual model plates. Made estimates of the computer-time/cost efficiency of the code for whole-earth grids, which were reasonable. Finally, converted an interactive graphical grid-designer program from Cartesian to spherical geometry to permit the beginning of serious modeling. For reasons of cost efficiency, models are isostatic, and do not consider the local effects of unsupported loads or bending stresses. The requirements are: (1) ability to represent rigid rotation on a sphere; (2) ability to represent a spatially uniform strain-rate tensor in the limit of small elements; and (3) continuity of velocity across all element boundaries. Author designed a 3-node triangle shell element which has two different sets of basis functions to represent (vector) velocity and all other (scalar) variables. Such elements can be shown to converge to the formulas for plane triangles in the limit of small size, but can also applied to cover any area smaller than a hemisphere. The difficult volume integrals involved in computing the stiffness of such elements are performed numerically using 7 Gauss integration points on the surface of the sphere, beneath each of which a vertical integral is performed using about 100 points.

Bird, Peter↗
NASA NTRS1988

Transonic aeroelastic analysis of launch vehicle configurations

A numerical study of the aeroelastic stability of typical launch vehicle configurations in transonic flight is performed. Recent computational fluid dynamics techniques are used to simulate the transonic aerodynamic flow fields, as opposed to relying on experimental data for the unsteady aerodynamic pressures. The flow solver is coupled to an appropriate structural representation of the vehicle. The aerodynamic formulation is based on the thin layer approximation to the Reynolds-Averaged Navier-Stokes equations, where the account for turbulent mixing is done by the two-layer Baldwin and Lomax algebraic eddy viscosity model. The structural-dynamic equations are developed considering free-free flexural vibration of an elongated beam with variable properties and are cast in modal form. Aeroelastic analyses are performed by integrating simultaneously in the two sets of equations. By tracing the growth or decay of a perturbed oscillation, the aeroelastic stability of a given constant configuration can be ascertained. The method is described in detail, and results that indicate its application are presented. Applications include some validation cases for the algorithm developed, as well as the study of configurations known to have presented flutter programs in the past.

Filgueirasdeazevedo, Joao Luiz↗