Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “wave 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 1,027 records · Page 57

Performance of high-order Godunov-type methods in simulations of astrophysical low Mach number flows

High-order Godunov methods for gas dynamics have become a standard tool for simulating different classes of astrophysical flows. Their accuracy is mostly determined by the spatial interpolant used to reconstruct the pair of Riemann states at cell interfaces and by the Riemann solver that computes the interface fluxes. In most Godunov-type methods, these two steps can be treated independently, so that many different schemes can in principle be built from the same numerical framework. Because astrophysical simulations often test out the limits of what is feasible with the computational resources available, it is essential to find the scheme that produces the numerical solution with the desired accuracy at the lowest computational cost. However, establishing the best combination of numerical options in a Godunov-type method to be used for simulating a complex hydrodynamic problem is a nontrivial task. In fact, formally more accurate schemes do not always outperform simpler and more diffusive methods, especially if sharp gradients are present in the flow. For this work, we used our fully compressible Seven-League Hydro (SLH) code to test the accuracy of six reconstruction methods and three approximate Riemann solvers on two- and three-dimensional (2D and 3D) problems involving subsonic flows only. We considered Mach numbers in the range from 10 −3 to 10 −1 , which are characteristic of many stellar and geophysical flows. In particular, we considered a well-posed, 2D, Kelvin–Helmholtz instability problem and a 3D turbulent convection zone that excites internal gravity waves in an overlying stable layer. Although the different combinations of numerical methods converge to the same solution with increasing grid resolution for most of the quantities analyzed here, we find that (i) there is a spread of almost four orders of magnitude in computational cost per fixed accuracy between the methods tested in this study, with the most performant method being a combination of a low-dissipation Riemann solver and a sextic reconstruction scheme; (ii) the low-dissipation solver always outperforms conventional Riemann solvers on a fixed grid when the reconstruction scheme is kept the same; (iii) in simulations of turbulent flows, increasing the order of spatial reconstruction reduces the characteristic dissipation length scale achieved on a given grid even if the overall scheme is only second order accurate; (iv) reconstruction methods based on slope-limiting techniques tend to generate artificial, high-frequency acoustic waves during the evolution of the flow; and (v) unlimited reconstruction methods introduce oscillations in the thermal stratification near the convective boundary, where the entropy gradient is steep.

79 ASTRONOMY AND ASTROPHYSICS↗

The collapse and three-dimensional explosion of three-dimensional massive-star supernova progenitor models

ABSTRACT The explosion outcome and diagnostics of core-collapse supernovae depend sensitively on the nature of the stellar progenitor, but most studies to date have focused exclusively on one-dimensional, spherically symmetric massive star progenitors. We present some of the first core-collapse supernovae simulations of three-dimensional massive star supernovae progenitors, a 12.5- and a 15-M⊙ model, evolved in three dimensions from collapse to bounce through explosion with the radiation-hydrodynamic code fornax. We compare the results using those starting from three-dimensional progenitors to three-dimensional simulations of spherically symmetric, one-dimensional progenitors of the same mass. We find that the models evolved in three dimensions during the final stages of massive star evolution are more prone to explosion. The turbulence arising in these multidimensional initial models serves as seed turbulence that promotes shock revival. Detection of gravitational waves and neutrinos signals could reveal signatures of pre-bounce turbulence.

79 ASTRONOMY AND ASTROPHYSICS↗

Thermal stability effects on the structure of the velocity field above an air-water interface

Mean velocity and turbulence measurements are described for turbulent flows above laboratory water waves, under various wind and thermal stratification conditions. Experimental results, when presented in the framework of Monin-Obukhov (1954) similarity theory, support local scaling based on evaluation of stratification effects at the same nondimensional distance from the mean water surface. Such scaling allows an extension of application of the above theory to the outer region of the boundary layer. Throughout the fully turbulent region, ratios of mean velocity gradients, eddy viscosities, and turbulence intensities under nonneutral and neutral conditions correlate well with the parameter z/Lambda (Lambda being a local Obukhov length and z the vertical coordinate of the mean air flow) and show good agreement with established field correlations. The influence of stratification on the wind-stress coefficient can be estimated from an empirical relationship in terms of its value under neutral conditions and a bulk Richardson number.

Papadimitrakis, Y. A.↗

A numerical study of a class of TVD schemes for compressible mixing layers

