Gyrokinetic Turbulence arXiv:0704.0044; 0806.1069 Alexander - - PowerPoint PPT Presentation

gyrokinetic turbulence
SMART_READER_LITE
LIVE PREVIEW

Gyrokinetic Turbulence arXiv:0704.0044; 0806.1069 Alexander - - PowerPoint PPT Presentation

Cracow 9.10.08 Gyrokinetic Turbulence arXiv:0704.0044; 0806.1069 Alexander Schekochihin Imperial College London Oxford w i t h Steve Cowley (UCLA Culham) Greg Howes (Berkeley Iowa) Bill Dorland, Tomo Tatsuno (Maryland) Michael Barnes


slide-1
SLIDE 1

Gyrokinetic Turbulence

Cracow 9.10.08 Alexander Schekochihin Imperial College London → Oxford w i t h Steve Cowley (UCLA → Culham) Greg Howes (Berkeley → Iowa) Bill Dorland, Tomo Tatsuno (Maryland) Michael Barnes (Maryland → Oxford) Eliot Quataert (Berkeley) Greg Hammett (Princeton) Gabriel Plunk (UCLA → Maryland) Reprints/references on http://www2.imperial.ac.uk/~aschekoc/ arXiv:0704.0044; 0806.1069

slide-2
SLIDE 2

Turbulence is Multiscale Disorder

[Image: Y. Kaneda et al., Earth Simulator, isovorticity surfaces, 40963]

slide-3
SLIDE 3

Turbulence is Multiscale Disorder

[Image: Y. Kaneda et al., Earth Simulator, isovorticity surfaces, 40963]

slide-4
SLIDE 4

Turbulence is Multiscale Disorder

[Image: Y. Kaneda et al., Earth Simulator, isovorticity surfaces, 40963]

slide-5
SLIDE 5

Turbulence is Multiscale Disorder

[Image: Y. Kaneda et al., Earth Simulator, isovorticity surfaces, 40963]

slide-6
SLIDE 6

Turbulence: A Nonlinear Route to Dissipation

energy injected E(k) k Inertial range energy dissipated energy transported ε

slide-7
SLIDE 7

Turbulence: A Nonlinear Route to Dissipation

energy injected E(k) k Inertial range energy dissipated energy transported If cascade is local, intermediate scales fill up Big whorls have little whorls That feed on their velocity, And little whorls have lesser whorls And so on to viscosity.

  • L. F. Richardson 1922

ε

slide-8
SLIDE 8

Turbulence: A Nonlinear Route to Dissipation

energy injected E(k) k Inertial range energy dissipated energy transported If cascade is local, intermediate scales fill up Big whorls have little whorls That feed on their velocity, And little whorls have lesser whorls And so on to viscosity.

  • L. F. Richardson 1922

k–5/3 ε

K41

slide-9
SLIDE 9

Plasma Turbulence: Analogous?

Turbulence in the solar wind

[Bale et al. 2005, PRL 94, 215002]

k–5/3 k–7/3 k–1/3

slide-10
SLIDE 10

Plasma Turbulence Extends to Collisionless Scales

Turbulence in the solar wind

[Bale et al. 2005, PRL 94, 215002]

k–5/3 k–7/3 k–1/3 λmfp ~ 108 km (~1 AU) ρi ~ 102 km

slide-11
SLIDE 11

Plasma Turbulence Extends to Collisionless Scales

Interstellar medium: “Great Power Law in the Sky” L ~ 1013 km (~100 pc) λmfp ~ 107 km ρi ~ 104 km k–5/3

[Armstrong et al. 1995, ApJ 443, 209]

slide-12
SLIDE 12

Plasma Turbulence Extends to Collisionless Scales

Intracluster (intergalactic) medium L ~ 1019 km (~1 Mpc) λmfp ~ 1016 km (~1 kpc) ρi ~ 104 km

Hydra A cluster [Vogt & Enßlin 2005, A&A 434, 67]

slide-13
SLIDE 13

