Preliminary Turbulence Results - CMSO

1 downloads 0 Views 5MB Size Report
Mar 11, 2014 - 1) Provide overview of plasma generation in SSX. 2) Demonstrate ability to conduct turbulence research in a laboratory plasma (two ways):.
Turbulence Analysis of an MHD wind-tunnel David Schaffner

Swarthmore College, NSF Center for Magnetic Self-Organization

with contributions from M. Brown, V. Lukin, A. Wan ‘15

CMSO General Meeting, Santa Fe, NM March 11 2014 Research supported by US DOE and NSF

Motivating Questions Is the statistical character of MHD turbulence universal?

How can laboratory plasmas be used to study MHD turbulence?

Example: Can laboratory turbulence be compared to solar wind? Space: Solar Wind*

Laboratory: SSX

fλp = 2MHz

*Cluster FGM and STAFF-SC Data: Sahraoui, PRL 2009

Objectives for this Talk

1) Provide overview of plasma generation in SSX 2) Demonstrate ability to conduct turbulence research in a laboratory plasma (two ways): A. Turbulence spectra, variance anisotropy B. Intermittency and scaling with helicity

3) Give overview of simulation efforts and preliminary comparisons

Image: Turbulent Channel Structures, Melissa Green Department of Mechanical and Aerospace Engineering, Princeton

OBJECTIVE 1: OVERVIEW OF PLASMA IN SSX

The SSX Laboratory 10kV/100kA Pulsed power 1mF banks

Cylindrical vacuum chamber (D = 0.5 m, L = 1 m)

Schematic + Simulation

High voltage plasma guns on each end

SSX Diagnostics measure B-field, density, flow, soft X-ray light, ion temperature, and gun I/V Copper Flux Conserver (Wind Tunnel) 16chan, 3 axis Bdot Probe 2-Side (Z-axis) Mach Probe

Off Midplane: HeNe Interferometer Gun Current/Voltage

Ion Doppler Spectrometer (IDS) Al/Zr Metal Soft X-ray Photodiode

Plasma gun produces spheromak; narrow cylinder boundary cause tilt and twist Initial State

Final State Turbulence

Taylor relaxation or selective decay of a plasma is a process where magnetic energy,

is minimized under the constraint that magnetic helicity,

is conserved.

Taylor State (Tilt and Twist) Double Helix Shape

Instability of field structure provides energy for cascade

Hall MHD simulation clearly illustrates plasma evolution

Initial State - Spheromak

Hall MHD simulation clearly illustrates plasma evolution

Field structure tilts

Hall MHD simulation clearly illustrates plasma evolution

Field structure twists

Hall MHD simulation clearly illustrates plasma evolution

Fields turbulent (Time frame for analysis)

Hall MHD simulation clearly illustrates plasma evolution

Fluctuations begin to decay

Hall MHD simulation clearly illustrates plasma evolution

Final State: Taylor State (double helix)

Typical Plasma Parameters Temp: = 25eV, = 10eV β = 0.1 Field: = 5kG  ρi = 0.1cm (tunnel diam = 15cm) Density: = 1x1015cm-3  λi = δi = c/ωpi = 0.54cm Velocity: VA > Cs > Axial Flow 256km/s > 31km/s > 20km/s Collisional Plasma – MFP ~ 0.2-3cm

OBJECTIVE 2: TURBULENCE RESEARCH ON SSX – SPECTRA AND ANISOTROPY

B-field spectra extracted from dB/dt signal

B-field spectra extracted from measured dB/dt time series using wavelet transform

Power spectrum

B-field spectra extracted from dB/dt signal

Ensemble Average:

x40+ shots

frequency

B-field spectra extracted from measured dB/dt time series using wavelet transform

time

Power spectrum

Spectra indicate power-law behavior Power-law scaling in Bfield, density and flow as function of frequency

Magnetic fluctuations steepen before flow fluctuations—suggest energy initially in fields, not flow fλp = 2-3MHz

Spectra indicate power-law behavior