At high Mach numbers the two-dimensional time-developing mixing layer develops shock waves, positioned around large-scale vortical structures. A suitable numerical method has to be able to capture the inherent instability of the flow, leading to the roll-up of vortices, and also must be able to capture shock waves when they develop. Standard schemes for low speed turbulent flows, for example spectral methods, rely on resolution of all flow-features and cannot handle shock waves, which become too thin at any realistic Reynolds number. The performance of a class of second-order explicit total variation diminishing (TVD) schemes on a compressible mixing layer problem was studied. The basic idea is to capture the physics of the flow correctly, by resolving down to the smallest turbulent length scales, without resorting to turbulence or sub-grid scale modeling, and at the same time capture shock waves without spurious oscillations. The present study indicates that TVD schemes can capture the shocks accurately when they form, but (without resorting to a finer grid) have poor accuracy in computing the vortex growth. The solution accuracy depends on the choice of limiter. However a larger number of grid points are in general required to resolve the correct vortex growth. The low accuracy in computing time-dependent problems containing shock waves as well as vortical structures is partly due to the inherent shock-capturing property of all TVD schemes. In order to capture shock waves without spurious oscillations these schemes reduce to first-order near extrema and indirectly produce clipping phenomena, leading to inaccuracy in the computation of vortex growth. Accurate simulation of unsteady turbulent fluid flows with shock waves will require further development of efficient, uniformly higher than second-order accurate, shock-capturing methods.

Sandham, N. D.↗

RNG in turbulence and modeling of bypass transition

Two projects are considered: the Renormalization Group (RNG) analysis of turbulence modeling, and the calculation of bypass transition through turbulence modeling. RNG is a process which eliminates small scales on the uneliminated large scales as the change in the transport properties. It is because of this property of RNG that it was previously suggested that RNG could be used as a model builder in turbulence modeling. The possibility is studied of constructing RNG based turbulence models, and to try to proceed to do the modeling through RNG in parallel with the classical approach. The numerical predictions made by RNG models and by classical models is compared against data from Direct Numerical Simulation. While in an environment with freestream turbulence, the transition initiated by the instability of the laminar boundary layer to Tollmien-Schlichting waves is found to be a bypass one in which turbulent spots are formed without T-S wave amplification. The formation is a random process, and flow within a turbulent spot is almost fully turbulent. This suggests the possibility of using turbulence modeling to describe and predict the bypass transition.

Yang, Zhigang↗

Turbulent resistivity, diffusion and heating

Experimental and theoretical studies are reported on ion acoustic and ion cyclotron turbulence and their roles in anomalous resistivity, viscosity, diffusion and heating and in the structure of collisionless electrostatic shocks. Resistance due to ion acoustic turbulence has been observed in experiments with a streaming cesium plasma in which electron current, potential rise due to turbulent resistivity, spectrum of unstable ion acoustic waves, and associated electron heating were all measured directly. Kinetic theory calculations for an expanding, unstable plasma, give results in agreement with the experiment. In a strong magnetic field, with T sub e/T sub i approximately 1 and current densities typical for present Tokomaks, the plasma is stable to ion acoustic but unstable to current driven electrostatic ion cyclotron waves. Relevant characteristics of these waves are calculated and it is shown that for ion, beta greater than m sub e/m sub i, the electromagnetic ion cyclotron wave has a lower instability threshold than the electrostatic one. However, when ion acoustic turbulence is present experiments with double plasma devices show rapid anomalous heating of an ion beam streaming through a plasma.

Fried, B. D.↗

Viscous effects in tube flow initiated by an expansion wave

An examination is made of the effect of the turbulent boundary layer on the uniformity of the flow initiated by a nonsteady expansion wave traveling in a long, constant-area tube. Closed-form expressions for the flow perturbations first are obtained from the linearized equations of motion. Measured static pressure histories for both lowand high-tube Mach numbers are used for comparison. Detailed discussions of the limitations of the linearized model and its extension to high-tube Mach numbers are given. Agreement is shown between a modified version of the linearized model and the experimental results. A model of viscous nonsteady pipe flow is used to explain the nonlinear pressure decay observed at very large length-to-diameter ratios.

Tong, K.-O.↗

A numerical study of fundamental shock noise mechanisms

