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.

SourceResearch connectionFrom yearThrough year

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

714 records · Page 40

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Density wavenumber spectrum measurements, synthetic diagnostic development, and tests of quasilinear turbulence modeling in the core of electron-heated DIII-D H-mode plasmas

Abstract Measurements of the turbulent density wavenumber spectrum, δ n ˆ e ( k ⊥ ) , using the Doppler Back-Scattering (DBS) diagnostic are reported from DIII-D H-mode plasmas with electron cyclotron heating as the only auxiliary heating method. These electron-heated plasmas have low collisionality, ν e ∗ < 1 , T e / T i > 1 , and zero injected torque—a regime expected to be relevant for future fusion devices. We probe density fluctuations in the core ( ρ ≈ 0.7) over a broad wavenumber range, 0.5 ⩽ k ⊥ ⩽ 16 cm −1 ( 0.1 ⩽ k ⊥ ρ s ⩽ 5 ), to characterize plasma instabilities and compare with theoretical predictions. We present a novel synthetic DBS diagnostic to relate the back-scattered power spectrum, P s ( k ⊥ ) —which is directly measured by DBS—to the underlying electron density fluctuation spectrum, δ n ˆ e ( k ⊥ ) . The synthetic DBS P s ( k ⊥ ) spectrum is calculated by combining the SCOTTY beam-tracing code with a model δ n ˆ e ( k ⊥ ) predicted either analytically or numerically. In this work we use the quasi-linear code Trapped Gyro-Landau Fluid (TGLF) to approximate the δ n ˆ e ( k ⊥ ) spectrum. We find that TGLF, using the experimental profiles, is capable of closely reproducing the DBS measurements. Both the DBS measurements and the TGLF-DBS synthetic diagnostic show a wavenumber spectrum with variable decay. The measurements show weak decay ( k −0.6 ) for k < 3.5 cm −1 , with k −2.6 at intermediate- k ( 3.5 ⩽ k ⩽ 8.5 cm −1 ), and rapid decay ( k −9.4 ) for k > 8.5 cm −1 . Scans of physics parameters using TGLF suggest that the normalized ∇ T e scale-length, R / L T e , is an important factor for distinguishing microturbulence regimes in these plasmas. A combination of DBS observations and TGLF simulations indicate that fluctuations remain peaked at ITG-scales (low k ) while R / L T e -driven TEM/ETG-type modes (intermediate/high k ) are marginally sub-dominant.

synthetic diagnostics↗

NASA NTRS1997

Shapes of Buoyant and Nonbuoyant Methane Laminar Jet Diffusion Flames

Laminar gas jet diffusion flames represent a fundamental combustion configuration. Their study has contributed to numerous advances in combustion, including the development of analytical and computational combustion tools. Laminar jet flames are pertinent also to turbulent flames by use of the laminar flamelet concept. Investigations into the shapes of noncoflowing microgravity laminar jet diffusion flames have primarily been pursued in the NASA Lewis 2.2-second drop tower, by Cochran and coworkers and by Bahadori and coworkers. These studies were generally conducted at atmospheric pressure; they involved soot-containing flames and reported luminosity lengths and widths instead of the flame-sheet dimensions which are of Greater value to theory evaluation and development. The seminal model of laminar diffusion flames is that of Burke and Schumann, who solved the conservation of momentum equation for a jet flame in a coflowing ambient by assuming the velocity of fuel, oxidizer and products to be constant throughout. Roper and coworkers improved upon this model by allowing for axial variations of velocity and found flame shape to be independent of coflow velocity. Roper's suggestion that flame height should be independent of gravity level is not supported by past or present observations. Other models have been presented by Klajn and Oppenheim, Markstein and De Ris, Villermaux and Durox, and Li et al. The common result of all these models (except in the buoyant regime) is that flame height is proportional to fuel mass flowrate, with flame width proving much more difficult to predict. Most existing flame models have been compared with shapes of flames containing soot, which is known to obscure the weak blue emission of flame sheets. The present work involves measurements of laminar gas jet diffusion flame shapes. Flame images have been obtained for buoyant and nonbuoyant methane flames burning in quiescent air at various fuel flow-rates, burner diameters and ambient pressures. Soot concentrations were minimized by selecting conditions at low flowrates and low ambient pressures; this allows identification of actual flame sheets associated with blue emissions of CH and CO2. The present modeling effort follows that of Roper and is useful in explaining many of the trends observed.

Sunderland, Peter B.↗

NASA NTRS1982

Propagation of wide bandwidth signals through strongly turbulent ionized media

Analytic and numerical techniques are presented which directly address the problem of propagation of wide bandwidth signals through random ionized media. This work is applicable to the problems of satellite communication and space based radar observation through a disturbed ionospheric propagation channel that would result from a high altitude chemical release or nuclear detonation. An analytic solution is presented for the two-position, two-frequency mutual coherence function for spherical wave propagation in the strong scattering limit. This solution is used to derive simple expressions for the power impulse response function and to determine its relationship to the important parameters of decorrelation distance and coherence bandwidth which describe the disturbed propagation channel. Results for mean time delay and time delay jitter are presented and compared to direct simulation results and to other theoretical calculations. A numerical/analytical solution to the parabolic wave equation is presented in the form of a multiple phase-screen (MPS) propagation simulation. In this very general solution technique, the random medium is divided into a finite number of layers. The field fluctuation through each layer is obtained by replacing the layer by a centrally located thin phase-screen, whose statistical properties are determined from the statistics of the electron-density irregularities. The waveform then propagates from phase-screen to phase-screen via an exact solution to the Fresnel-Kirchnoff equation. For wide bandwidth waveforms, numerical solutions are obtained at a number of discrete frequencies centered about the carrier and then time-domain solutions are obtained by Fourier transform techniques. Detailed results are presented for a number of interesting cases including propagation of a 20 MHz bandwidth signal through a finite barium cloud at a carrier frequency of 100 MHz. One of the major uses of the MPS propagation simulation has been to provide realizations of the received signal after propagation through a disturbed channel. The MPS simulation obtains a general solution of the parabolic wave equation under both weak and strong scattering conditions. A second technique to directly obtain realizations of wide bandwidth waveforms is presented here. This technique is limited to the case of strong scattering but requires only a fraction of the computer resources needed for MPS signal generation. Detailed comparisons of the two signal generation techniques are presented.

Random Media↗

NASA NTRS1994

Experimental studies of shock-wave/wall-jet interaction in hypersonic flow

Experimental studies have been conducted to examine slot film cooling effectiveness and the interaction between the cooling film and an incident planar shock wave in turbulent hypersonic flow. The experimental studies were conducted in the 48-inch shock tunnel at Calspan at a freestream Mach number of close to 6.4 and at a Reynolds number of 35 x 10(exp 6) based on the length of the model at the injection point. The Mach 2.3 planar wall jet was generated from 40 transverse nozzles (with heights of both 0.080 inch and 0.120 inch), producing a film that extended the full width of the model. The nozzles were operated at pressures and velocities close to matching the freestream, as well as at conditions where the nozzle flows were over- and under-expanded. A two-dimensional shock generator was used to generate oblique shocks that deflected the flow through total turnings of 11, 16, and 21 degrees; the flows impinged downstream of the nozzle exits. Detailed measurements of heat transfer and pressure were made both ahead and downstream of the injection station, with the greatest concentration of measurements in the regions of shock-wave/boundary layer interaction. The major objectives of these experimental studies were to explore the effectiveness of film cooling in the presence of regions of shock-wave/boundary layer interaction and, more specifically, to determine how boundary layer separation and the large recompression heating rates were modified by film cooling. Detailed distributions of heat transfer and pressure were obtained in the incident shock/wall-jet interaction region for a series of shock strengths and impingement positions for each of the two nozzle heights. Measurements were also made to examine the effects of nozzle lip thickness on cooling effectiveness. The major conclusion from these studies was that the effect of the cooling film could be readily dispersed by relatively weak incident shocks, so the peak heating in the recompression region was not significantly reduced by even the largest levels of film cooling. For the case studies in the absence of film cooling, the interaction regions were unseparated. However, adding film cooling resulted in regions of boundary layer separation induced in the film cooling layer -- the size of which regions first increased and then decreased with increased film cooling. Surprisingly, the size of the separated regions and the magnitude of the recompression heating were not strongly influenced by the thickness of the cooling film, nor by the point of shock impingement relative to the exit plane of the nozzles. The lip thickness was found to have little effect on cooling effectiveness. Measurements with and in the absence of shock interaction were compared with the results of earlier experimental studies and correlated in terms of the major parameters controlling these flows.

Holden, Michael S.↗

NASA NTRS1994

Experimental studies of shock-wave/wall-jet interaction in hypersonic flow, part A

Experimental studies have been conducted to examine slot film cooling effectiveness and the interaction between the cooling film and an incident planar shock wave in turbulent hypersonic flow. The experimental studies were conducted in the 48-inch shock tunnel at Calspan at a freestream Mach number of close to 6.4 and at a Reynolds number of 35 x 10(exp 6) based on the length of the model at the injection point. The Mach 2.3 planar wall jet was generated from 40 transverse nozzles (with heights of both 0.080 inch and 0.120 inch), producing a film that extended the full width of the model. The nozzles were operated at pressures and velocities close to matching the freestream, as well as at conditions where the nozzle flows were over- and under-expanded. A two-dimensional shock generator was used to generate oblique shocks that deflected the flow through total turnings of 11, 16, and 21 degrees; the flows impinged downstream of the nozzle exits. Detailed measurements of heat transfer and pressure were made both ahead and downstream of the injection station, with the greatest concentration of measurements in the regions of shock-wave/boundary layer interaction. The major objectives of these experimental studies were to explore the effectiveness of film cooling in the presence of regions of shock-wave/boundary layer interaction and, more specifically, to determine how boundary layer separation and the large recompression heating rates were modified by film cooling. Detailed distributions of heat transfer and pressure were obtained in the incident-shock/wall-jet interaction region for a series of shock strengths and impingement positions for each of the two nozzle heights. Measurements were also made to examine the effects of nozzle lip thickness on cooling effectiveness. The major conclusion from these studies was that the effect of the cooling film could be readily dispersed by relatively weak incident shocks, so the peak heating in the recompression region was not significantly reduced by even the largest levels of film cooling. For the case studies in the absence of film cooling, the interaction regions were unseparated. However, adding film cooling resulted in regions of boundary layer separation induced in the film cooling layer, the size of which regions first increased and then decreased with increased film cooling. Surprisingly, the size of the separated regions and the magnitude of the recompression heating were not strongly influenced by the thickness of the cooling film, nor by the point of shock impingement relative to the exit plane of the nozzles. The lip thickness was found to have little effect on cooling effectiveness. Measurements with and in the absence of shock interaction were compared with the results of earlier experimental studies and correlated in terms of the major parameters controlling these flows.

Holden, Michael S.↗

NASA NTRS2019

Nonmodal Growth of TravelingWaves on Blunt Cones at Hypersonic Speeds

The existing database of transition measurements in hypersonic ground facilities has established that, as the nosetip bluntness is increased, the onset of boundary layer transition over a circular cone at zero angle of attack shifts downstream. However, this trend is reversed at sufficiently large values of the nose Reynolds number, so that the transition onset location eventually moves upstream with a further increase in nose-tip bluntness. Because modal amplification is too weak to initiate transition at moderate-to-large bluntness values, nonmodal growth has been investigated as the potential basis for a physics-based model for the frustum transition. The present analysis investigates the nonmodal growth of traveling disturbances initiated within the nose-tip vicinity that peak within the entropy layer. Results show that, with increasing nose bluntness, both planar and oblique traveling disturbances experience appreciable energy amplification up to successively higher frequencies. For moderately blunt cones, the initial nonmmodal growth is followed by a partial decay that is more than overcome by an eventual, modal growth as Mack-mode waves. For larger bluntness values, the Mack-mode waves are not amplified anywhere upstream of the experimentally measured transition location, but the traveling modes still undergo a significant amount of nonmodal growth. This finding does not provide a definitive link between optimal growth and the onset of transition, but it is qualitatively consistent with the experimental observations that frustum transition in the absence of sufficient Mack-mode amplification implies a double peak in disturbance amplification and the appearance of transitional events above the boundary-layer edge.

Hypersonic speed↗

DOE OSTI2021

On closures for reduced order models - A spectrum of first-principle to machine-learned avenues

For over a century, reduced order models (ROMs) have been a fundamental discipline of theoretical fluid mechanics. Early examples include Galerkin models inspired by the Orr–Sommerfeld stability equation and numerous vortex models, of which the von Karman vortex street is one of the most prominent. Subsequent ROMs typically relied on first principles, like mathematical Galerkin models, weakly nonlinear stability theory, and two- and three-dimensional vortex models pioneered the data-driven proper orthogonal decomposition (POD) modeling. In early POD modeling, available data were used to build an optimal basis, which was then utilized in a classical Galerkin procedure to construct the ROM, but data have made a profound impact on ROMs beyond the Galerkin expansion. In this paper, we take a modest step and illustrate the impact of data-driven modeling on one significant ROM area. Specifically, we focus on ROM closures, which are correction terms that are added to the classical ROMs in order to model the effect of the discarded ROM modes in under-resolved simulations. Through simple examples, we illustrate the main modeling principles used to construct the classical ROMs, motivate and introduce modern ROM closures, and show how data-driven modeling, artificial intelligence, and machine learning have changed the standard ROM methodology over the last two decades. Finally, we outline our vision on how the state-of-the-art data-driven modeling can continue to reshape the field of reduced order modeling.

42 ENGINEERING↗

DOE OSTI2021

Self-consistent core-pedestal ITER scenario modeling

The purpose of these integrated ITER simulations is to identify dependencies that can impact the performance of ITER. The dependence of fusion power production, temperature and density pedestals, and core profiles on varying magnetic-q, edge density fueling strength, neoclassical transport, magnetic field strength, and alpha heating, are examined. It is noted that ion modes usually dominate the interior, while electron modes dominate the barrier. The electron mode without flow shear is substantially stronger. The q increase tends to cause more kinetic ballooning modes at the edge and the q decrease means that the peeling is getting worse, but the ballooning is getting better. The large edge q tends to provide hollow density while small edge q gives normal density profile. This is consistent with the observation that the density of particles is inversely proportional to the density of the plasma current. The rise in the source of the edge particle increases the density of the edge and the reaction rate close to the edge increases temperature fluxes, resulting in weak pedestal barriers. When neoclassical transport is turned off or the strength of B-field is increased, the pronounced edge barriers are identified. In addition, a distinct effect on the magnetohydrodynamic (MHD) modes is observed due to the increased B-field. The interesting aspect is that the barrier collapses about q = 2 as the stability of MHD has decreased to replace the barrier with a continuous slope but with roughly the same fusion power. This would actually be good for a reactor. The slope of the H-mode pedestal is found to be reduced due to the alpha heating. It is expected that the problem of edge-localized modes causing damage to the first wall will be eased. The effects of circular, elongated, and general geometry on plasma profiles are also compared and contrasted. The circular geometry has a hollow density profile and no edge or internal thermal transport barrier (ETB or ITB). However, both ETBs and ITBs can be found in the elongated geometry. In general geometry, there is a weaker ETB, comparable ITBs, but a higher edge particle transport barrier than in elongated geometry. Higher density near the edge, on the other hand, causes more wall erosion, so elongated geometry might be the best option.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS2022

Nonlinear Evolution of Instabilities in a Laminar Separation Bubble at a Hypersonic Mach Number

