Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “viscous effects”

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 217 records · Page 12

Extension of a three-dimensional viscous wing flow analysis

Three-dimensional unsteady viscous effects can significantly influence the performance of fixed and rotary wing aircraft. These effects are important in both flows about helicopter rotors in forward flight and flows about 3-D (swept and tapered) supercritical wings. A computational procedure for calculating such flow field is developed, and therefore would be of great value in the design process as well as in understanding the corresponding flow phenomena. The procedure is based upon an alternating direction technique employing the Linearized Block Implicit method for solving 3-D viscous flow problems. In order to demonstrate the viability of this method, 2-D and 3-D problems are computed. These include the flow over a 2-D NACA 0012 airfoil under steady and oscillating conditions, and the steady, skewed, 3-D flow on a flat plate. Although actual 3-D flows over wings were not obtained, the ground work was laid for considering such flows. The description of the computational procedure and results are given.

Weinberg, Bernard C.↗

An integrated aerodynamic/propulsion study for generic aero-space planes based on waverider concepts

Research efforts related to the development of a unified aerospace plane analysis based on waverider technology are summarized. Viscous effects on the forebodies of cone-derived waverider configurations were studied. A simple means for determining the average skin friction coefficient of laminar boundary layers was established. This was incorporated into a computer program that provides lift and drag coefficients and lift/drag ratio for on-design waveriders when the temperature and Reynolds number based on length are specified. An effort was made to carry out parabolized Navier-Stokes (PNS) calculations for cone-derived waveriders. When the viscous terms were turned off (in the Euler mode) computations for elliptic cone-derived waveriders could be carried out for a wide range of on-design and off-design situations. Work related to waveriders derived from power law shocks is described in some detail.

Emanuel, G.↗

Biowaste resistojet propellant system biological and functional analysis, task 3

Exhaust flow contamination aspects of the biowaste resistojet are studied by evaluating effects of operating pressure, temperature and composition. Biowaste propellant mixtures considered are comprised of: (1) The Sabatier reactor effluent; (2) the effluent of the cabin carbon dioxide molecular sieve; and (3) water and water vapor from various sources. Results show that plume shapes of resistojet thrusters in the 25 to 100 mlb range exhibit greater apex angles for a given density contour than a scaled inviscid jet. Operation at low thrust, low pressure and high temperature accentuates this pluming due to viscous effects in the nozzle flow. Since the biowaste resistojet effluent is traveling at high velocity in the plume away from the aircraft it is found to be a superior method of damping than the ambient venting.

Source record↗

Parametric Geometry, Structured Grid Generation, and Initial Design Study for REST-Class Hypersonic Inlets

Computational design and analysis of three-dimensional hypersonic inlets with shape transition has been a significant challenge due to the complex geometry and grid required for three-dimensional viscous flow calculations. Currently, the design process utilizes an inviscid design tool to produce initial inlet shapes by streamline tracing through an axisymmetric compression field. However, the shape is defined by a large number of points rather than a continuous surface and lacks important features such as blunt leading edges. Therefore, a design system has been developed to parametrically construct true CAD geometry and link the topology of a structured grid to the geometry. The Adaptive Modeling Language (AML) constitutes the underlying framework that is used to build the geometry and grid topology. Parameterization of the CAD geometry allows the inlet shapes produced by the inviscid design tool to be generated, but also allows a great deal of flexibility to modify the shape to account for three-dimensional viscous effects. By linking the grid topology to the parametric geometry, the GridPro grid generation software can be used efficiently to produce a smooth hexahedral multiblock grid. To demonstrate the new capability, a matrix of inlets were designed by varying four geometry parameters in the inviscid design tool. The goals of the initial design study were to explore inviscid design tool geometry variations with a three-dimensional analysis approach, demonstrate a solution rate which would enable the use of high-fidelity viscous three-dimensional CFD in future design efforts, process the results for important performance parameters, and perform a sample optimization.

Ferlemann, Paul G.↗

Reynolds Number Effects on the Richtmyer-Meshkov Instability

This presentation compares the results of two very different experimental studies of Richtmyer-Meshkov instability: shock tube experiments in which an air/SF6 interface is accelerated by a weak shock wave; and incompressible experiments in which a box containing two different density miscible liquids is impulsively accelerated by bouncing it off of a fixed coil spring. Both experiments are initiated with sinusoidal initial perturbations. The interface perturbation initially remains sinusoidal as it grows in amplitude, but eventually the interfacial vorticity concentrates into points, forming a row of line vortices of alternating sign. The Reynolds number based on vortex circulation ranges from 1,000 to 45,000 in these experiments. It is found that viscous effects have a large, quantifiable effect on the evolution of the individual vortices. The effects of viscosity on the overall perturbation amplitude, however, are small and will be compared to theory.

