Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numeric profiling”

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 541 records · Page 30

Viscous, resistive magnetohydrodynamic stability computed by spectral techniques

Expansions in Chebyshev polynomials are used to study the linear stability of one-dimensional magnetohydrodynamic quasi-equilibria, in the presence of finite resistivity and viscosity. The method is modeled on the one used by Orszag in accurate computation of solutions of the Orr-Sommerfeld equation. Two Reynolds-like numbers involving Alfven speeds, length scales, kinematic viscosity, and magnetic diffusivity govern the stability boundaries, which are determined by the geometric mean of the two Reynolds-like numbers. Marginal stability curves, growth rates versus Reynolds-like numbers, and growth rates versus parallel wave numbers are exhibited. A numerical result that appears general is that instability has been found to be associated with inflection points in the current profile, though no general analytical proof has emerged. It is possible that nonlinear subcritical three-dimensional instabilities may exist, similar to those in Poiseuille and Couette flow.

Dahlburg, R. B.↗

Far infrared remote sounding of stratospheric temperature and trace gas distributions

An evaluation is made of the potential usefulness of far-IR thermal emission measurements for simultaneous retrieval of temperature and constituent distributions. Attention is given to the design and performance of a high resolution Michelson interferometer which is currently used in balloon-borne far-IR stratospheric studies. Numerical results based on synthetic limb radiance data for model atmospheres are presented. Formal inversion techniques, together with a radiative transfer model that is based on line-by-line transmittance calculations, are used. Numerical results based on synthetic limb radiance data are given in order to assess the retrieval accuracies of the temperature profiles and the vertical mixing ratio profiles of O3, HCl, and OH.

Abbas, M. M.↗

Shock structure in non-circular jets

The shock-cell structure of supersonic jets with non-circular exit geometry is modeled using a linearized analysis. The model takes into account the finite thickness of the jet shear layer using realistic velocity and density profiles. The effects of the shear layer turbulence are included by incorporating eddy-viscosity terms. A finite-difference numerical method is used to solve the steady linearized equations of motion. A body-fitted coordinate system is used to describe the shear layer. The variation of the pressure fluctuation with downstream distance is given for circular jets and for an elliptic jet of aspect ratio 2.0. Comparisons with experimental data are made. Difficulties with the numerical technique are also discussed.

Morris, Philip J.↗

Heat Transfer to a Thin Solid Combustible in Flame Spreading at Microgravity

The heat transfer rate to a thin solid combustible from an attached diffusion flame, spreading across the surface of the combustible in a quiescent, microgravity environment, was determined from measurements made in the drop tower facility at NASA-Lewis Research Center. With first-order Arrhenius pyrolysis kinetics, the solid-phase mass and energy equations along with the measured spread rate and surface temperature profiles were used to calculate the net heat flux to the surface. Results of the measurements are compared to the numerical solution of the complete set of coupled differential equations that describes the temperature, species, and velocity fields in the gas and solid phases. The theory and experiment agree on the major qualitative features of the heat transfer. Some fundamental differences are attributed to the neglect of radiation in the theoretical model.

Bhattacharjee, S.↗

Integrated trajectory and control analysis for generic hypersonic vehicles

A tool which supports integrated assessment of air-breathing hypersonic vehicle trajectory management and control strategy has been developed by Hattis and Malchow (1991). Using a two-point boundary value problem solution technique, the tool can construct an entire near-minimum-fuel trajectory and desired control strategy from takeoff to orbit, while satisfying physically derived inequality constraints and while achieving efficient propulsive mode phasing. A strategy is also formulated to construct preliminary trajectory and control history representations with less computational burden than required for the overall flight profile assessment. A tabulated example hypersonic vehicle model has been used to demonstrate the integrated analysis methodology. To assure good numerical behavior when the algorithm uses tabulated model data, a multidimensional cubic spline curve smoothing routine has been developed. A dynamic-pressure-constrained, near-fuel-optimal trajectory was computed to provide examples of information that the methodology can provide. Previously unspecified propulsive discontinuities were located, flight regimes demanding rapid attitude changes were identified, available and closed-loop controller authority was ascertained, and inadequacies in vehicle and specific subsystem model representations were determined.

