CN103048522B - 常压下低温等离子体密度参数的诊别方法 - Google Patents
常压下低温等离子体密度参数的诊别方法 Download PDFInfo
- Publication number
- CN103048522B CN103048522B CN201310010841.6A CN201310010841A CN103048522B CN 103048522 B CN103048522 B CN 103048522B CN 201310010841 A CN201310010841 A CN 201310010841A CN 103048522 B CN103048522 B CN 103048522B
- Authority
- CN
- China
- Prior art keywords
- epsiv
- alpha
- partiald
- plasma density
- plasma
- 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.)
- Expired - Fee Related
Links
Landscapes
- Plasma Technology (AREA)
Abstract
常压下低温等离子体密度参数的诊别方法,属于等离子体领域,本发明为解决现有等离子体诊断方法在其测量过程中精确度低、效率差、对环境要求较高、使用寿命短、成本相对昂贵等一系列问题。本发明方法包括以下步骤:步骤一、用电压互感器检测放电通道内气体压降Ug(t);步骤二、用电流互感器检测为传导电流jc(x,t);步骤三、根据步骤一获取的放电通道内气体压降Ug(t)和传导电流jc(x,t)获取总的放电电流jT;步骤四、根据步骤三获取的总的放电电流jT获取等离子体密度n(t);步骤五、获取离散化等离子体密度nm+1(t),以实现对等离子体密度参数的诊别。
Description
技术领域
本发明涉及常压下低温等离子体密度参数的诊别方法,属于等离子体领域。
背景技术
近年来,基于介质阻挡放电(DBD)原理产生的低温等离子体在材料、微电子、化工、机械及环境保护等众多学科领域中有着广泛应用,且已形成一个崭新的工业——等离子体工业。例如在材料学科中,通常采用等离子体物理气相沉积(PVD)和增强化学气相沉积(PECVD)及等离子体源离子注入(PSII)等技术来制备一些发光、光电、微电子、耐蚀耐磨及超硬等新型多功能薄膜材料。利用等离子体进行材料表面改性,改善材料表面的可湿性、吸水性、可染性、粘着性和导电性等,更是常见的工业技术;在化工工业中,采用等离子聚合技术,可以印刷和制备出高分子薄膜材料;在微电子工业中的应用更加引人瞩目,目前微电子工业的全球销售额已达几千亿美元,其中三分之一以上微电子器件设备是采用低温等离子体技术生产的;在超大规模集成电路的生产工艺中,等离子体刻蚀技术可以实现高刻速率、高纵横比、高选择比、微观不均匀性小和低能量操作的刻蚀过程;在沉积无缺陷、附着力大的微电子薄膜材料以及微小器件晶片的清洗方面也显示出了巨大的优势。可以说,低温等离子体已与现代高新科技的发展紧密联系在一起,因而对低温等离子体及其应用诊断技术的研究也变得尤为重要。
由于等离子体的现象和行为复杂多变,使得人们常常无法对所需要求得的量进行直接测量,例如密度这个重要参数,人们常常只能通过等离子体在特定条件下所表现出来的某些现象来间接地求得这些参量,因此等离子体诊断是实现等离子体工业应用必不可少的方法手段。目前,常用的诊断技术主要包括探针诊断法、微波干涉诊断法、激光差拍法、光谱质谱诊断法等。当等离子体中包含波动、振荡和波时,探针方法的应用非常困难,有时甚至不能应用探针方法,且探针表面的杂技可能污染等离子体,这将使等离子体的I-V特性曲线发生变形,严重影响测量结果;微波法对等离子体的空间响应较差,动态范围较小;激光法为了能够测出散射信号而且有较小的统计误差必须采用大功率巨脉冲激光器作光源,采用灵敏度高、信噪比大、时间响应快的光探测器作接受器,不仅操作不便,且成本大大增加;光谱质谱极其复杂,较难精确解释,甚至在用作等离子体刻蚀工艺终点探测的分子谱线,有时并不清楚其来源,且其光学窗口上薄膜沉积或刻蚀能够大大改变或减弱光谱的信号。
发明内容
本发明目的是为了解决现有等离子体诊断方法在其测量过程中精确度低、效率差、对环境要求较高、使用寿命短、成本相对昂贵等一系列问题,提供了一种常压下低温等离子体密度参数的诊别方法。
本发明所述的常压下低温等离子体密度参数的诊别方法,该方法包括以下步骤:
步骤一、用电压互感器检测放电通道内气体压降Ug(t)
步骤二、用电流互感器检测为传导电流jc(x,t),x=0~dg,
式中:dg为放电通道长度,x为放电通道中某一位置的坐标,
步骤三、根据步骤一获取的放电通道内气体压降Ug(t)和传导电流jc(x,t),利用公式
获取总的放电电流jT,
式中:ε0为真空介电常数,εB为等离子体的相对介电常数,dB为所电极表面所覆盖介质厚度,
步骤四、根据步骤三获取的总的放电电流jT,利用公式
获取等离子体密度n(t),
步骤五、获取离散化等离子体密度nm+1(t),以实现对等离子体密度参数的诊别。
本发明的优点:本发明可根据外施参量的不同,得到不同条件下的密度结果,可将其应用于下述工程当中:
1.等离子体气象沉积技术,其辅助等离子体密度约为106~108cm-3;
2.等离子体刻蚀工艺,其等离子体密度约为108cm-3;
3.近空间再入飞行器地面模拟实验,其等离子体密度约为~1010cm-3;
4.辉光放电等离子体的研究工作,其等离子体密度为<1014cm-3;
对密度的测量实时性好,效率高,对环境要求低,使用寿命长,成本低。
附图说明
图1是本发明所述常压下低温等离子体密度参数的诊别方法的流程图;
图2是利用本发明方法计算得出的不同时刻放电空间电子数密度分布图;
图3是电子、离子及电场随时间分布图,其中虚线为电子密度;点划线为离子密度;实线为气隙间电场值;
图4是等离子体中电子密度随时间分布图;
图5是等离子体中负离子密度随时间分布图;
图6是等离子体中正离子密度随时间分布图;
图7是利用本发明方法计算结果与实验结果对比。其中“—·—”为计算结果;“O”为实验数据。
具体实施方式
具体实施方式一:下面结合图1至图7说明本实施方式,本实施方式所述常压下低温等离子体密度参数的诊别方法,该方法包括以下步骤:
步骤一、用电压互感器检测放电通道内气体压降Ug(t)
步骤二、用电流互感器检测为传导电流jc(x,t),x=0~dg,
式中:dg为放电通道长度,x为放电通道中某一位置的坐标,
步骤三、根据步骤一获取的放电通道内气体压降Ug(t)和传导电流jc(x,t),利用公式
获取总的放电电流jT,
式中:ε0为真空介电常数,εB为等离子体的相对介电常数,dB为所电极表面所覆盖介质厚度,
步骤四、根据步骤三获取的总的放电电流jT,利用公式
获取等离子体密度n(t),
步骤五、获取离散化等离子体密度nm+1(t),以实现对等离子体密度参数的诊别。
具体实施方式二:下面结合图2说明本实施方式,本实施方式对实施方式一作进一步说明,步骤五中获取离散化等离子体密度nm+1(t)的过程为:
以下离散过程中离散成k份,m=1,2,...k,k为正整数,
步骤51、电子、离子连续性方程差分方法离散化:
其中:S为描述带电粒子产生和复合的源项,
且通过S=αμe|E(x,t)|ne(x,t)-βne(x,t)ni(x,t)获取,
其中:α为汤森电离系数,β为电子与离子间的复合系数;
μe为电子迁移率,
μi为离子迁移率,
E(x,t)为电场强度;
ne(x,t)为电子密度,ni(x,t)为离子密度;
步骤52、带电粒子及亚稳态分布的输运用动量方程表示为:
其中:
De为电子扩散系数;Di为离子扩散系数;
则其离散化得到:
式中:
式中:εr为相对介电常数,且εr=ε0εB,
步骤53、根据步骤51和52获取离散化等离子体密度nm+1(t):
对所需诊断的等离子体采用一维流体力学模型,其带电粒子及亚稳态分布的输运由连续性方程及动量方程描述:
方程组(1)为连续性方程,方程组(2)为动量方程,其中ne、ni分别为电子、离子密度;且n=ne-ni,用je、ji分别为电子流、离子流密度;μe和De分别为电子迁移率和电子扩散系数;μi和Di分别为离子迁移率和离子扩散系数;E为电场强度;S为描述带电粒子产生和复合的源项,其具体表达式为:
S=αμe|E(x,t)|ne(x,t)-βne(x,t)ni(x,t) (3)
其中α为汤森电离系数;β为电子与离子间的复合系数。将方程组(1)中的两式相减可得:
其中ρ(x,t)为放电空间中的电荷密度;jc为放电过程中的传导电流。同时利用泊松方程来模拟放电通道中的电场强度:
其中ε0为真空介电常数。对上述泊松方程两端对时间求导,有:
同时将(4)式代入(6)式得到电流平衡方程:
即:
其中为放电过程中的位移电流;jT为放电过程中总放电电流。将(8)式两端对整个放电通道进行积分可得:
其中dg为放电通道长度;Ug为放电通道内气体压降。同时得到介质上的电流平衡方程为:
将(9)式与(10)式相加,可得:
从而可以得到总的放电电流jT,其表达式为:
εr为、dB为所电极表面所覆盖介质厚度,
由于等离子体中电子密度与离子密度为:
n=e[ji-(je+γji)]=ejT (13)
将公式(12)代入公式(13),得到等离子体密度随时间演化的表达式为:
则ne(t)和ni(t)的表达式按公式(14)类似,只替换相应参数即可;
由于上述方程很难得到解析解,因此对上述方程的求解是本发明的一大关键技术。将整个一维放电空间剖分为整结点及半结点,整结点上为电子、离子密度,而在半结点上为电场值,电子电流密度及离子电流密度。
将公式(15)和(16)进行离散,获取电子、离子连续性方程差分方法离散化:
电场的计算由电流平衡方程的离散化得到:
根据公式(17)至公式(21)获取离散化等离子体密度nm+1(t):
通过上述计算可在计算机上得到等离子体密度参量的实时控制。本文明可根据外施参量的不同,得到不同条件下的密度结果,可将其应用于下述工程当中:
1.等离子体气象沉积技术,其辅助等离子体密度约为106~108cm-3;
2.等离子体刻蚀工艺,其等离子体密度约为108cm-3;
3.近空间再入飞行器地面模拟实验,其等离子体密度约为~1010cm-3;
4.辉光放电等离子体的研究工作,其等离子体密度为<1014cm-3;
根据公式(22)获取了离散的等离子体的密度,使得我们可以实时监测当前参量是否适合于具体的工业应用。
Claims (2)
1.常压下低温等离子体密度参数的诊别方法,其特征在于,该方法包括以下步骤:
步骤一、用电压互感器检测放电通道内气体压降Ug(t)
步骤二、用电流互感器检测放电过程中的传导电流jc(x,t),x=0~dg,
式中:dg为放电通道长度,x为放电通道中某一位置的坐标,
步骤三、根据步骤一获取的放电通道内气体压降Ug(t)和传导电流jc(x,t),利用公式
获取总的放电电流jT,
式中:ε0为真空介电常数,εB为等离子体的相对介电常数,dB为所电极表面所覆盖介质厚度,
步骤四、根据步骤三获取的总的放电电流jT,利用公式
获取等离子体密度n(t),
步骤五、获取离散化等离子体密度nm+1(t),以实现对等离子体密度参数的诊别。
2.根据权利要求1所述常压下低温等离子体密度参数的诊别方法,其特征在于,步骤五中获取离散化等离子体密度nm+1(t)的过程为:
以下离散过程中离散成k份,m=1,2,…k,k为正整数,
步骤51、电子、离子连续性方程差分方法离散化:
其中:S为描述带电粒子产生和复合的源项,
且通过S=αμe|E(x,t)|ne(x,t)-βne(x,t)ni(x,t)获取,
其中:α为汤森电离系数,β为电子与离子间的复合系数;
μe为电子迁移率,
μi为离子迁移率,
E(x,t)为电场强度;
ne(x,t)为电子密度,ni(x,t)为离子密度;
步骤52、带电粒子及亚稳态分布的输运用动量方程表示为:
其中:
De为电子扩散系数;Di为离子扩散系数;
则其离散化得到:
式中:
式中:εr为相对介电常数,且εr=ε0εB,
步骤53、根据步骤51和52获取离散化等离子体密度nm+1(t):
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310010841.6A CN103048522B (zh) | 2013-01-11 | 2013-01-11 | 常压下低温等离子体密度参数的诊别方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310010841.6A CN103048522B (zh) | 2013-01-11 | 2013-01-11 | 常压下低温等离子体密度参数的诊别方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103048522A CN103048522A (zh) | 2013-04-17 |
CN103048522B true CN103048522B (zh) | 2015-03-11 |
Family
ID=48061234
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310010841.6A Expired - Fee Related CN103048522B (zh) | 2013-01-11 | 2013-01-11 | 常压下低温等离子体密度参数的诊别方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103048522B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111800932B (zh) * | 2020-08-07 | 2021-07-06 | 大连理工大学 | 一种等离子体放电过程模拟方法及系统 |
CN112949045B (zh) * | 2021-02-06 | 2022-10-14 | 四川大学 | 一种稳态微波等离子体的测算方法 |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5404984B2 (ja) * | 2003-04-24 | 2014-02-05 | 東京エレクトロン株式会社 | プラズマモニタリング方法、プラズマモニタリング装置及びプラズマ処理装置 |
JP2006186147A (ja) * | 2004-12-28 | 2006-07-13 | Konica Minolta Holdings Inc | シミュレーション装置及びシミュレーションプログラム |
JP2007011475A (ja) * | 2005-06-28 | 2007-01-18 | Nippon Hoso Kyokai <Nhk> | Pdp放電特性解析装置およびpdp放電特性解析プログラム |
JP4756272B2 (ja) * | 2006-07-28 | 2011-08-24 | 独立行政法人産業技術総合研究所 | マイクロ波プラズマ解析プログラム |
JP2008052530A (ja) * | 2006-08-25 | 2008-03-06 | Hamamatsu Metrix Kk | プラズマ粒子シミュレーション計算方法 |
CN201096521Y (zh) * | 2007-11-06 | 2008-08-06 | 南京理工大学 | 非接触式等离子体温度和电子密度测量装置 |
JP4982803B2 (ja) * | 2008-01-28 | 2012-07-25 | 独立行政法人産業技術総合研究所 | マイクロ波プラズマ解析プログラム |
CN102156001A (zh) * | 2011-03-17 | 2011-08-17 | 大连理工大学 | 一种射频放电等离子体自偏置探针诊断方法 |
-
2013
- 2013-01-11 CN CN201310010841.6A patent/CN103048522B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN103048522A (zh) | 2013-04-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Marjamäki et al. | Performance evaluation of the electrical low-pressure impactor (ELPI) | |
CN106354925B (zh) | 一种基于电晕放电空间电势分布的模拟方法 | |
CN103267723B (zh) | 基于场指纹法的金属管道、容器坑蚀检测方法 | |
CN101361176A (zh) | 等离子体诊断设备及方法 | |
CN103558531A (zh) | 基于pea空间电荷测量判断绝缘老化的方法 | |
CN103278540B (zh) | 一种对高压绝缘子等值盐密测量的方法 | |
CN103048522B (zh) | 常压下低温等离子体密度参数的诊别方法 | |
KR100902946B1 (ko) | 소프트 엑스레이 광이온화 하전기 | |
CN110688751A (zh) | 铂掺杂改性石墨炔传感器检测sf6的仿真方法 | |
Bui et al. | Pressure sensor based on bipolar discharge corona configuration | |
Liu et al. | Size-resolved chemical composition analysis of ions produced by a commercial soft X-ray aerosol neutralizer | |
Döring et al. | A novel approach for road surface wetness detection with planar capacitive sensors | |
CN103635004A (zh) | 一种等离子体中离子种类与数量密度分布的测量方法 | |
Wilkinson et al. | Evaluation of dielectric-barrier-discharge actuator substrate materials | |
Fanara et al. | Electrical characterization of atmospheric pressure arc plasmas: An overview | |
CN101408514B (zh) | 基于气体放电光谱分析的气体传感器及其检测气体的方法 | |
CN106290546A (zh) | 离子迁移谱仪 | |
CN104122492B (zh) | 一种预测半导体器件10年寿命对应的工作电压的方法 | |
Yuan et al. | Detection of decomposition products of SF6/air gas mixture by electron attachment mass spectrometry | |
CN105738037A (zh) | 一种等离子反应腔体的渗漏检测方法 | |
CN105466976B (zh) | 氢气传感器芯体用介质材料、氢气传感器芯体及其制备方法与应用 | |
Hou et al. | Development of a suitcase time‐of‐flight mass spectrometer for in situ fault diagnosis of SF6‐insulated switchgear by detection of decomposition products | |
CN106768034B (zh) | 一种测量真空罐体内沉积污染的方法 | |
CN104062350A (zh) | 一种用于检测铜金属蒸气浓度的多电极微型传感器 | |
CN103257284A (zh) | 一种基于激光热脉冲的有机聚合物薄膜空间电荷测量方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | 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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150311 Termination date: 20170111 |