Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical algorithms”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 883 records · Page 49
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Modeling flow at the nozzle of a solid rocket motor

The mechanical behavior of a rocket motor internal flow field results in a system of nonlinear partial differential equations which can be solved numerically. The accuracy and the convergence of the solution of the system of equations depends largely on how precisely the sharp gradients can be resolved. An adaptive grid generation scheme is incorporated into the computer algorithm to enhance the capability of numerical modeling. With this scheme, the grid is refined as the solution evolves. This scheme significantly improves the methodology of solving flow problems in rocket nozzle by putting the refinement part of grid generation into the computer algorithm.

Chow, Alan S.↗
NASA NTRS1992

Boundary layer study on nozzle wall at hypersonic velocities

The boundary layer on the wall of the Hypervelocity Tunnel 9 was investigated with pitot pressure and total temperature measurements. Experimental results are presented for standard and supercooled Mach 14 runs. The boundary layer data at supercooled conditions are compared to numerical predictions made with a Navier-Stokes algorithm including vibrational nonequilibrium and intermolecular force effects. For standard tunnel conditions, the numerical solutions agree well with experimental data. For the supercooled cases, the numerical code predicts the total temperature but overpredicts the pitot pressure.

Jones, Kenneth M.↗
NASA NTRS1993

A study of the compressible flow through a diffusing S-duct

Benchmark aerodynamic data are presented for compressible flow through a representative S-duct configuration. A numerical prediction of the S-duct flow field, obtained from a subsonic parabolized Navier-Stokes algorithm, is also shown. The experimental and numerical results are compared. Measurements of the three-dimensional velocity field, total pressures, and static pressures were obtained at five cross-sectional planes. Aerodynamic data were gathered with calibrated pneumatic probes. Surface static pressure and surface flow visualization data were also acquired. All reported tests were conducted with an inlet centerline Mach number of 0.6. The Reynolds number, based on the inlet centerline velocity and duct inlet diameter, was 2.6 x 10(exp 6). Thin inlet turbulent boundary layers existed. The collected data should be beneficial to aircraft inlet designers and the measurements are suitable for the validation of computational codes. The results show that a region of streamwise flow separation occurred within the duct. Details about the separated flow region, including mechanisms which drive this complicated flow phenomenon, are discussed. Results also indicate that the duct curvature induces strong pressure driven secondary flows. The cross flows evolve into counter-rotating vortices. These vortices convect low momentum fluid of the boundary layer toward the center of the duct, degrading both the uniformity and magnitude of the total pressure profile.

Wellborn, Steven R.↗
NASA NTRS2019

An Approach to Retrieve BRDF from Satellite and Airborne Measurements of Surface-Reflected Radiance Based on Decoupling of Atmospheric Radiative Transfer and Surface Reflection

Bi-directional Reflection Distribution Function (BRDF) defines anisotropy of the surface reflection. It is required to specify the boundary condition for radiative transfer (RT) modeling. Measurements of reflected radiance by satellite- and air-borne sensors provide information about anisotropy of surface reflection. Atmospheric correction needs to be performed to derive BRDF from the reflected radiance. Common approach for BRDF retrievals consists of the use of kernel-based BRDF and RT modeling that needs to be done anew at every step of the iterative process. The kernels’ weights are obtained by minimization of the difference between measured and modeled radiance. This study develops a new method of retrieving kernel-based BRDF that requires RT calculations to be done only once. The method employs the exact analytical expression of radiance at any atmospheric level through the solutions of two auxiliary atmosphere-only RT problems and the surface-reflected radiance at the surface level. The latter is related to BRDF and solutions of the auxiliary RT problems by a Fredholm integral equation of the second kind. The approach requires to perform RT calculations one time before the iterations. It can use observations taken at different atmospheric conditions assuming that surface conditions remain unchanged during the time span of observations. The algorithm accurately catches zero weights of the kernels that may be a concern if the number of kernels is greater than 3 in current mainstream approaches. The study presents numerical tests of the BRDF retrieval algorithm for various surface and atmospheric conditions.

Radkevich, Alexander↗
NASA NTRS1988

The use of multigrid techniques in the solution of the Elrod algorithm for a dynamically loaded journal bearing

A numerical solution to a theoretical model of vapor cavitation in a dynamically loaded journal bearing is developed, utilizing a multigrid iterative technique. The code is compared with a presently existing direct solution in terms of computational time and accuracy. The model is based on the Elrod algorithm, a control volume approach to the Reynolds equation which mimics the Jakobssen-Floberg and Olsson cavitation theory. Besides accounting for a moving cavitation boundary and conservation of mass at the boundary, it also conserves mass within the cavitated region via liquid striations. The mixed nature of the equations (elliptic in the full film zone and nonelliptic in the cavitated zone) coupled with the dynamic aspects of the problem create interesting difficulties for the present solution approach. Emphasis is placed on the methods found to eliminate solution instabilities. Excellent results are obtained for both accuracy and reduction of computational time.

