Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical computation”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Object oriented development of engineering software using CLIPS

Engineering applications involve numeric complexity and manipulations of a large amount of data. Traditionally, numeric computation has been the concern in developing an engineering software. As engineering application software became larger and more complex, management of resources such as data, rather than the numeric complexity, has become the major software design problem. Object oriented design and implementation methodologies can improve the reliability, flexibility, and maintainability of the resulting software; however, some tasks are better solved with the traditional procedural paradigm. The C Language Integrated Production System (CLIPS), with deffunction and defgeneric constructs, supports the procedural paradigm. The natural blending of object oriented and procedural paradigms has been cited as the reason for the popularity of the C++ language. The CLIPS Object Oriented Language's (COOL) object oriented features are more versatile than C++'s. A software design methodology based on object oriented and procedural approaches appropriate for engineering software, and to be implemented in CLIPS was outlined. A method for sensor placement for Space Station Freedom is being implemented in COOL as a sample problem.

Yoon, C. John↗
NASA NTRS1993

Computation of the sound generated by isotropic turbulence

The acoustic radiation from isotropic turbulence is computed numerically. A hybrid direct numerical simulation approach which combines direct numerical simulation (DNS) of the turbulent flow with the Lighthill acoustic analogy is utilized. It is demonstrated that the hybrid DNS method is a feasible approach to the computation of sound generated by turbulent flows. The acoustic efficiency in the simulation of isotropic turbulence appears to be substantially less than that in subsonic jet experiments. The dominant frequency of the computed acoustic pressure is found to be somewhat larger than the dominant frequency of the energy-containing scales of motion. The acoustic power in the simulations is proportional to epsilon (M(sub t))(exp 5) where epsilon is the turbulent dissipation rate and M(sub t) is the turbulent Mach number. This is in agreement with the analytical result of Proudman (1952), but the constant of proportionality is smaller than the analytical result. Two different methods of computing the acoustic power from the DNS data bases yielded consistent results.

Sarkar, S.↗
NASA NTRS1980

Flow oscillations of spike-tipped bodies

The oscillating flow field around a spike-tipped body has been simulated by the numerical solution of the 3-D compressible Navier-Stokes equations at a Mach number of 3.0 and a nominal Reynolds number of 7.87 x 10 to the 6th/m. Computations were performed using a vectorized 3-D Navier-Stokes program on the STAR 100 computer. Numerical solutions confirmed the experimental result that the self-sustained oscillation occurs within a limited range of the protruded spike length to shoulder height ratio. The numerical result predicted correctly the discrete frequency range as well as the rms pressure intensity. The detailed flow structure is also presented and discussed.

Shang, J. S.↗
NASA NTRS1992

The Osher scheme for non-equilibrium reacting flows

An extension of the Osher upwind scheme to nonequilibrium reacting flows is presented. Owing to the presence of source terms, the Riemann problem is no longer self-similar and therefore its approximate solution becomes tedious. With simplicity in mind, a linearized approach which avoids an iterative solution is used to define the intermediate states and sonic points. The source terms are treated explicitly. Numerical computations are presented to demonstrate the feasibility, efficiency and accuracy of the proposed method. The test problems include a ZND (Zeldovich-Neumann-Doring) detonation problem for which spurious numerical solutions which propagate at mesh speed have been observed on coarse grids. With the present method, a change of limiter causes the solution to change from the physically correct CJ detonation solution to the spurious weak detonation solution.

Suresh, Ambady↗
NASA NTRS1998

Performance Analysis and Portability of the PLUM Load Balancing System

The ability to dynamically adapt an unstructured mesh is a powerful tool for solving computational problems with evolving physical features; however, an efficient parallel implementation is rather difficult. To address this problem, we have developed PLUM, an automatic portable framework for performing adaptive numerical computations in a message-passing environment. PLUM requires that all data be globally redistributed after each mesh adaption to achieve load balance. We present an algorithm for minimizing this remapping overhead by guaranteeing an optimal processor reassignment. We also show that the data redistribution cost can be significantly reduced by applying our heuristic processor reassignment algorithm to the default mapping of the parallel partitioner. Portability is examined by comparing performance on a SP2, an Origin2000, and a T3E. Results show that PLUM can be successfully ported to different platforms without any code modifications.

