b2mn.dat_cf
Physics modules
*label (lblmn: character*60)
'D3D : no drifts : D+C : 152835.2500 (H-mode) : B2.5-Eirene'
*b2cmpa basic parameters
*b2cmpb boundary conditions
*cbregs
'south' 3 0.25 0.75 2.00
'north' 1 2.0
'west' 1 2.0
'east' 1 2.0
*b2cmpt transport coefficients
*cflim cflme cflmi cflmv ...
'-1' 0.300 0.300 0.5 0 0 0 0 0
*endphy
-2. RUN MODULES
'b2mndr_ntim' '5000'
'b2mndr_dtim' '5.0e-6'
'b2mndt_nstg0' '1'
'b2mndt_nstg1' '1'
'b2mndt_nstg2' '10'
'b2stbc_boundary_namelist' '1' #
'b2srdt_numerics_namelist' '1'
'b2tqna_transport_namelist' '1'
'b2tqna_inputfile' '1'
'b2mndr_na0eps' '1e14' # 所有网格,粒子初始密度,5.x,在iter中为b2mndr_na_new
'b2mndr_cpu' '1828000.0' # run CPU限制,置零无限制,其他,是上限。
'b2mndr_savecpu' '3600.0' # CPU时间间隔,保存plasmastate.***,中断备份
'b2mndr_eirene' '1'
'b2mndr_rescale_neutrals' '1.0' # 第一步中性密度乘子
'b2mndr_rescale_neutrals_sources' '1.0e-10' # 中性源乘子
-2. OUTPUT module
无 netcdf 和 balance_average
无 b2time.nc
无 tally
-3. GEOMETRY module
'b2mwti_jxi' '38'
'b2mwti_jxa' '54'
-4. 数值,物理,原子模块
-5. BOUND module
无 min_ares*
'b2stbc_sna0ep' '1.0e-10' # 引入网格的低等级源
'b2trcl_lluciani' '1' # physics 1 连接长度。2,反向兼容4.0。3,lluciani新的SPB格式输运系数
'b2news_potit' '100' # 电势方程计算的迭代最大次数。
# 增加开关
'b2news_potok' '1e-2'
无 recalculate
'eirene_repeat_first_call' '1' # 大于零,重复调用eirene。如果密度feedback gas puff打开,将会有用。
无extrap
DRIFT module
'b2news_facdrift_target' '0.00' # 抗磁漂移
'b2news_facdrift_start' '0.00'
'b2news_facdrift_inc' '1.01'
'b2news_facdrift_dec' '0.95'
'b2news_facExB_target' '0.00'
'b2news_facExB_start' '0.00'
'b2news_facExB_inc' '1.01'
'b2news_facExB_dec' '0.95' # EXB漂移
无粘滞漂移
'b2tfhe_vis_par' '0.0' # 平行粘滞驱动电流乘子
b2tfhe_vis_par一起使用(防止二次计算经典粘滞效应)
'b2trcl_lthf21' '0' # physics 非零,热传导模型通过21矩近似修改,仅用于多流体运行。
'b2tfhe_neutral' '0.0' # physics 非零打开离子——中性电流,需要和EXB和抗磁漂移配合使用
无 vis_q
-6. VISCOSITY TERM
无 addvis* 粘滞张量贡献
无hz
误sigp SPB压强梯度项
- 7. FLUX LIMIT
无 limit_to_dpa 压强驱动中性系数
无 b2tlc0z_alpha/gamma
无 lim_flux 电子热流&离子热流
无conductive_limit 离子热流流限制
-8. NEW manner of calculate drift
无b2fnb_drift 中性漂移选项
- About solver
无 9点矩阵计算
无 7点矩阵计算
无vis_per 垂直粘滞乘子
'b2ux5p_nltrsol' '0' # Output 输出矩阵计算器标签,越大输出等级越高
'b2ux5p_style' '2' # numerics 5点模板矩阵计算器
-10. SOME RECOMMENDED KEYS ABOUT DRIFT
无 zhadanov
无 pTe 新的电子热流计算
无 alfTeEh 旧的电子热流计算
无 phm0 离心力项乘子
无 lluciani 3 SPB格式lluciani系数
-11. SOME SWITCH FORM
无fchi 确定采用的离子流表达式
无 model_sig 计算反常系数是用中平面还是平均
无 dt_av 大于零计算平均
无 dt_av_continue 大于零计算continue
专属项
'b2trcl_lvis21' '0' # physics,非零,通过21-矩,只被用于多流体运算。不应该和
'b2sqel_artificial_radiation' '0.0' # 非零,把手动辐射损失项加入电子冷却率。art_rad取代等离子体中的杂质部分
'b2stbc_fchy_dia_coreonly' '1' # 置一,只在芯部使用新经典电流和抗磁热流对流边界条件。置零,引用到所有区域
'b2tlh0_alpha' '1.0' #中性离子热流传导流限制。经典流限制乘子。数值越大,流限制越小。
'b2tfnb_alpha' '1.0' # 指定对流中性流流限制参数。
'b2news_potitmin' '5' # 电势返程最小迭代次数
'b2stbr_neutrals_namelist' '1'
'b2stbc_feedback' '0' # 为边界条件开启feedback模式。
'b2stbc_diagno' '1' # 控制b2stbc的输出等级及其从属例程。
'b2stbc_solregno' '3' # b2ah.dat 和b2mn.dat中指定,SOL north 边界条件的下标,SN 对应3,DN对应7。
'eirene_dpc_fix' '1' # 置一,使用真实的粒子源。置二,粒子源置零。
'eirene_ionizing_core' '1' # 置一,被添加的流,直接修改边界网格,只有标准流边界条件被使用时工作;小于零,在边界条件为13时,边界条件使用绝对值。
'b2news_area_fix' '3' # 置零,恢复5.0,置一,计算极向速度。置二,使用平行连接区域的处理方式。置三,使用所有的平行接触区域。
b2ag.dat
b2tallies.nc存贮什么信息?
如何指定正场和反场?——线索(b2ag.dat)
*dimens (nx, ny, nx1, ny1; free format)
96 36 96 36
*param (param(0:99); free format)
-1.0, 99*0.0
'b2agfs_min_pitch' '0.01'
'b2agfs_geometry' 'F075075.004011-3.5s_v1.v001.geo'
'b2agfs_Bt_reversal' '0' # 置一,环向磁场翻转。在b2ag.dat中使用。
-3. b2.boundary
&BOUNDARY
nbc= 6,
bcchar= 'S', 'W', 'E', 'S', 'S', 'N',
bcpos= -1, -1, 96, -1, -1, 36,
bcstart= 24, -1, -1, -1, 72, -1,
bcend= 71, 36, 36, 23, 96, 96,
bcene= 8, 3, 3, 9, 9, 9,
bceni= 8, 3, 3, 9, 9, 9,
bcpot= 12, 11, 11, 2, 2, 2,
enepar(1,1)= 1.750E+06, 1.000 , 1.000 , 3.000E-02, 3.000E-02, 3.000E-02,
enipar(1,1)= 1.750E+06, 2.50 , 2.50 , 3.000E-02, 3.000E-02, 3.000E-02,
enipar(1,2)= 0.00 , 0.00 , 0.00 , 0.00 , 0.00 , 0.00 ,
potpar(1,1)= 0.00 , 0.00 , 0.00 , 0.00 , 0.00 , 0.00 ,
bccon(0,1)= 8, 1,
bccon(0,2)= 5, 14,
bccon(0,3)= 5, 14,
bccon(0,4)= 10, 10,
bccon(0,5)= 10, 10,
bccon(0,6)= 10, 10,
conpar(0,1,1)= 0.00 , 3.4000E+19,
conpar(0,2,1)= 0.00 , 1.0,
conpar(0,3,1)= 0.00 , 1.0,
conpar(0,4,1)=-1.000E-02, -1.0e-3,
conpar(0,5,1)=-1.000E-02, -1.0e-3,
conpar(0,6,1)=-1.000E-02, -1.0e-3,
bcmom(0,1)= 2, 2,
bcmom(0,2)= 13, 13,
bcmom(0,3)= 13, 13,
bcmom(0,4)= 2, 2,
bcmom(0,5)= 2, 2,
bcmom(0,6)= 2, 2,
mompar(0,1,1)= 0.00 , 0.00 ,
mompar(0,2,1)= 0.00 , 1.000 ,
mompar(0,3,1)= 0.00 , 1.000 ,
mompar(0,4,1)= 0.00 , 0.00 ,
mompar(0,5,1)= 0.00 , 0.00 ,
mompar(0,6,1)= 0.00 , 0.00 ,
mompar(0,2,2)= 0.00 , 0.00 ,
mompar(0,3,2)= 0.00 , 0.00 ,
lcbs= 1,
nniso= 0,
/
b2.boundary_use
&BOUNDARY
nbc= 6,
bcchar= 'S', 'S', 'S', 'N', 'W', 'E',
bcpos= -1, -1, -1, 36, -1, 96,
bcstart=24, 0, 72, 0, -1, -1,
bcend= 71, 23, 95, 95, 36, 36,
bcene= 8, 9, 9, 9, 3, 3,
bceni= 8, 9, 9, 9, 3, 3,
bcpot= 0, 0, 0, 0, 3, 3,
enepar(1,1)= 1.5e6, 0.01, 0.01, 0.01, 0.90, 0.90, enepar(2,2)=1e2, 1e2, 1e2,
enipar(1,1)= 1.5e6, 0.01, 0.01, 0.01, 2.50, 2.50, enepar(2,2)=1e2, 1e2, 1e2,
enipar(1,2)= 0.00 , 0.00 , 0.00 , 0.00 , 0.00 , 0.00 ,
potpar(1,1)= 0.00 , 0.00 , 0.00 , 0.00 , 0.00 , 0.00 , potpar(1,2)=6*0,
bccon(0,1)= 8, 1, 8, 8, 8, 8, 8, 8, 8,
bccon(0,2)= 10, 10, 10, 10, 10, 10, 10, 10, 10,
bccon(0,3)= 10, 10, 10, 10, 10, 10, 10, 10, 10,
bccon(0,4)= 10, 10, 10, 10, 10, 10, 10, 10, 10,
bccon(0,5)= 0, 3, 0, 3, 3, 3, 3, 3, 3,
bccon(0,6)= 0, 3, 0, 3, 3, 3, 3, 3, 3,
conpar(0,1,1)= 0.00 , 1.000E+19, 0.00 , 0.00 , 0.00 , 0.00 , 0.00 , 0.00 , 0.00 ,
conpar(0,2,1)=-1.000E-04, -1.000E-04, -1.000E-04, -1.000E-03, -1.000E-03, -1.000E-03, -1.000E-03, -1.000E-03, -1.000E-03,
conpar(0,3,1)=-1.000E-04, -1.000E-04, -1.000E-04, -1.000E-03, -1.000E-03, -1.000E-03, -1.000E-03, -1.000E-03, -1.000E-03,
conpar(0,4,1)=-1.000E-04, -1.000E-04, -1.000E-04, -3.000E-03, -3.000E-03, -3.000E-03, -3.000E-03, -3.000E-03, -3.000E-03,
conpar(0,5,1)= 0.000E-02, 0.000E-02, 0.000E-02, 0.000E-02, 0.000E-02, 0.000E-02, 0.000E-02, 0.000E-02, 0.000E-02,
conpar(0,6,1)= 0.000E-02, 0.000E-02, 0.000E-02, 0.000E-02, 0.000E-02, 0.000E-02, 0.000E-02, 0.000E-02, 0.000E-02,
bcmom(0,1)= 1, 4, 1, 4, 4, 4, 4, 4, 4,
bcmom(0,2)= 10, 10, 10, 10, 10, 10, 10, 10, 10,
bcmom(0,3)= 10, 10, 10, 10, 10, 10, 10, 10, 10,
bcmom(0,4)= 10, 10, 10, 10, 10, 10, 10, 10, 10,
bcmom(0,5)= 1, 3, 1, 3, 3, 3, 3, 3, 3,
bcmom(0,6)= 1, 3, 1, 3, 3, 3, 3, 3, 3,
mompar(0,1,1)= 0.00 , 0.00 , 0.00 , 0.00 , 0.00 , 0.00 , 0.00 , 0.00 , 0.00 ,
mompar(0,2,1)= 0.03 , 0.03 , 0.03 , 0.03 , 0.03 , 0.03 , 0.03 , 0.03 , 0.03 ,
mompar(0,3,1)= 0.03 , 0.03 , 0.03 , 0.03 , 0.03 , 0.03 , 0.03 , 0.03 , 0.03 ,
mompar(0,4,1)= 0.03 , 0.03 , 0.03 , 0.03 , 0.03 , 0.03 , 0.03 , 0.03 , 0.03 ,
mompar(0,5,1)= 0.00 , 1.00 , 0.00 , 1.00 , 1.00 , 1.00 , 1.00 , 1.00 , 1.00 ,
mompar(0,6,1)= 0.00 , 1.00 , 0.00 , 1.00 , 1.00 , 1.00 , 1.00 , 1.00 , 1.00 ,
mompar(0,1,2)= 1e05 , 1e05 , 1e05 , 1e05 , 1e05 , 1e05 , 1e05 , 1e05 , 1e05 ,
mompar(0,2,2)= 1e05 , 1e05 , 1e05 , 1e05 , 1e05 , 1e05 , 1e05 , 1e05 , 1e05 ,
mompar(0,3,2)= 1e05 , 1e05 , 1e05 , 1e05 , 1e05 , 1e05 , 1e05 , 1e05 , 1e05 ,
mompar(0,4,2)= 1e05 , 1e05 , 1e05 , 1e05 , 1e05 , 1e05 , 1e05 , 1e05 , 1e05 ,
mompar(0,5,2)= 0.00 , 0.00 , 0.00 , 0.00 , 0.00 , 0.00 , 0.00 , 0.00 , 0.00 ,
mompar(0,6,2)= 0.00 , 0.00 , 0.00 , 0.00 , 0.00 , 0.00 , 0.00 , 0.00 , 0.00 ,
lcbs= 1,
nniso= 0,
GAMMAI=1.666666667,GAMMAE=0.5,
LBNDUSR=f,
/
b2.neutrals
修改喷气速率
-
- b2mn.dat 整合(cfs and dyl)
按功能分块
整合 - b2.boundary.parameter保持
- b2ag.dat 修改
- b2.numerics.parameter speedup
- 尝试运行
网友评论