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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
NASA NTRS1983

Radio and X-ray observations of the radio halo source in A1367

Radio and X-ray observations of the 'radio halo' in A1367 are presented. The radio spectral index distribution is derived combining data at 610 and 1415 MHz taken with the Westerbork Synthesis Radio Telescope (WSRT). A high-resolution Very Large Array (VLA) map allows us to discriminate between the extended cluster emission and a triple radio source associated with a background stellar object (most probably a QSO) which appears projected onto it. X-ray observations taken with the high-resolution imager on board the Einstein satellite reveal an X-ray source associated with this background object. A second weak X-ray enhancement is found projected onto the extended radio component. Inverse Compton scattering of the synchrotron-emitting particles provides a possible interpretation of such a source. The presence of large-scale magnetic fields and particle acceleration due to the turbulent motion of galaxies through the intergalactic medium is discussed in relationship to the observed peculiarities of the galaxies found in this region.

Gavazzi, G.↗
NASA NTRS1988

Scaling laws for homogeneous turbulent shear flows in a rotating frame

The scaling properties of plane homogeneous turbulent shear flows in a rotating frame are examined mathematically by a direct analysis of the Navier-Stokes equations. It is proved that two such shear flows are dynamically similar if and only if their initial dimensionless energy spectrum E star (k star, 0), initial dimensionless shear rate SK sub 0/epsilon sub 0, initial Reynolds number K squared sub 0/nu epsilon sub 0, and the ration of the rotation rate to the shear rate omega/S are identical. Consequently, if universal equilibrium states exist, at high Reynolds numbers, they will only depend on the single parameter omega/S. The commonly assumed dependence of such equilibrium states on omega/S through the Richardson number Ri=-2(omega/S)(1-2 omega/S) is proven to be inconsistent with the full Navier-Stokes equations and to constitute no more than a weak approximation. To be more specific, Richardson number similarity is shown to only rigorously apply to certain low-order truncations of the Navier-Stokes equations (i.e., to certain second-order closure models) wherein closure is achieved at the second-moment level by assuming that the higher-order moments are a small perturbation of their isotropic states. The physical dependence of rotating turbulent shear flows on omega/S is discussed in detail along with the implications for turbulence modeling.

Speziale, Charles G.↗
NASA NTRS1989

Scaling laws for homogeneous turbulent shear flows in a rotating frame

The scaling properties of plane homogeneous turbulent shear flows in a rotating frame are examined mathematically by a direct analysis of the Navier-Stokes equations. It is proved that two such shear flows are dynamically similar if and only if their initial dimensionless energy spectrum E star (k star, 0), initial dimensionless shear rate SK sub 0/epsilon sub 0, initial Reynolds number K squared sub 0/nu epsilon sub 0, and the ration of the rotation rate to the shear rate omega/S are identical. Consequently, if universal equilibrium states exist, at high Reynolds numbers, they will only depend on the single parameter omega/S. The commonly assumed dependence of such equilibrium states on omega/S through the Richardson number Ri = -2(omega/S)(1-2 omega/S) is proven to be inconsistent with the full Navier-Stokes equations and to constitute no more than a weak approximation. To be more specific, Richardson number similarity is shown to only rigorously apply to certain low-order truncations of the Navier-Stokes equations (i.e., to certain second-order closure models) wherein closure is achieved at the second-moment level by assuming that the higher-order moments are a small perturbation of their isotropic states. The physical dependence of rotating turbulent shear flows on omega/S is discussed in detail along with the implications for turbulence modeling.

Speziale, Charles G.↗
DOE OSTI2022

High Performance Equilibrium Solvers for Integrated Magnetic Fusion Simulations

