SEARCH · Engineering Papers
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.
An experimental investigation of the effects of the geometry of body surface, Reynolds number, stream turbulence, and a roughness element (wire) on the reattachment of separated laminar boundary-layer flow on a bent flat plate is presented and discussed. The flow mechanisms determining reattachment of the boundary layer are analyzed and discussed.
Radar image simulation techniques are used to determine the character of SAR images of area-extensive targets, for the cases of flat underlying terrain and of moderate relief, with a view to the severity of elevation change effects on the detection and recognition of boundaries and shapes. The experiment, which demonstrated these effects on shapes and boundaries, was performed in order to establish Seasat-A SAR performance. Attention is given to the geometry/propagation effects in range perspective imaging that must be known for information extraction.
To study the azimuthal development of boundary-layer instabilities, a controlled, laser-generated perturbation was created in the freestream of the Boeing/U.S. Air Force Office of Scientific Research Mach 6 Quiet Tunnel. The freestream perturbation convected downstream in the wind tunnel to interact with a flared-cone model. The flared cone is a body of revolution bounded by a circular arc with a 3 m radius. Pressure transducers were used to measure a wave packet generated in the cone boundary layer by the freestream perturbation. Nine of these sensors formed three stations of azimuthal arrays and were used to determine the azimuthal variation of the wave packets in the boundary layer. The freestream laser-generated perturbation was positioned upstream of the model in three different configurations: along the centerline axis, offset from the centerline axis by 1.5 mm, and offset from the centerline axis by 3.0 mm. When the freestream perturbation was offset from the centerline of a flared cone with a 1.0 mm nose radius, a larger wave packet was generated on the side toward which the perturbation was offset. As a result, transition occurred earlier on that side. The offset perturbation did not have as large of an effect on the boundary layer of a nominally sharp flared cone.
This study describes the acceleration and displacement methods for use in the recovery of coupled system boundary forces. A simple two degree of freedom system has been used for illustration. The effect of the choice of method for use with indeterminate or over-constrained boundaries has been investigated. It has specifically looked at results from a simple two dimensional beam problem using both methods. Much work has been done on the effect of Craig-Bampton modal truncation system displacements and forces, however, little work has been done on system level modal truncation. The findings of this study indicate that the effect of this system level truncation is significant. This may be particularly true for the 35 Hz system cutoff frequency that is required by the space shuttle. From this study's findings, recommendations for areas of study with space shuttle payload systems are made.
Nonlinear oscillations and other motions of large axially symmetric liquid drops in zero gravity are studied numerically by a boundary-integral method. The effect of small viscosity is included in the computations by retaining first-order viscous terms in the normal stress boundary condition. This is accomplished by making use of a partial solution of the boundary-layer equations which describe the weak vortical surface layer. Small viscosity is found to have a relatively large effect on resonant mode coupling phenomena.
The determination of ignition or thermal explosion in an oxidizing porous body of material, as described by a dimensionless reaction-diffusion equation of the form .tu = .2u + .e-1/u over the bounded region O, is critically reexamined from a modern perspective using numerical methodologies. First, the classic stationary model is revisited to establish the proper reference frame for the steady-state solution space, and it is demonstrated how the resulting nonlinear two-point boundary value problem can be reexpressed as an initial value problem for a system of first-order differential equations, which may be readily solved using standard algorithms. Then, the numerical procedure is implemented and thoroughly validated against previous computational results based on sophisticated path-following techniques. Next, the transient nonstationary model is attacked, and the full nonlinear form of the reaction-diffusion equation, including a generalized convective boundary condition, is discretized and expressed as a system of linear algebraic equations. The numerical methodology is implemented as a computer algorithm, and validation computations are carried out as a prelude to a broad-ranging evaluation of the assembly problem and identification of the watershed critical initial temperature conditions for thermal ignition. This numerical methodology is then used as the basis for studying the relationship between the shape of the critical initial temperature distribution and the corresponding spatial moments of its energy content integral and an attempt to forge a fundamental conjecture governing this relation. Finally, the effects of dynamic boundary conditions on the classic storage problem are investigated and the groundwork is laid for the development of an approximate solution methodology based on adaptation of the standard stationary model.
A computational investigation of a pressure-gain combustor system for gas turbine applications is presented. The system consists of a valved pulse combustor and an ejector, housed within a shroud. The study focuses on two enhancements to previous models, related to the valve and ejector components. First, a new poppet inlet valve system is investigated, replacing the previously used reed valve configuration. Secondly, a new computational approach to approximating the effects of choked turbine inlet guide vanes present immediately downstream of the Ejector-Enhanced Resonant Pulse Combustor (EERPC) is investigated. Instead of specifying a back pressure at the EERPC exit boundary (as was done in previous studies) the new model adds a converging-diverging (CD) nozzle at the exit of the EERPC. The throat area of the CD nozzle can be adjusted to obtain the desired back pressure level and total mass flow rate. The results presented indicate that the new poppet valve configuration performs nearly as well as the original reed valve system, and that the addition of the CD nozzle is an effective method to approximate the exit boundary effects of a turbine present downstream of the EERPC. Furthermore, it is shown that the more acoustically reflective boundary imposed by a nozzle as compared to a constant pressure surface does not significantly affect operation or performance.
Utilizing a conceptual model for tropical convection and observational data for water vapor, the maintenance of the vertical distribution of the tropical tropospheric water vapor is discussed. While deep convection induces large-scale subsidence that constrains the turbulent downgradient mixing to within the convective boundary layer and effectively dries the troposphere through downward advection, it also pumps hydrometeors into the upper troposphere, whose subsequent evaporation appears to be the major source of moisture for the large-scale subsiding motion. The development of upper-level clouds and precipitation from these clouds may also act to dry the outflow, thus explaining the low relative humidity near the tropopause. A one-dimensional model is developed to simulate the mean vertical structure of water vapor in the tropical troposphere. It is also shown that the horizontal variation of water vapor in the tropical troposphere above the trade-wind boundary layer can be explained by the variation of a moisture source that is proportional to the amount of upper-level clouds. Implications for the nature of water vapor feedback in global warming are discussed.
The original aim of this work was to shed some light on the physics of turbulence over rough walls using large-eddy simulations and the suboptimal-control wall boundary conditions introduced by Nicoud et al. It was hoped that, if that algorithm was used to fit the mean velocity profile of the simulations to that of a rough-walled channel, instead of to a smooth one, the wall stresses introduced by the control algorithm would give some indication of what aspects of rough walls are most responsible for the modification of the flow in real turbulence. It was similarly expected that the structure of the resulting velocity fluctuations would share some of the characteristics of rough-walled flows, thus again suggesting what is intrinsic and what is accidental in the effect of geometric wall roughness. A secondary goal was to study the effect of 'unphysical' boundary conditions on the outside flow by observing how a relatively major change of the target velocity profile, and therefore presumably of the applied wall stresses, modifies properties such as the dominant length scales of the velocity fluctuations away from the wall. As will be seen below, this secondary goal grew more important during the course of the study, which was carried out during a short summer visit of the first two authors to the CTR. It became clear that there are open questions about the way in which the control algorithm models the boundary conditions, even for smooth walls, and that these questions make the physical interpretation of the results difficult. Considerable more work in that area seems to be needed before even relatively advanced large-eddy simulations, such as these, can be used to draw conclusions about the physics of wall-bounded turbulent flows. The numerical method is the same as in Nicoud et al. The modifications introduced in the original code are briefly described in section 2, but the original paper should be consulted for a full description of the algorithm. The results are presented in section 3 and summarized in section 4. The elementary properties of turbulence over rough walls which are used in the text have been taken from recent reviews such as Raupach et al. or Jimenez.
The second case for this workshop builds upon the isolated synthetic jet of Case 1 by adding a crossflow, with no streamwise pressure gradient, for the developing jet to interact with. Formally, Case 2 examines the interaction of a single, isolated, synthetic jet and a fully turbulent zero-pressure gradient boundary layer. The resulting flow has many of the characteristics that need to be modeled with fidelity if the results of the calculations are to serve as the basis for research and design with active flow control devices. These include the turbulence in the boundary layer, the time-evolution of the large vortical structure emanating from the jet orifice and its subsequent interaction with and distortion by the boundary layer turbulence, and the effect of the suction cycle on the boundary layer flow. In a synthetic jet, the flow through the orifice and out into the outer flowfield alternates between an exhaust and a suction cycle, driven by the contraction and expansion of a cavity internal to the actuator. In the present experiment, the volume changes in the internal cavity are accomplished by replacing one of the rigid walls of the cavity, the wall opposite the orifice exit, with a deformable wall. This flexible wall is driven by a bottom-mounted moveable piston. The piston is driven electro-mechanically. The synthetic jet issues into the external flow through a circular orifice. In the present experiment, this orifice has a diameter of 0.250 inches (6.35 mm). The flow is conceptually similar to that documented in Schaeffler [1]. To document the flow, several measurement techniques were utilized. The upstream boundary conditions (in-flow conditions), and several key phase-averaged velocity profiles were measured with a 3-component laser-Doppler velocimetry system. Phase-averaged velocity field measurements were made with both stereo digital particle image velocimetry and 2-D digital particle image velocimetry as the primary measurement system. Surface pressure measurements were made utilizing an electronically scanned pressure system.
The listings, user's instructions, sample inputs, and sample outputs of two computer programs which are especially useful in obtaining an approximate solution of the viscous flow over an arbitrary nonlifting three dimensional body are provided. The first program performs a potential flow solution by a well known panel method and readjusts this initial solution to account for the effects of the boundary layer displacement thickness, a nonuniform but unidirectional onset flow field, and the presence of air intakes and exhausts. The second program is effectually a geometry package which allows the user to change or refine the shape of a body to satisfy particular needs without a significant amount of human intervention. An effort to reduce the cruise drag of light aircraft through an analytical study of the contributions to the drag arising from the engine cowl shape and the foward fuselage area and also that resulting from the cooling air mass flowing through intake and exhaust sites on the nacelle is presented. The programs may be effectively used to determine the appropriate body modifications or flow port locations to reduce the cruise drag as well as to provide sufficient air flow for cooling the engine.
In this work, we demonstrate a nonlinear Hall effect due to the boundary spin accumulation in Pt films grown on Al 2 O 3 substrates. This Hall effect and the previously demonstrated Hanle magnetoresistance provide a complete picture of the spin-precession control of the spin and charge transport at the boundary of a spin-orbit coupled material, which we refer to as spin Hall Hanle effects (SHHE). We also show that the SHHE can be employed to measure the spin diffusion length, the spin Hall angle, and the spin relaxation time of heavy metal without the need of magnetic interface or the input from other measurements. The comprehensive demonstration of SHHE in such a simple system suggests they may be ubiquitous and needs to be considered for unraveling the spin and charge transport in more complex thin-film structures of spin-orbit coupled materials.
We report experimental observations on the effect of plasma boundary shaping towards balanced double-null (DN) configuration on the plasma performance in KSTAR. The transition from a single-null to a DN configuration resulted in improved plasma performance, manifested through changes in the pedestal region, decreased density, and core MHD activity variation. Specifically, the DN transition led to a wider and higher pedestal structure, accompanied by grassy edge-localized modes (ELMs) characteristics. The density decrease was a prerequisite for performance enhancement during DN shaping, increasing fast ion confinement. Optimizing the plasma near the core region was associated with the suppression of sawtooth instabilities and the occurrence of fishbone modes during the DN transition. Integrated modeling demonstrated that secondary effects of the DN shaping could increase core thermal energy confinement.
Regional surface fluxes are essential lower boundary conditions for large scale numerical weather and climate models and are the elements of global budgets of important trace gases. Surface properties affecting the exchange of heat, moisture, momentum and trace gases vary with length scales from one meter to hundreds of km. A classical difficulty is that fluxes have been measured directly only at points or along lines. The process of scaling up observations limited in space and/or time to represent larger areas was done by assigning properties to surface classes and combining estimated or calculated fluxes using an area weighted average. It is not clear that a simple area weighted average is sufficient to produce the large scale from the small scale, chiefly due to the effect of internal boundary layers, nor is it known how important the uncertainty is to large scale model outcomes. Simultaneous aircraft and tower data obtained in the relatively simple terrain of the western Alaska tundra were used to determine the extent to which surface type variation can be related to fluxes of heat, moisture, and other properties. Surface type was classified as lake or land with aircraft borne infrared thermometer, and flight level heat and moisture fluxes were related to surface type. The magnitude and variety of sampling errors inherent in eddy correlation flux estimation place limits on how well any flux can be known even in simple geometries.
Regional surface fluxes are essential lower boundary conditions for large scale numerical weather and climate models and are the elements of global budgets of important trace gases. Surface properties affecting the exchange of heat, moisture, momentum and trace gases vary with length scales from one meter to hundreds of km. A classical difficulty is that fluxes have been measured directly only at points or along lines. The process of scaling up observations limited in space and/or time to represent larger areas was done by assigning properties to surface classes and combining estimated or calculated fluxes using an area weighted average. It is not clear that a simple area weighted average is sufficient to produce the large scale from the small scale, chiefly due to the effect of internal boundary layers, nor is it known how important the uncertainty is to large scale model outcomes. Simultaneous aircraft and tower data obtained in the relatively simple terrain of the western Alaska tundra were used to determine the extent to which surface type variation can be related to fluxes of heat, moisture, and other properties. Surface type was classified as lake or land with aircraft borne infrared thermometer, and flight level heat and moisture fluxes were related to surface type. The magnitude and variety of sampling errors inherent in eddy correlation flux estimation place limits on how well any flux can be known even in simple geometries.
As the result of global, non-parallel flow stability analysis the single value of the disturbance growth-rate and respective frequency is obtained. This complex value characterizes the stability of the whole flow configuration and is not referred to any particular flow pattern. The global analysis assures that all the flow elements (wake, boundary and shear layer) are taken into account. The physical phenomena connected with the wake instability are properly reproduced by the global analysis. This enhances the investigations of instability of any 2-D flows, including ones in which the boundary layer instability effects are known to be of dominating importance. Assuming fully 2-D disturbance form, the global linear stability problem is formulated. The system of partial differential equations is solved for the eigenvalues and eigenvectors. The equations, written in the pure stream function formulation, are discretized via FDM using a curvilinear coordinate system. The complex eigenvalues and corresponding eigenvectors are evaluated by an iterative method. The investigations performed for various Reynolds numbers emphasize that the wake instability develops into the Karman vortex street. This phenomenon is shown to be connected with the first mode obtained from the non-parallel flow stability analysis. The higher modes are reflecting different physical phenomena as for example Tollmien-Schlichting waves, originating in the boundary layer and having the tendency to emerge as instabilities for the growing Reynolds number. The investigations are carried out for a circular cylinder, oblong ellipsis and airfoil. It is shown that the onset of the wake instability, the waves in the boundary layer, the shear layer instability are different solutions of the same eigenvalue problem, formulated using the non-parallel theory. The analysis offers large potential possibilities as the generalization of methods used till now for the stability analysis.
We consider a fourth order extension to MacCormack's scheme. The original extension was fourth order only for the inviscid terms but was second order for the viscous terms. We show how to modify the viscous terms so that the scheme is uniformly fourth order in the spatial derivatives. Applications are given to some boundary layer flows. In addition, for applications to shear flows the effect of the outflow boundary conditions are very important. We compare the accuracy of several of these different boundary conditions for both boundary layer and shear flows. Stretching at the outflow usually increases the oscillations in the numerical solution but the addition of a filtered sponge layer (with or without stretching) reduces such oscillations. The oscillations are generated by insufficient resolution of the shear layer. When the shear layer is sufficiently resolved then oscillations are not generated and there is less of a need for a nonreflecting boundary condition.
We consider a fourth order extension to MacCormack's scheme. The original extension was fourth order only for the inviscid terms but was second order for the viscous terms. We show how to modify the viscous terms so that the scheme is uniformly fourth order in the spatial derivatives. Applications are given to some boundary layer flows. In addition, for applications to shear flows the effect of the outflow boundary conditions are very important. We compare the accuracy of several of these different boundary conditions for both boundary layer and shear flows. Stretching at the outflow usually increases the oscillations in the numerical solution but the addition of a filtered sponge layer (with or without stretching) reduces such oscillations. The oscillations are generated by insufficient resolution of the shear layer. When the shear layer is sufficiently resolved then oscillations are not generated and there is less of a need for a nonreflecting boundary condition.