• Keine Ergebnisse gefunden

Plasma Facing Components (PFC)

N/A
N/A
Protected

Academic year: 2022

Aktie "Plasma Facing Components (PFC)"

Copied!
76
0
0

Wird geladen.... (Jetzt Volltext ansehen)

Volltext

(1)

S. W. Lisgo 1

with the PFC part based heavily on presentations by

F. Escourbiac 1 , T. Hirai 1

and, in general, helped by many conversations with

G. De Temmerman

1

and additional material provided by

S. Carpentier-Chouchana 1 , D. Coster 2 , A. Dvornova 1 , J. Gunn 3 , J. Harrison 5 , M. Kocan 1 , S. Panayotis 1 , R.A. Pitts 1 , J. Terry 6 ,

and the ITER Divertor Design Team 1,4

Plasma-Wall Interaction (PWI) and Plasma Facing Components (PFC) Plasma-Wall Interaction (PWI) and Plasma Facing Components (PFC)

The views and opinions expressed herein most certainly do not reflect those of the ITER Organization, at all

1

ITER Central Team, France

2

IPP Garching, Germany

3

IRFM/CEA, France

4

ITER Domestic Agencies, Earth

5

CCFE, UK

6

PSFC, MIT, USA

(2)

 Boundary plasma transport

 Plasma-material interaction

 Plasma facing components

(3)

The (ITER) tokamak (the “model tokamak” for this presentation)

(4)

The (ITER) tokamak

cryostat

(5)

The (ITER) tokamak

vacuum vessel port extensions:

diagnostics, heating, fueling, pumps, in-vessel component feeds, remote access

(6)

The (ITER) tokamak

(7)

The (ITER) tokamak

toroidal field coils

(8)

The (ITER) tokamak

vacuum vessel

blanket

divertor

central solenoid

(9)

The (ITER) tokamak

(10)

Magnetic field structure: toroidal field

(11)

Magnetic field structure: poloidal field

B poloidal  0.1 B toroidal

plasma current

(12)

Magnetic field structure: magnetic field line trajectory (q = 3)

B poloidal  0.1 B toroidal

(13)

Magnetic field structure: magnetic field line trajectory (q = 3)

MACHINE GEOMETRY MAGNETIC GEOMETRY

B poloidal  0.1 B toroidal

(14)

Magnetic field structure: x-point geometry directs plasma toward divertor

(15)

Magnetic field structure: “open” field lines intersect solid surfaces

(16)

Magnetic field structure: strong radial magnetic shear near the x-point

(17)

Plasma transport from the core to the boundary: turbulence

gravity

FLUID DYNAMICS: RALEIGH-TAYLOR INSTABILITY ( link)

(18)

Plasma transport from the core to the boundary: turbulence

C-MOD

gravity

EB

HIGH SPEED D a GAS PUFF IMAGING PLASMA DYNAMICS:

“INTERCHANGE” INSTABILITY

[J. Terry]

DEUTERIUM

GAS PUFF

(19)

Plasma transport from the core to the boundary: turbulence

C-MOD

1 cm

gravity

EB

HIGH SPEED D a GAS PUFF IMAGING PLASMA DYNAMICS:

“INTERCHANGE” INSTABILITY

[J. Terry]

fluctuations can also confine

plasma DEUTERIUM

GAS PUFF

(20)

Transport from the core to the boundary: ballooning transport

MAX = 4.310

19

ph m

-3

s

-1

INVERTED He +1 @ 468 nm IMAGE

5 mm mesh, 500500 image, shot 11849

MAST LIGHT FROM SINGLY CHARGED HELIUM

[S. Lisgo]

(21)

Main chamber recycling (MCR): plasma contact with main chamber surfaces

[A. Loarte 2007]

FAR-SOL FAR-SOL NEAR-SOL

NEAR-SOL

separat rix

RADIAL PLASMA DENSITY PROFILE

(22)

Main chamber recycling (MCR): plasma contact with main chamber surfaces

[A. Loarte 2007]

FAR-SOL FAR-SOL NEAR-SOL

NEAR-SOL

separat rix

RADIAL PLASMA DENSITY PROFILE

ESTIMATED RADIAL PLASMA VELOCITY

(23)

Turbulence suppression: H-mode MAST

15622, unfiltered

TURBULENT BOUNDARY ELMing H-MODE

(24)

Turbulence suppression: H-mode MAST