The results of this thesis demonstrate that direct numerical simulation can predict sound generation in unsteady aerodynamic flows containing shock waves. Shock waves can be significant sources of sound in high speed jet flows, on helicopter blades, and in supersonic combustion inlets. Direct computation of sound permits the prediction of noise levels in the preliminary design stage and can be used as a tool to focus experimental studies, thereby reducing cost and increasing the probability of a successfully quiet product in less time. This thesis reveals and investigates two mechanisms fundamental to sound generation by shocked flows: shock motion and shock deformation. Shock motion is modeled by the interaction of a sound wave with a shock. During the interaction, the shock wave begins to move and the sound pressure is amplified as the wave passes through the shock. The numerical approach presented in this thesis is validated by the comparison of results obtained in a quasi-one dimensional simulation with linear theory. Analysis of the perturbation energy demonstrated for the first time that acoustic energy is generated by the interaction. Shock deformation is investigated by the numerical simulation of a ring vortex interacting with a shock. This interaction models the passage of turbulent structures through the shock wave. The simulation demonstrates that both acoustic waves and contact surfaces are generated downstream during the interaction. Analysis demonstrates that the acoustic wave spreads cylindrically, that the sound intensity is highly directional, and that the sound pressure level increases significantly with increasing shock strength. The effect of shock strength on sound pressure level is consistent with experimental observations of shock noise, indicating that the interaction of a ring vortex with a shock wave correctly models a dominant mechanism of shock noise generation.

Meadows, Kristine R.↗

Critical Fluctuations in Beam-Plasma Systems and Solar Type III Radio Bursts

The type III radio bursts are the most intense radio emissions from the sun. In Figure 1, we present a typical type III burst observed by the STEREO spacecraft. Ginzburg and Zheleznyakov (1958) were the first to suggest that Langmuir waves excited by the solar flare accelerated electrons are the source of these bursts. The in situ detection of electron beams and Langmuir waves in association with type III bursts together with tracking of type III burst sources in the interplanetary medium confirmed this plasma hypothesis. However, the dynamics of type III electron beams and their interaction with the ambient plasma through the excited Langmuir waves is an unsolved problem. In order that the electron beams don't lose their energy by resonantly interacting with Langmuir waves, various nonlinear processes, such as the induced scattering of Langmuir waves off the background ions in weak turbulence regime and oscillating two-stream instability in strong turbulence regime were invoked. However, the observed electron density fluctuations do not allow Langmuir waves to grow to very high intensities as pointed out by several authors. As far as the conversion of Langmuir wave energy into electromagnetic energy is concerned, there is no consensus either for the fundamental or for the second harmonic emission. In this paper, we address these issues in terms of critical fluctuations near the boundary of the beam-plasma instability which is analogous to the anomalously growing fluctuations near the phase transition points called critical points. We show that this phenomenon can account for the survival of electron beams over large distances as well as the intensities of Langmuir waves excited during this phase transition regime are sufficient to explain the peak intensities of the fundamental as well as harmonic emissions

Fluctuations↗

Stochastic acceleration of electrons. I - Effects of collisions in solar flares

Stochastic acceleration of thermal electrons to nonrelativistic energies is studied under solar flare conditions. We show that, in turbulent regions, electron-whistler wave interactions can result in the acceleration of electrons in times comparable to or shorter than the Coulomb collision time. The kinetic equation describing the evolution of the electron energy distribution including stochastic acceleration by whistlers and energy loss via Coulomb interactions is solved for an initial thermal electron energy spectrum. In general, the shape of the resulting electron distributions are characterized by the energy E(c) where systematic energy gain by turbulence equals energy loss due to Coulomb collisions. For energies less than E(c), the spectra are steep (quasi-thermal) whereas above E(c), the spectra are power laws. We find that hard X-ray spectra computed using the electron distributions obtained from our numerical simulations are able to explain the complex spectral shapes and variations observed in impulsive hard X-ray bursts. In particular, we show that the gradual steepening observed by Lin et al. (1981) could be due to a systematic increase in the density of the plasma (due to evaporation) and the increasing importance of collisions instead of the appearance of a superhot thermal component.

Hamilton, Russell J.↗

On the linear approximation of gravity wave saturation in the mesosphere

Lindzen's model of gravity wave breaking is shown to be inconsistent with the process of convective adjustment and associated turbulent outbreak. The K-theory turbulent diffusion model used by Lindzen implies a spatially uniform turbulent field which is not in agreement with the fact that gravity wave saturation and the associated convection produce turbulence only in restricted zones. The Lindzen model may be corrected to some extent by taking the turbulent Prandtl number for a diffusion acting on the wave itself to be very large. The eddy diffusion coefficients computed by Lindzen then become a factor of 2 larger and eddy transports of heat and constituents by wave fields vanish to first order.

