SEARCH · Engineering Papers

Results for “Boundary layer 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 649 records · Page 36

All sources NASA NTRS DOE OSTI

NASA NTRS1992

Improved two-equation k-omega turbulence models for aerodynamic flows

Two new versions of the k-omega two-equation turbulence model will be presented. The new Baseline (BSL) model is designed to give results similar to those of the original k-omega model of Wilcox, but without its strong dependency on arbitrary freestream values. The BSL model is identical to the Wilcox model in the inner 50 percent of the boundary-layer but changes gradually to the high Reynolds number Jones-Launder k-epsilon model (in a k-omega formulation) towards the boundary-layer edge. The new model is also virtually identical to the Jones-Lauder model for free shear layers. The second version of the model is called Shear-Stress Transport (SST) model. It is based on the BSL model, but has the additional ability to account for the transport of the principal shear stress in adverse pressure gradient boundary-layers. The model is based on Bradshaw's assumption that the principal shear stress is proportional to the turbulent kinetic energy, which is introduced into the definition of the eddy-viscosity. Both models are tested for a large number of different flowfields. The results of the BSL model are similar to those of the original k-omega model, but without the undesirable freestream dependency. The predictions of the SST model are also independent of the freestream values and show excellent agreement with experimental data for adverse pressure gradient boundary-layer flows.

Menter, Florian R.↗

NASA NTRS1996

A Structured Grid Based Solution-Adaptive Technique for Complex Separated Flows

The objective of this work was to enhance the predictive capability of widely used computational fluid dynamic (CFD) codes through the use of solution adaptive gridding. Most problems of engineering interest involve multi-block grids and widely disparate length scales. Hence, it is desirable that the adaptive grid feature detection algorithm be developed to recognize flow structures of different type as well as differing intensity, and adequately address scaling and normalization across blocks. In order to study the accuracy and efficiency improvements due to the grid adaptation, it is necessary to quantify grid size and distribution requirements as well as computational times of non-adapted solutions. Flow fields about launch vehicles of practical interest often involve supersonic freestream conditions at angle of attack exhibiting large scale separate vortical flow, vortex-vortex and vortex-surface interactions, separated shear layers and multiple shocks of different intensity. In this work, a weight function and an associated mesh redistribution procedure is presented which detects and resolves these features without user intervention. Particular emphasis has been placed upon accurate resolution of expansion regions and boundary layers. Flow past a wedge at Mach=2.0 is used to illustrate the enhanced detection capabilities of this newly developed weight function.

Thornburg, Hugh↗

NASA NTRS1989

A block-corrected subdomain solution procedure for recirculating flow calculations

This paper describes a robust and efficient subdomain solution procedure for two-dimensional recirculating flows. The solution domain is divided into a number of overlapping subdomains, and a direct fully coupled solution is obtained for each subdomain using a sparse matrix form of LU decomposition. An effective parabolic block correction procedure, which calculates global corrections to the tentative solution by a marching technique similar to that used for boundary layer flows, is used to accelerate the convergence of the basic procedure. The use of effective block correction is found to be essential for the success of the subdomain approach on strongly recirculating flows. In a number of laminar two-dimensional flows, the new block-corrected method performed extremely well, rivaling the best direct methods in execution time, while requiring substantially less computer storage. The new method proved to be from two to ten times faster than conventional iterative methods, while requiring only a moderate increase in storage.

Braaten, M. E.↗

DOE OSTI2022

On the Response of a Herschel–Bulkley Fluid Due to a Moving Plate

In this paper, we study the boundary-layer flow of a Herschel–Bulkley fluid due to a moving plate; this problem has been experimentally investigated by others, where the fluid was assumed to be Carbopol, which has similar properties to cement. The computational fluid dynamics finite volume method from the open-source toolbox/library OpenFOAM is used on structured quad grids to solve the mass and the linear momentum conservation equations using the solver “overInterDyMFoam” customized with non-Newtonian viscosity libraries. The governing equations are solved numerically by using regularization methods in the context of the overset meshing technique. The results indicate that there is a good comparison between the experimental data and the simulations. The boundary layer thicknesses are predicted within the uncertainties of the measurements. The simulations indicate strong sensitivities to the rheological properties of the fluid.

36 MATERIALS SCIENCE↗

NASA NTRS1983

Comparison of predicted and experimental external heat transfer around a film cooled cylinder in crossflow

Calculations were made of the film cooling provided by rows of holes around the circumference of a cylinder in crossflow and the results were compared to experimental data. The calculations and experimental data were for conditions that simulate most of those that are typical of air cooled turbine vane leading edges. Injection was from a single and multiple rows of holes located at different angular locations from the stagnation line. The holes in the rows were angled normal to the flow direction and at a 25 degree angle to the cylinder wall. The calculations and experimental data were for several constant values of blowing ratios for all rows and for different blowing ratios for each row, representing a simulation of a common coolant plenum supply to multiple rows of holes. The calculations were made using a finite difference boundary layer code, STAN5. Contrary to initial expectations that injection would trip the boundary layer flow into the turbulent regime, the results indicated that the high free stream acceleration apparently kept the flow laminar for holes in the first 45 degrees past stagnation. The trend in Stanton number reduction due to coolant injection was predicted with generally good agreement at the lower blowing rates, but for multile rows of holes, agreement was poor beyond the first row.

Stepka, F. S.↗

NASA NTRS1938

Flow Observations with Tufts and Lampblack of the Stalling of Four Typical Airfoil Sections in the NACA Variable-density Tunnel

A preliminary investigation of the stalling processes of four typical airfoil sections was made over the critical range of the Reynolds Number. Motion pictures were taken of the movements of small silk tufts on the airfoil surface as the angle of attack increased through a range of angles including the stall. The boundary-layer flow also at certain angles of attack was indicated by the patterns formed by a suspension of lampblack in oil brushed onto the airfoil surface. These observations were analyzed together with corresponding force-test measurements to derive a picture of the stalling processes of airfoils.

Abbott, Ira H↗

NASA NTRS1997

Development of a Low-Reynolds Number, Nonlinear kappa-epsilon Model for the Reduced Navier-Stokes Equations

Previous work at NASA LeRC has shown that flow distortions in aircraft engine inlet ducts can be significantly reduced by mounting vortex generators, or small wing sections, on the inside surface of the engine inlet. The placement of the vortex generators is an important factor in obtaining the optimal effect over a wide operating envelope. In this regard, the only alternative to a long and expensive test program which would search out this optimal configuration is a good prediction procedure which could narrow the field of search. Such a procedure has been developed in collaboration with NASA LeRC, and results obtained by NASA personnel indicate that it shows considerable promise for predicting the viscous turbulent flow in engine inlet ducts in the presence of vortex generators. The prediction tool is a computer code which numerically solves the reduced Navier-Stokes equations and so is commonly referred to as RNS3D. Obvious deficiencies in RNS3D have been addressed in previous work. Primarily, it is known that the predictions of the mean velocity field of a turbulent boundary layer flow approaching separation are not in good agreement with data. It was suggested that the use of an algebraic mixing-length turbulence model in RNS3D is at least partly to blame for this. Additionally, the current turbulence model includes an assumption of isotropy which will ultimately fail to capture turbulence-driven secondary flow known to exist in noncircular ducts.

Boger, David A.↗

NASA NTRS2018

Nonlinear Evolution and Breakdown of Azimuthally Compact Crossflow Vortex Pattern over a Yawed Cone

Hypersonic boundary-layer flows over a circular cone at moderate angle of incidence can support strong crossflow instability in between the windward and leeward rays on the plane of symmetry. Due to the more efficient excitation of stationary crossflow vortices by surface roughness, a possible path to transition in such flows corresponds to rapid amplification of the high-frequency secondary instabilities of finite amplitude stationary crossflow vortices. In the present paper, the previous analyses of crossflow instability over a 7- degree half-angle, yawed circular cone in a Mach 6 free stream have been extended to the nonlinear evolution of azimuthally localized crossflow vortex packets and the amplification characteristics and nonlinear breakdown of high-frequency secondary instabilities associated with those packets. A comparison between plane marching PSE and direct Navier-Stokes simulations (DNS) reveals favorable agreement in regard to mode shapes, most amplified disturbance frequencies, and N-factor evolution. In contrast, the quasi-parallel predictions are found to result in severe underprediction of the N-factors. The direct numerical simulations also indicate that the breakdown of secondary instabilities in a 3D hypersonic boundary layer shares certain common features with the previous computations of crossflow transition over subsonic swept wings.

Vortex structure↗

NASA NTRS1977

Turbulence processes and simple closure schemes

The problem of closure in turbulence in the case of two-point correlations resides in the existence of two unknowns E and W, the energy spectrum function and the transfer function, respectively, in the spectrum equation. In the case of weak turbulence, W is negligible. In case of higher correlations, closure can be effective by neglecting the inertia term in the highest order term used. Specifying a certain number of spectra at an initial time is also a way of getting around the closure problem. A simple case of turbulent shear flow is then considered, where two-point correlation equations are used and the velocity is broken into mean and fluctuating components. This yields a differential equation for the energy spectrum, the three terms of which are the energy spectrum, production term and dissipation term. They are plotted for a particular time. Similar analyses and comparisons with experiment are made for pipe and boundary layer flows.