Niederhaus, Charles↗

Transonic Navier-Stokes computations of strake-generated vortex interactions for a fighter-like configuration

Transonic Euler/Navier-Stokes computations are accomplished for wing-body flow fields using a computer program called Transonic Navier-Stokes (TNS). The wing-body grids are generated using a program called ZONER, which subdivides a coarse grid about a fighter-like aircraft configuration into smaller zones, which are tailored to local grid requirements. These zones can be either finely clustered for capture of viscous effects, or coarsely clustered for inviscid portions of the flow field. Different equation sets may be solved in the different zone types. This modular approach also affords the opportunity to modify a local region of the grid without recomputing the global grid. This capability speeds up the design optimization process when quick modifications to the geometry definition are desired. The solution algorithm embodied in TNS is implicit, and is capable of capturing pressure gradients associated with shocks. The algebraic turbulence model employed has proven adequate for viscous interactions with moderate separation. Results confirm that the TNS program can successfully be used to simulate transonic viscous flows about complicated 3-D geometries.

Reznick, Steve↗

The Viscoplastic Behavior of Natural Hydrate‐Bearing Sandy‐Silts Under Uniaxial Strain Compression ( K 0 Loading)

Abstract The in‐situ stress state and geomechanical properties of hydrate‐bearing sediments impact hydrate formation and gas production strategies. We explore the uniaxial strain compression and stress evolution of natural hydrate‐bearing sandy‐silts from Green Canyon Block 955 in the deep‐water Gulf of Mexico. We performed constant rate of strain uniaxial strain experiments, interrupted by periods where we held the axial stress constant, to explore the vertical deformation and the evolution of the ratio of lateral to axial effective stress ( K 0 ) with time. The hydrate‐bearing sandy‐silt is stiffer and has a larger K 0 than the equivalent hydrate‐free sediment upon loading. During stress holds, the void ratio decreases sigmoidally with the log of time, and K 0 converges to isotropic conditions. We interpret that during loading, the hydrate bears the load and deforms. With time, the hydrate redistributes the load and K 0 increases. We used a viscoelastic model to describe the behavior. The model accurately captures deformation and K 0 trends but does not reproduce all the complex interactions of the hydrate with the porous skeleton. We anticipate that viscous effects within hydrate sediments will impact reservoir compression and stresses during production (hours to days), result in isotropic stress state over geological timescales, and explain the creeping movement in submarine landslides.

58 GEOSCIENCES↗

A Technique for Simulating Conditions of Walking and Performing Other Self-Locomotive Activities on the Lunar Terrain

One of the most important and probably the most interesting phase of a manned lunar mission will be the time the astronauts spend outside their vehicle on the moon's surface taking scientific measurements, exploring the surface features, surveying possible sites for a lunar base, inspecting their vehicle and preparing it for their return trip. Because the lunar gravity is only one-sixth that of the earth gravity, the explorers undoubtedly will have to adjust their accustomed methods of walking, climbing, jumping and performing other self-locomotive activities in order to carry out these various tasks. In as much as the over-all success of the lunar mission will depend to a large extent upon the self-reliance of the explorers, it will be necessary to have extensive knowledge of the effects of the moon's reduced gravity on the physical capabilities of man and of man's ability to adopt to the new environment prior to the planning and execution of the mission. At the present time there is a dearth of information on this subject due primarily to the lack of a practical technique for simulating the reduced gravity. Several techniques such as immersion in water and riding in an airplane flying a Keplerian trajectory have been used for zero-g or weightlessness studies to determine the physical capabilities of man but these techniques are limited in their usefulness either by restrictions imposed by the viscous effect of the water or by the short duration and small test area available in an airplane. Consequently, an effort was made at the NASA Langley Research Center to devise a new technique that would provide a realistic simlation of a reduced gravity for unlimited periods of time and allow freedom of movement over considerable distances. This paper concerns itself with a discussion of the newly developed simulation technique and a presentation of some preliminary results which were obtained utilizing a working model based on this scheme.

Hewes, Donald E.↗

Investigations related to the inviscid-viscous interaction in transonic flows about finite 3-D wings