Hattis, Philip D.↗

Supersonic boundary-layer flow turbulence modeling

Baldwin-Lomax and kappa-epsilon turbulence models were modified for use in Navier-Stokes numerical computations of Mach 2.9 supersonic turbulent boundary layer flows along compression ramps. The computational results of Reynolds shear stress profiles were compared with experimental data. The Baldwin-Lomax model was modified to account for the Reynolds shear stress amplification within the flow field. A hybrid kappa-epsilon model with viscous sublayer turbulence treatment was constructed to predict the Reynolds shear stress profiles within the entire flow field. These modified turbulence models were effective for the computations of the surface pressure and the skin friction factor variations along an 8 deg ramp surface. The hybrid kappa-epsilon model could improve the predictions of the Reynolds shear stress profile and the skin friction factor near the corner of a 16 deg ramp.

Wang, Chi-Rong↗

Buoyancy-Driven Heat Transfer During Application of a Thermal Gradient for the Study of Vapor Deposition at Low Pressure Using and Ideal Gas

A mathematical model has been developed to determine heat transfer during vapor deposition of source materials under a variety of orientations relative to gravitational accelerations. The model demonstrates that convection can occur at total pressures as low as 10-2 mm Hg. Through numerical computation, using physical material parameters of air, a series of time steps demonstrates the development of flow and temperature profiles during the course of vapor deposition. These computations show that in unit gravity vapor deposition occurs by transport through a fairly complicated circulating flow pattern when applying heat to the bottom of the vessel with parallel orientation with respect to the gravity vector. The model material parameters for air predict the effect of kinematic viscosity to be of the same order as thermal diffusivity, which is the case for Prandtl number approx. 1 fluids. Qualitative agreement between experiment and the model indicates that 6-(2-methyl-4-nitroanilino)-2,4-hexadiyn-l-ol (DAMNA) at these pressures indeed approximates an ideal gas at the experiment temperatures, and may validate the use of air physical constants. It is apparent that complicated nonuniform temperature distribution in the vapor could dramatically affect the homogeneity, orientation, and quality of deposited films. The experimental test i's a qualitative comparison of film thickness using ultraviolet-visible spectroscopy on films generated in appropriately oriented vapor deposition cells. In the case where heating of the reaction vessel occurs from the top, deposition of vapor does not normally occur by convection due to a stable stratified medium. When vapor deposition occurs in vessels heated at the bottom, but oriented relative to the gravity vector between these two extremes, horizontal thermal gradients induce a complex flow pattern. In the plane parallel to the tilt axis, the flow pattern is symmetrical and opposite in direction from that where the vessel is positioned vertically. The ground-based experiments are sufficient preliminary tests of theory and should be of significant interest regarding vapor deposited films in microgravity.

Frazier, D. O.↗

Improved Prediction of Losses with Large Eddy Simulation in a Low Pressure Turbine

We will present the wake profiles resulting from our previous computations for a range of FSTs. We will address key questions such as how the numerics, mesh density in addition to the FST level affect the wake prediction. We will utilize the spectral analysis and investigate the dynamic behavior of the wake as influenced by the flow variations in locations upstream.

LES↗

Interpreting Broad Double-Peaked Emission Lines in Active Galactic Nuclei

The principal objectives of this project were to probe the inner regions of active galactic nuclei and to test general relativity in the strong-field limit. The approach takes advantage of broad atomic line emission observed from material deep in the potential well of an active galactic nucleus which contains key information as to the physics of the system. Line profiles in a wide range of wavebands from optical to X-ray have provided compelling evidence of the existence of a relativistic accretion disk around a supermassive black hole in a number of galaxies. The simplest model posits a geometrically thin disk in Keplerian orbit, with general relativistic effects in evidence. This model is the point of departure for the proposed work. We developed a high-performance numerical code to calculate photon trajectories in a Schwarzschild or Kerr metric and implemented it on parallel supercomputers. This code includes a general purpose ray tracer that calculates line profiles, light curves, and other observable quantities for a wide variety of emitter configurations. The versatility comes from the fact that the ray tracing algorithm does not depend on any symmetries regarding emitter locations. The speed comes from parallel implementation which enables us to sample hitherto unattainable volumes of disk model parameter space. During the period 1 March 1997 through 28 February 1998, two papers, supported in whole or in part by this grant, were published in refereed journals. They are reproduced in their entirety in the next two sections of this report.

