Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “domination number”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

On the Behavior of Three-dimensional Wave Packets in Viscously Spreading Mixing Layers

We consider analytically the evolution of a three-dimensional wave packet generated by an impulsive source in a mixing layer. The base flow is assumed to be spreading due to viscous diffusion. The analysis is restricted to small disturbances (linearized theory). A suitable high-frequency ansatz is used to describe the packet; the key elements of this description are a complex phase and a wave action density. It is found that the product of this density and an infinitesimal material volume convecting at the local group velocity is not conserved: there is a continuous interaction between the base flow and the wave action. This interaction is determined by suitable mode-weighted averages of the second and fourth derivatives of the base-flow velocity profile. Although there is some tendency for the dominant wave number in the packet to shift from the most unstable value toward the neutral value, this shift is quite moderate. In practice, wave packets do not become locally neutral in a diverging base flow (as do instability modes), therefore, they are expected to grow more suddenly than pure instability modes and do not develop critical layers. The group velocity is complex; the full significance of this is realized by analytically continuing the equations for the phase and wave action into a complex domain. The implications of this analytic continuation are discussed vis-a-vis the secondary instabilities of the packet: very small scale perturbations on the phase can grow very rapidly initially, but saturate later because most of the energy in these perturbations is convected away by the group velocity. This remark, as well as the one regarding critical layers, has consequences for the nonlinear theories.

Balsa, Thomas F.↗
NASA NTRS2000

Convective Melting of Particles in Flow Under Microgravity Conditions

Study of melting of dispersed or packed solid particles in a fluid under gravity and microgravity conditions provides benchmark information for many engineering applications such as material processing, environmental assessment and protection and space fire protection. During such processes, packed or dispersed solid particles are interacting with fluid flow at above-melting temperatures. By unmasking the buoyancy effects in coupled flow, phase change and heat transport phenomena, a better understanding of melting rate in non-thermal equilibrium, convective conditions can be studied. A series of flight experiments were conducted onboard the NASA KC-135 Microgravity Research Airplane. The Particle Melting in Plow (PMF) module was designed to allow flow through the initially packed ice particles at controlled temperature and velocity. To achieve this, a close-loop flow system was designed. Video images were taken to record the visualization of the melting process, from which a time variation of packed particle thickness distribution at different times can be obtained by the image analysis method. The fluid temperature distribution within the melting zone is measured by thermocouples. An infrared camera was mounted from the top of the test section to record the ice-water thermal images at a given location. The results from thermal images yield local temperature variation between melting solid and liquid and local Stephan number. Typical results for a number of cases are presented. The mathematical model, describing mass, energy, and momentum balance equations for the liquid and solid phases, is presented. It is found that melting rate is influenced mainly by the ratio of Reynolds number (based on the initial particle diameter) to the Froud number, and Stephan number. At the absence of gravity, Froud number approaches zero, Reynolds number and Stephan number become dominant factors governing the melting rate. The numerically determined results are compared with the experimental ones. It is found that the discrepancy between the predicted and measured melting rate is largely due to the inaccuracy in the constitutive equations for effective thermal physical properties, such as effective thermal conductivity and diffusivity, and transport properties, such as particle interaction coefficient, and local heat transfer coefficient of particles.

Stewart, C.↗
NASA NTRS2005

Analysis of Hydrodynamics and Heat Transfer in a Thin Liquid Film Flowing over a Rotating Disk by Integral Method

An integral analysis of hydrodynamics and heat transfer in a thin liquid film flowing over a rotating disk surface is presented for both constant temperature and constant heat flux boundary conditions. The model is found to capture the correct trends of the liquid film thickness variation over the disk surface and compare reasonably well with experimental results over the range of Reynolds and Rossby numbers covering both inertia and rotation dominated regimes. Nusselt number variation over the disk surface shows two types of behavior. At low rotation rates, the Nusselt number exhibits a radial decay with Nusselt number magnitudes increasing with higher inlet Reynolds number for both constant wall temperature and heat flux cases. At high rotation rates, the Nusselt number profiles exhibit a peak whose location advances radially outward with increasing film Reynolds number or inertia. The results also compare favorably with the full numerical simulation results from an earlier study as well as with the reported experimental results.

Basu, S.↗
NASA NTRS2021

Experimental and Computational Investigation of Flow Boiling in Microgravity

