CN103617348B - 大气环境电子束等离子体参数获取方法 - Google Patents
大气环境电子束等离子体参数获取方法 Download PDFInfo
- Publication number
- CN103617348B CN103617348B CN201310568196.XA CN201310568196A CN103617348B CN 103617348 B CN103617348 B CN 103617348B CN 201310568196 A CN201310568196 A CN 201310568196A CN 103617348 B CN103617348 B CN 103617348B
- Authority
- CN
- China
- Prior art keywords
- electron beam
- thin layer
- temperature
- energy loss
- gas
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000010894 electron beam technology Methods 0.000 title claims abstract description 54
- 238000000034 method Methods 0.000 title claims abstract description 31
- 230000007613 environmental effect Effects 0.000 title claims abstract description 13
- 230000005540 biological transmission Effects 0.000 claims abstract description 15
- 238000012937 correction Methods 0.000 claims description 7
- 230000009471 action Effects 0.000 claims description 5
- 230000005611 electricity Effects 0.000 claims description 2
- 230000008569 process Effects 0.000 abstract description 6
- 230000008878 coupling Effects 0.000 abstract description 2
- 238000010168 coupling process Methods 0.000 abstract description 2
- 238000005859 coupling reaction Methods 0.000 abstract description 2
- 238000011160 research Methods 0.000 abstract description 2
- 239000002245 particle Substances 0.000 abstract 1
- 238000004088 simulation Methods 0.000 description 4
- 239000012530 fluid Substances 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 230000006978 adaptation Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000006353 environmental stress Effects 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 238000013316 zoning Methods 0.000 description 1
Landscapes
- Plasma Technology (AREA)
Abstract
一种大气环境电子束等离子体参数获取方法，其特点在于：采用多层网格剖分和蒙特卡罗方法处理电子束在大气中的输运和能量损失，考虑了能量损失对气体密度的影响，采用预估‑校正方法动态耦合电子束能量，采用动态扫描方法实现对高能粒子的跟踪。本发明解决了大气环境电子束等离子体参数获取中的多尺度和多物理场难题，并考虑能量损失对传输路径上气体密度的影响，实现了对大气环境电子束等离子体的准确研究。
Description
技术领域
本发明涉及一种电子束等离子体参数获取方法。
背景技术
中国期刊《物理学报》2006年第55卷第7期中，题为《电子束产生大尺度等离子体过程的数值模拟研究》中建立了一个四组分一维混合模型，对电子束注入大气产生大尺度等离子体的过程进行了数值模拟。该模型采用流体模型模拟电子束的传输过程，并以此计算电子束在气体中的能量损失，依据能量损失获得电子束的电离率，进而计算电子束等离子体的参数。这是一种模拟计算电子束等离子体参数的解析方法，可以快速获得等离子体参数，但其不足之处在于：流体模型采用简单的公式描述电子束的传输过程，忽略了电子束与气体原子、分子之间复杂的碰撞过程，且没有考虑由于碰撞导致的能量损失而引起的气体参数变化，因而造成计算得到的电子束传输距离和等离子体尺寸与实验值具有较大差异。
发明内容
为克服现有采用流体模型跟踪电子束的传输过程而导致获取的等离子体尺寸与实验值具有较大差异的技术缺陷，本发明提供一种大气环境电子束等离子体参数获取方法。
本发明的技术解决方案如下：
一种大气环境电子束等离子体参数获取方法，其特殊之处在于：包括以下步骤：
1】建立电子束的传输与碰撞模型；
1.1】利用蒙特卡罗方法建立电子束的传输与碰撞模型；
1.2】在步骤1.1】所建立的模型中采用多层网格剖分的方法将电子束与空气的作用区域划分为若个气体薄层，各个气体薄层的参数取决于气体压强与温度；
2】预估各个气体薄层的温度，并读取大气环境压强；
3】校准各个气体薄层的温度，具体步骤包括：
3.1】将步骤2】中预估的各气体薄层的温度及大气环境压强代入步骤1】所建立的传输与碰撞模型中，利用输入的预估温度数据初始化气体薄层的参数，利用该模型动态扫描跟踪电子束在气体薄层中的传输过程，并实时记录电子束的能量损失数据；
3.2】利用能量损失数据获取电子束作用区域的温度，并以此为校正温度数据取代前一次输入传输与碰撞模型中的温度数据，将其代入步骤1】所建立的传输与碰撞模型中，重新初始化气体薄层的参数；
利用该模型动态扫描跟踪电子束在气体薄层中的传输过程，再次记录电子束的能量损失数据；
3.3】将相邻两次的能量损失数据进行比较；
3.4】如果两者相对误差小于1％，则进行步骤4】，如果两者相差过大，则重复步骤3.2】-3.3】直至相邻两次的能量损失数据相对误差小于1％。
4】根据步骤3】所校正的能量损失数据或各个气体薄层的温度计算等离子体参数。
上述步骤1】中气体薄层的数量为1000个，且各气体薄层的参数取决于大气环境压强及电子束能量损失引起的温度变化。
上述步骤2】中各个气体薄层的温度采用一次预估温度数据。
上述步骤3】中记录各薄层中的电子能量损失，电子束传输1mm记录一次。
上述步骤3】中利用能量损失数据迭代计算校正温度数据，迭代计算次数为6-10次。
本发明相对于现有技术的优点是：
1、本发明通过采用蒙特卡罗方法，可以准确处理电子束与气体原子、分子的碰撞，大幅提高了参数获取的准确性。
2、本发明采用多层网格剖分对电子束的作用区域进行了气体薄层处理，动态耦合能量损失引起的温度变化数据，准确反映了电子束与气体作用的物理过程，实现了多物理场耦合。
3、本发明采用动态扫描方法跟踪电子束的高能电子，从而可以实时获得电子束任意输运距离处的能量损失。
4、本发明具有适应各种气体环境中电子束等离子体参数获取的特色，可为理论研究和工程研制提供有效手段。
附图说明
图1为本发明的流程图。
具体实施方式
本发明的大气环境电子束等离子体参数获取方法包括电子束的传输与碰撞模型的建立的步骤、利用该模型动态扫描跟踪电子束在气体薄层中的传输过程，记录电子束的能量损失数据的步骤、利用能量损失数据对气体温度进行校准的步骤、根据校准后的温度或能量损失数据获取等离子参数的步骤。
一个具体实例如图1所示，包括以下步骤：
(1)采用电子束加热模型，得到预估温度场，将温度数据输入到蒙特卡罗能量损失计算模块1，电子束作用区划分为1000个气体薄层，每个气体薄层厚度为1mm，利用输入的预估温度数据初始化气体薄层的参数(密度)；
(2)蒙特卡罗模型动态扫描跟踪1000万个电子束高能电子，记录电子束在传输时的能量损失数据；
(3)利用能量损失，采用温度校正模块2获取电子束的校正温度场；
(4)将校正温度场数据输入到蒙特卡罗能量损失计算模块1中，电子束作用区划分为1000个气体薄层，每个气体薄层厚度为1mm，利用输入的校正温度数据重新初始化气体薄层的参数(密度)；
(5)再次采用蒙特卡罗模型动态扫描跟踪1000万个电子束高能电子，记录电子束在传输时的能量损失数据；
(6)判断能量损失数据相对误差是否小于1％，若是，则进入电子束等离子体参数模块3；若否，则继续进入温度校正模块2校正温度场，重新获取能量损失数据；
(7)根据能量损失数据，采用电子束等离子体参数模块3获取电子束等离子体参数。
Claims (5)
1.一种大气环境电子束等离子体参数获取方法，其特征在于：包括以下步骤：
1】建立电子束的传输与碰撞模型；
1.1】利用蒙特卡罗方法建立电子束的传输与碰撞模型；
1.2】在步骤1.1】所建立的模型中采用多层网格剖分的方法将电子束与空气的作用区域划分为若个气体薄层，各个气体薄层的参数取决于气体压强与温度；
2】预估各个气体薄层的温度，并读取大气环境压强；
3】校准各个气体薄层的温度，具体步骤包括：
3.1】将步骤2】中预估的各气体薄层的温度及大气环境压强代入步骤1】
所建立的传输与碰撞模型中，利用输入的预估温度数据初始化气体薄层的参数，
利用该模型动态扫描跟踪电子束在气体薄层中的传输过程，并实时记录电子束的能量损失数据；
3.2】利用能量损失数据获取电子束作用区域的温度，并以此为校正温度数据取代前一次输入传输与碰撞模型中的温度数据，将其代入步骤1】所建立的传输与碰撞模型中，重新初始化气体薄层的参数；
利用该模型动态扫描跟踪电子束在气体薄层中的传输过程，再次记录电子束的能量损失数据；
3.3】将相邻两次的能量损失数据进行比较；
3.4】如果两者相对误差小于1％，则进行步骤4】，如果两者相差过大，则重复步骤3.2】-3.3】直至相邻两次的能量损失数据相对误差小于1％；
4】根据步骤3】所校正的能量损失数据或各个气体薄层的温度计算等离子体参数。
2.根据权利要求1所述的大气环境电子束等离子体参数获取方法，其特征在于：所述步骤1】中气体薄层的数量为1000个，且各气体薄层的参数取决于大气环境压强及电子束能量损失引起的温度变化。
3.根据权利要求1所述的大气环境电子束等离子体参数获取方法，其特征在于：所述步骤2】中各个气体薄层的温度采用一次预估温度数据。
4.根据权利要求1所述的大气环境电子束等离子体参数获取方法，其特征在于：所述步骤3】中记录各薄层中的电子能量损失，电子束传输1mm记录一次。
5.根据权利要求1所述的大气环境电子束等离子体参数获取方法，其特征在于：所述步骤3】中利用能量损失数据迭代计算校正温度数据，迭代计算次数为6-10次。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310568196.XA CN103617348B (zh) | 2013-11-13 | 2013-11-13 | 大气环境电子束等离子体参数获取方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201310568196.XA CN103617348B (zh) | 2013-11-13 | 2013-11-13 | 大气环境电子束等离子体参数获取方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN103617348A CN103617348A (zh) | 2014-03-05 |
| CN103617348B true CN103617348B (zh) | 2016-08-17 |
Family
ID=50168051
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201310568196.XA Active CN103617348B (zh) | 2013-11-13 | 2013-11-13 | 大气环境电子束等离子体参数获取方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103617348B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104866649B (zh) * | 2015-04-27 | 2017-09-26 | 北京计算科学研究中心 | 等离子体器件的光电转换效率多尺度建模和检测方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102033994A (zh) * | 2010-12-07 | 2011-04-27 | 北京航空航天大学 | 基于马尔可夫链蒙特卡罗的舵机可靠性仿真抽样方法 |
| CN102135628A (zh) * | 2010-12-31 | 2011-07-27 | 中国航天科技集团公司第六研究院第十一研究所 | 用于大气环境下电子束等离子体的诊断探针及诊断方法 |
| CN102507608A (zh) * | 2011-09-29 | 2012-06-20 | 西安空间无线电技术研究所 | 一种确定中低能电子非弹性散射的方法 |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5870687A (en) * | 1996-04-16 | 1999-02-09 | Majstorovic; Jovo Bozidar | Vehicle performance analyzer employing an accelerometer-based user interface |
| DE19725450A1 (de) * | 1997-06-16 | 1998-12-17 | Hoechst Schering Agrevo Gmbh | 4-Haloalkyl-3-heterocyclylpyridine und 4-Haloalkyl-5-heterocyclylpyrimidine, Verfahren zu ihrer Herstellung, sie enthaltende Mittel und ihre Verwendung als Schädlingsbekämpfungsmittel |
| KR20140057243A (ko) * | 2011-07-15 | 2014-05-12 | 오르보테크 엘티디. | 전자-빔 유도형 플라스마 프로브를 이용하는 전자 소자의 전기적 검사 |
-
2013
- 2013-11-13 CN CN201310568196.XA patent/CN103617348B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102033994A (zh) * | 2010-12-07 | 2011-04-27 | 北京航空航天大学 | 基于马尔可夫链蒙特卡罗的舵机可靠性仿真抽样方法 |
| CN102135628A (zh) * | 2010-12-31 | 2011-07-27 | 中国航天科技集团公司第六研究院第十一研究所 | 用于大气环境下电子束等离子体的诊断探针及诊断方法 |
| CN102507608A (zh) * | 2011-09-29 | 2012-06-20 | 西安空间无线电技术研究所 | 一种确定中低能电子非弹性散射的方法 |
Non-Patent Citations (3)
| Title |
|---|
| 低真空中高能电子束传输特性的MC模拟;栾希亭等;《核技术》;20101210;第33卷(第2期);全文 * |
| 大气环境电子束材料表面强化机理数值模拟;邓永锋;《火箭推进 》;20110415;第37卷(第2期);1-5 * |
| 电子束空气等离子体电子数密度分布模型研究;张志豪;《高电压技术 》;20130731;第39卷(第7期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN103617348A (zh) | 2014-03-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Kravtsov et al. | A new robust low-scatter X-ray mass indicator for clusters of galaxies | |
| Peiris et al. | Phenomenology of D-brane inflation with general speed of sound | |
| Chang | Marginal unscented Kalman filter for cross‐correlated process and observation noise at the same epoch | |
| Morzfeld et al. | Implicit particle filtering for models with partial noise, and an application to geomagnetic data assimilation | |
| CN105790258B (zh) | 基于正态Copula函数的拉丁超立方抽样法概率潮流计算方法 | |
| Paschalidis et al. | Head-on collisions of binary white dwarf-neutron stars: Simulations in full general relativity | |
| Desjacques et al. | Primordial Non‐Gaussianity in the Large‐Scale Structure of the Universe | |
| CN107765103A (zh) | 一种基于多传感器的复杂环境电磁态势反演方法 | |
| Boccaletti et al. | The semi-classical approximation at high temperature revisited | |
| Schianchi et al. | M 1 neutrino transport within the numerical-relativistic code BAM with application to low mass binary neutron star mergers | |
| CN103617348B (zh) | 大气环境电子束等离子体参数获取方法 | |
| CN118759514A (zh) | 一种多目标跟踪方法 | |
| Liu et al. | Recovery of high frequency wave fields from phase space–based measurements | |
| CN105183997B (zh) | 一种基于双层嵌套不确定性传播的热传导模型校准方法 | |
| CN111553110B (zh) | 一种基于蒙特卡罗模拟的核电磁脉冲电流源项计算方法 | |
| Pinheiro et al. | Efficient nonlinear data assimilation using synchronization in a particle filter | |
| Meng et al. | Joint localization algorithm of TDOA and FDOA considering the position error of underwater sensors | |
| Lei et al. | High‐precision hyperboloid location method using passive time‐difference‐of‐arrival measurements | |
| Chen et al. | Overview of cyber-physical temperature estimation in smart buildings: From modeling to measurements | |
| CN103322955B (zh) | 一种反向求解芯片单热源位置和面积的方法 | |
| CN105528735A (zh) | 基于量测风速与空间相关性的异常数据点的校正方法 | |
| Ferrari et al. | Influence of emitter-receiver number on measurement accuracy in acoustic pyrometry | |
| Fidler et al. | Suitable initial conditions for Newtonian simulations with massive neutrinos | |
| CN111277950B (zh) | 一种基于到达时差和到达频差的定位方法、装置及设备 | |
| Shang et al. | Approach of system error registration for two‐station coast radars for sea surface monitoring |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant |