Turbulence and Transport in Magnetically Confined Plasma Jens Juul - - PowerPoint PPT Presentation

turbulence and transport in magnetically confined plasma
SMART_READER_LITE
LIVE PREVIEW

Turbulence and Transport in Magnetically Confined Plasma Jens Juul - - PowerPoint PPT Presentation

Turbulence and Transport in Magnetically Confined Plasma Jens Juul Rasmussen Jens Madsen, Volker Naulin, Anders H. Nielsen Association EURATOM DTU, Department of Physics, Technical University of Denmark, DTU Ris Campus, DK-4000 Roskilde


slide-1
SLIDE 1

Department of Physics Technical University of Denmark

Turbulence and Transport in Magnetically Confined Plasma

Jens Juul Rasmussen

Jens Madsen, Volker Naulin, Anders H. Nielsen

Association EURATOM – DTU, Department of Physics, Technical University of Denmark, DTU Risø Campus, DK-4000 Roskilde

Odd Erik Garcia

Department of Physics and Technology, University of Tromsø, N-9037 Tromsø jjra@fysik.dtu.dk

slide-2
SLIDE 2

05-06-2012 2 DTU Physics SCT-2012

Outline/Motivation

  • Turbulence and transport in plasma – particularly in the

edge region of magnetically confined plasma

  • Results from numerical modelling with comparison with

experimental observations

  • In magnetically confined plasmas the anomalous –

turbulent - transport is the dominant mechanism for transport of particles and energy across the confining magnetic field – orders of magnitude higher than classical collisional transport

  • Understanding and predicting the transport is essential

for the viable operation of fusion power plants

slide-3
SLIDE 3

05-06-2012 3 DTU Physics SCT-2012

Tokamak configuration – edge/SOL

SOL Conditions at the edge determine plasma performance

slide-4
SLIDE 4

05-06-2012 4 DTU Physics SCT-2012

Turbulence and transport

  • large-amplitude, radially propagating blob-like structures of

particles and heat, generated close to the last closed flux surface

  • results in localized power loads at plasma facing components
  • lasting influence on the chamber wall and other plasma facing

components

  • strong demands on materials

Observed under a variety of conditions:

e.g., Zweben Phys. Fluids 28 974 (1985); Boedo et al PoP 10, 1670 (2003); Zweben et al. Nucl. Fus. 44, 134 (2004); Grulke et al PoP 13, 012306 (2006); Garcia et al. PPCF 48, L1 (2006); Garcia et al. PPCF 49, B47 (2007); Xu et al Nucl. Fus. 49, 092002 (2009); Nold et al.PPCF 52, 065005 (2010). ...

Reviews: Zweben et al Plasma Physics Control Fusion 49, S1 (2007); D’Ippolito et al. Phys Plasmas 18, 060501 (2011)

Cross field transport of particles and heat in magnetically confined plasmas is dominated by anomalous - turbulent - transport! In the edge/scrape-off-layer (SOL) region transport is strongly intermittent and characterized by:

slide-5
SLIDE 5

05-06-2012 5 DTU Physics SCT-2012

Typical density fluctuations in JET – Joint European Torus

  • G. Xu et al Nucl Fusion 49, 092002 (2009)

In the SOL Inside last closed flux surface

Intermittent density fluctuations

slide-6
SLIDE 6

05-06-2012 6 DTU Physics SCT-2012

Blob propagation, example from NSTX

2cm

2cm 208.853 ms

Separatrix Antenna limiter shadow

Shot 118152

Plasma blobs – field aligned filaments - detaches from confinement region and propagate into the SOL, where the plasma can flow towards material surfaces.

Zweben et al. Phys. Plasma 13, 056114 (2006); PPCF 49, S1 (2007). D’Ippolito et al. Phys. Plasmas 18, 060501 (2011) Maqueda et al. IEA Workshop Edge Transport in Fusion Plasmas, Sept. 11 – 13, 2006, Krakow, Poland

slide-7
SLIDE 7

05-06-2012 7 DTU Physics SCT-2012

Blob: parallel structure

Magnetic field line

Blobs are filaments stretched long the magnetic field lines

Quasi 2D dynamics

Grulke et al Phys. Plasma 13, 012306 (2006). Alcator C-Mod

slide-8
SLIDE 8

05-06-2012 8 DTU Physics SCT-2012

Propagation of a blob of plasma

Consider a localized perturbation of the plasma density/pressure at the

  • utboard side

Charge separation by gradB and curvature drifts sets up a vertical E- field, resulting in a radial ExB velocity Toroidal magnetic field curved and weaker on the outboard side Radially propagating filament

slide-9
SLIDE 9

05-06-2012 9 DTU Physics SCT-2012

Blob propagation

Initially a density perturbation and there is no flow; this arises by vertical charge polarization An isolated blob/filament accelerate and propagate a large radial distance Eventually the blob decelerates and disperses Density Vorticity

Garcia et al. Phys. Plasmas 12, 090701 (2005); 13, 082309 (2006)

slide-10
SLIDE 10

05-06-2012 10 DTU Physics SCT-2012

2D Equations for convection – fluid model

“effective gravity” Ω – vorticity φ- stream function θ- electron pressure (cold ions)

Rayleigh-Taylor type instability

“Minimal” model!

slide-11
SLIDE 11

05-06-2012 11 DTU Physics SCT-2012

Finite ion temperature – finite ion Larmor radius effect

