Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “pressure gradient”

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 1,297 records · Page 72
All sourcesNASA NTRSDOE OSTI
NASA NTRS1995

Boundary layer transition studies

A small-scale wind tunnel previously used for turbulent boundary layer experiments was modified for two sets of boundary layer transition studies. The first study concerns a laminar separation/turbulent reattachment. The pressure gradient and unit Reynolds number are the same as the fully turbulent flow of Spalart and Watmuff. Without the trip wire, a laminar layer asymptotes to a Falkner & Skan similarity solution in the FPG. Application of the APG causes the layer to separate and a highly turbulent and approximately 2D mean flow reattachment occurs downstream. In an effort to gain some physical insight into the flow processes a small impulsive disturbance was introduced at the C(sub p) minimum. The facility is totally automated and phase-averaged data are measured on a point-by-point basis using unprecedently large grids. The evolution of the disturbance has been tracked all the way into the reattachment region and beyond into the fully turbulent boundary layer. At first, the amplitude decays exponentially with streamwise distance in the APG region, where the layer remains attached, i.e. the layer is viscously stable. After separation, the rate of decay slows, and a point of minimum amplitude is reached where the contours of the wave packet exhibit dispersive characteristics. From this point, exponential growth of the amplitude of the disturbance is observed in the detached shear layer, i.e. the dominant instability mechanism is inviscid. A group of large-scale 3D vortex loops emerges in the vicinity of the reattachment. Remarkably, the second loop retains its identify far downstream in the turbulent boundary layer. The results provide a level of detail usually associated with CFD. Substantial modifications were made to the facility for the second study concerning disturbances generated by Suction Holes for laminar flow Control (LFC). The test section incorporates suction through interchangeable porous test surfaces. Detailed studies have been made using isolated holes in the impervious test plate that used to establish the Blasius base flow. The suction is perturbed harmonically and data are averaged on the basis of the phase of the disturbance, for conditions corresponding to strong suction and without suction. The technique was enhanced by using up to nine multiple probes to reduce the experimental run-time. In both cases, 3D contour surfaces in the vicinity of the hole show highly 3D TS waves which fan out in the spanwise direction forming bow-shaped waves downstream. The case without suction has proved useful for evaluating calculation methods. With suction, the perturbations on the centerline are much stronger and decay less rapidly, while the TS waves in the far field are similar to the case without suction. Downstream, the contour surfaces of the TS waves develop spanwise irregularities which eventually form into clumps. The spanwise clumping is evidence of a secondary instability that could be associated with suction vortices. Designers of porous surfaces use Goldsmith's Criterion to minimize cross-stream interactions. It is shown that partial TS wave cancellation is possible, depending on the hole spacing, disturbance frequency and free-stream velocity. New high-performance Constant Temperature Hot-Wire Anemometers were designed and built, based on a linear system theory analysis that can be extended to arbitrary order. The motivation was to achieve the highest possible frequency reponse while ensuring overall system stability. The performance is equal to or superior to commercially available instruments at about 10% of the cost. Details, such as fabrication drawings and a parts list, have been published to enable the instrument to be construced by others.

Watmuff, Jonathan H.↗
NASA NTRS1994

A mechanism to produce a dawn-dusk component of plasma flow during magnetic reconnection in the magnetotail

Magnetic reconnection between antiparallel field lines in the magnetotail is generally thought to produce plasma acceleration in the earthward-tailward direction. However, measurements of the plasma velocity in the magnetotail during substorm activity sometimes reveal a dawn-dusk component of plasma flow. In this paper, we show that a dawn-dusk component of plasma acceleration may be produced during reconnection if the neutral line is not perpendicular to the magnetic field. In this case, Magnetohydrodynamic (MHD) simulations show that reconnection between antiparallel field lines will initially produce plasma acceleration that is nearly parallel to the neutral line because the magnetic tension force is not opposed by a pressure gradient force in this direction. As the magnetic field topology evolves to a steady state, the plasma flow direction rotates until it is nearly parallel to the plane that initially contained the antiparallel magnetic field lines before reconnection (hereafter referred to as the initial field plane). However, the time required to reach a steady state (typically several hundred seconds in the magnetotail region) may be greater than the time during which the reconnection process is active. Consequently, bursts of plasma flow with a dawn-dusk component may occur in the magnetotail. The initial acceleration along the neutral line depends on the angle theta (sub B) between the neutral line and the initial field plane, with the largest burst of plasma flow along the neutral line occuring when theta (s ub B) = 45 degs.

