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

N-S Simulations of Crow-Type Instabilities in Vortex Wake

Three-dimensional Navier-Stokes simulations of the Crow instability of wake vortices are conducted using large-eddy simulations without background turbulence. Sinusoidal displacement has been specified as the initial perturbation for the vortex system. The results have shown that the minimum Crow instability wavelength is about one vortex spacing shorter than predicted by Crow's linear stability theory. The planar- standing-wave-angle value and the amplitude indifference behavior agree with Crow's analysis. Simulations with periodicity in the axial direction have indicated minor influence of axial flow on the Crow instability.

Zheng, Z. C.↗

Nosetip Bluntness Effects on Transition at Hypersonic Speeds: Experimental and Numerical Analysis

The existing database of transition measurements in hypersonic ground facilities has established that the onset of boundary layer transition over a circular cone at zero angle of attack shifts downstream as the nosetip bluntness is increased with respect to a sharp cone. However, this trend is reversed at sufficiently large values of the nosetip Reynolds number, so that the transition onset location eventually moves upstream with a further increase in nosetip bluntness. This transition reversal phenomenon, which cannot be explained on the basis of linear stability theory, was the focus of a collaborative investigation under the NATO STO group AVT-240 on Hypersonic Boundary-Layer Transition Prediction. The current paper provides an overview of that effort, which included wind tunnel measurements in three different facilities and theoretical analysis related to modal and nonmodal amplification of boundary layer disturbances. Because neither first and second-mode waves nor entropy-layer instabilities are found to be substantially amplified to initiate transition at large bluntness values, transient (i.e., nonmodal) disturbance growth has been investigated as the potential basis for a physics based model for the transition reversal phenomenon. Results of the transient growth analysis indicate that stationary disturbances that are initiated within the nosetip or in the vicinity of the juncture between the nosetip and the frustum can undergo relatively significant nonmodal amplification and that the maximum energy gain increases nonlinearly with the nose radius of the cone. This finding does not provide a definitive link between transient growth and the onset of transition, but it is qualitatively consistent with the experimental observations that frustum transition during the reversal regime was highly sensitive to wall roughness, and furthermore, was dominated by disturbances that originated near the nosetip. Furthermore, the present analysis shows significant nonmodal growth of traveling disturbances that peak within the entropy layer and could also play a role in the transition reversal phenomenon.

Freestream Mach Number↗

Field Effects of Buoyancy on Lean Premixed Turbulent Flames

The study of field effects of buoyancy on premixed turbulent flames is directed towards the advancement of turbulent combustion theory and the development of cleaner combustion technologies. Turbulent combustion is considered the most important unsolved problem in combustion science and laboratory studies of turbulence flame processes are vital to theoretical development. Although buoyancy is dominant in laboratory flames, most combustion models are not yet capable to consider buoyancy effects. This inconsistency has impeded the validation of theories and numerical simulations with experiments. Conversely, the understanding of buoyancy effects is far too limited to help develop buoyant flame models. Our research is also relevant to combustion technology because lean premixed combustion is a proven method to reduce the formation of oxides of nitrogen (NOx). In industrial lean premixed combustion systems, their operating conditions make them susceptible to buoyancy thus affecting heat distribution, emissions, stability, flashback and blowoff. But little knowledge is available to guide combustion engineers as to how to avoid or overcome these problems. Our hypothesis is that through its influence on the mean pressure field, buoyancy has direct and indirect effects on local flame/turbulence interactions. Although buoyancy acts on the hot products in the farfield the effect is also felt in the nearfield region upstream of the flame. These changes also influence the generation and dissipation of turbulent kinetic energy inside the flame brush and throughout the flowfield. Moreover, the plume of an open flame is unstable and the periodic fluctuations make additional contributions to flame front dynamics in the farfield. Therefore, processes such as flame wrinkling, flow acceleration due to heat release and flame- generated vorticity are all affected. Other global flame properties (e.g. flame stabilization limits and flame speed) may all be coupled to buoyancy. This problem poses major challenges to combustion modeling due to its need for a computation domain extending into the farfield and full specifications of upstream, wall and downstream boundary conditions.

Cheng, R. K.↗

The hydrodynamic role of fish squamosal integument as an analog of surfaces directly formed by turbulent flow. Report 1: Similarity between irregularities in squamosal integument and those on surfaces formed by river bed flow

The distribution of squamae on the fish body and that of the deposits in the bed of the river can be described by the same equation. The curves reflecting the relative elongation and stability of the body shape of the fish continue the curves showing the elongation of bank spit and stability of the bed of the river.

Kudryashov, A. F.↗

