Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “weak turbulence”

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 433 records · Page 24

A preliminary investigation of large-scale organized motions in a supersonic turbulent boundary layer on a curved surface

An experimental investigation concentrating on the change in the large-scale organized motions of a supersonic turbulent boundary layer subjected to a short region of concave surface curvature is presented. Two flows were studied, a relatively weak and a relatively strong perturbation. In both flows the freestream Mach number was 2.89 and the Reynolds number based on momentum thickness was 81,900. Hot-wire anemometry and schlieren photography provided the principle tools for this investigation. Hot-wire data were analyzed using spectral and correlation methods and conditional sampling techniques. Changes in structure angle and results of conditionally sampled velocity signals are discussed. The strong perturbation is found to increase the structure angle, and in both flows the length scales of the organized motions are altered.

Donovan, J. F.↗

Acoustic and aerodynamic study of a pusher-propeller aircraft model

An aerodynamic and acoustic study was made of a pusher-propeller aircraft model in the NASA-Ames 7 x 10 ft Wind Tunnel. The test section was changed to operate as an open jet. The 591 mm diameter unswept propeller was operated alone and in the wake of three empennages: an I tail, Y tail, and a V tail. The radiated noise and detailed wake properties were measured. Results indicate that the unsteady blade loading caused by the blade interactions with the wake mean velocity distribution had a strong effect on the harmonics of blade passage noise. The blade passage harmonics above the first were substantially increased in all horizontal directions by the empennage/propeller interaction. Directivity in the plane of the propeller was maximum perpendicular to the blade surface. Increasing the tail loading caused the propeller harmonics to increase 3 to 5 dB for an empennage/propeller spacing of 0.38 mean empennage chords. The interaction noise became weak as empennage propeller spacing was increased beyond 1.0 mean empennage chord lengths. Unlike the mean wake deficit, the wake turbulence had only a small effect on the propeller noise, that effect being a small increase in the broadband noise.

Soderman, Paul T.↗

Magnetic structure of slow shocks and the associated ion dissipation

The kinetic structure of slow magnetosonic shocks obtained from 1- and 2-D hybrid (fluid electrons, particle ions) simulations is discussed. The magnetic structure of the shock varies considerably with the Mach number. At low Mach numbers, the shock may be laminar with no associated waves, whereas, at high Mach numbers the shock may be unsteady with Alfvenic turbulence in the upstream and the downstream regions. The nature of ion dissipation also varies with the Mach number. At weak laminar shocks, ion dissipation consists of adiabatic heating and reflection of a small number of ions back into the upstream. At stronger shocks, particle scattering by Alfven/ion-cyclotron waves in the upstream contributes to the heating which continues well into the downstream region.

Omidi, N.↗

Background Turbulence Characterization of the Curved Duct Test Rig

Interest in applying acoustic treatment to nontraditional locations of turbofan engines has led to recent testing of acoustically treated airfoils in the NASA Langley Curved Duct Test Rig. Current testing focuses only on sound attenuation by the sample when exposed to an incident sound field driven by a loudspeaker array. However, there may be much to gain from an aeroacoustic study of the interaction of turbulence with treated airfoils, especially in the case of acoustically treated outlet guide vane designs. Thus, the intent of this study is to take the first step in assessing the aeroacoustic quality of the NASA Langley Curved Duct Test Rig (typically only used for grazing flow studies) by examining the flow quality. The background turbulence levels are measured using hot-wire anemometry just upstream of the test section. Measurements indicate that the turbulence intensity is less than 0.3% in the core region of the flow. The autospectral densities show cavity tones and weak vortex shedding present in the core. The vortex shedding is determined to be from the upstream total pressure probe used to determine flow speed. The cavity tones are found to be the fundamental and harmonic frequencies of a porous ceramic tubular acoustic liner sample that separates the acoustic drivers from the flow region. Recommendations for improving the aeroacoustic quality of the tunnel are provided, such as removing the upstream probe and redesigning the ceramic tubular liner. Future acoustic characterization of the background levels is also recommended to further understand the feasibility of aeroacoustic studies in CDTR.

Wind Tunnel Characterization↗

Simulation and Modeling of Hypersonic Turbulent Boundary Layers Subject to Favorable Pressure Gradients due to Streamline Curvature

