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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

A New Hybrid Mass-Flux/High-Order Turbulence Closure for Ocean Vertical Mixing

While various parameterizations of vertical turbulent fluxes at different levels of complexity have been proposed, each has its own limitations. For example, simple first-order closure schemes such as the K-Profile Parameterization (KPP) lack energetic constraints; two-equation models $k-ε$ like directly solve an equation for the turbulent kinetic energy but do not account for non-diffusive fluxes, and high-order closures that include the high-order transport terms are computationally expensive. To address these, we extend the Assumed-Distribution Higher-Order Closure (ADC) framework originally proposed for the atmospheric boundary layer and apply it to the ocean surface boundary layer. By assuming a probability distribution function relationship between the vertical velocity and tracers, all second-order and higher-order moments are exactly constructed and turbulence closure is achieved in the ADC scheme. In addition, this ADC parameterization has full energetic constraints and includes non-diffusive fluxes without the computational cost of a full higher-order closure scheme. We have tested the ADC scheme against a combination of large eddy simulation (LES), KPP, and $k-ε$ for surface buoyancy-driven convective mixing and found that the ADC scheme is robust with different vertical resolutions and compares well to the LES results.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Explosive dispersal of particles in high speed environments

In this paper, we present the results of the explosive dispersal of particles in high-speed environments. We carry out Euler–Lagrange numerical simulations of a source at quiescent ambient conditions as well as moving at Mach numbers of 3 and 6. Particle volume fractions of 0%, 1%, and 4.5% are presented. The detonation profile is computed with the Jones–Wilkins–Lee equation of state using a reactive burn model. Non-static cases provide a framework to consider the effect of a bow shock and pre-existing high-speed flow conditions on the dispersal process. We also compute averages of both static and dynamic pressures, as well as impulse density histories on virtual probe planes to characterize the momentum of the flow and particles that would deposit on a target. Results suggest that the presence of the particles can have a substantial effect on the pressure average of the virtual target planes.

Durant, B. A.↗
NASA NTRS1974

Westward equatorial electrojet during daytime hours

The phenomenon of the depression of the geomagnetic horizontal field during the daytime hours of magnetically quiet days at equatorial stations is described. These events are generally seen around 0700 and 1600 LT, being more frequent during the evening than the morning hours. The evening events are more frequent during periods of low solar activity and in the longitude region of weak equatorial electrojet currents. The latitudinal extent of the phenomenon is limited to the normal equatorial electrojet region, and on some occasions the phenomenon is not seen at both stations, separated by only a few hours in longitude. During such an event, the latitudinal profile of the geomagnetic vertical field across the equator is reversed, the ionospheric drift near the equator is reversed toward the east, the q type of sporadic E layer is completely absent, and the height of the peak ionization in the F2 region is decreased. It is suggested that these effects are caused by a narrow band of current flowing westward in the E region of the ionosphere and within the latitude region of the normal equatorial electrojet, due to the reversal of the east-west electrostatic field at low latitudes.

Rastogi, R. G.↗
NASA NTRS1974

Design and construction of a low-velocity boundary-layer temperature probe

This paper presents the results of a study performed to determine the probe geometry and thermal characteristics for a low-velocity boundary-layer temperature probe capable of eliminating conduction errors. Design curves for four common thermoelectric pairs are presented, and an experimental turbulent boundary-layer temperature profile is compared with the theoretical laminar sublayer equation to verify that the probe eliminates conduction errors.

Blackwell, B. F.↗
NASA NTRS1976

The proper weighting function for retrieving temperatures from satellite measured radiances

One class of methods for converting satellite measured radiances into atmospheric temperature profiles, involves a linearization of the radiative transfer equation: delta r = the sum of (W sub i) (delta T sub i) where (i=1...s) and where delta T sub i is the deviation of the temperature in layer i from that of a reference atmosphere, delta R is the difference in the radiance at satellite altitude from the corresponding radiance for the reference atmosphere, and W sub i is the discrete (or vector) form of the T-weighting (i.e., temperature weighting) function W(P), where P is pressure. The top layer of the atmosphere corresponds to i = 1, the bottom layer to i = s - 1, and i = s refers to the surface. Linearization in temperature (or some function of temperature) is at the heart of all linear or matrix methods. The weighting function that should be used is developed.

Arking, A.↗
NASA NTRS1978

Measurements in the turbulent boundary layer at constant pressure in subsonic and supersonic flow. Part 1: Mean flow