Containerless Ripple Turbulence

One of the longest standing unsolved problems in physics relates to the behavior of fluids that are driven far from equilibrium such as occurs when they become turbulent due to fast flow through a grid or tidal motions. In turbulent flows the distribution of vortex energy as a function of the inverse length scale [or wavenumber 'k'] of motion is proportional to 1/k(sup 5/3) which is the celebrated law of Kolmogorov. Although this law gives a good description of the average motion, fluctuations around the average are huge. This stands in contrast with thermally activated motion where large fluctuations around thermal equilibrium are highly unfavorable. The problem of turbulence is the problem of understanding why large fluctuations are so prevalent which is also called the problem of 'intermittency'. Turbulence is a remarkable problem in that its solution sits simultaneously at the forefront of physics, mathematics, engineering and computer science. A recent conference [March 2002] on 'Statistical Hydrodynamics' organized by the Los Alamos Laboratory Center for Nonlinear Studies brought together researchers in all of these fields. Although turbulence is generally thought to be described by the Navier-Stokes Equations of fluid mechanics the solution as well as its existence has eluded researchers for over 100 years. In fact proof of the existence of such a solution qualifies for a 1 M$ millennium prize. As part of our NASA funded research we have proposed building a bridge between vortex turbulence and wave turbulence. The latter occurs when high amplitude waves of various wavelengths are allowed to mutually interact in a fluid. In particular we have proposed measuring the interaction of ripples [capillary waves] that run around on the surface of a fluid sphere suspended in a microgravity environment. The problem of ripple turbulence poses similar mathematical challenges to the problem of vortex turbulence. The waves can have a high amplitude and a strong nonlinear interaction. Furthermore, the steady state distribution of energy again follows a Kolmogorov scaling law; in this case the ripple energy is distributed according to 1/k (sup 7/4). Again, in parallel with vortex turbulence ripple turbulence exhibits intermittency. The problem of ripple turbulence presents an experimental opportunity to generate data in a controlled, benchmarked system. In particular the surface of a sphere is an ideal environment to study ripple turbulence. Waves run around the sphere and interact with each other, and the effect of walls is eliminated. In microgravity this state can be realized for over 2 decades of frequency. Wave turbulence is a physically relevant problem in its own right. It has been studied on the surface of liquid hydrogen and its application to Alfven waves in space is a source of debate. Of course, application of wave turbulence perspectives to ocean waves has been a major success. The experiment which we plan to run in microgravity is conceptually straightforward. Ripples are excited on the surface of a spherical drop of fluid and then their amplitude is recorded with appropriate photography. A key challenge is posed by the need to stably position a 10cm diameter sphere of water in microgravity. Two methods are being developed. Orbitec is using controlled puffs of air from at least 6 independent directions to provided the positioning force. This approach has actually succeeded to position and stabilize a 4cm sphere during a KC 135 segment. Guigne International is using the radiation pressure of high frequency sound. These transducers have been organized into a device in the shape of a dodecahedron. This apparatus 'SPACE DRUMS' has already been approved for use for combustion synthesis experiments on the International Space Station. A key opportunity presented by the ripple turbulence data is its use in driving the development of codes to simulate its properties.

Putterman, Seth↗

New insights on wind turbine wakes from large‐eddy simulation: Wake contraction, dual nature, and temperature effects

Abstract Large‐eddy simulation (LES) has been adopted to study wind turbine wakes because it can capture fine‐scale details of turbulent wind flows and interactions with wind turbines. Here, we use the LES version of the Weather Research and Forecasting (WRF) model with an actuator disk model to gain insights on several wake effects that have been traditionally difficult to measure. The first finding is that the wake has a “dual nature,” meaning that the wind speed deficit behaves differently from the added turbulent kinetic energy (TKE) and the two are not co‐located in space. For example, the wind speed deficit peaks at hub height and reaches the ground within 8D (D is the rotor diameter), but added TKE peaks near the rotor tip and generally remains aloft. Second, temperature changes near the ground are driven by the added TKE in the rotor area and by atmospheric stability. The combination of these two factors determines the sign and intensity of the vertical heat flux divergence below the rotor, with convergence and warming associated with stable conditions and weak divergence and modest cooling with unstable conditions. Third, wakes do not expand indefinitely, as suggested by similarity theory applied to the wind speed deficit, but eventually stop expanding and actually contract, at different rates depending on atmospheric stability. The implication of these findings is that, in order to study wakes, it is not sufficient to focus on wind speed deficit alone, because TKE is also important and yet behaves differently from the wind speed deficit.

17 WIND ENERGY↗