The development of both convective stationary perturbation as well as global instabilities in the vicinity of a laminar separation bubble above an axisymmetric compression corner in a hypersonic flow is investigated using numerical simulations. The flow configuration of primary interest corresponds to the cone-cylinder-flare model used in experimental measurements in the Boeing/AFOSR Mach-6 Quiet Tunnel at Purdue University. For a flare angle of 10 degrees and a unit Reynolds number of 11.5 x 10 6 m -1 , their surface flow visualizations identified the presence of streamwise elongated thermal streaks near the reattachment position and the dominant azimuthal spacing between the streaks was determined to be approximately 10 degrees (i.e., an azimuthal mode number of 36). Previous linear stability analyses predicted that the amplification characteristics of small amplitude, unsteady, convective instabilities within this flow were consistent with the surface pressure fluctuations measured in the experiment. However, their accompanying analysis of global instabilities had found the separation bubble to be weakly unstable at the 10 degrees flare angle, with the most unstable global mode corresponding to a stationary disturbance with an azimuthal wavenumber of m ≈ 5, which was well below the measured wavenumber of m = 36. Besides confirming the presence of the global instability at these flow conditions, the present numerical simulations quantify the details of the stationary equilibrium state associated with the supercritical bifurcation resulting from the nonlinear saturation of the unstable global mode. Although velocity perturbations associated with the saturated global mode are dominated by the fundamental spanwise wavelength associated with the linear global instability, the surface heat flux downstream of reattachment is dominated by m = 36 in agreement with experimental measurements. However, the shorter wavelength peak ( m = 36) in the azimuthal spectrum is not manifested as streaks with an azimuthal spacing of 10 degrees. An additional simulation of the combined evolution of both unsteady convective and global instabilities indicates that the nonlinear interactions between a broadband spectrum of oblique Mack’s first mode instabilities generate stationary vortex modes with a higher dominant wavenumber ( m = 60) as that observed in the experiments. Additional analyses are performed to predict the global mode behavior at other flare angles and how it differs from the convective centrifugal instability due to concave streamline curvature over the compression corner.

Hypersonic flow↗

NASA NTRS2000

Experiments on the Richtmyer-Meshkov Instability of Incompressible Fluids

Richtmyer-Meshkov (R-M) instability occurs when two different density fluids are impulsively accelerated in the direction normal to their nearly planar interface. The instability causes small perturbations on the interface to grow and possibly become turbulent given the proper initial conditions. R-M instability is similar to the Rayleigh-Taylor (R-T) instability, which is generated when the two fluids undergo a constant acceleration. R-M instability is a fundamental fluid instability that is important to fields ranging from astrophysics to high-speed combustion. For example, R-M instability is currently the limiting factor in achieving a net positive yield with inertial confinement fusion. The experiments described here utilize a novel technique that circumvents many of the experimental difficulties previously limiting the study of the R-M instability. A Plexiglas tank contains two unequal density liquids and is gently oscillated horizontally to produce a controlled initial fluid interface shape. The tank is mounted to a sled on a high speed, low friction linear rail system, constraining the main motion to the vertical direction. The sled is released from an initial height and falls vertically until it bounces off of a movable spring, imparting an impulsive acceleration in the upward direction. As the sled travels up and down the rails, the spring retracts out of the way, allowing the instability to evolve in free-fall until impacting a shock absorber at the end of the rails. The impulsive acceleration provided to the system is measured by a piezoelectric accelerometer mounted on the tank, and a capacitive accelerometer measures the low-level drag of the bearings. Planar Laser-Induced Fluorescence is used for flow visualization, which uses an Argon ion laser to illuminate the flow and a CCD camera, mounted to the sled, to capture images of the interface. This experimental study investigates the instability of an interface between incompressible, miscible liquids with an initial sinusoidal perturbation. The amplitude of the disturbance during the experiment is measured and compared to theory. The results show good agreement (within 10%) with linear stability theory up to nondimensional amplitude ka = 0.7 (wavenumber x amplitude). These results hold true for an initial ka (before acceleration) of -0.7 less than ka less than -0.06, while the linear theory was developed for absolute value of ka much less than 1. In addition, a third order weakly nonlinear perturbation theory is shown to be accurate for amplitudes as large as ka = 1.3, even though the interface becomes double-valued at ka = 1.1. As time progresses, the vorticity on the interface concentrates, and the interface spirals around the alternating sign vortex centers to form a mushroom pattern. At higher Reynolds Number (based on circulation), an instability of the vortex cores has been observed. While time limitations of the apparatus prevent determination of a critical Reynolds Number, the lowest Reynolds Number this vortex instability has been observed at is 5000.

