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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 469 records · Page 26

All sources NASA NTRS DOE OSTI

DOE OSTI2025

A High-Order Discontinuous Galerkin Spectral Element Method for Compressible Reacting Flows

High-order methods have recently been shown to be an effective tool for high-fidelity flow computations like direct numerical simulations and large eddy simulations due to their strong balance between accuracy and computational cost. In this work, a high-order discontinuous Galerkin spectral element method (DGSEM) is developed to solve the chemically reactive Euler equations encountered in high-speed combustion. To handle the disparate length and time scales associated with these equations, we develop a novel method which combines the spectral accuracy of the SEM with the flexibility of DG approach. Thus, the framework is well suited to capture turbulence in smooth regions of the flow, while maintaining numerical stability in the presence of shocks. The numerical method is implemented within the spectral element solver Nek5000. Validation cases are conducted for both non-reactive and reactive discontinuous flows to demonstrate the solver capability. In particular, canonical one-dimensional and two-dimensional detonation simulations are performed and the high-order numerical results are validated against available literature data.

computational fluid dynamics (CFD)↗

NASA NTRS1984

A review and analysis of boundary layer transition data for turbine application

A symposium on transition in turbines was held at the NASA Lewis Research Center. One recommendation of the working groups was the collection of existing transition data to provide standard cases against which models could be tested. A number of data sets from the open literature that include heat transfer data in apparently transitional boundary layers, with particular application to the turbine environment, were reviewed and analyzed to extract transition information from the heat transfer data. The data sets reviewed cover a wide range of flow conditions, from low speed, flat plate tests to full scale turbine airfoils operating at simulated turbine engine conditions. The results indicate that free stream turbulence and pressure gradient have strong, and opposite, effects on the location of the start of transition and on the length of the transition zone.

Gaugler, R. E.↗

NASA NTRS1985

A review and analysis of boundary layer transition data for turbine application

Gaugler, R. E.↗

NASA NTRS1992

Theoretical and computational study of rotating isotropic turbulence

The various flow regimes that occur in rotating isotropic turbulence and their associated turbulence structure are studied. Direct numerical simulations are conducted for moderate and rapid rotation rates and comparisons are made with the predictions of generalized Eddy Damped Quasi-Normal Markovian (EDQNM) approximations. It is shown that at high rotation rates the nonlinear transfer terms remain small and that the development of the spectrum is through pure viscous decay. This shows that the effect of rapid rotation on the turbulence kinetic energy is through the shutting off of the production term in the dissipation rate equation. At moderate and relatively rapid rates, the rotation causes a discernible anisotropy to develop in the integral length scale with a mild trend toward a two-dimensionalization of the flow. As the flow decays, the Rossby number decreases, leading again to the shutting off of the nonlinear transfer terms. For extremely rapid rotation rates, the anisotropies in the integral length scales are small and direct numerical simulation results are in good agreement with rapid distortion theory.

Mansour, N. N.↗

NASA NTRS1995

Modeling the decay of energy containing eddies: A source of solar wind heating

To understand the solar wind heating and acceleration mechanisms one needs to understand the decay of energy containing eddies. With this goal in mind, attempts have been made to extend the fluid dynamic phenomenology of large scale quasi-equilibrium to the case of magnetohydrodynamics. Matthaeus et al. have proposed a model for the inhomogeneous transport and decay of five mean variables, namely, two mean square Elsasser variables z(exp 2) (sub +/-) their correlation lengths, and the difference between the kinetic and magnetic energies. We test the validity of this model in the simplified case of homogeneous turbulence simulated in a periodic box. We propose a class of models and show that they may fit the simulation satisfactorily. Analytic solutions of this class of model reveal their inherent properties and demonstrate the difficulties associated with finite cross helicity. It is noted that adjustments are required to make the simplest models, which are based upon isotropic turbulence, scale properly with respect to the strength of the mean magnetic field. This can be interpreted as due to anisotropic turbulence, which can be modelled by simple parameterization in the phenomenology.

Hossain, M.↗

DOE OSTI2024

Sheath constraints on turbulent magnetised plasmas

