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 505 records · Page 28

Ground-based sensors for the SR-71 sonic boom propagation experiment

This paper describes ground-level measurements of sonic boom signatures made as part of the SR-71 sonic boom propagation experiment recently completed at NASA Dryden Flight Research Center, Edwards, California. Ground level measurements were the final stage of this experiment which also included airborne measurements at near and intermediate distances from an SR-71 research aircraft. Three types of sensors were deployed to three station locations near the aircraft ground track. Pressure data were collected for flight conditions from Mach 1.25 to Mach 1.60 at altitudes from 30,000 to 48,000 ft. Ground-level measurement techniques, comparisons of data sets from different ground sensors, and sensor system strengths and weaknesses are discussed. The well-known N-wave structure dominated the sonic boom signatures generated by the SR-71 aircraft at most of these conditions. Variations in boom shape caused by atmospheric turbulence, focusing effects, or both were observed for several flights. Peak pressure and boom event duration showed some dependence on aircraft gross weight. The sonic boom signatures collected in this experiment are being compiled in a data base for distribution in support of the High Speed Research Program.

Norris, Stephen R.↗

Ground-Based Sensors for the SR-71 Sonic Boom Propagation Experiment

This paper describes ground-level measurements of sonic boom signatures made as part of the SR-71 sonic boom propagation experiment recently completed at NASA Dryden Flight Research Center, Edwards, California. Ground-level measurements were the final stage of this experiment which also included airborne measurements at near and intermediate distances from an SR-71 research aircraft. The types of sensors were deployed to three station locations near the aircraft ground track. Pressure data collected for flight conditions from Mach 1.25 to Mach 1.60 at altitudes from 30,000 to 48,000 ft. Ground-level measurement techniques, comparisons of data sets from different ground sensors, and sensor system strengths and weaknesses are discussed. The well-known N-wave structure dominated r sonic boom signatures generated by the SR-71 aircraft at most of these conditions. Variations in boom shape caused by atmospheric turbulence, focusing effects, or both, were observed for several flights. Peak pressure and boom event duration showed some dependence on aircraft gross weight. The sonic boom signatures collected in this experiment are being compiled in a data base for distribution in support of the High Speed Research Program.

Norris, Stephen R.↗

Current driven weak double layers

Double layers in plasmas can be created by different means. For example, a potential difference forms between two plasmas with different temperatures, in a plasma jet flowing along a converging magnetic field, in a quiescent plasma submitted to an external difference of potential, or in a turbulent plasma carrying an electric charge. The first three cases can be current-free, but not necessarily, although the numerical simulations were made under such conditions for the first two points. Apart from the third case, which is mainly of interest for laboratory experiments, these double layers are good candidates for accelerating the auroral electrons to the few kiloelectron volts observed.

Chanteur, Gerard↗

Direct numerical simulations of turbulent premixed cool flames: Global and local flame dynamics analysis

The cool flame dynamics, especially in turbulent flows, is of great interest for both practical application and fundamental research. Here, in this study, a series of direct numerical simulations of turbulent premixed n-C 7 H 16 /O 2 /O 3 /N 2 cool flames are performed, with the focus on the influence of turbulence intensity (u'/S L , where S L is the laminar flame speed) on the flame structure as well as the global and local cool flame dynamics. It is found that the cool flame front is considerably wrinkled by turbulence at high u'/S L , leading to significantly thickened turbulent cool flame brush and largely altered local reactivity compared with the reference laminar flame. However, the turbulent flame structure in the temperature space is found to be insensitive to u'/S L . Besides, with increasing u'/S L , the normalized turbulent cool flame speed (S L /S L ) is monotonically increased, attributed to substantial augmentation on the flame surface area (A T /A L ), while the stretching factor (I 0 ) remains almost constant and is smaller than 1. The underlying mechanisms for such variations are revealed through local flame dynamics analysis. Specifically, the local flame displacement speed S d is found to be strongly negatively correlated with flame curvature; meanwhile, such negative correlation and the probability distribution function (PDF) of flame curvature are barely influenced by u'/S L , leading to a weak dependence of I 0 on u'/S L . In contrast, the PDF of the tangential strain rate is found to span a much wider range and shift to the positive side as u'/S L increases, suggesting that the enhanced tangential strain rate is the main cause for the increase in surface area of the turbulent premixed cool flame. Finally, the influence of equivalence ratio on above findings is found to be insignificant, indicating that although the local reactivity of turbulent premixed cool flames is altered due to the differential diffusion, the resultant flame- stretch interaction is insensitive to the equivalence ratio. This study presents some unique cool flame dynamics that are distinct from hot flames, which can help improve the understanding and modeling of turbulent cool flames.

Cool flames↗

The role of breakup and coalescence in fine-scale bubble-induced turbulence. II. Kinematics

This second part of our research explores the kinematic aspect of fine-scale bubble-induced turbulence (BIT) to (i) present the effect of bubble breakup and coalescence and (ii) compare it against the universal kinematic fine-scale turbulence characteristics reported in the literature. To this end, we simulate a dilute bubbly system of 0.5% void fraction using two distinct numerical simulations. In the volume-of-fluid (VoF) simulation, bubbles undergo breakup and coalescence. In the immersed boundary method (IBM) simulation, however, they act as rigid spheres. We also perform a simulation of classical homogeneous isotropic turbulence (HIT). The first important outcome of this study is that BIT is radically different from HIT in terms of its kinematic fine-scale characteristics. In the vorticity-dominating regions, BIT exhibits a weak vortex stretching. This weak vortex stretching is due to (a) the intermediate strain-rate eigenvalues skewed weakly to positive and (b) the extensive strain-rate eigenvector aligning perpendicular to the vorticity vector. The BIT has, on average, not only a weak enstrophy production but also a weak strain production in strain-dominating regions. The weak strain production is due to (a) the presence of vortex stretching in highly strained fluid elements and (b) the absolute magnitude of compressive strain-rate eigenvalue being as close to the extensive strain-rate eigenvalue. Thus, none of the kinematic fine-scale HIT characteristics is noted for BIT. The second important conclusion is that bubble breakup and coalescence play little to no influence on the kinematics of fine-scale BIT as VoF and IBM simulations produce similar results.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Impact of plasma density/collisionality on divertor heat flux width

Both ASDEX-Upgrade (AUG) data and the generalized HD (GHD) model showed that the scrape-off width broadens as the density/collisionality increases [1, 2]. A series of BOUT++ transport simulations are performed to study the physics of the scaling characteristics of the divertor heat flux width vs density/collisionality via a plasma density scan with either fixed pressure profile or fixed temperature profile inside separatrix. Additionally, the simulations show that even in the drift dominated regime, the divertor heat flux width can be broadened due to the transition of the SOL residence time from the parallel particle flow time to the enhanced parallel conduction time as the collisionality/density increases as posited in the GHD model. In addition, the heat flux width is found to be proportional to the square root of ion mass for low collisionality while it has a weakly dependence on ion mass for high collisionality. Furthermore, our simulations show that as the density increases, the radial electric field (E r ) well shallows, which potentially weakens E r × B flow shear stabilization of turbulence at high density.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Theory of stellar coronae - An interpretation of X-ray emission from non-degenerate stellar sources

It is shown that the acoustic wave heating theory of stellar coronae explains neither Einstein Observatory coronae data nor previous UV and X-ray observations of the sun and other stars, on the evidence of data implying that magnetic fields, stellar rotation rates and convection zone parameters figure in the determination of coronal heating. Einstein Observatory results suggest that O-type star coronae are heated by the interaction of turbulent stellar winds with slowly-decaying primordial magnetic fields or by radiative instabilities in the flow. The apparent absence of coronae in Ap stars is due to the stability of atmospheres in which even weak convection is suppressed by the strong field. Dynamo action is implicated in some normal A-type stars and in F- and later-type dwarfs. Coronal characteristics of dMe and dM stars, close binaries, and K- and M-type giants are also considered.

Linsky, J. L.↗

Feasibility study of optical boundary layer transition detection method

A high sensitivity differential interferometer was developed to locate the region where the boundary layer flow undergoes transition from laminar to turbulent. Two laboratory experimental configurations were used to evaluate the performance of the interferometer: open shear layer, and low speed wind tunnel turbulent spot configuration. In each experiment, small temperature fluctuations were introduced as the signal source. Simultaneous cold wire measurements were compared with the interferometer data. The comparison shows that the interferometer is sensitive to very weak phase variations in the order of 0.001 the laser wavelength. An attempt to detect boundary layer transition over a flat plate at NASA-Langley Unitary Supersonic Wind Tunnel using the interferometer system was performed. The phase variations during boundary layer transition in the supersonic wind tunnel were beyond the minimum signal-to-noise level of the instrument.

Azzazy, M.↗

Laminar-Turbulent Transition Upstream of the Entropy-Layer Swallowing Location in Hypersonic Boundary Layers

Numerical and experimental studies have demonstrated that modal growth of planar Mack modes is responsible for laminar-turbulent transition on sharp cones at hypersonic speeds. However, the physical mechanisms that lead to transition onset upstream of the entropy-layer swallowing location over sufficiently blunt geometries are not well understood as yet. Modal amplification is too weak or nonexistent to initiate transition at moderate-to-large bluntness values. Nonmodal analysis shows that, with increasing nose bluntness, both planar and oblique traveling disturbances that peak within the entropy layer experience appreciable energy amplification. However, because of the relatively weak signature of the nonmodal traveling disturbances within the boundary-layer region, the route to transition onset subsequent to the nonmodal growth remains unclear. Thus, nonlinear parabolized stability equations (NPSE) and direct numerical simulations (DNS) have been used to investigate the potential transition mechanisms over a 7-degree blunt cone that was tested in the AFRL Mach-6 high-Reynoldsnumber facility. Computations are performed to separately follow the nonlinear development of two classes of inflow disturbances, namely, a pair of oblique traveling waves with equal but opposite angles with respect to the mean flow direction and a planar traveling wave. Results in both cases show an excellent agreement between the NPSE and DNS predictions, establishing that the NPSE is an accurate and efficient technique for predicting the nonlinear development for these particular nonmodal traveling disturbances. Computations reveal that the oblique mode interactions lead to the generation of stationary streaks inside the boundary layer that, in turn, facilitate the growth of a subharmonic sinuous disturbance. For relatively modest amplitudes of the inflow disturbance, the oblique-mode breakdown can lead to transition at the measured location of transition onset during the experiment. On the other hand, the nonlinear development of a planar traveling wave leads to the formation of inclined structures just above the boundary-layer edge and these structures are strongly reminiscent of the transitional events observed during blunt cone experiments by using schlieren flow visualizations.

Entropy↗

Direct simulation of compressible turbulence

Several direct simulations of 3-D homogeneous, compressible turbulence are presented with emphasis on the differences with incompressible turbulent simulations. A fully spectral collocation algorithm, periodic in all directions coupled with a 3rd order Runge-Kutta time discretization scheme is sufficient to produce well-resolved flows at Taylor Reynolds numbers below 40 on grids of 128x128x128. A Helmholtz decomposition of velocity is useful to differentiate between the purely compressible effects and those effects solely due to vorticity production. In the context of homogeneous flows, this decomposition in unique. Time-dependent energy and dissipation spectra of the compressible and solenoidal velocity components indicate the presence of localized small scale structures. These structures are strongly a function of the initial conditions. Researchers concentrate on a regime characterized by very small fluctuating Mach numbers Ma (on the order of 0.03) and density and temperature fluctuations much greater than sq Ma. This leads to a state in which more than 70 percent of the kinetic energy is contained in the so-called compressible component of the velocity. Furthermore, these conditions lead to the formation of curved weak shocks (or shocklets) which travel at approximately the sound speed across the physical domain. Various terms in the vorticity and divergence of velocity production equations are plotted versus time to gain some understanding of how small scales are actually formed. Possible links with Burger turbulence are examined. To visualize better the dynamics of the flow, new graphic visualization techniques have been developed. The 3-D structure of the shocks are visualized with the help of volume rendering algorithms developed in-house. A combination of stereographic projection and animation greatly increase the number of visual cues necessary to properly interpret the complex flow.

Zang, T. A.↗

Numerical Studies of Turbulent Supersonic Plane-Channel Flows

This two-part talk begins with an overview of the `turbulence problem’, the computational strategies used to address it, and the role of direct numerical simulation (DNS). The second part presents DNS results to determine the turbulent Prandtl number Pr t above cold (isothermal) and hot (adiabatic) walls in a family of low-supersonic channel flows. A range of mean temperature/density variations, corresponding to effective/edge Mach numbers between 1.1 to 2.2, and wall-variable-based Reynolds number Reτw from 73 to 3800, is considered. The adiabatic condition is a new feature of special interest. The value of Pr t away from the wall approaches 0.85 above both the isothermal and adiabatic walls. The variations of the near-wall Pr t profiles in both the present and previous, passive-scalar simulations collapse as a function of the semilocal y∗wall scaling proposed by Huang et al. (1995), with only a weak dependence on R eτw . This leads to a rather simple proposal for a model of heat transfer, attached to an eddy-viscosity model.

Direct numerical simulation↗

The impact of magnetic fields on momentum transport and saturation of shear-flow instability by stable modes

The Kelvin–Helmholtz (KH) instability of a shear layer with an initially uniform magnetic field in the direction of flow is studied in the framework of 2D incompressible magnetohydrodynamics with finite resistivity and viscosity using direct numerical simulations. The shear layer evolves freely, with no external forcing, and thus broadens in time as turbulent stresses transport momentum across it. As with hydrodynamic KH, the instability here features a conjugate stable mode for every unstable mode in the absence of dissipation. Stable modes are shown to transport momentum up its gradient, shrinking the layer width whenever they exceed unstable modes in amplitude. In simulations with weak magnetic fields, the linear instability is minimally affected by the field, but enhanced small-scale fluctuations relative to the hydrodynamic case are observed. These enhanced fluctuations coincide with increased energy dissipation and faster layer broadening, with these features more pronounced in simulations with stronger fields. These trends result from the magnetic field reducing the effects of stable modes relative to the transfer of energy to small scales. As field strength increases, stable modes become less excited, thus transporting less momentum against its gradient. Furthermore, the energy that would otherwise transfer back to the driving shear because of the stable modes is instead allowed to cascade to small scales, where it is lost to dissipation. Approximations of the turbulent state in terms of a reduced set of modes are explored. While the Reynolds stress is well-described using just two modes per wavenumber at large scales, the Maxwell stress is not.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Performance of a quantitative jet stream turbulence forecasting technique - The Specific CAT Risk (SCATR) index

The SCATR index is showing considerable promise as a tool for operational, nonconvective, high altitude turbulence forecasting. The results of nine case studies are summarized. Attention is then focused on one of these cases for which the data base is most complete. In light of these case studies, the strengths and weaknesses of the current formulation and suggested improvements are discussed. The SCATR index has applications both as a diagnostic tool and in operational forecasting.

Keller, J. L.↗

Black hole accretion disks - Electrodynamic coupling of accretion-disk coronae and the partitioning of soft and hard X-ray emission

It is demonstrated that the observed large ratio of hard to soft X-ray emission and the bimodel behavior of black hole accreting X-ray sources such as Cygnus X-1 can be described in terms of a magnetically structured accretion disk corona which is electrodynamically coupled to the disk turbulent motions while the disk is thermodynamically coupled to the corona as described by a feedback parameter delta. The observed ratio of hard to soft X-ray emission is independent of the disk thickness, and weakly dependent of the disk parameter alpha relating the disk viscous stresses to the total pressure. Observed values of the luminosity ratio point towards strong differences of the feedback of the low state compared to the high state, in the sense that low state means small feedback (delta less than 0.2) and high state means strong feedback delta of about 0.5.

Kuperus, M.↗

High sensitivity boundary layer transition detector

A high sensitivity differential interferometer has been developed to locate the region where the boundary layer flow changes from laminar to turbulent. Two experimental configurations have been used to evaluate the performance of the interferometer, open shear layer configuration and wind tunnel turbulent spot configuration. In each experiment small temperature fluctuations were introduced as the signal source. Simultaneous cold wire measurements have been compared with the interferometer data. The comparison shows that the interferometer is sensitive to very weak phase variations in the order of .001 the laser wavelength.

Azzazy, M.↗

Experimental and numerical study of swept ramp injection into a supersonic flowfield

Time-averaged measurements of pressure, temperature, velocity, and injectant mole fraction are presented using the planar laser-induced iodine fluorescence technique in the complex three-dimensional compressible flowfield around a swept ramp fuel injector. Within the range of thermodynamic conditions present in the test case studied, the technique's accuracy is estimated to be 4% for pressure, temperature, and velocity and 3% for injectant mole fraction. Comparisons with numerical simulations using the SPARK three-dimensional Navier-Stokes computer code with an algebraic turbulence model are made at the centerplane of the flowfield as well as on three crossflow planes downstream of the injector. Calculations and measurements are in good agreement throughout the flowfield, with deviations on the order of 5%; however, in specific regions, such as in the base of the ramp, deviations are larger. A weak asymmetry in the incoming flowfield appears to be amplified by boundary-layer separation occurring when the ramp-generated shock reflects off the tunnel walls. Ramp-generated vortices are weaker in the calculated results due to the effects of numerical viscosity in the vortex cores. This leads to less turning and mixing of the jet plume than observed in the experiments. The rate of decay of the maximum injectant mole fraction with streamwise distance is greater for the present ramp injection scheme than for previously measured transverse injection schemes. In this recirculation region at the base of the injector, laminar calculations show better agreement with the measurements than turbulent calculations.

Donohue, James M.↗

Lagrangian large eddy simulations via physics-informed machine learning

High-Reynolds number homogeneous isotropic turbulence (HIT) is fully described within the Navier–Stokes (NS) equations, which are notoriously difficult to solve numerically. Engineers, interested primarily in describing turbulence at a reduced range of resolved scales, have designed heuristics, known as large eddy simulation (LES). LES is described in terms of the temporally evolving Eulerian velocity field defined over a spatial grid with the mean-spacing correspondent to the resolved scale. This classic Eulerian LES depends on assumptions about effects of subgrid scales on the resolved scales. Here, we take an alternative approach and design LES heuristics stated in terms of Lagrangian particles moving with the flow. Our Lagrangian LES, thus L-LES, is described by equations generalizing the weakly compressible smoothed particle hydrodynamics formulation with extended parametric and functional freedom, which is then resolved via Machine Learning training on Lagrangian data from direct numerical simulations of the NS equations. The L-LES model includes physics-informed parameterization and functional form, by combining physics-based parameters and physics-inspired Neural Networks to describe the evolution of turbulence within the resolved range of scales. The subgrid-scale contributions are modeled separately with physical constraints to account for the effects from unresolved scales. We build the resulting model under the differentiable programming framework to facilitate efficient training. We experiment with loss functions of different types, including physics-informed ones accounting for statistics of Lagrangian particles. We show that our L-LES model is capable of reproducing Eulerian and unique Lagrangian turbulence structures and statistics over a range of turbulent Mach numbers.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

High-sensitivity density fluctuation detector

A high-sensitivity differential interferometer has been developed to detect small density fluctuations over an optical path length of the order of the boundary layer thickness near transition. Two experimental configurations have been used to evaluate the performance of the interferometer: an open shear-layer configuration and a wind-tunnel turbulent spot configuration. In each experiment small temperature fluctuations were introduced as the signal source. Simultaneous cold-wire measurements have been compared with the interferometer data. The comparison shows that the interferometer is sensitive to very weak phase variations of the order of 0.001 of the laser wavelength.

Azzazy, M.↗