Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “TKE”

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 19 records

Preprocessing for Eddy Dissipation Rate and TKE Profile Generation

The Aircraft Vortex Spacing System (AVOSS), a set of algorithms to determine aircraft spacing according to wake vortex behavior prediction, requires turbulence profiles to appropriately determine arrival and departure aircraft spacing. The ambient atmospheric turbulence profile must always be produced, even if the result is an arbitrary (canned) profile. The original turbulence profile code was generated By North Carolina State University and used in a non-real-time environment in the past. All the input parameters could be carefully selected and screened prior to input. Since this code must run in real-time using actual measurements in the field as input, it became imperative to begin a data checking and screening process as part of the real-time implementation. The process described herein is a step towards ensuring that the best possible turbulence profile is always provided to AVOSS. Data fill-ins, constant profiles and arbitrary profiles are used only as a last resort, but are essential to ensure uninterrupted application of AVOSS.

Zak, J. Allen↗

Interannual Variability of Seasonal Pblooms in the North Polar Atlantic in Response to Atmospheric Forcing

The year-to-year variability of timing, intensity, and spatial distribution of surface phytoplankton during spring-summer seasons is examined in the north polar region of the Atlantic using satellite-derived chlorophyll a concentration (Chl) over 7 years (1998 2004). Each year phytoplankton bloom differed in onset, temporal evolution, and intensity. This interannual variability is to a large degree controlled by local weather, as indicated by the analysis of coincident meteorological data. The atmospheric forcing was parameterized in terms of the generation rate of turbulent kinetic energy (TKE(sub RT)) supplied from the atmosphere to the ocean. Timing of the bloom is delayed in years with high TKE(sub RT) supplied to the ocean in March. In April, Chl (local and regionally averaged) and TKE(sub RT) are inversely related to one another. The late winter and early spring atmospheric conditions influence the seasonal Chl patterns. The seasonal (April-August) Chl correlates well with net heat flux, wind energy, and TKE(sub RT) in March and April, but the correlation can be negative or positive in different areas of the north polar Atlantic. The correlation between TKE(sub RT) and seasonal Chl is positive in the Greenland Gyre (higher seasonal Chl corresponds to higher TKE(sub RT) in March April) and negative in the regions of East Greenland and North Atlantic/West Spitsbergen Currents. The most likely explanation for the positive correlation between TKE(sub RT) and Chl within the Greenland Gyre is that higher TKE(sub RT) in spring increases the seasonal supply of nutrients into surface waters.

Stramska, Malgorzata↗

Flow Field and Acoustic Predictions for Three-Stream Jets

Computational fluid dynamics was used to analyze a three-stream nozzle parametric design space. The study varied bypass-to-core area ratio, tertiary-to-core area ratio and jet operating conditions. The flowfield solutions from the Reynolds-Averaged Navier-Stokes (RANS) code Overflow 2.2e were used to pre-screen experimental models for a future test in the Aero-Acoustic Propulsion Laboratory (AAPL) at the NASA Glenn Research Center (GRC). Flowfield solutions were considered in conjunction with the jet-noise-prediction code JeNo to screen the design concepts. A two-stream versus three-stream computation based on equal mass flow rates showed a reduction in peak turbulent kinetic energy (TKE) for the three-stream jet relative to that for the two-stream jet which resulted in reduced acoustic emission. Additional three-stream solutions were analyzed for salient flowfield features expected to impact farfield noise. As tertiary power settings were increased there was a corresponding near nozzle increase in shear rate that resulted in an increase in high frequency noise and a reduction in peak TKE. As tertiary-to-core area ratio was increased the tertiary potential core elongated and the peak TKE was reduced. The most noticeable change occurred as secondary-to-core area ratio was increased thickening the secondary potential core, elongating the primary potential core and reducing peak TKE. As forward flight Mach number was increased the jet plume region decreased and reduced peak TKE.

RANS CFD↗

Problems with the process partitioning theory of stratocumulus entrainment

Three different approaches to partitioning were proposed. Ball (1960), Lilly (1968), and Deardorff et al. (1969, 1974) considered the sign of the net buoyancy flux at each level. If the net flux is positive, it is counted as TKE producing; otherwise, it is counted as TKE consuming. This approach can be called Eulerian partitioning. The second approach can be called process partitioning. It is assumed that the various processes acting in concert each produce and consume the same energy as if they acted independently (Manins and Turner, 1978). The total rates of TKE production and consumption are obtained by summing the effects of all the forcing processes. The third approach is Lagrangian partitioning. Each air parcel is considered as either producing or consuming TKE, according to the sign of the product of its density and vertical velocity anomalies. Stage and Businger (1981a,b) have applied process partitioning to the cloud topped mixed layer. One of the most important processes influencing entrainment into such a layer is cloud top radiative cooling. The production and consumption due to entrainment and radiative cooling are thus closely related. Model results are sensitive to the choice of formulation.