Abstract A solid target in contact with a plasma charges (negatively) to reflect the more mobile species (electrons) and thus keep the bulk plasma quasineutral. To shield the bulk plasma from the charged target, there is an oppositely (positively) charged sheath with a sharp electrostatic potential variation on the Debye length scale λ D . In magnetised plasmas where the magnetic field is inclined at an oblique angle α with the target, some of the sheath potential variation occurs also on the ion sound gyroradius length scale ρ S cos ⁡ α , caused by finite ion gyro-orbit distortion and losses. We consider a collisionless and steady-state magnetised plasma sheath whose thickness l ms ∼ max ( λ D , ρ S cos ⁡ α ) is smaller than the characteristic length scale L of spatial fluctuations in the bulk plasma, such that the limit l ms / L → 0 is appropriate. Spatial structures are assumed to be magnetic field-aligned. In the case of small magnetic field angle α ∼ δ ≡ ρ S / L ≪ 1 , electric fields tangential to the target transport ions towards the target via E × B drifts at a rate comparable to the one from parallel streaming. A generalised form of the kinetic Bohm–Chodura criterion at the sheath entrance is derived by requiring that the sheath electric field have a monotonic spatial decay far from the target. The criterion depends on tangential gradients of potential and ion distribution function, with additional nontrivial conditions.

kinetic Bohm–Chodura condition↗

DOE OSTI2024

LLNL FESP Theory Highlights: March 2024

Predictions of heat load widths λ q based on particle orbits alone are very pessimistic. This paper shows that pedestal peeling-ballooning (P-B) MHD turbulence broadens the stable SOL by the transport, or spreading, of fluctuation energy from the pedestal. λ q is seen to increase with Γ ε , the fluctuation energy density flux. We elucidate the fundamental physics of the spreading process. Γ ε increases with pressure fluctuation correlation length. P-B turbulence is seen to be especially effective at spreading, on account of its large effective mixing length. Spreading is shown to be a multiscale process, which is enhanced by the synergy of large and small-scale modes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS2007

Effect of Temperature on Jet Velocity Spectra

Statistical jet noise prediction codes that accurately predict spectral directivity for both cold and hot jets are highly sought both in industry and academia. Their formulation, whether based upon manipulations of the Navier-Stokes equations or upon heuristic arguments, require substantial experimental observation of jet turbulence statistics. Unfortunately, the statistics of most interest involve the space-time correlation of flow quantities, especially velocity. Until the last 10 years, all turbulence statistics were made with single-point probes, such as hotwires or laser Doppler anemometry. Particle image velocimetry (PIV) brought many new insights with its ability to measure velocity fields over large regions of jets simultaneously; however, it could not measure velocity at rates higher than a few fields per second, making it unsuitable for obtaining temporal spectra and correlations. The development of time-resolved PIV, herein called TR-PIV, has removed this limitation, enabling measurement of velocity fields at high resolution in both space and time. In this paper, ground-breaking results from the application of TR-PIV to single-flow hot jets are used to explore the impact of heat on turbulent statistics of interest to jet noise models. First, a brief summary of validation studies is reported, undertaken to show that the new technique produces the same trusted results as hotwire at cold, low-speed jets. Second, velocity spectra from cold and hot jets are compared to see the effect of heat on the spectra. It is seen that heated jets possess 10 percent more turbulence intensity compared to the unheated jets with the same velocity. The spectral shapes, when normalized using Strouhal scaling, are insensitive to temperature if the stream-wise location is normalized relative to the potential core length. Similarly, second order velocity correlations, of interest in modeling of jet noise sources, are also insensitive to temperature as well.

Bridges, James E.↗

NASA NTRS2007

Offset Stream Technology Test-Summary of Results

Brown, Clifford A.↗

NASA NTRS2014

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.↗

NASA NTRS1990

Spanwise measurements of vertical components of atmospheric turbulence

Correlation and spectrum magnitude estimates are computed for vertical gust velocity measurements at the nose and wing tips of a NASA B-57B aircraft for six level flight, low speed and low altitude runs and are compared with those of the von Karman atmospheric turbulence model extended for spanwise relationships. The distance between the wing tips was 62.6 ft. Airspeeds ranged from about 330 to 400 ft/sec, heights above the ground ranged from near ground level to about 5250 ft. and gust velocity standard deviations ranged from 4.10 to 8.86 ft/sec. Integral scale lengths, determined by matching measured autocorrelation estimates with those of the model, ranged from 410 to 2050 ft. Digital signals derived from piezoelectric sensors provided continuous pressure and airspeed measurements. Some directional acceleration sensitivity of the sensors was eliminated by sensor orientation, and their performance was spectrally verified for the higher frequencies with supplemental onboard piezoresistive sensors. The model appeared to satisfactorily predict the trends of the measured cross-correlations and cross-spectrum magnitudes, particularly between the nose and wing tips. However, the measured magnitude estimates of the cross-spectra between the wing tips exceeded the predicted levels at the higher frequencies. Causes for the additional power across the wing tips were investigated. Vertical gust velocity components evaluated along and lateral to the flight path implied that the frozen-turbulence-field assumption is a suitable approximation.

