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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
NASA NTRS1984

Boundary-layer linear stability theory

Most fluid flows are turbulent rather than laminar and the reason for this was studied. One of the earliest explanations was that laminar flow is unstable, and the linear instability theory was first developed to explore this possibility. A series of early papers by Rayleigh produced many notable results concerning the instability of inviscid flows, such as the discovery of inflectional instability. Viscosity was commonly thought to act only to stabilize the flow, and flows with convex velocity profiles appeared to be stable. The investigations that led to a viscous theory of boundary layer instability was reported. The earliest application of linear stability theory to transition prediction calculated the amplitude ratio of the most amplified frequency as a function of Reynolds number for a Blasius boundary layer, and found that this quantity had values between five and nine at the observed Ret. The experiment of Schubauer and Skramstad (1947) completely reversed the prevailing option and fully vindicated the Gottingen proponents of the theory. This experiment demonstrated the existence of instability waves in a boundary layer, their connection with transition, and the quantitative description of their behavior by the theory of Tollmien and Schlichting. It is generally accepted that flow parameters such as pressure gradient, suction and heat transfer qualitatively affect transition in the manner predicted by the linear theory, and in particular that a flow predicted to be stable by the theory should remain laminar. The linear theory, in the form of the e9, or N-factor is today in routine use in engineering studies of laminar flow. The stability theory to boundary layers with pressure gradients and suction was applied. The only large body of numerical results for exact boundary layer solutions before the advent of the computer age by calculating the stability characteristics of the Falkner-Skan family of velocity profiles are given. When the digital computer reached a stage of development which permit the direct solution of the primary differential equations, numerical results were obtained from the linear theory during the next 10 years for many different boundary layer flows: three dimensional boundary layers; free convention boundary layers; compressible boundary layers; boundary layers on compliant walls; a recomputation of Falkner-Skan flows; unsteady boundary layers; and heated wall boundary layers.

Leslie M Mack↗
NASA NTRS1985

Diurnal tide in the equatorial middle atmosphere as seen in LIMS temperatures

The distribution of day-night temperature differences in the middle atmosphere determined by the Nimbus 7 LIMS experiment is described. Day-night differences maximize at and are approximately symmetric about the equator. Successive centers of opposite sign increase in amplitude with altitude, the pattern having a vertical wavelength of approximately 25 km. Profiles of rocket meridional wind at Kwajalein (8.7 deg N) and Ascension Island (8.0 deg S) taken near local noon and averaged over the LIMS data period, exhibit maxima which support the tidal interpretation of the equatorial temperature pattern. These characteristics are in general agreement with previous observational and theoretical results for the solar driven diurnal tide. Substantial time variations in amplitude and in location of the temperature maxima are observed. The diurnal tide near the equatorial stratopause appears to be influenced by the phase of the semiannual oscillation.

Hitchman, M. H.↗
NASA NTRS1990

Distributions of He(+) at middle and equatorial latitudes during solar maximum

The properties of the plasma composition in the topside ionosphere between 15 and 35 deg magnetic latitude are examined. Calculations are used to show that the distribution of neutral species at solar maximum, together with the appropriate ionization rates, can readily account for the dominance of He(+) at midlatitudes. The relative abundance of the atmospheric species is a sensitive function of local time and the associated evolution of the topside O(+) concentration profile. The existence of an ExB drift motion of the plasma is needed to explain the He(+) minimum at the dip equator. While the He(+) concentration there is produced in one day against the chemical loss process, at midlatitudes a large flux tube volume provides a reservoir in which the He(+) can accumulate each day.

Heelis, R. A.↗
NASA NTRS1991

A comprehensive magnetohydrodynamic model of the Venus ionosphere

The MHD Venus ionospheric model of Shinegawa and Cravens (1988) is improved here by including the energy equations for ions and electrons in a self-consistent manner. This new model reproduces observed electron density and magnetic field profiles very well, while the basic MHD process of the Venus ionosphere remain virtually unchanged. The results indicate that including energetics does not significantly alter the density and magnetic field profiles. Under unmagnetized conditions, heat fluxes for both ions and electrons must be imposed to reproduce the observed plasma temperature profiles. A heat source for the ions is probably present at higher altitudes in the magnetized ionosphere. Heating processes do not play a significant role in the dynamics at low altitudes.

