Survey
* Your assessment is very important for improving the workof artificial intelligence, which forms the content of this project
* Your assessment is very important for improving the workof artificial intelligence, which forms the content of this project
3rd IAEA TM and 11th IWS on ST Place: St.Petersburg Date: 04 Oct 2005 Model of filaments in plasma Nobuhiro Nishino Graduate school of Engineering Hiroshima University 1 Filaments with a wide-angle view RF antenna Filaments (and the vacuum vessel) with a wide-angle view from midplane port Center stack & gas puff RF antenna limiter 13500FPS R Usually filaments move across the magnetic field. It is clear that filaments moves across the LCFS at the initial discharge phase (Ohmic). 27000FPS Filaments hit the RF antenna and/or its limiter 13500FPS STW05 2 Multi-curtain and magnetic signal • • • • • Multi-curtain is correlated with the precursor. After ELM or blob goes through near X-point, ELM or blob is almost toroidal symmetry. Because the X-point has no poloidal field. two-dimensional multi-curtain Multi-curtain may show “ELM structure”. B.Davis R Reversed image STW05 3 Aim Filament-like edge turbulences called just “filaments” are seen in many ST and tokamak plasmas for fusion experiments. They should be related to the energy/particle confinement. Therefore, it is very important to understand the unknowns such as where the filament is forming, its radial and poloidal/toroidal extent and dynamics. What is the origin of the filament? What forces move the filament? Is it in infancy? ExB force or JxB force In experiment, a multi-curtain structure correlated to the magnetic oscillation signal. This result may show the magnetic reconnection occurs by the filament. The current may exist in the filament. Can single fluid MHD theory treat them? Single fluid MHD is very useful to treat real-sized plasmas. STW05 4 Scenario of this model (Overview) In this model, the origin of a filament is assumed to be non-homogenous (non-uniformity) heating in the same magnetic field line or flux. The origin is the hot region called “blob” due to the non-uniform heating there are the hot and cold regions in the same magnetic field line or flux. Where it is forming? Due to the temperature difference, thermal conduction parallel to the magnetic field by electrons occurs. Also thermal conduction perpendicular to the magnetic field by ions occurs simultaneously. Mainly parallel thermal conduction makes the blob into the filament-like structure. STW05 5 Scenario of this model (Overview) The current penetrate in the filament from its surroundings due to the magnetic diffusion. The change of current density makes two motions, such as expand shrink or pinch The magnetic diffusion and thermal conduction are the major role in the filament life. In this model, three important factors appears. thermal conduction parallel to the magnetic field thermal conduction perpendicular to the magnetic field magnetic diffusion STW05 6 Assumption of this model Initial plasma is the equilibrium state, then the non-uniform heating occurs. (e.g. Additional heating) The non-uniform heating process makes a “blob”. e.g. Non-uniformity of NBI heating is a few % to 10 few %. See the energy deposition profile of NB The “blob” is hotter region than that of the other region in the same magnetic field line or flux. The “blob” expands mainly along the magnetic field due to the thermal conduction. Te+DTe Te magnetic field line Te STW05 7 Aspect ratio of a filament Ratio of the thermal conductivities parallel and perpendicular to the magnetic field line are as follows. Using Lagrange system, the equation of thermal conduction is 3ne dT 2T 2 dt then 3ne L// 2 Dt// 2 // // 3.16 2 3ni a f Dt 2 neTe e me 2niTi mi i i where L// is a length of the filament along the field line, af is a length of the filament across the field line. Te=Ti, ne=ni are assumed due to the single fluid MHD STW05 8 Picture Initial condition Generation process af Additional heating thermal conduction L// magnetic field line STW05 9 Aspect ratio of a filament Let Dt //equal Dt . 2 3ne L// 2 3ni a f 2 // 2 af L // // n 2 n2 a f / L// 2 3 a f / L// 3.16e e i i B T BT 3/ 2 This value is very small at the edge parameters. - continued - 2 e.g. ne=5e18m-3, Te=20eV The ratio is 1.66e-4 Proportional to n, B-1, T-3/2 L// is not longer than 2pqR, then thermal conduction time should not excess 3ne (2p qR) 2 Dt// Dt 2 // Most likely L// is pqR. STW05 10 Evolution process of a filament Using single fluid MHD we can deduce a simple formula. dV p j B dt In steady state 0 p0 j0 B0 Then the non-uniform heating occurs. dV p p j0 j B0 B dt p nDT STW05 11 Estimation of the current density in a filament Magnetic diffusion time current penetration Using parallel circuit model the current density is estimated as at most (total plasma current may not change) 3 DTe j j0 0 1 j0 2 T e The current in the filament makes the filament motiom complicated. pinch p , B , j = j0(r) B major R Initial condition + blob or filament Hot blob in the outer region expand STW05 12 Evaluation of the jXB force of the filament Using the penetration current formula, we get j B0 j0 B 3DTe j0 B0 2Te 30 a f DTe 4Te j0 2 0 jp a f 2 0 a f j 30 a f DTe B j0 2p a f 2 4Te 90 a f DTe 2 2 jB j0 8e T Thus, j B0 is always largest. j B0 j0 B jB STW05 13 Evaluation of the jXB force of the filament -continued. Comparison with pressure gradient p DTe DT p0 e j0 B0 Te Te j B0 Thus, 3DTe j0 B0 2Te j B0 p JxB force due to the penetration current expand a part of the blob and also pinch the other part. Hot region expanding (low density region) + pinch (high density region) STW05 14 Picture Generation process Evolution process magnetic diffusion thermal conduction + pinch expand magnetic field line STW05 15 Movement depend on the initial figure of a “blob” In general a figure of blob is not spherical nor cylindrical symmetry. therefore, the movement due to the penetration current may be more complicated. pinch B new pinch region new expanding region expand The rotation may occur due to the conservation of momentum STW05 16 Gas Puff Imaging (GPI) provides various motions of filaments S.Zweben magnetic diffusion vs. thermal conduction Using the induction equation of the magnetic field B magnetic diffusion time is md 2 af 2 thermal conduction time is Dt 3ni a f 2 Therefore, if these times are equal, the generation rate of the filament may be the maximum. STW05 18 Birth location of a filament Let the magnetic diffusion time equal the thermal conduction time across the magnetic field to deduce the birth position of the filament. md Dt 0 a f 2 3n f a f 2 f 2 f n f Tf Bf 2 2 0 3 me 4 mi In H plasma, f=0.0124 In D plasma, f=0.00875 Therefore, the birth position of the filament is the edge of ST and/or tokamak plasmas. STW05 19 Evaluation of the filament velocity Vperp Using the equation of motion we get dV O( j B 0 ) dt V DTe j0 B0 DT p DT md e 0 md e md 8Te 8Te ni mi L p 8mi Lp Therefore, the filament moves rapidly with the increase temperature difference DTe. This velocity should be related to the filament generation rate, because the filament escape the heating region. STW05 20 Scaling law using this model The energy in a filament is estimated at DW f n f T f p a f 2 L// If the filament goes out the plasma within the magnetic diffusion time, then the power loss of the filament is 2 about 1 DW f n f T f p a f L// n f T f f p L// Pout ,1 fil 2 2 md 2 0 2 0 a f f where numerical factor ½ is due to the random walk of the filament. In general the energy confinement time is defined as follows. 0 Pinput W E Pinput 1 1 W E , NC E , fil STW05 21 Scaling law using this model - continued. In steady state, the power loss of the filaments may equal the total input power Pinput N Pout ,1 fil where N is related to the generation rate of the filaments. N may be related to the input power, because the generation rate of the filament depends on the temperature rise of the blob and the non-uniform efficiency. N OH POH NBI PNBI where each is the non-uniform efficiency of each heating method, respectively. In general these efficiencies are not equal. OH NBI In general OH addtional heating Then, only OH heating case N is proportional to Ptotal. In general N P , 0 1 total STW05 22 Scaling law using this model - continued. At last the energy confinement time is estimated as E , fil W N Pout ,1 fil N W n f T f f p L// 2 0 Let L// p Rq and B f Bt in above equation, and use f E , fil nT I PV Nn f 1.5 nT I PV P nf 1.5 , n f Tf Bf 2 2 0 3 me 4 mi 0 1 It seems L-mode scaling, even though this is not a dimensionless formula yet. 1 T E B F a2 B2 STW05 23 Conclusion A filament model using simple fluid MHD is proposed. In this model a “blob” mainly expand along the magnetic field line, and the blob becomes a filament-like shape. That is called the “filament”. In this model, the generation and extinction processes of the filament are decided by the magnetic diffusion and thermal conduction. According this model the scaling of the energy confinement time is also estimated. The scaling obtained using this model is similar to the Lmode scaling. STW05 24 Further problem of this model What condition determine the generation rate of filaments? Can we get completely L-mode scaling? dimensionless formula What is the L-H transition and H-mode? Can we get cheap nuclear fusion reactor? It is further question STW05 25 GPI Diagnostic setup in NSTX • Use re-entrant port and linear gas manifold. • Use He, D2, or Ar puffs. • Use beam-splitter and PMTs (100 kHz bandwidth) for discrete fast chords. Gas manifold GPI view Local magnetic field Side-viewing reentrant window S.Zweben Typical power deposition profile of NB in NSTX R.Bell STW05 27 Interpretation of multi-curtain Multi-curtain may be the magnetic surfaces, which are dragged by ELMs or blobs. • • • • • After ELM or blob goes through near X-point, ELM or blob is almost toroidal symmetry. Because the X-point has no poloidal field. two-dimensional multi-curtain Multi-curtain may show “ELM structure”. STW05 28