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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
NASA NTRS2020

Tuning Neural Network Models for Improved Prediction of Boundary Layer Transition

Boundary layer transition can strongly impact flight vehicle performance as it influences surface skin friction and aerodynamic heating, making accurate transition prediction a key to designing next generation aircraft. Artificial neural networks (ANNs) have shown promise toward predicting laminar-turbulent transition based on linear stability correlations. The computational efficiency of ANNs and the substantially reduced user involvement in relation to direct computations based on the linear stability theory (LST) makes them an attractive methodology for integrating the LST based correlations in computational fluid dynamics codes. Tollmien-Schlichting (TS) waves correspond to the dominant transition mechanism in 2D or weakly 3D subsonic boundary layers, such as those encountered in general aviation applications. Improvements to neural network model accuracy in predicting the amplification rates of TS instability waves have been investigated by leveraging recent machine learning developments in conjunction with surrogate optimization techniques and via suitable augmentation of the data used to train the networks. The optimized models trained on the modified dataset reduced the average transition location errors on different airfoils at several flow conditions by 51% of the original manually-tuned network’s errors on the same flow cases. The actual transition locations were derived from the Langley Stability and Transition Analysis Code (LASTRAC).

Machine Learning↗
NASA NTRSDate not supplied

An Observationally-Based Determination of the Arctic Sea Ice-Cloud Feedback Since 2000: Isolating the Arctic Cloud Response to Sea Ice Loss

Arctic sea ice responds to and drives Arctic climate change. The interactions between Arctic sea ice and clouds represent a mechanism through which sea ice can drive climate change. We composite active remote sensing satellite cloud properties for ice-free, marginal ice zone (MIZ), and ice-covered surfaces during MIZ crossing events to investigate the influence of the transition from an ice-covered to an ice-free surface on low-level clouds. We demonstrate that the event-based methodology controls for large-scale meteorological factors and isolates the sea ice effect on clouds. We find larger cloud fraction and total water content below ~1.5 km over ice-free relative to ice-covered surfaces during non-summer months, indicating a low-level cloud sensitivity to Arctic sea ice decline. During summer, results show larger cloud fraction and water content over ice-free surfaces, however the differences are statistically indistinguishable. Evidence is provided that atmospheric thermodynamic profile differences cause the cloud property differences, namely that ice-free footprints are warmer, moister, have more positive surface turbulent fluxes and are less stable than their ice-covered counterparts. Ice-free and ice-covered surface cloud property differences scale with surface temperature differences such that cloud property differences are only found in the presence of a surface temperature difference. We conclude that surface temperature differences modulate the cloud response to sea ice loss through influences on surface turbulent fluxes and lower tropospheric stability. The results imply a positive non-summer sea ice-cloud feedback and that up to a 0.02 cloud fraction and 0.05 g m 3 total water content increase in fall are due to the observed sea ice decline.

Patrick C Taylor↗
NASA NTRSDate not supplied

Estimating the Arctic Cloud-Sea Ice Feedback With Observations during the EOS Period

Arctic sea ice responds to and drives Arctic climate change. The interactions between Arctic sea ice and clouds represent a mechanism through which sea ice can drive climate change. We composite active remote sensing satellite cloud properties for ice-free, marginal ice zone (MIZ), and ice-covered surfaces during MIZ crossing events to investigate the influence of the transition from an ice-covered to an ice-free surface on low-level clouds. We demonstrate that the event-based methodology controls for large-scale meteorological factors and isolates the sea ice effect on clouds. We find larger cloud fraction and total water content below ~1.5 km over ice-free relative to ice-covered surfaces during non-summer months, indicating a low-level cloud sensitivity to Arctic sea ice decline. During summer, results show larger cloud fraction and water content over ice-free surfaces, however the differences are statistically indistinguishable. Evidence is provided that atmospheric thermodynamic profile differences cause the cloud property differences, namely that ice-free footprints are warmer, moister, have more positive surface turbulent fluxes and are less stable than their ice-covered counterparts. Ice-free and ice-covered surface cloud property differences scale with surface temperature differences such that cloud property differences are only found in the presence of a surface temperature difference. We conclude that surface temperature differences modulate the cloud response to sea ice loss through influences on surface turbulent fluxes and lower tropospheric stability. The results imply a positive non-summer sea ice-cloud feedback and that up to a 0.02 cloud fraction and 0.05 g m 3 total water content increase in fall are due to the observed sea ice decline.