If average flow is used to invoke Taylor Hypothesis: Break in B-field spectra near Doppler-shifted ion inertial length fλp = 2-3MHz

B-field spectra exhibit variance anisotropy Power of each component divided into perp/para using local B-field reference (at each timestep)

B(t) B(f)para

B(f)perp

Perp (red) Para (blue)

B-field spectra exhibit variance anisotropy Power of each component divided into perp/para using local B-field reference (at each timestep)

More perp power than parallel Ratio increases with frequency (up to ~1MHz)

Perp (red) Para (blue)

Ratio has power-law scaling ≈ 1/3 Low frequency ratio scales as frequency to the 1/3

Ratio peaks around 1MHz

Ratio shows isotropy at higher frequencies

Ratio shows some time dependence Black - Initial fluctuations (near isotropy)

Scaling increases from purple to blue to light blue

Spatial wavenumber spectra measured directly

Lower resolution spectra, but slopes are closer to a Kolmogorov -5/3 16 measurements along radius Scales between 8cm and 1cm

Frequency spectra steeper than wavenumber spectra Aligned by invoking Taylor Hypothesis with bulk flow of 20km/s Wavenumber spectra does not likely probe dissipation range Some effect causes frequency spectra to be steeper? Some type of dissipation? Not really a direct comparison: radial k versus axial freq

OBJECTIVE 2: TURBULENCE RESEARCH ON SSX – INTERMITTENCY/HELICITY SCALING

PDF constructed from increments in dB/dt

Probability Distribution Functions of Increments constructed by taking difference of dB/dt values at different time separations

PDF of increments have fat tails (intermittency): shows B-dot fluctuations are highly structured

Flatness increases with decreasing time scale

Flatness = normalized 4th moment of a PDF

Flatness = quantification of fat tails, departure from Gaussian  the larger the flatness, the more intermittency is observed

Power-law like scaling

Magnetic Helicity can be linearly modified

Mag Flux, Φ

Vary flux by varying strength of field produced by coil

Flatness (intermittency) scales with injected helicity Flatness increases at all timescales with helicity

Average flatness increases with helicity beyond 100 μWb2

Ion temperature bursts scale with helicity, but integrated SXR signal scales inversely Intermittency increases with helicity

Soft X-ray signal decreases with helicity Ti intermittency increases with helicity—mean value hold steady Possibly, trends suggest intermittency caused by current sheets—seen in solar wind turbulence simulation (Greco ApJ 09)

That is, helicity increases number of sites, but decreases their size

Shape/slopes of spectra generally unchanged with increasing helicity Spectra staggered for clarity

Fits (both high and low frequency) and breakpoints do not scale with changing helicity

COMPARISON TO HIFI SIMULATION

Similar turbulence statistics observed

Comparable slope of spectra observed

Similar wavenumber spectra observed

Similar turbulence statistics observed

Variance anisotropy observed

Intermittency increases with smaller scales

Concluding Remarks:

How well does lab plasma turbulence compare to solar wind turbulence? • B-field frequency spectra steeper than solar wind—flow fluctuations driven by magnetic? • Wavenumber spectra closer to Kolmogorov • Variance anisotropy is observed • Spectral break consistent with ion scale if Taylor Hypothesis invoked

• Bdot fluctuations exhibits intermittency—flatness decreases w/scale • Evidence for connection between current sheets/intermittency— found through lab helicity scaling

What is next?—Questions to Answer… • Is the Taylor Hypothesis Valid? – We plan to use multi-tipped probes aligned with the plasma flow to get a direct wavenumber spectrum to compare to the frequency spectrum w/Taylor Hypothesis

• What kind of modes are generated?—we plan to use a probe to measure the simultaneous Bfield and density fluctuations—some theories predict such correlation coupled with anisotropy can determine mode type (i.e. fast, slow, Alfven) [Klein 2012]

What is next?—Questions to Answer… • Does (system) size matter? – Comparison of experiment and simulation are favorable. We want to test how(if) spectra changes if volume is scaled up in simulation, in anticipation of possible larger experiment.

Thank you for your attention.