Direct numerical simulations (DNS) of favorable-pressure-gradient turbulent boundary layers are presented for a nominal freestream Mach number of 5, with the objective of assessing the limitations of the currently available Reynolds-averaged Navier-Stokes (RANS) models. The favorable pressure gradient is induced by the streamwise curvature of the two-dimensional, planar, convex measurement surface used during experiments at the Texas A&M University. The DNS data shows good comparison with the measured velocity profiles, strain rates, and some, but not all, of the Reynolds-stress components. The discrepancies between the predicted and the measured wall-normal as well as shear stress components are primarily attributed to the lower than actual values inferred from typical PIV measurements of turbulent boundary layers. The DNS data shows a zero or slightly negative Reynolds shear stress in the outer part of the boundary-layer, which is indicative of the decaying turbulent motion under a strong favorable pressure gradient. The DNS data is also compared with the results of RANS computations based on commonly used zero, one, and two equation eddy-viscosity models. The RANS models yield reasonable comparisons with the DNS-based skin friction under zero and weak pressure gradients, but significant discrepancies under a strong pressure gradient. The k-w SST model provided the best overall predictions of skin friction, except in the region where the flow transitions from a favorable to an adverse pressure gradient. While the RANS models examined herein also give good predictions of the Reynolds shear stress under a sufficiently weak pressure gradient, none of those models are able to appropriately capture the reduction in the Reynolds stresses when the flow was subjected to a strong pressure gradient. An a priori assessment of the turbulent heat-flux prediction based on the assumption of a constant turbulent Prandtl number with the DNS data shows that while the constant turbulent Prandtl number model is effective in predicting the wall-normal component of turbulent heat flux, it does not capture the turbulent heat transfer in the streamwise direction for all the pressure gradient cases. The failure of the constant turbulent Prandtl number model highlights a requirement for more advanced models of the turbulent heat flux.

turbulent boundary layers↗

The Evolution of Modulated Wavetrains Into Turbulent Spots

Experiment are being carried out to study the process by which th almost periodic disturbance waves generated naturally by the freestream evolve into turbulence. The boundary layer on a flat plate has been used for this study. The novelty of the approach is in the form of artificial excitation that is used. In this work the flow is excited artificially by deterministic white noise. The weak T-S wave created develops down stream, becomes nonlinear and blows up locally onto a highly distorted flow. These large local distortions of the mean flow allow very high frequency disturbances to grow and form into small turbulent spots. The spots arise from the excitation, and if the same noise sequence is repeated a spot will form at the same position and time instant relative to the excitation.

Gaster, M.↗

Simulations of a Turbulent Flow Subjected to Favorable and Adverse Pressure Gradients

This paper reports the results from a direct numerical simulation of an initially turbulent boundary layer passing over a wall-mounted “speed bump” geometry. The speed bump, represented in the form of a Gaussian distribution profile, generates a favorable pressure gradient region over the upstream half of the geometry, followed by an adverse pressure gradient over the downstream half. The boundary layer approaching the bump undergoes strong acceleration in the favorable pressure gradient region before experiencing incipient or very weak separation within the adverse pressure gradient region. These types of flows have proven to be particularly challenging to predict using lower-fidelity simulation tools based on various turbulence modeling approaches and warrant the use of the highest-fidelity simulation techniques. Simulation results are utilized to examine the key phenomena present in the flowfield, such as relaminarization/stabilization in the strong acceleration region succeeded by retransition to turbulence near the onset of adverse pressure gradient, incipient/weak separation, and development of internal layers where the sense of streamwise pressure gradient changes at the foot, apex and tail of the bump. The present direct numerical simulation is performed using a flow solver developed exclusively for graphics processing units, which is found to provide a significant speedup compared to an earlier solver optimized for central processing unit architectures.

Ali Uzun↗

Use, Assessment, and Improvement of the Loci-CHEM CFD Code for Simulation of Combustion in a Single Element GO2/GH2 Injector and Chamber

This document is a viewgraph presentation of a paper that documents a continuing effort at Marshall Space Flight Center (MSFC) to use, assess, and continually improve CFD codes to the point of material utility in the design of rocket engine combustion devices. This paper describes how the code is presently being used to simulate combustion in a single element combustion chamber with shear coaxial injectors using gaseous oxygen and gaseous hydrogen propellants. The ultimate purpose of the efforts documented is to assess and further improve the Loci-CHEM code and the implementation of it. Single element shear coaxial injectors were tested as part of the Staged Combustion Injector Technology (SCIT) program, where detailed chamber wall heat fluxes were measured. Data was taken over a range of chamber pressures for propellants injected at both ambient and elevated temperatures. Several test cases are simulated as part of the effort to demonstrate use of the Loci-CHEM CFD code and to enable us to make improvements in the code as needed. The simulations presented also include a grid independence study on hybrid grids. Several two-equation eddy viscosity low Reynolds number turbulence models are also evaluated as part of the study. All calculations are presented with a comparison to the experimental data. Weaknesses of the code relative to test data are discussed and continuing efforts to improve the code are presented.