Inviscid small-disturbance theory has been shown to predict three-dimensional transonic flows about finite wings reasonably well as long as viscous effects are negligible. In order to include these effects, the inviscid small-disturbance solution of Bailey and Ballhaus (1975) has been combined with a finite-difference solution for Prandtl's boundary-layer equations. This solution employs the conditionally stable Krause (1968) scheme, implicit in the direction normal to the wall, to cope with the domain-of-dependence problem that arises for reverse cross flow. To be consistent with the inviscid-flow solution, the boundary layer is computed in the representative wing planform plane which is transformed into rectangular shape in the computational domain. The flow has been assumed turbulent, and a scalar eddy-viscosity model is adopted. The interaction between inviscid and viscous flow is modeled with the help of the displacement surface which is added to the geometric wing shape. Sample distributions of displacement thickness for swept wings are presented for weak and strong interaction cases.

Kordulla, W.↗

Model of a Hollow Cathode Insert Plasma

A 2-D axisymmetric fluid model of the plasma in the insert region of a hollow cathode is presented. The level of sophistication included in the model is motivated in part by the need to determine quantitatively plasma fluxes to the emitter surface. The ultimate goal is to assess whether plasma effects can degrade the life of impregnated inserts beyond those documented throughout the 30-50 year history of vacuum cathode technologies. Results from simulations of a 1.2-cm diameter cathode operating at a discharge current of 25 A, and a gas flow rate of 5 sccm, suggest that approximately 10 A of electron current, and 3.5 A of ion current return to the emitter surface. The total emitted electron current computed by the model is about 35 A. Comparisons with plasma measurements suggest that anomalous heating of the plasma due to two-stream instabilities is possible near the orifice region. Solution to the heavy species energy equation, with classical transport and no viscous effects, predicts heavy species temperatures as high as 2640 K.

hollow cathodes↗

Effects of Atwood and Reynolds numbers on the evolution of buoyancy-driven homogeneous variable-density turbulence

The evolution of buoyancy-driven homogeneous variable-density turbulence (HVDT) at Atwood numbers up to 0.75 and large Reynolds numbers is studied by using high-resolution direct numerical simulations. To help understand the highly non-equilibrium nature of buoyancy-driven HVDT, the flow evolution is divided into four different regimes based on the behaviour of turbulent kinetic energy derivatives. The results show that each regime has a unique type of dependence on both Atwood and Reynolds numbers. It is found that the local statistics of the flow based on the flow composition are more sensitive to Atwood and Reynolds numbers compared to those based on the entire flow. It is also observed that, at higher Atwood numbers, different flow features reach their asymptotic Reynolds-number behaviour at different times. The energy spectrum defined based on the Favre fluctuations momentum has less large-scale contamination from viscous effects for variable-density flows with constant properties, compared to other forms used previously. The evolution of the energy spectrum highlights distinct dynamical features of the four flow regimes. Thus, the slope of the energy spectrum at intermediate to large scales evolves from -7/3 to -1 , as a function of the production-to-dissipation ratio. The classical Kolmogorov spectrum emerges at intermediate to high scales at the highest Reynolds numbers examined, after the turbulence starts to decay. Finally, the similarities and differences between buoyancy-driven HVDT and the more conventional stationary turbulence are discussed and new strategies and tools for analysis are proposed.

42 ENGINEERING↗

Earth Return Aerocapture for the TransHab/Ellipsled Vehicle

The current architecture being considered by NASA for a human Mars mission involves the use of an aerocapture procedure at Mars arrival and possibly upon Earth return. This technique would be used to decelerate the vehicles and insert them into their desired target orbits, thereby eliminating the need for propulsive orbital insertions. The crew may make the interplanetary journey in a large, inflatable habitat known as the TransHab. It has been proposed that upon Earth return, this habitat be captured into orbit for use on subsequent missions. In this case, the TransHab would be complimented with an aeroshell, which would protect it from heating during the atmospheric entry and provide the vehicle with aerodynamic lift. The aeroshell has been dubbed the "Ellipsled" because of its characteristic shape. This paper reports the results of a preliminary study of the aerocapture of the TransHab/Ellipsled vehicle upon Earth return. Undershoot and overshoot boundaries have been determined for a range of entry velocities, and the effects of variations in the atmospheric density profile, the vehicle deceleration limit, the maximum vehicle roll rate, the target orbit, and the vehicle ballistic coefficient have been examined. A simple, 180 degree roll maneuver was implemented in the undershoot trajectories to target the desired 407 km circular Earth orbit. A three-roll sequence was developed to target not only a specific orbital energy, but also a particular inclination, thereby decreasing propulsive inclination changes and post-aerocapture delta-V requirements. Results show that the TransHab/Ellipsled vehicle has a nominal corridor width of at least 0.7 degrees for entry speeds up to 14.0 km/s. Most trajectories were simulated using continuum flow aerodynamics, but the impact of high-altitude viscous effects was evaluated and found to be minimal. In addition, entry corridor comparisons have been made between the TransHab/Ellipsled and a modified Apollo capsule which is also being considered as the crew return vehicle; because of its slightly higher lift-to-drag ratio, the TransHab has a modest advantage with regard to corridor width. Stagnation-point heating rates and integrated heat loads were determined for a range of vehicle ballistic coefficients and entry velocities.