Chao, W. C.↗

Local effects of gravity wave propagation and saturation

In recent years, gravity waves were recognized to play a major role in the dynamics of the middle atmosphere. Perhaps the major effect of such motions are the reversal of the vertical shear of the mean zonal wind and the occurrence of a large turbulent diffusivity in the mesosphere due to gravity wave saturation. Yet, despite the importance of these gravity wave effects, the processes and the consequences of gravity wave propagation and saturation are only beginning to be understood in detail. The linear saturation theory predicts drag and turbulent diffusion due to saturating wave motion. This theory, however, fails to address a number of issues that are certain to be important for gravity wave propagation and saturation in the middle atmosphere. These issues, including wave transience, wave superposition, local convective adjustment, and nonlinearity, are discussed.

Fritts, D. C.↗

Solar Wind Turbulence Around Mars: Relation Between The Energy Cascade Rate And The Proton Cyclotron Waves Activity

The first estimation of the incompressible energy cascade rate at magnetohydrodynamic (MHD) scales in the plasma upstream of the Martian bow shock is obtained, making use of magnetic field and plasma observations provided by Mars Atmosphere and Volatile EvolutioN (MAVEN) over 600 orbits. In particular, the energy cascade rate is computed for events with and without proton cyclotron wave (PCW) activity, for time intervals when MAVEN was in the solar wind with no magnetic connection to the bow shock. It is shown that the nonlinear cascade of energy at the MHD scales is slightly amplified when PCWs are present in the plasma, around the Martian perihelion. In addition, the analysis of the normalized cross helicity and residual energy for the turbulent fluctuations shows the presence of Alfvenic and non-Alfvenic fluctuations in a magnetic dominant regime for the majority of the cases.

solar wind↗

Atmospheric turbulence chamber for optical transmission experiment Characterization by thermal method

Consideration is given to an atmospheric turbulence chamber designed for optical wave propagation experiments. The chamber consists of ten small electric heater/blowers with an aluminum foil screen and three screens of 2-mm aluminum wire meshes. Calculations are made of the temperature structure constant squared on the basis of temperature structure function measurements derived from a differential microthermocouple system. Values are presented for the refractive-index structure constant squared. The average wind velocity and temperature are found to be, respectively, 0.41 m/sec and 53 C. The inner and outer scales of turbulence are 5.0 mm and 6.5 cm. It is shown that the measured temperature structure function and the power spectrum of temperature fluctuations satisfy, respectively, the 2/3 and -5/3 power similarity laws in the inertial subrange. Possible chamber improvements are discussed.

Gamo, H.↗

Interaction of radio frequency waves with cylindrical density filaments: scattering and radiation pressure

The propagation of radio-frequency (RF) waves in tokamaks can be affected by filamentary structures, or blobs, that are present in the edge plasma and the scrape-off layer. The difference in the permittivity between the surrounding plasma and interior of a filament leads to reflection, refraction and diffraction of the waves. This, in turn, can affect the power flow into the core of the plasma and reduce the efficiency of heating and/or current generation. The scattering of RF waves, lower hybrid, helicon and ion cyclotron waves, by a single cylindrical filament, embedded in a background plasma, is studied using a full-wave analytical theory developed previously (Ram & Hizanidis, Phys. Plasmas , vol. 23, 2016, 022504). The theory assumes that the plasma in and around a filament is homogeneous and cold. A detailed scattering analysis reveals a variety of common features that exist among the three distinctly different RF waves. These common attributes can be inferred intuitively based on an examination of the cold plasma dispersion relation. The physical intuition is a useful step to understanding experimental observations on scattering, as well as results from simulations that include general forms of edge plasma turbulence. While a filament can affect the propagation of RF waves, the radiation force exerted by the waves can influence the filament. The force on a filament is determined using the Maxwell stress tensor. In 1905, Poynting was the first to evaluate and measure the radiation force on an interface separating two different dielectric media (Poynting, London Edinburgh Dublin Philos. Mag. J. Sci. , vol. 9, 1905, pp. 393–406). For ordinary light propagating in vacuum and incident on a glass surface, Poynting noted that the surface is ‘pulled’ towards the vacuum. In a magnetized cold plasma, there are two independent wave modes. Even if only one of these modes is excited by an RF antenna, a filament will couple power to the other mode: a consequence of electromagnetic boundary conditions. This facet of scattering has consequences on the radiation force that go beyond Poynting's seminal contribution. The direction of the force depends on the polarization of the incident wave and on the mode structure of the waves inside and in the vicinity of a filament. It can either pull the filament toward the RF source or push it away. For slow lower hybrid waves, filaments with densities greater than the ambient density are pulled in, while filaments with lower densities are pushed out, thereby enhancing the density in front of the antenna. In the case of fast helicon and ion cyclotron waves, the direction of the force depends on the plasma and wave parameters; in particular, on the ambient density. The radiation force, in all three frequency ranges, is large enough to affect the motion of a filament and could be measured experimentally. This also suggests the possibility of modifying the edge turbulence using RF waves.