Randall, D. A.↗

Mesoscale Simulations of a Florida Sea Breeze Using the PLACE Land Surface Model Coupled to a 1.5-Order Turbulence Parameterization

A sophisticated land-surface model, PLACE, the Parameterization for Land Atmospheric Convective Exchange, has been coupled to a 1.5-order turbulent kinetic energy (TKE) turbulence sub-model. Both have been incorporated into the Penn State/National Center for Atmospheric Research (PSU/NCAR) mesoscale model MM5. Such model improvements should have their greatest effect in conditions where surface contrasts dominate over dynamic processes, such as the simulation of warm-season, convective events. A validation study used the newly coupled model, MM5 TKE-PLACE, to simulate the evolution of Florida sea-breeze moist convection during the Convection and Precipitation Electrification Experiment (CaPE). Overall, eight simulations tested the sensitivity of the MM5 model to combinations of the new and default model physics, and initialization of soil moisture and temperature. The TKE-PLACE model produced more realistic surface sensible heat flux, lower biases for surface variables, more realistic rainfall, and cloud cover than the default model. Of the 8 simulations with different factors (i.e., model physics or initialization), TKE-PLACE compared very well when each simulation was ranked in terms of biases of the surface variables and rainfall, and percent and root mean square of cloud cover. A factor separation analysis showed that a successful simulation required the inclusion of a multi-layered, land surface soil vegetation model, realistic initial soil moisture, and higher order closure of the planetary boundary layer (PBL). These were needed to realistically model the effect of individual, joint, and synergistic contributions from the land surface and PBL on the CAPE sea-breeze, Lake Okeechobee lake breeze, and moist convection.

Lynn, Barry H.↗

Budgets of Reynolds Stresses and Turbulent Heat Flux for Hypersonic Turbulent Boundary Layers Subject to Pressure Gradients

In this paper, the budget terms in the transport equations of the Reynolds stresses, the turbulent kinetic energy (TKE), the internal energy, and the turbulent heat flux are computed using a direct numerical simulation (DNS) database of hypersonic turbulent boundary layers (TBLs) subject to favorable pressure gradients (FPG). Study of the Reynolds stress and TKE budgets shows that FPG has a damping effect on turbulence production, transport, and destruction. The semilocal scaling fails to collapse the Reynolds stress and TKE budget terms between the zero pressure gradient (ZPG) and the strong FPG cases, suggesting that the mechanical nonequilibrium due to FPG leads to changes in Reynolds-stress and TKE transport that cannot be captured by Morkovin’s hypothesis. Minimal change is seen in the internal energy budget in the presence of FPG, indicating a minimal impact due to FPG on the mean thermal field. While the near-wall transport of wall-normal turbulent heat flux is largely governed by the balance between the turbulent viscous-thermal transport term and the turbulent viscous-thermal dissipation term, the transport of the streamwise turbulent heat flux is influenced by contributions from almost all of the budget terms. For the strong FPG case in particular, production due to both the mean and the fluctuating strain rates begins to play a significant role in the near-wall transport of the streamwise turbulent heat flux, and their inclusion or modeling may be required to accurately predict turbulent heat fluxes for a strong FPG.

Reynolds stress↗

An integral turbulent kinetic energy analysis of free shear flows

Mixing of coaxial streams is analyzed by application of integral techniques. An integrated turbulent kinetic energy (TKE) equation is solved simultaneously with the integral equations for the mean flow. Normalized TKE profile shapes are obtained from incompressible jet and shear layer experiments and are assumed to be applicable to all free turbulent flows. The shear stress at the midpoint of the mixing zone is assumed to be directly proportional to the local TKE, and dissipation is treated with a generalization of the model developed for isotropic turbulence. Although the analysis was developed for ducted flows, constant-pressure flows were approximated with the duct much larger than the jet. The axisymmetric flows under consideration were predicted with reasonable accuracy. Fairly good results were also obtained for the fully developed two-dimensional shear layers, which were computed as thin layers at the boundary of a large circular jet.

Peters, C. E.↗

Turbulence modeling in hypersonic inlets