Plasma Turbulence Is Kinetic

Turbulence in the solar wind

[Bale et al. 2005, PRL 94, 215002]

k–5/3 Alfvén waves k–7/3 k–1/3 KAW

CASCADE CAS- CADE “DIS- SI- PA- TION”

  • What is cascading

in kinetic turbulence? (What is conserved?) What do the observed spectra tell us and how do we explain them?

  • Dissipation

(as usually understood) is “collisionless” (Landau damping) How does that heat particles? (ions, electrons, minority ions)

“Inertial range” “Dissipation range”

slide-14
SLIDE 14

Plasma Turbulence Ab Initio

slide-15
SLIDE 15

Plasma Turbulence Ab Initio

Work done

slide-16
SLIDE 16

Plasma Turbulence Ab Initio

Work done Entropy produced:

Boltzmann 1872

slide-17
SLIDE 17

Plasma Turbulence Ab Initio

Work done Entropy produced:

slide-18
SLIDE 18

Plasma Turbulence Ab Initio

Work done Entropy produced:

slide-19
SLIDE 19

Plasma Turbulence Ab Initio

Work done Heating: Fluctuation energy budget:

–TδS energy heating

slide-20
SLIDE 20

Plasma Turbulence: Generalised Energy Cascade

–TδS energy heating

Fowler 1968 Krommes & Hu 1994 Krommes 1999 Sugama et al. 1996 Hallatschek 2004 Howes et al. 2006 Schekochihin et al. 2007 Scott 2007

Generalised energy = free energy of the particles + fields

arXiv:0806.1069

slide-21
SLIDE 21

Plasma Turbulence: Generalised Energy Cascade

–TδS energy heating

Fowler 1968 Krommes & Hu 1994 Krommes 1999 Sugama et al. 1996 Hallatchek 2004 Howes et al. 2006 Schekochihin et al. 2007 Scott 2007

Generalised energy = free energy of the particles + fields

arXiv:0806.1069

Landau damping is a redistribution between e-m fluctuation energy and (negative) perturbed entropy (free energy). It was pointed out already by Landau 1946 that δfs does not decay: “ballistic response”

slide-22
SLIDE 22

Plasma Turbulence: Analogous to Fluid, But…

small scales in 3D physical space small scales in 6D phase space

energy heating arXiv:0806.1069 –TδS

slide-23
SLIDE 23

Plasma Turbulence: Analogous to Fluid, But…

energy heating

small scales in 6D phase space In gyrokinetic turbulence, the velocity-space and x-space cascades are intertwined, giving rise to a single phase-space cascade

arXiv:0806.1069 –TδS

slide-24
SLIDE 24

Plasma Turbulence: Analogous to Fluid, But…

energy heating arXiv:0806.1069 –TδS

SO, IDEA #1: GENERALISED ENERGY CASCADE THROUGH PHASE SPACE

slide-25
SLIDE 25

Critical Balance

IDEA #2: CRITICAL BALANCE

slide-26
SLIDE 26

Critical Balance

  • Strong nonlinearity:

Critical balance as a physical principle proposed for Alfvénic turbulence by Goldreich & Sridhar 1995 [ApJ 438, 763] More generally, one might argue that in a magnetised plasma, parallel linear propagation scale and perpendicular nonlinear interaction scale will adjust to each other and the turbulent cascade route will be determined by this principle

  • Strong anisotropy:

In magnetised plasma, confirmed by numerics (MHD) and observations (solar wind, ISM)

  • Weak turbulence drives itself into strong regime
  • 2D turbulence (“overstrong”) parallel-decorrelates

and returns to critical balance

slide-27
SLIDE 27

What Is Gyrokinetics?

  • Strong nonlinearity:

(critical balance as an ordering assumption)

  • Strong anisotropy: (this is the small parameter!)

[Howes et al. 2006, ApJ 651, 590]

slide-28
SLIDE 28

What Is Gyrokinetics?

  • Strong nonlinearity:

(critical balance as an ordering assumption)

  • Strong anisotropy: (this is the small parameter!)
  • Finite Larmor radius:

Low frequency GK ORDERING:

  • Weak collisions:

[Taylor & Hastie 1968, Plasma Phys. 10, 479; Rutherford & Frieman 1968, Phys. Fluids 11, 569; Catto 1977, Plasma Phys. 20, 719; Frieman & Chen 1982, Phys. Fluids 443, 209; for our derivation, notation, etc. see Howes et al. 2006, ApJ 651, 590]

slide-29
SLIDE 29

Particle dynamics can be averaged over the Larmor orbits and everything reduces to kinetics of Larmor rings centered at and interacting with the electromagnetic fluctuations.

Gyrokinetics: Kinetics of Larmor Rings

[Howes et al. 2006, ApJ 651, 590]

Catto transformation

  • nly two velocity variables,

i.e., 6D → 5D

slide-30
SLIDE 30

Particle dynamics can be averaged over the Larmor orbits and everything reduces to kinetics of Larmor rings centered at and interacting with the electromagnetic fluctuations.

Gyrokinetics: Kinetics of Larmor Rings

[Howes et al. 2006, ApJ 651, 590]

Catto transformation

+ Maxwell’s equations (quasineutrality and Ampère’s law)

slide-31
SLIDE 31

Gyrokinetics: Kinetics of Larmor Rings

[Howes et al. 2006, ApJ 651, 590]

+ Maxwell’s equations (quasineutrality and Ampère’s law)

Averaged gyrocentre drifts:

  • E×B0 drift
  • ∇B drift
  • motion along

perturbed fieldline Averaged wave-ring interaction

slide-32
SLIDE 32

Why is Gyrokinetics Valid?

KAW: k|| ~ k⊥

1/3

k–5/3 k–7/3 energy injected collisional (fluid) collisionless (kinetic) Alfvén waves: k|| ~ k⊥

2/3

GYROKINETICS

FLUID THEORY Because anisotropy makes frequencies low.

Cyclotron frequency

  • nly reached deep

in the dissipation range

arXiv:0704.0044

electron Landau damping ion Landau damping

slide-33
SLIDE 33

Why is Gyrokinetics Useful?

  • Because it is a simplifying

analytical step that is a natural staring point for further theory [Howes et al. 2006, ApJ 651, 690

Schekochihin et al., arXiv:0704.0044]

  • Because it reduces the

kinetic problem to 5D, making it numerically tractable (publicly available codes

developed in fusion research: e.g., GS2, GENE, GYRO…)

arXiv:0704.0044

slide-34
SLIDE 34

Why is Gyrokinetics Useful?

Alfvén-wave turbulence in the SW

[by Bale et al. 2005, PRL 94, 215002]

k–5/3 Alfvén waves k–7/3 k–1/3 KAW

  • Because it is a simplifying

analytical step that is a natural staring point for further theory [Howes et al. 2006, ApJ 651, 690

Schekochihin et al., arXiv:0704.0044]

  • Because it reduces the

kinetic problem to 5D, making it numerically tractable (publicly available codes

developed in fusion research: e.g., GS2, GENE, GYRO…)

slide-35
SLIDE 35
  • Because it is a simplifying

analytical step that is a natural staring point for further theory [Howes et al. 2006, ApJ 651, 690

Schekochihin et al., arXiv:0704.0044]

  • Because it reduces the

kinetic problem to 5D, making it numerically tractable (publicly available codes

created in fusion research: e.g., GS2, GENE, GYRO…)

Why is Gyrokinetics Useful?

Alfvén-wave turbulence using GS2 (by Greg Howes)

[Howes et al. 2008, PRL 100, 065004]

slide-36
SLIDE 36

Why is Gyrokinetics Useful?

Alfvén-wave turbulence in the SW

[by Bale et al. 2005, PRL 94, 215002]

k–5/3 Alfvén waves k–7/3 k–1/3 KAW

  • Because it is a simplifying