Oliker, Leonid↗
NASA NTRS1973

Thermal history and evolution of the moon.

In this work, theoretical lunar temperature models are computed taking into account different initial conditions to represent possible accretion models and various abundances of heat sources to correspond to different compositions. Differentiation and convection are simulated in the numerical computational scheme. Models of the thermal evolution of the moon that fit the chronology of igneous activity on the lunar surface, the stress history of the lunar lithosphere implied by the presence of mascons, and the surface concentrations of radioactive elements, involve extensive differentiation early in lunar history. This differentiation may be the result of rapid accretion and large-scale melting or of primary chemical layering during accretion. Differences in present-day temperatures for these two possibilities are significant only in the inner 1000 km of the moon and are not resolvable with presently available data.

Toksoz, M. N.↗
NASA NTRS1989

Efficient numerical method for computation of thermohydrodynamics of laminar lubricating films

The purpose of this paper is to describe an accurate, yet economical, method for computing temperature effects in laminar lubricating films in two dimensions. The procedure presented here is a sequel to one presented in Leeds in 1986 that was carried out for the one-dimensional case. Because of the marked dependence of lubricant viscosity on temperature, the effect of viscosity variation both across and along a lubricating film can dwarf other deviations from ideal constant-property lubrication. In practice, a thermohydrodynamics program will involve simultaneous solution of the film lubrication problem, together with heat conduction in a solid, complex structure. The extent of computation required makes economy in numerical processing of utmost importance. In pursuit of such economy, we here use techniques similar to those for Gaussian quadrature. We show that, for many purposes, the use of just two properly positioned temperatures (Lobatto points) characterizes well the transverse temperature distribution.

Elrod, Harold G.↗
NASA NTRS1989

Prediction of turbine rotor-stator interaction using Navier-Stokes methods

Flows in turbomachinery are generally complex and do not easily lend themselves to numerical computation. The flows are three-dimensional and inherently unsteady. Complicated blade geometries and flow phenomena such as separation and periodic transition from laminar to turbulent flow add to the numerical complexity. Nevertheless, the accurate numerical analysis of such flows is a problem of considerable interest and practical importance to the turbomachinery community. Much of the early work in turbomachinery flow prediction focussed on airfoil cascades. While such analyses of flows in isolated airfoil rows have helped improve understanding of the flow phenomena and have gained widespread acceptance in the industrial community as a design tool, they do not yield any information regarding the unsteady effects arising out of rotor-stator aerodynamic interaction. These interaction effects become increasingly important as the distance between successive stator and rotor rows is decreased. Thus, the need exists for analytical tools that treat the rotor and stator airfoils as a system and provide information regarding the magnitude and the impact of the unsteady effects. The focus a three-dimensional, time-accurate, thin-layer Navier-Stokes code that was recently developed to study rotor-stator interaction problems. A system of patched and overlaid grids that move relative to each other is used to discretize the flow field and the governing equations are integrated using a third-order upwind scheme set in an iterative, implicit framework. The code was used to simulate subsonic flow through an axial turbine configuration for which considerable experimental data exists. Grid refinement studies were also conducted as part of the code validation process. The current status of the research, along with planned future directions, are also discussed.

Madavan, Nateri K.↗
NASA NTRS2003

Study of Unsteady Flows with Concave Wall Effect

This paper presents computational fluid dynamic studies of the inlet turbulence and wall curvature effects on the flow steadiness at near wall surface locations in boundary layer flows. The time-stepping RANS numerical solver of the NASA Glenn-HT RANS code and a one-equation turbulence model, with a uniform inlet turbulence modeling level of the order of 10 percent of molecular viscosity, were used to perform the numerical computations. The approach was first calibrated for its predictabilities of friction factor, velocity, and temperature at near surface locations within a transitional boundary layer over concave wall. The approach was then used to predict the velocity and friction factor variations in a boundary layer recovering from concave curvature. As time iteration proceeded in the computations, the computed friction factors converged to their values from existing experiments. The computed friction factors, velocity, and static temperatures at near wall surface locations oscillated periodically in terms of time iteration steps and physical locations along the span-wise direction. At the upstream stations, the relationship among the normal and tangential velocities showed vortices effects on the velocity variations. Coherent vortices effect on the velocity components broke down at downstream stations. The computations also predicted the vortices effects on the velocity variations within a boundary layer flow developed along a concave wall surface with a downstream recovery flat wall surface. It was concluded that the computational approach might have the potential to analyze the flow steadiness in a turbine blade flow.