Muth, W. D.↗

An experimental investigation of a supersonic vortical flow

Although much research has been done on subsonic vortical flow, the current understanding of these flows remains limited. The effect is characterized of adding swirl to a supersonic jet. The motive is to study the enhancement of supersonic mixing in order to provide more efficient fuel injectors for supersonic combustion (scramjet) engines. The vortical flow was created by tangential injection into a swirl chamber ahead of a converging and/or diverging nozzle. The amount of swirl was varied by changing the number of tangential injection holes and with the removal of the end piece, the jet could be run without swirl. Shadowgraphy, conventional schlieren, and focusing schlieren were used to obtain a qualitative understanding of the jet flow structure. It was determined that an increase in swirl produced an increase in the shear layer growth. Pressure and temperature probes were used to obtain more flow data. The probe data compared favorably with the theoretical calculations, except in the viscous core where viscous effects were not considered negligible. These results verified that a supersonic vortical flow was being created with a maximum helix angle of 33 degs.

Levey, Brian S.↗

Wall interference and boundary simulation in a transonic wind tunnel with a discretely slotted test section

A computational simulation of a transonic wind tunnel test section with longitudinally slotted walls is developed and described herein. The nonlinear slot model includes dynamic pressure effects and a plenum pressure constraint, and each slot is treated individually. The solution is performed using a finite-difference method that solves an extended transonic small disturbance equation. The walls serve as the outer boundary conditions in the relaxation technique, and an interaction procedure is used at the slotted walls. Measured boundary pressures are not required to establish the wall conditions but are currently used to assess the accuracy of the simulation. This method can also calculate a free-air solution as well as solutions that employ the classical homogeneous wall conditions. The simulation is used to examine two commercial transport aircraft models at a supercritical Mach number for zero-lift and cruise conditions. Good agreement between measured and calculated wall pressures is obtained for the model geometries and flow conditions examined herein. Some localized disagreement is noted, which is attributed to improper simulation of viscous effects in the slots.

Al-Saadi, Jassim A.↗

Design of a Mach-3 Nozzle for TBCC Testing in the NASA LaRC 8-ft High Temperature Tunnel

A new nozzle is being constructed for the NASA Langley Research Center 8-Foot High Temperature Tunnel. The axisymmetric nozzle was designed with a Mach-3 exit flow for testing Turbine-Based Combined-Cycle engines at a Mach number in the vicinity of the transition from turbojet to ramjet operation. The nozzle contour was designed using the NASA Langley IMOCND computer program which solves the potential equation using the classical method of characteristics. To include viscous effects, the design procedure iterated the MOC contour generation with CFD Navier-Stokes calculations, adjusting MOC input parameters until target nozzle-exit conditions were achieved in the Navier-Stokes calculations. The design process was complicated by a requirement to use the final 29.5 inches of an existing 54.5-inch exit-diameter Mach-5 nozzle contour. This was accomplished by generating a Mach-3 contour that matched the radius of the Mach-5 contour at the match point and using a 3rd order polynomial to create a smooth transition between the two contours. During the final evaluation of the design it was realized that the throat diameter is more than half that of the upstream mixing chamber. This led to the concern that large vortical structures generated in the mixer would persist downstream, affecting nozzle-exit flow. This concern was addressed by analyzing the results of three-dimensional, viscous, numerical simulations of the entire flowfield, from the exit of the facility combustor to the nozzle exit. An analysis of the solution indicated that large scale structures do not pass through the throat and that both the total temperature and species (CO2) are well mixed in the mixer, providing uniform flow to the nozzle and subsequently the test cabin.

Gaffney, Richard L., Jr.↗

Computation of viscous transonic flow over porous airfoils

