边界条件选择
2020-11-21 本文已影响0人
锅炉工的自我修养
痛点
- 如何选择合适的边界条件
- 保证较小流出功率的基础上,加快运算速度
- 什么样的case,是稳态的case
- 不同边界条件的影响
python面向对象编程的基本原则
- 构想阶段——考虑步骤,不考虑实现(你能想到的事情,编程语言的库基本都帮你实现)
- 实施阶段——逐个完善各个函数的具体实现(易用、参数传递、加深考虑)
- 完善阶段——易读性、巧妙的实现方式(漂亮的代码是改出来的)
构想阶段
(1)输入文件
- b2mn.dat
- b2.boundary.parameters
- b2.neutrals.parameters
- b2.numerics.parameters
- b2.transport.parameters
- b2.transport.inputfile
- input.dat
(2)不同类型
- pure D + no chemical sputering、no pump 、no gas puff 、no drift
- instrisic impurity ——C、sputering、pump、no gas puff 、no drift
- extrisic impurity —— Ne、Ar 、no drift
- extrisic impurity —— Ne、Ar + drift (、dirft )
latex :what you see is what you want
simple latex
complicate latex
实施阶段
(2)不同类型: pass
(3)输入文件
<1>b2.boundary.parameters
- bcene/bceni
S_core:8——指定输入功率——2e6(W)
S_PFR:9——指定衰减长度——0.01m
N_SOL:9——指定衰减长度——0.01m
E_W_target:3——鞘层边界条件,能量传输系数——0.9
未知——enpar(2,2)——1e2
- bcene/bceni
- bcpot
no-target:0——不知道是啥——0
tatget:3——鞘层边界条件,potpar(,2)用于偏压
- bcpot
- 3.bccon
S_CORE:
——1——固定密度——5e19
其他——8——固定流密度——0
S_PFR:
全部——10——泄露因子边界条件——-原子 -0.01,离子 -0.001()
N_SOL:
原子——10——泄露因子——-0.01
离子——9——衰减长度——0.03m
WE_target:
原子——0——未知——0
离子——3——鞘层边界条件——0 - bcmom
S_CORE+PFR:
全为4——混合边界条件,指定速度——mompar(,1)描述平行速度0、mompar(,2)指定速度梯度。
N_SOL:
全为10——平行动量衰减长度——0.03m
WE——target:
原子——1——指定平行速度——mompar(,1)0
离子——3——鞘层边界条件,指定马赫数——mompar(,1)1.00指定马赫数值,mompar(,2)0控制是一起,还是单独的粒子。
- bcmom
好累休息一下,主要工作已经完成,python救我
<2>b2.neutrals.parameters
- WE靶板 、SSS芯部+PFR、N SOL、V 边缘
-
recyc
设为1,其他设为0
-
recyc
-
userfluxparm(1,1)
指定喷你速率
-
userfluxparm(1,1)
<3>b2.numerics.parameters
- 查看粒子种类是否对应,其他值全为1
<4>b2.transport.parameters
- dna transport coefficiention0.3
- hci/hce 热传导系数 1.0
<5>b2.transport.inputfile
- 模拟气球模,跑的贼慢
pass
- 模拟气球模,跑的贼慢
<6>input.dat
- eirene的输入模块
- 6a模块
开化学溅射
6a 参数设置
- 6a模块
- (i)12 2 2 ——打开化学溅射
- (ii)-2.6e-2——指定元素温度为300K(1eV=11600K)
- (iiii)2.0e-2——指定
pump 再循环系数
pump rate
其他 pass
<7>b2mn.dat
主函数,添加开关,选择不同的值,实现不同的功能。
-
使用 inputfile 确定输运系数
transport
-
-
指定内外中平面位置
midplane location
-
-
运算量控制
步数 $\times$ 步长
-
-
产生b2favere
b2favere
-
- feedback 对应开关
-
(i)
(a)feedback 开关(b)isfeedback 设置那种粒子
(c)nesepm_sol 设置边界密度
(d)nesepm_overshoot 设置波动范围
feedback -
(ii)增加gas puff 位置
在中平面增加喷气口,喷气速率 1.0e21 -
(iii)修改b2.boundary.parameters
boundary -
(iv)修改 b2.neutral.parameters 开关
neutral