Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “turbulence methods”

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 775 records · Page 43

What can the GRACE seasonal cycle tell us about lake-aquifer interactions?

Lake-aquifer interactions have been the subject of investigation and debate for decades. Traditional investigation methods include measurement of water flux across the groundwater-surface water interface, application of heat and environmental tracer methods, conducting numerical simulations of the water flow, and mass balance-based approaches. We first review, evaluate and describe the different approaches that have been applied for examining lake-aquifer interactions and propose an additional complementary approach. While the leakage from lakes and its slow movement through porous media (laminar flow) is well characterized by existing methods, modeling the rapid and turbulent flow through preferred pathways (faults/karst) remains a challenging task. To better understand the nature, and the full scale, of the lake aquifer interactions, and using Lake Nasser and the underlying Dakhla subbasin of the Nubian Sandstone Aquifer System (NSAS) in the Western Desert of Egypt (area: 0.66 x 10 6 km 2 ) as a test site, we examine, model, and correlate temporally and spatially, the variations in Gravity Recovery and Climate Experiment terrestrial water storage (GRACE TWS ) solutions, with precipitation, lake levels, area, and water volume. We review current GRACE applications in hydrology and present our novel approach. Findings include: firstly, large seasonal fluctuations (peak: Nov./Dec.; trough: July/Aug.) in Lake Nasser's surface water levels (average annual fluctuations: from 3.7 m to 7.5 m), area (average area: 3622 km 2 but up to 4530 km 2 ), and volume (annual average: 13.4 km 3 , and up to 34 km 3 ) are observed throughout years 2006 to 2015. These fluctuations are accompanied by an increase in GRACE TWS (average: 50 ± 13 mm/yr, up to 77 ± 18 mm/yr) over Lake Nasser and by a progression of a front of increasing GRACE TWS values (> 50 ± 13 mm) away from the lake reaching distances of up to 700 km some 3 to 5 months following peak lake level periods. The areas witnessing the seasonal increase in GRACE TWS display a progressive increase in phase and decrease in amplitude with distance from the lake. Secondly, the negligible precipitation over the Dakhla subbasin cannot account for the observed seasonal GRACE TWS patterns and neither can the leakage signal from Lake Nasser. Thirdly, overall similarities in the distribution of GRACE TWS seasonal spatial patterns are observed. During high lake level periods (e.g., Dec. 2007; Feb. 2015) additional, and more distant, areas from Lake Nasser saw increases in GRACE TWS and vice versa during low lake level periods (e.g., July 2006; Aug. 2010). These observations are consistent with Lake Nasser being the main source of modern recharge for the Dakhla subbasin and suggest a new conceptual model for the subbasin: a slow groundwater flow through a porous matrix and a rapid groundwater flow along a network of faults, fractures, and karst topography across the matrix. We suggest that the proposed conceptual model for the interaction between Lake Nasser and the Dakhla subbasin could be applicable to aquifers of similar geologic, climatic, and hydrologic settings worldwide and that approaches similar to those advocated here could be used to investigate the validity of this suggestion.

54 ENVIRONMENTAL SCIENCES↗

Alternative formulation of weak magnetohydrodynamic turbulence theory

In a recent paper, the weak turbulence theory for incompressible magnetohydrodynamics is formulated by employing the method customarily applied in the context of kinetic weak plasma turbulence theory. Such an approach simplified certain mathematical procedures including achieving the closure relationship. The formulation in the above-cited paper starts from the equations of incompressible magnetohydrodynamic (MHD) theory expressed via Elsasser variables. The derivation of nonlinear wave kinetic equation therein is obtained via a truncated solution at the second-order of iteration following the standard practice. In the present paper, the weak MHD turbulence theory is alternatively formulated by employing the pristine form of incompressible MHD equation rather than that expressed in terms of Elsasser fields. The perturbative expansion of the nonlinear momentum equation is carried out up to the third-order iteration rather than imposing the truncation at the second order. Furthermore, it is found that while the resulting wave kinetic equation is identical to that obtained in the previous paper cited above, the third-order nonlinear correction plays an essential role for properly calculating derived quantities such as the total and residual energies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The Effects of Channel Supplies on Overall Film-Cooling Effectiveness