Sleeper, Robert K.↗

DOE OSTI2021

Cloud droplet diffusional growth in homogeneous isotropic turbulence: bin microphysics versus Lagrangian super-droplet simulations

The increase in the spectral width of an initially monodisperse population of cloud droplets in homogeneous isotropic turbulence is investigated by applying a finite-difference fluid flow model combined with either Eulerian bin microphysics or a Lagrangian particle-based scheme. The turbulence is forced applying a variant of the so-called linear forcing method that maintains the mean turbulent kinetic energy (TKE) and the TKE partitioning between velocity components. The latter is important for maintaining the quasi-steady forcing of the supersaturation fluctuations that drive the increase in the spectral width. We apply a large computational domain (64 3 m 3 ), one of the domains considered in Thomas et al. (2020). The simulations apply 1 m grid length and are in the spirit of the implicit large eddy simulation (ILES), that is, with small-scale dissipation provided by the model numerics. This is in contrast to the scaled-up direct numerical simulation (DNS) applied in Thomas et al. (2020). Two TKE intensities and three different droplet concentrations are considered. Analytic solutions derived in Sardina et al. (2015), valid for the case when the turbulence integral timescale is much larger than the droplet phase relaxation timescale, are used to guide the comparison between the two microphysics simulation techniques. The Lagrangian approach reproduces the scalings relatively well. Representing the spectral width increase in time is more challenging for the bin microphysics because appropriately high resolution in the bin space is needed. The bin width of 0.5 µm is only sufficient for the lowest droplet concentration (26 cm -3 ). For the highest droplet concentration (650 cm -3 ), an order of magnitude smaller bin size is barely sufficient. The scalings are not expected to be valid for the lowest droplet concentration and the high-TKE case, and the two microphysics schemes represent similar departures. Finally, because the fluid flow is the same for all simulations featuring either low or high TKE, one can compare point-by-point simulation results. Such a comparison shows very close temperature and water vapor point-by-point values across the computational domain and larger differences between simulated mean droplet radii and spectral width. The latter are explained by fundamental differences in the two simulation methodologies, numerical diffusion in the Eulerian bin approach and a relatively small number of Lagrangian particles that are used in the particle-based microphysics.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

The role of ion-scale micro-turbulence in pedestal width of the DIII-D wide-pedestal QH mode

The low-edge rotation, intrinsically ELM-free, and improved confinement wide-pedestal quiescent H-mode (QH-mode), discovered in DIII-D tokamak, has pedestal widths exceeding the EPED-kinetic-ballooning mode (KBM) model scaling typically by at least 25%. Ion-scale ($k_yp_s$ < 1) microturbulence and its role in setting the pedestal structure is investigated using the radially local δ$f$ gyrokinetic code CGYRO. The electromagnetic trapped electron mode (TEM) is unstable at the pedestal top, while plasma beta (β$_e$) is ∼60% below the KBM onset threshold and the electron temperature gradient mode is found to be unstable in the peak gradient region. Nonlinear simulation reveals that the ion-scale turbulence could produce electron energy flux consistent with the flux inferred from power balance at the pedestal top, with a reasonable variation of the local shearing rate; and the local neoclassical transport from NEO is dominant over the simulated turbulent transport in the ion energy flux channel. The simulated ion-scale turbulence produces much lower electron energy flux than inferred from experiment in the pedestal peak gradient region. A correction to the EPED-KBM pedestal width scaling is obtained based on the two-dimensional scan of pedestal top plasma beta (β$_e$) and normalized electron density and temperature scale lengths,${a}$/${L_n}_e$, ${a}$/${L_T}_n$using CGYRO linear simulations. Mode transitions among TEM, micro-tearing mode, ion-temperature gradient mode and KBM, are observed in the 2D scan at the pedestal top. A fixed normalized growth rate for these drift-type modes is taken to determine the pedestal width scaling, which shows good consistency with the QH experimental database on pedestal heights and widths. The onset of KBM instabilities and the local E × B shear suppression criterion set the lower and upper limit for the pedestal width of standard QH-mode, wide-pedestal QH-mode and type-I ELMy H mode. A potentially higher and wider pedestal is expected from the new scaling of pedestal width. This work presents an improved understanding of the ion-scale micro-turbulence of wide-pedestal QH-mode and sheds light on a promising scenario for future reactors, including ITER and beyond.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Beta-dependent Properties of Solar Wind Current Sheets at 1 au