Modeling Arctic boundary layer cloud streets at grey-zone resolutions

To better understand how model resolution affects the formation of Arctic boundary layer clouds, we investigated the influence of grid spacing on simulating cloud streets that occurred near Utqiagvik (formerly Barrow), Alaska on 2 May 2013 and were observed by the Moderate Resolution Imaging Spectroradiometer (MODIS). The Weather Research and Forecasting model was used to simulate the clouds using nested domains with increasingly fine resolution ranging from a horizontal grid spacing of 27 km in the boundary-layer-parameterized mesoscale domain to a grid spacing of 0.111 km in the large-eddy-permitting domain. We investigated the model-simulated mesoscale environment, horizontal and vertical cloud structures, boundary layer stability, and cloud properties, all of which were subsequently used to interpret the observed roll-cloud case. Increasing model resolution led to a transition from a more buoyant boundary layer to a more shear-driven turbulent boundary layer. The clouds were stratiform-like in the mesoscale domain, but as the model resolution increased, roll-like structures, aligned along the wind field, appeared with ever smaller wavelengths. A stronger vertical water vapor gradient occurred above the cloud layers with decreasing grid spacing. With fixed model grid spacing at 0.333 km, changing the model configuration from a boundary layer parameterization to a large-eddy-permitting scheme produced a more shear-driven and less unstable environment, a stronger vertical water vapor gradient below the cloud layers, and the wavelengths of the rolls decreased slightly. In this study, only the large-eddy-permitting simulation with gird spacing of 0.111 km was sufficient to model the observed roll clouds.

Lai, Hui-Wen↗

Study of Stratus-Lowering Marine-Fog Events Observed During C-FOG

Two stratus lowering marine fog events observed on 28 September and 04 October 2018 during the Coastal Fog (C-Fog) field campaign that took place offshore of Eastern Canada during 1 September to 6 October 2018 are described. In-situ, profiling and remote sensing observations were made at selected land sites in eastern Newfoundland (NL) and Nova Scotia (NS) as well as aboard the research vessel (R/V) Hugh R. Sharp that cruised in adjoining coastal waters. Synoptic scale analysis showed that both fog episodes were an outcome of the interaction between synoptic-scale surface level low-pressure systems and a contiguous high-pressure system. At the same time, back trajectories revealed that the bulk of the fog layer is formed due to differential advection. Here, the diameter of the fog droplets at the surface gradually decreased from the centre of the fog layer to its leading/trailing edges. The bimodal fog droplet diameter distribution with peaks at 5-10 µm and 20-25 µm provided clues on droplet collision and coalescence processes. The observed difference between microphysical variables and droplet distribution between the two fog events and within same fog layer may have been governed by atmospheric boundary layer conditions (e.g., humidity conditions and turbulence) prevailed in the fog layer. Overall, it is concluded that the life cycle of observed stratus-lowering coastal-fog episodes is dependent on synoptic conditions as well as atmospheric boundary layer characteristics such as stability, cloud top cooling and entrainment.

54 ENVIRONMENTAL SCIENCES↗

Generation and development of small-amplitude disturbances in a laminar boundary layer in the presence of an acoustic field

A low-turbulence subsonic wind tunnel was used to study the influence of acoustic disturbances on the development of small sinusoidal oscillations (Tollmien-Schlichting waves) which constitute the initial phase of turbulent transition. It is found that acoustic waves propagating opposite to the flow generate vibrations of the model (plate) in the flow. Neither the plate vibrations nor the acoustic field itself have any appreciable influence on the stability of the laminar boundary layer. The influence of an acoustic field on laminar boundary layer disturbances is limited to the generation of Tollmien-Schlichting waves at the leading-edge of the plate.

Kachanov, Y. S.↗

Development of an explicit multigrid algorithm for quasi-three-dimensional viscous flows in turbo-machinery

A rapid quasi three-dimensional analysis was developed for blade-to-blade flows in turbomachinery. The analysis solves the unsteady Euler or thin layer Navier-Stokes equations in a body-fitted coordinate system. It accounts for the effects of rotation, radius change, and stream-surface thickness. The Baldwin-Lomax eddy-viscosity model is used for turbulent flows. The equations which are solved by a two-stage Runge-Kutta scheme made efficient by use of vectorization, a variable time-step, and a flux-based multigrid scheme, are described. A stability analysis is presented for the two-stage scheme. Results for a flat-plate model problem show the applicability of the method to axial, radial, and rotating geometries. Results for a centrifugal impeller and a radial diffuser show that the quasi three-dimensional viscous analysis can be a practical design tool.

Chima, R. V.↗

