Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “incompressible flow”

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 469 records · Page 26

Three-dimensional application of the Johnson-King turbulence model for a boundary-layer direct method

The Johnson-King turbulence model as extended to three-dimensional flows was evaluated using finite-difference boundary-layer direct method. Calculations were compared against the experimental data of the well-known Berg-Elsenaar incompressible flow over an infinite swept-wing. The Johnson-King model, which includes the nonequilibrium effects in a developing turbulent boundary-layer, was found to significantly improve the predictive quality of a direct boundary-layer method. The improvement was especially visible in the computations with increased three-dimensionality of the mean flow, larger integral parameters, and decreasing eddy-viscosity and shear stress magnitudes in the streamwise direction; all in better agreement with the experiment than simple mixing-length methods.

Kavsaoglu, Mehmet S.↗

A viscous/potential flow interaction analysis method for multi-element infinite swept wings, volume 1

An analysis method and computer program have been developed for the calculation of the viscosity dependent aerodynamic characteristics of multi-element infinite swept wings in incompressible flow. The wing configuration consisting at the most of a slat, a main element and double slotted flap is represented in the method by a large number of panels. The inviscid pressure distribution about a given configuration in the normal chord direction is determined using a two dimensional potential flow program employing a vortex lattice technique. The boundary layer development over each individual element of the high lift configuration is determined using either integral or finite difference boundary layer techniques. A source distribution is then determined as a function of the calculated boundary layer displacement thickness and pressure distributions. This source distribution is included in the second calculation of the potential flow about the configuration. Once the solution has converged (usually after 2-5 iterations between the potential flow and boundary layer calculations) lift, drag, and pitching moments can be determined as functions of Reynolds number.

Dvorak, F. A.↗

Two-dimensional compressible flow in turbomachines with conic flow surfaces

A general method of analysis is developed for two-dimensional, steady, compressible flow in stators or rotors of radial and mixed flow turbomachines with conic flow surfaces (surfaces of right circular cones generated by center line of flow passage in the axial-radial plane). The variables taken into account are: (1) tip speed of the rotor, (2) flow rate, (3) blade shape, (4) variation in passage height with radius, (5) number of blades, and (6) cone angle of the flow surface. Relaxation methods are used to solve the nonlinear differential equation for the stream function. Two numerical examples are presented; one for compressible and the other for incompressible flow in a centrifugal compressor with thin, straight blades. The results of these examples are given by plots of the streamlines, constant velocity-ratio lines, and constant pressure-ratio line.

Stanitz, John D↗

Field Sensitivity Analysis of Turbulence Model Parameters for Flow Over a Wing

Reynolds-averaged-Navier-Stokes (RANS) turbulence models are a critical tool in computational-fluid-dynamics simulations of aerodynamic systems, but simulation results can be highly sensitive to RANS-model parameter choices. Sensitivity analysis can be used to quantify these impacts, and the objective of this study is to demonstrate field sensitivity analysis with respect to ten parameters in the 2003 Menter shear-stress-transport (SST) turbulence model. Here, the analysis is demonstrated for an application relevant to wind energy, namely, flow over a NACA 0015 wing at 12 degree angle of attack and a Reynolds number of 1.5 x 10 6 . We quantify sensitivity using Sobol indices and the mean-squared gradient, which are estimated using polynomial chaos and active subspace models, respectively. Our results indicate that there are substantial spatial variations in parameter sensitivities, with different sets of most-sensitive parameters near the wing, as well as in the downstream wake, consistent with the physical interpretations of the turbulence model inputs. We show that, for this particular turbulence model and flow, simultaneous dimension reduction is possible across all quantities of interest, enabling efficient exploration of model outcomes. Ultimately, this analysis provides new insights into turbulence model parameter sensitivities in incompressible flows, and also demonstrates the implementation of field sensitivity analysis for applications relevant to aerodynamics simulations.

17 WIND ENERGY↗

Application of an invariant second-order closure model to compressible turbulent shear layers

A second-order closure model for two-dimensional, compressible shear flows is investigated using the invariant modeling technique developed by Donaldson. The invariant model parameters were originally selected by comparison of model predictions with critical experimental data on basic incompressible flows. Additional modeling for compressible flows has been introduced. A number of different shear flows such as the free shear layer, flat plate boundary layer and a simulation of the mixing region of a chemical laser are computed. The model is consistent with first-order closure turbulence models for equilibrium flows and is further capable of predicting nonequilibrium flows that cannot be correctly solved by eddy viscosity models. The results are in generally good agreement with experimental measurements but suggest the need for inclusion of a Mach number dependent model of the pressure diffusion terms in order to adequately represent high speed compressible flows.

Varma, A. K.↗

Progress in high-lift aerodynamic calculations