Patrick C. Taylor↗
NASA NTRS1991

Multigrid solution of compressible turbulent flow on unstructured meshes using a two-equation model

The system of equations consisting of the full Navier-Stokes equations and two turbulence equations was solved for in the steady state using a multigrid strategy on unstructured meshes. The flow equations and turbulence equations are solved in a loosely coupled manner. The flow equations are advanced in time using a multistage Runge-Kutta time stepping scheme with a stability bound local time step, while the turbulence equations are advanced in a point-implicit scheme with a time step which guarantees stability and positively. Low Reynolds number modifications to the original two equation model are incorporated in a manner which results in well behaved equations for arbitrarily small wall distances. A variety of aerodynamic flows are solved for, initializing all quantities with uniform freestream values, and resulting in rapid and uniform convergence rates for the flow and turbulence equations.

Mavriplis, D. J.↗
NASA NTRS1950

Turbulence and heat stratification

With applications to meteorology in view, the method of small oscillations is used to investigate the stability of a plane laminar boundary layer in which the density depends in a special way on the distance from the plate. It is found that for a constant Froude number the flow becomes more stable as the Richardson number increases until, at a critical value of the Richardson number that depends on the Froude number, the flow is stable for all disturbance wave lengths and Reynolds numbers. The calculated results are verified by experiment.

FLOW VELOCITY - DISTRIBUTION↗
NASA NTRS1957

Base Pressure at Supersonic Speeds on Two-dimensional Airfoils and on Bodies of Revolution with and Without Fins Having Turbulent Boundary Layers

An analysis has been made of available experimental data to show the effects of most of the variables that are more predominant in determining base pressure at supersonic speeds. The analysis covers base pressures for two-dimensional airfoils and for bodies of revolution with and without stabilizing fins and is restricted to turbulent boundary layers. The present status of available experimental information is summarized as are the existing methods for predicting base pressure. A simple semiempirical method is presented for estimating base pressure. For two-dimensional bases, this method stems from an analogy established between the base-pressure phenomena and the peak pressure rise associated with the separation of the boundary layer. An analysis made for axially symmetric flow indicates that the base pressure for bodies of revolution is subject to the same analogy. Based upon the methods presented, estimations are made of such effects as Mach number, angle of attack, boattailing, fineness ratio, and fins. These estimations give fair predictions of experimental results. (author)

LOVE EUGENE S↗
NASA NTRS1987

Prediction of transition on a flat plate under the influence of free-stream turbulence using low-Reynolds-number two-equation turbulence models

The prediction characteristics and capabilities of two popular two-equation low-Reynolds-number turbulence models (Jones, Launder 1972; Lam, Bremhorst 1981) have been evaluated with respect to the prediction of transition of a flat plate under the influence of free-stream turbulence. The sensitivity of the predictions to free-stream turbulence intensity, initial starting location of the calculation, and the assumed initial starting profiles for k and epsilon has been determined and presented. Although both models predict the correct qualitative characteristics of transition, they also exhibit significant quantitative deficiencies with regard to both the predicted location and the length of transition. A modification to the production term in the turbulent kinetic energy equation is proposed which is based on a simple stability criterion and correlated to the free-stream turbulence level. The modification becomes inactive in the fully turbulent regime, but is shown to improve both the qualitative and quantitative characteristics of the transition predictions.

Schmidt, Rodney C.↗
NASA NTRS1980

A piloted simulator investigation of static stability and stability/control augmentation effects on helicopter handling qualities for instrument approach