Hawkins, J. G.↗
NASA NTRS1995

Magnetohydrodynamic simulations of outflows from accretion disks

Magnetohydrodynamic simulations have been made of the formation of outflows from a Keplerian disk threaded by a magnetic field. The disk is treated as a boundary condition, where matter is ejected with Keplerian azimuthal speed and poloidal speed less than the slow magnetosonic velocity, and where boundary conditions on the magnetic field correspond to a highly conducting disk. Initially, the space above the disk, the corona, is filled with high specific entropy plasma in thermal equilibrium in the gravitational potential of the central object. The initial magnetic field is poloidal and is represented by a superposition of monopoles located below the plane of the disk. The rotation of the disk twists the initial poloidal magnetic field, and this twist propagates into the corona pushing and collimating matter into jetlike outflow in a cylindrical region. Matter outflowing from the disk flows and accelerates in the z-direction owing to both the magnetic and pressure gradient forces. The flow accelerates through the slow magnetosonic and Alfven surfaces and at larger distances through the fast magnetosonic surface. The flow velocity of the jet is approximately parallel to the z-axis, and the collimation results from the pinching force of the toroidal magnetic field. For a nonrotating disk no collimation is observed.

Ustyugova, G. V.↗
NASA NTRS1994

A hybrid simulation of contact discontinuity

Contact discontinuities in a collisionless plasma are studied by hybrid simulations, in which ions are treated as particles and electrons are considered as a fluid. It is demonstrated that contact discontinuity with a stable density ramp can exist in cases with a finite electron temperature. An electron pressure gradient is present across the contact discontinuity, leading to the presence of a parallel electric field and hence field-aligned potential increase (Delta Phi (sub parallel)) in the transition region. By reflecting ions at the discontinuity, this parallel electric potential peak reduces the interpenetration between hot and cold ions and maintains a stable density ramp across the contact discontinuity. The ratio of the field-aligned electric potential energy to ion thermal energy, e(Delta) Phi(sub parallel)/kT(sub i), is found to be an increasing function of T(sub e)/T(sub i), where T(sub e) and T(sub i) are respectively the electron and ion temperature.

Wu, B. H.↗
NASA NTRS1995

Observations in the vicinity of substorm onset: Implications for the substrom process

Multi-instrument data sets from the ground and satellites at both low and high altitude have provided new results concerning substorm onset and its source region in the magnetosphere. Twenty-six out of 37 substorm onset events showed evidence of azimuthally spaced auroral forms (AAFs) prior to the explosive poleward motion associated with optical substorm onset. AAFs can span 8 hours of local time prior to onset and generally propagate eastward in the morning sector. Onset itself is, however, more localized spanning only about 1 hour local time. AAF onset occur during time periods when the solar wind pressure is relatively high. AAFs brighten in conjunction with substorm onset leading to the conclusion that they are a growth phase activity casually related to substorm onset. Precursor activity associated with these AAFs is also seen near geosynchronous orbit altitude and examples show the relationship between the various instrumental definitions of substorm onset. The implied mode number (30 to 135) derived from this work is inconsistent with cavity mode resonances but is consistent with a modified flute/ballooning instability which requires azimuthal pressure gradients. The extended source region and the distance to the open-closed field line region constrain reconnection theory and local mechanisms for substorm onset. It is demonstrated that multiple onset substorms can exist for which localized dipolarizations and the Pi 2 occur simultaneously with tail stretching existing elsewhere. These pseudobreakups can be initiated by auroral streamers which originate at the most poleward set of arc systems and drift to the more equatorward main UV oval. Observations are presented of these AAFs in conjunction with low- and high-altitutde particle and magnetic field data. These place the activations at the interface between dipolar and taillike field lines probably near the peak in the cross-tail current. These onsets are put in the context of a new scenario for substorm morphology which employs individual modules which operate independently or couple together. This allows particular substorm events to be more accurately described and investigated.