The current work presents progress in the effort to numerically simulate the flow over high-lift aerodynamic components, namely, multi-element airfoils and wings in either a take-off or a landing configuration. The computational approach utilizes an incompressible flow solver and an overlaid chimera grid approach. A detailed grid resolution study is presented for flow over a three-element airfoil. Two turbulence models, a one-equation Baldwin-Barth model and a two equation k-omega model are compared. Excellent agreement with experiment is obtained for the lift coefficient at all angles of attack, including the prediction of maximum lift when using the two-equation model. Results for two other flap riggings are shown. Three-dimensional results are presented for a wing with a square wing-tip as a validation case. Grid generation and topology is discussed for computing the flow over a T-39 Sabreliner wing with flap deployed and the initial calculations for this geometry are presented.

Rogers, Stuart E.↗

Instability Wave-Streak Interactions in a Supersonic Boundary Layer

Optimal initial conditions for transient growth in a two-dimensional boundary layer flow correspond to stationary, counter-rotating vortices that subsequently develop into streamwise elongated streaks, which are characterized by an alternating pattern of low and high streamwise velocity. For incompressible flows, previous studies have shown that boundary layer modulation due to streaks below a threshold amplitude level can stabilize the Tollmien-Schlichting instability waves, resulting in a delay in the onset of laminar-turbulent transition. In the supersonic regime, the linearly, most-amplified waves become three-dimensional, corresponding to oblique, first-mode waves. This change in the character of dominant instabilities leads to an important change in the transition process, which is now dominated by oblique breakdown via nonlinear interactions between pairs of first-mode waves that propagate at equal but opposite angles with respect to the free stream. Because the oblique breakdown process is characterized by a rapid amplification of stationary streamwise streaks, artificial excitation of such streaks may be expected to promote transition in a supersonic boundary layer. Indeed, suppression of those streaks has been shown to delay the onset of transition in prior literature. Consistent with those findings, the present study shows that optimally growing stationary streaks indeed destabilize the first-mode waves, but only when the spanwise wavelength of the instability waves is equal to or smaller than twice the streak spacing. Transition in a benign disturbance environment typically involves first-mode waves with significantly longer spanwise wavelengths, and hence, these waves are stabilized by the optimal growth streaks. Thus, as long as the amplification factors for the destabilized, short wavelength instability waves remain below the threshold level for transition, a significant net stabilization is achieved, yielding a transition delay that is comparable to the length of the laminar region in the uncontrolled case.

Boundary layer transition↗

Advanced stability analysis for laminar flow control

Five classes of problems are addressed: (1) the extension of the SALLY stability analysis code to the full eighth order compressible stability equations for three dimensional boundary layer; (2) a comparison of methods for prediction of transition using SALLY for incompressible flows; (3) a study of instability and transition in rotating disk flows in which the effects of Coriolis forces and streamline curvature are included; (4) a new linear three dimensional instability mechanism that predicts Reynolds numbers for transition to turbulence in planar shear flows in good agreement with experiment; and (5) a study of the stability of finite amplitude disturbances in axisymmetric pipe flow showing the stability of this flow to all nonlinear axisymmetric disturbances.

Orszag, S. A.↗

Reynolds number effects on the turbulence field in compressible boundary layers

Detailed experiments were conducted in a zero pressure gradient, supersonic turbulent boundary layer, including measurements of the three components of velocity fluctuations and the turbulent shear stress, for Reynolds numbers ranging from 11.7 million to 105 million at a freestream Mach number of 2.3. The mean flow measurements established the fully developed and equilibrium nature of the boundary layer. Measurements of the turbulence field show that the vertical and transverse fluctuations are essentially equal throughout the boundary layer at all Reynolds numbers, a feature that is different from observations in incompressible flows. The data show that the boundary layer exhibits similarity in the turbulence profiles for the entire Reynolds number range and agrees with previous compressible and incompressible data using Morkovin's scaling to account for compressibility effects.

Acharya, M.↗

A quasi-three-dimensional blade surface boundary layer analysis for rotating blade rows

A quasi-three-dimensional approximation has been developed for a blade boundary layer which involves the calculation of the effect of nonzero pressure gradients, turbulent flow, and blade twist, but includes only a simple coupling between streamlines. The resulting set of equations is solved using Keller's box scheme. The solution scheme is checked against available incompressible flow solutions and then applied to a NASA low aspect ratio transonic compressor stage for which extensive experimental and computational data are available. It is found that the three-dimensional boundary layer separates significantly sooner and has a much larger influence on rotor performance than would be expected from a two-dimensional analysis.

Thompkins, W. T., Jr.↗

Transonic airfoil analysis and design in nonuniform flow

A nonuniform transonic airfoil code is developed for applications in analysis, inverse design and direct optimization involving an airfoil immersed in propfan slipstream. Problems concerning the numerical stability, convergence, divergence and solution oscillations are discussed. The code is validated by comparing with some known results in incompressible flow. A parametric investigation indicates that the airfoil lift-drag ratio can be increased by decreasing the thickness ratio. A better performance can be achieved if the airfoil is located below the slipstream center. Airfoil characteristics designed by the inverse method and a direct optimization are compared. The airfoil designed with the method of direct optimization exhibits better characteristics and achieves a gain of 22 percent in lift-drag ratio with a reduction of 4 percent in thickness.

