Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical approximations”

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 631 records · Page 35
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

Numerical investigation of supersonic turbulent boundary layers with high wall temperature

A direct numerical approach has been developed to simulate supersonic turbulent boundary layers. The mean flow quantities are obtained by solving the parabolized Reynolds-averaged Navier-Stokes equations (globally). Fluctuating quantities are computed locally with a temporal direct numerical simulation approach, in which nonparallel effects of boundary layers are partially modeled. Preliminary numerical results obtained at the free-stream Mach numbers 3, 4.5, and 6 with hot-wall conditions are presented. Approximately 5 million grid points are used in all three cases. The numerical results indicate that compressibility effects on turbulent kinetic energy, in terms of dilatational dissipation and pressure-dilatation correlation, are small. Due to the hot-wall conditions the results show significant low Reynolds number effects and large streamwise streaks. Further simulations with a bigger computational box or a cold-wall condition are desirable.

Guo, Y.↗
NASA NTRS1996

Wavelet Sparse Approximate Inverse Preconditioners

There is an increasing interest in using sparse approximate inverses as preconditioners for Krylov subspace iterative methods. Recent studies of Grote and Huckle and Chow and Saad also show that sparse approximate inverse preconditioner can be effective for a variety of matrices, e.g. Harwell-Boeing collections. Nonetheless a drawback is that it requires rapid decay of the inverse entries so that sparse approximate inverse is possible. However, for the class of matrices that, come from elliptic PDE problems, this assumption may not necessarily hold. Our main idea is to look for a basis, other than the standard one, such that a sparse representation of the inverse is feasible. A crucial observation is that the kind of matrices we are interested in typically have a piecewise smooth inverse. We exploit this fact, by applying wavelet techniques to construct a better sparse approximate inverse in the wavelet basis. We shall justify theoretically and numerically that our approach is effective for matrices with smooth inverse. We emphasize that in this paper we have only presented the idea of wavelet approximate inverses and demonstrated its potential but have not yet developed a highly refined and efficient algorithm.

Chan, Tony F.↗
NASA NTRS1984

Error analysis of analytic solutions for self-excited near-symmetric rigid bodies - A numerical study

Analytic error bounds are presented for the solutions of approximate models for self-excited near-symmetric rigid bodies. The error bounds are developed for analytic solutions to Euler's equations of motion. The results are applied to obtain a simplified analytic solution for Eulerian rates and angles. The results of a sample application of the range and error bound expressions for the case of the Galileo spacecraft experiencing transverse torques demonstrate the use of the bounds in analyses of rigid body spin change maneuvers.

Kia, T.↗
NASA NTRS2006

RAISHIN: A High-Resolution Three-Dimensional General Relativistic Magnetohydrodynamics Code

We have developed a new three-dimensional general relativistic magnetohydrodynamic (GRMHD) code, RAISHIN, using a conservative, high resolution shock-capturing scheme. The numerical fluxes are calculated using the Harten, Lax, & van Leer (HLL) approximate Riemann solver scheme. The flux-interpolated, constrained transport scheme is used to maintain a divergence-free magnetic field. In order to examine the numerical accuracy and the numerical efficiency, the code uses four different reconstruction methods: piecewise linear methods with Minmod and MC slope-limiter function, convex essentially non-oscillatory (CENO) method, and piecewise parabolic method (PPM) using multistep TVD Runge-Kutta time advance methods with second and third-order time accuracy. We describe code performance on an extensive set of test problems in both special and general relativity. Our new GRMHD code has proven to be accurate in second order and has successfully passed with all tests performed, including highly relativistic and magnetized cases in both special and general relativity.

Mizuno, Yosuke↗
NASA NTRS2003

An Approximate Dissipation Function for Large Strain Rubber Thermo-Mechanical Analyses