With the funding provided by this award, we developed numerical codes for the study of magnetically confined plasmas for fusion applications. Accordingly, our work can be divided into two separate categories: 1) the design and analysis of novel numerical methods providing high accuracy and high efficiency; 2) the study of the equilibrium and stability of magnetically confined plasmas with some of these numerical codes, as well as the study of the nature of the turbulent behavior which may arise in the presence of instabilities. We first developed new numerical schemes based on integral equation methods for the computation of steady-state magnetic configurations in fusion experiments, providing high accuracy for the magnetic field and its derivatives, which are required for stability and turbulence calculations. We employed different integral formulations depending on the application of interest: axisymmetric or non-axisymmetric equilibria, force-free or magnetohydrodynamic equilibria, fixed-boundary equilibria or free-boundary equilibria. While efficient, these methods do not yet apply to plasma boundaries which are not smooth, a situation which is fairly common in magnetic confinement experiments. To address this temporary weakness, we also constructed a new steady-state solver based on the Hybridizable Discontinuous Galerkin (HDG) method, which provides full geometric flexibility. In addition to these numerical tools focused on steady-states, we also contributed to the improvement of the speed and accuracy of codes simulating the plasma dynamics of fusion plasmas, by developing a novel velocity space representation for the efficient solution of kinetic equations, which most accurately describe the time evolution of hot plasmas in fusion experiments. Using the tools discussed above, we studied several questions pertaining to the equilibrium and stability of magnetically confined plasmas. In particular, we derived a new simple model for axisymmetric devices called tokamaks, to predict how elongated a fusion plasma can be before it becomes unstable and collapses. We also looked at the effect of the shape of the outer plasma surface on key properties of the steady-state magnetic configurations, and how these properties impact turbulence in fusion plasmas, and the corresponding transport of momentum. Likewise, we studied the role of large localized flows on the steady-state magnetic configurations, and how they may influence plasma stability and turbulence. Non-axisymmetric steady-state magnetic configurations are inherently more complex than axisymmetric steady-state configurations, and the subject of ongoing controversies regarding the regularity of the equations determining such steady-states, and their solutions. Implementing an existing NYU code in a new geometry, we studied the nature of the singularity of the solutions observed in the code, and methods to eliminate them. Our main conclusion is that by appropriately tailoring the plasma boundary, it is possible to eliminate the singularities otherwise appearing in our simulations, and to obtain steady-states which appear to be smooth. To gain further insights on incompletely understood turbulence phenomena, we proposed a new reduced model capturing most of these phenomena, which is simple enough to not require expensive numerical simulations on massive supercomputers to investigate them. We demonstrated the strong similarity between our simulations and published results obtained from computationally expensive simulations, and plan to rely on our reduced model to identify the key mechanisms determining the evolution and strength turbulent driven transport in fusion plasmas. Finally, we proposed a new framework for tokamak reactor design studies, enabling us to consider the relative merits of steady-state versus pulsed fusion reactors. We found that pulsed fusion reactors may benefit most from recent advances in magnet technology, and the availability of very high field magnets. As such, they may become more desirable than steady-state tokamak reactors for cost efficient electricity generation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Self-sustaining sound in collisionless, high- β plasma

Using analytical theory and hybrid-kinetic numerical simulations, we demonstrate that, in a collisionless plasma, long-wavelength ion-acoustic waves (IAWs) with amplitudes $δn/n_0≳2/β$ (where $β\gg1$ is the ratio of thermal to magnetic pressure) generate sufficient pressure anisotropy to destabilize the plasma to firehose and mirror instabilities. These kinetic instabilities grow rapidly to reduce the pressure anisotropy by pitch-angle scattering and trapping particles, respectively, thereby impeding the maintenance of Landau resonances that enable such waves’ otherwise potent collisionless damping. The result is wave dynamics that evince a weakly collisional plasma: the ion distribution function is near-Maxwellian, the field-parallel flow of heat resembles its Braginskii form (except in regions where large-amplitude magnetic mirrors strongly suppress particle transport), and the relations between various thermodynamic quantities are more ‘fluid-like’ than kinetic. A nonlinear fluctuation–dissipation relation for self-sustaining IAWs is obtained by solving a plasma-kinetic Langevin problem, which demonstrates suppressed damping, enhanced fluctuation levels and weakly collisional thermodynamics when IAWs with $δn/n_0≳2/β$ are stochastically driven. We investigate how our results depend upon the scale separation between the wavelength of the IAW and the Larmor radius of the ions, and discuss briefly their implications for our understanding of turbulence and transport in the solar wind and the intracluster medium of galaxy clusters.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Pressure–strain interaction. I. On compression, deformation, and implications for Pi-D

The pressure–strain interaction describes the rate per unit volume that energy is converted between bulk flow and thermal energy in neutral fluids or plasmas. The term has been written as a sum of the pressure dilatation and the collisionless analog of viscous heating referred to as Pi–D , which isolates the power density due to compressible and incompressible effects, respectively. It has been shown that Pi–D can be negative, which makes its identification as collisionless viscous heating troubling. We argue that an alternate decomposition of pressure–strain interaction can be useful for interpreting the underlying physics. Since Pi–D contains both normal deformation and shear deformation, we propose grouping the normal deformation with the pressure dilatation to describe the power density due to converging/diverging flows, with the balance describing the power density purely due to shear deformation. We then develop a kinetic theory interpretation of compression, normal deformation, and shear deformation. We use the results to determine the physical mechanisms that can make Pi–D negative. We argue that both decompositions can be useful for the study of energy conversion in weakly collisional or collisionless fluids and plasmas, and implications are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1977

Investigations related to the inviscid-viscous interaction in transonic flows about finite 3-D wings

Inviscid small-disturbance theory has been shown to predict three-dimensional transonic flows about finite wings reasonably well as long as viscous effects are negligible. In order to include these effects, the inviscid small-disturbance solution of Bailey and Ballhaus (1975) has been combined with a finite-difference solution for Prandtl's boundary-layer equations. This solution employs the conditionally stable Krause (1968) scheme, implicit in the direction normal to the wall, to cope with the domain-of-dependence problem that arises for reverse cross flow. To be consistent with the inviscid-flow solution, the boundary layer is computed in the representative wing planform plane which is transformed into rectangular shape in the computational domain. The flow has been assumed turbulent, and a scalar eddy-viscosity model is adopted. The interaction between inviscid and viscous flow is modeled with the help of the displacement surface which is added to the geometric wing shape. Sample distributions of displacement thickness for swept wings are presented for weak and strong interaction cases.

Kordulla, W.↗
NASA NTRS1985

Micron CO in the eclipse spectrum of Epsilon Aur

Two micrometer spectrum of epsilon Aur was observed at high resolution (lambda/delta lambda approx. 40000) with the KPNO 4m telescope and FTS approximately every 100 days during the recent eclipse. CO delta v=2 lines appeared in the spectrum only after mid-totality. No CO features were present JD 2445478 but weak absorption lines were present on 2445592. The CO strengthened in the spectrum and remained present until at least 2445812. CO was again absent after forth contact when the spectrum was observed on 2445912. The CO appears to originate in a cool (approx. 1000 K), turbulent region. Excitation temperatures, velocities, and column densities is presented.

Hinkle, K. H.↗
NASA NTRS1989

The use of dc-probe measurements in the middle atmosphere program

The middle atmosphere is weakly ionized, collision dominated plasma. Rocket-borne dc probes have proven effective in measuring conductivity and electron and ion concentrations in this plasma, and in some cases electron temperatures and ion composition were observed. During MAP, dc probes were also used successfully to study turbulence in the middle atmosphere, using the plasma as a passive tracer for the dynamical processes in the nonionized gas. Basic principles are reviewed of the dc-probe technique and examples are presented of results obtained during MAP.

Thrane, E. V.↗
NASA NTRS2000

Realistic Solar Surface Convection Simulations

We perform essentially parameter free simulations with realistic physics of convection near the solar surface. We summarize the physics that is included and compare the simulation results with observations. Excellent agreement is obtained for the depth of the convection zone, the p-mode frequencies, the p-mode excitation rate, the distribution of the emergent continuum intensity, and the profiles of weak photospheric lines. We describe how solar convection is nonlocal. It is driven from a thin surface thermal boundary layer where radiative cooling produces low entropy gas which forms the cores of the downdrafts in which most of the buoyancy work occurs. We show that turbulence and vorticity are mostly confined to the intergranular lanes and underlying downdrafts. Finally, we illustrate our current work on magneto-convection.

Stein, Robert F.↗
DOE OSTI2021

Dissipative magnetohydrodynamics for nonresistive relativistic plasmas: An implicit second-order flux-conservative formulation with stiff relaxation

Based on a 14-moment closure for nonresistive (general-) relativistic viscous plasmas, here we describe a new numerical scheme that is able to handle all first-order dissipative effects (heat conduction, bulk and shear viscosities), as well the anisotropies induced by the presence of magnetic fields. The latter is parametrized in terms of a thermal gyrofrequency or, equivalently, a thermal Larmor radius and allows to correctly capture the thermal Hall effect. By solving an extended Israel-Stewart-like system for the dissipative quantities that enforces algebraic constraints via stiff-relaxation, we are able to cast all first-order dissipative terms in flux-divergence form. This allows us to apply traditional high-resolution shock capturing methods to the equations, making the system suitable for the numerical study of highly turbulent flows. We present several numerical tests to assess the robustness of our numerical scheme in flat spacetime. The 14-moment closure can seamlessly interpolate between the highly collisional limit found in neutron star mergers, and the highly anisotropic limit of relativistic Braginskii magnetohydrodynamics appropriate for weakly collisional plasmas in black-hole accretion problems. We believe that this new formulation and numerical scheme will be useful for a broad class of relativistic magnetized flows.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1974

Magnetic fields in the sun

The observed properties of solar magnetic fields are reviewed, with particular reference to the complexities imposed on the field by motions of the highly conducting gas. Turbulent interactions between gas and field lead to heating or cooling of the gas according as the field energy density is less or greater than the maximum kinetic energy density in the convection zone. The field strength above which cooling sets in is 700 G. A weak solar dipole field may be primeval, but dynamo action is also important in generating new flux. The dynamo is probably not confined to the convection zone, but extends throughout most of the volume of the sun. Planetary tides appear to play a role in driving the dynamo.

Mullan, D. J.↗
NASA NTRS1993

An experimental study of the sources of fluctuating pressure loads beneath swept shock/boundary-layer interactions

An experimental research program providing basic knowledge and establishing a database on the fluctuating pressure loads produced on aerodynamic surfaces beneath three dimensional shock wave/boundary layer interactions is described. Such loads constitute a fundamental problem of critical concern to future supersonic and hypersonic flight vehicles. A turbulent boundary layer on a flat plate is subjected to interactions with swept planar shock waves generated by sharp fins at angle of attack. Fin angles from 10 to 20 deg at freestream Mach numbers of 3 and 4 produce a variety of interaction strengths from weak to very strong. Miniature Kulite pressure transducers flush-mounted in the flat plate are used to measure interaction-induced wall pressure fluctuations. The distributions of properties of the pressure fluctuations, such as their ring levels, amplitude distributions, and power spectra, are also determined. Measurements were made for the first time in the aft regions of these interactions, revealing fluctuating pressure levels as high as 160 dB. These fluctuations are dominated by low frequency (0-5 kHz) signals. The maximum ring levels in the interactions show an increasing trend with increasing interaction strength. On the other hand, the maximum ring levels in the forward portion of the interactions decrease linearly with increasing interaction sweep back. These ring pressure distributions and spectra are correlated with the features of the interaction flowfield. The unsteadiness of the off-surface flowfield is studied using a new, non-intrusive technique based on the shadow graph method. The results indicate that the entire lambda-shock structure generated by the interaction undergoes relatively low-frequency oscillations. Some regions where particularly strong fluctuations are generated were identified. Fluctuating pressure measurements are also made along the line of symmetry of an axisymmetric jet impinging upon a flat plate at an angle. This flow was chosen as a simple analog to the impinging jet region found in the rear portion of the shock wave/boundary layer interactions under study. It is found that a sharp peak in ring pressure level exists at or near the mean stagnation point. It is suggested that the phenomena responsible for this peak may be active in the swept interactions as well, and may cause the extremely high fluctuating pressures observed in the impinging jet region in the present experimental program.

Settles, G. S.↗
NASA NTRS2011

Frame-Transfer Gating Raman Spectroscopy for Time-Resolved Multiscalar Combustion Diagnostics

Accurate experimental measurement of spatially and temporally resolved variations in chemical composition (species concentrations) and temperature in turbulent flames is vital for characterizing the complex phenomena occurring in most practical combustion systems. These diagnostic measurements are called multiscalar because they are capable of acquiring multiple scalar quantities simultaneously. Multiscalar diagnostics also play a critical role in the area of computational code validation. In order to improve the design of combustion devices, computational codes for modeling turbulent combustion are often used to speed up and optimize the development process. The experimental validation of these codes is a critical step in accepting their predictions for engine performance in the absence of cost-prohibitive testing. One of the most critical aspects of setting up a time-resolved stimulated Raman scattering (SRS) diagnostic system is the temporal optical gating scheme. A short optical gate is necessary in order for weak SRS signals to be detected with a good signal- to-noise ratio (SNR) in the presence of strong background optical emissions. This time-synchronized optical gating is a classical problem even to other spectroscopic techniques such as laser-induced fluorescence (LIF) or laser-induced breakdown spectroscopy (LIBS). Traditionally, experimenters have had basically two options for gating: (1) an electronic means of gating using an image intensifier before the charge-coupled-device (CCD), or (2) a mechanical optical shutter (a rotary chopper/mechanical shutter combination). A new diagnostic technology has been developed at the NASA Glenn Research Center that utilizes a frame-transfer CCD sensor, in conjunction with a pulsed laser and multiplex optical fiber collection, to realize time-resolved Raman spectroscopy of turbulent flames that is free from optical background noise (interference). The technology permits not only shorter temporal optical gating (down to <1 s, in principle), but also higher optical throughput, thus resulting in a substantial increase in measurement SNR.

Nguyen, Quang-Viet↗
DOE OSTI2022

Impact of pedestal density gradient and collisionality on ELM dynamics

BOUT++ turbulence simulations are conducted to capture the underlying physics of small ELM characteristics achieved by increasing separatrix density via controlling strike points from vertical to horizontal divertor plates for three EAST discharges. BOUT++ linear simulations show that the most unstable modes change from high-n ideal ballooning modes to intermediate-n peeling–ballooning modes and eventually to peeling–ballooning stable plasmas in the pedestal. Nonlinear simulations show that the fluctuation is saturated at a high level for the lowest separatrix density. The ELM size decreases with increasing separatrix density, until the fraction of this energy lost during the ELM crash becomes less than 1% of the pedestal stored energy, leading to small ELMs. Simulations indicate that small ELMs can be triggered either by the marginally peeling–ballooning instability near the peak pressure gradient position inside the pedestal or by a local instability in the pedestal foot with a larger separatrix density gradient. The pedestal collisionality scan for type-I ELMs with steep pedestal density gradient shows that both linear growth rate and ELM size decrease with increasing collisionality. On the contrary, the pedestal collisionality and pedestal density width scan with a weak pedestal density gradient indicate small ELMs can either be triggered by a high-n ballooning mode or by a low-n peeling mode in a low collisionality region 0.04–0.1. We report the simulations indicate the weaker the linear unstable modes near marginal stability with small linear growth rate, the lower nonlinearly saturated fluctuation intensity and the smaller turbulence spreading from the linear unstable zone to stable zone in the nonlinear saturation phase, leading to small ELMs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