A ground simulator experiment was conducted on the Flight Simulator for Advanced Aircraft at Ames Research Center to investigate the influence of several static stability and stability/control augmentation design parameters on helicopter flying qualities during terminal area operations in instrument conditions. Effects of light turbulence were included. Two levels of static stability in each rotational axis (pitch, roll, yaw) were examined for a hingeless rotor configuration. The variations in pitch and roll were: (1) stable and (2) neutral static stability; in yaw there were two stable levels. Four types of stability/control augmentation were also examined for the lower level of static stability in each axis. This latter investigation covered three helicopter rotor types: hingeless, articulated, and teetering. Four pilots performed a total of 105 evaluations of these parameters for a representative VOR instrument approach task. Pilot rating results indicate the acceptability of neutral static stability longitudinally and laterally and the need for pitch-roll attitude augmentation to achieve a satisfactory system.

Lebacqz, J. V.↗
NASA NTRS1986

Jeans criterion in a turbulent medium

According to the classical Jeans analysis, all the molecular clouds of mass larger than a few 100 M(solar), size larger than about 1pc and kinetic temperature Tk less than 30K are gravitationally unstable. We have shown that in clouds supported by internal supersonic motions, local gravitational instabilities may appear within molecular clouds which are globally stable. The argument is threefold: (1) when the turbulent kinetic energy is included into the internal energy term, the virial equilibrium condition shows that molecular clouds such as those observed, which are gravitationally unstable according to the Jeans criterion, are indeed globally stable if supported by a turbulent velocity field of power spectrum steeper than 3; (2) 2D compressible hydrodynamical simulations show that a supersonic turbulent velocity field generates a turbulent pressure within clouds, the gradients of which stabilize the unstable scales (i.e., the largest scales and the cloud itself) against gravitational collapse; (3) an analysis similar to the Jeans approach but including the turbulent pressure gradient term, gives basically the same results as those given in (1). Clouds of mean density lower than a critical value are found to be stable even though more massive than their Jeans mass. In clouds of mean density larger than that critical value, the gravitational instability appears only over a range of scales smaller than the cloud size, the largest scales being stable. In practice, the observed mean densities are lower than this critical value: the observation of a small number of cores and stars of a few solar masses embedded in clouds of several hundred solar masses can only be understood in terms of small scale density fluctuations of large amplitude generated by the supersonic turbulence which would occasionally overtake the limit of gravitational stability.

Bonazzola, S.↗
DOE OSTI2025

Multi-scale Interaction Mechanism for Edge-Localized-Mode Suppression in the Tokamak Edge

A central challenge in fusion energy is reconciling the high-confinement mode required for reactor performance with the intense intermittent relaxation events it produces, known as edge-localized modes. These instabilities arise in the steep pressure pedestal at the plasma edge when magnetohydrodynamic thresholds are crossed, inflicting damaging heat loads on reactor components. Here, we show that multiscale interactions between microscopic turbulence and macroscopic magnetohydrodynamic modes provide encouraging prospects for self-organized edge-localized modes regulation. Using direct quantitative measurements of multiscale modes, eddy dynamics, and turbulent flux, we show that small-scale electron drift wave turbulence actively scatters the large-scale peeling-ballooning modes. This scattering decorrelates the pressure and velocity fields of the instability, so arresting its growth. Our modeling and theoretical analysis confirm this suppression mechanism is effective even when conventional linear stability thresholds are exceeded. This work establishes a nonlinear principle for edge-localized modes stability, revealing how ambient micro-turbulence can be leveraged to maintain a macro-stable, high-performance pedestal for future fusion reactors.

experimental nuclear physics↗
NASA NTRS1954

Wind-tunnel Tests of a Model of a Wingless Fin-controlled Missile to Obtain Static Stability and Control Characteristics Through a Range of Mach Numbers from 0.5 to 0.88