Development of an explicit multigrid algorithm for quasi-three-dimensional viscous flows in turbomachinery

A rapid quasi three-dimensional analysis was developed for blade-to-blade flows in turbomachinery. The analysis solves the unsteady Euler or thin layer Navier-Stokes equations in a body-fitted coordinate system. It accounts for the effects of rotation, radius change, and stream-surface thickness. The Baldwin-Lomax eddy-viscosity model is used for turbulent flows. The equations which are solved b a two-stage Runge-Kutta scheme made efficient by use of vectorization, a variable time-step, and a flux-based multigrid scheme, are described. A stability analysis is presented for the two-stage. Results for a flat-plate model problem show the applicability of the method to axial, radial, and rotating geometries. Results for a centrifugal impeller and a radial diffuser show that the quasi three-dimensional viscous analysis can be a practical design tool.

Chima, R. V.↗

Convective heat flow in space cryogenics plugs - Critical and moderate He II heat flux densities

Plug flow rates of entropy, heat and normal fluid in phase separators and in zero net mass flow systems are, to some extent, quite similar. A simplified analysis of critical conditions is presented in agreement with data trends. A critical temperature gradient arises on the basis of the He II two-fluid model at the stability limit constraining the thermohydrodynamics of the system. Thus, the question of critical thermodynamic fluctuations associated with nucleation versus the possibility of critical gradients in externally imposed parameters is answered in favor of the latter route toward turbulence. Furthermore, a similarity equation is presented which incorporates size dependent rates for moderate heat flow densities observed in experiments.

Yuan, S. W. K.↗

Does ion tearing exist?

Theoretical models of reconnection in the geomagnetic tail and its relationship to the collisionless tearing mode are examined analytically. The focus is on attempts to remove the stabilizing effects of electron compressibility in the quasi-neutral sheet. It is shown that this goal can be achieved by introducing spatial diffusion of electrons across the magnetic-flux surfaces, but not by wave-turbulent pitch-angle scattering (Coroniti, 1980) or nonadiabatic stochastic first-invariant diffusion (Buechner and Zelenyi, 1988). Hence the ion tearing mode does not exist, but the electron tearing mode does.

Pellat, R.↗

Alleviation of side force on tangent-ogive forebodies using passive porosity

An experimental investigation to determine the effectiveness of porosity for alleviating side forces on forebodies was conducted in the NASA Langley Research Center 7 x 10 ft High-Speed Wind Tunnel. The study consisted of a comparison of experimental force, moment, and surface pressure results obtained on a fineness ratio 5.0 tangent-ogive porous forebody with 0.020 in. hole diameter and 22 percent porosity with results obtained on a solid forebody. The forebodies were tested with cylindrical afterbodies. The solid forebody was tested with transition grit to simulate fully turbulent conditions and without transition grit to simulate free transition conditions. The extent of porosity on the forebody was varied to determine the extent of porosity needed to alleviate side forces. Static longitudinal and lateral-directional stability and surface pressure data were obtained at Mach numbers of 0.2, 0.5, and 0.8. The angle of attack range was from 5 to 45 deg and roll angles from -90 to 180 deg. The solid forebody exhibited large asymmetries at moderate to high angles of attack causing large side forces and yawing moments. The porous forebody exhibited no significant side forces or yawing moments at any angle of attack tested.

Bauer, Steven X. S.↗

Computational fluid mechanics

Two papers are included in this progress report. In the first, the compressible Navier-Stokes equations have been used to compute leading edge receptivity of boundary layers over parabolic cylinders. Natural receptivity at the leading edge was simulated and Tollmien-Schlichting waves were observed to develop in response to an acoustic disturbance, applied through the farfield boundary conditions. To facilitate comparison with previous work, all computations were carried out at a free stream Mach number of 0.3. The spatial and temporal behavior of the flowfields are calculated through the use of finite volume algorithms and Runge-Kutta integration. The results are dominated by strong decay of the Tollmien-Schlichting wave due to the presence of the mean flow favorable pressure gradient. The effects of numerical dissipation, forcing frequency, and nose radius are studied. The Strouhal number is shown to have the greatest effect on the unsteady results. In the second paper, a transition model for low-speed flows, previously developed by Young et al., which incorporates first-mode (Tollmien-Schlichting) disturbance information from linear stability theory has been extended to high-speed flow by incorporating the effects of second mode disturbances. The transition model is incorporated into a Reynolds-averaged Navier-Stokes solver with a one-equation turbulence model. Results using a variable turbulent Prandtl number approach demonstrate that the current model accurately reproduces available experimental data for first and second-mode dominated transitional flows. The performance of the present model shows significant improvement over previous transition modeling attempts.

