Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Turbulence Stabilization”

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 739 records · Page 41

Meteorological and aeronomical requirements for MST radar networks (keynote paper), part 1

Mesosphere - stratosphere - troposphere (MST) radar are phase coherent radars that measure the amplitude and Doppler shift of radio waves that are scattered back to the receiving antennas. For a monostatic system, the line-of-sight projection of the wind vector is obtained if one assumes that the atmospheric scatterers are being swept along with the wind velocity. The three-dimensional wind is then derived either by using multiple beams or by beam swinging. The turbulence intensity is derived either by measuring the backscattered power or by deriving the width of the autocorrelation function for the wind. Furthermore, some information on sharp changes in the atmospheric static stability (e.g., at the tropopause) can be obtained by looking for specular reflections. The discussion addresses the question of how these MST measurement capabilities can contribute to various meteorological and aeronomical research areas.

Geller, M. A.↗

Laminar-Turbulent Transition in High Mach-Number Flows

Accurate prediction of laminar-turbulent transition is critical for the design and characterization of high-speed flight systems. Presently, laminar-turbulent transition prediction in computational fluid dynamics (CFD) is accomplished using PDE-based models coupled to the governing equations of the flow system. Current industrial models suffer significant performance issues when analyzing high-speed systems because they were developed for lower speed flows. Current research focuses on multiple solutions to this problem. The coupling of physics-based stability analysis tools which produce quality transition predictions for high-speed flows to NASA’s OVERFLOW CFD solver is being pursued, as well as the implementation and improvement of emerging models developed specifically for high-speed flow.

High-Speed Flow↗

Logarithmic profile of temperature in sheared and unstably stratified atmospheric boundary layers

The impact of buoyancy on the mean velocity, temperature, and scalar concentration profiles in the lower atmosphere is typically investigated within the framework of Monin-Obukhov similarity theory (MOST). MOST is the theoretical foundation for parametrizing surface-atmosphere exchanges in nearly all weather, climate, and hydrological models.According to MOST, the classic logarithmic profiles of mean velocity and temperature break down as the buoyancy effects become important. However, recent studies on turbulent Rayleigh-Bénard convection and natural convection along vertical walls suggest that the mean temperature in the near-surface region still follows a logarithmic profile.Motivated by these new results, we study the mean potential temperature profile in sheared and unstably stratified atmospheric boundary layers using direct numerical simulations and field observations. We find that the mean potential temperature profile remains logarithmic across a wide range of stability parameters, which characterizes the relative importance of buoyancy versus shear effects. Compared to MOST, our results suggest that the buoyancy force does not modify the logarithmic nature of the mean potential temperature profile, but instead modulates its slope, which is no longer universal and differs from 1/κ,whereκisthe von Kármán constant. Furthermore, this study provides another perspective on scalar turbulence in the atmospheric boundary layer

58 GEOSCIENCES↗

Learning how structures form in drift-wave turbulence

Drift-wave turbulence produces anomalous transport via cross-correlations between fluctuations. This transport has profound implications for confinement, structure formation, and virtually all aspects of the non-linear turbulent dynamics. Herein, we use a data-driven method based on deep learning in order to study turbulent transport in the 2D Hasegawa–Wakatani system and infer a reduced mean-field model from numerical solution. In addition to the usual turbulent diffusion, we find an effect which couples the particle flux to the local gradient of vorticity, which tends to modulate the density profile. The direct coupling to the shear is relatively weak. In addition, the deep learning method finds a model for spontaneous zonal flow generation by negative viscosity, stabilized by non-linear and hyperviscous terms. We compare these results to analytic calculations using quasilinear theory and wave kinetics, finding qualitative agreement, though the calculations miss certain higher-order effects. A simplified, 1-D model for the evolution of the profile, flow, and intensity based on the deep learning results is solved numerically and compared to previous models for staircasing based on bistability. We see that the physics uncovered by the deep learning method provided simple explanations for the formation of zonal structures in the density, flow, and turbulence fields. We highlight the important role of symmetry in the deep learning method and speculate on the portability of the method to other applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A 7.5-Year Dataset of SSM/I-Derived Surface Turbulent Fluxes Over Global Oceans