A study is conducted to analyze the performance of different turbulence models when applied to flow through a Mach 7.4 hypersonic inlet. The analysis, which is two-dimensional, is done by comparing computational results from a Parabolized Navier-Stokes code and a full Navier-Stokes code, with experimental data. The McDonald-Camarata (MC) and Baldwin-Lomax (BL) models were the two zero-equation models used in the study. The Turbulent Kinetic Energy (TKE) model was chosen as a representative higher order model. The MC model, when run with transition of flow, provides a solution which compares excellently with the data. Transition has a first order effect on the overall solution provided by the code. The BL model predicts separation of flow in the inlet, which contradicts experimental findings. The TKE model does not perform any better than the MC and BL models, despite the fact that it is a higher order turbulence model. The BL and TKE models predict transition in the inlet at a location which is much earlier than observed in the experiment. This may be attributed to the empirical constants used to determine the point of transition.

Ng, W. F.↗

Effects of shock strength on shock turbulence interaction

Direct numerical simulation (DNS) and linear analysis (LIA) of isotropic turbulence interacting with a shock wave are performed for several upstream shock normal Mach numbers (M(sub 1)). Turbulence kinetic energy (TKE) is amplified across the shock wave, but this amplification tends to saturate beyond M(sub 1) = 3.0. TKE amplification and Reynolds stress anisotropy obtained in DNS are consistent with LIA predictions. Rapid evolution of TKE immediate downstream of the shock wave persists for all shock strengths and is attributed to the transfer between kinetic and potential modes of turbulence energy through acoustic fluctuations. Changes in energy spectra and various length scales across the shock wave are predicted by LIA, which is consistent with DNS results. Most turbulence length scales decrease across the shock. Dissipation length scale (rho-bar q(exp 3) / epsilon), however, increases slightly for shock waves with M(sub 1) less than 1.65. Fluctuations in thermodynamic variables behind the shock wave stay nearly isentropic for M(sub 1) less than 1.2 and deviate significantly from isentropy for the stronger shock waves due to large entropy fluctuation generated through the interaction.

Lee, Sangsan↗

Full-coverage film cooling. I - Three-dimensional measurements of turbulence structure. II - Prediction of the recovery-region hydrodynamics

Hydrodynamic measurements of turbulence structure were performed with a triaxial hot wire in the full coverage and the recovery regions following an array of injection holes under isothermal conditions at ambient temperature and pressure for blowing ratios of 0.9 and 0.4. High levels of turbulence kinetic energy (TKE) were determined for low blowing, and low TKE levels were found for the high blowing levels; in the recovery region, the flow can be represented by a model with an outer boundary layer and a 2-dimensional inner boundary layer. Recovery region hydrodynamics can be modelled by considering that a new boundary layer started to grow immediately after the end of blowing; the Prandtl mixing length distributions calculated from the values of mean velocity and turbulent shear stresses were consistent with the presence of a dual boundary layer structure in the recovery region. The program used here contains a one-equation model of turbulence, using turbulence kinetic energy with an algebraic mixing length; this 2-dimensional, finite difference program can predict the mean velocity and turbulence kinetic energy profiles based on initial values, boundary conditions, and a closure condition.

Yavuzkurt, S.↗

Buoyant production and consumption of turbulence kinetic energy in cloud-topped mixed layers

It is pointed out that studies of the entraining planetary boundary layer (PBL) have generally emphasized the role of buoyancy fluxes in driving entrainment. The buoyancy flux is proportional to the rate of conversion of the potential energy of the mean flow into the kinetic energy of the turbulence. It is not unusual for conversion to proceed in both directions simultaneously. This occurs, for instance, in both clear and cloudy convective mixed layers which are capped by inversions. A partitioning of the net conversion into positive parts, generating turbulence kinetic energy (TKE), and negative parts (TKE-consuming), would make it possible to include the positive part in the gross production rate, and closure would be achieved. Three different approaches to partitioning have been proposed. The present investigation is concerned with a comparison of the three partitioning theories. Particular attention is given to the cloud-topped mixed layer because in this case the differences between two partitioning approaches are most apparent.

Randall, D. A.↗

An Estimation of Turbulent Kinetic Energy and Energy Dissipation Rate Based on Atmospheric Boundary Layer Similarity Theory

