美文网首页
2020-02-24b2mn.dat.csf_formatted

2020-02-24b2mn.dat.csf_formatted

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

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 中性漂移选项
    1. 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

修改喷气速率


    1. b2mn.dat 整合(cfs and dyl)

    按功能分块
    整合

    1. b2.boundary.parameter保持
    1. b2ag.dat 修改
    1. b2.numerics.parameter speedup
    1. 尝试运行

相关文章

网友评论

      本文标题:2020-02-24b2mn.dat.csf_formatted

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