Experiments were carried out to test the accuracy of laser Doppler instrumentation for measurement of Reynolds stresses in turbulent boundary layers in supersonic flow. Two facilities were used to study flow at constant pressure. In one facility, data were obtained on a flat plate at M sub e = 0.1, with Re theta up to 8,000. In the other, data were obtained on an adiabatic nozzle wall at M sub e = 0.6, 0.8, 1.0, 1.3, and 2.2, with Re theta = 23,000 and 40,000. The mean flow as observed using Pitot tube, Preston tube, and floating element instrumentation is described. Emphasis is on the use of similarity laws with Van Driest scaling and on the inference of the shearing stress profile and the normal velocity component from the equations of mean motion. The experimental data are tabulated.

Collins, D. J.↗
NASA NTRS1984

Tension waves in tethered satellite cables

A one-degree-of-freedom simulation of the Tethered Satellite System (TSS) was programmed using a distributed system model of the tether based on the one-dimensional wave equation. This model represents the time varying tension profile along the tether as the sum of two traveling waves of tension moving in opposite directions. A control loop was devised which combines a deployment rate command with the measured tension at the deployer to produce a smooth, stable rate of deployment of the subsatellite. Simulation results show a buildup of periodic bursts of high frequency oscillation in tension. This report covers the mathematical modelling and simulation results and explains the reason for the observed oscillations. The design of a possible vibration damping device is discussed.

Lallman, F. J.↗
NASA NTRS1985

Surface Tensions and Their Variations with Temperature and Impurities

The surface tensions in this work were determined using the sessile drop technique. This method is based on a comparison of the profile of a liquid drop with the profile calculated by solving the Young-Laplace equation. The comparison can be made in several ways; the traditional Bashforth-Adams procedure was used in conjunction with recently calculated drop shape tables which virtually eliminate interpolation errors. Although previous study has found little difference in measurements with pure and oxygen doped silicon, there is other evidence suggesting that oxygen in dilute concentrations severely depresses the surface tension of silicon. The surface tension of liquid silicon in purified argon atmospheres was measured. A temperature coefficient near -0.28 mJ/square meters K was found. The experiments show a high sensitivity of the surface tension to what is believed are low concentrations of oxygen. Thus one cannot rule out some effect of low levels of oxygen in the results. However, the highest surface tension values obtained in conditions which minimized the residual oxygen pressure are in good agreement with a previous measurement in pure hydrogen. Therefore, depression of the surface tension by oxygen is insignificant in these measurements.

Hardy, S. C.↗
NASA NTRS1987

Removal of spurious modes encountered in solving stability problems by spectral methods

A technique based on the Galerkin approximation is developed to remove spurious roots arising when Chebyshev spectral methods are used to solve eigenvalue problems in hydrodynamic stability. The derivation of Galerkin-Chebyshev approximations is explained, and numerical results for the Orr-Sommerfeld equations of plane Poiseuille flow and a Blasius profile are presented in tables and compared with those obtained by the method of Zebib (1984). It is pointed out that the present method does not increase the size of the algebraic system to be solved.

Zebib, Abdelfattah↗
NASA NTRS1988

Semi-empirical models of the wind in cool supergiant stars

A self-consistent semi-empirical model for the wind of the supergiant in zeta Aurigae type systems is proposed. The damping of the Alfven waves which are assumed to drive the wind is derived from the observed velocity profile. Solution of the ionization balance and energy equation gives the temperature structure for given stellar magnetic field and wave flux. Physically acceptable solutions of the temperature structure place limits on the stellar magnetic field. A crude formula for a critical mass loss rate is derived. For a mass loss rate below the critical value the wind cannot be cool. Comparison between the observed and the critical mass loss rate suggests that the proposed theory may provide an explanation for the coronal dividing line in the Hertzsprung-Russell diagram. The physical explanation may be that the atmosphere has a cool wind, unless it is physically impossible to have one. Stars which cannot have a cool wind release their nonthermal energy in an outer atmosphere at coronal temperatures. It is possible that in the absence of a substantial stellar wind the magnetic field has less incentive to extend radially outward, and coronal loop structures may become more dominant.

Kuin, N. P. M.↗
NASA NTRS1989

Kinetic theory model for the flow of a simple gas from a two-dimensional nozzle

A system of nonlinear integral equations equivalent to the Krook kinetic equation for the steady state is the mathematical basis used to develop a computer code to model the flowfields for low-thrust two-dimensional nozzles. The method of characteristics was used to solve numerically by an iteration process the approximated Boltzmann equation for the number density, temperature, and velocity profiles of a simple gas as it exhausts into a vacuum. Results predict backscatter and show the effect of the inside wall boundary layer on the flowfields external to the nozzle.

Riley, B. R.↗
NASA NTRS1992

Lidar multiple scattering: Dependence on atmospheric parameters