Shinagawa, H.↗
NASA NTRS1992

One Year of Rayleigh Lidar Measurements at Toronto

The Rayleigh lidar technique is being used at York University to study middle atmospheric thermal structure with particular emphasis on the properties and influence of internal gravity waves. These measurements are now carried out on a routine basis (every clear night) in order to determine seasonal variability and, eventually, long term trends. The transmitter consists of a doubled Nd:YAG laser with an output of 0.4 Joules per pulse at 20Hz (8 Watts) which is expanded to obtain a divergence of 0.4 milliradians. Backscattered photons are collected by a half meter cassegrain telescope, filtered to a bandwidth of 1nm, detected by a PMT, and counted with a multi-channel scaling averager. Strong backscatter from lower altitudes is blocked using a rotating chopper (25000rpm) which triggers the laser. The range gate and hardware averaging interval are variable but typically 300m and 5 minutes. Relative density profiles are obtained by simply normalizing to a model (eg. CIRA) at some altitude. Absolute temperature is derived using the hydrostatic equation, ideal gas law, and a model pressure value at the altitude where the photocount statistical uncertainty is about 15 percent - typically 90km for the entire night's average. The lower limit is currently about 30-35km due the presence aerosol from the volcanic eruption of Mt. Pinatubo.

Whiteway, James A.↗
NASA NTRS1993

Quasi-periodic transverse plasma flow associated with an evolving MHD vortex street in the outer heliosphere

We study a transverse plasma flow induced by the evolution of a Karman vortex street using a Chebyshev-Fourier spectral algorithm to solve both the compressible Navier-Stokes and MHD equations. The evolving vortex street is formed by the nonlinear interaction of two vortex sheets initially in equilibrium. We study spatial profiles of the total plasma velocity, the density, the meridional flow angle and the location of sector boundaries and find generally good agreement with Voyager 2 measurements of quasi-periodic transverse flow in the outer heliosphere. The pressure pulses associated with the meridional flows in the simulation are too small, although they are correctly located, and this may be due to the lack of any 'warp' in the current sheet in this model. A strong, flow-aligned magnetic field, such as would occur in the inner heliosphere, is shown to lead to weak effects that would be masked by the background interplanetary turbulence. We also study the plasma and magnetic transport resulting from the meridional flow and find that deficits of magnetic quantities do occur near the ecliptic. While the effect is relatively small, it is in general agreement with the most recent analysis of 'flux deficit' in the outer heliosphere.

Siregar, Edouard↗
NASA NTRS1994

X-Ray Emission from a Simulated Cluster of Galaxies

Using the 1993 cluster simulation of Katz & White, we analyze the intracluster medium and investigate the accuracy of the standard hydrostatic method for determining cluster masses. We show that the simulated cluster gas is in hydrostatic equilibrium with a subsonic flow toward the center. Inside a radius of (approx.) 100 kpc, this flow is in a steady state. The cooling time is shorter than a Hubble time within the central 50 kpc. The flow rate is regulated by the gas sink in the middle of the cluster and the PdV work done as the gas flows in, verifying the standard cooling flow scenario. We simulate observations of the cluster using the instrument parameters of the EXOSAT ME detector and the Einstein IPC detector. Even though the intracluster gas is not isothermal, isothermal models of the cluster, excluding regions within 100 kpc of galaxies, fit the EXOSAT X-ray spectra as well as they fit real clusters. The X- ray surface brightness distribution is similar to that of real clusters, again excluding the galaxies. We simulate the procedure used to determine the masses of real clusters. We use the equation of hydrostatic equilibrium together with the temperature derived from an isothermal fit to the simulated EXOSAT spectrum and the density profile derived from a fit to the simulated IPC surface brightness profile to determine the mass. A comparison of the derived mass profile to the actual mass profile shows that errors of a factor of 2 are possible. If the actual temperature profile is used, the cluster mass is found to an accuracy of better than 25% within the virial radius.

Tsai, John C.↗
NASA NTRS2001

Numerical Study of Buoyancy and Different Diffusion Effects on the Structure and Dynamics of Triple Flames