15622, unfiltered

TURBULENT BOUNDARY ELMing H-MODE

(25)

[Eich 2013]

Power exhaust: peak steady-state heat-flux to the divertor targets

l q,|| for ITER currently predicted to be 1 mm, significantly smaller than the design value of 5 mm  reduced operational

space? [Kukushkin 2015]

S scaling? [Scarabosio 2015]

l int

l int  radial heat-flux decay length at target

l q,||

l q,||  radial decay length of heat-flux into divertor

S

S  heat-flux spreading due to divertor transport

l int  l q,|| + 1.64 S

(26)

Power exhaust: cross-field turbulent transport in the divertor (related to S)

LOCALIZED DIVERTOR TURBULENCE

X-POINT RADIAL MAGNETIC SHEAR

(27)

Power exhaust: ELM transient heat loads to PFCs

 Electron beam thermal cycling of ITER mono-block components indicate failure at

~0.5 MJ m -2 for both CFC and W [J. Linke et al JNM (2007)]

 Natural ELMs in ITER predicted to be 10-20 MJ m -2  ELM control required

(28)

Power exhaust: ELM control with Resonant Magnetic Perturbation (RMP) EMC3 SIMULATION OF DIII-D

MAST

He + , NO RMP

WITH RMP RMP

COILS

(29)

Tungsten influx during an ELM and core-boundary impurity coupling SOLPS-DIVIMP ELM SIMULATION

W DENSITY, SUM OVER ALL STATES

TARGET ELECTRON TEMPERATURE

[D. Coster 2015, S. Lisgo]

(30)

Impurities can be a very good thing: control actuator for dissipation

INFRARED, VIEW OF DIVERTOR FROM ABOVE [A. Kallenbach]

ASDEX UPDGADE

(31)

 Boundary plasma transport

 Plasma-material interaction

 Plasma facing components

(32)

Plasma-surface interaction in the divertor

ION FROM “UPSTREAM” ACCELERATING TOWARD THE TARGET

Angle between the magnetic field, B, and the surface is typically a few degrees

At 20 eV and 5 Tesla, the gryo-radius for a deuterium ion is

0.1 mm

(33)

The sheath

ELECTROSTATIC DEBYE SHEATH NEAR THE SURFACE

Electrostatic potential drop (acceleration) experienced by an ion passing through the sheath:

V sh  2 T i + 3 Z ion T e At T e = 20 eV and n e = 10 19 m -3 , the sheath thickness is

0.01 mm

(34)

Plasma recycling

PLASMA RECYCLING, FLUX AMPLIFICATION

(35)

Plasma recycling

PLASMA RECYCLING, FLUX AMPLIFICATION

ion flux  n e,upstream

when ionisation near the surface (local to the divertor)

becomes important

2

(36)

Charge exchange

CHARGE EXCHANGE (CX)

For hydrogen, the CX rate always larger than the ionisation rate

Important process for the transport of

energy across the

magnetic field, i.e. CX

(37)

Physical sputtering

PHYSCIAL SPUTTERING

4 𝑀 1 𝑀 2 (𝑀 1 + 𝑀 2 ) 2 Maximum energy

fraction transferred in a head-on collision:

= 100% for M 1 = M 2

and 4% for D  W

(38)

Physical sputtering

SPUTTERING YIELD PER INCIDENT DEUTERIUM ION [Stangeby pg 119]

D threshold T

C W

W self-sputtering yield > 1

(normal incidence)

implantation

(39)

Prompt re-deposition

GROSS AND NET EROSION, PROMPT RE-DEPOSITION

Net erosion can be much less than gross erosion

For example, 4 orders of magnitude lower for W sputtering during ELMs 

prompt re-deposition

[Chankin 2014]

(40)

Temperature gradients in the boundary

PARALLEL TEMPERATURE GRADIENTS, HIGH RECYCLING DIVERTOR

Energy dissipation in the boundary due to recycling and

impurity radiation can lead to strong temperature

gradients between

the core and divertor

(41)

Plasma “detachment”

“GAS TARGET” DIVERTOR, STRONG DISSPATION

For T e < 5 eV, plasma- surface interaction can be strongly reduced 

“detachment”

“Deep” detachment

can affect core plasma

performance

(42)

BALMER GAMMA EMISSION  T e < 1.5 eV

[C. Boswell 2001]

An example of a strongly dissipating divertor

C-MOD (SCALE 10)

ITER

(43)

Many open questions related to plasma-material interaction MANY ACTIVE AREAS OF RESEARCH

• mixed materials

• weakly bound atoms

• chemical erosion

• surface roughness

• incident angle

• plasma chemistry

• heavy ions

• plasma transients

• magnetic pre-sheath

(44)

ITER issue: tritium retention by co-deposition

 A 400 s Q=10 pulse will require ~50 g of T fuel, but the maximum mobilisable in- vessel T inventory is limited to 640 g (+180 g in pumps, +180 g uncertainty)

− nuclear safety issue (expensive too); tritium burn-up is only ~0.3% in ITER

T:Be trapped fraction

~1:20 @ 200 C

[G. De Temmerman, R. P. Doerner, et al., NF (2009)]

T:Be trapped fraction

~1:20 @ 200 C

[G. De Temmerman, R. P. Doerner, et al., NF (2009)]

T:Be <1:100 @ 300+ C T:Be <1:100 @ 300+ C

 T:Be depends sensitively on deposition rate, incoming particle energy, and surface

temperature  complex problem  on-going efforts to predict the level of T-

(45)

ITER issue: fluence

 Long-pulse, large size, and high density operation combine to give a significant increase in the ion fluence to the wall

No.

Pulses

Time in Diverted Phase (hours)

Outer Divertor Ion Fluence

JET 13466 40.5 ~5x10 27 [1]

ITER (Q=10)

No.

Pulses

Time in Diverted Phase (hours)

Outer Divertor Ion Fluence

JET 13466 40.5 ~5x10 27 [1]

ITER (Q=10) 1 0.15 ~1.5x10 27 [2]

DIVERTOR ION FLUX COMPARISON WITH JET (2000-2008, campaigns C1-C19)

[1] M.F. Stamp, CCFE, private communication

[2] SOLPS4.3 code result: run #1514, A. S. Kukushkin

 9 years of JET operation ≡ 3 ITER pulses @ Q=10 (~1.5 hours of real time)

− 3 decades of JET, or ~half the time humanity has spent on controlled fusion, in a morning

 The high resulting material turnover affects tritium retention, material mixing, layer

growth, and dust production

(46)

ITER issue: surface morphology: helium and tungsten [G. De Temmerman]

 Want to commission ITER systems during the non-nuclear phase (H and He plasmas), but H-mode may only be accessible with He plasmas

Nano-bubbles ‘Fuzz’

~ 600 - 700 K ~ 900 – 1900 K

500nm [M. Miyamoto et al, JNM, 415 (2011)]

> 2000 K

Large voids

1mm

 Helium insoluble in tungsten but can self-trap

 Clustering of helium to form

bubbles  material ‘swell’

(47)

 Boundary plasma transport

 Plasma-material interaction

 Plasma facing components

(48)

How can PFCs survive a steady-state plasma environment?

 Most tokamaks have been “short pulse” (~seconds)  PFCs heat up during a discharge and cool down before the next one  “inertial cooling”

− main thermo-physical property is the specific heat capacity of the material, C

p

(J Kg

-1

K

-1

)

− armour tiles held in place mechanically with bolts (C / Mo reliable up to ~2500 C)

 Long pulse (superconducting)  active cooling  complex PFCs

thermal conductivity, k (W m

-1

K

-1

), is the most important material property (a function of temperature and neutron fluence)

bonded components that are welded to the structure (bolts OK below 1-2 MW m

-2

)

HEAT

COOLANT (H 2 O) MAIN STRUCTURAL MATERIAL (316L(N)-IG SS)

PLASMA FACING MATERIAL, PFM (Be, W)

HIGH l STRUCTURAL MATERIAL (CuCrZr)

INTERLAYER (Cu)

(49)

The power of active cooling

 Boiling water in a paper cup…

− (warning: enthalpy of vapourization important here – not directly applicable to ITER PFCs)

[L. Jedral and his high school physics class, Ontario Public School Board]

(50)

Time-averaged heat flux design limit at the target: 10 MW m -2

 What does 10 MW m -2 look like?

 If there is no energy dissipation in the ITER boundary plasma (very unlikely), then for the ITER baseline assumption of l q,|| = 5 mm, there could be 40 MW m -2 at the

H EA T FL UX DE NSI TY (M W m -2 )

DISTANCE FROM CENTER OF ARC (mm)

