Search NASA⌕ Search

SEARCH · Search NASA

Results for “EQUATIONS - DRAG, PROFILE - BICKNELL”

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 145 records · Page 8

Integral method for the calculation of three-dimensional, laminar and turbulent boundary layers

The method for turbulent flows is a further development of an existing method; profile families with two parameters and a lag entrainment method replace the simple entrainment method and power profiles with one parameter. The method for laminar flows is a new development. Moment of momentum equations were used for the solution of the problem, the profile families were derived from similar solutions of boundary layer equations. Laminar and turbulent flows at the wings were calculated. The influence of wing tapering on the boundary layer development was shown. The turbulent boundary layer for a revolution ellipsoid is calculated for 0 deg and 10 deg incidence angles.

Stock, H. W.↗

Computation of turbulent boundary layer flows with an algebraic stress turbulence model

An algebraic stress turbulence model is presented, characterized by the following: (1) the eddy viscosity expression is derived from the Reynolds stress turbulence model; (2) the turbulent kinetic energy dissipation rate equation is improved by including a production range time scale; and (3) the diffusion coefficients for turbulence equations are adjusted so that the kinetic energy profile extends further into the free stream region found in most experimental data. The turbulent flow equations were solved using a finite element method. Examples include: fully developed channel flow, fully developed pipe flow, flat plate boundary layer flow, plane jet exhausting into a moving stream, circular jet exhausting into a moving stream, and wall jet flow. Computational results compare favorably with experimental data for most of the examples considered. Significantly improved results were obtained for the plane jet flow, the circular jet flow, and the wall jet flow; whereas the remainder are comparable to those obtained by finite difference methods using the standard kappa-epsilon turbulence model. The latter seems to be promising with further improvement of the expression for the eddy viscosity coefficient.

Kim, Sang-Wook↗

Transverse Mode Dynamics of VCSELs Undergoing Current Modulation

Transverse mode dynamics of a 20-micron-diameter vertical-cavity surface-emitting laser (VCSEL) undergoing gain switching by deep current modulation is studied numerically. The direct current (dc) level is set slightly below threshold and is modulated by a large alternating current (ac). The resulting optical pulse train and transverse-mode patterns are obtained numerically. The ac frequency is varied from 2.5 GHz to 10 GHz, and the ac amplitude is varied from one-half to four times that of the dc level. At high modulation frequencies, a regular pulse train is not generated unless the ac amplitude is large enough. At all modulation frequencies, the transverse spatial profile switches from single-mode to multiple-mode pattern as the ac pumping level is increased. Optical pulse widths vary in the range 5-30 ps. with the pulse width decreasing when either the frequency is increased or the ac amplitude is decreased. The numerical modeling uses an approximation form of the semiconductor Maxwell-Bloch equations. Temporal evolution of the spatial profiles of the laser (and of carrier density) is determined without any assumptions about the type or number of modes. Keywords: VCSELs, current modulation, gain switching, transverse mode dynamics, computational modeling

Goorjian, Peter M.↗

Computational simulation of transition to turbulence through inverse modeling

The present investigation has focused on a computational methodology for the fundamental case of transition in channel flow, in which recently published experimental data are utilized both as a stimulus and as a measure of merit of the method. The research has proceeded along three avenues in parallel. The first task has consisted of the development and verification of a computer code which calculates the mean evolution of flow in a channel similar to the one employed experimentally by Blair and Anderson. An analytical test case was created for the dual purposes of code verification and of highlighting the interactions between the Reynolds stress and the mean velocity profile. This test case generated a Reynolds stress by the residue in the momentum equation which is produced by a typical analytical velocity profile. By a substitution of this Reynolds stress into the appropriate code module, the correctness of the code may be verified, along with the accuracy of the computational method. The second task pursued has involved the development of a triple layer model for the Reynolds stress profile, which was suggested and derived from experimental velocity profiles. It is demonstrated that the innermost length scale is based on the local friction velocity, the intermediate layer corresponds to the usual logarithmic law of the wall region in which the normalized Reynolds stress is approximately unity, and the outermost layer is represented by a closed mathematical form depending explicitly on the velocity profile in the wake region. The third task was comprised of scrutiny of the excellent databases developed by Blair and others, and the planning of its incorporation into the transition analysis. These extensive measurements indicate that turbulent statistics in the transition regime may be considered to alternate between laminar and fully turbulent types, the proportions of which are quantified by a measured intermittency function.