Elphinstone, R. D.↗
NASA NTRS1995

Wing-section optimization for supersonic viscous flow

To improve the shape of a supersonic wing, an automated method that also includes higher fidelity to the flow physics is desirable. With this impetus, an aerodynamic optimization methodology incorporating thin-layer Navier-Stokes equations and sensitivity analysis had been previously developed. Prior to embarking upon the wind design task, the present investigation concentrated on testing the feasibility of the methodology, and the identification of adequate problem formulations, by defining two-dimensional, cost-effective test cases. Starting with two distinctly different initial airfoils, two independent shape optimizations resulted in shapes with similar features: slightly cambered, parabolic profiles with sharp leading- and trailing-edges. Secondly, the normal section to the subsonic portion of the leading edge, which had a high normal angle-of-attack, was considered. The optimization resulted in a shape with twist and camber which eliminated the adverse pressure gradient, hence, exploiting the leading-edge thrust. The wing section shapes obtained in all the test cases had the features predicted by previous studies. Therefore, it was concluded that the flowfield analyses and sensitivity coefficients were computed and fed to the present gradient-based optimizer correctly. Also, as a result of the present two-dimensional study, suggestions were made for the problem formulations which should contribute to an effective wing shape optimization.

Item, Cem C.↗
NASA NTRS1995

Wing-section optimization for supersonic viscous flows

The recent interest in High Speed Civil Transport (HSCT) has resulted in renewed research studies of optimized supersonic cruise transport configurations. Incorporation of flow viscosity effects in the design process of such a supersonic wing is currently under investigation. This may lead to more accurate problem formulations and, in turn, greater aerodynamic efficiency than can be obtained by the traditional, inviscid, linear theories. In this context, for a design code to be a candidate for a complex optimization problem, such as three-dimensional viscous supersonic wing design, it should be validated using simpler building-block shapes. To optimize the shape of a supersonic wing, an automated method that also includes higher fidelity to the flow physics is desirable. With this impetus, an aerodynamic optimization methodology incorporating Navier-Stokes equations and sensitivity analysis had been previously developed. Prior to embarking upon the wing design task, the present investigation concentrated on testing the flexibility of the methodology, and the identification of adequate problem formulations, by defining two-dimensional, cost-effective test cases. Starting with two distinctly different initial airfoils, two independent shape optimizations resulted in shapes with very similar features. Secondly, the normal section to the subsonic portion of the leading edge, which had a high normal angle-of-attack, was considered. The optimization resulted in a shape with twist and camber, which eliminated the adverse pressure gradient, hence, exploiting the leading-edge thrust. The wing section shapes obtained in all the test cases had the features predicted by previous studies. Therefore, it was concluded that the flowfield analyses and the sensitivity coefficients were computed and fed to the present gradient-based optimizer correctly. Also, as a result of the present two-dimensional study, suggestions were made for problem formulations which should contribute to an effective wing shape optimization.

Item, Cem C.↗
NASA NTRS1996

Implementation of a two-equation k-omega turbulence model in NPARC

The implementation of a two-equation k-omega turbulence model into the NPARC flow solver is described. Motivation for the selection of this model is given, major code modifications are outlined, new imputs to the code are described, and results are presented for several validation cases: an incompressible flow over a smooth flat plate, a subsonic diffuser flow, and a shock-induced separated flow. Comparison of results with the k-epsilon model indicate that the k-omega model predicts simple flows equally well whereas, for adverse pressure gradient flows, the k-omega model outperforms the other turbulence models in NPARC.

