2020-03-27input.dat 7A modules

作者: 锅炉工的自我修养 | 来源:发表于2020-03-27 22:15 被阅读0次

总结

1. 靶板eirene下标，默认值可能为25000
1. 对于C的追迹，发现测试粒子非常少
1. 通过手动添加，不成功？（化学溅射为0？）
1. 明白7A模块的结构
1. 即使深度脱靶，ITER的辐射值依然很低。（1.43MW，电子温度0.55eV）
1. sfd 辐射显著高于SD
1. 靶板夹角，导致粒子从内靶板向外办反射，同时，水平靶板导致外靶板在碳聚集在靶板附近

辐射统计

echo "write phys brna 0.0 rmma pot m* 0 0 sumz ev rm* -1.0 rm* brhe m- b2ra 0 0 sumz m+ b2br 0 0 sumz m+ 0 nreg m* sumx sumy 19 f.x surf" | b2plot ;mv b2pl.exe.dir/b2plot.write 2D_data/rad_test_w.last10

之前的公式

echo "comp b2ra b2br m+ abs 0 0 sumz neutrad abs 0 0 sumz
m+ sumx sumy jxa writ f.y" | b2plot

存在的问题

计算公式计算的结果，超过总的 $P_{SOL}$

trajectory

存在的问题

对于没有喷气，如何定位不同区域的trajectory
之前的source，不显示NPTS
定位显示图片的位置

image.png

添加靶板区域的source

1 .SOLPS5.0

对应的b2.neutral.parameters

&NEUTRALS
 NSTRAI=9,
 CRCSTRA= 'W','E','S','S','N','C','C','V','T',
 RCPOS=    -1, 96, -1, -1, 36, 0,   0, 0,  0,
 RCSTART=   0, 0,   0, 72,  0, 0,   0, 0,  0,
 RCEND=    35, 35, 23, 95, 95, 0,   0, 0,  0,

对应位置：IT,OT,PFR_inner,PFR_outer,North,Gas Puffing 1, Gas Puffing 2,Volume recombination

image.png

7A modules

     8
     0     0     0     0     0     0     0     0
 5.00000E-01
*   1 : Recycling - Target       1                                      
FFFF
 25000    -1     7     1
 1.00000E+00 0.00000E+00     0     0     0     0     0     0     0     0
FFFT
    -1
FFTFF
     1
     1     2     1     1    37     1     1     0     0
 1.00000E+00 1.04000E+02 1.00000E+00 0.00000E+00 2.00000E+03
    -1     0     0     0     0     0     0
 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
 1.00000E+00 9.00000E+01 0.00000E+00 0.00000E+00 0.00000E+00
*   2 : Recycling - Target       2                                      
FFFF
 25000    -1     7     1
 1.00000E+00 0.00000E+00     0     0     0     0     0     0     0     0
FFFT
    -1
FFTFF
     1
     2     2    99     1    37    99    99     0     0
 1.00000E+00 1.04000E+02 2.00000E+00 0.00000E+00 1.00000E+03
    -1     0     0     0     0     0     0
 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
 1.00000E+00 9.00000E+01 0.00000E+00 0.00000E+00 0.00000E+00
*   3 : Recycling - PFR side     1                                      
FFFF
  2000 10001     6     1
 1.00000E+00 0.00000E+00     0     0     0     0     0     0     0     0
FFFT
    -1
FFTFF
     1
     3     1     1     1     1     1    25     0     0
 1.00000E+00 1.40000E+02 3.00000E+00 0.00000E+00 2.11100E+03
    -1     0     0     0     0     0     0
 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
 1.00000E+00 9.00000E+01 0.00000E+00 0.00000E+00 0.00000E+00
*   4 : Recycling - PFR side     2                                      
FFFF
  2000 10001     6     1
 1.00000E+00 0.00000E+00     0     0     0     0     0     0     0     0
FFFT
    -1
FFTFF
     1
     4     1     1     1     1    75    99     0     0
 1.00000E+00 1.40000E+02 4.00000E+00 0.00000E+00 2.11100E+03
    -1     0     0     0     0     0     0
 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
 1.00000E+00 9.00000E+01 0.00000E+00 0.00000E+00 0.00000E+00
*   5 : Recycling - Outer side   1                                      
FFFF
  2000 10001     6     1
 1.00000E+00 0.00000E+00     0     0     0     0     0     0     0     0
FFFT
    -1
FFTFF
     1
     5     1    37    37    37     1    99     0     0
 1.00000E+00 1.40000E+02 5.00000E+00 0.00000E+00 1.11100E+03
    -1     0     0     0     0     0     0
 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
 1.00000E+00 9.00000E+01 0.00000E+00 0.00000E+00 0.00000E+00
*   6 :  Gas puffing source           1 : D2                            
FFFF
 10000    -1     1     1
 1.00000E+00 0.00000E+00     0     0     0     0     0     0     0     0
FTFF
     1
FFTFF
     1
     6     0    63     1     1     1     1     0     0
 1.00000E+00 1.20000E+02 0.00000E+00 0.00000E+00 1.00000E+03
     0     0     0     0     1     0     0
 0.00000E+00 0.00000E+00 7.96750E-01 2.33554E+01 0.00000E+00 0.00000E+00
 3.00000E-02 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
 1.00000E+00 9.00000E+01 0.00000E+00 0.00000E+00 0.00000E+00
*   7 :  Gas puffing source           2 : D2                            
FFFF
 10000  6001     1     1
 1.00000E+00 0.00000E+00     0     0     0     0     0     0     0     0
FTFF
     1
FFTFF
     1
     7     0   177     1     1     1     1     0     0
 1.00000E+00 1.02000E+02 0.00000E+00 0.00000E+00 1.00000E+03
     0     0     0     0     1     0     0
 1.52102E+02 1.57595E+02 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
 3.00000E-02 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
 1.00000E+00 9.00000E+01 0.00000E+00 0.00000E+00 0.00000E+00
*   8 : Volumetric recombination D                                      
FFFF
  7000  4001     6     1
 1.00000E+00 0.00000E+00     0     0     0     0     0     0     0     0
FFFT
     1
FFFTF
     1
     8     0     0     1    37     1    99     0     0
 1.00000E+00 1.40000E+02 1.00000E+00 0.00000E+00 0.00000E+00
     0     0     0     0     0     0     0
 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00

2. SOLPS-ITER

b2.neutrals.parameters

&NEUTRALS
 nstrai= 12, 
 crcstra= 'W', 'W', 'E', 'E', 'S', 'S', 'S', 'S', 'N', 'N', 'V', 'V',
 rcpos= -1, -1, 96, 96, -1, -1, -1, -1, 36, 36,  0,  0,  
 rcstart=  0,  0,  0,  0,  0,  0, 72, 72,  0,  0,  0,  0,  
 rcend= 35, 35, 35, 35, 23, 23, 95, 95, 95, 95,  0,  0,  
 species_start=  0,  2,  0,  2,  0,  2,  0,  2,  0,  2,  0,  2,  
 species_end=  1,  8,  1,  8,  1,  8,  1,  8,  1,  8,  1,  8,

IT,OT,PFR_inner,PFR_outer,North,Volume recombination for both D and C

3. 7A模块的书写逻辑

1. 指定初始源项的概率分布 $Q(i,$ r,t,v $)$
1. 初始分布状态：LOCATE(isp,r,t,v)
1. 五种类型的暂态初始源分布
- 点源
- 线源（Gas puffing）
- 表面源（target）
- 体复和源
- 统计源（“census source”）

前四种为非时间依赖项，第五种为时间依赖(给定时间点)

1. 因为输运方程式线性的，可以按顺序处理各种类型的测试源，然后线性叠加。
1. 源的强度通过指定FLUX(ISTRA)确定。

1. 可以将一层分成几个子层(sub-strata——NSRFSI)

*** 7. Data for primary sources, nstrai strata                          
     8
     0     0     0     0     0     0     0     0
 5.00000E-01

7. Data for primary sources
层数(ITER 为(WENSSV[C])粒子种类 )
不同层的文件读取类型
- 0-5:可用户例程修改（添加的源项，应该用0？）
- 6: 不从该处读取数据，（该层必须完整的被例程定义），在EIRENE+B2.5中耦合运行
- -1: 不修改该层（不调用 IF2COP(ISTRA)）

ALLOC: 分配不同层的时间占比 $(1-ALLOC)*NPTS+ALLOC*FLUX$
AMPTS：同时指定总CPU时间(NTCPU)和MC轨迹数（NPTS(ISTRA)）

*  11 : Volumetric recombination D     
FFFFF
  7000  4001     6     1     0    
 1.00000E+00 0.00000E+00     0     0     0     0     0     0     0     0    
FFFTF
     1    
FFFTF
     1    
    11     0     0     0     0     0     0     0     0    
 1.00000E+00 1.40000E+02 1.00000E+00 0.00000E+00 0.00000E+00
     0     0     0     0     0     0     0    
 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00

TXTSOU:提示文本
未使用，未使用，
NLSYNP：y轴的对齐文件
NLSYMT 环向对齐，ditto
缺省

NPTS：if NPTS>0; 最大测试粒子数；elseif 关闭该层；else 不限制；
NINITL：随机数种子：if NINITL>0;可重现；elseif 不初始化；else 时钟随机数种子
NEMODS:选择给定位置，粒子和能量粒子的能量分布（速度空间）
NAMODS：ditto，角度分布
- 1 沿着向量C
- 2 沿着向量-C
NMINPTS:每层测试粒子的最小值,默认为零

FLUX，SCALV: SCALV=0,源项单位为A
- FLUX 为原子流:FLUX/1.6E-19/2,for $H_{2}$ per second
- $SCALV\ne 0$ 设置中性离子密度描述一个特定网格
IVLSF: 1 volume avaerged ;2 surface averaged
ISCLS: 粒子下标选择标签
ISCLT: 记录数选择标签
ISCL1: if IVLSF=1, nx,ny,nt,nblock,nacell
ISCL2:
ISCL3:
ISCLB:
ISCLA:

粒子下标分布

NLATM: 原子源, ISP=IATM,初始权重 NPRTA(IATM),ITYP=1
NLMOL: 分子 ditto 2
NLION: 测试离子 ditto 3
NLPHOT:光子源，ditto 0
NLPLS: 主等离子体源 ditto 5

五个中只能有一个为T

NSPEZ:粒子源种类下标
- $\geq 1 \& \leq NATM, e.g.$ 指定粒子源，不随机取样

物理空间分布：

NLPNT 点
NLLNE 线
NLSRF surface
NLVOL 体
NLCNS 抽样源

NSRFSI：不同点的总数，该层的总数

image.png

INUM: 子层下标 11
INDIM: 指定源所在的表面类型
- 0 添加的表面：structure elements
- 1 标准surface x
- 2 ditto y
- 3 ditto z
- 4 多个elements组成的多边形
INSOR: 网格surface的标号 nx,ny,nz（surface），asurface
INGRDA: 定义源项在子区域的标准surface的变化
INGRDE:ditto

最后两个 multiple 3

image.png

SORWGT: 起始点的相对频率，归一化
SORLIM：=KLMN,抽样初始化，else :
- K:起始时间与分布下标
- L: Zsurface ditto
- M: y- ditto
- N: x ditto

M,N,L=0,相对坐标是从其他两个计算得到的
- [x] 有且只有一个为零
- [x] 如果INDIM=0,添加的表面，任何一个数都可以为0
- [x] 如果INDIM=1, x surface,N=0
- [x] 2 M=0 ditto
- [x] 3 L=0 ditto

L,M,N=1, $\delta$ 分布，（a+b）/2
L,M,N=2,均匀分布
L,M,N=3,指数衰减
L,M,N=4,分段函数，只能有一个为4
K=1, $\delta$ 分布
K=2,均匀分布

SORIND：if SORIND>100;SORIND=LMN,选择step函数
SOREXP: opt=3的衰减长度
SORIFL：前4个数可以更改表面初始的方向

image.png

NRSOR:点源
NPSOR：ditto
NTSOR：
NBSOR:
NASOR：
NISOR：

image.png

SORAD1~6
- [left~end]point for x,y,z [a,b]区间

distribution in velocity space

SORENI ： $T_e,T_I,V_{i}$ ,测试粒子的初始状态
SORENE
SORVDX
SORVDY
SORVDZ

SORCOS:
SORMAX: 角度分布的，依赖NAMODS值
SORCTX:非零，按照向量初始化
SORCTY:
SORCTZ:

不改变标签

image.png

*** 7. Data for primary sources, nstrai strata                          
    12
     0     0     0     0     0     0     0     0     0     0    -1    -1
 5.00000E-01
 *   1 : Recycling - Target       1 D                                      
FFFF
 25000    -1     7     1
 1.00000E+00 0.00000E+00     0     0     0     0     0     0     0     0
FFFT
    -1
FFTFF
     1
     1     2     1     1    37     1     1     0     0
 1.00000E+00 1.04000E+02 1.00000E+00 0.00000E+00 2.00000E+03
    -1     0     0     0     0     0     0
 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
 1.00000E+00 9.00000E+01 0.00000E+00 0.00000E+00 0.00000E+00
  *   2 : Recycling - Target       1 C                                      
FFFF
 25000    -1     7     1
 1.00000E+00 0.00000E+00     0     0     0     0     0     0     0     0
FFFT
    -1
FFTFF
     1
     2     2     1     1    37     1     1     0     0
 1.00000E+00 1.04000E+02 1.00000E+00 0.00000E+00 2.00000E+03
    -1     0     0     0     0     0     0
 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
 1.00000E+00 9.00000E+01 0.00000E+00 0.00000E+00 0.00000E+00
 *   3 : Recycling - Target       2 D                                      
FFFF
 25000    -1     7     1
 1.00000E+00 0.00000E+00     0     0     0     0     0     0     0     0
FFFT
    -1
FFTFF
     1
     3     2    99     1    37    99    99     0     0
 1.00000E+00 1.04000E+02 2.00000E+00 0.00000E+00 1.00000E+03
    -1     0     0     0     0     0     0
 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
 1.00000E+00 9.00000E+01 0.00000E+00 0.00000E+00 0.00000E+00
 *   4 : Recycling - Target       2 C                                      
FFFF
 25000    -1     7     1
 1.00000E+00 0.00000E+00     0     0     0     0     0     0     0     0
FFFT
    -1
FFTFF
     1
     4     2    99     1    37    99    99     0     0
 1.00000E+00 1.04000E+02 2.00000E+00 0.00000E+00 1.00000E+03
    -1     0     0     0     0     0     0
 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
 1.00000E+00 9.00000E+01 0.00000E+00 0.00000E+00 0.00000E+00
  *   5 : Recycling - PFR side     1 D                                      
FFFF
  2000 10001     6     1
 1.00000E+00 0.00000E+00     0     0     0     0     0     0     0     0
FFFT
    -1
FFTFF
     1
     5     1     1     1     1     1    25     0     0
 1.00000E+00 1.40000E+02 3.00000E+00 0.00000E+00 2.11100E+03
    -1     0     0     0     0     0     0
 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
 1.00000E+00 9.00000E+01 0.00000E+00 0.00000E+00 0.00000E+00
 *   6 : Recycling - PFR side     1 C                                      
FFFF
  2000 10001     6     1
 1.00000E+00 0.00000E+00     0     0     0     0     0     0     0     0
FFFT
    -1
FFTFF
     1
     6     1     1     1     1     1    25     0     0
 1.00000E+00 1.40000E+02 3.00000E+00 0.00000E+00 2.11100E+03
    -1     0     0     0     0     0     0
 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
 1.00000E+00 9.00000E+01 0.00000E+00 0.00000E+00 0.00000E+00
 *   7 : Recycling - PFR side     2 D                                      
FFFF
  2000 10001     6     1
 1.00000E+00 0.00000E+00     0     0     0     0     0     0     0     0
FFFT
    -1
FFTFF
     1
     7     1     1     1     1    75    99     0     0
 1.00000E+00 1.40000E+02 4.00000E+00 0.00000E+00 2.11100E+03
    -1     0     0     0     0     0     0
 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
 1.00000E+00 9.00000E+01 0.00000E+00 0.00000E+00 0.00000E+00
 *   8 : Recycling - PFR side     2 C                                      
FFFF
  2000 10001     6     1
 1.00000E+00 0.00000E+00     0     0     0     0     0     0     0     0
FFFT
    -1
FFTFF
     1
     8     1     1     1     1    75    99     0     0
 1.00000E+00 1.40000E+02 4.00000E+00 0.00000E+00 2.11100E+03
    -1     0     0     0     0     0     0
 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
 1.00000E+00 9.00000E+01 0.00000E+00 0.00000E+00 0.00000E+00
 *    9 : Recycling - Outer side   D                                      
FFFF
  2000 10001     6     1
 1.00000E+00 0.00000E+00     0     0     0     0     0     0     0     0
FFFT
    -1
FFTFF
     1
     9     1    37    37    37     1    99     0     0
 1.00000E+00 1.40000E+02 5.00000E+00 0.00000E+00 1.11100E+03
    -1     0     0     0     0     0     0
 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
 1.00000E+00 9.00000E+01 0.00000E+00 0.00000E+00 0.00000E+00
*   10 : Recycling - Outer side   C                                      
FFFF
  2000 10001     6     1
 1.00000E+00 0.00000E+00     0     0     0     0     0     0     0     0
FFFT
    -1
FFTFF
     1
    10     1    37    37    37     1    99     0     0
 1.00000E+00 1.40000E+02 5.00000E+00 0.00000E+00 1.11100E+03
    -1     0     0     0     0     0     0
 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
 1.00000E+00 9.00000E+01 0.00000E+00 0.00000E+00 0.00000E+00
*  11 : Volumetric recombination D                                      
FFFFF
  7000  4001     6     1     0
 1.00000E+00 0.00000E+00     0     0     0     0     0     0     0     0
FFFTF
     1
FFFTF
     1
    11     0     0     0     0     0     0     0     0
 1.00000E+00 1.40000E+02 1.00000E+00 0.00000E+00 0.00000E+00
     0     0     0     0     0     0     0
 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
*  12 : Volumetric recombination C                                      
FFFFF
  7000  4001     6     1     0
 1.00000E+00 0.00000E+00     0     0     0     0     0     0     0     0
FFFTF
     2
FFFTF
     1
    12     0     0     0     0     0     0     0     0
 1.00000E+00 1.40000E+02 2.00000E+00 0.00000E+00 0.00000E+00
     0     0     0     0     0     0     0
 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00

2020-03-27input.dat 7A modules

总结

辐射统计

trajectory

添加靶板区域的source

1 .SOLPS5.0

2. SOLPS-ITER

3. 7A模块的书写逻辑

不改变标签

相关文章

网友评论

延伸阅读

深度阅读

栏目导航

热点阅读