Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “turbulent length scale”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS1984

Effect of an oscillating flow direction on leading edge heat transfer

An experimental investigation was conducted to examine the effect of a periodic variation in the angle of attack on heat transfer at the leading edge of a gas turbine blade. A circular cylinder was used as a large-scale model of the leading edge region. The cylinder was placed in a wind tunnel and was oscillated rotationally about its axis. The incident flow Reynolds number and the Strouhal number of oscillation were chosen to model an actual turbine condition. Incident turbulence levels up to 4.9 percent were produced by grids placed upstream of the cylinder. The transfer rate was measured using a mass transfer technique and heat transfer rates inferred from the results. A direct comparison of the unsteady and steady results indicate that the effect is dependent on the Strouhal number, turbulence level, and the turbulence length scale, but that the largest observed effect was only a 10 percent augmentation at the nominal stagnation position.

Marziale, M. L.↗
NASA NTRS1976

Experimental and theoretical studies of subsonic fan noise

The noise generated by inlet turbulence impinging on a subsonic axial flow fan was studied as a function of tip speed, flow coefficient, and intensity and scale of turbulence was carried out. Both turbulence and far field acoustic measurements were made. The new elements introduced in the theoretical analysis were accounting for blade loading dependent noise mechanisms and consideration of anisotropic turbulence impinging on the rotor because of inlet flow contraction effects. Experimentally an unexplained increase of noise at about 1/2 and 1 1/2 times blade passsing frequency was observed at low flow coefficients even though there was no evidence of compressor surge. In the final version the theory does a fair job of predicting variations of noise with blade loading and tip speed. Alteration of inlet turbulence length scales produced some but not very pronounced changes in the far field PWL spectra. Some degree of eddy contraction and resulting anisotropy were essential to explain the concentration of energy around blade passing frequencies.

Mani, R.↗
DOE OSTI2024

DRDMannTurb: A Python package for scalable, data-driven synthetic turbulence

Synthetic turbulence models (STMs) are used in wind engineering to generate realistic flow fields and are employed as inputs to industrial wind simulations. Examples include prescribing inlet conditions in large eddy simulations that model loads on wind turbines and tall buildings. We are interested in STMs capable of generating fluctuations based on prescribed second-moment statistics since such models can simulate environmental conditions that closely resemble on-site observations. To this end, the widely used Mann model (see Mann, 1994, 1998) is the inspiration for DRDMannTurb. The Mann model is described by three physical parameters: a magnitude parameter influencing the global variance of the wind field and corresponding to the Kolmogorov constant multiplied by the rate of viscous dissipation of the turbulent kinetic energy to the two-thirds, αϵ 2/3 , a turbulence length scale parameter L, and a nondimensional parameter Γ related to the lifetime of the eddies. A number of studies, as well as international standards (e.g., those by the International Electrotechnical Commission (IEC)), include recommended values for these three parameters with the goal of standardizing wind simulations according to observed energy spectra. Yet, having only three parameters, the Mann model faces limitations in accurately representing the diversity of observable spectra. This Python package enables users to extend the Mann model and more accurately fit field measurements through flexible neural network models of the eddy lifetime function. Following Keith et al. (2021), we refer to this class of models as Deep Rapid Distortion (DRD) models. DRDMannTurb also includes a general module implementing an efficient method for synthetic turbulence generation based on a domain decomposition technique. This technique is also described in Keith et al. (2021).

17 WIND ENERGY↗
DOE OSTI2023

Wind Loading on Parabolic Trough Collectors: Wind and Structural Loads Measurements at an Operational Powerplant

Concentrated Solar Power (CSP) is a promising method for using solar power for electricity generation with thermal energy storage for industrial applications. Solar collectors constitute almost 1/3 of the total cost of the power plant. One of the primary drivers of unrealiability of these collectors is wind-driven loading of mirrors, support structures, and drives. To date, the design of the solar collector structures has relied on data from wind tunnels that do not adequately capture the dynamic effects observed at scale. NREL initiated a field measurement campaign at the operational Nevada Solar One (NSO) parabolic trough powerplant. At this plant, parabolic trough solar collectors track the sun from east to west in the course of a day and face varying wind loads depending on the wind properties and the angle of the troughs. The aim of the project is a detailed characterization of prevailing wind and turbulence conditions and resulting operational loads on parabolic troughs, providing insights into structural dynamic response, and generating a comprehensive wind-loading dataset for validating simulations of wind loading on collector structures. The measurements at NSO consist of sonic anemometers on masts at different heights to characterize the incoming flow and conditions at four trough rows at the edge of the trough field. In addition, a Doppler Lidar scans the horizontal plane above the troughs. The wind measurements at NSO have been continuously collecting data since October 2021 and are combined with structural load measurements that started in November 2022. The load measurements were installed on the same four outermost trough rows and include support structure bending moments, drive torque moments, dynamic accelerations of the spaceframe, mirror displacement, and tilt angles. Measurements are planned to continue until May 2023, providing a first-of-a-kind, high-resolution multi-month dataset of combined wind and load measurements. In this presentation, we show first results of the measurement campaign. Based on the measurements, we identified three main factors altering the flow over the parabolic troughs: Wind speed, wind direction, and the angle of the parabolic troughs. Most interactions between incoming wind and the trough field are observed when the wind blows perpendicular to the trough rows. In this case, the first rows experience the highest wind speed and block the downwind rows, creating conditions with decreased wind speed and enhanced turbulence within the trough field. Also, turbulent length scales are smaller after the first row. This leads to the highest static loads (bending and torque moments) at the first row but potentially increased dynamic loads within the field. We illustrate our findings with case studies focused on different wind conditions and will present fatigue analysis to highlight the impact of wind-driven loads on collector structures.

CSP collectors↗
DOE OSTI2024

The influence of the shock-to-reshock time on the Richtmyer–Meshkov instability in reshock

Experiments on the Richtmyer–Meshkov instability (RMI) in a dual driver vertical shock tube (DDVST) are described. An initially planar, stably stratified membraneless interface is formed by flowing air from above and sulfur hexafluoride from below the interface location using the method of Jones & Jacobs ( Phys. Fluids , vol. 9, issue 1997, 1997, pp. 3078–3085). A random three-dimensional, multi-modal initial perturbation is imposed by vertically oscillating the gas column to produce Faraday waves. The DDVST design generates two shock waves, one originating above and one below the interface, with these shocks having independently controllable strengths and interface arrival times. The shock waves have nominal strengths of $M_L=1.17$ and $M_H=1.18$ for the shock wave originating in the light and heavy gas, respectively, with these strengths chosen to result in arrested bulk interface motion following reshock. The influence of the length of the shock-to-reshock time, as well as the order of shock arrival, on the post-reshock RMI is examined. The mixing layer width grows according to $h\propto t^\theta$ , where $\theta _H=0.36\pm 0.018$ (95 %) and $\theta _L=0.38\pm 0.02$ (95 %) for heavy and light shock first experiments, respectively, indicating no strong dependence on the order of shock wave arrival. Volume integrated specific turbulent kinetic energy (TKE) in the mixing layer versus time is found to decay according to $E_{tot}/\bar {\rho }\propto t^p$ with $p_H=-0.823\pm 0.06$ (95 %) and $p_L=-1.061\pm 0.032$ (95 %) for heavy and light shock first experiments, respectively. Notably, the 95 % confidence intervals do not overlap. Analysis on the influence of the shock-to-reshock time on turbulent length scales, transition criteria, spectra and mixing layer anisotropy are also presented.

Ferguson, Kevin (ORCID:0000000166537482)↗
DOE OSTI2024

A method for examining ensemble averaging forms during the transition to turbulence in HED systems for application to RANS models