The viscous effects on transonic flow past an airfoil which contains a shallow cavity beneath a porous surface are studied numerically. The porous region occupies a small portion of the total airfoil surface, and is located near the shock. Both an interactive boundary layer (IBL) algorithm and a thin-layer Navier-Stokes (TLNS) algorithm have been modified for use in studying the outer flow, whereas a stream-function formulation has been used to model the inner flow in the small cavity. The coupling procedure at the porous surface is based on Darcy's law and on the assumption of a constant total presusre in the cavity. In addition, a modified Baldwin-Lomax turbulence model is used to consider the transpired turbulent boundary layer in the TLNS approach, and the Cebeci-Smith turbulence model is used in the IBL approach. According to the present analysis, a porous surface can reduce the wave drag appreciably, but it can also increase viscous losses. As has been observed experimentally, the numerical results indicate that the total drag is reduced at higher Mach numbers and increased at lower Mach numbers when the angles of attack are small. Furthermore, the streamline patterns of passive-shock and boundary-layer interaction are revealed in this study.

Chow, Chuen-Yen↗

Unsteady Flow Field in a Multistage Axial Flow Compressor

The flow field in a multistage compressor is three-dimensional, unsteady, and turbulent with substantial viscous effects. Some of the specific phenomena that has eluded designers include the effects of rotor-stator and rotor-rotor interactions and the physics of mixing of velocity, pressure, temperature and velocity fields. An attempt was made, to resolve experimentally, the unsteady pressure and temperature fields downstream of the second stator of a multistage axial flow compressor which will provide information on rotor-stator interaction effects and the nature of the unsteadiness in an embedded stator of a three stage axial flow compressor. Detailed area traverse measurements using pneumatic five hole probe, thermocouple probe, semi-conductor total pressure probe (Kulite) and an aspirating probe downstream of the second stator were conducted at the peak efficiency operating condition. The unsteady data was then reduced through an ensemble averaging technique which splits the signal into deterministic and unresolved components. Auto and cross correlation techniques were used to correlate the deterministic total temperature and velocity components (acquired using a slanted hot-film probe at the same measurement locations) and the gradients, distributions and relative weights of each of the terms of the average passage equation were then determined. Based on these measurements it was observed that the stator wakes, hub leakage flow region, casing endwall suction surface corner region, and the casing endwall region away from the blade surfaces were the regions of highest losses in total pressure, lowest efficiency and highest levels of unresolved unsteadiness. The deterministic unsteadiness was found to be high in the hub and casing endwall regions as well as on the pressure side of the stator wake. The spectral distribution of hot-wire and kulite voltages shows that at least eight harmonics of all three rotor blade passing frequencies are present at this measurement location. In addition to the basic three rotor blade passing frequencies (R1, R2 and R3) and their harmonics, various difference frequencies such as (2R1 -R2) and (2R3-R2) and their harmonics are also observed. These difference frequencies are due to viscous and potential interactions between rotors 1, 2 and 3 which are sensed by both the total pressure and aspirating probes at this location. Significant changes occur to the stator exit flow features with passage of the rotor upstream of the stator. Because of higher convection speeds of the rotor wake on the suction surface of the downstream stator than on the pressure side, the chopped rotor wake was found to be arriving at different times on either side of the stator wake. As the rotor passes across the stator.

Suryavamshi, N.↗

High pressure/high temperature multiphase simulations of dodecane injection to nitrogen: Application on ECN Spray-A

The present work investigates the complex phenomena associated with pressure/high temperature dodecane injection for the Engine Combustion Network (ECN) Spray-A case, employing more elaborate thermodynamic closures, to avoid well known deficiencies concerning density and speed of sound prediction using traditional cubic models. A tabulated thermodynamic approach is proposed here, based on log 10 (p)-T tables, providing very high accuracy across a large range of pressures, spanning from 0 to 2500 bar, with only a small number of interpolation points. The tabulation approach is directly extensible to any thermodynamic model, existing or to be developed in the future. Here NIST REFPROP properties are used, combined with PC-SAFT Vapor-Liquid-Equilibrium to identify the liquid in mixtures penetration, hence avoiding the use of an arbitrary threshold for mass fraction. Identified liquid and vapour penetration are compared against experimental data from the ECN database showing a good agreement, within approximately 3–8% for axial penetration of liquid, 2% for vapor axial penetration and within experimental uncertainty for radial distribution of mass fraction. Analysis of the vortex evolution indicates that driving mechanisms behind the jet break-up are vortex tilting/stretching, then baroclinic torque, leading to Rayleigh-Taylor instabilities, closely followed by vortex dilation and finally viscous effects.

42 ENGINEERING↗