Halpern, Jules↗

Turbulent Boundary Layers in Oscillating Flows. Part 1: an Experimental and Computational Study

An experimental-computational study of the behavior of turbulent boundary layers for oscillating air flows over a plane surface with a small favorable mean pressure gradient is described. Experimental studies were conducted for boundary layers generated on the test section wall of a facility that produces a flow with a mean free stream velocity and a superposed nearly-pure sinusoidal component over a wide range of frequency. Flow at a nominal mean free stream velocity of 50 m/s were studied at atmospheric pressure and temperature at selected axial positions over a 2 m test length for frequencies ranging from 4 to 29 Hz. Quantitative experimental results are presented for unsteady velocity profiles and longitudinal turbulence levels obtained from hot wire anemometer measurements at three axial positions. Mean velocity profiles for oscillating flows were found to exhibit only small deviations from corresponding steady flow profiles, while amplitudes and phase relationships exhibited a strong dependence on axial position and frequency. Since sinusoidal flows could be generated over a wide range of frequency, studies at fixed values of reduced frequency at different axial positions were studied. Results show that there is some utility in the use of reduced frequency to correlate unsteady velocity results. The turbulence level u' sub rms was observed to vary essentially sinusoidally around values close to those measured in steady flow. However, the amplitude of oscillation and phase relations for turbulence level were found to be strongly frequency dependent. Numerical predictions were obtained using an unsteady boundary layer computational code and the Cebeci-Smith and Glushko turbulence models. Predicted quantities related to unsteady velocity profiles exhibit fair agreement with experiment when the Cebeci-Smith turbulence model is used.

Cook, W. J.↗

Modeling of Convective-Stratiform Precipitation Processes: Sensitivity to Partitioning Methods and Numerical Advection Schemes

Six different convective-stratiform separation techniques, including a new technique that utilizes the ratio of vertical and terminal velocities, are compared and evaluated using two-dimensional numerical simulations of a tropical [Tropical Ocean Global Atmosphere Coupled Ocean-Atmosphere Response Experiment (TOGA COARE)] and midlatitude continental [Preliminary Regional Experiment for STORM-Central (PRESTORM)] squall line. The simulations are made using two different numerical advection schemes: 4th order and positive definite advection. Comparisons are made in terms of rainfall, cloud coverage, mass fluxes, apparent heating and moistening, mean hydrometeor profiles, CFADs (Contoured Frequency with Altitude Diagrams), microphysics, and latent heating retrieval. Overall, it was found that the different separation techniques produced results that qualitatively agreed. However, the quantitative differences were significant. Observational comparisons were unable to conclusively evaluate the performance of the techniques. Latent heating retrieval was shown to be sensitive to the use of separation technique mainly due to the stratiform region for methods that found very little stratiform rain. The midlatitude PRESTORM simulation was found to be nearly invariant with respect to advection type for most quantities while for TOGA COARE fourth order advection produced numerous shallow convective cores and positive definite advection fewer cells that were both broader and deeper penetrating above the freezing level.

Lang, Steve↗

Numerical solution of the two-dimensional Euler equations by second-order upwind difference schemes

Two time-level, five-point explicit and implicit upwind difference schemes based on the characteristic flux difference splitting concept have been developed for the two-dimensional Euler equations. The method is conservative, second-order accurate in time and space, and general coordinate systems are used to treat complex geometries. Nonlinear flux limiters are employed to yield oscillation free sharp shock profiles. Upstream interpolation is used to yield a class of higher-order upwind schemes which closely mimic the locally one-dimensional method of characteristics (with fixed time intervals) through operator splitting. Numerical results have been obtained for a plane shock reflection and for flow over a circular arc in a channel. Characteristics of upwind TVD schemes, as applied to two-dimensional flows with embedded shocks are discussed.

