240 likes | 389 Views
ITPA – CDBM H-mode threshold database IAEA TM on H-mode paper October 2007 - Naka. Yves Martin Centre de Recherches en Physique des Plasmas Association Euratom - Confédération Suisse Ecole Polytechnique Fédérale de Lausanne (EPFL) CH - 1015 Lausanne, S witzerland. Preliminary remarks.
E N D
ITPA – CDBMH-mode threshold databaseIAEA TM on H-mode paperOctober 2007 - Naka Yves Martin Centre de Recherches en Physique des Plasmas Association Euratom - Confédération Suisse Ecole Polytechnique Fédérale de Lausanne (EPFL) CH - 1015 Lausanne, Switzerland
Preliminary remarks • Constraints • IAEA TM papers to be published in JPCS (J.Phys.:Conf.Ser.) • ~30 papers (half presented papers) submitted • Delay: Oct 12th • Plan • Improved but similar coverage version to JPCS • Extended version to PPCF • Comments • Received comments TT, JDB, R.Sartori, AP, FR, JS, YA
Plan • Abstract / Introduction • Threshold power estimation • DB status • Data selection • Power law scaling • Predictions for ITER • Effect of CHS • Individual density dependence • Dependence in individual devices • Prediction for ITER • Density of the minimum threshold power • Search for minimum • Discussion • Access to an H-mode with good confinement properties • PoP in H-mode DB • High H factor close to threshold • Discussion - conclusion
Abstract / Introduction • Abstract • Too long • Introduction • Too long
Threshold power estimation • DB status • Data selection • Power law scaling • Predictions for ITER
DB status • New data • Changes in data • Distribution of data
Data selection • SELEC2007= If( • :SELDBsum == 1111111111 & • :PHASE_RED == "LH" & • :AUXHEAT != "NONE" & • :AUXHEAT != "EC" & • :CONFIGs == "SN" & • :PGASA == 2 & • :KAPPA > 1.2 & • :TOK != "NSTX" & • !( :TOK == "JFT2M" & :NEL < 25000000000000000000) & • !( :TOK == "JT60U" & :NEL < 15000000000000000000) & • !( :TOK == "CMOD" & :PHASE == "D") & • !( :TOK == "CMOD" & :PHASE == "DA") & • !Is Missing( :NEL) & • !Is Missing( :PLTH_MW), • 1, 0)
Data selection II • SELDBsum = Summary of SELDBxxx where x=1,2,4v1,4v2 If(!Is Missing( :SELDB2) & :SELDB2 > 0, :SELDB2, If(!Is Missing( :SELDB1) & :SELDB1 > 0, :SELDB1, If(!Is Missing( :SELDB4v2) & :SELDB4v2 > 0, :SELDB4v2, If(!Is Missing( :SELDB4v1) & :SELDB4v1 > 0, :SELDB4v1, 0)))) • CONFIGs = Simplified configuration • "SNT", "SN", "SNB", "SN", "SN(L)", "SN", • "SNL", "SN", "SNU", "SN", "SN(U)", "SN"
Data selection III • PHASE_RED = Reduction of PHASE • D LH DA LH DH LH EDA H EDAA H EDAL HL H H HA H HCMP H HELM H HELM? H HGELM H HGELMA H HGELMHA H HL HL HOHM HL HSELM H HSELMA H HSElMA H HSELML HL HSELMOHM HL L L LA L LA1 L LA2 L LA3 L LD LH LDH LH LDHSELM LH LDL L LDSELM LH LH LH LHLHL LH LHLHLA LH LHSELM LH OHA L OHM L OHMA L OHMD LH OHMDH LH OHMDSELM LH OHMH LH OHMHSELM LH OHMP L OHMPD LH", :PHASE || " is missing")
Data selection IV • Avoid low density region?SELDB should do it • Ohmic H-mode in CMOD do not have higher threshold
Recommended scaling • Scaled parameter: PLTH • Definition changed (PFLOSS included) • With a and R: • Pthresh = 2.15 n0.782+/-0.037 B0.772+/-0.031 a0.975+/-0.08 R0.999+/-0.101 • With S: • Pthresh = 0.049 n0.717+/-0.035 B0.803+/-0.032 S0.941+/-0.019 • Discuss Kadomtsev constraint • Discuss changes from earlier versions
Prediction for ITER • Specify: Predictions done for DD plasmas • Discuss DT predictions ?
Effect of CHS • Generalize discussion to low aspect ratio toks?
Individual density dependence • Include a figure showing the non-alignement of individual devices • Density dependence on individual devices at well populated field values • Prediction for ITER
Dependence in individual devices • Include uncertainties • Remove literature points • Keep as comment on additional uncertainty? • Compare with OK work in Ryter 1996 • JS not convinced by this effect
Prediction for ITER • Prediction with lowest RMS scaling (minor radius, open divertor for JT-60U): 135MW (at high density) • Include uncertainties
Density of the minimum threshold power • Search for minimum for individual devices at well populated field values • Interpretation of global picture
Search for minimum • Include examples (figure)
Discussion • Link between different density dependences and different minimum densities • Discuss neutral penetration / collisionality
Access to an H-mode with good confinement properties • PoP in H-mode DB for HSELM and HGELM • H factor vs PoP for HSELM and HGELM • Too weak
JET DIII-D Alcator C-mod ASDEX Upgrade 0 0 0 0 1 1 1 1 2 2 2 2 3 3 3 3 4 4 4 4 5 5 5 5 6 6 6 6 PoP in H-mode DB • Why PoP > 1 in most devices ? • Why PoP < 1 in Alcator C-Mod ?
ELMTYPE Count TYP-I-II TYPE-? TYPE-I TYPE-III TYPE-RF ... HGELM 31 28 1290 2 1 1352 HSELM 0 39 164 327 13 543 PHASE 31 67 1454 329 14 1895 High H factor close to threshold • Correlation PHASE / ELM type: 2291ts were included in the dataset (missing ELMTYPE data for ~400)
R.Sartori's comment • Type I ELM H factor ~ 1 (by construction!) • Type III ELM H factor = 0.8-0.9 x type I ELM H factor • @Threshold power (or just above): type III ELMs; more power to obtain type I ELMs • PtypeI = Y x Pthresh where 1.3<Y<2.5