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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Modeling Heat Transfer and Reaction Kinetics of Biomass in Pyrolysis Feeding Systems

Biomass feedstocks can be used to produce low carbon-intensity chemicals and liquid fuels by displacing fossil-fuel sources. Pyrolysis is the thermal decomposition of biomass in the absence of oxygen and can be used to produce crude biofuels. A typical feeding mechanism for pyrolysis reactors is by lock hopper followed by a horizontal auger feeder. Particle agglomeration and plugging of the auger screw are regularly occurring phenomena that is detrimental to the process since it results in the need for periodic inspection and cleaning. The underlying mechanisms leading to agglomeration and plugging are unclear, but a working hypothesis is that heat from the reactor raises the temperature of the biomass in the feeder to the point at which preliminary decomposition reactions occur and produce "sticky" products. Using principles of heat and mass transfer, we derived one-dimensional differential equations for heat flow from the pyrolysis reactor through the auger feeder. Temperature profiles were solved from the system of one-dimensional differential equations using the shooting method. Subsequently, we used the temperature profiles as a one-way coupling to a reaction model for the biomass in the feeder. Kinetic models were adopted from the literature and used to predict the formation of intermediates, products, char, extractives, and metaplastics as well as phase changes along the auger feeder. Species were identified that correlated with softening of the biomass, causing it to become sticky and plug the feeder. We will discuss the details of this work, our findings, and the direction of future work that includes fully coupled transport and reaction-kinetics simulations.

BIOMASS FUELS↗
NASA NTRS1976

A parametric analysis of visual approaches for helicopters

A flight investigation was conducted to determine the characteristic shapes of the altitude, ground speed, and deceleration profiles of visual approaches for helicopters. Two hundred thirty-six visual approaches were flown from nine sets of initial conditions with four types of helicopters. Mathematical relationships were developed that describe the characteristic visual deceleration profiles. These mathematical relationships were expanded to develop equations which define the corresponding nominal ground speed, pitch attitude, pitch rate, and pitch acceleration profiles. Results are applicable to improved helicopter handling qualities in terminal area operations.

Moen, G. C.↗
NASA NTRS1980

Variational description of the positive column with two-step ionization

The ionization balance in diffusion-dominated discharges may depend on both one- and two-step ionization processes. The Spenke diffusion equation describing such conditions is solved in this paper by the Rayleigh-Ritz variational method. Simple analytic approximations to the density profile, and the similarity relation between the equation parameter, and the discharge dimensions, are derived for planar and cylindrical geometry and compared with exact computations for certain limiting cases.

Crawford, F. W.↗
NASA NTRS1985

MAGSAT satellite magnetic anomaly map over South America

A scalar magnetic anomaly map was prepared for South America and adjacent marine areas directly from original MAGSAT orbits. Special problems associated with the separation of external field and crustal anomalies, and the reduction of data to a common altitude are addressed. External fields are manifested in a long-wavelength ring current effect, a medium-wavelength equatorial electrojet, and short-wavelength noise. The noise is reduced by selecting profiles from quiet periods (Kp or = 3), and the effect of the electrojet is minimized by drawing the data set from dawn profiles only. The ring current is corrected through the use of a standard equation, augmented by further digital band-pass filtering. Profiles thus filtered differ primarily in amplitude due solely to satellite altitude differences. These differences are normalized by an inversion of the profile data onto a grid of equivalent point dipoles, and recalculated at an altitude of 350 km. The low altitudes in the study area cause instability in the inversion, necessitating separate inversions of several sub-areas which are subsequently merged. Crustal anomalies reduced-to-the-pole exhibit marked correlations to known tectonic features.

Ridgway, J. R.↗
NASA NTRS1978

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.↗
NASA NTRS1986

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↗
NASA NTRS2000

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.↗
NASA NTRS1989

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↗
NASA NTRS1959

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.↗
NASA NTRS1972

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.↗
NASA NTRS1984

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.↗
NASA NTRS1989

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.↗
NASA NTRS1992

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.↗
NASA NTRS1978

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.↗
NASA NTRS1994

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.↗
NASA NTRS1986

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↗