Sepri, Paavo↗

Structure of Weak Shock Waves in a Monatomic Gas

The profiles and thicknesses of normal shock waves in argon at Mach numbers of 1.335, 1.454, 1.576, and 1-713 were determined experimentally by means of a free-molecule probe whose equilibrium temperature is related by kinetic theory to the local flow properties and their gradients. Comparisons were made between the experimental shock profiles and the theoretical profiles calculated from the Navier-Stokes equations, the Grad 13-moment equations, and the Burnett equations. New, very accurate numerical integrations of the Burnett equations were obtained for this purpose with results quite different from those found by Zoller, to whom the solution of this problem is frequently attributed. The experimental shock profiles were predicted with approximately equal success by the Navier-Stokes and Burnett theories, while the 13-moment method was definitely less satisfactory. A surprising feature of the theoretical results is the relatively small difference in predictions between the Navier-Stokes and Burnett theories in the present range of shock strengths and the contrastingly large difference between predictions of Burnett and the 13-moment theories. It is concluded that the Navier-Stokes equations are correct for weak shocks and that within the present shock strength range the Burnett equations make no improvement which merits the trouble of solving them. For shocks of noticeably greater strength, say with a shock Mach number of more than 2.5, it remains fundamentally doubtful that any of these theories can be correct.

Sherman, F. S.↗

Electron density measurements of an inhomogeneous plasma using millimeter wave Fabry-Perot interferometers.

The electron density of a laboratory generated inhomogeneous plasma was measured using 60-90, GHz flat-plate and 150 GHz semi-confocal Fabry-Perot interferometers. The plasma was a negative glow-type and had a measured parabolic electron density profile. An analytical comparison between the derived equations of electron density for an assumed homogeneous plasma slab and an inhomogeneous plasma with parabolic distribution showed that the slab approximation was valid for such a profile. Comparison of measured values of electron density with those using X-band and K-band interferometers for the same plasma indicated good agreement.

Kendall, B. M.↗

An analysis of temperature effect in a finite journal bearing with spatial tilt and viscous dissipation

The analysis presented herein deals with the evaluation of the pressure, velocity, and temperature profiles in a finite-length plane journal bearing. The geometry of the case under study consists of a spatially tilted shaft. The two-dimensional Reynolds equation accounts for the variation of the clearance gap h with x and z and is used to model the pressure field. The latter is solved for a variety of shaft tilt angles and then used to calculate the two-dimensional flow field. Finally, the flow field is used in the energy equation to solve for the film temperature profile, when the effect of viscous dissipation is taken into account.

Braun, M. J.↗

Comparison of Burnett, super-Burnett and Monte Carlo solutions for hypersonic shock structure

The continuum Navier-Stokes, Burnett, and super-Burnett equations are solved for one-dimensional shock structures in various monoatomic gases. Solutions for a hard sphere gas, argon, and a Maxwellian gas from Mach 1.3 to Mach 50 are obtained. A numerical method utilizing the complete time-dependent continuum equations and obtaining the steady-state shock structure by allowing the system to relax from arbitrary initial conditions is employed. Shock density, velocity, temperature, and entropy profiles are also obtained using the direct simulation Monte Carlo method, and these results are used as bases for comparison for continuum solution profiles. It is shown that the Burnett equations yield shock structure solutions in much closer agreement to both Monte Carlo and experimental results than do the Navier-Stokes equations. Solutions to the super-Burnett equations with coefficients as presently derived, however, are considered to be inferior to those of the Burnett equations.

Fiscko, Kurt A.↗

Wall functions for the kappa-epsilon turbulence model in generalized nonorthogonal curvilinear coordinates