Yoder, Dennis A.↗
NASA NTRS1995

Development of a two-equation turbulence model for hypersonic flows: Evaluation of a low Reynolds number correction to the Kappa - epsilon two equation compressible turbulence model - Volume 1

The objective of the current research is the development of an improved k-epsilon two-equation compressible turbulence model for turbulent boundary layer flows experiencing strong viscous-inviscid interactions. The development of an improved model is important in the design of hypersonic vehicles such as the National Aerospace Plane (NASP) and the High Speed Civil Transport (HSCT). Improvements have been made to the low Reynolds number functions in the eddy viscosity and dissipation of solenoidal dissipation of the k-epsilon turbulence mode. These corrections offer easily applicable modifications that may be utilized for more complex geometries. The low Reynolds number corrections are functions of the turbulent Reynolds number and are therefore independent of the coordinate system. The proposed model offers advantages over some current models which are based upon the physical distance from the wall, that modify the constants of the standard model, or that make more corrections than are necessary to the governing equations. The code has been developed to solve the Favre averaged, boundary layer equations for mass, momentum, energy, turbulence kinetic energy, and dissipation of solenoidal dissipation using Keller's box scheme and the Newton spatial marching method. The code has been validated by removing the turbulent terms and comparing the solution with the Blasius solution, and by comparing the turbulent solution with an existing k-epsilon model code using wall function boundary conditions. Excellent agreement is seen between the computed solution and the Blasius solution, and between the two codes. The model has been tested for both subsonic and supersonic flat-plate turbulent boundary layer flow by comparing the computed skin friction with the Van Driest II theory and the experimental data of Weighardt; by comparing the transformed velocity profile with the data of Weighardt, and the Law of the Wall and the Law of the Wake; and by comparing the computed results of an adverse pressure gradient with the experimental data of Fernando and Smits. Good agreement is obtained with the experimental correlations for all flow conditions.

Knight, Doyle D.↗
NASA NTRS1996

Control of Unstable Waves in Three Dimensional Boundary Layers

Stability experiments are conducted in the Arizona State University Unsteady Wind Tunnel on a 45 deg. swept airfoil. The pressure gradient is designed to provide purely crossflow-dominated transition; that is, the boundary layer is subcritical to Tollmien- Schlichting (T-S) disturbances. The airfoil surface is hand polished to a 0.25 micron rms finish. Under these conditions, stationary crossflow disturbances grow to nonuniform amplitude due to submicron surface irregularities near the leading edge. Spectral de-compositions isolate single-mode growth rates for the fundamental and harmonic disturbances. The measurements show early nonlinear growth causing amplitude saturation well before transition. Comparisons with nonlinear PSE calculations show excellent agreement in both the amplitude saturation and the disturbance mode shape.

Saric, William S.↗
NASA NTRS1995

Nonlinear Landau damping and Alfven wave dissipation

Nonlinear Landau damping has been often suggested to be the cause of the dissipation of Alfven waves in the solar wind as well as the mechanism for ion heating and selective preacceleration in solar flares. We discuss the viability of these processes in light of our theoretical and numerical results. We present one-dimensional hybrid plasma simulations of the nonlinear Landau damping of parallel Alfven waves. In this scenario, two Alfven waves nonresonantly combine to create second-order magnetic field pressure gradients, which then drive density fluctuations, which in turn drive a second-order longitudinal electric field. Under certain conditions, this electric field strongly interacts with the ambient ions via the Landau resonance which leads to a rapid dissipation of the Alfven wave energy. While there is a net flux of energy from the waves to the ions, one of the Alfven waves will grow if both have the same polarization. We compare damping and growth rates from plasma simulations with those predicted by Lee and Volk (1973), and also discuss the evolution of the ambient ion distribution. We then consider this nonlinear interaction in the presence of a spectrum of Alfven waves, and discuss the spectrum's influence on the growth or damping of a single wave. We also discuss the implications for wave dissipation and ion heating in the solar wind.