An investigation at medium to high subsonic speeds has been conducted in the Langley low-turbulence pressure tunnel to determine the static stability and control characteristics and to measure the fin normal forces and moments for a model of a wingless fin-controlled missile. The data were obtained at Reynolds number of 2.1 x 10(6) based on the missile maximum diameter or 17.7 x 10(6) based on missile length; this Reynolds number was found to be large enough to avoid any large scale effects between the test and the expected flight Reynolds number. With the horizontal-fin deflection limited to a maximum of 6 degrees, longitudinally stable and trimmed flight could not be maintained beyond an angle of attack of 17 degrees for a Mach number of 0.88 and beyond 20 degrees for a Mach number of 0.50 for any center-of-gravity location without the use of some auxiliary stability or control device such as jet vanes. Mach number had no appreciable effect on the center-of-pressure positions and only a slight effect on neutral-point position. There was a shift in neutral-point position of about 1 caliber as the angle of attack was varied through the range for which the neutral point could be determined. Yawing the model to angles of sideslip up to 7 degrees had little effect on the longitudinal stability at angles of attack up to 15 degrees; however, above 15 degrees, the effect of sideslip was destabilizing. With the vertical fins at a plus-or-minus 6 degree roll deflection, the rolling moment caused by yawing the model at high angles of attack could be trimmed out up to angles of sideslip of 6.5 degrees and an angle of attack of 26 degrees for a Mach number of 0.50; this range of sideslip angles was reduced to 3 degrees at a Mach number of 0.88. The data indicated that, at lower angles of attack, the trim range extended to higher angles of sideslip. The total normal-force and hinge-moment coefficients for both horizontal fins were slightly nonlinear with both angle-of-attack and fin deflection. The effect of Mach number was to reduce the slopes of the hinge-moment coefficient with angle of attack and deflection angle. In general, the effort of increasing the sideslip angle was to reduce the values of the fin normal-force and hinge-moment coefficients.

Burrows, Dale L↗
NASA NTRS1994

Modification of the Two-equation Turbulence Model in NPARC to a Chien Low Reynolds Number K-epsilon Formulation

This report documents the changes that were made to the two-equation k-epsilon turbulence model in the NPARC (National-PARC) code. The previous model based on the low Reynolds number model of Speziale, was replaced with the low Reynolds number k-epsilon model of Chien. The most significant difference was in the turbulent Prandtl numbers appearing in the diffusion terms of the k and epsilon transport equations. A new inflow boundary condition and stability enhancements were also implemented into the turbulence model within NPARC. The report provides the rationale for making the change to the Chien model, code modifications required, and comparisons of the performances of the new model with the previous k-epsilon model and algebraic models used most often in PARC/NPARC. The comparisons show that the Chien k-epsilon model installed here improves the capability of NPARC to calculate turbulent flows.

Georgiadis, Nicholas J.↗
NASA NTRS1997

Modeling of Turbulent Swirling Flows

Aircraft engine combustors generally involve turbulent swirling flows in order to enhance fuel-air mixing and flame stabilization. It has long been recognized that eddy viscosity turbulence models are unable to appropriately model swirling flows. Therefore, it has been suggested that, for the modeling of these flows, a second order closure scheme should be considered because of its ability in the modeling of rotational and curvature effects. However, this scheme will require solution of many complicated second moment transport equations (six Reynolds stresses plus other scalar fluxes and variances), which is a difficult task for any CFD implementations. Also, this scheme will require a large amount of computer resources for a general combustor swirling flow. This report is devoted to the development of a cubic Reynolds stress-strain model for turbulent swirling flows, and was inspired by the work of Launder's group at UMIST. Using this type of model, one only needs to solve two turbulence equations, one for the turbulent kinetic energy k and the other for the dissipation rate epsilon. The cubic model developed in this report is based on a general Reynolds stress-strain relationship. Two flows have been chosen for model evaluation. One is a fully developed rotating pipe flow, and the other is a more complex flow with swirl and recirculation.

Shih, Tsan-Hsing↗
DOE OSTI2025

Simulated meteorological impacts of offshore wind turbines and sensitivity to the amount of added turbulence kinetic energy