Woods, Claudia M.↗
NASA NTRS1979

Studies in astronomical time series analysis: Modeling random processes in the time domain

Random process models phased in the time domain are used to analyze astrophysical time series data produced by random processes. A moving average (MA) model represents the data as a sequence of pulses occurring randomly in time, with random amplitudes. An autoregressive (AR) model represents the correlations in the process in terms of a linear function of past values. The best AR model is determined from sampled data and transformed to an MA for interpretation. The randomness of the pulse amplitudes is maximized by a FORTRAN algorithm which is relatively stable numerically. Results of test cases are given to study the effects of adding noise and of different distributions for the pulse amplitudes. A preliminary analysis of the optical light curve of the quasar 3C 273 is given.

Scargle, J. D.↗
NASA NTRS1981

Simulation of large motions of nonuniform beams in orbit. I - The cantilever beam

An algorithm is developed for producing numerical simulations of large motions of a nonuniform cantilever beam in orbit. Special emphasis is given to the effective formulation of equations of motion and to the use of the finite element method to construct modal functions. Finite element methods are used to generate modal functions in such a way as to permit a particularly harmonious relationship to be established between the disciplines of rigid body dynamics and structural analysis. Sufficient information is provided to enable a reader to create, with relatively little effort, his own simulation program; simulation results are reported, both to provide check cases for other investigators and to illustrate certain important facets of the behavior of flexible spacecraft.

Levinson, D. A.↗
NASA NTRS1981

Simulation of large motions of nonuniform beams in orbit. II - The unrestrained beam

An algorithm is developed for the numerical simulation of large motions of an unrestrained, nonuniform beam in orbit. Explicit provisions are made for the accommodation of general, physically significant initial conditions and for the specification of control forces in practical terms. The evaluation of certain constants that appear in the equations of motion is considered, the finite element method being used for this purpose. The importance of initial conditions in the study of spacecraft motions is stressed.

Kane, T. R.↗
NASA NTRS1982

Measurement accuracies in band-limited extrapolation

The problem of numerical instability associated with extrapolation algorithms is addressed. An attempt is made to estimate the bounds for the acceptable errors and to place a ceiling on the measurement accuracy and computational accuracy needed for the extrapolation. It is shown that in band limited (or visible angle limited) extrapolation the larger effective aperture L' that can be realized from a finite aperture L by over sampling is a function of the accuracy of measurements. It is shown that for sampling in the interval L/b absolute value of xL, b1 the signal must be known within an error e sub N given by e sub N squared approximately = 1/4(2kL') cubed (e/8b L/L')(2kL') where L is the physical aperture, L' is the extrapolated aperture, and k = 2pi lambda.

Kritikos, H. N.↗
NASA NTRS1982

An implicit solution of the three-dimensional Navier-Stokes equations for an airfoil spanning a wind tunnel

An implicit finite-difference algorithm is developed for the numerical solution of the incompressible three dimensional Navier-Stokes equations in the non-conservative primitive-variable formulation. The flow field about an airfoil spanning a wind-tunnel is computed. The coordinate system is generated by an extension of the two dimensional body-fitted coordinate generation techniques of Thompson, as well as that of Sorenson, into three dimensions. Two dimensional grids are stacked along a spanwise coordinate defined by a simple analytical function. A Poisson pressure equation for advancing the pressure in time is arrived at by performing a divergence operation on the momentum equations. The pressure at each time-step is calculated on the assumption that continuity be unconditionally satisfied. An eddy viscosity coefficient, computed according to the algebraic turbulence formulation of Baldwin and Lomax, simulates the effects of turbulence.

Moitra, A.↗
NASA NTRS1984

Finite element analysis of steadily moving contact fields

By introducing a moving updated Lagrangian observer, this paper develops traveling finite elements with the capacity to handle the global response resulting from steadily moving contact fields. The generality of the results is such that large deformation kinematics and kinetics as well as the full compliment of inertial fields can be handled. To streamline the handling of nonlinear behavior, an elliptically constrained solution algorithm is also developed. Employing this algorithm, the results of several numerical benchmarking studies are presented which illustrate the capacity of the moving updated Lagrangian formulation as well as the potential effects of nonlinearity.

Padovan, J.↗
NASA NTRS1984

Estimation of discontinuous coefficients in parabolic systems: Applications to reservoir simulation