Jacobs, J.↗

NASA NTRS2022

Revolutionizing Our Understanding of Particle Energization in Space Plasmas Using on-Board Wave-Particle Correlator Instrumentation

A leap forward in our understanding of particle energization in plasmas throughout the heliosphere is essential to answer longstanding questions in heliophysics, including the heating of the solar corona, acceleration of the solar wind, and energization of particles that lead to observable phenomena, such as the Earth’s aurora. The low densities and high temperatures of typical heliospheric environments lead to weakly collisional plasma conditions. Under these conditions, the energization of particles occurs primarily through collisionless interactions between the electromagnetic fields and the individual plasma particles with energies characteristic of a particular interaction. To understand how the plasma heating and particle acceleration impacts the macroscopic evolution of the heliosphere, impacting phenomena such as extreme space weather, it is critical to understand these collisionless wave-particle interactions on the characteristic ion and electron kinetic timescales. Such understanding requires high-cadence measurements of both the electromagnetic fields and the three-dimensional particle velocity distributions. Although existing instrument technology enables these measurements, a major challenge to maximize the scientific return from these measurements is the limited amount of data that can be transmitted to the ground due to telemetry constraints. A valuable, but underutilized, approach to overcome this limitation is to compute on-board correlations of the maximum-cadence field and particle measurements to improve the sampling time by several orders of magnitude. Here we review the fundamentals of the innovative field-particle correlation technique, present a formulation of the technique that can be implemented as an on-board wave-particle correlator, and estimate results that can be achieved with existing instrumental capabilities for particle velocity distribution measurements.

plasma heating↗

DOE OSTI2021

Idealized Large-Eddy Simulations of Stratocumulus Advecting over Cold Water. Part I: Boundary Layer Decoupling

Abstract We explore the decoupling physics of a stratocumulus-topped boundary layer (STBL) moving over cooler water, a situation mimicking warm-air advection (WADV). We simulate an initially well-mixed STBL over a doubly periodic domain with the sea surface temperature decreasing linearly over time using the System for Atmospheric Modeling large-eddy model. Due to the surface cooling, the STBL becomes increasingly stably stratified, manifested as a near-surface temperature inversion topped by a well-mixed cloud-containing layer. Unlike the stably stratified STBL in cold-air advection (CADV) that is characterized by cumulus coupling, the stratocumulus deck in the WADV is unambiguously decoupled from the sea surface, manifested as weakly negative buoyancy flux throughout the subcloud layer. Without the influxes of buoyancy from the surface, the convective circulation in the well-mixed cloud-containing layer is driven by cloud-top radiative cooling. In such a regime, the downdrafts propel the circulation, in contrast to that in CADV regime for which the cumulus updrafts play a more determinant role. Such a contrast in convection regime explains the difference in many aspects of the STBLs including the entrainment rate, cloud homogeneity, vertical exchanges of heat and moisture, and lifetime of the stratocumulus deck, with the last being subject to a more thorough investigation in Part II. Finally, we investigate under what conditions a secondary stratus near the surface (or fog) can form in the WADV. We found that weaker subsidence favors the formation of fog whereas a more rapid surface cooling rate does not. Significant Statement The low-lying blanketlike clouds, called stratocumulus (Sc), reflect much incoming sunlight, substantially modulating Earth’s temperature. While much is known about how the Sc evolves when it moves over warmer water, few studies examine the opposite situation of Sc moving over colder water. We used a high-resolution numerical model to simulate such a case. When moving over cold water, the Sc becomes unambiguously decoupled from the water surface, distinctive from its warm counterpart in which the Sc interacts with the water surface via intermittent cauliflower-like clouds called cumulus clouds. Such decoupling influences many aspects of the Sc–sea surface system, which combine to alter the ability of the Sc to reflect sunlight, thereby influencing the climate. This work laid the foundation for future work that quantifies the contribution of such a decoupled Sc regime to Earth’s radiative budget and climate change.

54 ENVIRONMENTAL SCIENCES↗