On FIRE mode in KSTAR

We report on the status of Fast Ion Regulated Enhancement (FIRE) mode experiments in the Korea Superconducting Tokamak Advanced Research. This regime is being developed for high-performance, steady-state operation which features a stationary ion internal transport barrier, enabling a central ion temperature approaching 10 keV to be sustained for up to 50 s, without the need for delicate profile control and with no significant impurity accumulation. As its key novelty lies in the significant contribution of fast ions that stabilize core turbulence, the regime has been named FIRE mode. To achieve this regime, neutral beam injection is applied at moderate power levels near the L–H power threshold, while maintaining low plasma density to avoid the L–H transition. The scenario is typically established in diverted magnetic configurations. The core features of FIRE mode were investigated through power balance analysis and fluctuation measurements, revealing a clear transport bifurcation in the ion channel. At the plasma edge, FIRE mode occasionally exhibits I-mode characteristics, particularly in unfavorable magnetic null configurations with q 95 ∼ 4, including the presence of weakly coherent modes. In terms of MHD activity, sawtooth oscillations are observed but appear to be stabilized during the high-performance phase. Fast ion-driven Alfvénic eigenmodes (AEs), indicated by strong frequency chirping near 200 kHz, are also observed. Additionally, lower-frequency MHD activities, distinct from the fast ion-driven AEs, are present and have some influence on plasma performance. The enhancement of core confinement is primarily attributed to fast ion effects, which were evaluated with respect to dilution, alpha stabilization, and resonant interactions with turbulence using gyrokinetic analyses. Among these, the dilution effect was found to be the most dominant. The characteristics of FIRE mode were compared with those of other hot ion plasma scenarios, such as supershot and hot ion mode. While they share many similarities, FIRE mode is distinguished by the accessibility to conditions with T i ≈ T e , presence of I-mode edge features, and its long-duration sustainment. Predictive simulations of FIRE mode were performed using integrated transport modeling with TRIASSIC, employing the TGLF anomalous transport model. These simulations confirmed the critical role of fast ions in achieving this regime. The future prospect of FIRE mode for application in fusion reactors is discussed, with an emphasis on possibly extending the regime to higher density operation.

FIRE mode↗
NASA NTRS2008

The Aeroacoustics of Turbulent Flows

Aerodynamic noise prediction has been an important and challenging research area since James Lighthill first introduced his Acoustic Analogy Approach over fifty years ago. This talk attempts to provide a unified framework for the subsequent theoretical developments in this field. It assumes that there is no single approach that is optimal in all situations and uses the framework as a basis for discussing the strengths weaknesses of the various approaches to this topic. But the emphasis here will be on the important problem of predicting the noise from high speed air jets. Specific results will presented for round jets in the 0.5 to 1.4 Mach number range and compared with experimental data taken on the Glenn SHAR rig. It is demonstrated that nonparallel mean flow effects play an important role in predicting the noise at the supersonic Mach numbers. The results explain the failure of previous attempts based on the parallel flow Lilley model (which has served as the foundation for most jet noise analyses during past two decades).

Goldstein, M. E.↗
NASA NTRS1974

Magnetic fields in the sun

The observed properties of solar magnetic fields are reviewed, with particular reference to the complexities imposed on the field by motions of the highly conducting gas. Turbulent interactions between gas and field lead to heating or cooling of the gas according to whether the field energy density is less or greater than the maximum kinetic energy density in the convection zone. The field strength above which cooling sets in is 700 gauss. A weak solar dipole field may be primeval, but dynamo action is also important in generating new flux. The dynamo is probably not confined to the convection zone, but extends throughout most of the volume of the sun. Planetary tides appear to play a role in driving the dynamo.

Mullan, D. J.↗