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
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 |
JP2014521932A (ja) * | 2011-07-15 | 2014-08-28 | オーボテック リミテッド | 電子ビーム誘導プラズマプローブを用いた電子装置の電気検査 |
-
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 | |
Zhang et al. | Sparse target counting and localization in sensor networks based on compressive sensing | |
Fairbairn et al. | A comparison of 4DVar with ensemble data assimilation methods | |
Peiris et al. | Phenomenology of D-brane inflation with general speed of sound | |
CN106443627B (zh) | 基于Bayesian的K分布海杂波形状参数估计方法 | |
Chang | Marginal unscented Kalman filter for cross‐correlated process and observation noise at the same epoch | |
CN108229067B (zh) | 天线-激波管-天线等离子体电磁特性的仿真方法及仿真系统 | |
Morzfeld et al. | Implicit particle filtering for models with partial noise, and an application to geomagnetic data assimilation | |
Riva et al. | A methodology for the rigorous verification of particle-in-cell simulations | |
Boccaletti et al. | The semi-classical approximation at high temperature revisited | |
CN103617348B (zh) | 大气环境电子束等离子体参数获取方法 | |
Xiong et al. | Efficient bias reduction approach of time‐of‐flight‐based wireless localisation networks in NLOS states | |
Amiri et al. | Closed‐form positioning in MIMO radars with antenna location uncertainties | |
Liu et al. | Recovery of high frequency wave fields from phase space–based measurements | |
Schianchi et al. | M 1 neutrino transport within the numerical-relativistic code BAM with application to low mass binary neutron star mergers | |
Qu et al. | Recursive source localisation by time difference of arrival sensor networks with sensor position uncertainty | |
Gambino et al. | Comparison between time-of-flight measurements and numerical simulations for laser-generated plasmas | |
Pinheiro et al. | Efficient nonlinear data assimilation using synchronization in a particle filter | |
Lei et al. | High‐precision hyperboloid location method using passive time‐difference‐of‐arrival measurements | |
Vecherin et al. | Three-dimensional acoustic travel-time tomography of the atmosphere | |
Chen et al. | Overview of cyber-physical temperature estimation in smart buildings: From modeling to measurements | |
Fidler et al. | Suitable initial conditions for Newtonian simulations with massive neutrinos | |
Shang et al. | Approach of system error registration for two‐station coast radars for sea surface monitoring | |
Ferrari et al. | Influence of emitter-receiver number on measurement accuracy in acoustic pyrometry | |
Meng et al. | Joint localization algorithm of TDOA and FDOA considering the position error of underwater sensors |
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 |