SEARCH · Engineering Papers
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.
The flow over wave shaped surfaces is studied using a Navier Stokes solver. Detailed comparisons with theoretical results are presented, including the stability of a laminar flow over wavy surfaces. Drag characteristics of nonplanar surfaces are predicted using the Navier-Stokes solver. The secondary instabilities of wall bounded and free shear flows are also discussed.
Radiosonde and U-2 aircraft data from the Intertropical Convergence Zone experiment of July, 1977 were analyzed. The radiosonde data above 15 km was plotted in the form of time-height sections so that wave motions in the meteorological variables (temperature, zonal wind, and meridional wind) could be easily delineated. The results showed the presence of a planetary 16 day Kelvin wave. More significantly, large amplitude gravity waves with vertical wavelengths on the order of 3-6 km were found. Spectral analysis of frequencies showed peaks in meridional wind variance at 2-3 days and 5 days, with the shorter periods stronger at higher altitudes. In the time domain, two instances of organized wave groups propagating upward into the stratosphere were isolated. Associated with one of these wave groups was a region of very large vertical shear in which either turbulence or small scale (10 km wavelength or less) gravity wave activity had been found by the U-2 aircraft.
Supersonic inlets with external compression, having a good level performance at the critical operating point, exhibit a marked instability of the flow in some subcritical operation below a critical value of the capture mass flow ratio. This takes the form of severe oscillations of the shock system, commonly known as "buzz". The underlying purpose of this study is to indicate how Detached Eddy Simulation (DES) analysis of supersonic inlets will alter how we envision unsteady inlet aerodynamics, particularly inlet buzz. Presented in this paper is a discussion regarding the physical explanation underlying inlet buzz as indicated by DES analysis. It is the normal shock wave boundary layer separation along the spike surface which reduces the capture mass flow that is the controlling mechanism which determines the onset of inlet buzz, and it is the aerodynamic characteristics of a choked nozzle that provide the feedback mechanism that sustains the buzz cycle by imposing a fixed mean corrected inlet weight flow. Comparisons between the DES analysis of the Lockheed Martin Corporation (LMCO) N+2 inlet and schlieren photographs taken during the test of the Gulfstream Large Scale Low Boom (LSLB) inlet in the NASA 8x6 ft. Supersonic Wind Tunnel (SWT) show a strong similarity both in turbulent flow field structure and shock wave formation during the buzz cycle. This demonstrates the value of DES analysis for the design and understanding of supersonic inlets.
The ‘Acceleration, Reconnection, Turbulence and Electrodynamics of the Moon's Interaction with the Sun’ (ARTEMIS) mission provides a unique opportunity to study the structure of interplanetary shocks and the associated generation of plasma waves with frequencies between ~50-8000 Hz due to its long duration electric and magnetic field burst waveform captures. We compare wave properties and occurrence rates at 11 quasi-perpendicular interplanetary shocks with burst data within 10 minutes (~3200 proton gyroradii upstream, ~1900 downstream) of the shock ramp. A perturbed shock is defined as possessing a large amplitude whistler precursor in the quasi-static magnetic field with an amplitude greater than ⅓ the difference between the upstream and downstream average magnetic field magnitudes; laminar shocks lack these large precursors and have a smooth, step function-like transition. In addition to wave modes previously observed, including ion acoustic, whistler, and electrostatic solitary waves, waves in the ion acoustic frequency range that show rapid temporal frequency change are common. Three shocks had burst captures in the ramp; of these, the two laminar shocks with burst data in the ramp contained a wide range of large amplitude wave modes in the ramp whereas the one perturbed shock contained no such waves. Thus, energy dissipation through wave-particle interactions is more prominent in these two laminar shocks than in the perturbed shock. Based on observations from all 11 shocks, the The wave occurrence rates for laminar shocks are higher in the transition region, especially the ramp, than downstream. In contrast, perturbed shocks have approximately 2-3 times the wave occurrence rate downstream than laminar shocks.
KPP (Kolmogorov-Petrovsky-Piskunov) solutions of the reaction-diffusion equation have application in various physical phenomena occurring in biology, ecology, and reacting flows. In particular, these solutions are commonly used in turbulent combustion to scale turbulent burning velocities. Subject to certain conditions on reaction rate profile through the flame brush and turbulent diffusivity, this theory relates the turbulent burning velocity to the derivative of the reaction rate ($\tilde{ω}$) at the leading edge of the flame brush ($d\tilde{ω}$/$d\tilde{c}$ | $\tilde{c} =$0 ). Such waves are often referred to as “pulled fronts.” However, turbulent flames never actually satisfy the KPP conditions for a pulled front, as the turbulent flame brush, parametrized here by the thickness δ t , consists of an ensemble of laminar flamelets of thickness δ, where ε=δ/δ t $\ll$ 1 is very small, but nonzero, and $d\tilde{ω}$ /$d\tilde{c}$ tends to zero at the brush leading edge for high activation energy, combustion-type kinetics. Here, this paper analyzes these effects on KPP wave solutions, parametrized by ε=δ/δ t and Zeldovich number Ze focusing on whether turbulent flames retain their pulled front character and what the correction to the KPP wave speed is. Variational solutions of the reaction-diffusion equation show that the solution can be expanded in powers of 1/|lnε|. Both numerical and asymptotic results are presented, showing that the wave still exhibits pulled front solutions but with significant corrections to the KPP result. The leading order correction is of the form |lnε| –2 and independent of Ze. Higher order corrections are function of both ε and Ze. However, the dominant factor influencing the wave speed correction is due to the finite ε, with Ze exhibiting a weaker effect.
The role of zonal flow (ZF) in the turbulence saturation of trapped electron modes (TEM) in magnetically confined plasmas is revisited. Here, in this study, we examine ZF excitation and saturation mechanisms in TEM turbulence using detailed free-energy transfer diagnostics from nonlinear gyrokinetic simulations. When the ion channel becomes subdominant, electrons shift to transferring energy to the zonal component, which results in zonal flows always playing a significant though subdominant role, in the case of the temperature gradient driven TEM. The velocity-space structure of energy transfer reveals that trapped electrons directly couple to ZF. Regarding the saturation physics, the zonal flow advection, drift-wave–drift-wave interactions, and stable modes at the pump waves are compared. When ZF saturation is weak, drift-wave–drift-wave interactions play a larger role in maintaining turbulence saturation. The findings elucidate the roles of nonzonal to zonal flow coupling, stable modes, and nonzonal wave-wave interactions, and provide a basis for improving reduced models of turbulent transport.
Solutions are presented for the dispersion relation of waves propagating parallel to the ambient magnetic field, in a plasma having an arbitrary ratio of thermal to magnetic pressure, beta. Previous results applicable only to wave frequencies much lower than the ion cyclotron frequency are extended to include all frequencies lower than the electron cyclotron frequency. The cyclotron turnovers are found to occur at significantly lower frequencies in finite beta plasmas. Cyclotron damping is significant for small wavelength waves. These results show that in energetic particle scattering by electromagnetic turbulence in a plasma with beta value greater than 40, scattering by waves with wave frequencies greater than the ion cyclotron frequency can be neglected. This condition is satisfied in many astrophysical plasmas.
A two-dimensional computational code, PRLUS2D, which was developed for the reactive propulsive flows of ramjets and scramjets, was validated for two-dimensional shock-wave/turbulent-boundary-layer interactions. The problem of compression corners at supersonic speeds was solved using the RPLUS2D code. To validate the RPLUS2D code for hypersonic speeds, it was applied to a realistic hypersonic inlet geometry. Both the Baldwin-Lomax and the Chien two-equation turbulence models were used. Computational results showed that the RPLUS2D code compared very well with experimentally obtained data for supersonic compression corner flows, except in the case of large separated flows resulting from the interactions between the shock wave and turbulent boundary layer. The computational results compared well with the experiment results in a hypersonic NASA P8 inlet case, with the Chien two-equation turbulence model performing better than the Baldwin-Lomax model.
In this work, we are developing techniques for studying the Rayleigh–Taylor (RT) and Richtmyer–Meshkov (RM) instabilities in a planar geometry at high-energy-densities at the National Ignition Facility (NIF). In particular, through the improvement of experimental imaging quality, we are progressing toward the study of the turbulent regime of the mixing regions in capsule implosion experiments for inertial confinement fusion, which requires few micrometers resolution. Using 60 NIF beams, a solid shock tube is driven launching a shock wave that crosses the interface between a dense and a light material pre-machined in the target to obtain sinusoidal ripples, which results in RM and RT instabilities that are imaged using the NIF Crystal Backlighter Imager. High-quality images were obtained with a mean resolution of 7 μm and improved contrast. While the obtained resolution does not allow the observation of the smallest scale of the “turbulent” energy spectrum, the generated image encompasses 63% of the total flow energy, a 50% improvement over previous studies, which is observed for the first time a roll-up feature in a high energy density-type RT experiment.
The vector helium magnetometer on the International Cometary Explorer observed the magnetic fields induced by the interaction of comet Giacobini-Zinner with the solar wind. A magnetic tail was penetrated about 7800 kilometers downstream from the comet and was found to be 10,000 kilometers wide. It consisted of two lobes, containing oppositely directed fields with strengths up to 60 nanoteslas, separated by a plasma sheet about 1000 kilometers thick containing a thin current sheet. The magnetotail was enclosed in an extended ionosheath characterized by intense hydromagnetic turbulence and interplanetary fields draped around the comet. A distant bow wave, which may or may not have been a bow shock, was observed at both edges of the ionoshpeath. Weak turbulence was observed well upstream of the bow wave.
The relaminarization of the boundary layer over a flat plate in the shock tube was investigated by using the partially reflected shock wave technique. The flow Mach number was approximately 0.14, which corresponds to the inleft flow Mach number for the first row of vanes in a gas turbine. The thin film platinum heat gauges were used to measure the heat transfer rate and the Stanton number was calculated from the oscilloscope voltage traces. The Reynolds number was varied by changing the operation pressure of the shock tube and the values varied from 2.3 x 10 to the 4th to 5.3 x 10 to the 5th. For a Reynolds number range of 7 x 10 to the 4th to 3.5 x 10 to the 5th, the relaminarization of the boundary layer was observed. This phenomenon is due to the decay of the turbulence level in the flow as the reflected shock wave moves upstream from the flat plate. As the Reynolds number increased, the relaminarization was delayed and the delay was related to the turbulence generated by the reflected shock wave.
Freestream turbulence of weak but adjustable intensity incident upon a flat plate test model induced Tollmien-Schichting (TS) waves and wave packets. These were studied using a newly-developed sensor by which the packets could be followed from a station nearly as far forward as the minimum critical Reynolds number to the onset of transition. Conclusions reached from these studies were: (1) weak freestream turbulence induced TS wave packets in a Blasius boundary-layer, (2) packets gained strength during propagation through expansion of lateral and longitudinal scales, and through an increase in peak amplitude, (3) the average wave strength at stations of observation was not related to the freestream velocity fluctuation in a linear manner, and (4) packets typically evolved into turbulent spots.
In situ observations of turbulence spectra in space plasmas are usually interpreted as wavenumber spectra, assuming that the fluctuation frequency is negligible in the plasma flow frame. We explore the effects of nonzero frequency in the plasma flow frame on turbulence spectral observations. The finite frequency can be caused by either propagating waves or nonlinear broadening of nonpropagating structures. We show that the observed frequency spectrum can be modified by the nonzero frequency of turbulent fluctuations in several ways. Specifically, (i) frequency broadening results in a minor modification to the observed spectrum, primarily acting as a smoothing kernel of the spectrum near the spectral break, while the asymptotic spectral index remains unchanged; (ii) wave propagation can affect the observed spectral index for anisotropic turbulence. The effect is significant at low frequencies and weaker at high frequencies, leading to a "concave" shape of the observed perpendicular spectrum; (iii) the Doppler shift for forward- and backward-propagating Elsasser modes can result in a nonzero cross helicity for critical-balanced turbulence since the effect of the Doppler shift favors outward-propagating waves systematically, resulting in an observed imbalance. These results may have important implications for the interpretation of solar wind flows observed by Parker Solar Probe.
Structure formation in turbulence can be understood as an instability of “plasma” formed by fluctuations serving as effective particles. These “particles” are quantumlike in the sense that their wavelengths are non-negligible compared to the sizes of background coherent structures. The corresponding “kinetic equation” describes the Wigner matrix of the turbulent field, and the coherent structures serve as collective fields. This formalism is usually applied to manifestly quantumlike or scalar waves. Here, we show how to systematically extend it to more complex systems using compressible Navier–Stokes turbulence as an example. In this case, the fluctuation Hamiltonian is a five-dimensional matrix operator and diverse modulational modes are present. As an illustration, we calculate these modes for a sinusoidal shear flow and find two modulational instabilities. One of them is specific to supersonic flows, and the other one is a Kelvin–Helmholtz-type instability that is a generalization of the known zonostrophic instability. Our calculations are readily extendable to other types of turbulence, for example, magnetohydrodynamic turbulence in plasma.
Transition to turbulence in swept-wing flows has resisted correlation with linear theory because of its sensitivity to freestream conditions and 3-D roughness and because one of the principal instability modes quickly 'becomes nonlinear. In the face of such a formidable problem, two rather long-term fundamental efforts have been underway at DLR Gottinberg and Arizona State University that address swept-wing transition. These efforts have been recently reviewed by Bippes (1997) and Reibert and Saric (1997). Thus, the present work is a continuation of a series of studies on swept-wing boundary layers which have led to a better understanding of the transition process. In particular, we have taken advantage of the sensitivity to 3-D roughness and the modal nature of the instability in order to propose a particular control strategy. Complementing the two aforementioned reviews, general reviews of the swept-wing transition problem are found in Arnal (1997) and Kachanov (1996). Other recent reviews include Reshotko (t997), Crouch (1997), and Herbert (1997a,b). The failure of linear theory is discussed in Reed et al. (1996). The historical work is found in Reed and Sar-ic (1989). The basic idea is that the combination of sweep and chordwise pressure gradient within the boundary layer creates a velocity component perpendicular to the inviscid streamline. This crossflow profile is inflectional and exhibits both traveling and stationary unstable waves called crossflow vortices that are (approximately) aligned along the inviscid streamlines. Under conditions of low freestream turbulence levels, the dominant crossflow wave is stationary (Reibert and Saric t997) while moderate to high turbulence levels initiate dominant traveling waves (Dehle and Bippes 1996; Bippes 1997). 'Me mechanism is relatively insensitive to sound and 2-D surface roughness (Radeztsky et al. 1993) but very sensitive to 3-D roughness near the attachment line. We concentrate our work on low-turbulence freestream flows and stationary crossflow waves. Although the v' and w' components of the disturbances are very small, by convecting streamwise momentum in the wall-normal direction, they produce 0(l) changes inu'.Thus the mean flow is highly distorted with localized inflection points. Transition is then triggered by a high-frequency secondary instability of the distorted mean profile.
The influence of a positive pressure gradient on the boundary layer transition is studied. The mean velocity and turbulence profiles of four cases are examined. As the intensity of the pressure gradient is increased, the Reynolds number of the transition onset and the length of the transition region are reduced. The Tollmein-Schlichting waves disturb the laminar regime; the amplification of these waves is in good agreement with the stability theory. The three dimensional deformation of the waves leads finally to the appearance of turbulence. In the case of zero pressure gradient, the properties of the turbulent spots are studied by conditional sampling of the hot-wire signal; in the case of positive pressure gradient, the turbulence appears in a progressive manner and the turbulent spots are much more difficult to characterize.
The slow shocks which have been observed in the distant geomagnetic tail by ISEE 3 do not exhibit the large-amplitude magnetic field rotations which are theoretically expected. A possible explanation is that the rotation trailing wave train is resistively damped by strong lower hybrid drift turbulence excited in the shock front. Magnetic wave spectra for three slow shocks are constructed from ISEE 3 dc magnetometer and search coil experiments. The observed magnetic amplitudes near the lower hybrid frequency are below the level at which anomalous resistance would damp the slow shock magnetic rotations.
In this study, ensembles of experimental data are presented and utilized to compare and validate two models used in the simulation of variable-density (Atwood = 0.22), compressible turbulent mixing. Though models of this kind (Reynolds averaged NavierStokes and large-eddy simulations) have been validated extensively with more canonical flows in previous studies, here the present approach offers novelty in the complexity of the geometry, the ensemble-based validation, and the uniformity of the computational framework on which the models are tested. Moreover, all experimental and computational tasks were completed by the authors which has led to a tightly coupled experimental configuration with its “digital twin.” The experimental divergent-shock-tube facility and its data acquisition methods are described and replicated in simulation space. A 2D Euler model which neglects the turbulent mixing at the interface is optimized to experimental data using a Gaussian process. This model then serves as the basis for both the 2D RANS and 3D LES studies that make comparisons to the mixing-layer data from the experiment. A relatively simple RANS model is shown to produce good agreement with experimental data only at late flow development times. The LES ensembles generally show good agreement with experimental data but display sensitivity to the characterization of initial conditions. Resolution-dependent behavior is also observed for certain higher-order statistics of interest. Overall, the LES model successfully captures the effects of divergent geometry, compressibility, and combined nonlinear instabilities inherent to the problem. The successful prediction of mixing width and its growth rate highlight the existence of three distinct regimes in the development of the instability, each with similarities to previously studied instabilities.