Algorithms are developed to extract atmospheric boundary layer profiles for turbulence kinetic energy (TKE) and energy dissipation rate (EDR), with data from a meteorological tower as input. The profiles are based on similarity theory and scalings for the atmospheric boundary layer. The calculated profiles of EDR and TKE are required to match the observed values at 5 and 40 m. The algorithms are coded for operational use and yield plausible profiles over the diurnal variation of the atmospheric boundary layer.

Han, Jongil↗

Measurement of the Turbulence Kinetic Energy Budget of a Turbulent Planar Wake Flow in Pressure Gradients

Turbulence kinetic energy (TKE) is a very important quantity for turbulence modeling and the budget of this quantity in its transport equation can provide insight into the flow physics. Turbulence kinetic energy budget measurements were conducted for a symmetric turbulent wake flow subjected to constant zero, favorable and adverse pressure gradients in year-three of research effort. The purpose of this study is to clarify the flow physics issues underlying the demonstrated influence of pressure gradient on wake development and provide experimental support for turbulence modeling. To ensure the reliability of these notoriously difficult measurements, the experimental procedure was carefully designed on the basis of an uncertainty analysis. Four different approaches, based on an isotropic turbulence assumption, a locally axisymmetric homogeneous turbulence assumption, a semi-isotropy assumption and a forced balance of the TKE equation, were applied for the estimate of the dissipation term. The pressure transport term is obtained from a forced balance of the turbulence kinetic energy equation. This report will present the results of the turbulence kinetic energy budget measurement and discuss their implication on the development of strained turbulent wakes.

Liu, Xiao-Feng↗

Numerical Study of the Role of Shallow Convection in Moisture Transport and Climate

The objective of this investigation was to study the role of shallow convection on the regional water cycle of the Mississippi and Little Washita Basins of the Southern Great Plains (SGP) using a 3-D mesoscale model, the PSU/NCAR MM5. The underlying premise of the project was that current modeling of regional-scale climate and moisture cycles over the continents is deficient without adequate treatment of shallow convection. At the beginning of the study, it was hypothesized that an improved treatment of the regional water cycle can be achieved by using a 3-D mesoscale numerical model having high-quality parameterizations for the key physical processes controlling the water cycle. These included a detailed land-surface parameterization (the Parameterization for Land-Atmosphere-Cloud Exchange (PLACE) sub-model of Wetzel and Boone), an advanced boundary-layer parameterization (the 1.5-order turbulent kinetic energy (TKE) predictive scheme of Shafran et al.), and a more complete shallow convection parameterization (the hybrid-closure scheme of Deng et al.) than are available in most current models. PLACE is a product of researchers working at NASA's Goddard Space Flight Center in Greenbelt, MD. The TKE and shallow-convection schemes are the result of model development at Penn State. The long-range goal is to develop an integrated suite of physical sub-models that can be used for regional and perhaps global climate studies of the water budget. Therefore, the work plan focused on integrating, improving, and testing these parameterizations in the MM5 and applying them to study water-cycle processes over the SGP. These schemes have been tested extensively through the course of this study and the latter two have been improved significantly as a consequence.

Seaman, Nelson L.↗

Wake Vortex Algorithm Scoring Results

This report compares the performance of two models of trailing vortex evolution for which interaction with the ground is not a significant factor. One model uses eddy dissipation rate (EDR) and the other uses the kinetic energy of turbulence fluctuations (TKE) to represent the effect of turbulence. In other respects, the models are nearly identical. The models are evaluated by comparing their predictions of circulation decay, vertical descent, and lateral transport to observations for over four hundred cases from Memphis and Dallas/Fort Worth International Airports. These observations were obtained during deployments in support of NASA's Aircraft Vortex Spacing System (AVOSS). The results of the comparisons show that the EDR model usually performs slightly better than the TKE model.

Robins, R. E.↗

Interaction of a Boundary Layer with a Turbulent Wake