Offshore wind energy projects are currently in development off the east coast of the United States and may influence the local meteorology of the region. Wind power production and other commercial uses in this area are related to atmospheric conditions, and so it is important to understand how future wind plants may change the local meteorology. In the absence of measurements of potential wind plant impacts on meteorology, simulations offer the next-best possible insight into wake effects on boundary layer height, temperature, fluxes, and wind speeds. However, simulation tools that capture these effects offer multiple options for representing the amount of turbine-added turbulence that may impact assessments of micrometeorological effects. To explore this sensitivity, we compare 1 year of simulations from the Weather Research and Forecasting (WRF) model with and without wind plants incorporated, focusing on the lease area south of Massachusetts and Rhode Island. The simulations with wind plants are repeated to include both the maximum and minimum amounts of added turbulence to provide bounds on the potential impacts. We assess changes in wind speeds, 2 m temperature, surface heat flux, turbulence kinetic energy (TKE), and boundary layer height during different stability classifications and ambient wind speeds over the entire year and compare results for the degree of added turbulence in the wind plant simulations. Because the wake behavior may be a function of boundary layer stability, in this paper, we also present a machine learning algorithm to quantify the area and distance of the wake generated by the wind plant. This analysis enables us to identify the relationship between wake extent and boundary layer height. We find that hub-height wind speed is reduced within and downwind of the wind plant, with the strongest impacts occurring during stable conditions and faster wind speeds in region 3 of the turbine power curve, although impacts lessen as wind speeds increase past 15 m s−1. In contrast, wind speeds near the surface decrease when no turbine-added turbulence is included but can increase for stably stratified conditions when 100 % of possible TKE is included in the simulations. TKE increases at hub height in the simulations with added TKE for all stability classes, suggesting that atmospheric stability does not immediately modify the TKE generated by turbines. Negligible changes in hub-height TKE manifest in the simulations without the added TKE. At the surface, TKE increases in the simulations with maximum added turbulence only for unstable conditions. In the no-added-turbulence simulations, surface TKE decreases slightly in neutral and unstable simulations. Differences in 2 m temperatures and surface heat fluxes are small but vary considerably with atmospheric stability and the amount of added TKE. Boundary layer heights increase within the wind plant when turbine-added turbulence is included and decrease slightly downwind during stable conditions. In contrast, with no added turbulence, the boundary layer height is in general reduced in stable conditions with wind speeds less than 15 m s −1 and slightly increased in neutral conditions. Finally, shallower upwind boundary layer heights tend to correlate with larger wake areas and distances, though other factors likely also play a role in determining the extent of the wind plant wake. These simulation-based results provide a bound for micrometeorological impacts of wind plant wakes: simulations that couple the atmosphere to the ocean may reduce these impacts, and we await observational verification.

17 WIND ENERGY↗
NASA NTRS2000

The Surface Pressure Response of a NACA 0015 Airfoil Immersed in Grid Turbulence: Characteristics of the Turbulence - Volume 1

Two grids have been developed for the Virginia Tech 6 ft x 6 ft Stability wind tunnel for the purpose of generating homogeneous isotropic turbulent flows for the study of unsteady airfoil response. The first, a square bi-planar grid with a 12" mesh size and an open area ratio of 69.4%, was mounted in the wind tunnel contraction. The second grid, a metal weave with a 1.2 in. mesh size and an open area ratio of 68.2% was mounted in the tunnel test section. Detailed statistical and spectral measurements of the turbulence generated by the two grids are presented for wind tunnel free stream speeds of 10, 20, 30 and 40 m/s. These measurements show the flows to be closely homogeneous and isotropic. Both grids produce flows with a turbulence intensity of about 4% at the location planned for the airfoil leading edge. Turbulence produced by the large grid has an integral scale of some 3.2 inches here. Turbulence produced by the small grid is an order of magnitude smaller. For wavenumbers below the upper limit of the inertial subrange, the spectra and correlations measured with both grids at all speeds can be represented using the von Karman interpolation formula with a single velocity and length scale. The spectra maybe accurately represented over the entire wavenumber range by a modification of the von Karman interpolation formula that includes the effects of dissipation. These models are most accurate at the higher speeds (30 and 40 m/s).

Bereketab, Semere↗
DOE OSTI2023