Deissler, R. G.↗

NASA NTRS1982

Solution of 3-dimensional time-dependent viscous flows. Part 3: Application to turbulent and unsteady flows

A numerical scheme is developed for solving the time dependent, three dimensional compressible viscous flow equations to be used as an aid in the design of helicopter rotors. In order to further investigate the numerical procedure, the computer code developed to solve an approximate form of the three dimensional unsteady Navier-Stokes equations employing a linearized block implicit technique in conjunction with a QR operator scheme is tested. Results of calculations are presented for several two dimensional boundary layer flows including steady turbulent and unsteady laminar cases. A comparison of fourth order and second order solutions indicate that increased accuracy can be obtained without any significant increases in cost (run time). The results of the computations also indicate that the computer code can be applied to more complex flows such as those encountered on rotating airfoils. The geometry of a symmetric NACA four digit airfoil is considered and the appropriate geometrical properties are computed.

Weinberg, B. C.↗

NASA NTRS1983

The effect of a leading-edge slat on the dynamic stall of an oscillating airfoil

The dynamic stall characteristics of a slatted airfoil were investigated experimentally on a 2-ft-chord airfoil oscillating in pitch at M = 0.2 for a range of reduced frequency and mean angle of oscillation. The slat produced a flow that remained attached to the airfoil for angles well above those normally attained by the retreating blade of a helicopter during high speed flight. The dynamic stall vortex usually associated with these flight conditions was completely eliminated for all angles under 30 deg. Instantaneous surface pressure, lift, and pitching moment data are presented as a function of incidence throughout the oscillation cycle; a detailed analysis of instantaneous boundary-layer flow behavior for the various test conditions is included.

Carr, L. W.↗

NASA NTRS1959

Effects of Cross-Sectional Shape, Solidity, and Distribution of Heat-Transfer Coefficient on the Torsional Stiffness of Thin Wings Subjected to Aerodynamic Heating

A study has been made of the effects of varying the shape, solidity, and heat-transfer coefficient of thin wings with regard to their influence on the torsional-stiffness reduction induced by aerodynamic heating. The variations in airfoil shape include blunting, flattening, and combined blunting and flattening of a solid wing of symmetrical double-wedge cross section. Hollow double-wedge wings of constant skin thickness with and without internal webs also are considered. The effects of heat-transfer coefficients appropriate for laminar and turbulent flow are investigated in addition to a step transition along the chord from a lower to a higher constant value of heat-transfer coefficient. From the results given it is concluded that the flattening of a solid double wedge decreases the reduction in torsional stiffness while slight degrees of blunting increase the loss. The influence of chordwise variations in heat-transfer coefficient due to turbulent and laminar boundary-layer flow on the torsional stiffness of solid wings is negligible. The effect of a step transition in heat-transfer coefficient along the chord of a solid wing can, however, become appreciable. The torsional-stiffness reduction of multiweb and hollow double-wedge wings is substantially less than that calculated for a solid wing subjected to the same heating conditions.

Thomson, Robert G.↗

DOE OSTI2023

Analysis of an Arctic cold air outbreak during autumn and related air mass transformations forced by surface changes and advection in higher altitudes

During marine cold air outbreaks (MCAOs), cold and dry Arctic air masses are transported from the central Arctic southward across the closed sea ice and much warmer open oceans. They experience significant transformations including a rapid heating and moistening, often leading to cloud formation. While intense wintertime MCAOs have been analyzed widely, the air mass transformations during other seasons have been studied sparsely. We address this gap by investigating an MCAO case observed in September 2020. To study the transformation processes, we combine the fifth generation of atmospheric reanalyses of the global climate (ERA5), trajectory calculations, as well as shipborne and airborne measurements. In the central Arctic, observations acquired from aboard the research vessel (RV) Polarstern during the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition characterized the initial state of the air mass over closed sea ice. Trajectories indicated the pathway the air mass took from RV Polarstern southward to the Fram Strait. For the first 24 h of the southbound drift, the air masses remained quasi-stationary. Then, still 15 h ahead of the marginal sea ice zone, differential advection across the boundary layer flow introduced humidity and clouds at higher altitudes between 1.5 and 2.5 km. ERA5-derived temperature and humidity tendencies indicated complex vertical interactions. Radiative cloud-top cooling, entrainment, and turbulence were significantly reduced in the lower and enhanced in the upper advected cloud layer. Eventually, the lower cloud deck dissipated. After this confluence of 2 different air masses, observations gathered by Polar 5 in Fram Strait as part of the MOSAiC Airborne observations in the Central Arctic campaign revealed cloudy, moist layers throughout the lowest 3.5 km and an increasing boundary layer height. Comparing the initial with the final state 48 h later, the largest net heating of +8 K was found close to the surface, yet the largest net moistening of +2.5 g kg -1 at an altitude of 1 km, as the initial profile was exceptionally dry here. We conclude that the observed air mass transformations were driven by the surface changes from sea ice to open ocean but additionally strongly impacted by the differential advection of clouds and moisture across the near-surface MCAO flow.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS2015

