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 361 records · Page 20

Unveiling the Significance of Correlations in K-Space and Configuration Space for Drift Wave Turbulence in Tokamaks

Turbulence and transport phenomena play a crucial role in the confinement and stability of tokamak plasmas. Turbulent fluctuations in certain physical quantities, such as density or temperature fluctuations, can have a wide range of spatial scales, and understanding their correlation length is important for predicting and controlling the behavior of the plasma. The correlation length in the radial direction is identified as the critical length in real space. The dynamics in real space are of significant interest because transport in configuration space is primarily focused on them. When investigating transport caused by the E×B drift, the correlation length in real space represents the size of E×B whirls. It was numerically discovered that in drift wave turbulence, this length is inversely proportional to the normalized mode number of the fastest growing mode relative to the drift frequency. Considerable time was required before a proper analytical derivation of this condition was accomplished. Therefore, a connection has been established between phenomena occurring in real space and those occurring in k-space. Although accompanied by a turbulent spectrum in k-space with a substantial width, transport in real space is uniquely determined by the correlation length, allowing for accurate transport calculations through the dynamics of a single mode. Naturally, the dynamics are subject to nonlinear effects, with resonance broadening in frequency being the most significant nonlinear effect. Thus, mode number space is once again involved. Resonance broadening leads to the detuning of waves from particles, permitting a fluid treatment. It should be emphasized that the consideration here involves the total electric field, including the induction part, which becomes particularly important at higher beta plasmas.

Weiland, Jan (ORCID:0000000253496025)↗

Evaluation of Inflow Turbulence Methods in Large-Eddy Simulations of a Supersonic Boundary Layer

In the context of Large- Eddy Simulations (LES), Boundary-layer inflow turbulence is simulated using both the Synthetic Eddy Model (SEM) and Digital Filtering (DF). The effects of the projection error are investigated. The effect of the prescribed length scales on the adjustment region was found to be negligible for length scales less than one-tenth of the boundary-layer thickness. While it was conjectured that one method of the two might be more robust than the other, our results show that both the Digital Filtering Method and the Synthetic Eddy Method accurately replicate the boundary layer while successfully accounting for inflow turbulence.

Synthetic Eddy↗

Vertical Scales of Turbulence at the Mt. Wilson Observatory

The vertical scales of turbulence at the Mt. Wilson observatory are inferred from data from the University of California at Berkeley Infrared Spatial Interferometer (ISI), by modeling path length fluctuations observed in the interferometric paths to celestial objects and those in instrumental ground- based paths. The correlations between the stellar and ground-based path length fluctuations and the temporal statistic of those fluctuations are modeled on various time scales to constrain the vertical scales.

Mt. Wilson Observatory↗

Probing eddy size and its effective mixing length in stably stratified roughness sublayer flows

Stably stratified roughness sublayer flows are ubiquitous yet remain difficult to represent in models and to interpret using field experiments. Here, continuous high-frequency potential temperature profiles from the forest floor up to 6.5 times the canopy height observed with distributed temperature sensing (DTS) are used to link eddy topology to roughness sublayer stability correction functions and coupling between air layers within and above the canopy. The experiments are conducted at two forest stands classified as hydrodynamically sparse and dense. Near-continuous profiles of eddy sizes (length scales) and effective mixing lengths for heat are derived from the observed profiles using a novel conditional sampling approach. The approach utilizes potential temperature isoline fluctuations from a statically stable background state. The transport of potential temperature by an observed eddy is assumed to be conserved (adiabatic movement) and we assume that irreversible heat exchange between the eddy and the surrounding background occurs along the (vertical) periphery of the eddy. This assumption is analogous to Prandtl's mixing-length concept, where momentum is transported rapidly vertically and then equilibrated with the local mean velocity gradient. A distinct dependence of the derived length scales on background stratification, height above ground, and canopy characteristics emerges from the observed profiles. Implications of these findings for (1) the failure of Monin–Obukhov similarity in the roughness sublayer and (2) above-canopy flow coupling to the forest floor are examined. The findings have practical applications in terms of analysing similar DTS data sets with the proposed approach, modelling roughness sublayer flows, and interpreting nocturnal eddy covariance measurements above tall forested canopies.

54 ENVIRONMENTAL SCIENCES↗

Neoclassical transport and profile prediction in transport barriers

