1996R. Marchand (CPC)

Abstract

• 1. 自动产生曲线准正交网格的计算代码
• 2. 平衡产生考虑structure和边界
• 3. 支持实验与模拟对比

---

Method of solution

• 1. 通过参数化等高线（指定流函数，contourc for matlab）
• 2. 分布格点，通过松弛迭代计算临近格点位置

---

正交网格参考文献

• 1. 1982J. F. Thompson

https://books.google.co.jp/books?hl=zh-CN&lr=&id=ImaDT6ijKq4C&oi=fnd&pg=PR9&dq=J.F+Thompson+numerical+grid+generation+1982+pdf&ots=abWWJCtPkh&sig=UZN80nh_F6Gd-cTPHc2CUmrOufM#v=onepage&q&f=true

---

考虑的要素

• 1. XP（如何识别）
• 2. divertor plate（如何确定限制的位置）
• 3. separatrix 如何识别？？（如何让边界线连接在一起）
• 4. 计算分区

---

元素排序原则

• 1. 从上到下（相对于OMP）
• 2. 耳机排序：upper or lower
• 3. 拓扑：right or left

---

• 4. XP：
• 5. divertor plates
  • SN：靶板从右到左排序；
  • DN：上左，上右，下左，下右
  • 靶版需要特殊识别，通过标记为不同的结构
• 6. segments of separatrix
  • XP作为分割点
  • XP归属通过临近原则确定
  • 对于DDN，内内层闭环，需要分开成两半
• 7. region
  • r1:连接靶版的外侧区域，对于DDN是inner separatrix右侧的区域
  • r2:远离右侧靶版的区域
    image.png

image.png
image.png
image.png

---

methodology

1.步骤

• 1. 识别磁场位型
• 2. 找到 separtrix 和 boundaries，分区
• 3. 指定初始flux的值
• 4. 把separatrix通过格点进行划分
• 5. 将格点从一个flux映射到另一个表面

---

2. 依赖条件

• 1. 找到靶板和限制器位置

• 2. 获得对应位置的极向磁场

     
     image.png
• 3. 极向磁场通过flux计算
  • 求gradient
  • need to be done
• 4. structure:(r,z (mm))

---

3. 定义磁场位型

• 1. 计算磁通的导数，找打导数为零的点
  • OP：磁轴
  • XP
  • 用户手动选择XP和OP（多个导数名零点）
  • 对于两个零点，特别是两个零点距离很近的情况
• 2. spearatrix， boundaries以及区域划分
• 通过XP对segments进行划分，确定位置
• 计算通过XP的等高线（知道XP位置的磁通值）

如何处理，通过XP的等高线不在XP交叉的问题？？

• 3. 子边界线计算
  • 从XP到与靶板的交叉点
  • 同时定义偏滤器plates
  • 挡板是边界线从XP远离方向的第一个交叉点
• 4. 寻找outer boundaries of SOL and PFR
  • 判断磁通是否与初始靶板以为的元素交叉，最后一个不与其他structure交叉的位置boundary
  • 沿着板步进

如何判断断交叉
https://blog.csdn.net/s77108887/article/details/7911953?utm_medium=distribute.pc_relevant.none-task-blog-BlogCommendFromMachineLearnPai2-17.nonecase&depth_1-utm_source=distribute.pc_relevant.none-task-blog-BlogCommendFromMachineLearnPai2-17.nonecase
https://blog.csdn.net/u010429167/article/details/73187854

• done

如何找到跨越靶板的等高线部分？

如何获得排好序的登高线数据？

查看fortran 源码

• 1. 数组：f(a,b)

---

需求

• 1. btf
• 2. r=eps
• 3. 全磁场

---

• 5. 如何确定芯部厚度
  • XP towads OP方向的厚度
  • CDN 上XP
  • DDN inner XP

---

如何获得MATLAB contour的顺序

如何获得登高线数据的排列顺序

---

• 6. 边界线上的格点分布

     
     藏宝图