Wang, Chi R.↗
NASA NTRS1980

On the interactions between large-scale structure and fine-grained turbulence in a free shear flow. III - A numerical solution

The results of a numerical computation of the interactions between the horizontally periodic monochromatic component of a large-scale coherent structure and the fine-grained turbulence in a mixing layer are presented. In the numerical calculations, the appropriate dependent variable is one which comprises both the mean and the large-scale coherent structure. The dynamical equations obtained for such a total coherent structure quantity are identical to the unsteady equations for the mean quantities in the Reynolds sense, except that the fine-grained turbulent stresses are interpreted as being conditionally averaged. It is shown how three dimensional fine-grained turbulence can be produced indirectly from the two dimensional large-scale structure via the isotropizing process of the approximated pressure-strain correlation.

Gatski, T. B.↗
NASA NTRS1967

Computer used to program numerically controlled milling machine

Computer program automatically directs a numerically controlled milling machine through a series of cutting and trimming actions. It accepts engineering data points, passes smooth curve segments through the points, breaks the resulting curves into a series of closely spaced points, and transforms these points into the form required by the mechanism.

Harris, T. C.↗
DOE OSTI2025

Tropical amplitudes for colored Lagrangians

Recently a new formulation for scattering amplitudes in Tr(Φ 3 ) theory has been given based on simple combinatorial ideas in the space of kinematic data. This allows all-loop integrated amplitudes to be expressed as “curve integrals” defined using tropical building blocks — the “headlight functions”. This paper shows how the formulation extends to the amplitudes of more general Lagrangians. We will present a number of different ways of introducing tropical “numerator functions” that allow us to describe general Lagrangian interactions. The simplest family of these “tropical numerators” computes the amplitudes of interesting Lagrangians with infinitely many interactions. We also describe methods for tropically formulating the amplitudes for general Lagrangians. One uses a variant of “Wick contraction” to glue together numerator factors for general interaction vertices. Another uses a natural characterization of polygons on surfaces to give a novel combinatorial description of all possible diagrams associated with arbitrary valence interactions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1992

An efficient HZETRN (a galactic cosmic ray transport code)

An accurate and efficient engineering code for analyzing the shielding requirements against the high-energy galactic heavy ions is needed. The HZETRN is a deterministic code developed at Langley Research Center that is constantly under improvement both in physics and numerical computation and is targeted for such use. One problem area connected with the space-marching technique used in this code is the propagation of the local truncation error. By improving the numerical algorithms for interpolation, integration, and grid distribution formula, the efficiency of the code is increased by a factor of eight as the number of energy grid points is reduced. The numerical accuracy of better than 2 percent for a shield thickness of 150 g/cm(exp 2) is found when a 45 point energy grid is used. The propagating step size, which is related to the perturbation theory, is also reevaluated.

Shinn, Judy L.↗
NASA NTRS2014

APC: A New Code for Atmospheric Polarization Computations

A new polarized radiative transfer code Atmospheric Polarization Computations (APC) is described. The code is based on separation of the diffuse light field into anisotropic and smooth (regular) parts. The anisotropic part is computed analytically. The smooth regular part is computed numerically using the discrete ordinates method. Vertical stratification of the atmosphere, common types of bidirectional surface reflection and scattering by spherical particles or spheroids are included. A particular consideration is given to computation of the bidirectional polarization distribution function (BPDF) of the waved ocean surface.

surface reflection↗
NASA NTRS1990

The route to chaos for the Kuramoto-Sivashinsky equation