analytical step that is a natural staring point for further theory [Howes et al. 2006, ApJ 651, 690

Schekochihin et al., arXiv:0704.0044]

  • Because it reduces the

kinetic problem to 5D, making it numerically tractable (publicly available codes

created in fusion research: e.g., GS2, GENE, GYRO…)

slide-37
SLIDE 37

Kinetics vs. Fluid Models: What Is New?

k–5/3 Alfvén waves k–7/3 k–1/3 KAW

CASCADE CAS- CADE “DIS- SI- PA- TION”

  • What is cascading

in kinetic turbulence? (What is conserved?) What do the observed spectra tell us and how do we explain them?

  • Dissipation

(as usually understood) is “collisionless” (Landau damping) How does that heat particles? (ions, electrons, minority ions)

“Inertial range” “Dissipation range”

Alfvén-wave turbulence in the SW

[by Bale et al. 2005, PRL 94, 215002]

slide-38
SLIDE 38

Gyrokinetics: Kinetics of Larmor Rings

+ Maxwell’s equations (quasineutrality and Ampère’s law)

[Howes et al. 2006, ApJ 651, 590]

SO, IDEA #3: GYROAVEARGED KINETIC THEORY AT LOW FREQUENCIES

  • Only two velocity variables, i.e., 6D → 5D
  • All high-frequency stuff averaged out
slide-39
SLIDE 39

Generalised Energy in Gyrokinetics

energy heating –TδS

+ Maxwell’s equations (quasineutrality and Ampère’s law)

[Howes et al. 2006, ApJ 651, 590]

slide-40
SLIDE 40

Generalised Energy in Gyrokinetics

arXiv:0704.0044 energy –TδS

+ Maxwell’s equations (quasineutrality and Ampère’s law)

[Howes et al. 2006, ApJ 651, 590]

slide-41
SLIDE 41

The Grand Kinetic Cascade

k–5/3 k–7/3 energy injected Alfvén waves KAW electron Landau damping ion Landau damping

arXiv:0704.0044

MHD KINETICS

slide-42
SLIDE 42

The Grand Kinetic Cascade

k–5/3 k–7/3

slow waves entropy fluctuations Alfvén waves

KAW Alfvén waves slow waves entropy mode

arXiv:0704.0044

ion Landau damping electron Landau damping

MHD KINETICS

energy injected

slide-43
SLIDE 43

The Grand Kinetic Cascade

k–5/3 k–7/3

compressive fluctuations Alfvén waves

KAW Alfvén waves compressive fluctuations

arXiv:0704.0044

ion Landau damping electron Landau damping energy injected

slide-44
SLIDE 44

The Grand Kinetic Cascade

k–5/3 k–7/3 KAW all modes mixed

arXiv:0704.0044

ion Landau damping electron Landau damping energy injected

slide-45
SLIDE 45

The Grand Kinetic Cascade

k–5/3 k–7/3

kinetic Alfvén waves entropy cascade

arXiv:0704.0044

ion Landau damping electron Landau damping energy injected

slide-46
SLIDE 46

The Grand Kinetic Cascade

k–5/3 k–7/3

arXiv:0704.0044

ion Landau damping electron Landau damping Dissipated by collisions ion heating electron heating So the cascade split at ion gyroscale determines relative heating of the species energy injected

slide-47
SLIDE 47

The Grand Kinetic Cascade

k–5/3 k–7/3 energy injected ion Landau damping electron Landau damping Dissipated by collisions ion heating electron heating So the cascade split at ion gyroscale determines relative heating of the species

SO, IDEA #4: DECOUPLING-RECOUPLING OF SUBCASCADES → HEATING

slide-48
SLIDE 48

Ion Gyroscale Transition: GK DNS by G. Howes

Alfvén-wave turbulence in the solar wind [by Bale et al. 2005, PRL 94, 215002]

k–5/3 Alfvén waves k–7/3 k–1/3 KAW

Alfvén-wave turbulence using GS2 [by Howes et al. 2008, PRL 100, 065004]