Backscatter lidars are very useful tools in gaining information about the atmosphere. Inversion of backscatter signals can be used to retrieve the spatial distribution of clouds and aerosols, e.g., vertical extensions of clouds and aerosol layers, detection of 'subvisible' cirrus, etc., very accurately compared to passive radiometers. To retrieve extinction coefficient profiles, in general, the so-called lidar equation is solved. This equation describes a laser pulse which propagates through an atmosphere with transmission exp(- integral of sigma(z')dz'), which is scattered singly under 180 degrees towards a receiver, and which is attenuated on its way back again by exp(- integral of sigma(z')dz'). For inversion, further information on optical properties of the scatterers are required (extinction to backscatter ratio). However, it is evident that, at least for optically thick clouds and/or large receiver's field-of-views, this approach may be erroneous because of multiply scattered photons. Several studies were performed to simulate multiple scattering in case of collimated beams or pulses including Monte Carlo, transport, and diffusion theory. Unfortunately, most of the approaches require an enormous amount of computer time. Thus, it is meaningful to elaborate approximate methods, and one such method is discussed. Our method is based on Bissonnett (1988). We have improved the numerical simulation of limited field-of-view detectors, and compared the results to Monte Carlo calculations and an exact model of Altmann (1989).

Wiegner, Matthias↗
NASA NTRS1992

A comparison of the reduced and approximate systems for the time dependent computation of the polar wind and multiconstituent stellar winds

The paper derives the 'reduced' system of equations commonly used to describe the time evolution of the polar wind and multiconstituent stellar winds from the equations for a multispecies plasma with known temperature profiles by assuming that the electron thermal speed approaches infinity. The reduced system is proved to have unbounded growth near the sonic point of the protons for many of the standard parameter cases. For the same parameter cases, the unmodified system exhibits growth in some of the Fourier modes, but this growth is bounded. An alternate system (the 'approximate' system) in which the electron thermal speed is slowed down is introduced. The approximate system retains the mathematical behavior of the unmodified system and can be shown to accurately describe the smooth solutions of the unmodified system. Other advantages of the approximate system over the reduced system are discussed.

Browning, G. L.↗
NASA NTRS1991

Kinetic theory model for the flow of a simple gas from a three-dimensional axisymmetric nozzle

A system of nonlinear integral equations equivalent to the Krook kinetic equations for the steady state is the mathematical basis used to develop a computer code to model the flowfields for low-thrust three-dimensional axisymmetric nozzles. The method of characteristics is used to solve numerically by an iteration process the approximated Boltzmann equation for the number density, temperature, and velocity profiles of a simple gas as it expands into a vacuum. Results predict backscatter and show the effect of the nozzle wall boundary layer on the external flowfields.

Riley, B. R.↗
NASA NTRS1996

Symmetries in Turbulent Boundary Layer Flows

The objective is the development of a new theory which enables the algorithmic computation of all self-similar mean velocity profiles. The theory is based on Liegroup analysis and unifies a large set of self-similar solutions for the mean velocity of stationary parallel turbulent shear flows. The results include the logarithmic law of the wall, an algebraic law, the viscous sublayer, the linear region in the middle of a Couette flow and in the middle of a rotating channel flow, and a new exponential mean velocity profile not previously reported. Experimental results taken in the outer parts of a high Reynolds number flat-plate boundary layer, strongly support the exponential profile. From experimental as well as from DNS data of a turbulent channel flow the algebraic scaling law could be confirmed in both the center region and in the near wall region. In the case of the logarithmic law of the wall, the scaling with the wall distance arises as a result of the analysis and has not been assumed in the derivation. The crucial part of the derivation of all the different mean velocity profiles is to consider the invariance of the equation for the velocity fluctuations at the same time as the invariance of the equation for the velocity product equations. The latter is the dyad product of the velocity fluctuations with the equation for the velocity fluctuations. It has been proven that all the invariant solutions are also consistent with similarity of all velocity moment equations up to any arbitrary order.

Oberlack, M.↗
NASA NTRS2004

Chemical models for simulating single-walled nanotube production in arc vaporization and laser ablation processes

Chemical kinetic models for the nucleation and growth of clusters and single-walled carbon nanotube (SWNT) growth are developed for numerical simulations of the production of SWNTs. Two models that involve evaporation and condensation of carbon and metal catalysts, a full model involving all carbon clusters up to C80, and a reduced model are discussed. The full model is based on a fullerene model, but nickel and carbon/nickel cluster reactions are added to form SWNTs from soot and fullerenes. The full model has a large number of species--so large that to incorporate them into a flow field computation for simulating laser ablation and arc processes requires that they be simplified. The model is reduced by defining large clusters that represent many various sized clusters. Comparisons are given between these models for cases that may be applicable to arc and laser ablation production. Solutions to the system of chemical rate equations of these models for a ramped temperature profile show that production of various species, including SWNTs, agree to within about 50% for a fast ramp, and within 10% for a slower temperature decay time.

Nanotubes, Carbon/chemistry↗