Strong gradient regions in tokamaks, such as the pedestal or internal transport barriers, are regions of reduced turbulence where neoclassical transport can play a dominant role. However, standard neoclassical transport theory assumes that the gradient length scales of density, temperature, and potential are of the order of the system size. In the pedestal, gradient length scales are much shorter and are measured to be of the order of the ion poloidal gyroradius. We present an extension of neoclassical theory that is applicable in transport barriers of large aspect ratio tokamaks. We show that particle and momentum transport are connected in such a way that a source of parallel momentum can drive a significant neoclassical ion particle flux. In strong gradient regions, density, electric potential, mean parallel flow, and ion temperature are shown to no longer be flux functions. Instead, they have a small but important poloidally varying piece that modifies the transport equations to lowest order. This introduces a nonlinearity in the transport problem through the coupling with quasineutrality that yields multiple co-existing solutions when solving for the plasma profiles. The different solutions could be connected to low and high transport states and jumps between solutions could provide a new neoclassical explanation for H-L back-transitions.

Trinczek, Silvia [Princeton Plasma Physics Laborat↗

Modeling of Turbulent Flow in Electromagnetically Levitated Metal Droplets

This article details an effort to improve the understanding and prediction of turbulent flow inside a droplet of molten metal levitated in an electromagnetic field. It is shown that the flow field in a test case, a nickel droplet levitated under microgravity conditions, is in the transitional regime between laminar and turbulent flow. Past research efforts have used laminar, enhanced viscosity, and k-epsilon turbulence models to describe the flow. The method highlighted in our study is the renormalization group (RNG) algorithm. We show that an accurate description of the turbulent eddy viscosity is critical in order to obtain realistic velocity fields, and that the turbulent eddy viscosity cannot be uniform in levitated droplets. The RNG method does not impose isotropic length or time scales on the flow field, thus allowing such nonuniform features to be captured. A number of other materials processing applications exhibit similarly complex flow characteristics, such as highly recirculating, transitional, and free surface flows, for which this modeling approach may prove useful.

Berry, S.↗

Transport of Solar Wind Fluctuations: A Two-Component Model

We present a new model for the transport of solar wind fluctuations which treats them as two interacting incompressible components: quasi-two-dimensional turbulence and a wave-like piece. Quantities solved for include the energy, cross helicity, and characteristic transverse length scale of each component, plus the proton temperature. The development of the model is outlined and numerical solutions are compared with spacecraft observations. Compared to previous single-component models, this new model incorporates a more physically realistic treatment of fluctuations induced by pickup ions and yields improved agreement with observed values of the correlation length, while maintaining good observational accord with the energy, cross helicity, and temperature.

Oughton, S.↗

Calculation of the Heat Transfer Coefficient in the Outer Body for a Rotational Detonation

Unsteady heat transfer characterization on the combustion surfaces of Rotational Detonation Engines (RDE) is not well understood. It is generally thought that the complex nature of the unsteady, reacting, compressible fluid flow inside the combustion anulus of the RDE causes the convective heat transfer coefficient to be significantly higher than it is in other applications. Empirical models that have been used to analyze this strictly apply to steady flows where dimensionless groups can be employed. The reacting flow fields that are characteristic of RDEs are inherently compressible, three dimensional, unsteady, and turbulent, having properties that change by orders of magnitude throughout the flow field. They will therefore contain multiple length scales and time scales operating everywhere in the flow during all times. Because of this it is not likely that the RDE flow fields will lend themselves to explanation using simple dimensionless parameters. The dimensionless groups have meaning only in situations where length scales and time scales are singular and well defined. In spite of this it may be possible to get a relatively good idea of what the convection heat transfer coefficient is. In this work a numerical study is performed where the inside wall surface temperature distribution in the RDE outer body is systematically changed over a given range that would be characteristic of the start-up flows inside an RDE. For each case, temperature distributions inside the outer containment wall of the RDE was calculated and compared with experimental data. The closet match can then be used to directly calculate the convection heat transfer coefficient on the inside surface of the RDE.

VanOsdol, John↗

PARTICLE RESOLVED DNS STUDY OF TURBULENCE EFFECTS ON HYPORHEIC MIXING IN RANDOMLY PACKED SEDIMENT BEDS

Pore-resolved direct numerical simulations (DNS) are used to investigate the interactions between stream-water flow turbulence and groundwater flow through a porous sediment bed in the hyporheic zone. Two permeability Reynolds numbers (2.56 and 5.17), representative of aquatic systems and representing ratio of permeability to viscous length scales, were simulated to understand its influence on the momentum exchange at the sediment-water interface (SWI). A doubleaveraging methodology is used to compute the Reynolds stresses, form-induced stresses, and pressure fluctuations. It is observed that both shear layer and turbulent shear stress penetration increases with ReK. Reynolds and form-induced bed-normal stresses increase with ReK. The peak values of the form-induced stresses for the lower (2.56) and higher (5.17) ReK happen within the top layer of the sediment bed. The sum of turbulent and form-induced pressure fluctuations, analyzed at their respective zero-displacement planes, are statistically similar and can be well approximated by a t location-scale distribution fit providing with a model that could potentially be used to impose boundary conditions at the SWI in reach scale simulations.

Karra, Shashank↗

Estimating and Evaluating Roughness Length and Displacement Height in Heterogeneous Urban Environments

The roughness length (z 0 ) and displacement height (z d ) are essential surface-layer parameters in numerical models (e.g., weather, climate, wall-modeled LES, etc.). This work evaluates the consistency of z 0 and z d estimates from morphometric and anemometric methods using data from two eddy-covariance flux towers (AmeriFlux US-INg and US-INc) in Indianapolis, IN. Results show inconsistencies in estimated z 0 and z d values depending on the chosen method. The two evaluated anemometric methods estimate non-physical values of z d when compared to roughness elements surrounding both towers. Additionally, predictions of mean wind speed using surface-layer similarity theory with morphometric estimates exhibit a bias during near-neutral and stable conditions relative to observations. The overestimation of mean wind speed by surface layer similarity theory is consistent with previous observational and modeling studies in urban areas, suggesting that the application of similarity theories to urban environments may have limitations. Differentiation of vegetation from built structures appears to impact morphometric z 0 and z d estimates, particularly where vegetation is abundant; however, it has little impact on correcting biases in the similarity theory. Specifically, we find that existing similarity theories using morphometric estimates underestimate integral velocity and length scales, and the degree of underestimation depends on the stability conditions. Accounting for the degree of anisotropy in surface-layer turbulence helps reduce the biases between similarity theories and observations during unstable conditions, but not in near-neutral cases. Future work is needed to identify the cause of such biases for near-neutral conditions.

Aerodynamic roughness length↗

Scaling of Turbulent Viscosity and Resistivity: Extracting a Scale-dependent Turbulent Magnetic Prandtl Number

Turbulent viscosity νt and resistivity ηt are perhaps the simplest models for turbulent transport of angular momentum and magnetic fields, respectively. The associated turbulent magnetic Prandtl number Pr t ≡ ν t /η t has been well recognized to determine the final magnetic configuration of accretion disks. Here, we present an approach to determining these "effective transport" coefficients acting at different length scales using coarse-graining and recent results on decoupled kinetic and magnetic energy cascades. By analyzing the kinetic and magnetic energy cascades from a suite of high-resolution simulations, we show that our definitions of ν t , η t , and Pr t have power-law scalings in the "decoupled range." We observe that Pr t ≈ 1–2 at the smallest inertial-inductive scales, increasing to ≈5 at the largest scales. Furthermore, based on physical considerations, our analysis suggests that Pr t has to become scale independent and of order unity in the decoupled range at sufficiently high Reynolds numbers (or grid resolution) and that the power-law scaling exponents of velocity and magnetic spectra become equal. In addition to implications for astrophysical systems, the scale-dependent turbulent transport coefficients offer a guide for large-eddy simulation modeling.

79 ASTRONOMY AND ASTROPHYSICS↗

Calculation of free turbulent mixing by interaction approach.

The applicability of Bradshaw's interaction hypothesis to two-dimensional free shear flows was investigated. According to it, flows with velocity extrema may be considered to consist of several interacting layers. The hypothesis leads to a new expression for the shear stress which removes the usual restriction that shear stress vanishes at the velocity extremum. The approach is based on kinetic energy and the length scale equations. The compressible flow equations are simplified by restriction to low Mach numbers, and the range of their applicability is discussed. The empirical functions of the turbulence model are found here to be correlated with the spreading rate of the shear layer. The analysis demonstrates that the interaction hypothesis is a workable concept.

Morel, T.↗

Split-Step Simulations to Assess the Effects of Atmospheric Boundary Layer Turbulence on the Dose Variability of N-Waves and Shaped Booms

The effects of atmospheric boundary layer turbulence on the loudness variability of a sonic boom N-wave and shaped boom are examined with split-step simulations. The shaped boom is representative of a design iteration of the NASA X-59 aircraft. Inhomogeneous atmospheric boundary layer turbulence is generated in the computational domain by a Fourier synthesis method. The N-wave and shaped boom are propagated through turbulent fields representing eight different convection levels measured at the NASA Kennedy Space Center and the NASA Armstrong Flight Research Center. Probability density functions of the formation of caustic regions along the propagation direction are computed from the N-wave results, and a parameter to collapse the caustic PDFs that accounts for both fluctuation intensities and length scales is proposed. Statistical results concerning loudness metric variability are presented, and the standard deviations of several metrics are shown to collapse across different convection levels of turbulence for small nondimensional propagation distances. The loudness metric distributions are observed to be well approximated by a normal distribution for a given range of propagation distances, and become increasingly skewed as distance increases. A model function for the dose variability is proposed, and the function parameters are found to be related to the convection level of the turbulence. The model for the dose variability distribution is compared to simulation data that were not used to find the regression parameters of the model. At several nondimensional propagation distances, agreement is observed between the model and the simulation data. These results indicate that the model may be suitable for providing quick estimates of noise dose variability in the primary carpet region across a wide range of atmospheric boundary layer conditions.

Sonic boom↗

Low-subsonic aerodynamic characteristics of a shuttle-orbiter configuration designed for reduced length

An investigation has been made in a low-turbulence pressure tunnel to determine the low-subsonic aerodynamic characteristics of a 0.01875-scale model of a potential shuttle orbiter. The design has the rocket engines mounted in fairings on either side of the body on top of the wing. The wing had a leading-edge sweep of 50 and a trailing-edge sweep of minus 4. configurations investigated included engine-mounted twin dorsal tails at various rollout angles, a body-mounted center-line vertical tail, cylindrical and boattailed afterbody, and elevon and rudder at several deflections.

Ware, G. M.↗

Measurements of the imploding plasma sheath in triple-nozzle gas-puff z pinches

Gas-puff z-pinch implosions are characterized by the formation of a dense annular plasma shell, the sheath, that is driven to the axis by magnetic forces and therefore subject to the magneto-Rayleigh–Taylor instability. Here, the conditions within these sheaths are measured on the 1-MA COBRA generator at Cornell University for various gas species and initial fill densities. The gas-puff loads are initialized by a 7 cm diameter triple-nozzle gas valve assembly with concentric outer and inner annular nozzles and a central gas jet. Thomson scattering and laser interferometry provide spatially resolved flow, temperature, and electron density profiles midway through the implosion, while extreme ultraviolet pinhole cameras record the evolution of the plasma column and photoconducting diodes measure x-ray emission. Analysis of the scattering spectra includes a means of discriminating between thermal and non-thermal broadening to test for the presence of hydrodynamic turbulence. Two types of sheath profiles are observed, those with sharp discontinuities at the leading edge and those with smooth gradients. In both cases, non-thermal broadening is generally peaked at the front of the sheath and exhibits a characteristic decay length that roughly scales with the sheath ion mean free path. Here, we demonstrate that this non-thermal broadening term is inconsistent with laminar velocity gradients and is more consistent with dissipative turbulence driven by unstable plasma waves in a collisionless shock. The resulting differences in sheath profile are then set by the sheath ion collisionality in a manner consistent with recent 1D kinetic simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Comet-solar wind interaction - Dynamical length scales and models

ICE magnetometer measurements at Comet Giacobini-Zinner and model simulations of comet-solar wind interactions are analyzed. The magnetometer data reveal the existence of intense hydromagnetic turbulence, a draping of the magnetic field lines to form a magnetotail, a weak shock, and a magnetic barrier region in the magnetosphere. The global models of the comet-solar wind interaction are described. The observed data and models are compared and good correlation is displayed.

Mendis, D. A.↗

Assessment of turbulent models for scramjet flowfields

The behavior of several turbulence models applied to the prediction of scramjet combustor flows is described. These models include the basic two equation model, the multiple dissipation length scale variant of the two equation model, and the algebraic stress model (ASM). Predictions were made of planar backward facing step flows and axisymmetric sudden expansion flows using each of these approaches. The formulation of each of these models are discussed, and the application of the different approaches to supersonic flows is described. A modified version of the ASM is found to provide the best prediction of the planar backward facing step flow in the region near the recirculation zone, while the basic ASM provides the best results downstream of the recirculation. Aspects of the interaction of numerica modeling and turbulences modeling as they affect the assessment of turbulence models are discussed.

Sindir, M. M.↗

Inverse energy cascades in three-dimensional turbulence

Fully three-dimensional magnetohydrodynamic (MHD) turbulence at large kinetic and low magnetic Reynolds numbers is considered in the presence of a strong uniform magnetic field. It is shown by numerical simulation of a model of MHD that the energy inverse cascades to longer length scales when the interaction parameter is large. While the steady-state dynamics of the driven problem is three-dimensional in character, the behavior has resemblance to two-dimensional hydrodynamics. These results have implications in turbulence theory, MHD power generator, planetary dynamos, and fusion reactor blanket design.

Hossain, Murshed↗