Abstract Cooling components in the hot section of a gas turbine are essential to component durability. Common methods of cooling include rib turbulators in internal passages and film cooling on external surfaces. The holes that produce the film cooling are fed from the internal channels often containing ribs. Consequently, there is an interdependence of internal heat transfer and external film cooling. The purpose of this study was to obtain a better understanding of the interaction of ribs and film cooling. To quantify the cooling performance, the surface temperatures were measured, from which overall effectiveness was calculated. For the experiments, additively manufactured test coupons were made of Inconel 718 to match engine Biot numbers. These test coupons had internal feed channels with and without ribs and had both cylindrical holes and meter-diffuser-shaped holes with 15 deg lateral expansion angles and a 1 deg forward expansion angle. A single rectangular channel was one type of feed channel. The other type of feed channel was individual circular channels, with each circular channel supplying an individual film-cooling hole. The experimental results showed that the circular individual channels have 80% higher baseline overall effectiveness than the single rectangular channels without any film cooling. Ribbed turbulators without film cooling also increased the overall effectiveness by 21% for single rectangular channels and by 29% for the circular individual channels compared to the respective non-ribbed channels. While the film cooling increased the overall effectiveness of all geometries, the single rectangular channels had increased overall effectiveness levels by up to twice that of the no film-cooling case. On average, the single rectangular channels had an 80% improvement from film cooling, whereas the individual channel feeds, on average, had only a 50% improvement, given their high baseline effectiveness levels.

Engineering↗

kokkos-fft: A shared-memory FFT for the Kokkos ecosystem

kokkos-fft provides a unified, performance-portable interface for Fast Fourier Transforms (FFTs) within the Kokkos ecosystem (C. Trott et al., 2021). It seamlessly integrates with leading local FFT libraries including FFTW, cuFFT, rocFFT, and oneMKL. Designed for simplicity and efficiency, kokkos-fft offers a user experience akin to numpy.fft for in-place and out-of-place transforms, while leveraging the raw speed of vendor-optimized libraries. A demonstration solving 2D Hasegawa-Wakatani turbulence with the Fourier spectral method illustrates how kokkos-fft can deliver significant speedups over Python-based alternatives without drastically increasing code complexity, empowering researchers to perform high-performance FFTs simply and effectively.

97 MATHEMATICS AND COMPUTING↗

Local correlation-based transition models for high-Reynolds-number wind-turbine airfoils

Abstract. Modern wind-turbine airfoil design requires robust performance predictions for varying thicknesses, shapes, and appropriate Reynolds numbers. The airfoils of current large offshore wind turbines operate with chord-based Reynolds numbers in the range of 3–15 million. Turbulence transition in the airfoil boundary layer is known to play an important role in the aerodynamics of these airfoils near the design operating point. While the lack of prediction of lift stall through Reynolds-averaged Navier–Stokes (RANS) computational fluid dynamics (CFD) is well known, airfoil design using CFD requires the accurate prediction of the glide ratio (L/D) in the linear portion of the lift polar. The prediction of the drag bucket and the glide ratio is greatly affected by the choice of the transition model in RANS CFD of airfoils. We present the performance of two existing local correlation-based transition models – one-equation model (γ− SA) and two-equation model (γ-Reθt‾- SA) coupled with the Spalart–Allmaras (SA) RANS turbulence model – for offshore wind-turbine airfoils operating at a high Reynolds number. We compare the predictions of the two transition models with available experimental and CFD data in the literature in the Reynolds number range of 3–15 million including the AVATAR project measurements of the DU00-W-212 airfoil. Both transition models predict a larger L/D compared to fully turbulent results at all Reynolds numbers. The two models exhibit similar behavior at Reynolds numbers around 3 million. However, at higher Reynolds numbers, the one-equation model fails to predict the natural transition behavior due to early transition onset. The two-equation transition model predicts the aerodynamic coefficients for airfoils of various thickness at higher Reynolds numbers up to 15 million more accurately compared to the one-equation model. As a result, the two-equation model predictions are more comparable to the predictions from eN transition model. However, a limitation of this model is observed at very high Reynolds numbers of around 12–15 million where the predictions are very sensitive to the inflow turbulent intensity. The combination of the two-equation transition model coupled with the Spalart–Allmaras (SA) RANS turbulence model is a good method for performance prediction of modern wind-turbine airfoils using CFD.

17 WIND ENERGY↗

Review of the Air Force summer study program on the integration of wind tunnels and computers

The present state of computational fluid dynamics and its impact on the design cycle and computer requirements for future developments in this field were explored. The increase in productivity and efficiency which experimental facilities can achieve by a close integration with computers was investigated together with possible improvements in simulation quality of wind tunnels in conjunction with computer control. Research experiments are outlined to provide a better understanding of the physics of fluid flow and to assist in the modeling of these phenomena for computational methods, with primary emphasis on turbulent flows.