Vinas, Adolfo F.↗
NASA NTRS1995

The tilts of corotating interaction regions at mid heliographic latitudes

Ulysses observations of corotating interaction regions (CIRs) at mid heliographic latitudes have shown that the flow downstream of the forward shock (or wave) on the leading edge of a CIR generally turns northward and westward, while the flow downstream of the reverse wave on the trailing edge generally turns southward and eastward. These systematic flow deflections are a natural consequence of large scale pressure gradients associated with the CIRs, and indicate that the forward waves tend to propagate toward and across the equator with increasing heliocentric distance, while the reverse waves tend to propagate toward the pole. Recent determinations of CIR shock normals using the Ulysses magnetic field data appear to confirm these plasma results (Burton, private communication). Numerical simulations indicate that these effects (which imply that CIRs are systematically tilted in the north-south direction at mid latitudes) are a natural consequence of the tilt of the solar magnetic dipole axis relative to the solar rotation axis. The present work utilizes a variety of techniques to analyze the flow deflections observed within CIRs from which we can infer the overall orientations of the CIRs and the speeds and directions of propagation of the waves. Where possible, the observations are quantitatively compared with the results of 3-dimensional MHD simulations.

Riley, Pete↗
NASA NTRS1995

Large-eddy simulations with wall models

The near-wall viscous and buffer regions of wall-bounded flows generally require a large expenditure of computational resources to be resolved adequately, even in large-eddy simulation (LES). Often as much as 50% of the grid points in a computational domain are devoted to these regions. The dense grids that this implies also generally require small time steps for numerical stability and/or accuracy. It is commonly assumed that the inner wall layers are near equilibrium, so that the standard logarithmic law can be applied as the boundary condition for the wall stress well away from the wall, for example, in the logarithmic region, obviating the need to expend large amounts of grid points and computational time in this region. This approach is commonly employed in LES of planetary boundary layers, and it has also been used for some simple engineering flows. In order to calculate accurately a wall-bounded flow with coarse wall resolution, one requires the wall stress as a boundary condition. The goal of this work is to determine the extent to which equilibrium and boundary layer assumptions are valid in the near-wall regions, to develop models for the inner layer based on such assumptions, and to test these modeling ideas in some relatively simple flows with different pressure gradients, such as channel flow and flow over a backward-facing step. Ultimately, models that perform adequately in these situations will be applied to more complex flow configurations, such as an airfoil.

Cabot, W.↗
NASA NTRS1995

Large-eddy simulation of flow around an airfoil on a structured mesh

The diversity of flow characteristics encountered in a flow over an airfoil near maximum lift taxes the presently available statistical turbulence models. This work describes our first attempt to apply the technique of large-eddy simulation to a flow of aeronautical interest. The challenge for this simulation comes from the high Reynolds number of the flow as well as the variety of flow regimes encountered, including a thin laminar boundary layer at the nose, transition, boundary layer growth under adverse pressure gradient, incipient separation near the trailing edge, and merging of two shear layers at the trailing edge. The flow configuration chosen is a NACA 4412 airfoil near maximum lift. The corresponding angle of attack was determined independently by Wadcock (1987) and Hastings & Williams (1984, 1987) to be close to 12 deg. The simulation matches the chord Reynolds number U(sub infinity)c/v = 1.64 x 10(exp 6) of Wadcock's experiment.

Kaltenbach, Hans-Jakob↗
NASA NTRS1995

Analysis of the two-point velocity correlations in turbulent boundary layer flows