This study explores use of Computational Fluid Dynamics (CFD) to predict near-saturated flow boiling of FC-72 in microgravity. The computational method employs transient analysis to predict detailed interfacial behavior and heat transfer characteristics along a rectangular channel heated along two opposite walls. Predicted results are validated against experimental temperature measurements and high-speed video images captured during a series of parabolic aircraft maneuvers for three sets of operating conditions which include variations of both mass velocity and wall heat flux. The computational method is based on the multi-phase volume of fluid (VOF) model, which is combined with appropriate phase change and turbulence models, and accounts for both shear-lift force on bubbles and conjugate heat transfer along the heating walls. A key advantage of the CFD method is ability to capture details that are very difficult to measure experimentally, including detailed spatial variations of bubble shape, void fraction, mixture fluid temperature, liquid velocity, and vapor velocity, results for which are presented for each of the three test cases. Different flow regimes predicted along the heated length exhibit a number of dominant mechanisms including bubble nucleation, bubble growth, coalescence, vapor blankets, interfacial waviness, and residual liquid sub-layer, all of which agree well with experiment. Vapor velocity is shown to increase appreciably along the heated length because of increased void fraction, while liquid velocity experiences large fluctuations. Non-equilibrium effects are accentuated with increasing mass velocity, contributing minor deviations of fluid temperature from simulations compared to those predicted by the analytical method. Predicted wall temperature is fairly uniform in the middle of the heated length but increases in the entrance region, due to sensible heat transfer in the subcooled liquid, and decreases toward the exit, mostly because of flow acceleration resulting from increased void fraction.

two-phase cooling; flow boiling; CFD; microgravity↗
NASA NTRS1986

Evidence of a saturated gravity-wave spectrum throughout the atmosphere

The view adapted here is that the dominant mesoscale motions are due to internal gravity waves and show that previous and new vertical wave number spectra of horizontal winds are consistent with the notion of a saturation limit on wave amplitudes. It is also proposed that, at any height, only those vertical wave numbers m less than m sub asterisk are at saturation amplitudes, where m sub asterisk is the vertical wave number of the dominant energy-containing scale. Wave numbers m less than m sub asterisk are unsaturated, but experience growth with height due to the decrease of atmospheric density. The result is a saturated spectrum of gravity waves with both m sub asterisk decreasing and wave energy increasing with height. This saturation theory is consistent with a variety of atmospheric spectral observations and provides a basis for the notion of a universal spectrum of atmospheric gravity waves.

Fritts, D. C.↗
NASA NTRS1975

Estimation of transonic aircraft aerodynamics to high angles of attack

A recently developed method for estimating transonic aircraft aerodynamics to high angles of attack which is especially suited to computerized aircraft design and optimization studies is described. Following theoretical and empirical guidelines, explicit nonlinear equations are formulated for subsonic, transonic, and supersonic speeds; new algorithms are used for assessing compressibility effects and for analyzing transonic, shock-dominated flows which adhere to Laitone's limit Mach-number criterion. Viscous dominated flows are not analyzed. The transonic influence of viscosity is indirectly accommodated by assigning designated inputs rather than extracted solutions for the chordwise locations of shock and separation. The method is extended to complete aircraft configurations by including accounts for nose lift, the wing downwash field, and tail lift. Several comparisons of experiment and estimate are included.

Axelson, J. A.↗
NASA NTRS2001

Warm-Core Intensification Through Horizontal Eddy Heat Transports into the Eye

A simulation of Hurricane Bob (1991) using the PSU/NCAR MM5 mesoscale model with a finest mesh spacing of 1.3 km is used to diagnose the heat budget of the hurricane. Heat budget terms, including latent and radiative heating, boundary layer forcing, and advection terms were output directly from the model for a 6-h period with 2-min frequency. Previous studies of warm core formation have emphasized the warming associated with gentle subsidence within the eye. The simulation of Hurricane Bob confirms subsidence warming as a major factor for eye warming, but also shows a significant contribution from horizontal advective terms. When averaged over the area of the eye, subsidence is found to strongly warm the mid-troposphere (2-9 km) while horizontal advection warms the mid to upper troposphere (5-13 km) with about equal magnitude. Partitioning of the horizontal advective terms into azimuthal mean and eddy components shows that the mean radial circulation does not, as expected, generally contribute to this warming, but that it is produced almost entirely by the horizontal eddy transport of heat into the eye. A further breakdown of the eddy components into azimuthal wave numbers 1, 2, and higher indicates that the warming is dominated by wave number 1 asymmetries, with smaller coming from higher wave numbers. Warming by horizontal eddy transport is consistent with idealized modeling of vortex Rossby waves and work is in progress to identify and clarify the role of vortex Rossby waves in warm-core intensification in both the full-physics model and idealized models.

