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 415 records · Page 23

Dead Zone Accretion Flows in Protostellar Disks

Planets form inside protostellar disks in a dead zone where the electrical resistivity of the gas is too high for magnetic forces to drive turbulence. We show that much of the dead zone nevertheless is active and flows toward the star while smooth, large-scale magnetic fields transfer the orbital angular momentum radially outward. Stellar X-ray and radionuclide ionization sustain a weak coupling of the dead zone gas to the magnetic fields, despite the rapid recombination of free charges on dust grains. Net radial magnetic fields are generated in the magnetorotational turbulence in the electrically conducting top and bottom surface layers of the disk, and reach the midplane by ohmic diffusion. A toroidal component to the fields is produced near the midplane by the orbital shear. The process is similar to the magnetization of the solar tachocline. The result is a laminar, magnetically driven accretion flow in the region where the planets form.

solar system↗

Interactive phenomena in supersonic jet mixing problems. II - Numerical studies

Applications of the numerical techniques discussed in Part I of this article are presented. Fundamental wave/shear-layer interaction calculations are performed exhibiting the waves produced by high-speed turbulent mixing processes and the interactions of an expansion fan and shock wave with a developing shear layer. Calculations for an underexpanded supersonic jet into a supersonic external stream at different velocity ratios are presented which delineate betweeen weakly interactive and strongly interactive processes in the jet near field. Calculations for the multiple-cell shock structure of an underexpanded jet into still air exhibit the damping of wave intensities by the turbulence and compare quite favorably with available data.

Dash, S. M.↗

Pressure–strain interaction. II. Decomposition in magnetic field-aligned coordinates

In weakly collisional and collisionless magnetized plasmas, the pressure–strain interaction describes the rate of conversion between bulk flow and thermal energy density. In this study, we derive an analytical expression for the pressure–strain interaction in a coordinate system with an axis aligned with the local magnetic field. The result is eight groups of terms corresponding to different physical mechanisms that can contribute to the pressure–strain interaction. We provide a physical description of each term. Here, the results are immediately of interest to weakly collisional and collisionless magnetized plasmas and the fundamental processes that happen therein, including magnetic reconnection, magnetized plasma turbulence, and collisionless shocks. The terms in the field-aligned coordinate decomposition are likely accessible to measurement with satellite observations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Magnetic shear effect on plasma transport at T e / T i ~ 1 through electron cyclotron heating in DIII-D plasmas

The effect of magnetic shear on plasma transport for an electron to ion temperature ratio ( T e / T i ) near unity has been explored in DIII-D utilizing electron cyclotron heating (ECH). Previous reports showed that significant confinement degradation occurred at T e / T i ~ 1 in positive shear (PS) plasmas in DIII-D, whereas reduced confinement degradation was observed in negative central shear (NCS) plasmas. In this study, plasma transport in weak magnetic shear (WS) plasmas with ECH is investigated and compared with that in NCS and PS plasmas. Here the magnetic shears ($\hat s$) are $\hat s$ > 0.5, ~0 and <-0.1 in the core region (ρ~ 0.3–0.4) of PS, WS and NCS plasmas, respectively, and flat or negative inside ρ~ 0.4 in the WS and NCS plasmas. Weak magnetic shear is found to be effective in minimizing degradation of ion thermal confinement as T e / T i increases through ECH application, and an improved confinement factor of H 98y2 ~ 1.2 is maintained, similar to NCS plasmas. At T e / T i ~ 1, the ion thermal diffusivity around an internal transport barrier decreases when changing the magnetic shear from positive to weak or negative shear. Also, reduced local particle and momentum transport was indicated by steeper density and toroidal rotation profiles in the weak and negative shear regimes. Linear gyrokinetic simulations predict little change in growth rates of low-k turbulence with ECH application in the WS and NCS plasmas, which is consistent with the transport and profile analyses.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Transonic Flow Past a Symmetrical Airfoil-- Inviscid and Turbulent Flow Properties

Flowfield measurements are presented for a symmetrical NACA 64A010 airfoil section at transonic conditions. Measurements were obtained for three angles of attack with the freestream Mach number fixed at O.S. The cases studied included a weak shock-wave/boundary-layer interaction, an interaction of medium strength with mild separation, and an interaction of sufficient strength to produce a shock-induced stall situation. Two nonintrusive optical techniques, lower velocimetry and holographic interferometry, were used to characterize the flows. The results include Mach number contours and flow angle distributions in the inviscid flow regions, and turbulent flow properties, including the turbulent Reynolds stresses, of the upper surface viscous layers, and of the near-wake. The turbulent flow measurements reveal that the turbulence fluctuations attain equilibrium with the local mean flow much faster than previously expected.

