美文网首页
2020-03-27input.dat 7A modules

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,\leq 0抽样初始化,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 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

    总结 靶板eirene下标,默认值可能为25000 对于C的追迹,发现测试粒子非常少 通过手动添加,不成功?(化学...

  • 赞美

    7A上下为主勤劳!

  • 7A

    [太阳][玫瑰]家庭教育每日分享:家庭教育中的“7A”原则 “7A”原则,是由美国著名婚姻与子女教育专家约瑟·D·...

  • 2.4.4 modules配置

    2.4.4 modules配置 问题一:modules配置? resolve.modules配置 Webpack ...

  • 嘟嘟莱雅

    嘟嘟莱雅是7A优品自营的一款品牌。

  • A优品创始人王正:社交电商的下半场

    7A优品——美好生活如7而至 长青(中国)7A优品创始人王正 2019(第二届)中国社交电商新零售产业峰会长青(中...

  • C++ Modules 初窥

    C++ Modules Modules 被认为是C++20的BigFour特性(Modules,Concept,R...

  • AS出现Unsupported modules detected

    AS出现Unsupported modules detected项目中碰到Unsupported modules ...

  • 7A优品

    7A优品是由长青(中国)日用品有限公司斥资1亿元人民币与广州美佳慧科技有限公司联合打造的新型轻社交电商平台,致力于...

  • Modules

    module_using_sys.py module.using_name.py mymodule.py mymo...

网友评论

      本文标题:2020-03-27input.dat 7A modules

      本文链接:https://www.haomeiwen.com/subject/jwukuhtx.html