Numerical simulations of 2-D triple flames under gravity force have been implemented to identify the effects of gravity on triple flame structure and propagation properties and to understand the mechanisms of instabilities resulting from both heat release and buoyancy effects. A wide range of gravity conditions, heat release, and mixing widths for a scalar mixing layer are computed for downward-propagating (in the same direction with the gravity vector) and upward-propagating (in the opposite direction of the gravity vector) triple flames. Results of numerical simulations show that gravity strongly affects the triple flame speed through its contribution to the overall flow field. A simple analytical model for the triple flame speed, which accounts for both buoyancy and heat release, is developed. Comparisons of the proposed model with the numerical results for a wide range of gravity, heat release and mixing width conditions, yield very good agreement. The analysis shows that under neutral diffusion, downward propagation reduces the triple flame speed, while upward propagation enhances it. For the former condition, a critical Froude number may be evaluated, which corresponds to a vanishing triple flame speed. Downward-propagating triple flames at relatively strong gravity effects have exhibited instabilities. These instabilities are generated without any artificial forcing of the flow. Instead disturbances are initiated by minute round-off errors in the numerical simulations, and subsequently amplified by instabilities. A linear stability analysis on mean profiles of stable triple flame configurations have been performed to identify the most amplified frequency in spatially developed flows. The eigenfunction equations obtained from the linearized disturbance equations are solved using the shooting method. The linear stability analysis yields reasonably good agreements with the observed frequencies of the unstable triple flames. The frequencies and amplitudes of disturbances increase with the magnitude of the gravity vector. Moreover, disturbances appear to be most amplified just downstream of the premixed branches. The effects of mixing width and differential diffusion are investigated and their roles on the flame stability are studied.

Chen, Jyh-Yuan↗
NASA NTRS1998

Simulation of Ultra-Small MOSFETs Using a 2-D Quantum-Corrected Drift-Diffusion Model

We describe an electronic transport model and an implementation approach that respond to the challenges of device modeling for gigascale integration. We use the density-gradient (DG) transport model, which adds tunneling and quantum smoothing of carrier density profiles to the drift-diffusion model. We present the current implementation of the DG model in PROPHET, a partial differential equation solver developed by Lucent Technologies. This implementation approach permits rapid development and enhancement of models, as well as run-time modifications and model switching. We show that even in typical bulk transport devices such as P-N diodes and BJTs, DG quantum effects can significantly modify the I-V characteristics. Quantum effects are shown to be even more significant in small, surface transport devices, such as sub-0.1 micron MOSFETs. In thin-oxide MOS capacitors, we find that quantum effects may reduce gate capacitance by 25% or more. The inclusion of quantum effects in simulations dramatically improves the match between C-V simulations and measurements. Significant quantum corrections also occur in the I-V characteristics of short-channel MOSFETs due to the gate capacitance correction.

Biegel, Bryan A.↗
NASA NTRS2015

An Analytical Solution for Yaw Maneuver Optimization on the International Space Station and Other Orbiting Space Vehicles

This paper suggests a new method for optimizing yaw maneuvers on the International Space Station (ISS). Yaw rotations are the most common large maneuvers on the ISS often used for docking and undocking operations, as well as for other activities. When maneuver optimization is used, large maneuvers, which were performed on thrusters, could be performed either using control moment gyroscopes (CMG), or with significantly reduced thruster firings. Maneuver optimization helps to save expensive propellant and reduce structural loads ‐ an important factor for the ISS service life. In addition, optimized maneuvers reduce contamination of the critical elements of the vehicle structure, such as solar arrays. This paper presents an analytical solution for optimizing yaw attitude maneuvers. Equations describing pitch and roll motion needed to counteract the major torques during a yaw maneuver are obtained. A yaw rate profile is proposed. Also the paper describes the physical basis of the suggested optimization approach. In the obtained optimized case, the torques are significantly reduced. This torque reduction was compared to the existing optimization method which utilizes the computational solution. It was shown that the attitude profiles and the torque reduction have a good match for these two methods of optimization. The simulations using the ISS flight software showed similar propellant consumption for both methods. The analytical solution proposed in this paper has major benefits with respect to computational approach. In contrast to the current computational solution, which only can be calculated on the ground, the analytical solution does not require extensive computational resources, and can be implemented in the onboard software, thus, making the maneuver execution automatic. The automatic maneuver significantly simplifies the operations and, if necessary, allows to perform a maneuver without communication with the ground. It also reduces the probability of command errors. The suggested analytical solution provides a new method of maneuver optimization which is less complicated, automatic and more universal. A maneuver optimization approach, presented in this paper, can be used not only for the ISS, but for other orbiting space vehicles.