[A. B. Murphy CSIRO]

HEAT FLUX TO THE WORKPIECE FOR AN ARC WELDER

[http://www.youtube.com/watch?v=Ng-urnWhBR8&hd=1]

40 MW m -2

10 MW m -2

(51)

Divertor cassettes

 54 cassettes, 8.7 tonnes each, “easily” removable

− 6 months to replace

(52)

Divertor cassettes

(53)

Divertor cassettes

(54)

Divertor cassettes

(55)

Vertical target and plasma facing unit (PFU) configuration

PLASMA FACING UNIT (PFU)

water flowing in a tube

CuCrZr

 A water cooled tube made from a high thermal

conductivity metal can remove heat very efficiently

(56)

Vertical target and plasma facing unit (PFU) configuration

PLASMA FACING UNIT (PFU)

 Increase tube thickness for greater strength and to maximize lifetime (tokamak plasma contact will erode the tube)

CuCrZr

(57)

Vertical target and plasma facing unit (PFU) configuration

PLASMA FACING UNIT (PFU)

W (6-8 mm thick) Cu

 A good structural metal may not be compatible with plasma interaction, so armour the tube with an

erosion resistant and high melting point material such as tungsten (note: tungsten is a poor structural material at tube temperatures appropriate for water cooling)

CuCrZr

 A layer of copper (soft) compensates for the thermal

expansion mismatch between the armour and tube

(58)

Vertical target and plasma facing unit (PFU) configuration

PLASMA FACING UNIT (PFU)

Cu

 Castellate to reduce thermomechanical stress in the armour and on the bond with the tube, as well as electromagnetic forces during plasma disruptions

W (6-8 mm thick)

CuCrZr

(59)

Vertical target and plasma facing unit (PFU) configuration

“VAPOUR MATTRESS”

• heat transfer dramatically reduced

• “critical heat flux” (CHF)  imminent failure of component

“NUCLEATE BOILING”

• onset at T sat

• heat transfer coefficient goes up by a factor 3-8

[SOLPS, A. Kukushkin]

[RACLETTE, A. Sashala Naik]

(60)

Vertical target and plasma facing unit (PFU) configuration

T wall “VAPOUR

MATTRESS”

• heat transfer dramatically reduced

• “critical heat flux” (CHF) 

(61)

Vertical target and plasma facing unit (PFU) configuration

T wall “VAPOUR

MATTRESS”

• heat transfer dramatically reduced

• “critical heat flux” (CHF)  imminent failure of component

 The CHF is a function of the overall system design and can be

raised by increasing the coolant velocity + pressure, the PFC

complexity, and by sub-cooling  more money

(62)

Vertical target and plasma facing unit (PFU) configuration

 The CHF is a function of the overall system design and can be raised by increasing the coolant velocity + pressure, the PFC complexity, and by sub-cooling  more money

“SWIRL TUBE”

(63)

 The CHF is a function of the overall system design and can be raised by increasing the coolant velocity + pressure, the PFC complexity, and by sub-cooling  more money

Vertical target and plasma facing unit (PFU) configuration

ICHF VERSUS COOLANT FLOW VELOCITY

[A.R. Raffray, 1999]

(64)

 The CHF is a function of the overall system design and can be raised by increasing the coolant velocity + pressure, the PFC complexity, and by sub-cooling  more money

Vertical target and plasma facing unit (PFU) configuration

ICHF VERSUS COOLANT FLOW VELOCITY

[A.R. Raffray, 1999]

SCREW TUBE SCREW TUBE SCREW TUBE

Screw tube not chosen despite higher CHF  concerns that the grooves will seed crack formation in the tube

[A.R. Raffray, Fus. Eng. and Design 45 (1999) 377–407]

(65)

Vertical target and plasma facing unit (PFU) configuration

SMOOTH TUBE

SWIRL TUBE WITH TWISTED TAPE

HIGHER CHF

CuCrZr Cu

W

(66)

The dome

W FLAT TILE

ARMOUR

HYPER-

VAPOTRON

(67)

Hypervapotron

NUCLEATE BOILING IMAGES (low T)

[P. Chen, PhD thesis, 2007]

 Good CHF (to ~20 MW m -2 so far) with little loss in pressure  turbulence mixing of coolant across the square channels

− limitation perhaps from structural integrity at present (thermal gradients), rather than CHF

[J. Milnes, PhD thesis, 2011]

EDDY SIMULATIONS

[O. Aybay, Fluid Mechanics and Multidisciplinary Design, 2010]

FLOW

(68)

Exposed “leading edges” can receive the full parallel heat-flux

 Magnetic field lines can penetrate gaps between

monoblocks, allowing the plasma to heat the “leading edges”

of the PFCs  power is concentrated onto a smaller area, resulting in strong local heating

a

2D 2L

(69)

Exposed “leading edges” can receive the full parallel heat-flux

 Magnetic field lines can penetrate gaps between

monoblocks, allowing the plasma to heat the “leading edges”

of the PFCs  power is concentrated onto a smaller area, resulting in strong local heating

a

2D 2L

(70)

Exposed “leading edges” can receive the full parallel heat-flux

 Magnetic field lines can penetrate gaps between

monoblocks, allowing the plasma to heat the “leading edges”

of the PFCs  power is concentrated onto a smaller area, resulting in strong local heating

a

2D 2L

 The shallow angle between the magnetic field lines and the

surface, due to the large toroidal component of the field and

the vertical target geometry, means that small mono-block

(71)

Misalignments accounted for in the manufacturing specifications

 Geometrical tolerances are specified considering operational requirements at global and local scales

− all monoblocks aligned within 0.5 mm range

− steps of neighbouring blocks smaller than 0.3 mm

 To achieve the required tolerance during manufacture:

− identify poor tolerance operation: joining (welding) and assembly

− recover by higher tolerance operation (machining) and custom

machining of tolerance compensators

(72)

 Monoblock edges must be protected in an all-W divertor (C is more “forgiving”  “plasma polishing”)

− straight chamfer at 0.45º angle protects 0.3 mm with factor 2 margin 

~23% loss of wetted area for a = 3º incident angle

Protecting leading edges with monoblock shaping

 Extensive and detailed field line tracing studies have been performed that take into account variation of magnetic

a

2D 2L

h

x

b

(73)

ANSALDO W THICKNESS 6.0 mm PLANSEE W THICKNESS 7.5 mm

ITER Divertor Test Facility (IDTF) for high heat-flux testing of PFUs

 Comprehensive small-scale prototype testing

− “self-castellation” observed with some samples  does not affect qualification process

(74)

Advanced high-heat flux PFC technologies: helium

 He cooling  single phase coolant  10 MW m -2 system demonstrated [KIT]

− high temperature targets possible  increased efficiency of electricity production

− a lot of power may be required to move the helium…

KIT He COOLED “THIMBLE” (HEMJ)

(75)

Advanced high-heat flux PFC technologies: liquid lithium targets

 Can’t melt what’s already melted!  but early days…

− lots of potential issues, other metals being considered

NSTX LIQUID LITHIUM DIVERTOR (LLD)

(76)

Summary

 There are many active areas of PWI research, i.e. open questions

 The plasma boundary is complex and dynamic

 Plasma physics is an integral part of the PWI solution for high high-performance, high duty-cycle fusion devices

 A lot of recent focus has been provided by acute ITER needs, i.e. design guidance and nuclear regulation

 High-performance and robust active cooling technologies are required for the next- step in fusion energy development

 ITER PFC designs are based on conventional technologies and are well advanced

Referenzen

ÄHNLICHE DOKUMENTE

(Upper) Residual gravity anomaly over the Luusika region with outlines of gravity (solid) and magnetic (dashed) models (elliptic pipes).. Gravity profiles (measured and calculated)

M.. Proposed models, algorithms.. c ) Mapping of economic damage. d ) Calculation of fields of pollutant maximum concentrations under normal un- favorable

b) diapiric emplacement (taken from Clemens et al., 1997) 1) Tail of the diapir reaches deep into the crust; 2) Surrounding rocks and sediments are deflected upwards due to

We perform numerical simulations of the stochastic Landau-Lifshitz-Gilbert equation of motion for an atomistic spin model and calculate the magnon accumulation in linear

FIGURE 14: EXCESS DEMAND FUNCTION WITH A TARIFF.

Field and laboratory analysis of the runout characteristics of hillslope debris flows in

We have seen in the previous section that in the DLCQ we are left with all states which have positive momentum in the 11 direction. We motivated that for Matrix theory it is useful

Al- together, Lewin’s Berlin research group conceived a theoretical and empirical framework through whose prism such broad psychological categories as action (Handlung),