Chang, J. F.↗

Numerical study of juncture flows

The present paper describes a computational study of laminar/turbulent and subsonic/supersonic horseshoe vortex systems generated by a cylindrical protuberance mounted on a flat plate. Various vortex structures are predicted and discussed. Low-speed laminar juncture flows are computed to determine the Reynolds number effect with the same incoming boundary-layer thickness. For a low subsonic laminar flow, the number of vortex arrays increases with the Reynolds number, in agreement with both experimental and numerical observations. Qualitative comparisons are made along with the computations, experimental observations, and analytical work. For incompressible flow, the relationships among pressure extrema, vorticity, and singular points in flow structure are discussed. A parametric study of the effect of the free-stream Mach number on the flow structure for laminar flow is conducted. The juncture flow when the incoming flow is turbulent and supersonic is computed.

Chen, Chung-Lung↗

Amplitude-dependent neutral modes in compressible boundary layer flows

The ideas of Benney and Bergeron (1969) and Davies (1970) on nonlinear critical layers are extended, and some new nonlinear neutral modes are computed for compressible boundary layer flow. A special case of the work is when the generalized inflexion point criterion holds. Neutral modes are found for a range of phase-speeds, dependent on the Mach number, and the properties of these are discussed. As in the linear case when the flow is relatively supersonic, multiple neutral modes exist. The behavior of the neutral amplitude in some limiting cases is also considered, and it is found that the results are significantly different from that in incompressible flow when the flow is locally supersonic.

Gajjar, J. S. B.↗

Amplitude-dependent neutral modes in compressible boundary layer flows

The ideas of Benney and Bergeron (1969) and Davies (1969) on nonlinear critical layers are extended, and some new nonlinear neutral modes are computed for compressible boundary layer flow. A special case of the work is when the generalized inflexion point criterion holds. Neutral modes are found for a range of phase-speeds, dependent on the Mach number, and the properties of these are discussed. As in the linear case when the flow is relatively supersonic, multiple neutral modes exist. The behavior of the neutral amplitude in some limiting cases is also considered, and it is found that the results are significantly different from that in incompressible flow when the flow is locally supersonic.

Gajjar, J. S. B.↗

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.↗

Sound Generation by Flow over a Two-Dimensional Cavity

Sound generated by flow over a cavity at a Mach number of 0.1 and a Reynolds number based on cavity length of 5000 is calculated. The computation utilizes a two part technique where the time-dependent incompressible flow is first obtained and then a second calculation is performed for the compressible aspects of the flow. This second calculation utilizes a grid and numerical scheme designed for resolution of acoustic waves. The cavity flowfield is observed to oscillate quite regularly at the Strouhal number of 0.58 which produces an acoustic source of the same frequency. Time histories, spectra, and directivity of the sound radiation are computed.

Hardin, Jay C.↗

Two-dimensional compressible flow in centrifugal compressors with straight blades

Six numerical examples are presented for steady, two-dimensional, compressible, nonviscous flow in centrifugal compressors with thin straight blades, the center lines of which generate the surface of a right circular cone when rotated about the axis of the compressor. A seventh example is presented for incompressible flow. The solutions were obtained in a region of the compressors, including the impeller tip, that was considered to be unaffected by the diffuser vanes or by the impeller-inlet configuration. Each solution applies to radial and mixed flow compressors with various cone angles but with the same angle between blades on the conic flow surface. The solution also apply to radial and mixed flow turbines with the rotation and the flow direction reversed. The effects of variations in the following parameters were investigated: (1) flow rate, (2) impeller-tip speed, (3) variation of passage height with radius, and (4) angle between blades on conic flow surface. The numerical results are presented in plots of the streamlines and constant Mach number lines. Correlation equations are developed whereby the flow conditions in any impeller with straight blades can be determined (in the region investigated by this analysis) for all operating conditions.

Stanitz, John D↗

Finite element analyses of plane thermal entry-length flows

It is pointed out that the design of actively cooled engine and airframe structures for hypersonic flight requires strong interaction between fluid, thermal, and structural analyses. Research programs are being conducted to develop an integrated thermal-structural analysis capability based upon the finite element method. In connection with these programs, a finite element engineering approach has been developed based upon a number of assumptions customarily used in practical heat transfer analysis. An evaluation is conducted of the finite element engineering approach for entry-length flows. The evaluation makes use of an analysis of plane thermal entry-length flows by a finite element approach utilizing mean fluid temperatures with a convection coefficient. Another analysis is based on a finite element continuum approach, taking into account the rigorous momentum, mass, and energy equations of viscous, incompressible flow. The engineering approach was found to give excellent agreement with the rigorous continuum approach for Peclet number values not less than 100.

Thornton, E. A.↗