Dobrinskaya, Tatiana↗
NASA NTRS1999

Reusable Launch Vehicle Control in Multiple Time Scale Sliding Modes

A reusable launch vehicle control problem during ascent is addressed via multiple-time scaled continuous sliding mode control. The proposed sliding mode controller utilizes a two-loop structure and provides robust, de-coupled tracking of both orientation angle command profiles and angular rate command profiles in the presence of bounded external disturbances and plant uncertainties. Sliding mode control causes the angular rate and orientation angle tracking error dynamics to be constrained to linear, de-coupled, homogeneous, and vector valued differential equations with desired eigenvalues placement. The dual-time scale sliding mode controller was designed for the X-33 technology demonstration sub-orbital launch vehicle in the launch mode. 6DOF simulation results show that the designed controller provides robust, accurate, de-coupled tracking of the orientation angle command profiles in presence of external disturbances and vehicle inertia uncertainties. It creates possibility to operate the X-33 vehicle in an aircraft-like mode with reduced pre-launch adjustment of the control system.

Shtessel, Yuri↗
NASA NTRS2005

Navier-Stokes Computations With One-Equation Turbulence Model for Flows Along Concave Wall Surfaces

This report presents the use of a time-marching three-dimensional compressible Navier-Stokes equation numerical solver with a one-equation turbulence model to simulate the flow fields developed along concave wall surfaces without and with a downstream extension flat wall surface. The 3-D Navier- Stokes numerical solver came from the NASA Glenn-HT code. The one-equation turbulence model was derived from the Spalart and Allmaras model. The computational approach was first calibrated with the computations of the velocity and Reynolds shear stress profiles of a steady flat plate boundary layer flow. The computational approach was then used to simulate developing boundary layer flows along concave wall surfaces without and with a downstream extension wall. The author investigated the computational results of surface friction factors, near surface velocity components, near wall temperatures, and a turbulent shear stress component in terms of turbulence modeling, computational mesh configurations, inlet turbulence level, and time iteration step. The computational results were compared with existing measurements of skin friction factors, velocity components, and shear stresses of the developing boundary layer flows. With a fine computational mesh and a one-equation model, the computational approach could predict accurately the skin friction factors, near surface velocity and temperature, and shear stress within the flows. The computed velocity components and shear stresses also showed the vortices effect on the velocity variations over a concave wall. The computed eddy viscosities at the near wall locations were also compared with the results from a two equation turbulence modeling technique. The inlet turbulence length scale was found to have little effect on the eddy viscosities at locations near the concave wall surface. The eddy viscosities, from the one-equation and two-equation modeling, were comparable at most stream-wise stations. The present one-equation turbulence model is an effective approach for turbulence modeling in the near solid wall surface region of flow over a concave wall.

Wang, Chi R.↗
NASA NTRS1994

Molecular rotational line profiles from oxygen-rich red giant winds

We have developed a radiative transfer model of the dust and gas envelopes around late-type stars. The gas kinetic temperature for each star is calculated by solving equations of motion and the energy balance simultaneously. The main processes include viscous heating and adiabatic and radiative cooling. Heating is dominated by viscosity as the grains stream outward through the gas, with some contribution in oxygen-rich stars by near-IR pumping of H2O followed by collisional de-excitation in the inner envelope. For O-rich stars, rotational H2O cooling is a dominant mechanism in the middle part of the envelope, with CO cooling being less significant. We have applied our model to three well-studied oxygen-rich red giant stars. The three stars cover a wide range of mass-loss rates, and hence they have different temperature structures. The derived temperature structures are used in calculating CO line profiles for these objects. Comparison of the dust and gas mass-loss rates suggests that mass-loss rates are not constant during the asymptotic giant branch phase. In particular, the results show that the low CO 1-0 antenna temperatures of OH/IR stars reflect an earlier phase of much lower mass-loss rate.

Justtanont, K.↗
NASA NTRS1978

Effect of anomalous transport coefficients on the thermal structure of the storm time auroral ionosphere

