PASSKEy输入文件的问题（原飞书群问题迁移，20210315）

admin

朱老师您好：
根据您上次的建议，我正尝试使用您开发的PASSKEy程序。在输入文件设置方面，我遇到如下问题想向您请教。
PASSKEy.inp

• 在&reactions_plasma中输入反应方程时，反应中的三体m应该怎么表达？
• 反应焓变应该怎么计算或者获取？
• 反应速率的输入应该用什么单位？
  

朱：

• 手工把对应的中性组分填进去，有多少中性成分就写多少个反应
• 从文献获取，或者计算能级差减去断键的能量
• 国际标准单位，m6/s,m3/s或者/s
  

Coefficients.input

• 我们有一台工作站是双cpu，总共28核56线程，&MPIOMP中的num_cores和num_processes是不是应该分别设为56和2呢？
  

朱：

• num_cores是线程数，填56
• num_processes其实是节点数，一台工作站的话填0或1就好了，MPI没开的时候这个选项没用
  

udf_geometry.f90

• f()和g()是什么参数，我没有理解为什么要用单元中一条边的一个点的y坐标减去另一个点x坐标的函数。
  

PASSKEy程序中是怎么界定neutral gas和plasma的？example文件中，在定义material时，先把所有单元定义成neutral gas后又定义成了plasma，其实neutral gas是被plasma覆盖掉了的吧。
   3. cell()%Y和cell()%boundary_cell()%X分别指的是哪里的Y和X坐标？在定义driven electrode和ground electrode的时候，example中给出的电极坐标范围分别是用cell()%Y和cell()%boundary_cell()%X来约束的。

   4. 如果要定义SDBD的外形，按照程序中的思路，是不是应该先把全域定义plasma，然后定义电介质层，最后给出接电极和接地极的范围呢？

朱：

• 这是之前对应Kulikovsky流注几何设置的两个参考函数，用来划定针尖电极的，可以忽略不计，只是做参考的
• Neutral gas是只计算泊松方程、流体方程、燃烧等的，plasma部分不仅计算这些，还计算等离子体输运方程和反应方程。案例中plasma全部覆盖了neutral gas，其实可以不用，这样计算速度会大大提升。
• cell()%Y是指网格中心点的Y坐标，cell(i)%boundary_cell(j)%X是指网格i的第j条边中心位置的X坐标，网格的边的编号是从底部开始逆时针从1到4，如下图所示：
  

几何形状约束方法比较自由，只要材料覆盖到就行了。


4. 应该先定义一个neutral gas区域，再定义一个比较小的等离子体区域，然后介质，电极。参照2017年的PSST给出的示意图。

hylyl

pumpkin计算等离子体反应路径的贡献损失率，有人会用这个软件吗？