Physics↗

Interaction of radio frequency waves with cylindrical density filaments -- scattering and radiation pressure

The propagation of radio-frequency (RF) waves in tokamaks can be affected by filamentary structures, or blobs, that are present in the edge plasma and the scrape-off layer. The difference in the permittivity between the surrounding plasma and interior of a filament leads to reflection, refraction and diffraction of the waves. This, in turn, can affect the power flow into the core of the plasma and reduce the efficiency of heating and/or current generation. The scattering of RF waves, lower hybrid, helicon and ion cyclotron waves, by a single cylindrical filament, embedded in a background plasma, is studied using a full-wave analytical theory developed previously (Ram & Hizanidis, Phys. Plasmas, vol. 23, 2016, 022504). The theory assumes that the plasma in and around a filament is homogeneous and cold. A detailed scattering analysis reveals a variety of common features that exist among the three distinctly different RF waves. These common attributes can be inferred intuitively based on an examination of the cold plasma dispersion relation. The physical intuition is a useful step to understanding experimental observations on scattering, as well as results from simulations that include general forms of edge plasma turbulence. While a filament can affect the propagation of RF waves, the radiation force exerted by the waves can influence the filament. The force on a filament is determined using the Maxwell stress tensor. In 1905, Poynting was the first to evaluate and measure the radiation force on an interface separating two different dielectric media (Poynting, London Edinburgh Dublin Philos. Mag. J. Sci., vol. 9, 1905, pp. 393–406). For ordinary light propagating in vacuum and incident on a glass surface, Poynting noted that the surface is ‘pulled’ towards the vacuum. In a magnetized cold plasma, there are two independent wave modes. Even if only one of these modes is excited by an RF antenna, a filament will couple power to the other mode: a consequence of electromagnetic boundary conditions. This facet of scattering has consequences on the radiation force that go beyond Poynting's seminal contribution. The direction of the force depends on the polarization of the incident wave and on the mode structure of the waves inside and in the vicinity of a filament. It can either pull the filament toward the RF source or push it away. For slow lower hybrid waves, filaments with densities greater than the ambient density are pulled in, while filaments with lower densities are pushed out, thereby enhancing the density in front of the antenna. In the case of fast helicon and ion cyclotron waves, the direction of the force depends on the plasma and wave parameters; in particular, on the ambient density. The radiation force, in all three frequency ranges, is large enough to affect the motion of a filament and could be measured experimentally. This also suggests the possibility of modifying the edge turbulence using RF waves.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Determination of atmospheric structure function by using a single coherent detector

The wave structure function associated with atmospheric turbulence is important in the effective design of telescopes for coherent optical receivers. This paper describes a method of measurement of the wave structure function by employing an already existing optical receiver. The aperture of this receiver's telescope is masked by an opaque screen with two small holes. It is then shown that when the incident beam obeys Gaussian statistics, the receiver's A.M. output can be simply related to the wave structure function of the incident beam.

Lang, R. H.↗

Space plasma simulations; Proceedings of the Second International School for Space Simulations, Kapaa, HI, February 4-15, 1985. Parts 1 & 2

Space plasma simulations, observations, and theories are discussed. Papers are presented on the capabilities of various types of simulation codes and simulation models. Consideration is given to plasma waves in the earth's magnetotail, outer planet magnetosphere, geospace, and the auroral and polar cap regions. Topics discussed include space plasma turbulent dissipation, the kinetics of plasma waves, wave-particle interactions, whistler mode propagation, global energy regulation, and auroral arc formation.

Ashour-Abdalla, M.↗