Marschner, B. W.↗

Large-eddy simulation of a turbulent mixing layer

The three dimensional, time dependent (incompressible) vorticity equations were used to simulate numerically the decay of isotropic box turbulence and time developing mixing layers. The vorticity equations were spatially filtered to define the large scale turbulence field, and the subgrid scale turbulence was modeled. A general method was developed to show numerical conservation of momentum, vorticity, and energy. The terms that arise from filtering the equations were treated (for both periodic boundary conditions and no stress boundary conditions) in a fast and accurate way by using fast Fourier transforms. Use of vorticity as the principal variable is shown to produce results equivalent to those obtained by use of the primitive variable equations.

Mansour, N. N.↗

Calculations of transonic flow about an airfoil in a wind tunnel

A combined experimental and numerical study was performed to include wind-tunnel wall interference effects in calculations for airfoil flows at transonic speeds. Pressure-survey-tube and laser-Doppler velocimeter measurements were made in the flow field about an airfoil in the 2- by 2-Foot Transonic Wind Tunnel at Ames Research Center. The results were then used as boundary data in a Navier-Stokes code modified by incorporating a pressure condition on the upper and lower computational boundaries. Comparison of calculated results and experimental data obtained from the surface of the airfoil indicates that the pressure-boundary condition is particularly effective in moving the shock to a position near that observed experimentally when the flow remains attached. For flows with large separation, shock position and viscous-layer properties are not well predicted, principally because of the inadequacies of the algebraic turbulence models employed with the method.

King, L. S.↗

Comments on settling chamber design for quiet, blowdown wind tunnels

Transfer of an existing continous circuit supersonic wind tunnel to Langley and its operation there as a blowdown tunnel is planned. Flow disturbance requirements in the supply section and methods for reducing the high level broad band acoustic disturbances present in typical blowdown tunnels are reviewed. Based on recent data and the analysis of two blowdown facilities at Langley, methods for reducing the total turbulence levels in the settling chamber, including both acoustic and vorticity modes, to less than one percent are recommended. The pertinent design details of the damping screens and honeycomb and the recommended minimum pressure drop across the porous components providing the required two orders of magnitude attenuation of acoustic noise levels are given. A suggestion for the support structure of these high pressure drop porous components is offered.

Beckwith, I. E.↗

Boundary layer transition effects on flow separation around V/STOL engine inlets at high incidence

Numerical methods for calculating laminar and turbulent boundary layer development around vertical-short take off and landing engine inlets at high incidence angles are investigated. Various transition models were compared and evaluated in calculations off flow separation bound inside the inlet. Results of the transition effects on the boundary layer characteristics at onset of separation for two types of engine inlet geometries are presented. Some of the numerical results are compared with existing wind-tunnel test data for scaled inlet models to demonstrate the effects of transition models in the numerical scheme. The effects of transition modeling on the boundary layer development are illustrated for typical engine operating conditions.

Chou, D. C.↗

An interplanetary magnetic field ensemble at 1 AU

A method for calculation ensemble averages from magnetic field data is described. A data set comprising approximately 16 months of nearly continuous ISEE-3 magnetic field data is used in this study. Individual subintervals of this data, ranging from 15 hours to 15.6 days comprise the ensemble. The sole condition for including each subinterval in the averages is the degree to which it represents a weakly time-stationary process. Averages obtained by this method are appropriate for a turbulence description of the interplanetary medium. The ensemble average correlation length obtained from all subintervals is found to be 4.9 x 10 to the 11th cm. The average value of the variances of the magnetic field components are in the approximate ratio 8:9:10, where the third component is the local mean field direction. The correlation lengths and variances are found to have a systematic variation with subinterval duration, reflecting the important role of low-frequency fluctuations in the interplanetary medium.

Matthaeus, W. H.↗

SSME Turbopump Turbine Computations