This paper discusses a strategy to initialize a two-dimensional (2D) Reynolds-averaged Navier–Stokes model [LANL's Besnard–Harlow–Rauenzahn (BHR) model] in order to describe an unsteady transitional Richtmyer–Meshkov (RM)-induced flow observed in on-going high-energy-density ensemble experiments performed on the OMEGA-EP facility. The experiments consist of a nominal single-mode perturbation (initial amplitude a 0 ≈ 10 and wavelength $λ$ = 100μm) with target-to-target variations in the surface roughness subjected to the RM instability with delayed Rayleigh–Taylor in a heavy-to-light configuration. Our strategy leverages high-resolution three-dimensional (3D) implicit large eddy simulations (ILES) simulations to initialize BHR-relevant parameters and subsequently validate the 2D BHR results against the 3D ILES simulations. A suite of five 3D ILES simulations corresponding to five experimental target profiles is undertaken to generate an ensemble dataset. Using ensemble averages from the 3D simulations to initialize the turbulent kinetic energy in the BHR model ( K 0 ) demonstrates the ability of the model to predict the time evolution of the interface as well as the density-specific-volume covariance, b . To quantify the sensitivity of the BHR results to the choice of K 0 and the initial turbulent length scale, S 0 , we execute a parameter sweep spanning four orders of magnitude for both S 0 and K 0 , generating a parameter space consisting of 26 simulations. The Pearson's correlation coefficient is used as a measure of discrepancy between the 2D BHR and 3D ILES simulations and reveals that the ranges 8≲S 0 ≲20 μm and 10 9 ≲K 0 ≲10 10 cm 2 /s 2 produce predictions that agree best with the 3D ILES results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Direct Numerical Simulation of the Flow through a Randomly Packed Pebble Bed

The proposition for molten salt and high-temperature gas-cooled reactors has increased the focus on the dynamics and physics in randomly packed pebble beds. Research is being conducted on the validity of these designs as a possible contestant for the fourth-generation nuclear power systems. A detailed understanding of the coolant flow behavior is required in order to ensure proper cooling of the reactor core during normal and accident conditions. In order to increase the understanding of the flow through these complex geometries and enhance the accuracy of lower-fidelity modeling, high-fidelity approaches such as direct numerical simulation (DNS) can be utilized. Nek5000, a spectral-element computational fluid dynamics (CFD) code, was used to develop DNS fluid flow data. The flow domain consisted of 147 pebbles enclosed by a bounding wall. In the work presented, the Reynolds numbers ranged from 430 to 1050 based on the pebble diameter and inlet velocity. Characteristics of the flow domain such as volume averaged porosity, axial porosity, and radial porosity were studied and compared with correlations available in the literature. Friction factors from the DNS results for all Reynolds numbers were compared with correlations in the literature. The first- and second-order statistics show good agreement with the available experimental data. Turbulence length scales were analyzed in the flow. Reynolds stress anisotropy was characterized by utilizing invariant analysis. Overall, the results of the analysis in this study provide deeper understanding of the flow behavior and the effect of the wall in packed beds.

Yildiz, Mustafa Alper↗
NASA NTRS1974

Determination of critical nondimensional parameters in aircraft dynamic response to random input

The critical parameters of subsonic jet aircraft response in a random atmospheric environment are determined. Equations of motion are presented for semirigid aircraft with a flexible primary airfoil. However, the analysis is easily extendable to include additional appendage flexibility. The analysis establishes the mechanical admittance values for pitching, plunging, and the first mode effects from wing elastic bending and torsion. Nondimensional parameters are established which allow the representation of all subsonic jet transport aircraft with one nondimensional model. The critical parameters for random forcing are found to be aircraft relative mass, reduced natural and forcing frequencies, and Mach number. Turbulence scale lengths are found to be directly related to the critical values of reduced forcing frequency. Results are given for subsonic craft traveling at constant altitude. Specific values of admittance functions are tabulated at Mach numbers of 0.2, 0.5, and 0.7. The relative mass range covers all aircraft currently in operation.

Hillard, S. E.↗
NASA NTRS1978

Analysis of turbulent free jet hydrogen-air diffusion flames with finite chemical reaction rates

The nonequilibrium flow field resulting from the turbulent mixing and combustion of a supersonic axisymmetric hydrogen jet in a supersonic parallel coflowing air stream is analyzed. Effective turbulent transport properties are determined using the (K-epsilon) model. The finite-rate chemistry model considers eight reactions between six chemical species, H, O, H2O, OH, O2, and H2. The governing set of nonlinear partial differential equations is solved by an implicit finite-difference procedure. Radial distributions are obtained at two downstream locations of variables such as turbulent kinetic energy, turbulent dissipation rate, turbulent scale length, and viscosity. The results show that these variables attain peak values at the axis of symmetry. Computed distributions of velocity, temperature, and mass fraction are also given. A direct analytical approach to account for the effect of species concentration fluctuations on the mean production rate of species (the phenomenon of unmixedness) is also presented. However, the use of the method does not seem justified in view of the excessive computer time required to solve the resulting system of equations.

Sislian, J. P.↗
NASA NTRS1980

Flow measurements in a model swirl combustor

Flow measurements in a model combustor composed of two confined coaxial swirling jets under noncombusting conditions are presented. Mean flow results are obtained for five different flow conditions to determine the effect of outer swirl on the recirculation zone (which is used for flame stabilization under combustion conditions). As the outer swirl magnitude is first decreased from maximum counter-swirl to zero and then increased to give co-swirl conditions, the size and the reverse flow velocity in the recirculation zone diminish. Detailed time mean and fluctuating flow measurements are obtained for a co-swirl and a counter-swirl condition with a directional pitot probe and hot-wire anenometry. For these two cases, recirculation zone occurs only with counter-swirl, near the exit of the inner jet. The recirculation zone is in the form of a one celled toroidal vortex having very low swirl velocity. Axial development of fluctuation levels, energy dissipation rates and turbulence length scales are described for the two flow conditions. Spectral analysis reveals periodic oscillations in both flows. The oscillations originate from the inner jet and up to 4 harmonics are observed. The fundamental frequencies are comparable to the rotational frequencies of the inner jets at the exit which are approximately under solid body rotation. The significance of the present results for the combustion process is discussed.

Vu, B. T.↗
NASA NTRS1980

On the calculation of turbulent heat transport downstream from an abrupt pipe expansion

A numerical study is reported of flow and heat transfer in the separated flow region created by an abrupt pipe expansion. Computations employed an adaptation of the TEACH-2E computer program with the standard model of turbulence. Emphasis is given to the simulation, from both a physical and numerical viewpoint, of the region in the immediate vicinity of the wall where turbulent transport gives way to molecular conduction and diffusion. Wall resistance laws or wall functions used to bridge this near-wall region are based on the idea that, beyond the viscous sublayer, the turbulent length scale is universal, increasing linearly with distance from the wall. Predictions of expermental data for a diameter ratio of 0.54 show generally encouraging agreement with experiment. At a diameter of 0.43 different trends are discernible between measurement and calculation though this appears to be due to effects unconnected with the wall region studied.

Chieng, C. C.↗
NASA NTRS1983

Determination of the properties of magnetic turbulence in radio sources

The transport of polarized synchrotron radiation in a source possessing a highly irregular magnetic field, as proposed by Laing has been considered. The transport equation has been solved in a special case, relating the observable correlation functions in the Stokes parameters Q and U to the correlation functions of magnetic field and plasma density in the source. A rough application of the results to observations of the radio galaxy 3C 166 indicates that the turbulent scale length may be a few percent of the lobe size.

Spangler, S. R.↗
NASA NTRS1983

Turbulence and surface heat transfer near the stagnation point of a circular cylinder in turbulent flow

A turbulent boundary layer flow analysis of the momentum and thermal flow fields near the forward stagnation point due to a circular cylinder in turbulent cross flow is presented. Turbulence modeling length scale, anisotropic turbulence initial profiles and boundary conditions were identified as functions of the cross flow turbulence intensity and the boundary layer flow far field velocity. These parameters were used in a numerical computational procedure to calculate the mean velocity, mean temperature, and turbulence double correlation profiles within the flow field. The effects of the cross flow turbulence on the stagnation region momentum and thermal flow fields were investigated. This analysis predicted the existing measurements of the stagnation region mean velocity and surface heat transfer rate with cross flow Reynolds number and turbulence intensity less than 250,000 and 0.05, respectively.

Wang, C. R.↗
NASA NTRS1987

Two- and three-dimensional viscous computations of a hypersonic inlet flow

The three-dimensional parabolized Navier-Stokes code has been used to investigate the flow through a Mach 7.4 inlet. A two-dimensional parametric study of grid resolution, turbulence modeling and effect of gamma has been done and compared with experimental results. The results show that mesh resolution of the shock waves, real gas effects and turbulence length scaling are very important to get accurate results for hypersonic inlet flows. In addition a three-dimensional calculation of the Mach 7.4 inlet has been done on a straight sideplate configuration. The results show that the glancing shock/boundary layer interaction phenomena causes significant three-dimensional flow in the inlet.

Kunik, William G.↗
NASA NTRS1987

Two- and three-dimensional viscous computations of a hypersonic inlet flow

The three-dimensional parabolized Navier-Stokes code has been used to investigate the flow through a Mach 7.4 inlet. A two-dimensional parametric study of grid resolution, turbulence modeling and effect of gamma has been done and compared with experimental results. The results show that mesh resolution of the shock waves, real gas effects and turbulence length scaling are very important to get accurate results for hypersonic inlet flows. In addition a three-dimensional calculation of the Mach 7.4 inlet has been done on a straight sideplate configuration. The results show that the glancing shock/boundary layer interaction phenomena causes significant three-dimensional flow in the inlet.

Kunik, William G.↗
NASA NTRS1988

Simulations using a 3D-moment turbulence closure in a two-dimensional cloud model - Comparisons with observations

A method for deep cumulus cloud simulations is presented, in which the plumes or clouds are simulated using a two-dimensional axisymmetric numerical cloud model with third-moment turbulence closure. This method is based on the assumption that the only uncertainty in the turbulence closure is in the turbulent length scale specification. Therefore, the length scale is varied, and, using different length scales, the simulated plumes are compared to the observed ones, making it possible to optimize the length scale prescription.

Krueger, Steven K.↗
NASA NTRS1990

On the possibility of an alpha-sq omega-type dynamo in a thin layer inside the sun

If the solar dynamo operates in a thin layer of 10,000-km thickness at the interface between the convection zone and the radiative core, using the facts that the dynamo should have a period of 22 years and a half-wavelength of 40 deg in the theta-direction, it is possible to impose restrictions on the values which various dynamo parameters are allowed to have. It is pointed out that the dynamo should be of alpha-sq omega nature, and kinematical calculations are presented for free dynamo waves and for dynamos in thin rectangular slabs with appropriate boundary conditions. An alpha-sq omega dynamo is expected to produce a significant poloidal field which does not leak to the solar surface. It is found that the turbulent diffusity eta and alpha-coefficient are restricted to values within about a factor of 10, the median values being eta of about 10 to the 10th sq cm/sec and alpha of about 10 cm/sec. On the basis of mixing length theory, it is pointed out that such values imply a reasonable turbulent velocity of the order 30 m/s, but rather small turbulent length scales like 300 km.

Choudhuri, Arnab Rai↗
NASA NTRS1991

Direct numerical simulation and analysis of shock turbulence interaction

Two kinds of linear analysis, rapid distortion theory (RDT) and linear interaction analysis (LIA), were used to investigate the effects of a shock wave on turbulence. Direct numerical simulations of two-dimensional isotropic turbulence interaction with a normal shock were also performed. The results from RDT and LIA are in good agreement for weak shock waves, where the effects of shock front curvature and shock front unsteadiness are not significant in producing vorticity. The linear analyses predict wavenumber-dependent amplification of the upstream one-dimensional energy spectrum, leading to turbulence scale length scale decrease through the interaction. Instantaneous vorticity fields show that vortical structures are enhanced while they are compressed in the shock normal direction. Entrophy amplfication through the shock wave compares favorably with the results of linear analyses.

Lee, Sangsan↗