Thermodynamic and kinematic drivers of atmospheric boundary layer stability in the central Arctic during the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC)

Observations collected during the Multidisciplinary drifting Observatory for the Study of Arctic Climate (MOSAiC) provide a detailed description of the impact of thermodynamic and kinematic forcings on atmospheric boundary layer (ABL) stability in the central Arctic. This study reveals that the Arctic ABL is stable and near-neutral with similar frequencies, and strong stability is the most persistent of all stability regimes. MOSAiC radiosonde observations, in conjunction with observations from additional measurement platforms, including a 10 m meteorological tower, ceilometer, microwave radiometer, and radiation station, provide insight into the relationships between atmospheric stability and various atmospheric thermodynamic and kinematic forcings of ABL turbulence and how these relationships differ by season. We found that stronger stability largely occurs in low-wind (i.e., wind speeds are slow), low-radiation (i.e., surface radiative fluxes are minimal) environments; a very shallow mixed ABL forms in low-wind, high-radiation environments; weak stability occurs in high-wind, moderate-radiation environments; and a near-neutral ABL forms in high-wind, high-radiation environments. Surface pressure (a proxy for synoptic staging) partially explains the observed wind speeds for different stability regimes. Cloud frequency and atmospheric moisture contribute to the observed surface radiation budget. Unique to summer, stronger stability may also form when moist air is advected from over the warmer open ocean to over the colder sea ice surface, which decouples the colder near-surface atmosphere from the advected layer, and is identifiable through observations of fog and atmospheric moisture.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Dynamic Subgrid Turbulence Modeling for Shallow Cumulus Convection Simulations beyond LES Resolutions

A scale-dependent dynamic Smagorinsky model is implemented in the Met Office/NERC Cloud (MONC) model using two averaging flavors, along Lagrangian pathlines and local moving averages. The dynamic approaches were compared against the conventional Smagorinsky–Lilly scheme in simulating the diurnal cycle of shallow cumulus convection. The simulations spanned from the LES to the near-gray-zone and gray-zone resolutions and revealed the adaptability of the dynamic model across the scales and different stability regimes. The dynamic model can produce a scale- and stability-dependent profile of the subfilter turbulence length scale across the chosen resolution range. At gray-zone resolutions the adaptive length scales can better represent the early precloud boundary layer leading to temperature and moisture profiles closer to the LES compared to the standard Smagorinsky. As a result, the initialization and general representation of the cloud field in the dynamic model is in good agreement with the LES. In contrast, the standard Smagorinsky produces a less well-mixed boundary layer, which fails to ventilate moisture from the boundary layer, resulting in the delayed spinup of the cloud layer. Moreover, strong downgradient diffusion controls the turbulent transport of scalars in the cloud layer. However, the dynamic approaches rely on the resolved field to account for nonlocal transports, leading to overenergetic structures when the boundary layer is fully developed and the Lagrangian model is used. Introducing the local averaging version of the model or adopting a new Lagrangian time scale provides stronger dissipation without significantly affecting model behavior.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Organizing bacterial vortex lattices by periodic obstacle arrays

Recent experiments have shown that the complex spatio-temporal vortex structures emerging in active fluids are susceptible to weak geometrical constraints. This observation poses the fundamental question of how boundary effects stabilize a highly ordered pattern from seemingly turbulent motion. Here we show, by a combination of continuum theory and experiments on a bacterial suspension, how artificial obstacles guide the flow profile and reorganize topological defects, which enables the design of bacterial vortex lattices with tunable properties. To this end, the continuum model is extended by appropriate boundary conditions. Beyond the stabilization of square and hexagonal lattices, we also provide a striking example of a chiral, antiferromagnetic lattice exhibiting a net rotational flow, which is induced by arranging the obstacles in a Kagome-like array. Recent experiments showed that weak geometrical constraints can organize topological defects in turbulent bacterial suspensions. Here, the authors use a continuum model to study the connection between symmetry and stability of emergent vortex patterns and the geometry of constraining pillar arrays.

59 BASIC BIOLOGICAL SCIENCES↗