The general objective of the present work is to explore the use of Rapid Distortion Theory (RDT) in analysis of the two-point statistics of the log-layer. RDT is applicable only to unsteady flows where the non-linear turbulence-turbulence interaction can be neglected in comparison to linear turbulence-mean interactions. Here we propose to use RDT to examine the structure of the large energy-containing scales and their interaction with the mean flow in the log-region. The contents of the work are twofold: First, two-point analysis methods will be used to derive the law-of-the-wall for the special case of zero mean pressure gradient. The basic assumptions needed are one-dimensionality in the mean flow and homogeneity of the fluctuations. It will be shown that a formal solution of the two-point correlation equation can be obtained as a power series in the von Karman constant, known to be on the order of 0.4. In the second part, a detailed analysis of the two-point correlation function in the log-layer will be given. The fundamental set of equations and a functional relation for the two-point correlation function will be derived. An asymptotic expansion procedure will be used in the log-layer to match Kolmogorov's universal range and the one-point correlations to the inviscid outer region valid for large correlation distances.

Oberlack, M.↗
NASA NTRS1996

Comparative Study Of Four Models Of Turbulence

Report presents comparative study of four popular eddy-viscosity models of turbulence. Computations reported for three different adverse pressure-gradient flowfields. Detailed comparison of numerical results and experimental data given. Following models tested: Baldwin-Lomax, Johnson-King, Baldwin-Barth, and Wilcox.

Menter, Florian R.↗
NASA NTRS1996

Turbulence Modeling and Computation of Turbine Aerodynamics and Heat Transfer

The objective of the present research is to develop improved turbulence models for the computation of complex flows through turbomachinery passages, including the effects of streamline curvature, heat transfer and secondary flows. Advanced turbulence models are crucial for accurate prediction of rocket engine flows, due to existance of very large extra strain rates, such as strong streamline curvature. Numerical simulation of the turbulent flows in strongly curved ducts, including two 180-deg ducts, one 90-deg duct and a strongly concave curved turbulent boundary layer have been carried out with Reynolds stress models (RSM) and algebraic Reynolds stress models (ARSM). An improved near-wall pressure-strain correlation has been developed for capturing the anisotropy of turbulence in the concave region. A comparative study of two modes of transition in gas turbine, the by-pass transition and the separation-induced transition, has been carried out with several representative low-Reynolds number (LRN) k-epsilon models. Effects of blade surface pressure gradient, freestream turbulence and Reynolds number on the blade boundary layer development, and particularly the inception of transition are examined in detail. The present study indicates that the turbine blade transition, in the presence of high freestream turbulence, is predicted well with LRN k-epsilon models employed. The three-dimensional Navier-Stokes procedure developed by the present authors has been used to compute the three-dimensional viscous flow through the turbine nozzle passage of a single stage turbine. A low Reynolds number k-epsilon model and a zonal k-epsilon/ARSM (algebraic Reynolds stress model) are utilized for turbulence closure. An assessment of the performance of the turbulence models has been carried out. The two models are found to provide similar predictions for the mean flow parameters, although slight improvement in the prediction of some secondary flow quantities has been obtained by the ARSM model. It's found that the wake profiles inside the endwall boundary layers are predicted better than those near the mid-span.

Lakshminarayana, B.↗
NASA NTRS1993

Effect of sound on boundary layer stability

Experiments are conducted in the Arizona State University Unsteady Wind Tunnel with a zero-pressure-gradient flat-plate model that has a 67:1 elliptical leading edge. Boundary-layer measurements are made of the streamwise fluctuating-velocity component in order to identify the amplified T-S waves that are forced by downstream-travelling, sound waves. Measurements are taken with circular 3-D roughness elements placed at the Branch 1 neutral stability point for the frequency under consideration, and then with the roughness element downstream of Branch 1. These roughness elements have a principal chord dimension equal to 2(lambda)(sub TS)/pi, of the T-S waves under study and are 'stacked' in order to resemble a Gaussian height distribution. Measurements taken just downstream of the roughness (with leading-edge T-S waves, surface roughness T-S waves, instrumentation sting vibrations and the Stokes wave subtracted) show the generation of 3-D-T-S waves, but not in the characteristic heart-shaped disturbance field predicted by 3-D asymptotic theory. Maximum disturbance amplitudes are found on the roughness centerline. However, some near-field characteristics predicted by numerical modelling are observed.

Saric, William S.↗