Implicit Large-Eddy Simulations of Zero-Pressure Gradient, Turbulent Boundary Layer

A set of direct simulations of zero-pressure gradient, turbulent boundary layer flows are conducted using various span widths (62-630 wall units), to document their influence on the generated turbulence. The FDL3DI code that solves compressible Navier-Stokes equations using high-order compact-difference scheme and filter, with the standard recycling/rescaling method of turbulence generation, is used. Results are analyzed at two different Re values (500 and 1,400), and compared with spectral DNS data. They show that a minimum span width is required for the mere initiation of numerical turbulence. Narrower domains ((is) less than 100 w.u.) result in relaminarization. Wider spans ((is) greater than 600 w.u.) are required for the turbulent statistics to match reference DNS. The upper-wall boundary condition for this setup spawns marginal deviations in the mean velocity and Reynolds stress profiles, particularly in the buffer region.

Turblent boundary layer↗

DOE OSTI2024

Aerodynamic and production comparison of wind farms with downwind versus conventional upwind turbines

Ever-increasing turbine scales and their associated logistical challenges have reignited questions about the performance of downwind rotor configurations. A particular potential benefit of downwind rotor configurations is the farm-scale power increase that may be conferred by tilt-driven downward wake entrainment and associated wake recovery. In this work, a comprehensive aerodynamic analysis is carried out to understand the mechanisms for wake entrainment and recovery across a spectrum of velocity and inflow alignment conditions on a small, structured farm in order to understand the impact of downwind rotors on farm production. The results show that the benefits demonstrated previously in the literature for downwind-rotor farms in aligned flows are fragile, and, outside of strong farm/flow alignment conditions, power production benefits for small farms with downwind rotor configurations are significantly if not completely mitigated by misalignment effects. The work indicates that farm-scale benefits for downwind rotors must be realized either from large-scale entrainment benefits, with more exotic farm arrangements that can take advantage of the aerodynamic effects, or from beneficial fatigue impacts from entrainment of less turbulent outer boundary layer flows.

17 WIND ENERGY↗

NASA NTRS1987

Experimental study of low Reynolds number nozzles

High-performance electrothermal thrusters operate in a low nozzle-throat Reynolds number regime. Under these conditions, the flow boundary layer occupies a large volume inside the nozzle, contributing to large viscous losses. Four nozzles (conical, bell, trumpet, and modified trumpet) and a sharp-edged orifice were evaluated over a Reynolds number range of 500 to 9000 with unheated nitrogen and hydrogen. The nozzles showed significant decreases in specific impulse efficiency with decreasing Reynolds number. At Reynolds numbers less than 1000, all four nozzles were probably filled with a large boundary layer. The discharge coefficient decreased with Reynolds number in the same manner as the specific impulse efficiency. The bell and modified trumpet nozzles had discharge coefficients 4 to 8 percent higher than those of the cone or trumpet nozzles. The Two-Dimensional Kinetics (TDK) nozzle analysis computer program was used to predict nozzle performance. The results were then compared to the experimental results in order to determine the accuracy of the program within this flow regime.

Grisnik, Stanley P.↗

NASA NTRS1987

Experimental study of low Reynolds number nozzles

Grisnik, Stanley P.↗

NASA NTRS1975

Turbulent flow between a rotating disk and a parallel wall

A scalar effective viscosity method is used to calculate the incompressible flow in the gap between two infinite, parallel smooth disks, one of them rotating and the other stationary. The flow is calculated for two values of the gap Reynolds number omega(s squared)/nu, namely 296 and 2852, where omega is the disk angular velocity, s the gap width, and nu the kinematic viscosity of the fluid; and over a range of radial Reynolds number omega(r squared)/nu from 0 to 10 million, where r is the radius from the axis of rotation. A numerical shooting technique is used to match the two boundary-layer flows arising on the disk and wall. Gradual transition is assumed for radial Reynolds numbers from 160,000 to 250,000. The results are in good agreement with Daily and Nece's measurements of velocity profiles and disk friction drag for enclosed rotating disks, and at sufficiently large radius the local skin friction approaches that of turbulent Couette flow.

Cooper, P.↗