slide-49
SLIDE 49

Ion Gyroscale Transition: GK DNS by G. Howes

Alfvén-wave turbulence in the solar wind [by Bale et al. 2005, PRL 94, 215002]

k–5/3 Alfvén waves k–7/3 k–1/3 KAW

Alfvén-wave turbulence using GS2 [by Howes et al. 2008, PRL 100, 065004]

slide-50
SLIDE 50

Main Points So Far

  • IDEA #1: Kinetic turbulence is a generalised energy cascade

in phase space towards collisional scales

  • IDEA #2: Cascade is anisotropic and critically balanced

(linear parallel propagation scale = nonlinear perpendicular interaction scale)

  • IDEA #3: Can be described by gyrokinetics — gyroangle

averaged low frequency kinetics of Larmor rings

  • IDEA #4: Cascade splits into various non-energy-exchanging

channels in different ways, depending on scales (some of these described by fluid/hybrid models); mixing and resplitting of these subcascades at ion gyroscale determines relative heating of the two species Details are in these preprints: arXiv:0704.0044, 0806.1069

slide-51
SLIDE 51

Further Topics

  • Alfvénic turbulence and passive compressive fluctuations

in the inertial range

  • Energetic minority ions and their heating
  • Kinetic Alfvén wave turbulence in the “dissipation range”
  • Entropy cascade in phase space and nonlinear phase mixing
  • Pressure anisotropies and resulting instabilities
  • Magnetogenesis
  • The answer to the general question about life, universe,

and everything… Details are in these preprints: arXiv:0704.0044, 0806.1069

slide-52
SLIDE 52

Further Topics

  • Alfvénic turbulence and passive compressive fluctuations

in the inertial range

  • Energetic minority ions and their heating
  • Kinetic Alfvén wave turbulence in the “dissipation range”
  • Entropy cascade in phase space and nonlinear phase mixing
  • Pressure anisotropies and resulting instabilities
  • Magnetogenesis
  • The answer to the general question about life, universe,

and everything… Details are in these preprints: arXiv:0704.0044, 0806.1069

slide-53
SLIDE 53

Kinetic Reduced MHD

k–5/3 k–7/3 collisional (fluid) collisionless (kinetic)

GYROKINETICS

FLUID THEORY

magnetised ions isothermal electrons

arXiv:0704.0044

ion Landau damping electron Landau damping energy injected

slide-54
SLIDE 54

KRMHD: Alfvén Waves

  • Alfvénic fluctuations and

rigourously satisfy Reduced MHD Equations:

[Strauss 1976, Phys. Fluids 19, 134] [Schekochihin et al., arXiv:0704.0044

  • cf. Higdon 1984, ApJ 285, 109; Lithwick & Goldreich 2001, ApJ 562, 279]
slide-55
SLIDE 55

KRMHD: Alfvén Waves

  • Alfvénic fluctuations and

rigourously satisfy Reduced MHD Equations:

[Kadomtsev & Pogutse 1974,

  • Sov. Phys. JETP 38, 283]

[Schekochihin et al., arXiv:0704.0044

  • cf. Higdon 1984, ApJ 285, 109; Lithwick & Goldreich 2001, ApJ 562, 279]
slide-56
SLIDE 56

KRMHD: Alfvén Waves

  • Alfvénic fluctuations and

rigourously satisfy Reduced MHD Equations:

[Kadomtsev & Pogutse 1974,

  • Sov. Phys. JETP 38, 283

Strauss 1976, Phys. Fluids 19, 134]

  • Alfvén-wave cascade is indifferent to collisions and damped
  • nly at the ion gyroscale
  • The GS95 theory describes this part of the turbulence
  • Alfvén waves are decoupled from density and magnetic-field-strength

fluctuations (slow waves and entropy mode in the fluid limit)

[Schekochihin et al., arXiv:0704.0044

  • cf. Higdon 1984, ApJ 285, 109; Lithwick & Goldreich 2001, ApJ 562, 279]
slide-57
SLIDE 57

KRMHD: Alfvén Waves

  • Alfvénic fluctuations and

rigourously satisfy Reduced MHD Equations:

[Kadomtsev & Pogutse 1974,

  • Sov. Phys. JETP 38, 283

Strauss 1976, Phys. Fluids 19, 134] [Schekochihin et al., arXiv:0704.0044

  • cf. Higdon 1984, ApJ 285, 109; Lithwick & Goldreich 2001, ApJ 562, 279]

SO, IDEA #5: DECOUPLED RMHD ALFVENIC CASCADE IN THE INERTIAL RANGE

slide-58
SLIDE 58

ISM: Density Fluctuations

Electron-density fluctuations in the interstellar medium [Armstrong et al. 1995, ApJ 443, 209]

k–5/3

slide-59
SLIDE 59

ISM: Density Fluctuations

Electron-density fluctuations in the interstellar medium [Armstrong et al. 1995, ApJ 443, 209]

k–5/3 “Great Power Law In the Sky”

… coined by Steve Spangler

slide-60
SLIDE 60

SW: Density and Field-Strength Fluctuations

[Bershadskii & Sreenivasan 2004, PRL 93, 064501] Spectrum of magnetic-field strength in the solar wind at ~1 AU (1998) Density fluctuations in the solar wind at ~1 AU (31 Aug. 1981) [Celnikier, Muschietti & Goldman1987, A&A 181, 138]

k–5/3

FLR: density mode mixing with Alfvén waves

slide-61
SLIDE 61

KRMHD: Density and Magnetic-Field Strength

Density and field-strength fluctuations are passively mixed by Alfvén waves require kinetic description: our expansion gives

Maxwellian equilibrium

KRMHD

[Schekochihin et al., arXiv:0704.0044

  • cf. Higdon 1984, ApJ 285, 109; Lithwick & Goldreich 2001, ApJ 562, 279]
slide-62
SLIDE 62

KRMHD: Density and Magnetic-Field Strength

require kinetic description: our expansion gives In the Lagrangian frame of the Alfvén waves…

[Schekochihin et al., arXiv:0704.0044]

slide-63
SLIDE 63

KRMHD: Density and Magnetic-Field Strength

require kinetic description: our expansion gives In the Lagrangian frame of the Alfvén waves… equation is linear!

[Schekochihin et al., arXiv:0704.0044]

slide-64
SLIDE 64

KRMHD: Density and Magnetic-Field Strength

require kinetic description: our expansion gives In the Lagrangian frame of the Alfvén waves… equation is linear! No refinement of scale along perturbed magnetic field (but there is along the guide field, i.e. kz grows)

slide-65
SLIDE 65

Collisionless Damping

require kinetic description: our expansion gives equation is linear!

[Barnes 1966, Phys. Fluids 9, 1483]

time to be cascaded in k⊥ by Alfvén waves, for which

Cascades of density and field strength fluctuations are undamped above ion gyroscale

… but parallel cascade might be induced due to dissipation [Lithwick & Goldreich 2001, ApJ 562, 279]

slide-66
SLIDE 66

Damping of Cascades

k–5/3 k–7/3

If they have a parallel cascade, density and field strength are damped Alfvén waves Landau damped via conversion into density/field-strength fluctuations arXiv:0704.0044

ion Landau damping electron Landau damping Barnes damping KAW energy injected

slide-67
SLIDE 67

Damping of Cascades

k–5/3 k–7/3

If their parallel cascade is inefficient, density and field strength are only weakly damped above ρi Alfvén waves Landau damped via conversion into density/field-strength fluctuations arXiv:0704.0044

ion Landau damping electron Landau damping KAW energy injected

slide-68
SLIDE 68

Damping of Cascades

k–5/3 k–7/3

If their parallel cascade is inefficient, density and field strength are only weakly damped above ρi Alfvén waves Landau damped via conversion into density/field-strength fluctuations arXiv:0704.0044

ion Landau damping electron Landau damping KAW energy injected