Westra, Douglas G.↗

Analysis of Unsteady Flow Phenomena: Shock - Vortex and Shock - Boundary Layer Interactions

The interaction of shock waves with vortices has received much attention in the past, mainly because shock-vortex interaction closely models the interaction of a shock wave with the coherent structures of a turbulent flow-field, and is a key feature in the broad-band shock noise for supersonic jets in off-project conditions. Chu and Kovasznay have shown that a weak disturbance in a viscous heat conducting fluid can be decomposed as the sum of three basic modes, namely acoustic, vortical and entropy mode; the interaction of any of these modes with a shock wave gives rise to all three disturbance modes downstream of the shock. The vortical mode is important since it constitutes the basis of the coherent structures that have been observed to dominate turbulence for low- to moderate-flow speed. Hollingsworth et al. have experimentally investigated the interaction of a cylindrical shock-induced starting vortex with a plane normal shock, and have shown that the interaction generates a cylindrical acoustic pulse that exhibits a quadrupolar structure consisting of four alternate compression and expansion regions centered around the transmitted vortex. The investigations of Hollingsworth and Richards have been extended by Dosanjh and Weeks that have analyzed the interaction of a columnar spiral vortex with a normal shock wave, thus obtaining quantitative measurements and confirming the generation of a progressive cylindrical wavefront of alternate compression-expansion nature. Naumann and Hermanns' have experimentally addressed the non-linear aspects of shock-vortex interaction, and have shown that the interaction causes both a diffraction and a reflection of the shock with a pattern consisting of either a regular-or a Mach-reflection depending on the shock and the vortex strengths. An attempt to theoretically explain the production of sound from the shock-vortex interaction was carried out by Ribner. Pao and Salas have numerically studied two-dimensional shock-vortex interactions.

Grasso, Francesco↗

GPU-enabled extreme-scale turbulence simulations: Fourier pseudo-spectral algorithms at the exascale using OpenMP offloading

Fourier pseudo-spectral methods for nonlinear partial differential equations are of wide interest in many areas of advanced computational science, including direct numerical simulation of three-dimensional (3-D) turbulence governed by the Navier-Stokes equations in fluid dynamics. This paper presents a new capability for simulating turbulence at a new record resolution up to 35 trillion grid points, on the world's first exascale computer, Frontier, comprising AMD MI250x GPUs with HPE's Slingshot interconnect and operated by the US Department of Energy's Oak Ridge Leadership Computing Facility (OLCF). Key programming strategies designed to take maximum advantage of the machine architecture involve performing almost all computations on the GPU which has the same memory capacity as the CPU, performing all-to-all communication among sets of parallel processes directly on the GPU, and targeting GPUs efficiently using OpenMP offloading for intensive number-crunching including 1-D Fast Fourier Transforms (FFT) performed using AMD ROCm library calls. With 99% of computing power on Frontier being on the GPU, leaving the CPU idle leads to a net performance gain via avoiding the overhead of data movement between host and device except when needed for some I/O purposes. Memory footprint including the size of communication buffers for MPI_ALLTOALL is managed carefully to maximize the largest problem size possible for a given node count. Detailed performance data including separate contributions from different categories of operations to the elapsed wall time per step are reported for five grid resolutions, from 2048 3 on a single node to 32768 3 on 4096 or 8192 nodes out of 9408 on the system. Both 1D and 2D domain decompositions which divide a 3D periodic domain into slabs and pencils respectively are implemented. The present code suite (labeled by the acronym GESTS, GPUs for Extreme Scale Turbulence Simulations) achieves a figure of merit (in grid points per second) exceeding goals set in the Center for Accelerated Application Readiness (CAAR) program for Frontier. The performance attained is highly favorable in both weak scaling and strong scaling, with notable departures only for 2048 3 where communication is entirely intra-node, and for 32768 3 , where a challenge due to small message sizes does arise. Communication performance is addressed further using a lightweight test code that performs all-to-all communication in a manner matching the full turbulence simulation code. Performance at large problem sizes is affected by both small message size due to high node counts as well as dragonfly network topology features on the machine, but is consistent with official expectations of sustained performance on Frontier. Overall, although not perfect, the scalability achieved at the extreme problem size of 32768 3 (and up to 8192 nodes — which corresponds to hardware rated at just under 1 exaflop/sec of theoretical peak computational performance) is arguably better than the scalability observed using prior state-of-the-art algorithms on Frontier's predecessor machine (Summit) at OLCF. New science results for the study of intermittency in turbulence enabled by this code and its extensions are to be reported separately in the near future.

3D fast Fourier transform↗

Edge radiated temperature fluctuations across confinement regime transitions in favorable and unfavorable drift configurations at ASDEX Upgrade

Abstract The physics behind the difference in heating power threshold for the low to high confinement (L–H) transition between favorable ion B × ∇ B drift (toward active X -point) and unfavorable ion B × ∇ B drift (away from active X -point) magnetic configurations is a longstanding open question. In this work we show differences in edge ( ρ p o l = 0.95 –1.0) turbulence in matched heating power favorable and unfavorable magnetic configuration discharges using a correlation electron cyclotron emission diagnostic to measure radiated temperature ( T r a d ) fluctuations. During power matched L-mode phases, the edge T r a d fluctuations have higher amplitude in the unfavorable configuration plasma as compared to the favorable configuration plasma. The lower fluctuation amplitude in favorable configuration is observed concomitant with increased E × B shear, as compared to the unfavorable case. Although this turbulence is different in amplitude, it is similar in radial and frequency structure between favorable and unfavorable configurations and grows to the weakly coherent mode during the I-mode phase of the unfavorable configuration discharge. A broadband pedestal fluctuation is seen during an inter-ELM H-mode period of the favorable configuration discharge.

Physics↗

Strong scintillations in astrophysics. I - The Markov approximation, its validity and application to angular broadening

The Markov approximation to the propagation of waves in an extended, irregular medium is discussed in an astrophysical context. A new derivation is presented which is simple and which shows that the assumption of Gaussian statistics used by previous authors is irrelevant. We discuss the relevance of the approximation and show that it may apply in many situations of interest, including interstellar scintillations of pulsar signals. The approximation does not require the assumption of weak scattering or Gaussian correlation functions. The Markov equation for the angular spectrum is particularly simple, and solutions are discussed for typical turbulence spectra. It is found that the equation for the angular spectrum is very nearly that used by previous authors, and the present discussion shows that these results are much more general than previously thought. A possible observational test for distinguishing between Gaussian and power-law interstellar density spectra is discussed.

Lee, L. C.↗

Particle acceleration by electromagnetic ion cyclotron turbulence

The LF EM-turbulence which furnishes energy for the acceleration of ions in various regions of the earth's magnetosphere efficiently accomplishes its transfer of energy from waves to particles through ion cyclotron resonance (ICR) with the left-hand polarized component of the turbulence; the result of this interaction is a heating of the particle distribution. A general theoretical treatment of ICR heating in a weakly inhomogeneous magnetic geometry is presented, en route to a more detailed examination of auroral ion conics' formation. A substantial simplification of the analysis of the altitude-asymptotic form of the conic distribution is obtained via the similarity transformation introduced into the properties of the electric field spectral density and the earth's dipolar magnetic field.

Crew, G. B.↗

Interaction of the sonic boom with atmospheric turbulence

Theoretical research has been carried out to study the effect of free-stream turbulence on sonic boom pressure fields. A new transonic small-disturbance model to analyze the interactions of random disturbances with a weak shock has been developed. The model equation has an extended form of the classic small-disturbance equation for unsteady transonic aerodynamics. An alternative approach shows that the pressure field may be described by an equation that has an extended form of the classic nonlinear acoustics equation that describes the propagation of sound beams with narrow angular spectrum. The model shows that diffraction effects, nonlinear steepening effects, focusing and caustic effects and random induced vorticity fluctuations interact simultaneously to determine the development of the shock wave in space and time and the pressure field behind it. A finite-difference algorithm to solve the mixed-type elliptic-hyperbolic flows around the shock wave has also been developed. Numerical calculations of shock wave interactions with various deterministic and random fluctuations will be presented in a future report.

Rusak, Z.↗

Interaction of the sonic boom with atmospheric turbulence

Theoretical research was carried out to study the effect of free-stream turbulence on sonic boom pressure fields. A new transonic small-disturbance model to analyze the interactions of random disturbances with a weak shock was developed. The model equation has an extended form of the classic small-disturbance equation for unsteady transonic aerodynamics. An alternative approach shows that the pressure field may be described by an equation that has an extended form of the classic nonlinear acoustics equation that describes the propagation of sound beams with narrow angular spectrum. The model shows that diffraction effects, nonlinear steepening effects, focusing and caustic effects and random induced vorticity fluctuations interact simultaneously to determine the development of the shock wave in space and time and the pressure field behind it. A finite-difference algorithm to solve the mixed type elliptic-hyperbolic flows around the shock wave was also developed. Numerical calculations of shock wave interactions with various deterministic and random fluctuations will be presented in a future report.