Braun, Scott A.↗
NASA NTRS2001

Warm-Core Intensification of a Hurricane Through Horizontal Eddy Heat Transports Inside the Eye

A simulation of Hurricane Bob (1991) using the PSU/NCAR MM5 mesoscale model with a finest mesh spacing of 1.3 km is used to diagnose the heat budget of the hurricane. Heat budget terms, including latent and radiative heating, boundary layer forcing, and advection terms were output directly from the model for a 6-h period with 2-min frequency. Previous studies of warm core formation have emphasized the warming associated with gentle subsidence within the eye. The simulation of Hurricane Bob also identifies subsidence warming as a major factor for eye warming, but also shows a significant contribution from horizontal advective terms. When averaged over the area of the eye, excluding the eyewall (at least in an azimuthal mean sense), subsidence is found to strongly warm the mid-troposphere (2-9 km) while horizontal advection warms the mid to upper troposphere (5-13 km) with about equal magnitude. Partitioning of the horizontal advective terms into azimuthal mean and eddy components shows that the mean radial circulation cannot, as expected, generally contribute to this warming, but that it is produced almost entirely by the horizontal eddy transport of heat into the eye. A further breakdown of the eddy components into azimuthal wave numbers 1, 2, and higher indicates that the warming is dominated by wave number 1 asymmetries, with smaller contributions coming from higher wave numbers. Warming by horizontal eddy transport is consistent with idealized modeling of vortex Rossby waves and work is in progress to identify and clarify the role of vortex Rossby waves in warm-core intensification in both the full-physics model and idealized models.

Braun, Scott A.↗
NASA NTRS2010

Characteristics of the Dust-Plasma Interaction Near Enceladus' South Pole

