1

原研における TOPICS を用いた輸送研究の現状
白井浩

日本原子力研究所
1.
2.
3.
4.

那珂研究所

TOPICS 概要
原研における輸送研究
輸送モデル
統合コードと TOPICS

平成15年7月31日, 8月1日 京都大学工学部原子核工学専攻会議室

1. TOPICS 概要

2(TOkamak Prediction and Interpretation Code System)
機能 : 2 次元 MHD 平衡 + 1 次元輸送 (粒子輸送、エネルギー輸送、磁場拡散)
解析手法 : 定常解析（実験解析 ← JT-60 データ読み込み）
非定常解析（シミュレーション、実験解析 (JT-60 時系列データ)）
図形出力 : Graphic library : GDlib + Eagle（原研のオリジナル）
定常解析 : 計算終了後自動的に図形出力
非定常解析 : ステップごと・時間間隔ごとにバイナリファイル出力
別のソフト Topplt で図形出力
ファイル出力 : UFILE 形式 (1時刻、多時刻)、テキスト (CSV) 形式、バイコンパイル : Fortran 90
実行環境 : IBM SP3, Compaq(Dec) alfa WS, PC Linux, SGI Origin
移植性 : 実験解析用ファイル入力 (JT-60 実験データはバイナリで読み込み) と図形出力ルーチン以外は基本的に修正なしで可能

3TOPICS の基本式
粒子保存則
∂n i
∂
=– 1V 'Γ i + S i∂tV ' ∂ρΓi = –∇ρ

2D i ∇n i + n i v iエネルギー保存則3 ∂ n T = – 1 ∂ V 'Q + P Gain + P EQ – P Loss
jjjjV ' ∂ρ
2 ∂t j j
Qj = –∇ρn j χ j ∇T j + 3 T j Γ j22

磁場拡散∂B pη ||∂=–∂t∂ρ T R –2(j = e, i, He, imp)

j ||⋅B – j e||xt ⋅Bj ||⋅B =T = RB φTρ
V'2R–2∂∂ρ

∇V 2R

R–2Bp
ρ新古典理論
Hirshman&Sigmar
粘性係数・摩擦係数

OFMC code

MHD現象

異常輸送

Sawtooth, Tearing,Ballooning

経験則, ηi, CDBM

MI method,
NCLASS

運動量源

4TOPICS構成図

輸送係数

新古典輸送
C&H, H&H,
NCLASS

MImethod,α加熱
NBI/NBCD

主プラズマ

Stix,

熱・粒子輸送ECH/ECCD

1D transport

Bootstrap電流

1or2DFPRay tracing,
Relativistic FP

電流拡散
1D diffusion

径電場

LHH/LHCD

検討中

MHD平衡

Bonolicode*

Grad-Shafranov

新古典抵抗
H&H model

ペレット
ExB pellet model

中性粒子
SOL/Divertor プラズマ‖(Cs) & ⊥(Bohm) transport

Monte-Carlo

不純物
1D transport
IMPACT code

放射損失

* ACCOME に組み込んでいるのものを移植予定

52. 原研における輸送研究
最外殻磁気面の同定(First Boundary Identification Code)*(MHD 平衡と分布)− 定常輸送解析 −2(境界条件)

Y(ψ ) = 1 – ψ 3n + A(4) + A(2) ψ n – ψ 3n+ A(3) ψ 2n – ψ 3n

* 辻、他、JAERI-M 86-006 (1986).
** T. Takeda, et al., JAERI-M 8042(1979).

-1-22Te (keV)

2Te64200.5r/a

15Ti

8400034R (m)128ne40Ti (keV)

6ne (1019 m-3)

接線磁気プローブ
垂直磁気プローブ
ワンターンループ
円環電流フィラメント

フィッティング関数は数種類ある

1Z (m)

dP = C A(1) Y(ψ )dψ
1 F dF = C 1 – A(1) R 2 Y(ψ )pµ 0 dψψ – ψ (0)ψn =
ψ (1) – ψ (0)平衡 DBに登録

最外殻磁気面の同定(SELENE Code)**00.5r/a1

00.5r/a1

6TOPICS の定常輸送解析結果出力例
E29728, Ip=1.5 MA, Bt=3.5 T (負磁気シア)
JOULE = 0.3831
NBI-E = 1.8069
NBI-I = 6.3547
RAD = -2.5040
RFEH = 0.0000
RFIH = 0.0000
ECH = 0.0000
TOT-IN= 8.5447

イオン密度(高速成分)

イオン温度

10 1

ジュール加熱

0.0温度

4.0

電子温度

0.40.6RHO

0.81.0-1.00.20.40.6
RHO

ブートスト
ラップ電流

0.00.20.81.0Q
J//
BSC
R FC D
NBCD
ECCD

安全係数

2.0
0.0

0.0

0.2

11.0

Ele.
Ion

6.0

8.0

0.8

Ele.
Ion

放射損失

NBI 電子加熱*10
1.00.40.6
RH O

NBIE
NBII
JOULRAD
ECH
RFE
R FI

NBI イオン加熱

Q , J/JMAX
0.0

0.00.210.00.0

DD = 0.000D+00
TH = 0.000D+00
TB = 0.000D+00
BB = 0.000D+00
DT = 0.000D+00
TH = 0.000D+00
TB = 0.000D+00
T-FRAC= 5.000D-01

プラズマ電流0.40.6RHO電流0.8

イオン
熱拡散係数

輸送

電子
熱拡散係数

10 -22.0

イオン密度(熱化成分)

DATE=03-07-18

TAUE = 0.4538
TAUP = 0.3000
ZEF-SP= 0.0000
ZEF-NC= 0.0000
CKI-NC= 0.0000
TI0-NC= 0.0000
NC-FAC= 1.0000

パワー

POW (MW/m**3)

3.04.0

密度6.680sec

KAI (m**2/s)
10 -1
10 0Ele.
Ion
fast

電子密度1.0

DEN (10**20/m**3)
TEM (keV)*10 -1
0.0 1.0 2.0 3.0 4.0 5.0 6.0 7.06.680sec / EQDB_RV :-1*106.
05.0E029728 case012 tokrdrv_ofmc :DISTRIBUTION OF PLASMA
BETA-T= 0.6578
AV-DE = 0.2332
BETA-J= 0.9123
AV-TE = 4.6076
TCUR = 1.4971
AV-TI = 6.2651
TBSC = 0.9247
WS-E = 1.6953
QAXI = 5.4176
WS-I = 1.4761
QSUR = 6.4695
WS-B = 0.7062
VL = 0.0000
TOT-WS= 3.8776
VL-NC = -0.2045
ZEFF = 2.8000

1.00.00.2

0.40.6
RHO0.81.0TOPICS 非定常輸送解析（高イオン温度 H モード） 7ne (1019 m-3)

Ti (keV)Te (keV)
402063047.0 107.0
6.56.526.0e
ti m

5.50.50.0

1.005.50.00.5r (m)1004(1) (2)

(3)

100
30051L-mode(3)0time (s)70.5(2)0.1
61.0s)

10(1)55.5
0.0e(mti

実験データをもとに多時刻分布フィッティングデータを用意して、フォッカープランク方程式から NBI パワーの時間変化を求め、輸送係数、中性子発生率などを評価する。

Bt = 4.0 T

080

6.0

r (m)

χe (m2s-1)(MW)(arb.)(MJ) (1019 m-3)

neWdia
PNBI

1.0s)

r (m)

Ip = 2.1 MADαdiv(eti m

6.07.06.5

χi (m2s-1)0(s)2080.2 0.4 0.6 0.8r/a1(1)
1 Hot Ion
L-mode(2)(3)
Hot IonH-mode0.1
00.2 0.4 0.6 0.8r/a1

8

シミュレーション例 (新古典テアリングモード)

林 伸彦 (解析室)JT-60SC のパラメータを想定

Ip=3.0 MA, BT=3.8 T, R=2.9 m, a=0.85 m, PNBI=30 MW, ne=4.2×1019 m-3,
Te=6.4 keV, Ti= 7.5 keV・拡散係数は経験則
Di = χ e =χi= 0.06(1 + 2 ρ 2 )(1 + PNBI 0.5 ) m2 / s2・NBI 加熱分布と放射損失分布を固定
・磁気島内の拡散係数を増加して密度・温度の平坦化
・ブートストラップ電流で電流分布が変化
・NTM の最初のトリガーは∆'>0NBI30 MW

02002468 10 12 14

時刻 (s)

140t=3s

2
1.57s0ρs
全電流

1∆'(0)∆'βNj (MA/m2)

2.5∆'(W)ρsブートストラップ磁気島位置

0.5 電流0-8000ρρ

磁気島幅 W

0.2

10

4812
時刻 (s)

9

シミュレーション例（新古典テアリングモード続き）
8s (NBI)

6 12s(NBI)
4 14s (OH)
200t=2s

161212(OH)ρ2.584001

Ti (keV)

816Te (keV)ρ10.60ρρ1ρ143210145jBS (MA/m2)

j (MA/m2)
08006qne (1019m-3)

100

ρ100

径電場の評価 (径電場シアと輸送係数の軽減)
運動量バランスの式

m an adua= e a n a E + u a×B – ∇Pa – ∇⋅Π a + F a 1 + M adt
摩擦テンソル

熱運動量バランスの式

10

摩擦力

運動量源

qaema ∂= a q a×B – 5 n a∇T a – ∇⋅Θ a + F a2 + Q a∂t T a2Ta
粘性
熱摩擦力テンソル

熱運動量源

定常状態 (左辺 = 0) を仮定して磁力線に平行な成分を磁気面平均B ⋅ F a 1 = B ⋅ ∇⋅Π a – B ⋅ M a – e a n a B ⋅ E
B ⋅ F a2 = B ⋅ ∇⋅Θ a – B ⋅ Q a
7×7 の行列式にまとめる
L U || = M U || + V – S – εすなわち U || = L – M–1M V – L–M–1L 規格化摩擦行列
M 規格化粘性行列
U || 磁力線方向の流速に関するベクトル
V 熱力学的力
S 運動量・熱運動量源
ε 磁力線方向の電場に関するベクトル

S– L–M–1ε

M. Kikuchi and M. Azumi
PPCF, 37, 1215 (1995)

11

径電場分布の評価 (続き)･

磁気面平均された磁力線に平行な不純物速度 (不純物の実験データを使うため)
7Bu ||I =

ΣL–M

j=1LM
Vj
Sj

–1MI,jVj – L – M

ブート
ストラップ項–1B•M j – e j n j L – MI,j–1I,j

誘導項

外部駆動項

B•E j

(無視)

摩擦行列
･ 運動量
･ 熱運動量j=1:
電子
j=5 : 電子
粘性行列
=2j 熱化イオン:
j=6 : 熱化イオン
熱力学的力
j=3 : 不純物
j=7 不純物:=4j 高速イオン:
運動量・熱運動量源 (OFMC code)･ 不純物のトロイダル速度Bφ

荷電交換u Iφ =
2
再結合分光
B

Er = – dΦ = 4 1dρΣ α* I,jj=1u*I φ = u I φ +BφB

Bu ||I – 1 –B2

B 2φB2d PI
R d Φ + e 1ndψI I dψ

(ポロイダル回転は新古典理論に従う)

43dP jdT jdψ11*–Ru I φ –α* I,j e n d ρ + α* I,j+4 e d ρ2 dρj jjj=1j =1RB φΣ

7–1B•ML –MΣ
2 j =1I ,jj

α i ,j = M – L

Σ–1Mi ,jα* I ,j = α I ,j –B 2 / B 2φ – 1 δI ,j

12

定常輸送解析（内部輸送障壁と熱拡散係数, 径電場）
イオン温度 (keV)

内部輸送障壁

10径電場(kV/m)

00-25-50
10径電場シアが小さい

イオン
熱拡散係数
新古典拡散係数(MI method)

0.1
00.5r/a1510
内部輸送障壁500-25-5010(輸送障壁なし)

1(箱型内部輸送障壁)

熱拡散係数(m2/s)

20

熱拡散係数
(m2/s)

径電場(kV/m)

イオン温度 (keV)(パラボラ型内部輸送障壁)

径電場
シアが大きい(輸送障壁なし)

1

新古典拡散係数(MI method)

0.1
100.5r/a

イオン熱拡散係数1径電場・フローの分布の時間発展シミュレーション (予定)
運動量バランスの式・熱運動量バランスの式を解く

13→ 実験におけるフロー分布 (不純物のトロイダル回転) の時間変化との比較
ITB の時間発展との consistency (Er 分布)、輸送モデリングへの寄与(1) 時間微分項を残す
運動量・熱運動量バランスの式の左辺の m an a d u a⋅Bdt

1
および m a ∂∂t T a q a⋅B

10(CO 入射82u)

6
不純物

4電子
水素

2000.5r/a1

運動量源(10 -2 N/m3)

運動量源(10 -2 N/m3)(2) 運動量源の評価
0電子-1

水素-2-3-4

現在 M a は OFMC の定常計算で評価

不純物(CTR 入射 2u)00.5r/aフォッカープランク方程式で評価

1(3) 異常粘性項 (運動量に関する異常拡散係数)を考慮
22 1ano右辺に m an a µ ano
∇u
⋅
Bmµ∇
q ⋅B
および
aaa qaTa aanoµ ano~
χを仮定
aaを追加

3. 輸送モデル14

弱い磁気シアで輸送係数が減少する輸送モデルx e = x ae no + C e3x Ni CE (s) =x i = x ai no + C i3x Ni C11+α[s-∆∆s])α[s+∆∆s])
1 + exp(-α1 + exp(α12x = C e1(1 + C e2ρ ) HHa(y,2) Ee(s)τE2a2anox i = C i1(1 + C i2ρ ) HH(y,2) Ei(s)τE2–1Eβ (s) = 1 + exp α β (s + ∆ β )
–1

∆s
= 0.5

2.50.8E (s)

anoe

1.0

0.6

2.00.41

1.50.2+ 1 + exp –α β (s – ∆ β )α=100-4

D = Cd xe-3-2-101234shear sαe/i: = 5 (熱拡散係数の減少の度合い)∆e/i: = 1 (ITB の脚の位置)

Cd: = 0.1 (ne の ITB を再現するため)10.8E (s)χiNC : 新古典イオン熱拡散係数
Ce/i 1 : τEth を τEHH(y,2) に合わせる
Ce/i 2 : = 4 (周辺での熱拡散係数の増加)
Ce/i 3 : = 1

α=10.60.420.243∆s=15・平均的なχの絶対値は・s → 0 で異常輸送

ITER 閉じ込め則から与える→ 00-4-3-2-1

K. Ushigusa, et al. Proc. 17th IAEA 1998, IAEA-CN-69/FTP12 (1999).

0shear s

123415輸送モデルと JT-60U の輸送係数との比較
E32611
Ip
[MA]

100

αDα[A.U.]

10Wdia[MJ]χmodelχeχi215050PNB[MW]E32611(7.7s)

ダブルバリアプラズマ

準定常状態
050135[1019m-3]
60.1

00.20.4

r/a0.60.8

1.0

C1=0.5, C2=10, α =1.2, △ =2.5
比較的強い ITB を持つプラズマで係数の値を設定

7Time[s]91113[keV]
At t=7.7s
10
Ip=1.5MA
8
Bt=3.5T4
R=3.6mne6Ti
a=0.84m
4q95=3.8
Te2q
PNB=10.5MW2
τE=0.33s
00QDTeq=0.12
0 0.2 0.4 0.6 0.8 1.0 0 0.2 0.4 0.6 0.8 1.0
ρρ16ITB 輸送モデルの改良・牛草モデルは qmin に ITB の脚があり、s の値のみで χ の低減を決めている→ 箱型 (強い) ITB や電子系の ITB の場合は良さそうだが、弱い ITB の場合は一考が必要
・牛草モデルの表式ではχの値を調整する係数が多く、汎用的に使えない・これらの係数を JT-60 実験結果を基にプラズマパラメータで置き換えて、ITB の領域から輸送モデルの改良を試みる(base )
χj = χ jω E xB T e
ρ*1–fγ L , T i , s, α, ρ * ,...

ITB

輸送障壁 (ITB, ETB) 以外の領域の輸送としては、従来より用いている ITER 閉込め則ベースの経験的熱拡散係数や、Weiland ηi [1], MMM (Multi-mode mod[2] を当面使用し、これに ITB の実験結果を相乗する[1] H. Nordman, et al., Nucl. Fusion, 37, 413 (1997).
[2] J. Kinsey and G. Bateman, Phys. Plasmas, 3, 3344 (1996).

ITB 輸送モデル改良について

17

ITB の脚の位置：∆e , ∆i 、 ITB の広がりと強さ：αe , αi → プラズマパラメータ化
・ITB 形成に関して電子系 (s が効く) とイオン系 (シアフローが効く) を分離
・ITB の脚位置は∇Ti の値如何で ρmin より内側にも外側にもなる(坂本氏 '00 IAEA)→ 正の s のある値の所で脚の位置が決まる (∇Ti が小さい場合) → ∆i をシフト
qinteger と ITB 形成の関連性 (JET) は?
・ITB の幅 ∝1/Bp (強い ITB の場合)、弱い ITB ほど幅が広くなる？→ αe , αi の値
・ITB の強さ (∇P の大きさ) (← 磁気シア、外部運動量、ECH) → αe , αi の値
・ITB の加熱パワー閾値をどう考える？
(local な熱流束そのもの?, ∇P?, α パラメータ?)
・イオン系は ωExB が陽に入った表式が必要 ?（但し ∝ {1-ωExB/γL} で良いのか）・Ti ITB と Te ITB は互いに影響を及ぼし合うことを考慮→ χe と χi の表式の一部をクロスさせる
更に・ITB における熱ピンチ・粒子ピンチの有無
・ETB と ITB の関連性 (ダブルバリアプラズマの形成・維持条件)

更に

18

ベースになる輸送モデルも改良
・JT-60 過渡応答特性実験 (今年の冬) の解析結果を輸送モデルに反映
ガスパフモジュレーション
ECH モジュレーション
運動量入力モジュレーション
理論シミュレーションからのフィードバック (NEXT)
ジャイロ粒子・流体シミュレーション・揺らぎの自己組織化 (乱流 ⇄ 帯状流)・異なるスケールの揺らぎ間の複合的な相互作用↓(1) 電子系・イオン系輸送の関連性(2) 輸送障壁形成・維持 (揺らぎの自己抑制) の条件、輸送の制御条件のモデル化
輸送コードで使用可能な定式化が必要

19

4. 統合コードと TOPICS
原研の解析コード群(整備中のものも含む)

統合コード

JT-60実験

平衡・輸送をベースにした
磁気核融合プラズマ時間発展解析コード・実験データとの比較による最適化
・核燃焼プラズマ予測
・運転シナリオの最適化

燃焼模擬実験、長時間帰還制御等

解析コード群の開発・統合化

高エネルギー粒子
挙動解析コード(OFMC、他)

周辺プラズマモデルコードコアプラズマ
輸送解析コード（TOPICS）(TOPICS)

原研 NEXT 研究会 (平成 15
核燃焼プラズマコード構想

年 3

MHD 解析コード(ERATO、TSC 他)

月 4 日)

ITER,

JT-60SC,

コア・ダイ
バータ結合
コード

JT-60U

ダイバータプラズ
マ輸送解析コード(SONIC)

プラズマを対象

TOPICS と TASK のベンチマーク
ベンチマーク
テス
トTASK

20

JT-60

(定常解析・時間発展

UFILE(テキスト形式データ)

解析 実験データベース)

I T PATOPICS

ITPA トピカル物理グループ
・輸送及び内部障壁の物理・閉じ込めデータベースとモデリング(・定常運転と高エネルギー粒子)

JT-60 実験テーマ班「ITB のモデリングサブジェクト」
・輸送モデルの比較・検証
・輸送モデルの構築 (精度の向上)
・輸送機構の解明

UFILE の例 (ITER 物理 R&D 分布 DB 用、テキスト形式)21j 34487 JT60 2 0 6;-SHOT #- F(X) DATA -uwrite- 18-APL-00(E34487, t=4.6s, Ti)
19-DEC-02;-SHOT DATE- UFILES ASCII FILE SYSTEM
0;-NUMBER OF ASSOCIATED SCALAR QUANTITIESRHO;-INDEPENDENT VARIABLE LABEL: XTIMESECONDS ;-INDEPENDENT VARIABLE LABEL: YTIeV;-DEPENDENT VARIABLE LABEL4;-PROC CODE- 0:RAW 1:AVG 2:SM 3:AVG+SM
51;-# OF PTS- X, F(X) DATA FOLLOW:
1;-# OF Y PTS- X,Y,F(X,Y) DATA FOLLOW:
0.000000E+00 2.000000E-02 4.000000E-02 6.000000E-02 8.000000E-02 1.000000E-01
1.200000E-01 1.400000E-01 1.600000E-01 1.800000E-01 2.000000E-01 2.200000E-01
2.400000E-01 2.600000E-01 2.800000E-01 3.000000E-01 3.200000E-01 3.400000E-01
3.600000E-01 3.800000E-01 4.000000E-01 4.200000E-01 4.400000E-01 4.600000E-01
4.800000E-01 5.000000E-01 5.200000E-01 5.400000E-01 5.600000E-01 5.800000E-01
6.000000E-01 6.200000E-01 6.400000E-01 6.600000E-01 6.800000E-01 7.000000E-01
7.200000E-01 7.400000E-01 7.600000E-01 7.800000E-01 8.000000E-01 8.200000E-01
8.400000E-01 8.600000E-01 8.800000E-01 9.000000E-01 9.200000E-01 9.400000E-01
9.600000E-01 9.800000E-01 1.000000E+00
4.600000E+00
1.874068E+04 1.870750E+04 1.860800E+04 1.844235E+04 1.821087E+04 1.791393E+041.755144E+04 1.712592E+04 1.664338E+04 1.610725E+04 1.552086E+04 1.488848E+04
1.421486E+04 1.350508E+04 1.276465E+04 1.200010E+04 1.121801E+04 1.042509E+049.629133E+03 8.838983E+03 8.063406E+03 7.311970E+03 6.595521E+03 5.925367E+035.314926E+03 4.778537E+03 4.333539E+03 4.001851E+03 3.821413E+03 3.685258E+033.523527E+03 3.366604E+03 3.215970E+03 3.071439E+03 2.933008E+03 2.800604E+032.674175E+03 2.553689E+03 2.439055E+03 2.330246E+03 2.227179E+03 2.129805E+032.038052E+03 1.951848E+03 1.871130E+03 1.795830E+03 1.725874E+03 1.661193E+031.601667E+03 1.547420E+03 1.500000E+03;----END-OF-DATA-----------------COMMENTS:-----------TOPICS からの UFILE 出力変数一覧
小半径定義
小半径の対応 (RHOPHI)

平衡関係
プラズマ体積 (VOL)
メッシュゾーンごとの体積 (DVOL)
プラズマ断面積 (AREA)
プラズマ表面積 (SUF)
大半径 (RMAJ)
小半径 (RMIN)
楕円度 (ELL)
三角度 (TRG)
メトリック <¦∇ρ¦> (GR1)
メトリック <¦∇ρ¦^2> (GR2)

実験データ
電子密度 (NE)
高速イオン密度 (BDENS)
電子温度 (TE)
イオン温度 (TI)
有効電荷数 (ZE)
トロイダル回転速度 (VTCX)
トロイダル回転周波数 (WT)

圧力電流電場
プラズマ圧力（バルク）(PTH)
プラズマ圧力（高速イオン）(PBEAM)
プラズマ電流 (CUR)
ブートストラップ電流 (CURBS)
安全係数（計算値） (Q)
安全係数（MHD 平衡データ） (QQ)
周回電圧 (V)
新古典プラズマ抵抗 (ETA̲NC)
径電場 (ER)ExB シア率 (EBSHEAR)

source項

NBI による粒子源 (SBTH)
リサイクリングによる粒子源 (SWTOT)
ジュール加熱パワー (POH)
NBI 加熱パワー（電子）(PBE)
NBI 加熱パワー（全イオン）(PBI)
放射損失パワー (PRAD)
運動量源（電子）(MELE)
運動量源（水素イオン）(MDEU)
運動量源（不純物イオン）(MIMP)
運動量源（全イオン）(MION)

22

輸送関係
電子熱拡散係数 (CONDE)
イオン熱拡散係数 (CONDI)
伝導熱流束（電子）(QCNDE)
伝導熱流束（イオン）(QCNDI)
対流熱流束（電子）(QCNVE)
対流熱流束（イオン）(QCNVI)
エネルギー等分配熱流束 (QIE)
放射損失熱流束 (QRAD)
荷電交換損失熱流束 (QCX)

23

進行状況

(ITB モデリングサブジェクト・ITPA 閉じ込め DB とモデリング・・の輸送研究)• 輸送モデルとしてCDBM を導入、JT-60 データでの検証• TOPICS/TASK の非定常輸送シミュレーションのベンチマークテスト• プラズマ流／回転の効果を取り入れた輸送モデルの構築• ITB 形成に関する過去のデータの解析• ITB 準定常維持条件に関する過去のデータの解析• JT-60 の過去の過渡応答実験のデータの見直し• JT-60 過渡応答実験 (12 月)
NBI 運動量入力変調、ECH 加熱パワー変調、ECH 局所加熱、ペレット・ガス注入密度摂動→ 解析が終わったら ITPA 分布データベースへの登録

更に•CDBM 以外の ITPAの輸送モデルの検証•

輸送理論 (粒子シミュレーションなどの結果)のモデル化・定式化

24

今後の課題
・輸送モデルの構築(輸送係数だけでなく、輸送障壁の形成・維持条件、フロー・径電場の時間発展・・)→ 核燃焼プラズマ予測・最適運転シナリオの研究
・MHD 不安定性の効果のプラズマ性能 (輸送)への反映

(高β核燃焼プラズマ)
アルファ粒子 → アルヴェン固有モードなども・TOPICS/TASK のほかのベンチマーク、解析結果の共有