We utilize plasma turbulence in interplanetary coronal mass ejections as a natural laboratory to analyze more than 4600 current sheets observed over a broad range of electron and proton beta values, 10 −2 ≲ β e ≲ 10 and 10 −3 ≲ β p ≲ 10. The analysis shows that the shear angle Δθ and the thickness λ CS of the current sheets depend on both electron and proton beta. Specifically, the dependencies on electron beta are Δθ ≈ 12.6° $β^{0.41}_{e}$ and, λ CS /λ p ≈ 1.9 $β^{0.26}_{e}$ where λ p is the local proton inertial length. We argue that the beta dependence of the shear angle is an intrinsic feature of solar wind turbulence arising from the natural correlation between turbulence intensity and plasma beta. In contrast, the beta dependence of current sheet scales is likely a general characteristic of plasma turbulence. According to recent theory, current sheets formed in turbulence are disrupted by the electron tearing instability once their thickness is smaller than a critical scale. We demonstrate that normalizing the thickness of current sheets by this critical scale eliminates the dependence on beta across the entire range of plasma beta. This observation not only supports the theory suggesting that current sheets mediate the transition from the inertial to the kinetic cascade through the electron tearing instability, but also demonstrates the plausibility of similar processes in astrophysical plasma. The presented results also indicate that the theory may be applicable beyond its formal limits, m e /m i ≪ β ≪ 1.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1977

A study of key features of random atmospheric disturbance models for the approach flight phase

An analysis and brief simulator experiment were performed to identify and classify important features of random turbulence for the landing approach flight phase. The analysis of various wind models was carried out within the context of the longitudinal closed-loop pilot/vehicle system. The analysis demonstrated the relative importance of atmospheric disturbance scale lengths, horizontal versus vertical gust components, decreasing altitude, and spectral forms of disturbances versus the pilot/vehicle system. Among certain competing wind models, the analysis predicted no significant difference in pilot performance. This was confirmed by a moving base simulator experiment which evaluated the two most extreme models. A number of conclusions were reached: attitude constrained equations do provide a simple but effective approach to describing the closed-loop pilot/vehicle. At low altitudes the horizontal gust component dominates pilot/vehicle performance.

Heffley, R. K.↗

NASA NTRS1994

Turbulent diffusion and concentration of chondrules in the protoplanetary nebula

Turbulence is known to possess structure on many scales. The largest or integral scale L is usually taken to be the largest dimension of the system. Two different turbulent regimes are likely to be of importance. Early stage convective turbulence in a hot nebula probably extends over the entire nebula scale height, and probably has typical eddy frequency comparable to the orbit frequency. At a later stage, midplane shear layer turbulence can be driven by the presence of a differently rotating, densely settled particle layer, with different length and timescales.

Cuzzi, J. N.↗

NASA NTRS1988

Turbulent transport observed just above the Amazon forest

Observations of turbulent heat, moisture, and momentum transport made at two levels, approximately 5 and 10 m above the Amazon rain forest canopy, are presented. Data acquired at 10 Hz included variances and some mixed third moments of vertical velocity, temperature, and humidity. Two features of the data appear to question the displacement height hypothesis: (1) the characteristic dissipation length scale in the near-canopy layer varied between 20 m in stable conditions to approximately 150 m during afternoon convective conditions, generally larger scales than would be expected, and (2) no appreciable differences in dissipation scales was seen at two observed levels. Heat budgets on selected days show that frequent periods with negative heat flux concurrent with continuing positive moisture flux occur in early afternoon, and this is believed to indicate the patchy nature of canopy-atmosphere coupling. Vertical velocity skewness was observed to be negative on three successive days and exhibited a sharp positive gradient. Time series of some of the terms in the turbulence budgets of vertical velocity variance and kinematic heat flux are presented.

Fitzjarrald, David R.↗

NASA NTRS2006

e(sup +/-) Pair Loading and the Origin of the Upstream Magnetic Field in GRB Shocks

We investigate here the effects of plasma instabilities driven by rapid e(sup +/-) pair cascades, which arise in the environment of GRB sources as a result of back-scattering of a seed fraction of their original spectrum. The injection of e(sup +/-) pairs induces strong streaming motions in the ambient medium. One therefore expects the pair-enriched medium ahead of the forward shock to be strongly sheared on length scales comparable to the radiation front thickness. Using three-dimensional particle-in-cell simulations, we show that plasma instabilities driven by these streaming e(sup +/-) pairs are responsible for the excitation of near-equipartition, turbulent magnetic fields. Our results reveal the importance of the electromagnetic filamentation instability in ensuring an effective coupling between e(sup +/-) pairs and ions, and may help explain the origin of large upstream fields in GRB shocks.

Ramirez-Ruiz, Enrico↗