Mechanically induced viscoelastic dissipation is difficult to compute. When the constitutive model is defined by history integrals, the formula for dissipation is a double convolution integral. Since double convolution integrals are difficult to approximate, coupled thermo-mechanical analyses of highly viscous rubber-like materials cannot be made with most commercial finite element software. In this study, we present a method to approximate the dissipation for history integral constitutive models that represent Maxwell-like materials without approximating the double convolution integral. The method requires that the total stress can be separated into elastic and viscous components, and that the relaxation form of the constitutive law is defined with a Prony series. Numerical data is provided to demonstrate the limitations of this approximate method for determining dissipation. Rubber cylinders with imbedded steel disks and with an imbedded steel ball are dynamically loaded, and the nonuniform heating within the cylinders is computed.

Johnson, Arthur R.↗
NASA NTRS2018

The Magellan Evolution of Galaxies Spectroscopic and Ultraviolet Reference Atlas (MegaSaura). II. Stacked Spectra

We stack the rest-frame ultraviolet spectra of N = 14 highly magnified gravitationally lensed galaxies at redshifts 1.6 < z < 3.6. The resulting new composite spans 900 < lambda rest < 3000 Angstrom, with a peak signal-to-noise ratio of 103 per spectral resolution element (approximately 100 km s (exp -1)). It is the highest signalto-noise ratio, highest spectral resolution composite spectrum of z approximately 2-3 galaxies yet published. The composite reveals numerous weak nebular emission lines and stellar photospheric absorption lines that can serve as new physical diagnostics, particularly at high redshift with the James Webb Space Telescope (JWST). We report equivalent widths to aid in proposing for and interpreting JWST spectra. We examine the velocity profiles of strong absorption features in the composite, and in a matched composite of z approximately 0 COS/HST galaxy spectra. We find remarkable similarity in the velocity profiles at z approximately 0 and z approximately 2, suggesting that similar physical processes control the outflows across cosmic time. While the maximum outflow velocity depends strongly on ionization potential, the absorptionweighted mean velocity does not. As such, the bulk of the high-ionization absorption traces the low-ionization gas, with an additional blueshifted absorption tail extending to at least approximately 2000 km s (exp -1). We interpret this tail as arising from the stellar wind and photospheres of massive stars. Starburst99 models are able to replicate this high-velocity absorption tail. However, these theoretical models poorly reproduce several of the photospheric absorption features, indicating that improvements are needed to match observational constraints on the massive stellar content of star-forming galaxies at z approximately 2. We publicly release our composite spectra.

Rigby, J. R.↗
NASA NTRS1996

Intergalactic Hydrogen Clouds at Low Redshift: Connections to Voids and Dwarf Galaxies