The surface turbulent fluxes of momentum, latent heat, and sensible heat over global oceans are essential to weather, climate and ocean problems. Wind stress is the major forcing for driving the oceanic circulation, while Evaporation is a key component of hydrological cycle and surface heat budget. We have produced a 7.5-year (July 1987-December 1994) dataset of daily, individual monthly-mean and climatological (1988-94) monthly-mean surface turbulent fluxes over the global oceans from measurements of the Special Sensor Microwave/Imager (SSM/I) on board the US Defense Meteorological Satellite Program F8, F10, and F11 satellites. It has a spatial resolution of 2.0x2.5 latitude-longitude. Daily turbulent fluxes are derived from daily data of SSM/I surface winds and specific humidity, National Centers for Environmental Prediction (NCEP) sea surface temperatures, and European Centre for Medium-Range Weather Forecasts (ECMWF) air-sea temperature differences, using a stability-dependent bulk scheme. The retrieved instantaneous surface air humidity (with a 25-km resolution) IS found to be generally accurate as compared to the collocated radiosonde observations over global oceans. The surface wind speed and specific humidity (latent heat flux) derived from the F10 SSM/I are found to be -encrally smaller (larger) than those retrieved from the F11 SSM/I. The F11 SSM/I appears to have slightly better retrieval accuracy for surface wind speed and humidity as compared to the F10 SSM/I. This difference may be due to the orbital drift of the F10 satellite. The daily wind stresses and latent heat fluxes retrieved from F10 and F11 SSM/Is show useful accuracy as verified against the research quality in si -neasurerrients (IMET buoy, RV Moana Wave, and RV Wecoma) in the western Pacific warm pool during the TOGA COARE Intensive observing period (November 1992-February 1993). The 1988-94 seasonal-mean turbulent fluxes and input variables derived from FS and F11 SSM/Is show reasonable patterns related to seasonal variations of atmospheric general circulation. This dataset of SSM/I-derived turbulent fluxes is useful for climate studies, forcing of ocean models, and validation of coupled ocean-atmosphere global models and can be accessed through the NASA/GSFC Distributed Active Archive Center.

Chou, Shu-Hsien↗

Stability Analysis of the Flow Over a Swept Forward-Facing Step Using PIV Base Flows

Step excrescences are a type of surface imperfection encountered on swept wings of commercial aircraft, commonly because accessibility requirements prevent creating the wing’s surface from a single panel. If the step height is too large, super-critical, the flow can undergo an early transition to turbulence, which can be highly detrimental for the performance of the wing. In the case of a forward-facing step on a swept wing in a low-disturbance environment, stationary crossflow vortices can develop a significant amplitude upstream of the step and hence dominate the structure of the boundary-layer flow over the step. The main goal of the present investigation is to illuminate the path to transition supported by the fascinatingly complex flow field in the direct downstream vicinity of a step with a super critical height. The high-resolution, stereographic Particle Image Velocimetry (PIV) measurement dataset presently available for this flow field provides a complete description of the laminar flow for the execution of BiGlobal stability analysis in a plane parallel to the step. Although the notorious sensitivity of stability results to the description of the base flow demands a very careful uncertainty analysis of those results, it is argued that this very fact can be leveraged to produce new insight into the supported perturbation dynamics. In performing the analysis, several unsteady mode families are discovered that display the explosive perturbation expected for early transition to be induced. In considering domain widths equal to an integer-multiple of the incident crossflow-vortex wavelength and analyzing an extent of 5 crossflow-vortex wavelengths parallel to the step, it is found that the stability results converge while increasing the domain width. It is demonstrated, moreover, that the results for the wider domains can be approximated by appropriately averaging the results on neighboring single-crossflow-vortex-wavelength domains covering the same region. Besides being useful for computational purposes, this observed property suggests interpreting the instability mechanism as a distorted primary mechanism rather than a “proper” secondary mechanism. This follows in the context of the secondary in-stability analysis of three-dimensional boundary layers, because the secondary mechanism is usually characterized by being localized in a pocket of strong shear, while the distorted primary mechanism typically has an infinite support in the direction parallel to the step. Even though the growth rates are found to be sensitive to the interrogation-window size inherent to the PIV post-processing procedure, the spatial structure of the eigenfunctions is found to be relatively insensitive. Lastly, the spatial structure of the eigenfunctions corresponding to all velocity components are matched with the shape functions determined by computing the Spectral Proper Orthogonal Decomposition (SPOD) of a time-resolved measurement of the perturbation content.

forward-facing step↗

Monte Carlo turbulence simulation using rational approximations to von Karman spectra

Turbulence simulation is computationally much simpler using rational spectra, but turbulence falls off as f exp -5/3 in frequency ranges of interest to aircraft response and as predicted by von Karman's model. Rational approximations to von Karman spectra should satisfy three requirements: (1) the rational spectra should provide a good approximation to the von Karman spectra in the frequency range of interest; (2) for stability, the resulting rational transfer function should have all its poles in the left half-plane; and (3) at high frequencies, the rational spectra must fall off as an integer power of frequency, and since the -2 power is closest to the -5/3 power, the rational approximation should roll off as the -2 power at high frequencies. Rational approximations to von Karman spectra that satisfy these three criteria are presented, along with spectra from simulated turbulence. Agreement between the spectra of the simulated turbulence and von Karman spectra is excellent.

Campbell, C. W.↗

Measuring the Effects of Wake Steering With Nacelle-Mounted Scanning Lidar

In large wind plants, wakes from upwind turbines affect downwind turbines by reducing wind speeds and increasing turbulence. Wake steering may mitigate this effect by deflecting the wakes of upwind turbines away from downwind units. Here, we characterize the impact of wakes from yawed turbines at a commercial-scale wind plant under varying atmospheric and turbine operating conditions. Six months of nacelle-based lidar measurements were collected as part of a field campaign in 2019-2020 in the northern US Great Plains to test the effectiveness of wake steering. We separate these lidar scans by atmospheric stability and turbine operating condition to summarize how yawed wind turbine wakes vary with these input parameters in the atmospheric boundary layer. We summarize the impact of wake steering on various wake characteristics including velocity deficit, wake width, and wake center as retrieved from these lidar data. Yawed wakes have significantly different centerline characteristics compared to unyawed wakes, with large regions of meander in the mid wake region. Yawed wakes are also deflected farther than unyawed wakes in less turbulent conditions with velocity deficits persisting further downstream in stable atmospheric conditions. Overall, yawed wind turbine wakes are larger and wider in lower wind speed environments, and yawed wakes are deflected farther in less turbulent conditions, suggesting that wake steering is most effective in stable atmospheric stratification.

17 WIND ENERGY↗

Stability of compressible Taylor-Couette flow

The objectives of this paper are to: (1) develop both analytical and numerical tools that can be used to predict the onset of instability and subsequently to simulate the transition process by which the originally laminar flow evolves into a turbulent flow; and (2) conduct the preliminary investigations with the purpose of understanding the mechanisms of the vortical structures of the compressible flow between tow concentric cylinders.

Kao, K.↗

Conservative discontinuous Galerkin scheme of a gyro-averaged Dougherty collision operator

Here, a conservative discontinuous Galerkin scheme for a non-linear Dougherty collision operator in full-f long-wavelength gyrokinetics is presented. Analytically this model operator has the advective-diffusive form of Fokker-Planck operators, it has a non-decreasing entropy functional, and conserves particles, momentum and energy. Discretely these conservative properties are maintained exactly as well, independent of numerical resolution. In this work the phase space discretization is performed using a novel version of the discontinuous Galerkin scheme, carefully constructed using concepts of weak equality and recovery. Discrete time advancement is carried out with an explicit time-stepping algorithm, whose stability limits we explore. The formulation and implementation within the long-wavelength gyrokinetic solver of Gkeyll are validated with relaxation tests, collisional Landau-damping benchmarks and the study of 5D gyrokinetic turbulence on helical, open field lines.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

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

Satellite-Based Surface Heat Budgets and Sea Surface Temperature Tendency in the Tropical Eastern Indian and Western Pacific Oceans for the 1997/98 El Nino and 1998/99 La Nina

The 1997/98 is a strong El Nino warm event, while the 1998/99 is a moderate La Nina cold event. We have investigated surface heat budgets and sea surface temperature (SST) tendency for these two events in the tropical western Pacific and eastern Indian Oceans using satellite-retrieved surface radiative and turbulent fluxes. The radiative fluxes are taken from the Goddard Satellite-retrieved Surface Radiation Budget (GSSRB), derived from radiance measurements of the Japanese Geostationary Meteorological Satellite 5. The GSSRB covers the domain 40 deg S - 4 deg N, 90 deg E-17 deg W and a period from October 1997 to December 2000. The spatial resolution is 0.5 deg x 0.5 deg lat-long and the temporal resolution is 1 day. The turbulent fluxes are taken from Version 2 of the Goddard Satellite-based Surface Turbulent Fluxes (GSSTF-2). The GSSTF-2 has a spatial resolution of 1 deg x 1 deg lat-long over global Oceans and a temporal resolution of 1 day covering the period July 1987-December 2000. Daily turbulent fluxes are derived from the S S M (Special Sensor Microwave/Imager) surface wind and surface air humidity, and the SST and 2-m air temperature of the NCEP/NCAR reanalysis, using a stability-dependent bulk flux algorithm. The changes of surface heat budgets, SST and tendency, cloudiness, wind speed, and zonal wind stress of the 1997/98 El Nino relative to the1998/99 La Nina for the northern winter and spring seasons are analyzed. The relative changes of surface heat budgets and SST tendency of the two events are quite different between the tropical eastern Indian and western Pacific Oceans. For the tropical western Pacific, reduced solar heating (more clouds) is generally associated with decreased evaporative cooling (weaker winds), and vise versa. The changes in evaporative cooling over-compensate that of solar heating and dominate the spatial variability of the changes in net surface heating. Both solar heating and evaporative cooling offset each other to reduce interannual variability of net surface heating. In addition, the area of increased SST tendency is larger than that of increased net surface heating, due to less solar radiation penetration through the bottom of deeper ocean mixed layer (stronger winds). For the tropical eastern Indian Ocean, enhanced solar heating (less clouds) is generally associated with reduced evaporative cooling (weaker winds). Both solar heating and evaporative cooling reinforce each other to increase interannual variability of net surface heating. In addition, the area of increased SST tendency is smaller than that of increased net surface heating in the southern domain. The relative changes in wind and zonal wind stress indicate more solar radiation penetration through the ocean mixed layer and more northward heat transport by Ocean current from the south to the north Indian Ocean for the El Nino than for the La Nina.

Chou, Shu-Hsien↗

Studying Turbulence Using Numerical Simulation Databases - X Proceedings of the 2004 Summer Program

This Proceedings volume contains 32 papers that span a wide range of topics that reflect the ubiquity of turbulence. The papers have been divided into six groups: 1) Solar Simulations; 2) Magnetohydrodynamics (MHD); 3) Large Eddy Simulation (LES) and Numerical Simulations; 4) Reynolds Averaged Navier Stokes (RANS) Modeling and Simulations; 5) Stability and Acoustics; 6) Combustion and Multi-Phase Flow.

Moin, Parviz↗

Turbulence model reduction by deep learning

