CN107992684A - 一种时变等离子体等效分层介质模型建模方法 - Google Patents
一种时变等离子体等效分层介质模型建模方法 Download PDFInfo
- Publication number
- CN107992684A CN107992684A CN201711269808.XA CN201711269808A CN107992684A CN 107992684 A CN107992684 A CN 107992684A CN 201711269808 A CN201711269808 A CN 201711269808A CN 107992684 A CN107992684 A CN 107992684A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- plasma
- msup
- time
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Aviation & Aerospace Engineering (AREA)
- Automation & Control Theory (AREA)
- Computer Graphics (AREA)
- Software Systems (AREA)
- Plasma Technology (AREA)
Abstract
本发明公开了一种时变等离子体等效分层介质模型建模方法,包含以下过程:采用等值面提取算法对高超声速目标绕流流场数据进行分析,获取等离子体等电子数密度廓面信息。利用所述等离子体等电子数密度廓面将临近空间超高声速目标绕流流场区域划分为边界层与普通层,并对各层等离子体所对应的特征频率、介电常数进行求解,建立等离子体的等效分层介质模型。本发明具有实现对等离子体动态时变特性的描述,简化了后续电磁特性求解运算,极大的扩展了适用范围的优点。
Description
技术领域
本发明涉及电磁散射建模领域,特别涉及一种时变等离子体等效分层介质模型建模方法。
背景技术
当超高速目标运行于稀薄大气层内时,由于强烈的激波压缩和粘性摩擦作用,其周围的温度迅速升高,使空气发生离解和电离,形成包覆目标的激波等离子体。再入目标包覆的等离子体是不均匀的弯曲等离子体层,同时,其边界层存在湍流现象,由此引起的等离子鞘套参数具有时变特性。电磁波在这种时变非线性等离子体中传播时,会呈现出频率漂移、非线性幅度衰减、电磁波断裂等非线性现象。因此,针对激波等离子体包覆目标的散射传输建模需求,建立逼真的时变等离子体等效分层介质模型,具有重要意义。目前,对于均匀等离子体等效介质模型的研究已日趋完善,而对于稀薄大气层内超高速目标激发的时变非均匀等离子体等效分层介质模型建模方法则有待发展。
对现有技术进行了国内外数据库的检索,在国家发明专利“一种电磁环境体数据等值面提取方法”(申请号:201310258955.2)中提出了一种高效的电磁环境等值面提取方法。该方法通过将电磁环境数据按照八叉树的结构进行存储,通过遍历八叉树,找出含有等值面的体素,实现对等值面的提取,但是无法对以不规则点云形式存储的电磁环境体数据进行等值面提取。
在国家发明专利“高超声速飞行器等离子体鞘套与电磁波相互作用预测方法”(申请号:201611033297.7)中提出了一种基于等离子体鞘套分层模型的等离子体电磁传播预测方法。该方法实质是通过平行分层结构模拟非均匀等离子体的空间分布,并通过传输矩阵方法实现等离子体透射的计算,但是平行分层结构无法实现对等离子体分布的逼真模拟,将引入很大误差。
在国家发明专利“一种基于分层结构的动态等离子鞘套电子密度建模方法”(申请号:201710052765.3)中公开了一种基于分层结构的动态等离子鞘套电子密度建模方法。该方法通过时空电子密度矩阵来实现对动态等离子体鞘套参数的存储,对内存需求较大,后续电磁传输计算复杂。
2014年《IEEE Transactions on Plasma Science》期刊中公开的文献《Effectsof Reentry Plasma Sheath on the Polarization Properties of Obliquely IncidentEM waves》介绍了一种非均匀等离子体等效介质模型建模方法。该方法可将再入飞行器等离子体鞘套等效为多层平行分层的介质,但是在等离子体分布情况比较复杂时,平行分层处理会引入很大误差。
2015年《真空科学与技术学报》期刊中公开的论文《大面积等离子体片分层现象的实验研究》介绍了一种用等离子体片分层模型等效非均匀等离子体分布的方法。该方法实质仍是采用多层平行分层模型来实现对非均匀等离子体模拟,将引入很大误差。
针对稀薄大气层内目标高超声速运动所激发的等离子体鞘套,为了对其电磁特性进行建模需建立其等效介质模型,目前多采用多层平行分层介质来实现对非均匀分布的等离子体鞘套的模拟。该方法虽然简化了后期电磁特性仿真的难度,但是采用的平行分层模型属于一种近似,与真实的等离子体鞘套相比存在很大误差。因此,现有技术方案的适用性受到很大限制。
发明内容
本发明的目的是提供一种时变等离子体等效分层介质模型建模方法,通过Matching Cubes算法对高超声速目标等离子体流场数据的等电子数密度廓面进行提取,利用分层的有耗介质体来等效时变非均匀等离子体鞘套对电磁波的影响,采用等效分层介质对非均匀等离子体的分布特性进行描述;通过将等离子体参数表征为时间的指数级数形式,实现对等离子体动态时变特性的描述,简化了后续电磁特性求解运算,扩展适用范围的目的。
为了实现以上目的,本发明通过以下技术方案实现:
一种时变等离子体等效分层介质模型建模方法,包含以下过程:采用等值面提取算法对高超声速目标绕流流场数据进行分析,获取等离子体等电子数密度廓面信息。利用所述等离子体等电子数密度廓面将临近空间超高声速目标绕流流场区域划分为边界层与普通层,并对各层等离子体所对应的特征频率、介电常数进行求解,建立等离子体的等效分层介质模型。
优选地,所述等值面提取算法为Marching Cubes算法,通过逐个处理数据场中的立方体,分类出与等值面相交的立方体,采用插值方法计算出等值面与立方体边的交点;根据立方体每一顶点与等值面的相对位置,将等值面与立方体边的交点连接生成等值面。通过对点云信息数据中所有立方体元素的遍历,获得所述等值面在每个立方体元素中的分布情况,连接各三角面元获得所述等电子数密度廓面信息。
优选地,所述点云信息数据包括:压强、碰撞频率、温度、电子数密度、各类气体组分密度在空间的分布,所述点云信息数据为通过采用热化学非平衡热环境的N-S/DSMC耦合算法对稀薄大气层内高超声速目标绕流流场的建模得到。
优选地,临近空间超高声速目标绕流流场中的普通层区域的等离子体相对介电常数为:
式中:ω为电磁波频率;ωp为等离子振荡角频率,qe为电子电荷,ε0为真空介电常数,me为电子质量,Ne为电子数密度(cm-3);ν为电子与中性粒子和其它粒子的碰撞频率。
优选地,在纯空气的情况下,所述碰撞频率ν与温度、压力的经验关系式为:
ν=5.82×1012T-1/2P
式中,T为温度,2000K<T<6000K,P为压强;
对带有烧蚀组分的混合气体,碰撞频率ν为:
式中,na为所有粒子数密度,Xi为混合气体中每个组分的质量百分率,Qi为每一组分电子的矩传输截面,k为波耳兹曼常数;
k=β-iα
式中,折射率衰减率c为光速,ω为电磁波频率,k0为自由空间波数。
优选地,对于位于绕流流场边界层的时变等离子体,等离子体的电子数密度Ne表示为指数形式级数,其形式如下
式中,cn为级数的系数,s为与等离子体振荡频率相关的系数;
边界层的时变等离子体的特征频率为
式中,为常数,其中e、me分别为电子电荷和电子质量;ε0为真空介电常数;
根据时变等离子体的特征频率,边界层的复介电常数为
式中,b为与时间无关的复数,Ne(t)为边界层的电子数密度,ω为电磁波频率,ν为电子与中性粒子和其它粒子的碰撞频率。
本发明与现有技术相比具有以下优点:
本发明通过Matching Cubes算法对高超声速目标等离子体流场数据的等电子数密度廓面进行提取,利用分层的有耗介质体来等效时变非均匀等离子体鞘套对电磁波的影响,采用等效分层介质对非均匀等离子体的分布特性进行描述;通过将等离子体参数表征为时间的指数级数形式,实现对等离子体动态时变特性的描述,简化了后续电磁特性求解运算,极大的扩展了适用范围。
附图说明
图1为根据轴对称性进行截取后的钝锥体及绕流流场区域示意图,图中最内层为钝锥体模型,其球头半径Rn=0.08m,钝锥总长L=0.3m,半锥角θ=8°,外层绕流流场被划分为14个分区;
图2为飞行速度10Ma,高度65km的钝锥体目标的等效分层介质模型示意图,图中内外两层曲面对应的电子数密度分别为1E12cm-3和6E9cm-3;
图3为图2中所示外层等离子体的等效介电常数的实部随频率的变化曲线;
图4为图2中所示外层等离子体的等效介电常数的虚部随频率的变化曲线;
图5为本发明一种时变等离子体等效分层介质模型建模方法的流程图。
具体实施方式
以下结合附图,通过详细说明一个较佳的具体实施例,对本发明做进一步阐述。
如图5所示,本发明一种时变等离子体等效分层介质模型建模方法,包含以下过程:
步骤S1、等离子体等电子数密度廓面提取。利用等值面提取算法对高超声速目标绕流流场数据进行分析,建立等离子体等电子数密度廓面模型。
所述步骤S1进一步包含以下过程:
高超声速目标及其周围等离子体流场构成目标头身部区域,其电磁散射模型相当于一个有耗非均匀介质体对电磁波的散射。
通过采用热化学非平衡热环境的N-S/DSMC耦合算法实现稀薄大气层内高超声速目标绕流流场的建模,获取包含压强、碰撞频率、温度、电子数密度、各类气体组分密度等参数在空间分布的点云信息数据。
提取高超声速目标绕流流场体数据中电子数密度分布的最大值与最小值,将其N等分,对于各分界值调用等值面提取的Marching Cubes(MC)算法,获取等电子数密度廓面信息。利用等电子数密度廓面将临近空间超高速目标绕流流场区域进行划分,从而获得等离子体的等效分层介质模型。
通过高超声速目标流场建模求解获得的临近空间目标等离子体点云信息数据是一种内部含有信息的三维体数据,并且没有三角面片或者边缘等几何信息。
MC算法是一种体素单元内等值面抽取算法。这一算法所处理的一般是三维正交的数据场,可以表示为
Fi,j,k=F(xi,yj,zk)(i=1,...,Nx,j=1,...,Ny,k=1,...,Nz) (1)
式中,Fi,j,k是在点(i,j,k)处的取值,F(xi,yj,zk)是一个以离散坐标为自变量的函数。
MC算法中的立方体体素是一逻辑上的立方体,由相邻层上的各四个像素组成立方体上的八个顶点。算法以扫描线方式逐个处理数据场中每一个立方体体素,求出每一体素内包含的等值面,由此生成整个数据场的等值面。
MC算法的基本原理是逐个处理数据场中的立方体,分类出与等值面相交的立方体,采用插值方法计算出等值面与立方体边的交点。根据立方体每一顶点与等值面的相对位置,将等值面与立方体边的交点按一定的方式连接生成等值面,作为等值面在该立方体内的一个逼近表示。因此,MC算法中的两个主要计算是:关于体素中由三角片逼近的等值面的计算,以及关于三角片各顶点法向量的计算。
MC算法的基本假设是沿着立方体的边,数据场呈现连续线性变化,也就是讲,如一条边的两个顶点分别大于小于等值面的值,则在该边上必有也仅有一点是这条边与等值面的交点。确定立方体体素中等值面的分布是该算法的基础。
交点的求法是采用线性插值计算法,如式(2)所示。
X=i+(V-Fi,j,k)/(Fi+1,j,k-Fi,j,k) (2)
式中,Fi,j,k为点(i,j,k)处的取值;Fi+1,j,k为点(i+1,j,k)处的取值;V是当前求解等值面的值。
由于函数F(x,y,z)的梯度垂直于等值面,因此数据场Fi,j,k在等值面点P的梯度可以作为P点的法向量。由于是三维正交数据场,可以采用中心差分计算出立方体顶点的梯度,再用线性插值计算P点的梯度值,计算公式如下
式中,(a,b,c)是立方体体素的大小,函数f(x,y,z)是在点(x,y,z)处的取值,(gx,gy,gz)是点(x,y,z)处的梯度。
该算法执行时,每次扫描其中两层图像,构造这两层之间的立方体体素。首先对立方体顶点进行分类,由分类顶点建立该立方体在检索分类表中的索引,索引分类表中对应的等值面分布模式。由线性插值计算出三角片顶点的位置和梯度值,将立方体上的交点按等值面片连接模式连接成三角片或多边形,最后使用图形API绘制出来。其中根据立方体中八个顶点与等值面的包含关系,不同的索引号共有256种,排除对称性后,等值面在立方体素中的分布情况,可以概括为十五种情况。
通过对点云数据中所有立方体元素的遍历,获得等值面在每个立方体元素中的分布情况,连接这些三角面元即可获得最终等离子体数据中等电子数密度的廓面信息。
步骤S2、等离子体等效介质模型建模。对各层等离子体所对应的特征频率、介电常数等参数进行求解,建立等离子体的等效分层介质模型。
所述步骤S2进一步包含以下过程:等离子体相对介电常数是一复数,即
式中:ω为电磁波频率;ωp为等离子振荡角频率,qe为电子电荷,ε0为真空介电常数,me为电子质量,Ne为电子数密度(cm-3);ν为电子与中性粒子和其它粒子的碰撞频率。纯空气中碰撞频率与温度、压力的经验关系式是ν=5.82×1012T-1/2P,式中T为温度,2000K<T<6000K,P为压强。
对带有烧蚀组分的混合气体,碰撞频率可由气体运动论得到na为所有粒子数密度,Xi为混合气体中每个组分的质量百分率,Qi为每一组分电子的矩传输截面,k为波耳兹曼常数。
等离子中电磁波的波矢是一复数,即k=β-iα,折射率衰减率c为光速,ω为电磁波频率,k0为自由空间波数。
美国相关机构对RAM-3飞行器的风洞试验研究指出,在同一固定高度,电子密度随厚度逐渐变化,在0mm~8mm厚度处,电子密度梯度较大,通常将其视作边界层。等离子体鞘套的边界层存在湍流现象,湍流变化频率最高可达数十千赫兹,并且频率越大抖动值越小。
将等离子体绕流流场区域分为边界层和普通层。等离子绕流流场在边界层湍流脉动会引起的等离子体鞘套的时变特性,对于位于绕流流场边界层的时变等离子体,等离子体的电子数密度Ne是一个随时间变化的函数,变化频率可达到数十千赫兹,可以表示为指数形式级数,其形式如下
式中,cn为级数的系数,s为与等离子体振荡频率相关的系数。
边界层的介电常数等特性也相应具有时变特性。
针对等离子体鞘套边界层的时变现象,时变等离子体的特征频率可以表示为
式中,为常数,其中e、me分别为电子电荷和电子质量;ε0为真空介电常数。
根据时变等离子体的特征频率,可求得其对应的复介电常数
式中,b为与时间无关的复数,Ne(t)为边界层的电子数密度,ω为电磁波频率,ν为电子与中性粒子和其它粒子的碰撞频率。
采用上述方法可基于临近空间超高声速目标绕流流场信息数据获取等离子体的等效分层介质模型。
本发明的一个实施例为以稀薄大气层内高超声速飞行钝锥体目标绕流流场的等效分层介质模型建模,包含以下过程:
图1为根据轴对称性进行截取后的钝锥体及绕流流场区域示意图,图中最内层为钝锥体模型,其球头半径Rn=0.08m,钝锥总长L=0.3m,半锥角θ=8°,外层绕流流场被划分为14个分区1~14;首先,利用上文所述的MC算法对该外层绕流流场的14个分区1~14的电子数密度分布数据进行分析,提取等电子数密度廓面模型,该等电子数密度廓面模型如图2所示,当钝锥体模型在高度65km,以速度为10Ma进行飞行时,钝锥体模型的等效分层介质模型的内外两层曲面对应的电子数密度分别为1E12cm-3和6E9cm-3。之后对各层等离子体所对应的等效介质参数进行求解,并通过在边界层引入随时间呈指数级数形式变化的等离子体参数,实现对时变等离子体的描述。结合图3与图4所示,外层等离子体的等效介电常数的实部与虚部随频率的变化曲线,随着频率的升高等离子体的等效介电常数实部逐渐向1逼近,虚部向0逼近,等离子体的等效介电特性逐渐接近空气,对电磁波传输的影响逐渐减弱。
尽管本发明的内容已经通过上述优选实施例作了详细介绍,但应当认识到上述的描述不应被认为是对本发明的限制。在本领域技术人员阅读了上述内容后,对于本发明的多种修改和替代都将是显而易见的。因此,本发明的保护范围应由所附的权利要求来限定。
Claims (6)
1.一种时变等离子体等效分层介质模型建模方法,其特征在于,包含以下过程:
采用等值面提取算法对高超声速目标绕流流场数据进行分析,获取等离子体等电子数密度廓面信息;
利用所述等离子体等电子数密度廓面将临近空间超高声速目标绕流流场区域划分为边界层与普通层,并对各层等离子体所对应的特征频率、介电常数进行求解,建立等离子体的等效分层介质模型。
2.如权利要求1所述的时变等离子体等效分层介质模型建模方法,其特征在于,所述等值面提取算法为Marching Cubes算法,
通过逐个处理数据场中的立方体,分类出与等值面相交的立方体,采用插值方法计算出等值面与立方体边的交点;根据立方体每一顶点与等值面的相对位置,将等值面与立方体边的交点连接生成等值面;
通过对点云信息数据中所有立方体元素的遍历,获得所述等值面在每个立方体元素中的分布情况,连接各三角面元获得所述等电子数密度廓面信息。
3.如权利要求2所述的时变等离子体等效分层介质模型建模方法,其特征在于,所述点云信息数据包括:压强、碰撞频率、温度、电子数密度、各类气体组分密度在空间的分布,所述点云信息数据为通过采用热化学非平衡热环境的N-S/DSMC耦合算法对稀薄大气层内高超声速目标绕流流场的建模得到。
4.如权利要求1所述的时变等离子体等效分层介质模型建模方法,其特征在于,临近空间超高声速目标绕流流场中的普通层区域的等离子体相对介电常数为:
<mrow>
<msub>
<mi>&epsiv;</mi>
<mi>r</mi>
</msub>
<mo>=</mo>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<msubsup>
<mi>&omega;</mi>
<mi>p</mi>
<mn>2</mn>
</msubsup>
<mrow>
<msup>
<mi>&omega;</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mi>v</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>-</mo>
<mi>i</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mi>v</mi>
<mi>&omega;</mi>
</mfrac>
<mfrac>
<msubsup>
<mi>&omega;</mi>
<mi>p</mi>
<mn>2</mn>
</msubsup>
<mrow>
<msup>
<mi>&omega;</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mi>v</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
</mrow>
式中:ω为电磁波频率;ωp为等离子振荡角频率,qe为电子电荷,ε0为真空介电常数,me为电子质量,Ne为电子数密度(cm-3);ν为电子与中性粒子和其它粒子的碰撞频率。
5.如权利要求4所述的时变等离子体等效分层介质模型建模方法,其特征在于,在纯空气的情况下,所述碰撞频率ν与温度、压力的经验关系式为:
ν=5.82×1012T-1/2P
式中,T为温度,2000K<T<6000K,P为压强;
对带有烧蚀组分的混合气体,碰撞频率ν为:
<mrow>
<mi>v</mi>
<mo>=</mo>
<msub>
<mi>n</mi>
<mi>a</mi>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<mn>3</mn>
<mi>k</mi>
<mi>T</mi>
</mrow>
<msub>
<mi>m</mi>
<mi>e</mi>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<mrow>
<mn>1</mn>
<mo>/</mo>
<mn>2</mn>
</mrow>
</msup>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msub>
<mi>X</mi>
<mi>i</mi>
</msub>
<msub>
<mi>Q</mi>
<mi>i</mi>
</msub>
</mrow>
式中,na为所有粒子数密度,Xi为混合气体中每个组分的质量百分率,Qi为每一组分电子的矩传输截面,k为波耳兹曼常数;
k=β-iα
式中,折射率衰减率c为光速,ω为电磁波频率,k0为自由空间波数。
6.如权利要求1所述的时变等离子体等效分层介质模型建模方法,其特征在于,对于位于绕流流场边界层的时变等离子体,等离子体的电子数密度Ne表示为指数形式级数,其形式如下
<mrow>
<msub>
<mi>N</mi>
<mi>e</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munder>
<mo>&Sigma;</mo>
<mi>n</mi>
</munder>
<msub>
<mi>c</mi>
<mi>n</mi>
</msub>
<msup>
<mi>e</mi>
<mrow>
<mi>j</mi>
<mi>n</mi>
<mi>s</mi>
<mi>t</mi>
</mrow>
</msup>
</mrow>
式中,cn为级数的系数,s为与等离子体振荡频率相关的系数;
边界层的时变等离子体的特征频率为
<mrow>
<msub>
<mi>&omega;</mi>
<mi>p</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msqrt>
<mrow>
<mo>(</mo>
<msub>
<mi>N</mi>
<mi>e</mi>
</msub>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<msup>
<mi>e</mi>
<mn>2</mn>
</msup>
<mo>)</mo>
<mo>/</mo>
<mo>(</mo>
<msub>
<mi>m</mi>
<mi>e</mi>
</msub>
<msub>
<mi>&epsiv;</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
</msqrt>
<mo>&ap;</mo>
<mi>a</mi>
<msqrt>
<mrow>
<msub>
<mi>N</mi>
<mi>e</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
</msqrt>
</mrow>
式中,为常数,其中e、me分别为电子电荷和电子质量;ε0为真空介电常数;
根据时变等离子体的特征频率,边界层的复介电常数为
<mrow>
<msub>
<mover>
<mi>&epsiv;</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>r</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<mrow>
<msup>
<msub>
<mi>&omega;</mi>
<mi>p</mi>
</msub>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msup>
<mi>&omega;</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mi>&upsi;</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>+</mo>
<mi>i</mi>
<mfrac>
<mrow>
<msup>
<msub>
<mi>&omega;</mi>
<mi>p</mi>
</msub>
<mn>2</mn>
</msup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mrow>
<mo>(</mo>
<mi>&upsi;</mi>
<mo>/</mo>
<mi>&omega;</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msup>
<mi>&omega;</mi>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mi>&upsi;</mi>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>=</mo>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>bN</mi>
<mi>e</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
</mrow>
式中,b为与时间无关的复数,Ne(t)为边界层的电子数密度,ω为电磁波频率,ν为电子与中性粒子和其它粒子的碰撞频率。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711269808.XA CN107992684B (zh) | 2017-12-05 | 2017-12-05 | 一种时变等离子体等效分层介质模型建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711269808.XA CN107992684B (zh) | 2017-12-05 | 2017-12-05 | 一种时变等离子体等效分层介质模型建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107992684A true CN107992684A (zh) | 2018-05-04 |
CN107992684B CN107992684B (zh) | 2021-01-19 |
Family
ID=62035799
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711269808.XA Active CN107992684B (zh) | 2017-12-05 | 2017-12-05 | 一种时变等离子体等效分层介质模型建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107992684B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110308654A (zh) * | 2019-07-01 | 2019-10-08 | 西安电子科技大学 | 基于模型预测控制的电感耦合等离子体产生系统的整定方法 |
CN111182705A (zh) * | 2020-01-03 | 2020-05-19 | 西安电子科技大学 | 基于自动编码器的时变等离子体诊断方法及诊断系统 |
CN111368380A (zh) * | 2018-12-24 | 2020-07-03 | 中国空气动力研究与发展中心超高速空气动力研究所 | 一种用于n-s/dsmc耦合算法的区域边界优化方法 |
CN112153796A (zh) * | 2020-09-27 | 2020-12-29 | 西安电子科技大学 | 时变等离子体电子密度抖动频率的测量方法 |
CN112966401A (zh) * | 2021-05-12 | 2021-06-15 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种热化学非平衡多级气体模型自适应算法 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1933328A (zh) * | 2005-09-13 | 2007-03-21 | 精工爱普生株式会社 | 脉冲产生电路和使用该电路的设备及信息传送方法 |
CN102405404A (zh) * | 2009-02-02 | 2012-04-04 | 行星排放管理公司 | 监测温室气体通量的各系统的系统 |
CN103295269A (zh) * | 2013-06-26 | 2013-09-11 | 电子科技大学 | 一种电磁环境体数据等值面提取方法 |
CN104561139A (zh) * | 2013-10-18 | 2015-04-29 | 中国石油化工股份有限公司 | 一种提高微生物终细胞密度并缩短培养时间的方法 |
CN105260507A (zh) * | 2015-09-16 | 2016-01-20 | 西安电子科技大学 | 一种快速的动态等离子鞘套电波传播计算方法 |
CN105653747A (zh) * | 2014-11-14 | 2016-06-08 | 南京理工大学 | 超高速飞行器共形亚网格电磁散射特性分析的仿真方法 |
WO2016112464A1 (en) * | 2015-01-15 | 2016-07-21 | General Fusion Inc. | Apparatus and method for generating a vortex cavity in a rotating fluid |
CN106934095A (zh) * | 2017-01-24 | 2017-07-07 | 西安电子科技大学 | 一种基于分层结构的动态等离子鞘套电子密度建模方法 |
CN107367371A (zh) * | 2017-07-14 | 2017-11-21 | 北京机电工程研究所 | 非定常空化流场可视化试验数据处理方法 |
-
2017
- 2017-12-05 CN CN201711269808.XA patent/CN107992684B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1933328A (zh) * | 2005-09-13 | 2007-03-21 | 精工爱普生株式会社 | 脉冲产生电路和使用该电路的设备及信息传送方法 |
CN102405404A (zh) * | 2009-02-02 | 2012-04-04 | 行星排放管理公司 | 监测温室气体通量的各系统的系统 |
CN103295269A (zh) * | 2013-06-26 | 2013-09-11 | 电子科技大学 | 一种电磁环境体数据等值面提取方法 |
CN104561139A (zh) * | 2013-10-18 | 2015-04-29 | 中国石油化工股份有限公司 | 一种提高微生物终细胞密度并缩短培养时间的方法 |
CN105653747A (zh) * | 2014-11-14 | 2016-06-08 | 南京理工大学 | 超高速飞行器共形亚网格电磁散射特性分析的仿真方法 |
WO2016112464A1 (en) * | 2015-01-15 | 2016-07-21 | General Fusion Inc. | Apparatus and method for generating a vortex cavity in a rotating fluid |
CN105260507A (zh) * | 2015-09-16 | 2016-01-20 | 西安电子科技大学 | 一种快速的动态等离子鞘套电波传播计算方法 |
CN106934095A (zh) * | 2017-01-24 | 2017-07-07 | 西安电子科技大学 | 一种基于分层结构的动态等离子鞘套电子密度建模方法 |
CN107367371A (zh) * | 2017-07-14 | 2017-11-21 | 北京机电工程研究所 | 非定常空化流场可视化试验数据处理方法 |
Non-Patent Citations (2)
Title |
---|
HUIYUAN ZHANG ET AL.: "Bi-static Angular Glint Calculation on Complex Targets in Near-Regions via Multilevel Fast Multipole Algorithm", 《2016 PROGRESS IN ELECTROMAGNETIC RESEARCH SYMPOSIUM (PIERS)》 * |
张粲宇,等: "非平行分层等离子体鞘套电波传播的复射线方法", 《电波科学学报》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111368380A (zh) * | 2018-12-24 | 2020-07-03 | 中国空气动力研究与发展中心超高速空气动力研究所 | 一种用于n-s/dsmc耦合算法的区域边界优化方法 |
CN111368380B (zh) * | 2018-12-24 | 2022-07-26 | 中国空气动力研究与发展中心超高速空气动力研究所 | 一种用于n-s/dsmc耦合算法的区域边界优化方法 |
CN110308654A (zh) * | 2019-07-01 | 2019-10-08 | 西安电子科技大学 | 基于模型预测控制的电感耦合等离子体产生系统的整定方法 |
CN110308654B (zh) * | 2019-07-01 | 2020-06-05 | 西安电子科技大学 | 基于模型预测控制的电感耦合等离子体产生系统的整定方法 |
CN111182705A (zh) * | 2020-01-03 | 2020-05-19 | 西安电子科技大学 | 基于自动编码器的时变等离子体诊断方法及诊断系统 |
CN112153796A (zh) * | 2020-09-27 | 2020-12-29 | 西安电子科技大学 | 时变等离子体电子密度抖动频率的测量方法 |
CN112966401A (zh) * | 2021-05-12 | 2021-06-15 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种热化学非平衡多级气体模型自适应算法 |
CN112966401B (zh) * | 2021-05-12 | 2021-07-16 | 中国空气动力研究与发展中心计算空气动力研究所 | 一种热化学非平衡多级气体模型自适应方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107992684B (zh) | 2021-01-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107992684A (zh) | 一种时变等离子体等效分层介质模型建模方法 | |
CN107942309B (zh) | 一种稀薄大气层内超高速目标电磁散射快速计算方法 | |
CN102129523B (zh) | 基于mda和mlssm的分析复杂目标电磁散射的方法 | |
CN106096267B (zh) | 一种腔体电磁散射特性快速计算方法 | |
CN109765538B (zh) | 非均匀介质目标体的散射场确定方法 | |
Zheng et al. | An efficient method to reduce the numerical dispersion in the ADI-FDTD | |
CN103425864A (zh) | 应用于金属复杂非均匀媒质混合目标的电磁散射分析方法 | |
CN106446470A (zh) | 一种高效并行的非均匀介质频域有限差分方法 | |
Lu et al. | A study on zoning coating method of absorbing materials for stealth aircraft | |
CN102254347A (zh) | 电磁环境体数据三维可视化方法 | |
CN104778286B (zh) | 掠海飞行器电磁散射特性快速仿真方法 | |
CN104915326A (zh) | 基于等效原理的区域分解阶数步进时域积分方法 | |
CN104346488B (zh) | 电大复杂外形金属目标混合建模及电磁散射快速仿真方法 | |
Tian et al. | Improved hybrid algorithm for composite scattering from multiple 3D objects above a 2D random dielectric rough surface | |
CN108090255B (zh) | 一种计算等离子体覆盖目标电磁散射的高低频混合方法 | |
Aghabarati et al. | Algebraic multigrid combined with domain decomposition for the finite element analysis of large scattering problems | |
CN106156394B (zh) | 基于显式差分格式的电磁特性提取方法 | |
CN114398716A (zh) | 一种热力电三场耦合下天线罩电性能全波数值计算方法 | |
Yuan et al. | Electromagnetic scattering by arbitrarily shaped stratified anisotropic media using the equivalent dipole moment method | |
CN111931353A (zh) | 一种应用于仿真fss结构的散射场求解方法 | |
Zhao et al. | Sparse matrix canonical grid method for three-dimension rough surface | |
CN104657527B (zh) | 超高速薄涂敷隐身飞行目标的电磁散射分析方法 | |
Minjie et al. | An efficient adaptive space partitioning algorithm for electromagnetic scattering calculation of complex 3D models | |
Gu et al. | Parallel Integral Equation Based Non-overlapping DDM for Fast Solving Electromagnetic Scattering Problems with Changeable Parts | |
Liu et al. | High-order dgtd for solving em scattering from hypersonic aircraft with plasma sheath |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |