Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “shear flow”

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 505 records · Page 28

Analysis of the leading edge effects on the boundary layer transition

A general theory of boundary layer control by surface heating is presented. Some analytical results for a simplified model, i.e., the optimal control of temperature fluctuations in a shear flow are described. The results may provide a clue to the effectiveness of the active feedback control of a boundary layer flow by wall heating. In a practical situation, the feedback control may not be feasible from the instrumentational point of view. In this case the vibrational control introduced in systems science can provide a useful alternative. This principle is briefly explained and applied to the control of an unstable wavepacket in a parallel shear flow.

Chow, Pao-Liu↗

Perturbations of the Richardson number field by gravity waves

An analytic solution is presented for a stratified fluid of arbitrary constant Richardson number. By computer aided analysis the perturbation fields, including that of the Richardson number can be calculated. The results of the linear analytic model were compared with nonlinear simulations, leading to the following conclusions: (1) the perturbations in the Richardson number field, when small, are produced primarily by the perturbations of the shear; (2) perturbations of in the Richardson number field, even when small, are not symmetric, the increase being significantly larger than the decrease (the linear analytic solution and the nonlinear simulations both confirm this result); (3) as the perturbations grow, this asymmetry increases, but more so in the nonlinear simulations than in the linear analysis; (4) for large perturbations of the shear flow, the static stability, as represented by N2, is the dominating mechanism, becoming zero or negative, and producing convective overturning; and (5) the convectional measure of linearity in lee wave theory, NH/U, is no longer the critical parameter (it is suggested that (H/u sub 0) (du sub 0/dz) takes on this role in a shearing flow).

Wurtele, M. G.↗

Mitigation of resistive drift wave and ion temperature gradient instabilities by velocity shear

The effects of velocity shear on the resistive drift wave instability in the non-adiabatic limit and the toroidal ion temperature gradient instability are investigated for a plasma of inhomogeneous density/ion temperature, respectively. For the resistive drift wave, we find that the instability growth rate decreases monotonically with increasing magnitude of shear, but we find that complete stabilization is impossible. For the ion temperature gradient instability, we find that the standard WKB approximation is insufficient to describe the full behavior of the instability, and that an analysis of the localized eigenmode problem reveals two separate unstable solutions, which the WKB approximation does not predict. The impact of flow shear on these two new unstable solutions is discussed. In both resistive drift wave and ion temperature gradient instabilities, the sheared flow causes a shifting, tilting, and sharpening of the electrostatic potential eddies.

Physics↗

Time dependent turbulence modeling and analytical theories of turbulence

By simplifying the direct interaction approximation (DIA) for turbulent shear flow, time dependent formulas are derived for the Reynolds stresses which can be included in two equation models. The Green's function is treated phenomenologically, however, following Smith and Yakhot, we insist on the short and long time limits required by DIA. For small strain rates, perturbative evaluation of the correlation function yields a time dependent theory which includes normal stress effects in simple shear flows. From this standpoint, the phenomenological Launder-Reece-Rodi model is obtained by replacing the Green's function by its long time limit. Eddy damping corrections to short time behavior initiate too quickly in this model; in contrast, the present theory exhibits strong suppression of eddy damping at short times. A time dependent theory for large strain rates is proposed in which large scales are governed by rapid distortion theory while small scales are governed by Kolmogorov inertial range dynamics. At short times and large strain rates, the theory closely matches rapid distortion theory, but at long times it relaxes to an eddy damping model.

Rubinstein, R.↗

The mean logarithm emerges with self-similar energy balance.

The attached eddy hypothesis of Townsend (The Structure of Turbulent Shear Flow, 1956, Cambridge University Press) states that the logarithmic mean velocity admits self-similar energy-containing eddies which scale with the distance from the wall. Over the past decade, there has been a significant amount of evidence supporting the hypothesis, placing it to be the central platform for the statistical description of the general organisation of coherent structures in wall-bounded turbulent shear flows. Nevertheless, the most fundamental question, namely why the hypothesis has to be true, has remained unanswered over many decades. Under the assumption that the integral length scale is proportional to the distance from the wall y , in the present study we analytically demonstrate that the mean velocity is a logarithmic function of y if and only if the energy balance at the integral length scale is self-similar with respect to y , providing a theoretical basis for the attached eddy hypothesis. The analysis is subsequently verified with the data from a direct numerical simulation of incompressible channel flow at the friction Reynolds number Reτ≃5200 (Lee & Moser, J. Fluid Mech., vol. 774, 2015, pp. 395–415).

42 ENGINEERING↗

Coarse-Grained Simulations of Polyrotaxane Hydrogels under Quiescent and Shear Conditions

Cyclodextrin-based polyrotaxanes (PR) form hydrogels in water through cyclodextrin (CD) aggregation and crystallization. These networks can break under shear flows, making them versatile platforms for extrusion-based 3D printing. To optimize the material properties of 3D-printed PR gels, a microscopic understanding of the structural evolution during 3D printing is necessary. Here, we employ coarse-grained (CG) simulations to reveal the PR assembly process at the molecular level. Our simulations reproduce the experimental crystal morphologies of PR at varying concentrations and chain lengths under quiescent conditions. Using nonequilibrium simulations, we show shear flow ruptures crystalline domain connectivity in PR gels and stacks the lamellae in the gradient direction, allowing the materials to flow during 3D printing. After the cessation of flow, the anisotropic crystal alignment and absence of available dangling PR diminish the intercrystal connectivity. The printed materials are therefore mechanically weaker than the pristine hydrogels, in agreement with experimental results. Nonetheless, by relating the microscopic structural evolution with viscoelastic properties of PR gels and solutions, we elucidate how flow conditions and sample composition affect the 3D printing performance of PR hydrogels.

Smith, Cameron D. [Dartmouth College, Hanover, NH ↗

Approach to the origin of turbulence on the basis of two-point kinetic theory

Equations for the fluctuation correlation in an incompressible shear flow are derived on the basis of kinetic theory, utilizing the two-point distribution function which obeys the BBGKY hierarchy equation truncated with the hypothesis of 'ternary' molecular chaos. The step from the molecular to the hydrodynamic description is accomplished by a moment expansion which is a two-point version of the thirteen-moment method, and which leads to a series of correlation equations, viz., the two-point counterparts of the continuity equation, the Navier-Stokes equation, etc. For almost parallel shearing flows the two-point equation is separable and reduces to two Orr-Sommerfeld equations with different physical implications.

Tsuge, S.↗

The instability of sheared liquid layers

The critical values of the Reynolds number (Rc) above which unstable travelling waves exist in a steady shear flow are calculated with the inclusion of a term in the normal stress boundary condition omitted by Miles (1960). Results show a value of Rc = 34.2, compared to Miles' value of 203, for the case of a zero-surface tension interface. The critical wavenumber, critical phase speed, and Rc are calculated as functions of a nondimensional surface tension, and the mechanism of the instability is investigated by examining the disturbance energy equation. A long wave instability is found at small values of the Reynolds number when the shear flow has a parabolic velocity profile. In addition, these results are extended to larger values of the wavenumber by numerical methods, and the relation between this long wave instability and the profile curvature is examined.

Smith, M. K.↗

Local isotropy in distorted turbulent boundary layers at high Reynolds number

This is a report on the continuation of our experimental investigations of the hypothesis of local isotropy in shear flows. This hypothesis, which states that at sufficiently high Reynolds numbers the small-scale structures of turbulent motions are independent of large-scale structures and mean deformations, has been used in theoretical studies of turbulence and computational methods such as large-eddy simulation. Since Kolmogorov proposed his theory, there have been many experiments, conducted in wakes, jets, mixing layers, a tidal channel, and atmospheric and laboratory boundary layers, in which attempts have been made to verify - or refute - the local-isotropy hypothesis. However, a review of the literature over the last five decades indicated that, despite all these experiments in shear flows, there was no consensus in the scientific community regarding this hypothesis, and, therefore, it seemed worthwhile to undertake a fresh experimental investigation into this question.

Saddoughi, Seyed G.↗

Cross phases of temperature-gradient-driven turbulence as a model basis for I -mode particle transport

As a model for understanding the type of transport behavior characteristic of the tokamak I mode, cross-phase physics for particle-transport is studied analytically for turbulence dominated by either ion-temperature-gradient (ITG) or electron-temperature-gradient (ETG) instability. I mode is a transport-barrier regime of reduced thermal transport but essentially unaffected particle transport. It is assumed that ITG turbulence applies to the baseline L mode, ETG to I mode, and that E × B flow shear is stronger in I mode, lowering all fluxes. In ITG turbulence, particle transport is governed by trapped electrons. Sensitivity to collisions produces the well-known temperature-gradient-driven pinch that offsets density-gradient-driven outward diffusion, weakening particle transport in L mode. In ETG turbulence, nonadiabatic ions are collisionless. Nonzero transport requires an ion spectrum feature whose magnetic-drift resonance supplies the necessary cross phase. If frequencies of order the ion diamagnetic drift frequency dominate the ion part of the spectrum, as would occur with weakly unstable ITG turbulence, all components of the particle transport are outward and can offset flow-shear-induced flux reductions to produce a flux that is similar to the ITG L-mode particle flux. Nonlinear frequencies are potentially relevant and discussed in relation to I mode.

Physics↗

The structure of a passive scalar field with a uniform mean gradient in rapidly sheared homogeneous turbulent flow

The effect of an arbitrarily oriented mean passive scalar gradient on one-point passive scalar statistics is studied in homogeneous turbulent shear flow in the limit of rapid shearing. By neglecting the nonlinear inertial transfer to small scales an analytical solution for individual Fourier modes is obtained for the case of unity Prandtl number. This solution is used to compute the development of one-point statistics for the velocity and scalar fields in the inviscid limit. Comparisons to direct numerical simulations of the full nonlinear equations for the same flow show that, in addition to describing the early time response to the imposed shear, the linear solution gives reasonable estimates of several correlation coefficients for the developed shear flow.

Rogers, Michael M.↗

Structure Formation Through Magnetohydrodynamical Instabilities in Primordial Disks

The shear flow instabilities under the presence of magnetic fields in the primordial disk can greatly facilitate the formation of density structures that serve as seeds prior to the onset of the gravitational Jeans instability. We evaluate the effects of the Parker, magnetorotational and kinematic dynamo instabilities by comparing the properties of these instabilities. We calculate the mass spectra of coagulated density structures by the above mechanism in the radial direction for an axisymmetric magnetohydrodynamic (MHD) torus equilibrium and power density profile models. Our local three-dimensional MHD simulation indicates that the coupling of the Parker and magnetorotational instability creates spiral arms and gas blobs in an accretion disk, reinforcing the theory and model. Such a mechanism for the early structure formation may be tested in a laboratory. The recent progress in experiments involving shear flows in rotating tokamak, field reversed configuration (FRC) and laser plasmas may become a key element to advance in nonlinear studies.

Noguchi, Koichi↗

Simulation of Li-pellet triggered ELMs in EAST with an impurity model implemented under BOUT++ framework

A simple impurity model has been developed under the BOUT++ framework to investigate the Li-pellet triggered edge localized mode (ELM) in EAST configuration. The present impurity model decouples the ion pressure enhancement effect (IPEE) (which affects radial E X B flow shear and gyro-viscosity of deuterium ions and is implicitly included in our previous work) to make the simulation of Li pellet injection more physically reasonable. In addition, the present impurity model also includes the impurity equilibrium effect (IEE) (which induces modifications on vorticity and gyro-viscosity). The simulation results show that without IPEE a turbulent ELM induced by Li pellet is triggered by multiple peeling-ballooning modes (PBMs) rather than by a single dominant mode, and the simultaneous growth of multiple modes is conducive to reducing the time for pedestal entering the energy loss state. When a turbulent ELM occurs, the nonlinearly dominant modes undergo a secondary fast growth during fast crash phase. The finding explains the secondary increase of D α emission observed in the DIII-D Li injection experiments. During the evolution of a turbulent ELM, the continuous pedestal collapse dominates the pedestal energy loss in turbulent transport phase; however, in saturation phase, the stabilizing effects on PBMs (n ≠ 0) by n = 0 mode and radial E X B flow shear within the steepest gradient region are more prominent. It is also found that with IEE and without IPEE, the size threshold of ELM triggering and the magnitude of ELM sizes both show a good agreement with EAST and DIII-D experimental observations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

On sublayer streaks

Turbulent sublayer streaks are studied with the aid of a simplified theoretical model. In this the nonlinear activity is assumed to be intermittent and to act locally in space during a very short initial time so as to set up the initial conditions for the subsequent linear and inviscid evolution of the resulting three-dimensional flow disturbance. The mean shear flow is taken as a parallel one and a correction for the long-tern effects of viscosity is applied. A model for the initial nonlinear phase is chosen to represent the local Reynolds stresses that would be produced by a patch of local inflectional instability. The streamwise dimension of the resulting eddy is found to grow linearly with time in accordance with the algebraic instability mechanism (Landahl 1980). The associated Reynolds shear stress is expressible in a simple manner in terms of the liftup of the fluid elements and is suggestive of an algebraic-type Reynolds stress model similar to, but not identical to, that of Prandtl's (1925) mixing-length theory.

Landahl, M. T.↗

Direct measurements of the 3D plasma velocity in single-helical-axis RFP plasmas

The first local velocity measurements of helical equilibrium plasmas in the Reversed Field Pinch (RFP) Single Helical Axis (SHAx) state using a Charge Exchange Recombination Spectroscopy (CHERS) diagnostic are presented. Measurements show strong axisymmetric and non-axisymmetric flow, with n ¼ 5 components of flow related to the ðm; nÞ¼ð1; 5Þ dominant magnetic mode on the order of the axisymmetric flow in certain regions of the plasma, as well as significant n > 5 flow. Flow measurements are compared with NIMROD simulations of visco-resistive, single-fluid MHD in toroidal and cylindrical geometries with limited axial periodicity. Both measurements and the simulation with toroidal geometry show stronger inboard flows relative to the outboard flows, which is attributed to the toroidal geometry of the device. In the experiment, the n ¼ 5 component of flow is phase shifted from the reconnection-like flow pattern observed in the single-fluid simulations, possibly due to decoupling of the ion and electron fluids over much of the plasma. Finally, the strength of the helical angular flow shear relative to the critical shear necessary to disrupt nonlinear coupling between tearing modes is calculated around the helical magnetic axis. The shear in the measured flow is on the order of the theoretical critical threshold needed to nonlinearly decouple modes, but the measurement uncertainty in the gradient of the flow is large.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Assessment of Model Hub Height Wind Speed Performance Using DOE Lidar Buoy Data

As the offshore wind industry continues to gain momentum both in the U.S. and globally, reanalysis products remain essential tools for the lifecycle of a wind farm. From project siting and production estimates to determining construction and maintenance windows, reanalysis products play an important role. These regional and global long-term models also serve as inputs to higher resolution models. With such wide-ranging impacts on the offshore wind industry, the need for validation of reanalysis products is strong, particularly at typical turbine hub heights. Pacific Northwest National Laboratory (PNNL) operates two AXYS WindSentinelTM buoys for the U.S. Department of Energy (DOE) in order to collect meteorological and oceanographic data, including hub height wind speed, in areas of interest for offshore wind development. This work uses DOE buoy observations off the coasts of New Jersey and Virginia to evaluate the performance of Modern Era Retrospective Analysis for Research and Applications-2 (MERRA-2), the North American Regional Reanalysis (NARR), the European Center for Medium-Range Weather Forecasts (ECMWF) Reanalysis 5th generation (ERA5), and the National Oceanic and Atmospheric Administration (NOAA) assimilation system Rapid Refresh (RAP). Biases and degrees of correspondence are determined for each reanalysis product in order to provide insights on the performance and uncertainty for long-term wind resource characterization. With the baseline metrics determined, this work proceeds with an investigation into the sources of large deviations between the modeled and observed hub height wind speeds during the East Coast lidar buoy deployments. Consistent sources of reanalysis model error include stably stratified flow conditions, high wind shear, flow parallel to the coastline, and tropical and winter storms. For the near shore New Jersey location, offshore winds (i.e. blowing from land to water) tend to be associated with model overestimation of observed hub height wind speed, while winds parallel to the coastline are correlated with model underestimation of wind speed. Similarly, near the Virginia buoy, during southerly winds (parallel to the coastline) models underestimate the observed wind speed due to coastal upwelling. Reanalysis model wind speeds underestimate the hub height winds with increasing stable atmospheric stratification. Reanalysis model biases are seen when the observed wind speeds approach or exceed typical turbine cut-out speeds (around 20 m s-1). Large model errors during winter and tropical storm events, high wind speed events, and ramp events are often the result of wind speed timing and magnitude offsets. Overall, ERA5 provides the most successful representation of observed offshore hub height wind speeds at the U.S. East Coast buoy locations and is therefore best suited as an input boundary condition for model case studies. ERA5 performs well for a variety of atmospheric phenomena, including storms and sea breezes, however RAP, the next most successful reanalysis model, performs best in capturing the frequency of ramp events.

17 WIND ENERGY↗

Higher-level simulations of turbulent flows

The fundamentals of large eddy simulation are considered and the approaches to it are compared. Subgrid scale models and the development of models for the Reynolds-averaged equations are discussed as well as the use of full simulation in testing these models. Numerical methods used in simulating large eddies, the simulation of homogeneous flows, and results from full and large scale eddy simulations of such flows are examined. Free shear flows are considered with emphasis on the mixing layer and wake simulation. Wall-bounded flow (channel flow) and recent work on the boundary layer are also discussed. Applications of large eddy simulation and full simulation in meteorological and environmental contexts are included along with a look at the direction in which work is proceeding and what can be expected from higher-level simulation in the future.

Ferziger, J. H.↗