A k-epsilon turbulence model suitable for compressible flow, including the new wall function formulation, has been incorporated into an existing compressible Reynolds-averaged Navier-Stokes code, F3D. The low Reynolds number k-epsilon model of Chien (1982) was added for comparison with the present method. A number of features were added to the F3D code including improved far-field boundary conditions and viscous terms in the streamwise direction. A series of computations of increasing complexity was run to test the effectiveness of the new formulation. Flow over a flat plate was computed by using both orthogonal and nonorthogonal grids, and the friction coefficients and velocity profiles compared with a semi-empirical equation. Flow over a body of revolution at zero angle of attack was then computed to test the method's ability to handle flow over a curved surface. Friction coefficients and velocity profiles were compared to test data. All models gave good results on a relatively fine grid, but only the wall function formulation was effective with coarser grids. Finally, in order to demonstrate the method's ability to handle complex flow fields, separated flow over a prolate spheroid at angle of attack was computed, and results were compared to test data. The results were also compared to a k-epsilon model by Kim and Patel (1991), in which one equation model patched in at the wall was employed. Both models gave reasonable solutions, but improvement is required for accurate prediction of friction coefficients in the separated regions.

Sondak, D. L.↗

Effect of wall temperature on a supersonic turbulent boundary layer

Measurement of mean flow profiles in a fully developed Mach 3 turbulent boundary layer with negligible pressure gradient is reported. Data were acquired at several streamwise locations for wall-to-total temperature ratios of 0.94, 0.71, and 0.54. The results demonstrate that the velocity defect formulation of the law-of-the-wake, which successfully correlates compressible, adiabatic boundary layers, is also valid for nonadiabatic flows. It is also shown that for adiabatic walls, the linear Crocco relation between total temperature and velocity does not provide a valid test of the nature of the boundary-layer flow for practical cases where the Prandtl number departs from unity. Finally, the turbulent shear stress, mixing length, and eddy viscosity were extracted from the 'time-averaged' conservation equations using the measured mean flow profiles and found to be insensitive to wall temperature. In particular, the latter properties are in good agreement with earlier compressible, adiabatic correlations of turbulent transport properties.

Laderman, A. J.↗

Mapping the dark matter in the NGC 5044 group with ROSAT: Evidence for a nearly homogeneous cooling flow with a cooling wake

The NGC 5044 group of galaxies was observed by the ROSAT Position Sensitive Proportional Counter (PSPC) for 30 ks during its reduced pointed phase (1991 July). Due to the relatively cool gas temperature in the group (kT = 0.98 +/- 0.02 keV) and the excellent photon statistics (65,000 net counts), we are able to determine precisely a number of fundamental properties of the group within 250 kpc of the central galaxy. In particular, we present model-independent measurements of the total gravitating mass, the temperature and abundance profiles of the gas, and the mass accretion rate. Between 60 and 250 kpc, the gas is nearly isothermal with T varies as r(exp (-0.13 +/- 0.03)). The total gravitating mass of the group can be unambiguously determined from the observed density and temperature profiles of the gas using the equation of hydrostatic equilibrium. Within 250 kpc, the gravitating mass is 1.6 x 10(exp 13) solar mass, yielding a mass-to-light ratio of 130 solar mass/solar luminosity. The baryons (gas and stars) comprise 12% of the total mass within this radius. At small radii, the temperature clearly increases outward and attains a maximum value at 60 kpc. The positive temperature gradient in the center of the group confirms the existence of a cooling flow. The cooling flow region extends well beyond the temperature maximum with a cooling radius between 100 and 150 kpc. There are two distinct regions in the cooling flow separated by the temperature maximum. In the outer region, the gas is nearly isothermal with a unifor m Fe abundance of approximately 80% solar, the flow is nearly homogeneous with dot-M= 20 to 25 solar mass/year, the X-ray contours are spherically symmetric, and rho(sub gas) varies as r(exp -1.6). In the inner region, the temperature profile has a positive gradient, the mass accretion rate decreases rapidly inward, the gas density profile is steeper, and the X-ray image shows some substrucutre. NGC 5044 is offset from the centroid of the outer X-ray contours indicating that the central galaxy may have a residual velocity with respect to the center of the group potential. There is also a linear X-ray feature with an extent of approximately 30 kpc with one end coincident with NGC 5044. The X-ray emission from this feature is softer than the ambient gas. We interpret this feature as a 'cooling wake' formed by the accreting gas as it is gravitationally focused into the wake of NGC 5044. One of the most surprising results of our PSPC observation is the discovery of a nearly homogeneous cooling flow. Prior results concerning the mass accretion profile in cooling flows indicate that dot-M varies as r. This relation implies that significant mass deposition occurs at large radii which generates an inhomogeneous flow. The mass accretion rate in the NGC 5044 group is essentially a constant beyond 40 kpc (well within the cooling radius). Significant mass deposition (a declining dot-M) does not commence until the gas accretes to within 40 kpc of the group center where the radiative cooling time is approximately equals 10(exp 9) year. Th is radius also corresponds to the temperature maximum, the break in gas density profile, and the onset of structure in the X-ray image. A Hubble constant of H(sub 0) = 50 km/sec/Mpc is used throughout the paper.