SO, IDEA #6: PASSIVE COMPRESSIVE MODES IN THE INERTIAL RANGE WITH NO PARALLEL CASCADE?

slide-69
SLIDE 69

Electron Reduced MHD

k–5/3 k–7/3

arXiv:0704.0044

ion Landau damping electron Landau damping energy injected

Boltzmann ions magnetised electrons (still isothermal)

slide-70
SLIDE 70

Electron Reduced MHD

This is the anisotropic version of EMHD [Kingsep et al. 1990,

  • Rev. Plasma Phys. 16, 243],

which is derived (for βI >>1) by assuming magnetic field frozen into electron fluid and doing a RMHD-style anisotropic expansion:

Start with GK, consider the scales such that

arXiv:0704.0044

slide-71
SLIDE 71

Kinetic Alfvén Waves

Linear wave solutions:

Eigenfunctions:

  • Critical balance + constant flux argument à la K41/GS95 give

spectrum of magnetic field with anisotropy

  • There is a cascade of KAW,
  • Electric field has spectrum:

Start with GK, consider the scales such that

arXiv:0704.0044 [Biskamp et al. 1996, PRL 76, 1264; Cho & Lazarian 2004, ApJ 615, L41]

slide-72
SLIDE 72

Kinetic Alfvén Waves

Linear wave solutions:

Eigenfunctions:

Start with GK, consider the scales such that

arXiv:0704.0044

SO, IDEA #7: CRITICALLY BALANCED KAW CASCADE IN THE DISSIPATION RANGE

slide-73
SLIDE 73

Dissipation Range of the SW: KAW?

Magnetic- and electric-field fluctuations in the solar wind at ~1 AU (19 Feb. 2002) [Bale et al. 2005, PRL 94, 215002]

k–5/3 Alfvén waves k–7/3 k–1/3 KAW

slide-74
SLIDE 74

Dissipation Range of the SW: No KAW?

Magnetic-field fluctuations in the solar wind at ~1 AU (19 Feb. 2002) [Leamon et al. 1998, JGR 103, 4775]

slide-75
SLIDE 75

Dissipation Range of the SW: ???

Spectral indices in the inertial and dissipation ranges [Smith et al. 2006, ApJ 645, L85]

slide-76
SLIDE 76

Nonlinear Perpendicular Phase Mixing

arXiv:0806.1069

IDEA #8: DUAL (ION) ENTROPY CASCADE IN VELOCITY AND POSITION SPACE

slide-77
SLIDE 77

Nonlinear Perpendicular Phase Mixing

This comes from gyroaveraging

NB: In fluid models (like EMHD) these fluctuations are invisible Low-frequency electrostatic fluctuations

arXiv:0806.1069

slide-78
SLIDE 78

Nonlinear Perpendicular Phase Mixing

  • Potential mixes hi via this term,

so hi developes small (perpendicular) scales in the gyrocenter space: Low-frequency electrostatic fluctuations

arXiv:0806.1069

slide-79
SLIDE 79

Nonlinear Perpendicular Phase Mixing

  • Potential mixes hi via this term,

so hi developes small (perpendicular) scales in the gyrocenter space:

  • Two values of the gyroaveraged potential

come from spatially decorrelated fluctuations if

°

  • Low-frequency

electrostatic fluctuations

arXiv:0806.1069 [The perpendicular nonlinear phase-mixing mechanism was anticipated in the work of Dorland & Hammett 1993]

slide-80
SLIDE 80

Entropy Cascade

  • Electrostatic fluctuations come from ion-entropy fluctuations:
  • Entropy is conserved, so use const-flux argument:
  • Nonlinear decorrelation time:

Low-frequency electrostatic fluctuations

arXiv:0806.1069

slide-81
SLIDE 81

Entropy Cascade

  • Electrostatic fluctuations come from ion-entropy fluctuations:
  • Entropy is conserved, so use const-flux argument:
  • Nonlinear decorrelation time:

Low-frequency electrostatic fluctuations

arXiv:0806.1069