A two-dimensional viscous code was developed to be used in the prediction of the flow in the SSME high-pressure turbopump blade passages. The rotor viscous code (RVC) employs a four-step Runge-Kutta scheme to solve the two-dimensional, thin-layer Navier-Stokes equations. The Baldwin-Lomax eddy-viscosity model is used for these turbulent flow calculations. A viable method was developed to use the relative exit conditions from an upstream blade row as the inlet conditions to the next blade row. The blade loading diagrams are compared with the meridional values obtained from an in-house quasithree-dimensional inviscid code. Periodic boundary conditions are imposed on a body-fitted C-grid computed by using the GRAPE GRids about Airfoils using Poisson's Equation (GRAPE) code. Total pressure, total temperature, and flow angle are specified at the inlet. The upstream-running Riemann invariant is extrapolated from the interior. Static pressure is specified at the exit such that mass flow is conserved from blade row to blade row, and the conservative variables are extrapolated from the interior. For viscous flows the noslip condition is imposed at the wall. The normal momentum equation gives the pressure at the wall. The density at the wall is obtained from the wall total temperature.

Jorgenson, P. G. E.↗

Boundary layers on oscillating airfoils

A two-point finite-difference unsteady laminar and turbulent boundary-layer computational method has been used to investigate the properties of the flow around an airfoil (NACA 0012) oscillating through angles of attack up to 18 degrees, and for reduced frequencies of 0.01 and 0.20. The unsteady potential flow was determined using the method of Geissler. The influence of transition location on stall behavior has been investigated, using both experimentally determined transition information, and transition located at the pressure peak; the results show the need for viscous-inviscid interaction in future computations of such flows.

Carr, L. W.↗

An interplanetary magnetic field ensemble at 1 AU

A method for calculation ensemble averages from magnetic field data is described. A data set comprising approximately 16 months of nearly continuous ISEE-3 magnetic field data is used in this study. Individual subintervals of this data, ranging from 15 hours to 15.6 days comprise the ensemble. The sole condition for including each subinterval in the averages is the degree to shich it represents a weakly time-stationary process. Averages obtained by this method are appropriate for a turbulence description of the interplanetary medium. The ensemble average correlation length obtained from all subintervals is found to be 4.9 x 10 to the 11th cm. The average value of the variances of the magnetic field components are in the approximate ratio 8:9:10, where the third component is the local mean field direction. The correlation lengths and variances are found to have a systematic variation with subinterval duration, reflecting the important role of low-frequency fluctuations in the interplanetary medium.

Matthaeus, W. H.↗

Further results of natural laminar flow flight test experiments

Flight test experiments were conducted to measure the extent and nature of natural laminar flow on a smoothed test region of a swept-wing business jet wing. Surface hot film anemometry and sublimating chemicals were used for transition detection. Surface pressure distributions were measured using pressure belts. Engine noise was monitored by a microphone attached to the wing surface to study possible acoustic effects on stability of the laminar boundary layer. Side-slip conditions were flown to simulate changes in effective wing sweep. Flight instrumentation and ground data analysis techniques and a method for measuring intermittency of turbulence are described. Correlation was obtained between the hot film gage signals and chemicals for transition detection. Cross-flow vortices were observed for some flight conditions. Results of spectral and statistical analysis of the hot film signals for various flight test conditions are presented.

Wentz, W. H., Jr.↗

Heat and mass transfer in flames

Heat- and mass-transfer processes in turbulent diffusion flames are discussed, considering turbulent mixing and the structure of single-phase flames, drop processes in spray flames, and nonluminous and luminous flame radiation. Interactions between turbulence and other phenomena are emphasized, concentrating on past work of the author and his associates. The conserved-scalar formalism, along with the laminar-flamelet approximation, is shown to provide reasonable estimates of the structure of gas flames, with modest levels of empiricism. Extending this approach to spray flames has highlighted the importance of drop/turbulence interactions; e.g., turbulent dispersion of drops, modification of turbulence by drops, etc. Stochastic methods being developed to treat these phenomena are yielding encouraging results.

Faeth, G. M.↗

Method for laminar boundary layer transition visualization in flight

Disclosed is a method of visualizing laminar to turbulent boundary layer transition, shock location, and laminar separation bubbles around a test surface. A liquid crystal coating is formulated using an unencapsulated liquid crystal operable in a temperature bandwidth compatible with the temperature environment around the test surface. The liquid crystal coating is applied to the test surface, which is preferably pretreated by painting with a flat, black paint to achieve a deep matte coating, after which the surface is subjected to a liquid or gas flow. Color change in the liquid crystal coating is produced in response to differences in relative shear stress within the boundary layer around the test surface. The novelty of this invention resides in the use of liquid crystals which are sensitive to shear stress to show aerodynamic phenomena such as a boundary layer transition, shock location, and laminar separation bubbles around a test surface.

Bruce J Holmes↗