Rusak, Zvi↗

The microscopic state of the solar wind: Links between composition, velocity distributions, waves and turbulence

An overview is given of the microscopic state of the solar wind with emphasis on recent Ulysses high-latitude observations and previous Helios in-ecliptic observations. The possible links between composition, ionization state. velocity distribution functions of electrons, protons and heavy ions. kinetic plasma waves and MHD-scale turbulence are elaborated. Emphasis is placed on a connection of interplanetary kinetic-scale phenomena with their generating microscopic processes in the corona. The fast streams seem to consist of mesoscale pressure-balanced plasma filaments and magnetic flux tubes, reminiscent of the supergranular-size structures building the open corona, from which copious Alfven waves emanate. The wind from the magnetically structured and active corona shows developed compressive turbulence and considerable abundance and ionization state variations. Some modelling attempts to explain the observed element fractionation are briefly reviewed. The causes of the nonthermal particle features. such as proton-ion differential streaming, ion beams. temperature anisotropies, and skewed distributions associated with collisionless heat conduction, are ultimately to be searched in the fact, that the corona is never quiescent but fundamentally variable in space and time. Consequently, the radial evolution of the internal state of the wind resembles at all latitudes a complicated relaxation process, in the course of which the free (in comparison with LTE conditions) particle kinetic energy is converted into plasma waves and MHD turbulence on a wide range of scales. This leads to intermittent wave-particle interactions and unsteady anomalous transport, mixed with the weak effects of the rare Coulomb collisions. Spherical expansion and large-scale inhomogeneity forces the wind to attain microscopically a complex internal state of dynamic statistical equilibrium.

Marsch, E.↗

Insights of Boundary Layer Turbulence Over the Complex Terrain of Central Himalaya from GVAX Field Campaign

Limited observations hinder understanding of turbulent characteristics in mountainous terrain resulting from heating or cooling of slopes, wind, vertical motions, and heat or moisture advection, which disperse aerosols and other pollutants over the region. In this study, the 1290 MHz radar wind profiler data are utilized to compute the boundary layer height (BLH), the refractive index structure constant (C n 2 ), and the energy dissipation rate (ε) over the central Himalayan site for the period of November 2011 to March 2012, from the intense Ganges Valley Aerosol Experiment (GVAX) field measurements. The radar wind profiler (RWP) based estimation of BLH and ε is validated against the radiosonde, representing the effectiveness of the datasets for further investigation. The strong seasonal variation of log C n 2 and log ε, with average values of ≈ -12 m -2/3 and -2 m 2 s -3 , respectively, is associated with the mountain-induced local circulations and stability in the atmospheric boundary layer. The weak stratification during weak flow is found to be responsible for deep mixing, particularly in the nocturnal boundary layer in spring. Furthermore, the level of cloud cover significantly impacts the strength of turbulence, with the highest cloud cover resulting in a substantial increase in log C n 2 (approximately -11 m -2/3 ) due to intense updraft and downdraft motions compared to clear skies. Additionally, the distribution of aerosol loading across the site, coupled with the behavior of BLH, atmospheric stability, and orographic-induced circulations, implies distinctive seasonal mechanisms for transporting aerosols toward the mountains. This study offers valuable insights into the diurnal and seasonal patterns of turbulent mixing and the mechanisms behind the transport of pollutants through boundary layer processes over the region.

54 ENVIRONMENTAL SCIENCES↗

Oscillation-convection coupling: Cause of supergranulation

A semianalytic method is derived for dealing simultaneously with large numbers of linear stellar oscillation modes trapped in a cavity (a shell) of fluid which is rotating and convecting. A simple generalization of mixing-length theory shows how convection is modulated by weak rotational effects and by the horizontal wind fields of linear r-mode oscillations. The modulated convection is then used to compute the energy lost to turbulent viscosity by a family of nondegenerate oscillations. Viscosity terms of fourth degree in the wind shear can be included if they are a perturbation affecting only a small portion of the r-mode. Viscous energy loss strenghthens convection in a narrow layer near the base of the H and He ionization zone. In the Sun, this layer is about 7 Mm thick and centered at 0.932 of a solar radius where convection cells have a typical size of about 20 Mm and a lifetime of 0.3 Ms, both similar to what is observed in supergranules. If the rms velocity of r-modes at the surface exceeds 5 m/s, then energy is deposited inside the Sun at a sufficient rate to power the supergranulation and impose on it a weak latitude dependence.

Wolff, C.↗