A defining problem of turbulence theory is to produce a predictive model for turbulent fluxes. These have profound implications for virtually all aspects of the turbulence dynamics. In magnetic confinement devices, drift-wave turbulence produces anomalous fluxes via cross-correlations between fluctuations. In this work, we introduce an alternative, data-driven method for parametrizing these fluxes. The method uses deep supervised learning to infer a reduced mean-field model from a set of numerical simulations. We apply the method to a simple drift-wave turbulence system and find a significant new effect which couples the particle flux to the local gradient of vorticity. Notably, here, this effect is much stronger than the oft-invoked shear suppression effect. We also recover the result via a simple calculation. The vorticity gradient effect tends to modulate the density profile. In addition, our method recovers a model for spontaneous zonal flow generation by negative viscosity, stabilized by nonlinear and hyperviscous terms. We highlight the important role of symmetry to implementation of the new method.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A computational examination of directional stability for smooth and chined forebodies at high-alpha

Computational Fluid Dynamics (CFD) has been used to study aircraft forebody flowfields at low-speed, angle-of-attack conditions with sideslip. The purpose is to define forebody geometries which provide good directional stability characteristics under these conditions. The flows over the experimentally investigated F-5A forebody and chine type configuration, previously computed by the authors, were recomputed with better grid topology and resolution. The results were obtained using a modified version of CFL3D (developed at NASA Langley) to solve either the Euler equations or the Reynolds equations employing the Baldwin-Lomax turbulence model with the Degani-Schiff modification to account for massive crossflow separation. Based on the results, it is concluded that current CFD methods can be used to investigate the aerodynamic characteristics of forebodies to achieve desirable high angle-of-attack characteristics. An analytically defined generic forebody model is described, and a parametric study of various forebody shapes was then conducted to determine which shapes promote a positive contribution to directional stability at high angle-of-attack. An unconventional approach for presenting the results is used to illustrate how the positive contribution arises. Based on the results of this initial parametric study, some guidelines for aerodynamic design to promote positive directional stability are presented.

Ravi, Ramakrishnan↗

Transition in Turbine Flows

We have further developed our capabilities to analyze transition in turbine boundary layers from first principles by integrating the nonlinear parabolized stability equations (PSE) with improved initial and boundary conditions. With modified iteration schemes, we are able to proceed deeper into the transition region where skin friction coefficient and heat transfer coefficient significantly increase. Initial and boundary conditions at elevated turbulence levels can be derived by receptivity analysis. Test runs for ERCOFTAC test case T3A at 2.4\% turbulence level provide results in good agreement with the experimental data. The sharper minimum of the skin coefficient also shown by DNS results is likely due to the missing intermittency. The method has been applied to various experimentally studied turbine blades (UTRC, VKI, Zierke, Langston, Hippensteele, and others). The PSE results, though physically reasonable, do not agree as well as expected with the experimental findings. We have, therefore, performed an extensive search for the reasons of the seemingly systematic deviations. A first source of uncertainty has been found in the often insufficient documentation of the experiments (e.g. on blockage by end-wall boundary layers). However, variation of the relevant parameters does not lead to more satisfactory agreement. A second reason has been found in the "standard procedure" which considers a 2D flow at midspan and uses a panel code and subsequent boundary-layer code to obtain the laminar basic flow for the transition analysis. Comparison with the pressure distribution obtained with a 3D design code (RVC3D) shows significant three-dimensionality of the flow (e.g. in the UTRC experiments). The spanwise variation has been neglected in our original PSE code. To overcome this problem, we have developed the PSE/3D for fully 3D boundary layers to account for streamwise and spanwise variations. Since the design code does not provide the boundary-layer flow with sufficient resolution, we have generated the Euler solution and employed a 3D boundary-layer code to obtain the viscous basic flow. Although only the linear stability level of PSE/3D has been implemented so far, the discrepancies with the experiments change but do not disappear. We still find deviations between the computed and experimental variations of C(sub f), and St along the blade for laminar flow. The main reason can be seen by comparing the solution of the boundary-layer code with the viscous results of the design code. The conventional boundary-layer solution exhibits an asymptotic behavior appropriate in external aerodynamics but does not match the steep gradients of the inviscid flow through the passage and consequently provides biased results for C(sub f), and St. An attempt is currently being made to correct this deficiency. Before attempting to perform the transition analysis for the viscous flow provided by the design code, we have analyzed the implementation and "best possible" results. Code and results exhibit flaws that may negatively affect the design and are intolerable for transition analysis. Therefore, we have decided to develop a new code to obtain a reliable basis for stability and transition studies. We expect to report improved results by the time of the meeting.