The results of extensive numerical experiments of the spatially periodic initial value problem for the Kuramoto-Sivashinsky equation. This paper is concerned with the asymptotic nonlinear dynamics at the dissipation parameter decreases and spatio-temporal chaos sets in. To this end the initial condition is taken to be the same for all numerical experiments (a single sine wave is used) and the large time evolution of the system is followed numerically. Numerous computations were performed to establish the existence of windows, in parameter space, in which the solution has the following characteristics as the viscosity is decreased: a steady fully modal attractor to a steady bimodal attractor to another steady fully modal attractor to a steady trimodal attractor to a periodic attractor, to another steady fully modal attractor, to another periodic attractor, to a steady tetramodal attractor, to another periodic attractor having a full sequence of period-doublings (in parameter space) to chaos. Numerous solutions are presented which provide conclusive evidence of the period-doubling cascades which precede chaos for this infinite-dimensional dynamical system. These results permit a computation of the length of subwindows which in turn provide an estimate for their successive ratios as the cascade develops. A calculation based on the numerical results is also presented to show that the period doubling sequences found here for the Kuramoto-Sivashinsky equation, are in complete agreement with Feigenbaum's universal constant of 4,669201609... . Some preliminary work shows several other windows following the first chaotic one including periodic, chaotic, and a steady octamodal window; however, the windows shrink significantly in size to enable concrete quantitative conclusions to be made.

Papageorgiou, Demetrios T.↗
NASA NTRS1991

The route to chaos for the Kuramoto-Sivashinsky equation

The results of extensive numerical experiments of the spatially periodic initial value problem for the Kuramoto-Sivashinsky equation. This paper is concerned with the asymptotic nonlinear dynamics at the dissipation parameter decreases and spatio-temporal chaos sets in. To this end the initial condition is taken to be the same for all numerical experiments (a single sine wave is used) and the large time evolution of the system is followed numerically. Numerous computations were performed to establish the existence of windows, in parameter space, in which the solution has the following characteristics as the viscosity is decreased: a steady fully modal attractor to a steady bimodal attractor to another steady fully modal attractor to a steady trimodal attractor to a periodic attractor, to another steady fully modal attractor, to another periodic attractor, to a steady tetramodal attractor, to another periodic attractor having a full sequence of period-doublings (in parameter space) to chaos. Numerous solutions are presented which provide conclusive evidence of the period-doubling cascades which precede chaos for this infinite-dimensional dynamical system. These results permit a computation of the length of subwindows which in turn provide an estimate for their successive ratios as the cascade develops. A calculation based on the numerical results is also presented to show that the period doubling sequences found here for the Kuramoto-Sivashinsky equation, are in complete agreement with Feigenbaum's universal constant of 4,669201609 .... Some preliminary work shows several other windows following the first chaotic one including periodic, chaotic, and a steady octamodal window; however, the windows shrink significantly in size to enable concrete quantitative conclusions to be made.

Papageorgiou, Demetrios T.↗
NASA NTRS1986

A simple and fast method for computing the relativistic Compton Scattering Kernel for radiative transfer

The Klein-Nishina differential cross section averaged over a relativistic Maxwellian electron distribution is analytically reduced to a single integral, which can then be rapidly evaluated in a variety of ways. A particularly fast method for numerically computing this single integral is presented. This is, to the authors' knowledge, the first correct computation of the Compton scattering kernel.

Kershaw, David S.↗
DOE OSTI2023

An analytic approach to quasinormal modes for coupled linear systems

Quasinormal modes describe the ringdown of compact objects deformed by small perturbations. In generic theories of gravity that extend General Relativity, the linearized dynamics of these perturbations is described by a system of coupled linear differential equations of second order. We first show, under general assumptions, that such a system can be brought to a Schrödinger-like form. We then devise an analytic approximation scheme to compute the spectrum of quasinormal modes. We validate our approach using a toy model with a controllable mixing parameter ε and showing that the analytic approximation for the fundamental mode agrees with the numerical computation when the approximation is justified. The accuracy of the analytic approximation is at the (sub-) percent level for the real part and at the level of a few percent for the imaginary part, even when ε is of order one. Our approximation scheme can be seen as an extension of the approach of Schutz and Will [1] to the case of coupled systems of equations, although our approach is not phrased in terms of a WKB analysis, and offers a new viewpoint even in the case of a single equation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