We present RPWS Langmuir probe data from the third Enceladus flyby (E3) showing (he presence of dusty plasma near Enceladus' South Pole. There is a sharp rise in both the electron and ion number densities when the spacecraft traverses through Enceladus plume. The ion density near Enceladus is found to increase abruptly from about 10(exp 2) cm (exp -3) before the closest approach to 10(exp 5) cm (exp -3) just 30 s after the closest approach, an amount two orders of magnitude higher than the electron density. Assuming that the inconsistency between the electron and ion number densities is due to the presence of dust particles that are collecting the missing electron charges, we present dusty plasma characteristics down to sub-micron particle sizes. By assuming a differential dust number density for a range in dust sizes and by making use of Langmuir probe data, the dust densities for certain lower limits in dust size distribution were estimated. In order to achieve the dust densities of micrometer and larger sized grains comparable to the ones reported in the literature. we show that the power law size distribution must hold down to at least 0.03 micron such that the total differential number density is dominated by the smallest sub-micron sized grains. The total dust number density in Enceladus' plume is of the order of l0(exp 2) cm(exp -3) reducing to 1 cm(exp -3) in the E- ring. The dust density for micrometer and larger sized grains is estimated to be about 10(exp -4) cm(exp -3) in the plume while it is about 10(exp -6) - 10(exp -7) cm(exp -3) in the E-ring. Dust charge for micron sized grains is estimated to be about eight thousand electron charges reducing to below one hundred electron charges for 0.03 micron sized grains. The effective dusty plasma Debye length is estimated and compared with intergrain distance as well as the electron Debye length. The maximum dust charging time of 1.4 h is found for 0.03 11mmicron sized grains just 1 min before the closest approach. The charging time decreases substantially in the plume where it is only a fraction of a second for 1 micron sized grains, 1 s for 0.l micron sized grains and about 10 s for 0.03 micron sized grains.

Shafiq, Muhammad↗
NASA NTRS1998

The NASA Low-Pressure Turbine Flow Physics Program

An overview of the NASA Lewis Low-Pressure Turbine (LPT) Flow Physics Program will be presented. The program was established in response to the aero-engine industry's need for improved LPT efficiency and designs. Modern jet engines have four to seven LPT stages, significantly contributing to engine weight. In addition, there is a significant efficiency degradation between takeoff and cruise conditions, of up to 2 points. Reducing the weight and part count of the LPT and minimizing the efficiency degradation will translate into fuel savings. Accurate prediction methods of LPT flows and losses are needed to accomplish those improvements. The flow in LPT passages is at low Reynolds number, and is dominated by interplay of three basic mechanisms: transition, separation and wake interaction. The affecting parameters traditionally considered are Reynolds number, freestream turbulence intensity, wake frequency parameter, and the pressure distribution (loading). Three-dimensional effects and additional parameters, particularly turbulence characteristics like length scales, spectra and other statistics, as well as wake turbulence intensity and properties also play a role. The flow of most interest is on the suction surface, where large losses are generated as the flow tends to separate at the low Reynolds numbers. Ignoring wakes, a common flow scenario, there is laminar separation, followed by transition on the separation bubble and turbulent reattachment. If transition starts earlier the separation will be eliminated and the boundary layer will be attached leading to the well known bypass transition issues. In contrast, transition over a separation bubble is closer to free shear layer transition and was not investigated as well, particularly in the turbine environment. Unsteadiness created by wakes complicates the picture. Wakes induce earlier transition, and the calmed regions trailing the induced turbulent spots can delay or eliminate separation via shear stress modification. Three-dimensional flow physics and geometry will have strong effects. Altogether a very complex and challenging problem emerges. The objective of the program is to provide improved models and physical understanding of the complex flow, which are essential for accurate prediction of flow and losses in the LPT. Experimental, computational and analytical work as complementing and augmenting approaches are used. The program involves industry, universities and research institutes, and other government laboratories. It is characterized by strong interaction among participants, quick dissemination of results, and responsiveness to industry's needs. The presentation will describe the work elements. Highlighting some activities in progress are experiments on simulated blade suction surface in low-speed wind tunnels, on curved wall, and on a flat-plate, both with pressure gradient. In the area of computation, assessment of existing models is performed using RANS (Reynolds Averaged Navier Stokes) simulations. Laminar flow DNS was completed. Analytical studies of instability and receptivity in attached and separated flows were started. In the near future the program is moving to include wake effects and development of improved modeling. Experimental work in preparation stages are: (1) Addition of wakes to the curved tunnel experiment; (2) Low-speed rotating rig experiment on GE90 engine LPT; and (3) Transonic cascade. In the area of computation, it is expected to move from model assessment towards development of improved models. In addition, a new project of Large Eddy Simulation (LES) of LPT is to begin and will provide numerical data bases. It is planned to implement the emerging improved models in a multistage turbomachinery code and to validate against the GE90 engine LPT.

Ashpis, David E.↗
NASA NTRS1990

Schlieren studies of compressibility effects on dynamic stall of airfoils in transient pitching motion

Compressibility effects on the flowfield of an airfoil executing rapid transient pitching motion from 0 - 60 degrees over a wide range of Mach numbers and pitching rates were studied using a stroboscopic schlieren flow visualization technique. The studies have led to the first direct experimental documentation of multiple shocks on the airfoil upper surface flow for certain conditions. Also, at low Mach numbers, additional coherent vortical structures were found to be present along with the dynamic stall vortex, whereas at higher Mach numbers, the flow was dominated by a single vortex. The delineating Mach number for significant compressibility effects was 0.3 and the dynamic stall process was accelerated by increasing the Mach number above that value. Increasing the pitch rate monotonically delayed stall to angles of attack as large as 27 degrees.

Chandrasekhara, M. S.↗
NASA NTRS1990

Evolution of the total ozone field during the breakdown of the Antarctic circumpolar vortex

Nine years of total ozone measurements from the TOMS on Nimbus 7 are used to study the evolution of the Southern Hemisphere total ozone field during the breakdown of the Antarctic circumpolar vortex. The TOMS data provide detailed maps of the morphology of the ozone field and reliable estimates of the vertically integrated meridional transport of ozone during the springtime period when the breakdown occurs (September, October, November). In estimating the ozone transport, chemical effects, including those thought to be responsible for the Antarctic ozone hole, are neglected. This approximation appears to be valid for time scales of a few days to a week. On this time scale, local ozone changes are primarily due to transport. Planetary-scale waves, especially zonal wave numbers 1 and 2 dominate the eddy variance and ozone transport. Wave number 1 is quasistationary, while wave number 2 is eastward moving with a period of about 10 days.

Bowman, Kenneth P.↗
NASA NTRS1986

Numerical simulation of boundary layers. Part 2: Ribbon-induced transition in Blasius flow

The early three-dimensional stages of transition in Blasius boundary layers are studied by numerical solution of the Navier-Stokes equations. A finite-amplitude two-dimensional wave and random low-amplitude three-dimensional disturbances are introduced. Rapid amplification of the three-dimensional components is observed and leads to transition. For intermediate amplitudes of the two-dimensional wave the breakdown is of subharmonic type, and the dominant spanwise wave number increases with the amplitude. For high amplitudes the energy of the fundamental mode is comparable to the energy of the subharmonic mode, but never dominates it; the breakdown is of mixed type. Visualizations, energy histories, and spectra are presented. The sensitivity of the results to various physical and numerical parameters is studied. Agreement with experimental and theoretical results is discussed.

Spalart, P.↗
NASA NTRS2000

Global Lightning Variations Caused by Changes in Thunderstorm Flash Rate and by Changes in Number of Thunderstorms

Global lightning activity is highly variable on many time scales. This variability is attributable to changes in the flash rate per thunderstorm, the number of thunderstorms, or a combination. The TRMM Mission offers lightning observations from the Optical Transient Detector (OTD) and the Lightning Imaging Sensor (LIS) in space. Both are used to examine the response of these parameters to thermodynamic forcing of deep convection on the diurnal and annual time scales. On both time scales, the changes in the number of storms dominate the variations in total lightning activity. On the diurnal time scale, the mean flash rate appears to vary with cloud buoyancy, peaking in early afternoon and declining in late afternoon, but the contribution of number of thunderstorms is 2-3 times greater that the mean storm flash rate. On the annual time scale, almost all of the total lightning response is due to changes in the number of storms, with a negligible contribution from flash rate. Evidence is presented that the LIS/OTD 'area' is a meaningful objective identifier for a thunderstorm.

Williams, E. K.↗
NASA NTRS1993

Efficiency and reliability enhancements in propulsion flowfield modeling

The implementation of traditional CFD algorithms in practical propulsion related flowfields often leads to dramatic reductions in efficiency and/or robustness. The present research is directed at understanding the reasons for this deterioration and finding methods to circumvent it. Work to date has focussed on low Mach number regions, viscous dominated regions, and high grid aspect ratios. Time derivative preconditioning, improved definition of the local time stepping, and appropriate application of boundary conditions are employed to decrease the required time to obtain a solution, while maintaining accuracy. A number of cases having features typical of rocket engine flowfields are computed to demonstrate the improvement over conventional methods. These cases include laminar and turbulent high Reynolds number flat plate boundary layers, flow over a backward-facing step, a diffusion flame, and wall heat-flux calculations in a turbulent converging-diverging nozzle. Results from these cases show convergence that is virtually independent of the local Mach number and the grid aspect ratio, which translates to a convergence speed-up of up to several orders of magnitude over conventional algorithms. Current emphasis is in extending these results to three-dimensional flows with highly stretched grids.

Buelow, Philip E. O.↗
NASA NTRS1990

Relaxation processes in a turbulent compressible magnetofluid

The compressible extensions of time asymptotic relaxation states of incompressible two-dimensional magnetohydrodynamic turbulence are studied. A polytropic equation of state is used with viscous and resistive dissipation. The incompressible case is known to allow three distinct time asymptotic types of behavior: magnetic energy dominated relaxation, kinetic energy dominated relaxation, and cross helicity dominated relaxation. At low Mach numbers the incompressible scenario is reproducible from the compressible simulations, and compressibility plays only a secondary role. At moderate, but still subsonic, Mach numbers the distinct incompressible processes are still recognizable, but strong compressibility features dominate the high-wave-number regime of several simulations. In particular, the magnetic and kinetic energy dominated simulations display regions of strong acoustic turbulence near the dissipation scale.

Ghosh, S.↗
NASA NTRS1988

On the two-dimensional modelling of horizontal CVD reactors and its limitations

In order to clarify under what conditions two-dimensional theoretical models can realistically account for the actual growth rate distribution in horizontal CVD reactors, a systematic numerical study was performed, and the results were compared with published experimental data for the MOCVD of gallium arsenide. Three models of increasing complexity were studied. It is shown that two-dimensional models can produce realistic predictions only for reactors with large width-to-height (aspect) ratios, that are operated at subcritical Rayleigh number,i.e., under dominant forced flow conditions. But even then, depending on the carrier gas, thermal (Soret) diffusion must be included in the model. Furthermore, velocity corrections for finite aspect ratios must be made, and buoyancy effects can be significant in the entrance region. On the other hand, the results are relatively insensitive to axial diffusion, temperature dependence of transport properties (versus averaged properties), and an assumed Poiseuille velocity profile (versus a rigorous solution of the velocity field). Yet, boundary layer approximations are prone to give unrealistic results at the low Reynolds numbers typical for CVD operations.

Ouazzani, Jalil↗