Herbert, Thorwald↗

Two-Phase Turbulence Statistics from High Fidelity Dispersed Droplet Flow Simulations in a Pressurized Water Reactor (PWR) Sub-Channel with Mixing Vanes

In the dispersed flow film boiling regime (DFFB), which exists under post-LOCA (loss-of-coolant accident) conditions in pressurized water reactors (PWRs), there is a complex interplay between droplet dynamics and turbulence in the surrounding steam. Experiments have accredited particular significance to droplet collision with the spacer-grids and mixing vane structures and their consequent positive feedback to the heat transfer recorded in the immediate downstream vicinity. Enabled by high-performance computing (HPC) systems and a massively parallel finite element-based flow solver—PHASTA (Parallel Hierarchic Adaptive Stabilized Transient Analysis)—this work presents high fidelity interface capturing, two-phase, adiabatic simulations in a PWR sub-channel with spacer grids and mixing vanes. Selected flow conditions for the simulations are informed by the experimental data found in the literature, including the steam Reynolds number and collision Weber number (Wec={40,80}), and are characteristic of the DFFB regime. Data were collected from the simulations at an unprecedented resolution, which provides detailed insights into the continuous phase turbulence statistics, highlighting the effects of the presence of droplets and the comparative effect of different Weber numbers on turbulence in the surrounding steam. Further, axial evolution of droplet dynamics was analyzed through cross-sectionally averaged quantities, including droplet volume, surface area and Sauter mean diameter (SMD). The downstream SMD values agree well with the existing empirical correlations for the selected range of Wec. The high-resolution data repository from the simulations herein is expected to be of significance to guide model development for system-level thermal hydraulic codes.

Saini, Nadish↗

Hypersonic Boundary-Layer Instabilities over Ogive-Cylinder Models

Computational investigations of an ogive-cylinder geometry with varying nosetips at zero degrees angle of attack are presented. The model geometry and conditions are selected to match experiments conducted in the Air Force Research Laboratory (AFRL) Mach 6 Ludwieg Tube. Five nosetips of interest were selected for the computational studies herein: two sharp ogives, two blunt ogives, and one hemispherical nosetip. Computations are performed at a freestream Reynolds number of 7.01 × 10 6 m -1 . Each sharp and blunt tip ogives had a 14 and 28 degree version for the tip angles. A cylindrical section follows the nosetip, resulting in a meter long model such as in the experiments. The laminar flow solutions are analyzed. The boundary-layer-edge properties and the velocity and temperature profiles are compared across streamwise locations aft of the ogive-cylinder junction. The blunt nosetips induce an entropy layer that envelopes the boundary-layer profiles. Modal stability analysis identifies most amplified frequencies corresponding to Mack’s second modes that agree with experimental results for the sharp tips. Similar to the experimental measurements based on wall-mounted pressure sensors, no unstable modes are found for the blunt models. Nonmodal analysis revealed a broadband set of disturbances present for the blunter tips, in agreement with experimental observations. Flow perturbation contours of most amplified planar and oblique disturbances are shown to qualitatively match wind tunnel Schlieren images, with a switch from rope-like to elongated structures, i.e., from high frequency Mack’s second modes to low frequency Mack’s first modes, as the nosetip angle is increased for the sharp tip, and from boundary-layer to entropy-layer disturbances as the bluntness is increased.

Boundary layer transition↗