Important when blob size σ< 10ρT ion Larmor radius Ions then responds to a potential averaged over the fast gyro motion: ρT resulting in a different ExB velocity than the electrons. leads to the derivation of the Gyro- fluid model – “lowest order kinetic effects”

slide-12
SLIDE 12

05-06-2012 12 DTU Physics SCT-2012

Madsen et al Phys. Plasma 18, 112504 (2011)

Gyro-fluid simulation of blob propagation

Ti/T

e = 0

Ti/T

e = 3

Propagation of density blob with finite ion temperature effects – compact density blob – like experiments -

slide-13
SLIDE 13

05-06-2012 13 DTU Physics SCT-2012

Finite ion temperature blob propagation

Thermal energy IC carried by blob

Finite ion temperature structures travel further and live longer. Broaden energy deposition (good), but may run into vessel wall (bad).

Madsen et al Phys. Plasma 18, 112504 (2011)

Blob size σ<5 ρi ion Larmor radius : concentrated blobs

slide-14
SLIDE 14

05-06-2012 14 DTU Physics SCT-2012

Blobs in turbulence – the ESEL model

A self-consistent description of fluctuations and intermittent transport in the edge/SOL by employing the ESEL ( Edge SOL Electrostatic) model for interchange dynamics. 2D model perpendicular to magnetic field – dynamics/losses along magnetic field are parameterized

  • include separate plasma production ``edge'' and loss region

``SOL'',

  • allow self-consistent flows and profile relaxations,
  • profiles and fluctuations are NOT separated,
  • conserve particles and energy in collective dynamics.

Results agree very well with experimental observations! Being applied at many laboratories.

Garcia, Naulin, Nielsen, Rasmussen, PRL 92 165003 (2004); Phys. Plasma 12, 090701 (2005); Physica ScriptaT122, 89 (2006).

slide-15
SLIDE 15

05-06-2012 15 DTU Physics SCT-2012

Slab domain at

  • utboard midplane
  • Local slab 2D geometry, (x,y)
  • Including edge and SOL
  • Global model with self-consistent profiles

B 

B 

Model domain: ESEL

slide-16
SLIDE 16

05-06-2012 16 DTU Physics SCT-2012

Energy evolution and transfer

slide-17
SLIDE 17

05-06-2012 17 DTU Physics SCT-2012

Spatial structure during burst - developing blob structure

Particle density (left) and vorticity (right) during a burst (Δt = 500) Blob like-structure in plasma density and dipole structure in vorticity

slide-18
SLIDE 18

05-06-2012 18 DTU Physics SCT-2012

Single point PDF, density fluctuations

Scaled probability density distribution functions, PDF, of density fluctuations at Pi . i > 2 – inside SOL - exponential tails, indicating strong blob structures. – Fit a Gamma distribution i >2 Skewness, S, and flatness, F, factors for the density fluctuations

slide-19
SLIDE 19

05-06-2012 19 DTU Physics SCT-2012

Particle density flux

Re-scaled PDF of particle density flux measured at the probes, Pi. Exponential tails: flux dominated by strong bursts. Re-scaled PDF of the turbulent radial ExB-velocity, vx, recorded at the probes Pi .

Turbulent particle density flux : Γ = n vx = n vExB

Transport is NOT diffusive

slide-20
SLIDE 20

05-06-2012 20 DTU Physics SCT-2012

Density fluctuations statistics and wave form

Conditionally averaged density wave form in far SOL Characteristics of blob propagation and fits well simulation results Rescaled PDFs of density fluctuations in far SOL Skewed to the positive side. Fits well simulation results.

Garcia et al. PPCF 48, L1 (2006); Nucl. Fusion 47, 667 (2007)

Direct comparison with experimental results from the TCV- Tokamak, Lausanne: excellent quantitative agreement

slide-21
SLIDE 21

05-06-2012 21 DTU Physics SCT-2012

Particle density flux statistics

Rescaled PDF of particle flux in far SOL at TCV Tokamak. Almost independent of ne. Flux dominated by strong bursts and agreement with simulation results. Exponential tail -- mean value only contain limited information

Transport is NOT diffusive! No simple parameterization in terms of an effective diffusivity and a convection velocity

Garcia et al., Nucl. Fusion 47, 667 (2007)

Particle density flux: Γ = n vx = n vExB

slide-22
SLIDE 22

05-06-2012 22 DTU Physics SCT-2012

Scatter plot for the flux-gradient relation.

Parametrization of density flux?

Transport cannot be parameterized by an effective diffusivity and a convection velocity

Garcia et al., J. Nucl. Mater. 263-265, 575 (2007) Naulin, J. Nucl. Mater. 263-265, 24 (2007)

TCV ESEL

Transport modelling: linear combination of convection and diffusion:

slide-23
SLIDE 23

05-06-2012 23 DTU Physics SCT-2012

Summary

  • Transport of particles and heat in the edge region of

magnetically confined plasmas is dominated by large- amplitude, radially propagating blob-like structures of particles and heat

  • No parameterisation of transport; not Fickian diffusion.

The flux is given by its PDF, characterized by “fat” tails

  • Experimental results well described by the ESEL-model
  • Blob structures give rise to localized power loads on

plasma facing components

  • Strong demands on the materials and sets engineering

limits to power plant operation

  • For ITER, events with power loads of several MW/cm2 are

expected and control will be essential

slide-24
SLIDE 24

05-06-2012 24 DTU Physics SCT-2012

Thanks for your attention