By analyzing an observed storm time auroral electron temperature profile it is shown that anomalous transport effects strongly influence the thermal structure of the disturbed auroral ionosphere. Such anomalous transport effects are a consequence of plasma turbulence, the existence of which has been established by a large number of observations in the auroral ionosphere. The electron and composite ion energy equations are solved with anomalous electron thermal conductivity and parallel electrical resistivity coefficients. The solutions are parameterized with respect to a phenomenological altitude-dependent anomaly coefficient A and are compared with an observed storm time electron temperature profile above Chatanika. The calculated temperature profile for the classical case (A = 1) disagrees considerably with the measured profile over most of the altitude range up to 450 km. It is shown that an anomaly coefficient with a sharp peak of the order of 10,000 centered around the F2 peak is consistent with observations.

Fontheim, E. G.↗
DOE OSTI2024

Investigation of core transport changes in DIII-D discharges with off-axis T e profile peaks

Abstract DIII-D discharges that transition to H-mode solely with off-axis electron cyclotron heating (ECH) often exhibit strong off-axis peaking of electron temperature profiles at the heating location. Electron heat transport properties near these off-axis temperature peaks have been studied using modulated ECH. The Fourier analyzed electron temperature data have been used to infer electron thermal diffusivity. Comparisons with numerical solutions of the time-dependent electron thermal equation find that the data are consistent with a narrow region with electron diffusivity χ e an order of magnitude lower than the average value across the plasma, suggesting an electron internal transport barrier (ITB) near the ECH heating location. Detailed profile analysis and equilibrium reconstructions suggest that the formation of these ITBs are correlated with off-axis values of the safety factor q being near 1. Furthermore, the ECH driven H-mode discharges demonstrate more rapid electron heating rate near the ECH deposition location than L-mode discharges with higher auxiliary ECH heating power. Additional modeling attributes this difference to the modification of electron heat transport in the core at the L-H transition, which also sustains the off-axis electron temperature peaks.

Physics↗
DOE OSTI2026

Uncertainty quantification for Joule heating processes in fibrous pore-resolved media

Joule heating (JH) is an energy-efficient and sustainable technique for heating materials. Its application for industrial heating, particularly, has been gaining attention due to its potential for increasing the yield of various chemical products. The process involves the use of heating elements (materials that are highly conductive electrically and thermally) to heat up other materials or substances. These conductors, however, can exhbit varying degrees of uncertainty due to non-linearities in their temperature-dependent properties, which could result in variable material behavior. In this work, we carry out uncertainty quantification (UQ) at the pore scale to describe the uncertainty of such materials. In so doing, we applied the non-intrusive polynomial chaos expansion (PCE) technique to quantify the uncertainty within the system. The steady state Joule heating equation was solved numerically at the pore scale mimicking conditions within a heating chamber for propane dehydrogenation, and various electro-thermal profiles were obtained. We also examined the effect of the number of sampling points (20 – 100) and order of the PCE coefficients (2 – 5) on the accuracy of the temperature evaluations. The results were then benchmarked with the standard Monte Carlo (MC) method. The average temperature of the 4th-order global PCE showed good agreement with the MC results (which were positively skewed). Orders greater than 4 gave an underestimation of the temperatures while predictions for the peak temperature improved as the number of sampling points increased.

Fagbemi, Samuel [ORNL] (ORCID:0000000236995025)↗
DOE OSTI2024

Simultaneous measurement of surface velocity and plasma density with interferometric velocimetry

The apparent velocity measured by an interferometric surface velocimeter is a function of both the surface velocity and the time derivative of the refractive index along the measurement path. We employed this dual sensitivity to simultaneously measure km/s surface velocities and 10 18 cm -3 average plasma densities with combined VISAR (velocity interferometer system for any reflector) and PDV (photonic Doppler velocimetry) measurements in experiments performed on the Z Pulsed Power Facility. We detail the governing equations, associated assumptions, and analysis specifics and show that the surface velocity can be extracted without knowledge of the specific plasma density profile.

47 OTHER INSTRUMENTATION↗
DOE OSTI2025

Atmospheric Profile Builder

APBuilder is a software tool designed to generate atmospheric profiles for use with AC2Dr by downloading weather model data and transforming it into 1D or 2D binary profiles. These binary profiles serve as required inputs for AC2Dr, an LLNL-developed, open source, two-dimensional numerical solver for the acoustic wave equation based on the finite difference method.

Kim, Keehoon [Lawrence Livermore National Laborato↗