Hassan, H. A.↗

Numerical simulation of supersonic boundary layer transition

The present contribution reviews some of the recent progress obtained at our group in the direct numerical simulation (DNS) of compressible boundary layer transition. Elements of the different simulation approaches and numerical techniques employed are surveyed. Temporal and spatial simulations, as well as comparisons with results obtained from Parabolized Stability Equations, are discussed. DNS results are given for flat plate boundary layers in the Mach number range 1.6 to 4.5. A temporal DNS at Mach 4.5 has been continued through breakdown all the way to the turbulent stage. In addition results obtained with a recently developed extended temporal DNS approach are presented, which takes into account some nonparallel effects of a growing boundary layer. Results from this approach are quite close to those of spatial DNS, while preserving the efficiency of the temporal DNS.

Guo, Y.↗

Novel Epoxy Particulate Composites for Mitigation of Insect Residue Adhesion on Future Aircraft Surfaces

Drag is reduced significantly for airflow over surfaces when laminar flow can be maintained over greater chord lengths, the distance from the leading edge of an airfoil.1 However, surface imperfections, such as chipped paint, scratches, and events that change topography on a microscopic scale can introduce airflow instabilities resulting in premature transition to turbulent flow.1 Although many of these surface imperfections can be avoided with proper maintenance, advanced materials, and advanced manufacturing practices, topographical surface anomalies arising during flight from insect impacts cannot be controlled and can influence laminar flow stability. Practical solutions to this operational challenge need to be developed for future aircraft to have full advantage of laminar flow designs that improve fuel efficiency.2 Researchers have investigated various methods to mitigate insect residue adhesion for decades.3 Although several techniques have demonstrated efficacy including mechanical scrapers, active liquid discharge systems, and sacrificial paper coatings, they have not been commercially implemented due to increased manufacturing and operational complexity, environmental impact, and weight penalties. Coatings offer a simple route for passive insect residue adhesion prevention without many of the challenges associated with maintenance of laminar flow.4 In our previous work, we determined that most commercially available materials were not effective at insect residue adhesion.5 We also identified improvements when both surface energy could be controlled by surface modifying agents and the topography could be altered through the use of micron-sized and nanometer-sized filler materials.6 In this work, these general principles were applied to an epoxy system to evaluate the behavior of the surface modifying agent, a fluorinated alkyl ether oligomer, on surface energy and insect residue adhesion properties.

Wohl, Christopher J.↗

On the role of ion temperature gradient turbulence in driving ion thermal transport in neutral beam injection-heated L-mode plasmas in a superconducting tokamak

In this paper, we report a detailed experimental study of the role of ion-temperature-gradient driven turbulence in driving ion thermal transport in neutral beam injection (NBI)-heated L-mode plasmas in Experimental Advanced Superconducting Tokamak (EAST) (Wan et al 2000 Nucl. Fusion 40 1057) coupled with linear and nonlinear gyrokinetic simulations. Significant ion-scale turbulence k ⟂ < cm −1 (k ⟂ < 1.5, where k ⟂ where is the perpendicular wavenumber and ρ s is the ion gyroradius calculated using local electron temperature T e ), measured with a microwave reflectometer, is observed in the plasma core at for r/a ≈ 0.265 and 0.5 (where r is half the diameter of the closed flux surface at a given radial position, and a is half the diameter of the last closed flux surface). Local linear stability analysis with experimental equilibrium quantities at these two radial locations using the GS2 gyrokinetic code shows that the most unstable ion-scale micro-instability is the ion temperature gradient (ITG) mode. Since the computed maximum ion-scale linear growth rates are larger than the local Waltz-Miller E × B shearing rate (Waltz and Miller 1999 Phys. Plasmas 6 4265), the E × B shear is unable to suppress ion-scale turbulence, consistent with the experimental observation of ion-scale turbulence. Experimental ion and electron thermal transports, calculated with power balance analysis, are both anomalous with ion thermal transport being the dominant channel. Nonlinear gyrokinetic simulations using the electrostatic global particle-in-cell gyrokinetic tokamak simulation (GTS) code (Weixing Wang et al 2010 Phys. Plasmas 17 072511) show good quantitative agreement between predicted ion thermal transport and that from the power balance analysis. These results demonstrate for the first time the important role of ITG turbulence in driving thermal transport in NBI-heated L-mode plasmas in EAST, showing the effectiveness of GTS code in simulating electrostatic turbulence in EAST and the need of suppressing ITG turbulence to improve energy confinement in EAST.

Sun, P. J. [Chinese Academy of Sciences (CAS), Hef↗