Yang, J. Y.↗

SINDA, Systems Improved Numerical Differencing Analyzer

Computer program has been written to analyze group of 100-node areas and then provide for summation of any number of 100-node areas to obtain temperature profile. SINDA program options offer user variety of methods for solution of thermal analog modes presented in network format.

Fink, L. C.↗

Near field antenna patterns of obstructed Cassegrainian telescopes

The near field antenna pattern of an obstructed optical transmitter, such as a Cassegrainian telescope, has been computed using the Fresnel approximation to the scalar diffraction integral. A number of cases have been considered. These include uniformly illuminated, obstructed circular apertures, and apertures coherently illuminated with a divergent gaussian beam wave. Other beam profiles that may be derived by energy redistribution of the laser beam have also been considered. The techniques used for the numerical evaluation of the Fresnel integral are discussed and results for the optical systems to be used on the High Altitude Aircraft Optical Communications Experiment are presented.

Webb, W. E.↗

Jupiter radiation belt models (July 1974)

Flux profiles which were derived from data returned by Pioneer 10 during Jupiter encounter, form the basis for a new set of numerical models for the energy spectra of electrons and protons in Jupiter's inner magnetosphere.

Divine, N.↗

OSO-7 observations of solar X-rays in the energy range 10-100 keV

Data on 123 hard X-ray bursts observed by the satellite OSO-7 between Oct. 10, 1971 and June 6, 1972 are described and evaluated. Typical duration of a burst is 100 sec. Average spectral indices lie between 3.5 and 5.5 for two-thirds of the 123 bursts, with a median of 4.6. In some events a soft-hard-soft pattern is observed, but there are numerous examples in which the spectrum softens continuously throughout the burst. The mean shape of the hard X-ray time profile as measured by the full width at half maximum does not depend on burst amplitude; nor does the spectral hardness correlate with the flux. The distribution of burst peak fluxes and the observation of large soft X-ray bursts without accompanying hard X rays suggest the existence of a distinct class of solar flares which emit only soft X rays. No center-to-limb variation was found in the frequency of occurrence of bursts or in the fraction with a nonthermal component. Estimates of the energy in the form of nonthermal electrons and in the flare plasma derived from these data indicate that the total amounts in each are comparable.

Datlowe, D. W.↗

Measurements of solar magnetic fields by Fourier transform techniques. II - Saturated and blended lines

Fourier techniques have been exhaustively calibrated using Unno's (1956) results for the absorption profile of a simple Zeeman triplet. If a simple transformation is applied to the normalized line depths, then magnetic-field strengths and inclination angles can be measured very accurately from noisy saturated line profiles. Systematic errors caused by saturation effects can be estimated and reduced by varying one parameter. When a significant fraction of the line profile is unsplit and unpolarized, large errors may be made in measurements of low fields, unless the line is sufficiently weak. For a weak line, a vertical field of 1600 gauss can be measured to 10% accuracy even when 70% of the line profile is stray light. These stray-light errors are troublesome in measuring fields of gaps and pores but not sunspots. Numerical results of the error analysis are presented graphically.

Tarbell, T. D.↗

Turbulent wake/shock interaction at Mach 6

The paper presents measurements of mean pitot, static pressure and total temperature profiles for a two-dimensional turbulent wake in the vicinity and downstream of intersection with a shock/expansion wave system. The wake centerline Mach number, velocity, and total temperature distributions indicated that the shock increased turbulent mixing, while the expansion field that followed the shock decreased mixing. The overall effect of the shock/expansion wave interaction was to increase mixing for shock angles 15 deg and less for the shock generator geometry used in the present study. Wake growth rates determined from the total temperature profiles confirmed the increased turbulent mixing effect of the shock. These data could be used to validate nonquilibrium turbulence modeling and numerical solution of the time-averaged Navier-Stokes equations.

Walsh, M. J.↗