Spline based techniques for estimating spatially varying parameters that appear in parabolic distributed systems (typical of those found in reservoir simulation problems) are presented. The problem of determining discontinuous coefficients, estimating both the functional shape and points of discontinuity for such parameters is discussed. Convergence results and a summary of numerical performance of the resulting algorithms are given.

Lamm, P. D.↗
NASA NTRS1984

Intersection of parametric surfaces using lookup tables

When primitive structures in the form of parametric surfaces are combined and modified interactively to form complex intersecting surfaces, it becomes important to find the curves of intersection. One must distinguish between finding the shape of the intersection curve, which may only be useful for display purposes, and finding an accurate mathematical representation of the curve, which is important for any meaningful geometric modeling, analysis, design, or manufacturing involving the intersection. The intersection curve between two or more parametric surfaces is important in a variety of computer-aided design and manufacture areas. A few examples are shape design, analysis of groins, design of fillets, and computation of numerically controlled tooling paths. The algorithm presented here provides a mathematical representation of the intersection curve to a specified accuracy. It also provides the database that can simplify operations such as hidden-surface removal, surface rendering, profile identification, and interference or clearance computations.

Hanna, S. L.↗
NASA NTRS1984

Retrieval of vertical moisture profiles with microwave radiometry

A statistical correlation technique is applied to the retrieval of vertical moisture profiles from downlooking radiometric measurements of atmospheric radiation at microwave wavelengths. Only an optimum subset of available radiometer channels is selected for estimating water vapor at specific pressure levels. To test its validity the algorithm was applied, in a numerical experiment, to 50 independent tropical radiosondes over a sea surface. It was also used to retrieve continuous sequences of atmospheric moisture profiles from a set of data obtained with a 4-channel microwave radiometer carried aboard an aircraft over a land surface.

Kakar, R. K.↗
NASA NTRS1984

High temperature thermomechanical analysis of ceramic coatings

This paper investigates the thermomechanical response of ceramically coated metal parts in elevated thermal environments. This is made possible through the development of an improved finite element algorithm that enables the efficient and stable solution of the inherently nonlinear elastic-creep (inelastic) type thermomechanical field equations associated with high temperature. Based on the improved algorithm, the results of several numerical experiments are presented. These illustrate the significant influence of inelastic behavior in generating residual stress fields.

Padovan, J.↗
NASA NTRS1984

A statistical correlation method for the retrieval of atmospheric moisture profiles by microwave radiometry

A statistical correlation technique is applied to the retrieval of vertical moisture profiles under clear-sky conditions from down-looking radiometric measurements of atmospheric radiation at microwave wavelengths. For a given set of channels, the method selects the optimum radiometric channels for estimating water vapor at specific pressure levels between the surface and 300 mb. The water vapor mixing ratio at these pressure levels is then calculated from a linear combination of the selected channel brightness temperatures. To test its validity the algorithm was applied, in a numerical experiment, to fifty independent tropical radiosondes. The rms absolute deviation of the estimated moisture profiles from the actual profiles was comparable to that obtained using an iterative retrieval method reported earlier. The statistical method, however, requires several orders of magnitude less computer time than the iterative method; it is suitable for high speed processing of large amounts of data.

Kakar, R. K.↗
NASA NTRS1985

Generalized formulation of a class of explicit and implicit TVD schemes

A one parameter family of second order explicit and implicit total variation diminishing (TVD) schemes is reformulated so that a simpler and wider group of limiters is included. The resulting scheme can be viewed as a symmetrical algorithm with a variety of numerical dissipation terms that are designed for weak solutions of hyperbolic problems. This is a generalization of Roe and Davis's recent works to a wider class of symmetric schemes other than Lax-Wendroff. The main properties of the present class of schemes are that they can be implicit, and when steady state calculations are sought, the numerical solution is independent of the time step.

Yee, H. C.↗
NASA NTRS1985

CAG12 - A CSCM based procedure for flow of an equilibrium chemically reacting gas

The Conservative Supra Characteristic Method (CSCM), an implicit upwind Navier-Stokes algorithm, is extended to the numerical simulation of flows in chemical equilibrium. The resulting computer code known as Chemistry and Gasdynamics Implicit - Version 2 (CAG12) is described. First-order accurate results are presented for inviscid and viscous Mach 20 flows of air past a hemisphere-cylinder. The solution procedure captures the bow shock in a chemically reacting gas, a technique that is needed for simulating high altitude, rarefied flows. In an initial effort to validate the code, the inviscid results are compared with published gasdynamic and chemistry solutions and satisfactorily agreement is obtained.

Green, M. J.↗