slide-82
SLIDE 82

Entropy Cascade

We get the following set of scaling relations:

arXiv:0806.1069

Low-frequency electrostatic fluctuations

slide-83
SLIDE 83

Entropy Cascade: GK 4D DNS by T. Tatsuno

arXiv:0806.1069

2562×722

slide-84
SLIDE 84

Entropy Cascade: GK 4D DNS by T. Tatsuno

arXiv:0806.1069

2562×722

Similar (density) spectra also reported in 3D ITG/ETG tokamak flux-tube GK simulations by Görler & Jenko (2008)

slide-85
SLIDE 85

Entropy Cascade: GK 4D DNS by T. Tatsuno

Distribution function develops small-scale structure in velocity space

arXiv:0806.1069

slide-86
SLIDE 86

Entropy Cascade: GK 4D DNS by T. Tatsuno

Distribution function develops small-scale structure in velocity space

arXiv:0806.1069

slide-87
SLIDE 87

Entropy Cascade: GK 4D DNS by T. Tatsuno

Distribution function develops small-scale structure in velocity space

arXiv:0806.1069

slide-88
SLIDE 88

Entropy Cascade: GK 4D DNS by T. Tatsuno

Distribution function develops small-scale structure in velocity space

arXiv:0806.1069

  • G. Plunk has developed a “kinematics of phase-space turbulence”

to quantify perpendicular velocity-space structure via Hankel transforms and derived scaling relations à la K41

slide-89
SLIDE 89

Entropy Cascade: GK 4D DNS by T. Tatsuno

Distribution function develops small-scale structure in velocity space

arXiv:0806.1069

  • G. Plunk has developed a “kinematics of phase-space turbulence”

to quantify perpendicular velocity-space structure via Hankel transforms and derived scaling relations à la K41

slide-90
SLIDE 90

Phase-Space Cutoff

arXiv:0806.1069

Distribution function develops small-scale structure in velocity space

slide-91
SLIDE 91

Phase-Space Cutoff

Distribution function develops small-scale structure in velocity space

arXiv:0806.1069 characteristic time at the ion gyroscale

slide-92
SLIDE 92

Phase-Space Cutoff

Distribution function develops small-scale structure in velocity space

arXiv:0806.1069

Do–3/5

characteristic time at the ion gyroscale

Dorland Number x- and v-space resolution are related

  • cf. kcL ~ Re3/4 in Kolmogorov fluid turbulence
slide-93
SLIDE 93

Linear Parallel Phase Mixing

Parallel phase mixing is due to the “ballistic response”:

arXiv:0806.1069

if linear propagation time ~ nonlinear decorrelation time (“critical balance”) So the nonlinear perpendicular phase mixing dominates after t ~ τλ

slide-94
SLIDE 94

Dissipation Range With and Without KAW

With KAW Without KAW High-frequency, electromagnetic, fluid-like (EMHD) Low-frequency, electrostatic, purely kinetic (GK ions)

slide-95
SLIDE 95

Dissipation Range With and Without KAW

k–5/3 Alfvén waves k–7/3 k–1/3 KAW [Leamon et al. 1998, JGR 103, 4775] [Bale et al. 2005, PRL 94, 215002]

With KAW Without KAW High-frequency, electromagnetic, fluid-like (EMHD) Low-frequency, electrostatic, purely kinetic (GK ions)

slide-96
SLIDE 96

Dissipation Range of the Solar Wind

Variable spectral index in the dissipation range may be due to superposition of KAW and no KAW cascades With KAW Without KAW

[Smith et al. 2006, ApJ 645, L85] arXiv:0704.0044

slide-97
SLIDE 97

Dissipation Range of the Solar Wind

Variable spectral index in the dissipation range may be due to superposition of KAW and no KAW cascades With KAW Without KAW

[Smith et al. 2006, ApJ 645, L85] arXiv:0704.0044

SO, THIS WAS IDEA #9: WE MAY BE OBSERVING THE ENTROPY CASCADE