David, Laurence P.↗

Time-dependent profiles of lava flows

This work investigates the combined effects of a time-dependent effusion rate and a spatially varying viscosity on the thickness profile of a flowing lava. Relatively simple governing equations are formulated as a primitive dynamical model for these two influences on lava flow morphology. Free boundary solutions for the profile of the flow and its extent are presented and analyzed for different boundary conditions at the source of the flow and several models for the spatial dependence of the viscosity. The results presented here suggest that time dependence in the flow depth at the source and the form of the viscosity variation are significant influences on the morphology and dimensions of lava flows. Moreover, this analysis implies that time-dependent source conditions may contribute to the scatter in plots of eruption rate versus flow length and disparities between field and laboratory estimates of lava viscosity.

Baloga, Stephen↗

Absorption line profiles in a moving atmosphere - A single scattering linear perturbation theory

An integral equation is derived which linearly relates Doppler perturbations in the spectrum of atmospheric absorption features to the wind system which creates them. The perturbation theory is developed using a single scattering model, which is validated against a multiple scattering calculation. The nature and basic properties of the kernels in the integral equation are examined. It is concluded that the kernels are well behaved and that wind velocity profiles can be recovered using standard inversion techniques.

Hays, P. B.↗

Estimates of equatorial upwelling between 140 deg and 110 deg W during 1984

The equation of continuity is used to estimate profiles of vertical velocity between 25 and 120 m from moored current measurements in arrays nested within the triangle with vertices located at 1 deg 30 arcmins S, 140 deg W and along the equator at 140 deg and 110 deg W during December 1983 through March 1984 and May-September 1984. All directions of the 4- to 5-month mean values were upward. The ensemble averaged mean vertical velocity was 0.000022 m/s. Upwelling speeds decreased eastward.

Halpern, David↗

Statistical versus nonstatistical temperature inversion methods

Vertical temperature profiles are derived from radiation measurements by inverting the integral equation of radiative transfer. Because of the nonuniqueness of the solution, the particular temperature profile obtained depends on the numerical inversion technique used and the type of auxiliary information incorporated in the solution. The choice of an inversion algorithm depends on many factors; including the speed and size of computer, the availability of representative statistics, and the accuracy of initial data. Results are presented for a numerical study comparing two contrasting inversion methods: the statistical-matrix inversion method and the nonstatistical-iterative method. These were found to be the most applicable to the problem of determining atmospheric temperature profiles. Tradeoffs between the two methods are discussed.

Smith, W. L.↗

Aerosol studies in mid-latitude coastal environments in Australia

The results of the evaluation of several inversion procedures that were used to select one which provides the most accurate atmospheric extinction profiles for small aerosol extinction coefficients (that often predominate in the maritime airmass) are presented. Height profiles of atmospheric extinction calculated by a two component atmospheric solution to the LIDAR equation will be compared with corresponding in-situ extinction profiles based on the size distribution profiles obtained in Western Australia. Values of the aerosol backscatter to extinction ratio obtained from multi-angle LIDAR measurements will be used in this solution.

Young, S. A.↗