美文网首页
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