The objective of this grant was to study the transition mechanisms on a flat-plate boundary layer interacting with the wake of a bluff body. This is a simplified configuration presented and designed to exemplify the phenomena that occur in multi-element airfoils, in which the wake of an upstream element impinges on a downstream one. Some experimental data is available for this configuration at various Reynolds numbers. The first task carried out was the implementation and validation of the immersed-boundary method. This was achieved by performing calculations of the flow over a cylinder at low and moderate Reynolds numbers. The low-Reynolds number results are discussed, which is enclosed as Appendix A. The high-Reynolds number results are presented in a paper in preparation for the Journal of Fluid Mechanics. We performed calculations of the wake-boundary-layer interaction at two Reynolds numbers, Re approximately equal to 385 and 1155. The first case is discussed and a comparison of the two calculations is reported. The simulations indicate that at the lower Reynolds number the boundary layer is buffeted by the unsteady Karman vortex street shed by the cylinder. This is shown: long streaky structures appear in the boundary layer in correspondence of the three-dimensionalities in the rollers. The fluctuations, however, cannot be self-sustained due to the low Reynolds-number, and the flow does not reach a turbulent state within the computational domain. In contrast, in the higher Reynolds-number case, boundary-layer fluctuations persist after the wake has decayed (due, in part, to the higher values of the local Reynolds number Re achieved in this case); some evidence could be observed that a self-sustaining turbulence generation cycle was beginning to be established. A third simulation was subsequently carried out at a higher Reynolds number, Re=3900. This calculation gave results similar to those of the Re=l155 case. Turbulence was established at fairly low Reynolds number, as a consequence of the high level of the free-stream perturbation. An instantaneous flow visualization for that case is shown. A detailed examination of flow statistics in the transitional and turbulent regions, including the evolution of the turbulent kinetic energy (TKE) budget and frequency spectra showed the formation and evolution of turbulent spots characteristic of the bypass transition mechanism. It was also observed that the turbulent eddies achieved an equilibrium, fully developed turbulent states first, as evidenced by the early agreement achieved by the terms in the TKE budget with those observed in turbulent flows. Once a turbulent Reynolds stress profile had been established, the velocity profile began to resemble a turbulent one, first in the inner region and later in the outer region of the wall layer. An extensive comparison of the three cases, including budgets, mean velocity and Reynolds stress profiles and flow visualization, is included. The results obtained are also presented.

Piomelli, Ugo↗

Influence of Satellite-Based Heterogeneous Vegetation Momentum Roughness on Mesoscale Model Dynamics During IHOP 2002

The sensitivity of mesoscale weather prediction model to a vegetation roughness initialization is investigated for the south central United States. Three different roughness databases are employed: i) a control or standard lookup table roughness that is a function only of land cover type, ii) a spatially heterogeneous roughness database previously derived using a physically based procedure and MODIS imagery, and iii) a MODIS climatologic roughness database that possesses the same spatial heterogeneity as (i) but with mean land class values from (ii). The model used is the Weather Research and Forecast Model (WRF) coupled to the Community Land Model within the Land Information System (LIS). For each simulation, a statistical comparison is made between modeled results and ground observations from meteorological stations within the Oklahoma mesonet and surrounding region during IHOP20O2. A sensitivity analysis on the impact the MODIS-based roughness fields is also made through a time-series intercomparison of temperature bias, probability of detection (POD), average wind speed, boundary layer height, and turbulent kinetic energy (TKE) the results that, for the current replacement of the standard land-cover type based roughness values with the satellite-derived fields statistically improves model performance for most of the observed variables. Further, the satellite-based roughness enhances the surface wind speed, PBL height and TKE production on the order of 3 to l0 percent, with a lesser effect over grassland and cropland domains, and the greater effect over mixed land cover domains

Jasinski, Michael↗

Numerical Model Sensitivity to Heterogeneous Satellite Derived Vegetation Roughness

The sensitivity of a mesoscale weather prediction model to a 1 km satellite-based vegetation roughness initialization is investigated for a domain within the south central United States. Three different roughness databases are employed: i) a control or standard lookup table roughness that is a function only of land cover type, ii) a spatially heterogeneous roughness database, specific to the domain, that was previously derived using a physically based procedure and Moderate Resolution Imaging Spectroradiometer (MODIS) imagery, and iii) a MODIS climatologic roughness database that like (i) is a function only of land cover type, but possesses domain specific mean values from (ii). The model used is the Weather Research and Forecast Model (WRF) coupled to the Community Land Model within the Land Information System (LIS). For each simulation, a statistical comparison is made between modeled results and ground observations within a domain including Oklahoma, Eastern Arkansas, and Northwest Louisiana during a 4-day period within IHOP 2002. Sensitivity analysis compares the impact the three roughness initializations on time-series temperature, precipitation probability of detection (POD), average wind speed, boundary layer height, and turbulent kinetic energy (TKE). Overall, the results indicate that, for the current investigation, replacement of the standard look-up table values with the satellite-derived values statistically improves model performance for most observed variables. Such natural roughness heterogeneity enhances the surface wind speed, PBL height and TKE production up to 10 percent, with a lesser effect over grassland, and greater effect over mixed land cover domains.

Jasinski, Michael↗