We provide new post-COSTAR data on one sightline (Mrk 421) and updated data from another (I Zw 1) from our Hubble Space Telescope (HST) survey of intergalactic Ly(alpha) clouds located along sightlines to four bright quasars passing through well-mapped galaxy voids (16000 km/s pathlength) and superclusters (18000 km/s). We report two more definite detections of low-redshift Ly(alpha) clouds in voids: one at 3047 km/s (heliocentric) toward Mrk 421 and a second just beyond the Local Supercluster at 2861 km/s toward I Zw 1, confirming our earlier discovery of Ly(alpha) absorption clouds in voids (Stocke et al., ApJ, 451, 24). We have now identified ten definite and one probable low-redshift neutral hydrogen absorption clouds toward four targets, a frequency of approximately one absorber every 3400 km/s above 10(exp 12.7/sq cm column density. Of these ten absorption systems, three lie within voids; the probable absorber also lies in a void. Thus, the tendency of Ly(alpha) absorbers to 'avoid the voids' is not as clear as we found previously. If the Ly(alpha) clouds are approximated as homogeneous spheres of 100 kpc radius, their masses are approximately 10(exp 9)solar mass (about 0.01 times that of bright L* galaxies) and they are 40 times more numerous, comparable to the density of dwarf galaxies and of low-mass halos in numerical CDM simulations. The Ly(alpha) clouds contribute a fraction Omega(sub cl)approximately equals 0.003/h(sub 75) to the closure density of the universe, comparable to that of luminous matter. These clouds probably require a substantial amount of nonbaryonic dark matter for gravitational binding. They may represent extended haloes of low-mass protogalaxies which have not experienced significant star formation or low-mass dwarf galaxies whose star formation ceased long ago, but blew out significant gaseous material.

Shull, J. Michael↗
NASA NTRS1972

A discrete variable approximation to minimum weight panel designs subject to a supersonic flutter-speed constraint.

A numerical technique is presented that determines the minimum-weight thickness distribution for a one-dimensional, simply supported panel subjected to a parallel supersonic flow on one side. An aeroelastic eigenvalue, which characterizes flutter speed, is held constant. The governing differential equations are approximated by sets of difference equations adjoined to the weight function via a penalty function. A conjugate gradient method then solves the resulting sequence of unconstrained minimization problems. Numerical results are obtained without imposing a minimum thickness constraint; the effect of constant inplane tensile and compressive stresses is also investigated.

Pierson, B. L.↗
NASA NTRS1973

The particle resonance in spiral galaxies - Nonlinear effects.

A theory is developed to account for the nonlinear effects, near the particle resonance, found by numerical integration. There are four equilibrium points in the rotating frame of reference: two of them unstable, at the minima of potential (L1, L2), and two stable, at the maxima of potential (L4, L5). Many particles are trapped in librating orbits around L4, L5. Using the lowest-order terms of a 'third integral' of motion near L4, L5, the behavior of the trapped orbits is determined, and the approximate theoretical results are compared with the orbits found numerically by computer. An estimate of the trapped mass is given. It amounts to about 30% of the total mass between 9 and 11 kpc when the particle resonance is at r sub s = 10 kpc and the maximum force due to the spiral is 4% of the axisymmetric force. Some effects due to these mass concentrations are found numerically and theoretically.

Contopoulos, G.↗
NASA NTRS1989

Markov reliability models for digital flight control systems

The reliability of digital flight control systems can often be accurately predicted using Markov chain models. The cost of numerical solution depends on a model's size and stiffness. Acyclic Markov models, a useful special case, are particularly amenable to efficient numerical solution. Even in the general case, instantaneous coverage approximation allows the reduction of some cyclic models to more readily solvable acyclic models. After considering the solution of single-phase models, the discussion is extended to phased-mission models. Phased-mission reliability models are classified based on the state restoration behavior that occurs between mission phases. As an economical approach for the solution of such models, the mean failure rate solution method is introduced. A numerical example is used to show the influence of fault-model parameters and interphase behavior on system unreliability.

Mcgough, John↗
NASA NTRS1961

Laminar Boundary Layer Behind a Strong Shock Moving into Air

The laminar wall boundary layer behind a strong shock advancing into stationary air has been determined. Numerical results have been obtained for shock Mach numbers up to 14 using real gas values for density and viscosity and assuming Prandtl and Lewis numbers of 0.72 and 1, respectively. The numerical results for shear and heat transfer agree, within 4 percent, with a previously presented approximate analytical expression for these quantities. A slight modification of this expression results in agreement with the numerical data to within 2.5 percent. Analytical expressions for boundary-layer thickness and displacement thickness, correct to within 4 percent for the present data, have also been obtained.

Mirels, Harold↗
NASA NTRS1971

Pressurized expulsion of nonisothermal single-phase cryogen

The performance of single-phase storage and expulsion systems as affected by temperature variations within the stored cryogen which are generated during heat transfer is discussed. Peculiar operating responses are indicated by spontaneous changes in fluid pressure which accompany g level changes, increased heater surface temperature, and durations of pressure cycles which differ considerably from that which is computed for an isothermal cryogen. The nonisothermal characteristics are predicted with a numerical model which includes the simultaneous solution of the time dependent conservation equations of mass, energy, and momentum in two space dimensions of Cartesian coordinates for boundary conditions which approximate those of the flight cryogenic system. The methodology of the numerical method and some comparisons between the predictions and the Apollo 12 flight data are included.

Forester, C. K.↗
NASA NTRS1971

Reverse and forward slow shocks in the solar wind.

Probable reverse and forward slow shocks were found in plasma and magnetic-field data from Pioneer 6. The shocks were oblique and weak (Mach of approximately 1.2). No reverse fast shocks were found. Numerous (50) other discontinuities were found, most of which are probably tangential, since all the plasma parameters changed across most of them.

Burlaga, L. F.↗
NASA NTRS1973

Round-off errors in cutting plane algorithms based on the revised simplex procedure

This report statistically analyzes computational round-off errors associated with the cutting plane approach to solving linear integer programming problems. Cutting plane methods require that the inverse of a sequence of matrices be computed. The problem basically reduces to one of minimizing round-off errors in the sequence of inverses. Two procedures for minimizing this problem are presented, and their influence on error accumulation is statistically analyzed. One procedure employs a very small tolerance factor to round computed values to zero. The other procedure is a numerical analysis technique for reinverting or improving the approximate inverse of a matrix. The results indicated that round-off accumulation can be effectively minimized by employing a tolerance factor which reflects the number of significant digits carried for each calculation and by applying the reinversion procedure once to each computed inverse. If 18 significant digits plus an exponent are carried for each variable during computations, then a tolerance value of 0.1 x 10 to the minus 12th power is reasonable.

Moore, J. E.↗
NASA NTRS1974

Post-flight differential correction analysis using Vinti's spheroidal method for the small astronomy satellite orbit

The results of an intensive analysis of a differential orbit improvement method utilizing observational data for a 550-kilometer altitude, near-circular, near-equatorial satellite orbit are presented. Observations of the Small Astronomy Satellite (SAS-1) are in the form of direction cosines as measured at two ground interferometer tracking stations near the Equator during the first 22 orbital revolutions (approximately 37 hours) after launch of the spacecraft. Numerical results, in both tabular and graphical form, are displayed for numerous iterated fittings of various observational arcs by differential correction of the orbital elements. Parameters varied in these comparative cases include the time duration of the observational data block, the number of pairs of direction cosine data and the number of tracking station passes included in the solution, the distribution of such passes between the two available tracking stations, and the acceptance criterion for the observational residuals in the least squares fitting procedure. It was found that three observational pairs of direction cosine data, the minimum number possible for a uniquely determined solution in theory, are sufficient to promote covergence to an accurate solution, if properly selected.

Walden, H.↗
NASA NTRS1974

On laminar boundary-layer separation

Iterative finite-difference techniques are developed for integrating the boundary-layer equations, without approximation, through a region of reversed flow. The numerical procedures are used to calculate incompressible laminar separated flows and to investigate the conditions for regular behavior at the point of separation. Regular flows are shown to be characterized by an integrable saddle-type singularity that makes it difficult to obtain numerical solutions which pass continuously into the separated region. The singularity is removed and continuous solutions ensured by specifying the wall shear distribution and computing the pressure gradient as part of the solution. Calculated results are presented for a number of separated flows and the accuracy of the method is verified.

Klineberg, J. M.↗
NASA NTRS1979

A group-velocity criterion for breakdown of vortex flow: An application to measured inlet profiles

Vortex flows exhibiting breakdown in a slightly divergent duct were measured. The slowly varying vortex flow field downstream of the entrance and upstream of the breakdown region is obtained numerically by using the inviscid quasi-cylindrical approximation. In these calculations, the Faler and Lebovich's experimental data were used as the starting conditions at the entrance of the duct. The group velocity of wave propagation for the axisymmetric mode (n = 0) and the asymmetric modes (n = + or - 1 and n = + or - 2) are calculated for the entrance conditions. For the theoretically predicted slowly varying flow field downstream of the entrance, the wave characteristics of the n = 0 and n = + or - 1 modes are presented. It was concluded that the flows which subsequently undergo vortex breakdown are all predicted to be supercritical and stable to infinitesimal inviscid disturbances, including the axially symmetric as well as the nonsymmetric perturbations.

Tsai, C. Y.↗