Johnson D. A.↗

Multifrequency studies of bright radio supernova remnants. 2: W49B

We report radio observations of the supernova remnant W49B using the Very Large Array (VLA) of the NRAO at 90, 20, and 6 cm. This work continues the study of the properties of young, bright supernova remnants (SNRs) begun with 3C 391 (Moffett & Reynolds 1994). Here we present high-resolution images of total intensity, polarization, and spectral index of W49B. In large-scale morphology it is basically a box-shaped remnant expanding into an apparently inhomogeneous medium, but we also find small-scale structures, arcs and filaments, which lie interior to the outer shell in projection. The actual spatial position of the filaments is unknown, though we suspect that they are in the remnant interior rather than on the front or back face of the blast wave. In addition, their distribution in two dimensions suggests the projection of a helical structure. The shell edge, sharply rising in brightness for at least 3/4 of its circumference, still remains unresolved at our highest resolution (4 sec). This for our assumed distance of 8 kpc, the width of the region in which the emission appears is less than 0.16 pc, indicating short mean free paths for shock-accelerated electrons and high levels of MHD turbulence presumably causing the scattering. We find no polarized flux at 90 or 20 cm, with 3 sigma upper limits in polarized intensity of 36 mJy at 90 cm (45 sec resolution) and 7.2 mJy at 20 cm (5 sec resolution), or 22 micro-Jy/sq arcsec and 370 micro-jy/sq arcsec, respectively, for any emission extended on those scales. Polarized flux is present at 6 cm, but at a very low mean polarized fraction (total polarized flux divided by total flux) of 0.44% +/- 0.06%, which, as for 3C 391, is much lower than typical for bright SNRs at this frequency. The morphology in polarized intensity is poorly correlated with that in total intensity. We see excursions in the polarized fraction up to at least 10% in a few locations, but even where polarization is seen, typical levels are a few percent. Tangled or disordered magnetic fields in the emitting region of the radio shell may be responsible for depolarizing the radio synchrotron radiation, but some form of internal Faraday depolarization may also occur. We estimate the foreground Faraday rotation measure to be about -450 rad/sq m, similar to that found for 3C 391, which is also the Galactic plane and just about as distant. Spectral index images created from the total intensity images show that the spectral index across W49B is constant to within about Delta alpha approximately 0.1 in bright regions. This result weakly supports a common origin of the radio-emitting electrons, as in the blast wave, rather than in inhomogeneous turbulent regions of differing properties due to the stochastic (second-order Fermi) acceleration process. Variations at the level of Delta alpha approximately 0.1 are seen, but their significance is doubtful. New observations at 90 cm, using experimental three-dimensional imaging technology, may improve on this limit.

Moffett, David A.↗

Turbulent Compressible Convection with Rotation: Flow Structure and Evolution - Part 1

The effects of Coriolis forces on compressible convection are studied using three-dimensional numerical simulations carried out within a local modified f-plane model. The physics is simplified by considering a perfect gas occupying a rectilinear domain placed tangentially to a rotating sphere at various latitudes, through which a destabilizing heat flux is driven. The resulting convection is considered for a range of Rayleigh, Taylor, and Prandtl (and thus Rossby) numbers, evaluating conditions where the influence of rotation is both weak and strong. Given the computational demands of these high-resolution simulations, the parameter space is explored sparsely to ascertain the differences between laminar and turbulent rotating convection. The first paper in this series examines the effects of rotation on the flow structure within the convection, its evolution, and some consequences for mixing. Subsequent papers consider the large-scale mean shear flows that are generated by the convection, and the effects of rotation on the convective energetics and transport properties. It is found here that the structure of rotating turbulent convection is similar to earlier nonrotating studies, with a laminar, cellular surface network disguising a fully turbulent interior punctuated by vertically coherent structures. However, the temporal signature of the surface flows is modified by inertial motions to yield new cellular evolution patterns and an overall increase in the mobility of the network. The turbulent convection contains vortex tubes of many scales, including large-scale coherent structures spanning the full vertical extent of the domain involving multiple density scale heights. Remarkably, such structures align with the rotation vector via the influence of Coriolis forces on turbulent motions, in contrast with the zonal tilting of streamlines found in laminar flows. Such novel turbulent mechanisms alter the correlations which drive mean shearing flows and affect the convective transport properties. In contrast to this large-scale anisotropy, small-scale vortex tubes at greater depths are randomly orientated by the rotational mixing of momentum, leading to an increased degree of isotropy on the medium to small scales of motion there. Rotation also influences the thermodynamic mixing properties of the convection. In particular, interaction of the larger coherent vortices causes a loss of correlation between the vertical velocity and the temperature leaving a mean stratification which is not isentropic.

Brummell, Nicholas H.↗

Tethered balloon-borne profile measurements of atmospheric properties in the cloudy atmospheric boundary layer over the Arctic sea ice during MOSAiC: Overview and first results

The tethered balloon-borne measurement system BELUGA (Balloon-bornE moduLar Utility for profilinG the lower Atmosphere) was deployed over the Arctic sea ice for 4 weeks in summer 2020 as part of the Multidisciplinary drifting Observatory for the Study of Arctic Climate expedition. Using BELUGA, vertical profiles of dynamic, thermodynamic, aerosol particle, cloud, radiation, and turbulence properties were measured from the ground up to a height of 1,500 m. BELUGA was operated during an anomalously warm period with frequent liquid water clouds and variable sea ice conditions. Three case studies of liquid water phase, single-layer clouds observed on 3 days (July 13, 23, and 24, 2020) are discussed to show the potential of the collected data set to comprehensively investigate cloud properties determining cloud evolution in the inner Arctic over sea ice. Simulated back-trajectories show that the observed clouds have evolved within 3 different air masses (“aged Arctic,” “advected over sea ice,” and “advected over open ocean”), which left distinct fingerprints in the cloud properties. Strong cloud top radiative cooling rates agree with simulated results of previous studies. The weak warming at cloud base is mostly driven by the vertical temperature profile between the surface and cloud base. In-cloud turbulence induced by the cloud top cooling was similar in strength compared to former studies. From the extent of the mixing layer, it is speculated that the overall cloud cooling is stronger and thus faster in the warm oceanic air mass. Larger aerosol particle number concentrations and larger sizes were observed in the air mass advected over the sea ice and in the air mass advected over the open ocean.

Arctic↗

Development of Renormalization Group Analysis of Turbulence

The renormalization group (RG) procedure for nonlinear, dissipative systems is now quite standard, and its applications to the problem of hydrodynamic turbulence are becoming well known. In summary, the RG method isolates self similar behavior and provides a systematic procedure to describe scale invariant dynamics in terms of large scale variables only. The parameterization of the small scales in a self consistent manner has important implications for sub-grid modeling. This paper develops the homogeneous, isotropic turbulence and addresses the meaning and consequence of epsilon-expansion. The theory is then extended to include a weak mean flow and application of the RG method to a sequence of models is shown to converge to the Navier-Stokes equations.

L M Smith↗

Energy Transfer and Triadic Interactions in Compressible Turbulence

Using a two-point closure theory, the Eddy-Damped-Quasi-Normal-Markovian (EDQNM) approximation, we have investigated the energy transfer process and triadic interactions of compressible turbulence. In order to analyze the compressible mode directly, the Helmholtz decomposition is used. The following issues were addressed: (1) What is the mechanism of energy exchange between the solenoidal and compressible modes, and (2) Is there an energy cascade in the compressible energy transfer process? It is concluded that the compressible energy is transferred locally from the solenoidal part to the compressible part. It is also found that there is an energy cascade of the compressible mode for high turbulent Mach number (M(sub t) greater than or equal to 0.5). Since we assume that the compressibility is weak, the magnitude of the compressible (radiative or cascade) transfer is much smaller than that of solenoidal cascade. These results are further confirmed by studying the triadic energy transfer function, the most fundamental building block of the energy transfer.

Bataille, F.↗

Evaluation of Turbulence-Model Performance in Jet Flows

The importance of reducing jet noise in both commercial and military aircraft applications has made jet acoustics a significant area of research. A technique for jet noise prediction commonly employed in practice is the MGB approach, based on the Lighthill acoustic analogy. This technique requires as aerodynamic input mean flow quantities and turbulence quantities like the kinetic energy and the dissipation. The purpose of the present paper is to assess existing capabilities for predicting these aerodynamic inputs. Two modern Navier-Stokes flow solvers, coupled with several modern turbulence models, are evaluated by comparison with experiment for their ability to predict mean flow properties in a supersonic jet plume. Potential weaknesses are identified for further investigation. Another comparison with similar intent is discussed by Barber et al. The ultimate goal of this research is to develop a reliable flow solver applicable to the low-noise, propulsion-efficient, nozzle exhaust systems being developed in NASA focused programs. These programs address a broad range of complex nozzle geometries operating in high temperature, compressible, flows. Seiner et al. previously discussed the jet configuration examined here. This convergent-divergent nozzle with an exit diameter of 3.6 inches was designed for an exhaust Mach number of 2.0 and a total temperature of 1680 F. The acoustic and aerodynamic data reported by Seiner et al. covered a range of jet total temperatures from 104 F to 2200 F at the fully-expanded nozzle pressure ratio. The aerodynamic data included centerline mean velocity and total temperature profiles. Computations were performed independently with two computational fluid dynamics (CFD) codes, ISAAC and PAB3D. Turbulence models employed include the k-epsilon model, the Gatski-Speziale algebraic-stress model and the Girimaji model, with and without the Sarkar compressibility correction. Centerline values of mean velocity and mean temperature are compared with experimental data.

Woodruff, S. L.↗

Assessment of Turbulent Shock-Boundary Layer Interaction Computations Using the OVERFLOW Code

The performance of two popular turbulence models, the Spalart-Allmaras model and Menter s SST model, and one relatively new model, Olsen & Coakley s Lag model, are evaluated using the OVERFLOWcode. Turbulent shock-boundary layer interaction predictions are evaluated with three different experimental datasets: a series of 2D compression ramps at Mach 2.87, a series of 2D compression ramps at Mach 2.94, and an axisymmetric coneflare at Mach 11. The experimental datasets include flows with no separation, moderate separation, and significant separation, and use several different experimental measurement techniques (including laser doppler velocimetry (LDV), pitot-probe measurement, inclined hot-wire probe measurement, preston tube skin friction measurement, and surface pressure measurement). Additionally, the OVERFLOW solutions are compared to the solutions of a second CFD code, DPLR. The predictions for weak shock-boundary layer interactions are in reasonable agreement with the experimental data. For strong shock-boundary layer interactions, all of the turbulence models overpredict the separation size and fail to predict the correct skin friction recovery distribution. In most cases, surface pressure predictions show too much upstream influence, however including the tunnel side-wall boundary layers in the computation improves the separation predictions.

Oliver, A. B.↗

Transonic flow about a two-dimensional airfoil inviscid and turbulent flow properties.

Flow-field measurements are presented for a symmetrical NACA 64A010 airfoil section at transonic conditions. Measurements were obtained for three angles of attack with the free-stream Mach number fixed at 0.8. The cases studied included a weak shock wave/boundary layer interaction, an interaction of medium strength with mild separation, and an interaction of sufficient strength to produce a shock-induced stall situation. Two nonintrusive optical techniques, laser velocimetry and holographic interferometry, were used to characterize the flows. The results include Mach number contours and flow angle distributions in the inviscid flow regions, and turbulent flow properties, including the turbulent Reynolds stresses, of the upper surface viscous layers, and of the near-wake. The turbulent flow measurements reveal that the turbulence fluctuations attain equilibrium with the local mean flow much faster than previously expected.

Johnson, D. A.↗

Impact of Cloud-Base Turbulence on CCN Activation: CCN Distribution

Abstract Following our previous investigation of the turbulence impact on cloud-base single-size CCN activation, this study considers a similar problem assuming CCN size distribution obtained from field measurements. The total CNN concentration is taken as either 200 cm −3 to represent clean conditions, or as 2000 cm −3 to represent polluted conditions. CCN is assumed to be sodium chloride. The CCN activation in the rising nonturbulent adiabatic parcel is contrasted with the activation within a rising adiabatic parcel filled with inertial-range homogeneous isotropic turbulence. The turbulent parcel of 64 3 m 3 and the turbulent kinetic energy dissipation rate of 10 −3 m −2 s −3 are used in most of the simulations. Results for a range of mean parcel ascent rates, between 0.125 and 8 m s −1 , are discussed. Overall, the adiabatic turbulent parcel simulations show results consistent with the adiabatic nonturbulent parcel, with higher activated CCN concentrations for stronger parcel ascent rates. The key difference is a blurriness of the separation between dry CCN size bins featuring activated and nonactivated (haze) CCN, especially for weak mean ascent rates. The blurriness comes from CCN getting activated and subsequently deactivated in the fluctuating supersaturation field, instead of all becoming cloud droplets above the cloud base. This leads to significantly larger spectral widths in turbulent parcel simulations compared to the nonturbulent parcel when activation is completed. Modeling results are discussed in the context of the impact of turbulent fluctuations on CCN activation documented in laboratory experiments using the Pi chamber.

Meteorology & Atmospheric Sciences↗

Stability of superthermal strahl electrons in the solar wind

We present a kinetic stability analysis of the solar wind electron distribution function consisting of the Maxwellian core and the magnetic-field aligned strahl, a superthermal electron beam propagating away from the sun. We use an electron strahl distribution function obtained as a solution of a weakly collisional drift-kinetic equation, representative of a strahl affected by Coulomb collisions but unadulterated by possible broadening from turbulence. This distribution function is essentially non-Maxwellian and varies with the heliospheric distance. The stability analysis is performed with the Vlasov–Maxwell linear solver leopard. We find that depending on the heliospheric distance, the core-strahl electron distribution becomes unstable with respect to sunward-propagating kinetic-Alfvén, magnetosonic, and whistler modes, in a broad range of propagation angles. The wavenumbers of the unstable modes are close to the ion inertial scales, and the radial distances at which the instabilities first appear are on the order of 1 au. However, we have not detected any instabilities driven by resonant wave interactions with the superthermal strahl electrons. Instead, the observed instabilities are triggered by a relative drift between the electron and ion cores necessary to maintain zero electric current in the solar wind frame (ion frame). Contrary to strahl distributions modelled by shifted Maxwellians, the electron strahl obtained as a solution of the kinetic equation is stable. Our results are consistent with the previous studies based on a more restricted solution for the electron strahl.

79 ASTRONOMY AND ASTROPHYSICS↗

Closures for multicomponent reacting flows based on dispersion analysis

This work presents algebraic closure models associated with advective transport and nonlinear reactions in a Reynolds-averaged Navier-Stokes context for a system of species subject to binary reactions and transport by advection and diffusion. Expanding upon analysis originally developed for non-reactive transport in the context of Taylor dispersion of scalars, this work extends the modified gradient diffusion model explicated by Peters [N. Peters, Turbulent Combustion, Cambridge Monographs on Mechanics (Cambridge University Press, Cambridge, 2000)] and based on work by Corrsin [S. Corrsin, The reactant concentration spectrum in turbulent mixing with a first-order reaction, J. Fluid Mech. 11, 407 (1961)] beyond single-component transport phenomena and involving nonlinear reactions. The presented model forms, from this weakly nonlinear extension of the original dispersion theory, lead to an analytic expression for the eddy diffusivity matrix that explicitly captures the influence of the reaction kinetics on the closure operators. Furthermore, we demonstrate that the derived model form directly translates between flow topologies through a priori and a posteriori testing of a binary species system subject to homogeneous isotropic turbulence. Using two- and three-dimensional direct numerical simulations involving laminar and turbulent flows, it is shown that this framework improves prediction of mean quantities compared to previous results. Lastly, the presented model form, collapses to the earlier gradient diffusion and its modified version derived by Corrsin in the limits of nonreactive species and linear reactions, respectively.

42 ENGINEERING↗

Inverse transonic airfoil design including viscous interaction

A numerical technique was developed for the analysis of specified transonic airfoils or for the design of airfoils having a prescribed pressure distribution, including the effect of weak viscous interaction. The method uses the full potential equation, a stretched Cartesian coordinate system, and the Nash-MacDonald turbulent boundary layer method. Comparisons with experimental data for typical transonic airfoils show excellent agreement. An example shows the application of the method to design a thick aft-cambered airfoil, and the effects of viscous interaction on its performance are discussed.

Carlson, L. A.↗

Numerical simulation of viscous transonic airfoil flows

Numerical simulations of transonic airfoil flows using the Reynolds-averaged Navier-Stokes equations and various turbulence models are presented and compared with experimental data. Three different airfoils were investigated under varying flow conditions ranging from subcritical unseparated flows to supercritical separated flows. The turbulence models investigated consisted of three zero-equation models and one two-equation model. For unseparated flows involving weak viscous-inviscid interactions, the four models were comparable in their agreement with experiment. For separated flows involving strong viscous-inviscid interactions, the nonequilibrium zero-equation model of Johnson and King gave the best overall agreement with experiment.

Coakley, Thomas J.↗