Growth of turbulent magnetic fields.
Weak random initial magnetic field evolution in highly conducting isotropically turbulent fluid, noting exact initial growth expression of magnetic energy spectrum
SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Weak random initial magnetic field evolution in highly conducting isotropically turbulent fluid, noting exact initial growth expression of magnetic energy spectrum
This is the final technical report for the project "Collaborative Research: Unveiling the Interplay Between Magnetic Reconnection and Turbulence — Theory and Modeling." This project investigated the interplay between two fundamental plasma processes — turbulence and magnetic reconnection — and the associated heating and particle acceleration. In particular, we focused on a new regime where the plasmoid instability mediates the turbulent energy cascade. In this new regime, reconnecting current sheets are disrupted by the growth of plasmoids/flux ropes on time scales shorter than the typical turnover times of eddies, speeding up the energy cascade and steepening the turbulence energy spectrum. This project addressed the following outstanding open questions: (1) What are the fundamental differences between plasmoid‐mediated turbulence cascade in three dimensions (3D) and two dimensions (2D)? (2) What are the essential features of plasmoid‐mediated turbulence in weakly collisional or collisionless regimes beyond resistive MHD? (3) How does the interplay of turbulence and reconnection regulate the energy release and dissipation in space and astrophysical plasmas? To address these questions, we have carried out three investigations: (1) Three-dimensional plasmoid-mediated turbulence energy cascade; (2) plasmoid-mediated reconnection and turbulence in three-dimensional Hall MHD; (3) reconnection and heating in coronal loops. This technical report summarizes the key findings of these investigations.
The general topic of atmospheric turbulence effects on sonic boom propagation is addressed with especial emphasis on taking proper and efficient account of the contributions of the portion oi the turbulence that is associated with extremely high wavenumber components. The recent work reported by Bart Lipkens in his doctoral thesis is reexamined to determine whether the good agreement between his measured rise times with the 1971 theory of the author is fortuitous. It is argued that Lipken's estimate of the distance to the first caustic was a gross overestimate because of the use of a sound speed correlation function shaped like a gaussian curve. In particular, it is argued that the expected distance to the first caustic varies with the kinematic viscosity nu and the energy epsilon dissipated per unit mass per unit time, and the sound speed c as : d(sub first caustic) = nu(exp 7/12) c(exp 2/3)/ epsilon(exp 5/12)(nu x epsilon/c(exp 4))(exp a), where the exponent a is greater than -7/12 and can be argued to be either O or 1/24. In any event, the surprising aspect of the relationship is that it actually goes to zero as the viscosity goes to zero with s held constant. It is argued that the apparent overabundance of caustics can be grossly reduced by a general computational and analytical perspective that partitions the turbulence into two parts, divided by a wavenumber k(sub c). Wavenumbers higher than kc correspond to small-scale turbulence, and the associated turbulence can be taken into account by a renormalization of the ambient sound speed so that the result has a small frequency dependence that results from a spatial averaging over of the smaller-scale turbulent fluctuations. Selection of k(sub c). can be made so large that only a very small number of caustics are encountered if one adopts the premise that the frequency dispersion of pulses is caused by that part of the turbulence spectrum which lies in the inertial range originally predicted by Kolmogoroff. The acoustic propagating wave's dispersion relation has the acoustic wavenumber being of the form k = (omega/c) + F(omega), where c is a spatially averaged sound speed and where, for mechanical turbulence, the extra term F(omega) must depend on only the angular frequency omega, the sound speed c, and the turbulent energy dissipation epsilon per unit fluid mass and per unit time. If the turbulence is weak, then the quantity F(omega) has to be of second order in the portions of the turbulent fluid velocity in the inertial range, so, following Kolmogoroff's reasoning, it must vary with epsilon as epsilon(exp 2/3). Simple dimensional analysis then reveals that F(omega) is K epsilon(exp 2/3) c(exp -7/3) omega(exp l/3), K being a universal dimensionless complex constant.
An experimental study of the generation of Tollmien-Schlichting waves and wave packets in a flat-plate boundary-layer by weak freestream turbulence has been conducted with the intent of clarifying receptivity mechanisms. Emphasis was placed upon the properties of such waves at stations as far forward as the minimum critical Reynolds number. It was found that alteration of the flow about the leading edge, due either to an asymmetry associated with lift, or due to a change of the fineness ratio of the leading edge, altered the T-S wave amplitude at early stations. The subsequent growth of the waves proceeded faster than expected according to certain stability theory results. Speculation regarding receptivity mechanisms is made.
We investigate the spectral properties of the electromagnetic fluctuations of sub-ion scale turbulence in weakly collisional, low-beta plasmas using a two-field isothermal gyrofluid model. The numerical results strongly support a description of the turbulence as a critically balanced Kolmogorov-like cascade of kinetic Alfvén wave fluctuations, as amended by previous studies to include intermittency effects. The measured universal index of the energy spectra from systems with different flux-unfreezing mechanisms excludes the role of tearing mediation in determining the spectra. The fluctuations remain isotropic in the plane perpendicular to the strong background magnetic fields as they cascade to smaller scales, which explains the absence of tearing mediation. The calculation of high-order, multipoint structure functions of magnetic fluctuations suggests that the intermittent structures have a quasi-2D, sheet-type morphology. In conclusion, these results are useful for explaining recent observations of the spectrum and structure of magnetic and density fluctuations in the solar wind at sub-proton scales, and are relevant for modelling the energy dissipation in a broad range of astrophysical systems.
An important issue for shaped minimized sonic booms is whether turbulence-induced distortions will adversely affect the benefits gained by shaping. This question was considerably simplified by two recent results. The first is the finding that the loudness of sonic booms is well quantified by loudness. The second is that loudness of a shaped boom is dominated by the shock waves. The issue is now the effect of turbulence on weak (1 psf or less) sonic booms. Since it is clear that molecular relaxation effects have a significant effect on shock structure and loudness, turbulence effects must be examined in conjunction with relaxation-thickened shocks. This analysis must be directed toward loudness calculations and include all pertinent mechanisms.
The interaction of the drift wave (DW) turbulence and zonal flow (ZF) is investigated with the modified Hasegawa–Mima equation taking into account the backreaction of ZF velocity on DW turbulence. It is shown that the y-averaged enstrophy of DW turbulence and the velocity of ZF are intrinsically related. By utilizing this feature, a nonlinear stage of DW modulational instability is considered within the framework of the wave kinetic equation. It is shown that in this approximation, the nonlinear stage of the modulational instability results in the collapsing solutions, accompanied by the “wave breaking” phenomenon. Numerical simulations based on the Hasegawa–Mima equation show that for a weak DW turbulence, Φ ̃ = ( e φ ̃ / T e ) ( L n / ρ s ) ⪝ 1 , the collapsing-like features on both ZF and y-averaged enstrophy of DW turbulence decay in time and then re-emerge again at different locations. For the case of a strong DW turbulence, Φ ̃ > 1 , where nonlinear interactions of DW harmonics dominate, stable spatial structures of ZF and y-averaged enstrophy of DW turbulence emerge.
Abstract Predictions of heat load widths λ q based on particle orbits alone are very pessimistic. This paper shows that pedestal peeling-ballooning (P-B) magnetohydrodynamic (MHD) turbulence broadens the stable scrape-off layer (SOL) by the transport, or spreading, of fluctuation energy from the pedestal. λ q is seen to increase with Γ ε , the fluctuation energy density flux. We elucidate the fundamental physics of the spreading process. Γ ε increases with pressure fluctuation correlation length. P-B turbulence is seen to be especially effective at spreading, on account of its large effective mixing length. Spreading is shown to be a multiscale process, which is enhanced by the synergy of large and small-scale modes. Pressure fluctuation skewness correlates well with the spreading flux–with the zero crossing of skewness and Γ ε spatially coincident–suggesting the role of coherent fluctuation structures and the presence of intermittency in λ q broadening. λ q ∼ B p − 1 scaling persists for the broadened SOL. We show that the spreading flux increases for increasing pedestal pressure gradient ∇ P 0 and for decreasing pedestal collisionality υ ped ∗ . This trend is due to the dominance of peeling modes for large ∇ P 0 and low υ ped ∗ . Ultimately, we see that a state of weak MHD turbulence, as for small ELMs, is very attractive for heat load management. Our findings have transformative implications for future fusion reactor designs and call for experimental investigations to validate the observed trends.
Weak scattering, wave optics results are given for the first-order log-amplitude and phase, and for the second-order phase and angle of arrival, of radiation propagating through a turbulent planetary atmosphere. The analytical results include the effects of coupling to the inhomogeneous background upon which the turbulence is superimposed. Such coupling significantly affects signal characteristics in stellar and spacecraft occultation experiments, in contrast to situations involving much shorter paths, as in terrestrial propagation experiments.
We propose a conservative Galerkin scheme for the quasilinear model in three-dimensional momentum space and three-dimensional spectral space, with cylindrical symmetry. We construct an unconditionally conservative weak form and use a discretization that preserves conservation properties independent of the wave emission probability. The discrete operators, combined with a consistent quadrature rule, preserve all the conservation laws rigorously. The proposed scheme is quite general: it works for both relativistic and non-relativistic systems, for both magnetized and unmagnetized plasmas, and even for problems with time-dependent dispersion relations. We represent the particle distribution by continuous basis functions and use discontinuous basis functions for the wave spectral energy density, which enables the application of a positivity-preserving technique. We adopt the marching simplex algorithm, designed initially for computer graphics, for numerical integration on the resonance manifold. Furthermore, the numerical examples with a bump-on-tail initial configuration show how the unstable waves produce strong momentum space diffusion.
We describe the interaction of parallel-propagating Alfvén waves with ion-acoustic waves and other Alfvén waves, in magnetized, high- β collisionless plasmas. This is accomplished through a combination of analytical theory and numerical fluid simulations of the Chew–Goldberger–Low (CGL) magnetohydrodynamic (MHD) equations closed by Landau-fluid heat fluxes. An asymptotic ordering is employed to simplify the CGL-MHD equations and derive solutions for the deformation of an Alfvén wave that results from its interaction with the pressure anisotropy generated either by an ion-acoustic wave or another, larger-amplitude Alfvén wave. The difference in time scales of acoustic and Alfvénic fluctuations at high- β means that interactions that are local in wavenumber space yield little modification to either mode within the time it takes the acoustic wave to Landau damp away. Instead, order-unity changes in the amplitude of Alfvénic fluctuations can result after interacting with frequency-matched acoustic waves. Additionally, we show that the propagation speed of an Alfvén-wave packet in an otherwise homogeneous background is a function of its self-generated pressure anisotropy. This allows for the eventual interaction of separate co-propagating Alfvén-wave packets of differing amplitudes. The results of the CGL-MHD simulations agree well with these predictions, suggesting that theoretical models relying on the interaction of these modes should be reconsidered in certain astrophysical environments. Applications of these results to weak Alfvénic turbulence and to the interaction between the compressive and Alfvénic cascades in strong, collisionless turbulence are also discussed.
Biomass-burning smoke drives large uncertainty in climate projections of the Earth's radiative balance. This is due to the chemical and physical evolution of smoke and its impact on clouds and radiation. Here we focus on the southeastern Atlantic region and its inflow of African biomass-burning smoke during August 2017. We evaluate smoke properties and processes in two coupled earth-system models, the Energy Exascale Earth System Model (E3SM) and Community Earth System Model (CESM). These are compared against in situ aircraft observations from two field campaigns, ObseRvations of Aerosols above CLouds and their intEractionS (ORACLES) and CLoud–Aerosol–Radiation Interaction and Forcing: Year 2017 (CLARIFY-2017). Observations reveal an increase and subsequent decrease in smoke mean diameter, and a steady decrease in the mass ratio of organic aerosol (OA) to black carbon aerosol (BC) (OA:BC) over 4–12 days of aging, neither captured by the base models. Implementation of a photolytic loss scheme for secondary organic aerosol (SOA)—as a proxy for other heterogeneous volatilization chemistry—and a ∼1-day conversion for primary OA to SOA significantly improves the representation of this loss. In the boundary layer, both models show dimethyl sulfide driving a large increase in the sulfate aerosol mass fraction from the free troposphere, which is consistent with observations. Finally, models tend to underpredict cloud droplet number concentration partially due to weak modeled turbulent updraft strength, and model performance improves when the parameterized turbulent updraft strength is increased substantially. These results are expected to provide insights into future model development to reduce climate model uncertainties.
Abstract Using the assumption of a weak normalized turbulent viscosity, usually valid in practice, the modulated-transport model (Stoltzfus-Dueck 2012 Phys. Plasmas 19 055908) is generalized to allow the turbulent transport coefficient to vary in an arbitrary way on radial and poloidal position. The new approach clarifies the physical interpretation of the earlier results and significantly simplifies the calculation, via a boundary-layer asymptotic method. Rigorous detailed appendices verify the result of the simple boundary-layer calculation, also demonstrating that it achieves the claimed order of accuracy and providing a concrete prediction for the strong plasma flows in the immediate vicinity of the last closed flux surface. The new formulas are used to predict plasma rotation at the core-edge boundary, in cases with and without externally applied torque. Dimensional formulas and extensive discussion are provided, to support experimental application of the new model.
ABSTRACT Our aim in this work is to identify and explain the necessary conditions required for an energetic explosion of a Chandrasekhar-mass white dwarf. We construct and analyse weakly compressible turbulence models with nuclear burning effects for carbon/oxygen plasma at a density expected for the deflagration-to-detonation transition (DDT) to occur. We observe the formation of carbon deflagrations and transient carbon detonations at early times. As turbulence becomes increasingly inhomogeneous, sustained carbon detonations are initiated by the Zel’dovich reactivity gradient mechanism. The fuel is suitably preconditioned by the action of compressive turbulent modes with wavelength comparable to the size of resolved turbulent eddies; no acoustic wave is involved in this process. Oxygen detonations are initiated, aided either by reactivity gradients or by collisions of carbon detonations. The observed evolutionary time-scales are found to be sufficiently short for the above process to occur in the expanding, centrally ignited massive white dwarf. The inhomogeneous conditions produced prior to the DDT might be of consequence for the chemical composition of the outer ejecta regions of Type Ia supernovae from the single degenerate channel, and offer the potential for validation of the proposed model.
A kinetic equation for interplanetary cosmic rays is set up with the aid of weak-plasma-turbulence theory for an idealized radially symmetric model of the interplanetary magnetic field. As a starting point, this treatment invokes the Vlasov equation instead of the traditional Fokker-Planck equation. Quasi-linear theory is applied to obtain a momentum diffusion equation for the heliocentric frame of reference which describes the interaction of cosmic rays with convecting magnetic irregularities in the solar-wind plasma. Under restricted conditions, the well-known equation of solar modulation can be obtained from this kinetic equation.
The various mechanisms which generate broadband noise on a range of rotors studied include load fluctuations due to inflow turbulence, due to turbulent boundary layers passing the blades' trailing edges, and due to tip vortex formation. Existing analyses are used and extensions to them are developed to make more accurate predictions of rotor noise spectra and to determine which mechanisms are important in which circumstances. Calculations based on the various prediction methods in existing experiments were compared. The present analyses are adequate to predict the spectra from a wide variety of experiments on fans, full scale and model scale helicopter rotors, wind turbines, and propellers to within about 5 to 10 dB. Better knowledge of the inflow turbulence improves the accuracy of the predictions. Results indicate that inflow turbulence noise depends strongly on ambient conditions and dominates at low frequencies. Trailing edge noise and tip vortex noise are important at higher frequencies if inflow turbulence is weak. Boundary layer trailing edge noise, important, for large sized rotors, increases slowly with angle of attack but not as rapidly as tip vortex noise.
A three-dimensional laser Doppler anemometer (3D-LDA) was used in a wind tunnel to measure a jet in a crossflow. Measurements were made in the vicinity of a 5-cm-diam jet which issued normally into a 10.65 m/sec wind tunnel crossflow; the velocity ratio Vjet/Vinf was 8. Detailed lateral surveys were made at two elevations (z = cm and 2 cm); both elevations were within the region affected by the boundary layer on the plate. The results are believed to provide reliable velocity field information in the boundary layer of the jet in a crossflow. Turbulence information also is available and believed to be roughly correct, although it may be subject to broadening effects for the lower values of turbulence. A weak vortex pair was observed in the wake at the plate surface. This structure existed in the boundary layer and built confidence because the 3D-LDA was, indeed, able to resolve fine detail in the wake. The capabilities of the 3D-LDA not only allow the making of the velocity surveys, but can be utilized to follow mean streamlines in the flow.
The various source mechanisms which generate broadband noise on a range of rotors are reviewed. Analyses of these mechanisms are presented and compared to existing experimental data. The sources considered are load fluctuations due to inflow turbulence, due to turbulent blade boundary layers passing the trailing edge, and due to tip vortex formation turbulence. Vortex shedding noise due to laminar boundary layers and blunt trailing edges is not considered in detail as it can be avoided in most cases. Present analyses are adequate to predict the spectra from a wide variety of experiments on fans, helicopter rotors, and wind turbines to within about 5 to 10 dB. Better knowledge of the inflow turbulence improves the accuracy of the predictions. Inflow turbulence noise depends strongly on ambient conditions and dominates at low frequencies. Trailing edge and tip vortex noise are important at higher frequencies if inflow turbulence is weak. Boundary layer trailing edge noise increases slowly with angle of attack but not as rapidly as tip vortex formation noise. Tip noise can be important at high angles of attack for wide chord, square edge tips.