raybet雷竞app|雷竞技官网下载|雷电竞下载苹果

AMP

数学物理的发展

1687 - 9139 1687 - 9120

Hindawi

10.1155 / 2017/6473234

6473234

研究文章

研究有效方法解决萨哈方程数值模拟的毛细管放电

http://orcid.org/0000 - 0002 - 8797 - 4144 刘

Hengzhu

http://orcid.org/0000 - 0001 - 7714 - 9928 吴

程

一个

Fengjiang

² 廖

沙沙村

² Daripa

Prabir

^{1

电机工程学院

北京理工学院

北京100081年

中国

bit.edu.cn} ^{2

爆炸科学与技术国家重点实验室

北京理工学院

北京100081年

中国

bit.edu.cn}

2017年

20. 11 2017年

2017年 14 06 2017年 17 08年 2017年 22 10 2017年 20. 11 2017年

这是一个开放的文章在知识共享归属许可下发布的,它允许无限制的使用,分布和繁殖在任何媒介,提供最初的工作是正确的引用。

数值模拟的毛细管放电的电热(ET)或ET-chemical(等)发射系统,需要解决萨哈方程来确定组成的等离子体传输系数。这将导致一个相对长的仿真运行时相比,传统的计算流体动力学问题。摘要相对困难的多维问题解决萨哈方程转化为一维问题的模拟毛细管放电通过构造迭代方程对温度。的介绍,以确保收敛系数迭代方程。为了提高计算效率,这个系数是进一步优化和动态设置迭代初始值的方法。数进行模拟试验研究这两种方法的性能。结果表明,仿真运行时可以显著降低本文中提出的方法。 1。介绍</t我tle> <p>电热(ET)或ET-chemical(等)发射技术可以实现出口速度超过2公里/秒(<gydF4y2Baxref ref-type="bibr" rid="B1"> 1</gydF4y2Baxref>]。超速的传统发射火药驱动的技术几乎可以实现。等的关键物理过程或等发射系统,为这一领域的研究总是集中在毛细管放电所产生的高压和高温等离子体去除管作为工作介质的墙壁或点燃推进剂(<gydF4y2Baxref ref-type="bibr" rid="B2"> 2</gydF4y2Baxref>]。在过去的几十年中,不同程度的简化数学模型,建立了基于欧拉方程来研究这个物理过程(<gydF4y2Baxref ref-type="bibr" rid="B3"> 3</gydF4y2Baxref>- - - - - -<gydF4y2Baxref ref-type="bibr" rid="B9"> 9</gydF4y2Baxref>]。作为一个可以找到这些模型,特征是这些司机术语或传输系数,如电导率和烧蚀率,都是有关等离子体组成,所以,需要确定在每个计算等离子体的组成细胞首先解决数学模型之前下一个时间步的等离子体特性。</pgydF4y2Ba> <p>等或等发射技术,典型的温度和压力范围内的毛细管放电产生的等离子体10000 - 40000 K和100 MPa,分别为(<gydF4y2Baxref ref-type="bibr" rid="B10"> 10</gydF4y2Baxref>]。在这样的极端条件下,等离子体中认为是局部热力学平衡(LTE)和熔化的材料完全分离成基本的成分,部分或全部电离(<gydF4y2Baxref ref-type="bibr" rid="B11"> 11</gydF4y2Baxref>]。这些假设,等离子体的成分可以被解决萨哈方程决定。萨哈方程定义的关系在两个相邻离子电离阶段。它通常是一个多维的问题如果直接解决它。如前所述,扎格卢勒(<gydF4y2Baxref ref-type="bibr" rid="B12"> 12</gydF4y2Baxref>),这将是很难找到解决方案。工作的Trayner和Głowacki [<gydF4y2Baxref ref-type="bibr" rid="B13"> 13</gydF4y2Baxref>]情况理想的萨哈方程与单一元素等离子体物种首先新配方为一维方程,然后用一些解决数值算法对分法等。之后,扎格卢勒扩展这种方法的情况不理想的等离子体与多元素物种,同时保持一个维数的方法<gydF4y2Baxref ref-type="bibr" rid="B12"> 12</gydF4y2Baxref>]。的方法Trayner Głowacki和扎格卢勒是突出的,因为他们的方法使解决萨哈方程单独的工作变得非常容易。</pgydF4y2Ba> <p>毛细管放电等中使用的模拟或等发射系统,采用的方法的问题Trayner Głowacki或扎格卢勒是那些谈话方程数学模型的毛细管放电不能直接给温度而应作为输入数据为解决萨哈方程。的系列研究Zoler et al。<gydF4y2Baxref ref-type="bibr" rid="B5"> 5</gydF4y2Baxref>,<gydF4y2Baxref ref-type="bibr" rid="B14"> 14</gydF4y2Baxref>),能源对话方程是通过变量代换使它能够直接给温度。然而,夫妻偏微分项也会引入模型,将导致额外的计算量计算这些术语。在其他作品<gydF4y2Baxref ref-type="bibr" rid="B4"> 4</gydF4y2Baxref>,<gydF4y2Baxref ref-type="bibr" rid="B9"> 9</gydF4y2Baxref>,<gydF4y2Baxref ref-type="bibr" rid="B15"> 15</gydF4y2Baxref>),萨哈方程直接解决了在仿真与特定的能量方程。在这篇文章中,这种直接的理想是保留。方法的基础上Trayner Głowacki和扎格卢勒,多维问题的解决萨哈方程模拟毛细管放电转化为一维问题通过修改特定能量方程为迭代方程对温度。</pgydF4y2Ba> <p>这种方法使解决萨哈方程的工作变得相当简单的模拟毛细管放电;然而,仍然需要解决迭代方程数值算法在每一个细胞。由于计算量,毛细管放电通常需要长时间的仿真与传统的计算流体动力学问题。特别是在长管的情况下可能联合的退出毛细管作为桶或阴极或其他用途<gydF4y2Baxref ref-type="bibr" rid="B4"> 4</gydF4y2Baxref>,<gydF4y2Baxref ref-type="bibr" rid="B16"> 16</gydF4y2Baxref>,<gydF4y2Baxref ref-type="bibr" rid="B17"> 17</gydF4y2Baxref>),更大的计算区域将进一步提高仿真运行时。本文介绍了两种方法更好的计算效率包括优化迭代方程中引入系数和动态设置迭代初始值。加速机制与这些方法也是数值计算研究。</pgydF4y2Ba> </sec> <sec id="sec2"> <title>2。萨哈方程</t我tle> <p>在仿真研究等的毛细管放电或等发射技术,萨哈方程的一种广泛使用的形式是由(<gydF4y2Baxref ref-type="bibr" rid="B15"> 15</gydF4y2Baxref>]<d我年代p- - - - - -for米ul一个> <mml:math display="block" xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M1"> <mml:mtable> <mml:mlabeledtr id="EEq1"> <mml:mtd> <mml:mtext> (1)</米米l:mtext> </mml:mtd> <mml:mtd> <mml:mfrac> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mi> 我</米米l:mi> <mml:mo> +</米米l:mo> <mml:mn mathvariant="normal"> 1</米米l:mn> </mml:mrow> </mml:msub> </mml:mrow> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mi> 我</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:mfrac> <mml:mo> =</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> K</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mi> 我</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> ,</米米l:mo> </mml:mtd> </mml:mlabeledtr> </mml:mtable> </mml:math> </disp-formula>在哪里<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M2"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mi> 我</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>的比例吗<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M3"> <mml:mrow> <mml:mi> 我</米米l:mi> </mml:mrow> </mml:math> </inline-formula>倍电离离子的元素种类<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M4"> <mml:mrow> <mml:mi> j</米米l:mi> </mml:mrow> </mml:math> </inline-formula>总重粒子(<gydF4y2Baxref ref-type="bibr" rid="B12"> 12</gydF4y2Baxref>]。显然,确定等离子体的成分通过求解萨哈方程实际上意味着确定<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M5"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mi> 我</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>通过求解(<gydF4y2Baxref ref-type="disp-formula" rid="EEq1"> 1</gydF4y2Baxref>)。的功能<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M6"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> K</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mi> 我</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>是由(<gydF4y2Baxref ref-type="bibr" rid="B12"> 12</gydF4y2Baxref>,<gydF4y2Baxref ref-type="bibr" rid="B15"> 15</gydF4y2Baxref>,<gydF4y2Baxref ref-type="bibr" rid="B18"> 18</gydF4y2Baxref>]<d我年代p- - - - - -for米ul一个> <mml:math display="block" xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M7"> <mml:mtable> <mml:mlabeledtr id="EEq2"> <mml:mtd> <mml:mtext> (2)</米米l:mtext> </mml:mtd> <mml:mtd> <mml:msub> <mml:mrow> <mml:mi> K</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mi> 我</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:mfrac> <mml:mrow> <mml:mn mathvariant="normal"> 1</米米l:mn> </mml:mrow> <mml:mrow> <mml:mi> n</米米l:mi> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> e</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:mfrac> <mml:mfrac> <mml:mrow> <mml:mn mathvariant="normal"> 2</米米l:mn> <mml:msub> <mml:mrow> <mml:mi> U</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mi> 我</米米l:mi> <mml:mo> +</米米l:mo> <mml:mn mathvariant="normal"> 1</米米l:mn> </mml:mrow> </mml:msub> </mml:mrow> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> U</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mi> 我</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:mfrac> <mml:msup> <mml:mrow> <mml:mfenced separators="|"> <mml:mrow> <mml:mfrac> <mml:mrow> <mml:mn mathvariant="normal"> 2</米米l:mn> <mml:mi> π</米米l:mi> <mml:msub> <mml:mrow> <mml:mi> 米</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> e</米米l:mi> </mml:mrow> </mml:msub> <mml:mi> k</米米l:mi> <mml:mi> T</米米l:mi> </mml:mrow> <mml:mrow> <mml:msup> <mml:mrow> <mml:mi> h</米米l:mi> </mml:mrow> <mml:mrow> <mml:mn mathvariant="normal"> 2</米米l:mn> </mml:mrow> </mml:msup> </mml:mrow> </mml:mfrac> </mml:mrow> </mml:mfenced> </mml:mrow> <mml:mrow> <mml:mrow> <mml:mrow> <mml:mn mathvariant="normal"> 3</米米l:mn> </mml:mrow> <mml:mo> /</米米l:mo> <mml:mrow> <mml:mn mathvariant="normal"> 2</米米l:mn> </mml:mrow> </mml:mrow> </mml:mrow> </mml:msup> <mml:mrow> <mml:mrow> <mml:mi mathvariant="normal"> 经验值</米米l:mi> </mml:mrow> <mml:mo> ⁡</米米l:mo> <mml:mrow> <mml:mfenced open="[" close="]" separators="|"> <mml:mrow> <mml:mfrac> <mml:mrow> <mml:mo> - - - - - -</米米l:mo> <mml:mfenced separators="|"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> 我</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mi> 我</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> - - - - - -</米米l:mo> <mml:mi> Δ</米米l:mi> <mml:msub> <mml:mrow> <mml:mi> 我</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mi> 我</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:mfenced> </mml:mrow> <mml:mrow> <mml:mi> k</米米l:mi> <mml:mi> T</米米l:mi> </mml:mrow> </mml:mfrac> </mml:mrow> </mml:mfenced> </mml:mrow> </mml:mrow> <mml:mo> ,</米米l:mo> </mml:mtd> </mml:mlabeledtr> </mml:mtable> </mml:math> </disp-formula>在哪里<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M8"> <mml:mrow> <mml:mi> n</米米l:mi> </mml:mrow> </mml:math> </inline-formula>是总重粒子的数密度,<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M9"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> e</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>电子比总重粒子,<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M10"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> U</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mi> 我</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>表示的配分函数依赖于温度<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M11"> <mml:mrow> <mml:mi> T</米米l:mi> </mml:mrow> </mml:math> </inline-formula>,<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M12"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> 米</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> e</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>电子的质量,<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M13"> <mml:mrow> <mml:mi> k</米米l:mi> </mml:mrow> </mml:math> </inline-formula>是玻尔兹曼常数,<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M14"> <mml:mrow> <mml:mi> h</米米l:mi> </mml:mrow> </mml:math> </inline-formula>普朗克常数,<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M15"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> 我</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mi> 我</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>电离势。的值为<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M16"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> 我</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mi> 我</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>和必要的光谱数据用于计算<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M17"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> U</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mi> 我</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>,一个可以在网上找到他们从数据库(<gydF4y2Baxref ref-type="bibr" rid="B19"> 19</gydF4y2Baxref>]。<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M18"> <mml:mi> Δ</米米l:mi> <mml:msub> <mml:mrow> <mml:mi> 我</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mi> 我</米米l:mi> </mml:mrow> </mml:msub> </mml:math> </inline-formula>介绍了电离势的降低,因为等离子体的非理想的效果。<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M19"> <mml:mi> Δ</米米l:mi> <mml:msub> <mml:mrow> <mml:mi> 我</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mi> 我</米米l:mi> </mml:mrow> </mml:msub> </mml:math> </inline-formula>是由(<gydF4y2Baxref ref-type="bibr" rid="B4"> 4</gydF4y2Baxref>,<gydF4y2Baxref ref-type="bibr" rid="B12"> 12</gydF4y2Baxref>]<d我年代p- - - - - -for米ul一个> <mml:math display="block" xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M20"> <mml:mtable> <mml:mlabeledtr id="EEq3"> <mml:mtd> <mml:mtext> (3)</米米l:mtext> </mml:mtd> <mml:mtd> <mml:mi> Δ</米米l:mi> <mml:msub> <mml:mrow> <mml:mi> 我</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mi> 我</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:mfrac> <mml:mrow> <mml:mfenced separators="|"> <mml:mrow> <mml:mi> 我</米米l:mi> <mml:mo> +</米米l:mo> <mml:mn mathvariant="normal"> 1</米米l:mn> </mml:mrow> </mml:mfenced> <mml:msup> <mml:mrow> <mml:mi> e</米米l:mi> </mml:mrow> <mml:mrow> <mml:mn mathvariant="normal"> 2</米米l:mn> </mml:mrow> </mml:msup> </mml:mrow> <mml:mrow> <mml:mn mathvariant="normal"> 4</米米l:mn> <mml:mi> π</米米l:mi> <mml:msub> <mml:mrow> <mml:mi> ε</米米l:mi> </mml:mrow> <mml:mrow> <mml:mn mathvariant="normal"> 0</米米l:mn> </mml:mrow> </mml:msub> <mml:mfenced separators="|"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> λ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> D</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> +</米米l:mo> <mml:mrow> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> Λ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> B</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> <mml:mo> /</米米l:mo> <mml:mrow> <mml:mn mathvariant="normal"> 8</米米l:mn> </mml:mrow> </mml:mrow> </mml:mrow> </mml:mfenced> </mml:mrow> </mml:mfrac> <mml:mo> ,</米米l:mo> </mml:mtd> </mml:mlabeledtr> </mml:mtable> </mml:math> </disp-formula>在哪里<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M21"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ε</米米l:mi> </mml:mrow> <mml:mrow> <mml:mn> 0</米米l:mn> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>的介电常数是自由的空间,<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M22"> <mml:mrow> <mml:mi> e</米米l:mi> </mml:mrow> </mml:math> </inline-formula>是电子电荷,<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M23"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> Λ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> B</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>是德布罗意波长,<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M24"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> λ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> D</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>是标准的德拜长度可表示为<d我年代p- - - - - -for米ul一个> <mml:math display="block" xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M25"> <mml:mtable> <mml:mlabeledtr id="EEq4"> <mml:mtd> <mml:mtext> (4)</米米l:mtext> </mml:mtd> <mml:mtd> <mml:msub> <mml:mrow> <mml:mi> λ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> D</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:msqrt> <mml:mfrac> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ε</米米l:mi> </mml:mrow> <mml:mrow> <mml:mn mathvariant="normal"> 0</米米l:mn> </mml:mrow> </mml:msub> <mml:mi> k</米米l:mi> <mml:mi> T</米米l:mi> </mml:mrow> <mml:mrow> <mml:msup> <mml:mrow> <mml:mi> e</米米l:mi> </mml:mrow> <mml:mrow> <mml:mn mathvariant="normal"> 2</米米l:mn> </mml:mrow> </mml:msup> <mml:mi> n</米米l:mi> <mml:mfenced separators="|"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> e</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> +</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 一个</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:mfenced> </mml:mrow> </mml:mfrac> </mml:msqrt> <mml:mo> ,</米米l:mo> </mml:mtd> </mml:mlabeledtr> </mml:mtable> </mml:math> </disp-formula>在哪里<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M26"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 一个</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>被定义为<d我年代p- - - - - -for米ul一个> <mml:math display="block" xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M27"> <mml:mtable> <mml:mlabeledtr id="EEq5"> <mml:mtd> <mml:mtext> (5)</米米l:mtext> </mml:mtd> <mml:mtd> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 一个</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:mrow> <mml:munder> <mml:mstyle displaystyle="true"> <mml:mo stretchy="false"> ∑</米米l:mo> </mml:mstyle> <mml:mrow> <mml:mi> j</米米l:mi> </mml:mrow> </mml:munder> <mml:mrow> <mml:mrow> <mml:munderover> <mml:mstyle displaystyle="true"> <mml:mo stretchy="false"> ∑</米米l:mo> </mml:mstyle> <mml:mrow> <mml:mi> 我</米米l:mi> <mml:mo> =</米米l:mo> <mml:mn mathvariant="normal"> 1</米米l:mn> </mml:mrow> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> z</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:munderover> <mml:mrow> <mml:msup> <mml:mrow> <mml:mi> 我</米米l:mi> </mml:mrow> <mml:mrow> <mml:mn mathvariant="normal"> 2</米米l:mn> </mml:mrow> </mml:msup> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mi> 我</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:mrow> </mml:mrow> </mml:mrow> <mml:mo> ,</米米l:mo> </mml:mtd> </mml:mlabeledtr> </mml:mtable> </mml:math> </disp-formula>在哪里<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M28"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> z</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>代表的最大允许元素的电离阶段物种<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M29"> <mml:mrow> <mml:mi> j</米米l:mi> </mml:mrow> </mml:math> </inline-formula>。<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M30"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> z</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>可以估计的上限温度在一个特定的仿真案例。两个补充方程求解给出了萨哈方程<d我年代p- - - - - -for米ul一个> <mml:math display="block" xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M31"> <mml:mtable> <mml:mlabeledtr id="EEq6"> <mml:mtd> <mml:mtext> (6)</米米l:mtext> </mml:mtd> <mml:mtd> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> e</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:mrow> <mml:munder> <mml:mstyle displaystyle="true"> <mml:mo stretchy="false"> ∑</米米l:mo> </mml:mstyle> <mml:mrow> <mml:mi> j</米米l:mi> </mml:mrow> </mml:munder> <mml:mrow> <mml:mrow> <mml:munderover> <mml:mstyle displaystyle="true"> <mml:mo stretchy="false"> ∑</米米l:mo> </mml:mstyle> <mml:mrow> <mml:mi> 我</米米l:mi> <mml:mo> =</米米l:mo> <mml:mn mathvariant="normal"> 1</米米l:mn> </mml:mrow> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> z</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:munderover> <mml:mrow> <mml:mi> 我</米米l:mi> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mi> 我</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:mrow> </mml:mrow> </mml:mrow> <mml:mo> ,</米米l:mo> </mml:mtd> </mml:mlabeledtr> <mml:mlabeledtr id="EEq7"> <mml:mtd> <mml:mtext> (7)</米米l:mtext> </mml:mtd> <mml:mtd> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:mrow> <mml:munderover> <mml:mstyle displaystyle="true"> <mml:mo stretchy="false"> ∑</米米l:mo> </mml:mstyle> <mml:mrow> <mml:mi> 我</米米l:mi> <mml:mo> =</米米l:mo> <mml:mn mathvariant="normal"> 0</米米l:mn> </mml:mrow> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> z</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:munderover> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mi> 我</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:mrow> <mml:mo> ,</米米l:mo> </mml:mtd> </mml:mlabeledtr> </mml:mtable> </mml:math> </disp-formula>(在哪里<gydF4y2Baxref ref-type="disp-formula" rid="EEq6"> 6</gydF4y2Baxref>)是等离子体的电中性条件。<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M32"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>元素种类的比例吗<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M33"> <mml:mrow> <mml:mi> j</米米l:mi> </mml:mrow> </mml:math> </inline-formula>总重粒子。<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M34"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>可以得到等离子体的密度与材料的化学分子式切除从墙上的毛细管。为解决萨哈方程表示为(<gydF4y2Baxref ref-type="disp-formula" rid="EEq1"> 1</gydF4y2Baxref>)、等离子体密度、数密度<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M35"> <mml:mrow> <mml:mi> n</米米l:mi> </mml:mrow> </mml:math> </inline-formula>或质量密度<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M36"> <mml:mrow> <mml:mi> ρ</米米l:mi> </mml:mrow> </mml:math> </inline-formula>,<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M37"> <mml:mrow> <mml:mi> T</米米l:mi> </mml:mrow> </mml:math> </inline-formula>需要的输入数据。上面的方程是用来说明解决萨哈方程的方法。更多细节关于萨哈方程,他们可以找到一个在<gydF4y2Baxref ref-type="bibr" rid="B4"> 4</gydF4y2Baxref>,<gydF4y2Baxref ref-type="bibr" rid="B12"> 12</gydF4y2Baxref>,<gydF4y2Baxref ref-type="bibr" rid="B15"> 15</gydF4y2Baxref>]。</pgydF4y2Ba> </sec> <sec id="sec3"> <title>3所示。单独解决萨哈方程的方法</t我tle> <p>下面列出的方法解决萨哈方程分别来自Trayner和Głowacki和扎格卢勒<gydF4y2Baxref ref-type="bibr" rid="B12"> 12</gydF4y2Baxref>,<gydF4y2Baxref ref-type="bibr" rid="B13"> 13</gydF4y2Baxref>]:元素的物种<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M38"> <mml:mrow> <mml:mi> j</米米l:mi> </mml:mrow> </mml:math> </inline-formula>,(<gydF4y2Baxref ref-type="disp-formula" rid="EEq1"> 1</gydF4y2Baxref>)可以表示为<d我年代p- - - - - -for米ul一个> <mml:math display="block" xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M39"> <mml:mtable> <mml:mlabeledtr id="EEq8"> <mml:mtd rowspan="4"> <mml:mtext> (8)</米米l:mtext> </mml:mtd> <mml:mtd> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mn mathvariant="normal"> 1</米米l:mn> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mn mathvariant="normal"> 0</米米l:mn> </mml:mrow> </mml:msub> <mml:mo> ·</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> K</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mn> 0</米米l:mn> </mml:mrow> </mml:msub> <mml:mo> ,</米米l:mo> </mml:mtd> </mml:mlabeledtr> <mml:mtr> <mml:mtd> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mn mathvariant="normal"> 2</米米l:mn> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mn mathvariant="normal"> 1</米米l:mn> </mml:mrow> </mml:msub> <mml:mo> ·</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> K</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mn> 1</米米l:mn> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mn> 0</米米l:mn> </mml:mrow> </mml:msub> <mml:mo> ·</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> K</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mn> 0</米米l:mn> </mml:mrow> </mml:msub> <mml:mo> ·</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> K</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mn> 1</米米l:mn> </mml:mrow> </mml:msub> <mml:mo> ,</米米l:mo> </mml:mtd> </mml:mtr> <mml:mtr> <mml:mtd> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mn> 3</米米l:mn> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mn> 2</米米l:mn> </mml:mrow> </mml:msub> <mml:mo> ·</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> K</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mn> 2</米米l:mn> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mn> 0</米米l:mn> </mml:mrow> </mml:msub> <mml:mo> ·</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> K</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mn> 0</米米l:mn> </mml:mrow> </mml:msub> <mml:mo> ·</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> K</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mn> 1</米米l:mn> </mml:mrow> </mml:msub> <mml:mo> ·</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> K</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mn> 2</米米l:mn> </mml:mrow> </mml:msub> <mml:mo> ,</米米l:mo> </mml:mtd> </mml:mtr> <mml:mtr> <mml:mtd> <mml:mo> ⋮</米米l:mo> <mml:mo> 。</米米l:mo> </mml:mtd> </mml:mtr> </mml:mtable> </mml:math> </disp-formula>也就是说,<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M40"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mi> 我</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula> <inline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M41"> <mml:mrow> <mml:mfenced separators="|"> <mml:mrow> <mml:mi> 我</米米l:mi> <mml:mo> =</米米l:mo> <mml:mn> 1、2</米米l:mn> <mml:mo> ,</米米l:mo> <mml:mo> …</米米l:mo> <mml:mo> ,</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> z</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:mfenced> </mml:mrow> </mml:math> </inline-formula>可以表示为<d我年代p- - - - - -for米ul一个> <mml:math display="block" xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M42"> <mml:mtable> <mml:mlabeledtr id="EEq9"> <mml:mtd> <mml:mtext> (9)</米米l:mtext> </mml:mtd> <mml:mtd> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mi> 我</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mn mathvariant="normal"> 0</米米l:mn> </mml:mrow> </mml:msub> <mml:mrow> <mml:munderover> <mml:mstyle displaystyle="true"> <mml:mo stretchy="false"> ∏</米米l:mo> </mml:mstyle> <mml:mrow> <mml:mi> k</米米l:mi> <mml:mo> =</米米l:mo> <mml:mn mathvariant="normal"> 0</米米l:mn> </mml:mrow> <mml:mrow> <mml:mi> 我</米米l:mi> <mml:mo> - - - - - -</米米l:mo> <mml:mn mathvariant="normal"> 1</米米l:mn> </mml:mrow> </mml:munderover> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> K</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mi> k</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:mrow> <mml:mo> 。</米米l:mo> </mml:mtd> </mml:mlabeledtr> </mml:mtable> </mml:math> </disp-formula>从(<gydF4y2Baxref ref-type="disp-formula" rid="EEq7"> 7</gydF4y2Baxref>),<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M43"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mn> 0</米米l:mn> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>可以表示为<d我年代p- - - - - -for米ul一个> <mml:math display="block" xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M44"> <mml:mtable> <mml:mlabeledtr id="EEq10"> <mml:mtd> <mml:mtext> (10)</米米l:mtext> </mml:mtd> <mml:mtd> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mn> 0</米米l:mn> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> - - - - - -</米米l:mo> <mml:mrow> <mml:munderover> <mml:mstyle displaystyle="true"> <mml:mo stretchy="false"> ∑</米米l:mo> </mml:mstyle> <mml:mrow> <mml:mi> 我</米米l:mi> <mml:mo> =</米米l:mo> <mml:mn mathvariant="normal"> 1</米米l:mn> </mml:mrow> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> z</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:munderover> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mi> 我</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:mrow> <mml:mo> 。</米米l:mo> </mml:mtd> </mml:mlabeledtr> </mml:mtable> </mml:math> </disp-formula>用(<gydF4y2Baxref ref-type="disp-formula" rid="EEq9"> 9</gydF4y2Baxref>)(<gydF4y2Baxref ref-type="disp-formula" rid="EEq10"> 10</gydF4y2Baxref>),一些重组,其中一个<d我年代p- - - - - -for米ul一个> <mml:math display="block" xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M45"> <mml:mtable> <mml:mlabeledtr id="EEq11"> <mml:mtd> <mml:mtext> (11)</米米l:mtext> </mml:mtd> <mml:mtd> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mn mathvariant="normal"> 0</米米l:mn> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:mfrac> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> <mml:mrow> <mml:mfenced separators="|"> <mml:mrow> <mml:mn mathvariant="normal"> 1</米米l:mn> <mml:mo> +</米米l:mo> <mml:mrow> <mml:msubsup> <mml:mo stretchy="false"> ∑</米米l:mo> <mml:mrow> <mml:mi> 我</米米l:mi> <mml:mo> =</米米l:mo> <mml:mn mathvariant="normal"> 1</米米l:mn> </mml:mrow> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> z</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:msubsup> <mml:mrow> <mml:mrow> <mml:msubsup> <mml:mo stretchy="false"> ∏</米米l:mo> <mml:mrow> <mml:mi> k</米米l:mi> <mml:mo> =</米米l:mo> <mml:mn> 0</米米l:mn> </mml:mrow> <mml:mrow> <mml:mi> 我</米米l:mi> <mml:mo> - - - - - -</米米l:mo> <mml:mn> 1</米米l:mn> </mml:mrow> </mml:msubsup> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> K</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mi> k</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:mrow> </mml:mrow> </mml:mrow> </mml:mrow> </mml:mfenced> </mml:mrow> </mml:mfrac> <mml:mo> 。</米米l:mo> </mml:mtd> </mml:mlabeledtr> </mml:mtable> </mml:math> </disp-formula>首先,对于一个更简单的情况<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M46"> <mml:mi> Δ</米米l:mi> <mml:msub> <mml:mrow> <mml:mi> 我</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mi> 我</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:mn> 0</米米l:mn> </mml:math> </inline-formula>理想的等离子体,<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M47"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> e</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>在<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M48"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> K</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mi> k</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>作为一个可以找到从未知(<gydF4y2Baxref ref-type="disp-formula" rid="EEq2"> 2</gydF4y2Baxref>),所以,<gydF4y2Baxref ref-type="disp-formula" rid="EEq11"> 11</gydF4y2Baxref>)可以表示为<d我年代p- - - - - -for米ul一个> <mml:math display="block" xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M49"> <mml:mtable> <mml:mlabeledtr id="EEq12"> <mml:mtd> <mml:mtext> (12)</米米l:mtext> </mml:mtd> <mml:mtd> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mn mathvariant="normal"> 0</米米l:mn> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> f</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mn> 0</米米l:mn> </mml:mrow> </mml:msub> <mml:mfenced separators="|"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> e</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:mfenced> <mml:mo> 。</米米l:mo> </mml:mtd> </mml:mlabeledtr> </mml:mtable> </mml:math> </disp-formula>方程(<gydF4y2Baxref ref-type="disp-formula" rid="EEq12"> 12</gydF4y2Baxref>)表明,<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M50"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mn> 0</米米l:mn> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>可以明确表达的吗<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M51"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> e</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>。用(<gydF4y2Baxref ref-type="disp-formula" rid="EEq12"> 12</gydF4y2Baxref>)(<gydF4y2Baxref ref-type="disp-formula" rid="EEq9"> 9</gydF4y2Baxref>),一个人<d我年代p- - - - - -for米ul一个> <mml:math display="block" xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M52"> <mml:mtable> <mml:mlabeledtr id="EEq13"> <mml:mtd> <mml:mtext> (13)</米米l:mtext> </mml:mtd> <mml:mtd> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mi> 我</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> f</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mn> 0</米米l:mn> </mml:mrow> </mml:msub> <mml:mfenced separators="|"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> e</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:mfenced> <mml:mrow> <mml:munderover> <mml:mstyle displaystyle="true"> <mml:mo stretchy="false"> ∏</米米l:mo> </mml:mstyle> <mml:mrow> <mml:mi> k</米米l:mi> <mml:mo> =</米米l:mo> <mml:mn mathvariant="normal"> 0</米米l:mn> </mml:mrow> <mml:mrow> <mml:mi> 我</米米l:mi> <mml:mo> - - - - - -</米米l:mo> <mml:mn mathvariant="normal"> 1</米米l:mn> </mml:mrow> </mml:munderover> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> K</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mi> k</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:mrow> <mml:mo> ,</米米l:mo> </mml:mtd> </mml:mlabeledtr> </mml:mtable> </mml:math> </disp-formula>在哪里<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M53"> <mml:mi> 我</米米l:mi> <mml:mo> =</米米l:mo> <mml:mn> 1、2</米米l:mn> <mml:mo> ,</米米l:mo> <mml:mo> …</米米l:mo> <mml:mo> ,</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> z</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> </mml:mrow> </mml:msub> </mml:math> </inline-formula>。同样,作为(<gydF4y2Baxref ref-type="disp-formula" rid="EEq11"> 11</gydF4y2Baxref>),只有<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M54"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> e</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>在右侧(<gydF4y2Baxref ref-type="disp-formula" rid="EEq13"> 13</gydF4y2Baxref>)是未知的表明<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M55"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mi> 我</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula> <inline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M56"> <mml:mrow> <mml:mfenced separators="|"> <mml:mrow> <mml:mi> 我</米米l:mi> <mml:mo> =</米米l:mo> <mml:mn> 1、2</米米l:mn> <mml:mo> ,</米米l:mo> <mml:mo> …</米米l:mo> <mml:mo> ,</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> z</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:mfenced> </mml:mrow> </mml:math> </inline-formula>也可以表达明确的吗<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M57"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> e</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>。然后,在一起(<gydF4y2Baxref ref-type="disp-formula" rid="EEq12"> 12</gydF4y2Baxref>),所有<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M58"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mi> 我</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula> <inline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M59"> <mml:mrow> <mml:mfenced separators="|"> <mml:mrow> <mml:mi> 我</米米l:mi> <mml:mo> =</米米l:mo> <mml:mn> 0,2</米米l:mn> <mml:mo> ,</米米l:mo> <mml:mo> …</米米l:mo> <mml:mo> ,</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> z</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:mfenced> </mml:mrow> </mml:math> </inline-formula>元素的种类<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M60"> <mml:mrow> <mml:mi> j</米米l:mi> </mml:mrow> </mml:math> </inline-formula>可以明确表达的吗<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M61"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> e</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>。也就是说,<d我年代p- - - - - -for米ul一个> <mml:math display="block" xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M62"> <mml:mtable> <mml:mlabeledtr id="EEq14"> <mml:mtd> <mml:mtext> (14)</米米l:mtext> </mml:mtd> <mml:mtd> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mi> 我</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> f</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mi> 我</米米l:mi> </mml:mrow> </mml:msub> <mml:mfenced separators="|"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> e</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:mfenced> <mml:mo> 。</米米l:mo> </mml:mtd> </mml:mlabeledtr> </mml:mtable> </mml:math> </disp-formula>自(<gydF4y2Baxref ref-type="disp-formula" rid="EEq14"> 14</gydF4y2Baxref>)对任何元素实际上是有效的物种<我t一个l我c> j</我t一个l我c>被认为是在模拟的情况下,(<gydF4y2Baxref ref-type="disp-formula" rid="EEq14"> 14</gydF4y2Baxref>)也意味着<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M63"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mi> 我</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula> <inline-formula> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M64"> <mml:mrow> <mml:mfenced separators="|"> <mml:mrow> <mml:mi> 我</米米l:mi> <mml:mo> =</米米l:mo> <mml:mn> 0,2</米米l:mn> <mml:mo> ,</米米l:mo> <mml:mo> …</米米l:mo> <mml:mo> ,</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> z</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:mfenced> </mml:mrow> </mml:math> </inline-formula>每个元素的种类<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M65"> <mml:mrow> <mml:mi> j</米米l:mi> </mml:mrow> </mml:math> </inline-formula>考虑到仿真可以明确的表达<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M66"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> e</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>。然后,用(<gydF4y2Baxref ref-type="disp-formula" rid="EEq14"> 14</gydF4y2Baxref>)(<gydF4y2Baxref ref-type="disp-formula" rid="EEq6"> 6</gydF4y2Baxref>),可以得到迭代方程<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M67"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> e</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>。也就是说,<d我年代p- - - - - -for米ul一个> <mml:math display="block" xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M68"> <mml:mtable> <mml:mlabeledtr id="EEq15"> <mml:mtd> <mml:mtext> (15)</米米l:mtext> </mml:mtd> <mml:mtd> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> e</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:mrow> <mml:munder> <mml:mstyle displaystyle="true"> <mml:mo stretchy="false"> ∑</米米l:mo> </mml:mstyle> <mml:mrow> <mml:mi> j</米米l:mi> </mml:mrow> </mml:munder> <mml:mrow> <mml:mrow> <mml:munderover> <mml:mstyle displaystyle="true"> <mml:mo stretchy="false"> ∑</米米l:mo> </mml:mstyle> <mml:mrow> <mml:mi> 我</米米l:mi> <mml:mo> =</米米l:mo> <mml:mn mathvariant="normal"> 1</米米l:mn> </mml:mrow> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> z</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:munderover> <mml:mrow> <mml:mi> 我</米米l:mi> <mml:msub> <mml:mrow> <mml:mi> f</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mi> 我</米米l:mi> </mml:mrow> </mml:msub> <mml:mfenced separators="|"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> e</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:mfenced> </mml:mrow> </mml:mrow> </mml:mrow> </mml:mrow> <mml:mo> 。</米米l:mo> </mml:mtd> </mml:mlabeledtr> </mml:mtable> </mml:math> </disp-formula>本文简单的定点迭代法(<gydF4y2Baxref ref-type="bibr" rid="B20"> 20.</gydF4y2Baxref>)是用来解决迭代方程(<gydF4y2Baxref ref-type="disp-formula" rid="EEq15"> 15</gydF4y2Baxref>)。过程如下:<l是t> <list-item> <label>(我)</l一个bel> </list-item> </list></p> <p>集<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M69"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> e</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>迭代初始值。</pgydF4y2Ba> <list-item> <label>(2)</l一个bel> <p>计算(<gydF4y2Baxref ref-type="disp-formula" rid="EEq11"> 11</gydF4y2Baxref>)获得<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M70"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mn> 0</米米l:mn> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>。</pgydF4y2Ba> </list-item> <list-item> <label>(3)</l一个bel> <p>计算(<gydF4y2Baxref ref-type="disp-formula" rid="EEq9"> 9</gydF4y2Baxref>)获得<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M71"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mi> 我</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>。</pgydF4y2Ba> </list-item> <list-item> <label>(iv)</l一个bel> <p>重复步骤(2),(3)对所有元素的物种被认为是在模拟的情况下。</pgydF4y2Ba> </list-item> <list-item> <label>(v)</l一个bel> <p>计算(<gydF4y2Baxref ref-type="disp-formula" rid="EEq6"> 6</gydF4y2Baxref>)获得一个新的<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M72"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> e</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>表示为<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M73"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> e</米米l:mi> <mml:mo> _</米米l:mo> <mml:mi mathvariant="normal"> r</米米l:mi> <mml:mi mathvariant="normal"> e</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>。</pgydF4y2Ba> </list-item> <list-item> <label>(vi)</l一个bel> <p>集<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M74"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> e</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>是<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M75"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> e</米米l:mi> <mml:mo> _</米米l:mo> <mml:mi mathvariant="normal"> r</米米l:mi> <mml:mi mathvariant="normal"> e</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>重复步骤(2),直到满足终止条件(v)。</pgydF4y2Ba> </list-item> <p></p> <p>作为一个可以找到从程序(i) - (vi),<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M76"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mi> 我</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>同时可以确定吗<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M77"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> e</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>。也就是说,的情况下<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M78"> <mml:mi> Δ</米米l:mi> <mml:msub> <mml:mrow> <mml:mi> 我</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mi> 我</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:mn> 0</米米l:mn> </mml:math> </inline-formula>,一个相对困难的多维问题降低到一维问题。对于目前的情况<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M79"> <mml:mi> Δ</米米l:mi> <mml:msub> <mml:mrow> <mml:mi> 我</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mi> 我</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> ≠</米米l:mo> <mml:mn> 0</米米l:mn> </mml:math> </inline-formula>我们可以看到从(<gydF4y2Baxref ref-type="disp-formula" rid="EEq2"> 2</gydF4y2Baxref>)- (<gydF4y2Baxref ref-type="disp-formula" rid="EEq5"> 5</gydF4y2Baxref>),<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M80"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 一个</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>不能明确的表达出来<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M81"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> e</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>除非<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M82"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> z</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>是1,所以这两个吗<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M83"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> e</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>和<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M84"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 一个</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>是未知的<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M85"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> K</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mi> k</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>。在这种情况下,(<gydF4y2Baxref ref-type="disp-formula" rid="EEq12"> 12</gydF4y2Baxref>)和(<gydF4y2Baxref ref-type="disp-formula" rid="EEq14"> 14</gydF4y2Baxref>)成为<d我年代p- - - - - -for米ul一个> <mml:math display="block" xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M86"> <mml:mtable> <mml:mlabeledtr id="EEq16"> <mml:mtd rowspan="2"> <mml:mtext> (16)</米米l:mtext> </mml:mtd> <mml:mtd> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mn mathvariant="normal"> 0</米米l:mn> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> g</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mn> 0</米米l:mn> </mml:mrow> </mml:msub> <mml:mfenced separators="|"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> e</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> ,</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 一个</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:mfenced> <mml:mo> ,</米米l:mo> </mml:mtd> </mml:mlabeledtr> <mml:mtr> <mml:mtd> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mi> 我</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> g</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mi> 我</米米l:mi> </mml:mrow> </mml:msub> <mml:mfenced separators="|"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> e</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> ,</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 一个</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:mfenced> <mml:mo> 。</米米l:mo> </mml:mtd> </mml:mtr> </mml:mtable> </mml:math> </disp-formula>然后,(<gydF4y2Baxref ref-type="disp-formula" rid="EEq15"> 15</gydF4y2Baxref>)成为<d我年代p- - - - - -for米ul一个> <mml:math display="block" xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M87"> <mml:mtable> <mml:mlabeledtr id="EEq18"> <mml:mtd> <mml:mtext> (17)</米米l:mtext> </mml:mtd> <mml:mtd> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> e</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:mrow> <mml:munder> <mml:mstyle displaystyle="true"> <mml:mo stretchy="false"> ∑</米米l:mo> </mml:mstyle> <mml:mrow> <mml:mi> j</米米l:mi> </mml:mrow> </mml:munder> <mml:mrow> <mml:mrow> <mml:munderover> <mml:mstyle displaystyle="true"> <mml:mo stretchy="false"> ∑</米米l:mo> </mml:mstyle> <mml:mrow> <mml:mi> 我</米米l:mi> <mml:mo> =</米米l:mo> <mml:mn mathvariant="normal"> 1</米米l:mn> </mml:mrow> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> z</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:munderover> <mml:mrow> <mml:mi> 我</米米l:mi> <mml:msub> <mml:mrow> <mml:mi> g</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mi> 我</米米l:mi> </mml:mrow> </mml:msub> <mml:mfenced separators="|"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> e</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> ,</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 一个</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:mfenced> </mml:mrow> </mml:mrow> </mml:mrow> </mml:mrow> <mml:mo> 。</米米l:mo> </mml:mtd> </mml:mlabeledtr> </mml:mtable> </mml:math> </disp-formula>显然,(<gydF4y2Baxref ref-type="disp-formula" rid="EEq18"> 17</gydF4y2Baxref>)无法解决程序(我)——直接(vi)。然而,如果<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M88"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 一个</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>人为地将是一个经验值,<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M89"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> e</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>和<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M90"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mi> 我</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>可能仍然得到解决(<gydF4y2Baxref ref-type="disp-formula" rid="EEq18"> 17</gydF4y2Baxref>)程序(i) - (vi)作为解决(<gydF4y2Baxref ref-type="disp-formula" rid="EEq15"> 15</gydF4y2Baxref>)。然后,一个新的<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M91"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 一个</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>可以通过计算(<gydF4y2Baxref ref-type="disp-formula" rid="EEq5"> 5</gydF4y2Baxref>)。这实际上是一个迭代的过程<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M92"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 一个</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>可以表示为<d我年代p- - - - - -for米ul一个> <mml:math display="block" xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M93"> <mml:mtable> <mml:mlabeledtr id="EEq19"> <mml:mtd> <mml:mtext> (18)</米米l:mtext> </mml:mtd> <mml:mtd> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 一个</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:mover accent="false"> <mml:mrow> <mml:mi> g</米米l:mi> </mml:mrow> <mml:mo> ¯</米米l:mo> </mml:mover> <mml:mfenced separators="|"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 一个</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:mfenced> <mml:mo> 。</米米l:mo> </mml:mtd> </mml:mlabeledtr> </mml:mtable> </mml:math> </disp-formula>酒吧在<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M94"> <mml:mrow> <mml:mi> g</米米l:mi> </mml:mrow> </mml:math> </inline-formula>指的是数值函数的性质(<gydF4y2Baxref ref-type="bibr" rid="B12"> 12</gydF4y2Baxref>]。解决的过程(<gydF4y2Baxref ref-type="disp-formula" rid="EEq19"> 18</gydF4y2Baxref>)如下:<l是t> <list-item> <label>(一)</l一个bel> </list-item> </list></p> <p>集<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M95"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 一个</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>迭代初始值。</pgydF4y2Ba> <list-item> <label>(b)</l一个bel> <p>解决(<gydF4y2Baxref ref-type="disp-formula" rid="EEq18"> 17</gydF4y2Baxref>)程序(i) - (vi)获得<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M96"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mi> 我</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>。</pgydF4y2Ba> </list-item> <list-item> <label>(c)</l一个bel> <p>计算(<gydF4y2Baxref ref-type="disp-formula" rid="EEq5"> 5</gydF4y2Baxref>)获得一个新的<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M97"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 一个</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>表示为<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M98"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 一个</米米l:mi> <mml:mo> _</米米l:mo> <mml:mi mathvariant="normal"> r</米米l:mi> <mml:mi mathvariant="normal"> e</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>。</pgydF4y2Ba> </list-item> <list-item> <label>(d)</l一个bel> <p>集<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M99"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 一个</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>是<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M100"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 一个</米米l:mi> <mml:mo> _</米米l:mo> <mml:mi mathvariant="normal"> r</米米l:mi> <mml:mi mathvariant="normal"> e</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>和重复步骤(b) - (c),直到满足终止条件。</pgydF4y2Ba> </list-item> <p></p> <p>扎格卢勒如前所述,上面的方法将相对困难的多维问题转化为一维问题<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M101"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 一个</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>它很容易找到解决任意精度(<gydF4y2Baxref ref-type="bibr" rid="B12"> 12</gydF4y2Baxref>]。</pgydF4y2Ba> </sec> <sec id="sec4"> <title>4所示。解决萨哈方程模拟的毛细管放电</t我tle> <p> <disp-formula> <mml:math display="block" xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M102"> <mml:mtable> <mml:mlabeledtr id="EEq20"> <mml:mtd> <mml:mtext> (19)</米米l:mtext> </mml:mtd> <mml:mtd> <mml:mfrac> <mml:mrow> <mml:mo> ∂</米米l:mo> <mml:mi> ρ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mo> ∂</米米l:mo> <mml:mi> t</米米l:mi> </mml:mrow> </mml:mfrac> <mml:mo> +</米米l:mo> <mml:mi> w</米米l:mi> <mml:mfrac> <mml:mrow> <mml:mo> ∂</米米l:mo> <mml:mi> ρ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mo> ∂</米米l:mo> <mml:mi> z</米米l:mi> </mml:mrow> </mml:mfrac> <mml:mo> +</米米l:mo> <mml:mi> ρ</米米l:mi> <mml:mfrac> <mml:mrow> <mml:mo> ∂</米米l:mo> <mml:mi> w</米米l:mi> </mml:mrow> <mml:mrow> <mml:mo> ∂</米米l:mo> <mml:mi> z</米米l:mi> </mml:mrow> </mml:mfrac> <mml:mo> =</米米l:mo> <mml:msubsup> <mml:mrow> <mml:mi> ρ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> α</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> ′</米米l:mi> </mml:mrow> </mml:msubsup> <mml:mo> ,</米米l:mo> </mml:mtd> </mml:mlabeledtr> <mml:mlabeledtr id="EEq21"> <mml:mtd> <mml:mtext> (20)</米米l:mtext> </mml:mtd> <mml:mtd> <mml:mi> ρ</米米l:mi> <mml:mfrac> <mml:mrow> <mml:mo> ∂</米米l:mo> <mml:mi> w</米米l:mi> </mml:mrow> <mml:mrow> <mml:mo> ∂</米米l:mo> <mml:mi> t</米米l:mi> </mml:mrow> </mml:mfrac> <mml:mo> +</米米l:mo> <mml:mi> ρ</米米l:mi> <mml:mi> w</米米l:mi> <mml:mfrac> <mml:mrow> <mml:mo> ∂</米米l:mo> <mml:mi> w</米米l:mi> </mml:mrow> <mml:mrow> <mml:mo> ∂</米米l:mo> <mml:mi> z</米米l:mi> </mml:mrow> </mml:mfrac> <mml:mo> +</米米l:mo> <mml:mfrac> <mml:mrow> <mml:mo> ∂</米米l:mo> <mml:mi> p</米米l:mi> </mml:mrow> <mml:mrow> <mml:mo> ∂</米米l:mo> <mml:mi> z</米米l:mi> </mml:mrow> </mml:mfrac> <mml:mo> =</米米l:mo> <mml:mo> - - - - - -</米米l:mo> <mml:mi> w</米米l:mi> <mml:msubsup> <mml:mrow> <mml:mi> ρ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> α</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> ′</米米l:mi> </mml:mrow> </mml:msubsup> <mml:mo> ,</米米l:mo> </mml:mtd> </mml:mlabeledtr> <mml:mlabeledtr id="EEq22"> <mml:mtd> <mml:mtext> (21)</米米l:mtext> </mml:mtd> <mml:mtd> <mml:mi> ρ</米米l:mi> <mml:mfrac> <mml:mrow> <mml:mo> ∂</米米l:mo> <mml:mi> E</米米l:mi> </mml:mrow> <mml:mrow> <mml:mo> ∂</米米l:mo> <mml:mi> t</米米l:mi> </mml:mrow> </mml:mfrac> <mml:mo> +</米米l:mo> <mml:mi> ρ</米米l:mi> <mml:mi> w</米米l:mi> <mml:mfrac> <mml:mrow> <mml:mo> ∂</米米l:mo> <mml:mi> E</米米l:mi> </mml:mrow> <mml:mrow> <mml:mo> ∂</米米l:mo> <mml:mi> z</米米l:mi> </mml:mrow> </mml:mfrac> <mml:mo> +</米米l:mo> <mml:mi> p</米米l:mi> <mml:mfrac> <mml:mrow> <mml:mo> ∂</米米l:mo> <mml:mi> w</米米l:mi> </mml:mrow> <mml:mrow> <mml:mo> ∂</米米l:mo> <mml:mi> z</米米l:mi> </mml:mrow> </mml:mfrac> <mml:mo> =</米米l:mo> <mml:mfrac> <mml:mrow> <mml:msup> <mml:mrow> <mml:mi> J</米米l:mi> </mml:mrow> <mml:mrow> <mml:mn> 2</米米l:mn> </mml:mrow> </mml:msup> </mml:mrow> <mml:mrow> <mml:mi> σ</米米l:mi> </mml:mrow> </mml:mfrac> <mml:mo> - - - - - -</米米l:mo> <mml:msubsup> <mml:mrow> <mml:mi> ρ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> α</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> ′</米米l:mi> </mml:mrow> </mml:msubsup> <mml:mfenced separators="|"> <mml:mrow> <mml:mi> E</米米l:mi> <mml:mo> - - - - - -</米米l:mo> <mml:mfrac> <mml:mrow> <mml:msup> <mml:mrow> <mml:mi> w</米米l:mi> </mml:mrow> <mml:mrow> <mml:mn mathvariant="normal"> 2</米米l:mn> </mml:mrow> </mml:msup> </mml:mrow> <mml:mrow> <mml:mn mathvariant="normal"> 2</米米l:mn> </mml:mrow> </mml:mfrac> </mml:mrow> </mml:mfenced> <mml:mo> ,</米米l:mo> </mml:mtd> </mml:mlabeledtr> <mml:mlabeledtr id="EEq23"> <mml:mtd> <mml:mtext> (22)</米米l:mtext> </mml:mtd> <mml:mtd> <mml:mi> E</米米l:mi> <mml:mo> =</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> E</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> t</米米l:mi> <mml:mi> r</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> +</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> E</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> e</米米l:mi> <mml:mi> l</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> +</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> E</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> c</米米l:mi> <mml:mi> h</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> +</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> E</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> v</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> +</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> E</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> d</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> ,</米米l:mo> </mml:mtd> </mml:mlabeledtr> <mml:mlabeledtr id="EEq24"> <mml:mtd> <mml:mtext> (23)</米米l:mtext> </mml:mtd> <mml:mtd> <mml:mi> p</米米l:mi> <mml:mo> =</米米l:mo> <mml:mfenced separators="|"> <mml:mrow> <mml:mi> n</米米l:mi> <mml:mo> +</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> n</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> e</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:mfenced> <mml:mi> k</米米l:mi> <mml:mi> T</米米l:mi> <mml:mo> 。</米米l:mo> </mml:mtd> </mml:mlabeledtr> </mml:mtable> </mml:math> </disp-formula></p> <p>方程(<gydF4y2Baxref ref-type="disp-formula" rid="EEq20"> 19</gydF4y2Baxref>)- (<gydF4y2Baxref ref-type="disp-formula" rid="EEq24"> 23</gydF4y2Baxref>)是时间的一维数学模型用来模拟毛细管放电等或等发射技术(<gydF4y2Baxref ref-type="bibr" rid="B4"> 4</gydF4y2Baxref>,<gydF4y2Baxref ref-type="bibr" rid="B11"> 11</gydF4y2Baxref>,<gydF4y2Baxref ref-type="bibr" rid="B21"> 21</gydF4y2Baxref>]。在这组模型,<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M103"> <mml:mrow> <mml:mi> t</米米l:mi> </mml:mrow> </mml:math> </inline-formula>代表时间和<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M104"> <mml:mrow> <mml:mi> z</米米l:mi> </mml:mrow> </mml:math> </inline-formula>代表轴坐标,<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M105"> <mml:mrow> <mml:mi> ρ</米米l:mi> </mml:mrow> </mml:math> </inline-formula>质量密度,<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M106"> <mml:mrow> <mml:mi> w</米米l:mi> </mml:mrow> </mml:math> </inline-formula>等离子体流的速度,<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M107"> <mml:mrow> <mml:msubsup> <mml:mrow> <mml:mi> ρ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> α</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> ′</米米l:mi> </mml:mrow> </mml:msubsup> </mml:mrow> </mml:math> </inline-formula>烧蚀率的增加导致吗<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M108"> <mml:mrow> <mml:mi> ρ</米米l:mi> </mml:mrow> </mml:math> </inline-formula>,<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M109"> <mml:mrow> <mml:mi> p</米米l:mi> </mml:mrow> </mml:math> </inline-formula>的压力,<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M110"> <mml:mrow> <mml:mi> E</米米l:mi> </mml:mrow> </mml:math> </inline-formula>具体的内部能量(您),<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M111"> <mml:mrow> <mml:mi> σ</米米l:mi> </mml:mrow> </mml:math> </inline-formula>是电导率,<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M112"> <mml:msup> <mml:mrow> <mml:mi> J</米米l:mi> </mml:mrow> <mml:mrow> <mml:mn> 2</米米l:mn> </mml:mrow> </mml:msup> <mml:mo> /</米米l:mo> <mml:mi> σ</米米l:mi> </mml:math> </inline-formula>是输入能量焦耳加热的形式。<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M113"> <mml:mrow> <mml:mi> σ</米米l:mi> </mml:mrow> </mml:math> </inline-formula>确定基于电子碰撞中性原子和离子的碰撞频率与离子评估模型Zollweg和利伯曼模型(<gydF4y2Baxref ref-type="bibr" rid="B22"> 22</gydF4y2Baxref>]。详细的理论或计算的方法<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M114"> <mml:mrow> <mml:mi> σ</米米l:mi> </mml:mrow> </mml:math> </inline-formula>以及传输系数等<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M115"> <mml:mrow> <mml:msubsup> <mml:mrow> <mml:mi> ρ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> α</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> ′</米米l:mi> </mml:mrow> </mml:msubsup> </mml:mrow> </mml:math> </inline-formula>和这些条款您方程(<gydF4y2Baxref ref-type="disp-formula" rid="EEq23"> 22</gydF4y2Baxref>从(),他们可以找到一个<gydF4y2Baxref ref-type="bibr" rid="B4"> 4</gydF4y2Baxref>,<gydF4y2Baxref ref-type="bibr" rid="B15"> 15</gydF4y2Baxref>]。如前所述,这些传输系数都有关<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M116"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mi> 我</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>,所以需要解决萨哈方程在每个计算单元来确定<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M117"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mi> 我</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>首先,然后,(<gydF4y2Baxref ref-type="disp-formula" rid="EEq20"> 19</gydF4y2Baxref>)- (<gydF4y2Baxref ref-type="disp-formula" rid="EEq24"> 23</gydF4y2Baxref>)可以得到解决,下一个时间步的等离子体流特性。然而,(<gydF4y2Baxref ref-type="disp-formula" rid="EEq20"> 19</gydF4y2Baxref>)- (<gydF4y2Baxref ref-type="disp-formula" rid="EEq22"> 21</gydF4y2Baxref>)不能给<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M118"> <mml:mrow> <mml:mi> T</米米l:mi> </mml:mrow> </mml:math> </inline-formula>直接而<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M119"> <mml:mrow> <mml:mi> T</米米l:mi> </mml:mrow> </mml:math> </inline-formula>需要的输入数据为解决萨哈方程。同样,作为建设(<gydF4y2Baxref ref-type="disp-formula" rid="EEq18"> 17</gydF4y2Baxref>新出现。<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M120"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 一个</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>对于非理想的等离子体前面提到的,该方法Trayner Głowacki和扎格卢勒可以进一步扩展通过构造迭代方程<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M121"> <mml:mrow> <mml:mi> T</米米l:mi> </mml:mrow> </mml:math> </inline-formula>同时保持方法的一个维度。的迭代方程<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M122"> <mml:mrow> <mml:mi> T</米米l:mi> </mml:mrow> </mml:math> </inline-formula>可以很容易地获得通过修改(<gydF4y2Baxref ref-type="disp-formula" rid="EEq23"> 22</gydF4y2Baxref>)<d我年代p- - - - - -for米ul一个> <mml:math display="block" xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M123"> <mml:mtable> <mml:mlabeledtr id="EEq25"> <mml:mtd> <mml:mtext> (24)</米米l:mtext> </mml:mtd> <mml:mtd> <mml:mi> T</米米l:mi> <mml:mo> =</米米l:mo> <mml:mfenced separators="|"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> E</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi mathvariant="normal"> 我</米米l:mi> <mml:mi mathvariant="normal"> n</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> - - - - - -</米米l:mo> <mml:mi> E</米米l:mi> <mml:mfenced separators="|"> <mml:mrow> <mml:mi> T</米米l:mi> <mml:mo> ,</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> ρ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi mathvariant="normal"> 我</米米l:mi> <mml:mi mathvariant="normal"> n</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:mfenced> </mml:mrow> </mml:mfenced> <mml:mo> ·</米米l:mo> <mml:mi> δ</米米l:mi> <mml:mo> +</米米l:mo> <mml:mi> T</米米l:mi> <mml:mo> ,</米米l:mo> </mml:mtd> </mml:mlabeledtr> </mml:mtable> </mml:math> </disp-formula>在哪里<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M124"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> E</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi mathvariant="normal"> 我</米米l:mi> <mml:mi mathvariant="normal"> n</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>和<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M125"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ρ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi mathvariant="normal"> 我</米米l:mi> <mml:mi mathvariant="normal"> n</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>代表您和密度取决于(<gydF4y2Baxref ref-type="disp-formula" rid="EEq20"> 19</gydF4y2Baxref>)- (<gydF4y2Baxref ref-type="disp-formula" rid="EEq22"> 21</gydF4y2Baxref>)。这两个可能被视为解决输入数据(<gydF4y2Baxref ref-type="disp-formula" rid="EEq25"> 24</gydF4y2Baxref>)。<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M126"> <mml:mi> E</米米l:mi> <mml:mfenced separators="|"> <mml:mrow> <mml:mi> T</米米l:mi> <mml:mo> ,</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> ρ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi mathvariant="normal"> 我</米米l:mi> <mml:mi mathvariant="normal"> n</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:mfenced> </mml:math> </inline-formula>代表由计算结果(<gydF4y2Baxref ref-type="disp-formula" rid="EEq23"> 22</gydF4y2Baxref>)。解决的过程(<gydF4y2Baxref ref-type="disp-formula" rid="EEq25"> 24</gydF4y2Baxref>)如下:<l是t> <list-item> <label>(一)</l一个bel> </list-item> </list></p> <p>集<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M127"> <mml:mrow> <mml:mi> T</米米l:mi> </mml:mrow> </mml:math> </inline-formula>迭代初始值。</pgydF4y2Ba> <list-item> <label>(B)</l一个bel> <p>解决(<gydF4y2Baxref ref-type="disp-formula" rid="EEq19"> 18</gydF4y2Baxref>)程序(一)——(d)获得<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M128"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> j</米米l:mi> <mml:mi> 我</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>。</pgydF4y2Ba> </list-item> <list-item> <label>(C)</l一个bel> <p>计算(<gydF4y2Baxref ref-type="disp-formula" rid="EEq23"> 22</gydF4y2Baxref>)获得<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M129"> <mml:mi> E</米米l:mi> <mml:mfenced separators="|"> <mml:mrow> <mml:mi> T</米米l:mi> <mml:mo> ,</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> ρ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi mathvariant="normal"> 我</米米l:mi> <mml:mi mathvariant="normal"> n</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:mfenced> </mml:math> </inline-formula>。</pgydF4y2Ba> </list-item> <list-item> <label>(D)</l一个bel> <p>计算(<gydF4y2Baxref ref-type="disp-formula" rid="EEq25"> 24</gydF4y2Baxref>)获得一个新的<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M130"> <mml:mrow> <mml:mi> T</米米l:mi> </mml:mrow> </mml:math> </inline-formula>表示为<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M131"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> T</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi mathvariant="normal"> r</米米l:mi> <mml:mi mathvariant="normal"> e</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>。</pgydF4y2Ba> </list-item> <list-item> <label>(E)</l一个bel> <p>集<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M132"> <mml:mrow> <mml:mi> T</米米l:mi> </mml:mrow> </mml:math> </inline-formula>是<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M133"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> T</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi mathvariant="normal"> r</米米l:mi> <mml:mi mathvariant="normal"> e</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>和重复步骤(B) - (C),直到满足终止条件。</pgydF4y2Ba> </list-item> <p></p> <p>通过上述方法,相对困难的多维问题解决萨哈方程模拟的毛细管放电转化为一维问题<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M134"> <mml:mrow> <mml:mi> T</米米l:mi> </mml:mrow> </mml:math> </inline-formula>。在(<gydF4y2Baxref ref-type="disp-formula" rid="EEq25"> 24</gydF4y2Baxref>),一个系数<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M135"> <mml:mrow> <mml:mi> δ</米米l:mi> </mml:mrow> </mml:math> </inline-formula>介绍了。它是用来使迭代方程(<gydF4y2Baxref ref-type="disp-formula" rid="EEq25"> 24</gydF4y2Baxref>)满足收敛条件。(<gydF4y2Baxref ref-type="disp-formula" rid="EEq25"> 24</gydF4y2Baxref>),这种收敛条件可以表示为(<gydF4y2Baxref ref-type="bibr" rid="B20"> 20.</gydF4y2Baxref>]<d我年代p- - - - - -for米ul一个> <mml:math display="block" xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M136"> <mml:mtable> <mml:mlabeledtr id="EEq26"> <mml:mtd> <mml:mtext> (25)</米米l:mtext> </mml:mtd> <mml:mtd> <mml:mfenced open="|" close="|" separators="|"> <mml:mrow> <mml:msup> <mml:mrow> <mml:mi> h</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> ′</米米l:mi> </mml:mrow> </mml:msup> <mml:mfenced separators="|"> <mml:mrow> <mml:mi> T</米米l:mi> </mml:mrow> </mml:mfenced> </mml:mrow> </mml:mfenced> <mml:mo> <</米米l:mo> <mml:mn mathvariant="normal"> 1</米米l:mn> <mml:mo> ,</米米l:mo> </mml:mtd> </mml:mlabeledtr> </mml:mtable> </mml:math> </disp-formula>在哪里<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M137"> <mml:mi> h</米米l:mi> <mml:mo stretchy="false"> (</米米l:mo> <mml:mi> T</米米l:mi> <mml:mo stretchy="false"> )</米米l:mo> </mml:math> </inline-formula>代表的右边(<gydF4y2Baxref ref-type="disp-formula" rid="EEq25"> 24</gydF4y2Baxref>)。不等式(<gydF4y2Baxref ref-type="disp-formula" rid="EEq26"> 25</gydF4y2Baxref>),一个可以获得的范围<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M138"> <mml:mrow> <mml:mi> δ</米米l:mi> </mml:mrow> </mml:math> </inline-formula>,这是<d我年代p- - - - - -for米ul一个> <mml:math display="block" xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M139"> <mml:mtable> <mml:mlabeledtr id="EEq27"> <mml:mtd> <mml:mtext> (26)</米米l:mtext> </mml:mtd> <mml:mtd> <mml:mn mathvariant="normal"> 0</米米l:mn> <mml:mo> <</米米l:mo> <mml:mi> δ</米米l:mi> <mml:mo> <</米米l:mo> <mml:mfrac> <mml:mrow> <mml:mn mathvariant="normal"> 2</米米l:mn> </mml:mrow> <mml:mrow> <mml:mrow> <mml:mrow> <mml:mo> ∂</米米l:mo> <mml:mi> E</米米l:mi> </mml:mrow> <mml:mo> /</米米l:mo> <mml:mrow> <mml:mo> ∂</米米l:mo> <mml:mi> T</米米l:mi> </mml:mrow> </mml:mrow> </mml:mrow> </mml:mfrac> <mml:mo> 。</米米l:mo> </mml:mtd> </mml:mlabeledtr> </mml:mtable> </mml:math> </disp-formula>在毛细管放电或等发射技术,的大小顺序<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M140"> <mml:mo> ∂</米米l:mo> <mml:mi> E</米米l:mi> <mml:mo> /</米米l:mo> <mml:mo> ∂</米米l:mo> <mml:mi> T</米米l:mi> </mml:math> </inline-formula>通常是<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M141"> <mml:mn> 1</米米l:mn> <mml:mo> ×</米米l:mo> <mml:msup> <mml:mrow> <mml:mn> 10</米米l:mn> </mml:mrow> <mml:mrow> <mml:mn> 4</米米l:mn> </mml:mrow> </mml:msup> </mml:math> </inline-formula>(国际标准单位)。很明显,(<gydF4y2Baxref ref-type="disp-formula" rid="EEq25"> 24</gydF4y2Baxref>)如果不能满足收敛条件<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M142"> <mml:mrow> <mml:mi> δ</米米l:mi> </mml:mrow> </mml:math> </inline-formula>不是介绍或设置为1。的特定值<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M143"> <mml:mrow> <mml:mi> δ</米米l:mi> </mml:mrow> </mml:math> </inline-formula>在一个细胞,一个可以设置为最大的计算效率。不等式(<gydF4y2Baxref ref-type="disp-formula" rid="EEq26"> 25</gydF4y2Baxref>),可以发现,如果<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M144"> <mml:mrow> <mml:mi> δ</米米l:mi> </mml:mrow> </mml:math> </inline-formula>使<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M145"> <mml:msup> <mml:mrow> <mml:mi> h</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> ′</米米l:mi> </mml:mrow> </mml:msup> <mml:mfenced separators="|"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> T</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> t</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:mfenced> <mml:mo> =</米米l:mo> <mml:mn> 0</米米l:mn> </mml:math> </inline-formula>,在那里<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M146"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> T</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> t</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>代表了真正的根细胞中,它可能需要至少迭代找到这个细胞的根。也就是说,<d我年代p- - - - - -for米ul一个> <mml:math display="block" xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M147"> <mml:mtable> <mml:mlabeledtr id="EEq28"> <mml:mtd> <mml:mtext> (27)</米米l:mtext> </mml:mtd> <mml:mtd> <mml:mi> δ</米米l:mi> <mml:mo> =</米米l:mo> <mml:msubsup> <mml:mrow> <mml:mfenced separators="|"> <mml:mrow> <mml:mfrac> <mml:mrow> <mml:mo> ∂</米米l:mo> <mml:mi> E</米米l:mi> </mml:mrow> <mml:mrow> <mml:mo> ∂</米米l:mo> <mml:mi> T</米米l:mi> </mml:mrow> </mml:mfrac> </mml:mrow> </mml:mfenced> </mml:mrow> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> T</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> t</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> <mml:mrow> <mml:mo> - - - - - -</米米l:mo> <mml:mn mathvariant="normal"> 1</米米l:mn> </mml:mrow> </mml:msubsup> <mml:mo> 。</米米l:mo> </mml:mtd> </mml:mlabeledtr> </mml:mtable> </mml:math> </disp-formula>理想的情况是,如果<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M148"> <mml:mi> h</米米l:mi> <mml:mo stretchy="false"> (</米米l:mo> <mml:mi> T</米米l:mi> <mml:mo stretchy="false"> )</米米l:mo> </mml:math> </inline-formula>是一个线性函数,<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M149"> <mml:mrow> <mml:mi> δ</米米l:mi> </mml:mrow> </mml:math> </inline-formula>由(<gydF4y2Baxref ref-type="disp-formula" rid="EEq28"> 27</gydF4y2Baxref>),只有一次的迭代将需要找到解决方案和收敛速度将是无穷大。对于非线性函数,如果<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M150"> <mml:mrow> <mml:mi> δ</米米l:mi> </mml:mrow> </mml:math> </inline-formula>使<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M151"> <mml:msup> <mml:mrow> <mml:mi> h</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> ′</米米l:mi> </mml:mrow> </mml:msup> <mml:mfenced separators="|"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> T</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> t</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:mfenced> <mml:mo> =</米米l:mo> <mml:mn> 0</米米l:mn> </mml:math> </inline-formula>在一个细胞,这意味着,<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M152"> <mml:mrow> <mml:mi> T</米米l:mi> </mml:mrow> </mml:math> </inline-formula>接近<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M153"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> T</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> t</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>,<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M154"> <mml:msup> <mml:mrow> <mml:mi> h</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> ′</米米l:mi> </mml:mrow> </mml:msup> <mml:mfenced separators="|"> <mml:mrow> <mml:mi> T</米米l:mi> </mml:mrow> </mml:mfenced> </mml:math> </inline-formula>会越来越接近<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M155"> <mml:msup> <mml:mrow> <mml:mi> h</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> ′</米米l:mi> </mml:mrow> </mml:msup> <mml:mfenced separators="|"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> T</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> t</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:mfenced> </mml:math> </inline-formula>0,这个过程的收敛速度也会增加。它更像是一个加速收敛速度的影响。与牛顿迭代法收敛的平方,简单的定点迭代方法本身是线性收敛的<gydF4y2Baxref ref-type="bibr" rid="B20"> 20.</gydF4y2Baxref>),但仍可能获得一些通过优化加速度的影响<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M156"> <mml:mrow> <mml:mi> δ</米米l:mi> </mml:mrow> </mml:math> </inline-formula>。也就是说,如果<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M157"> <mml:mrow> <mml:mi> δ</米米l:mi> </mml:mrow> </mml:math> </inline-formula>将更接近<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M158"> <mml:mrow> <mml:msup> <mml:mrow> <mml:mfenced separators="|"> <mml:mrow> <mml:mrow> <mml:mrow> <mml:mo> ∂</米米l:mo> <mml:mi> E</米米l:mi> </mml:mrow> <mml:mo> /</米米l:mo> <mml:mrow> <mml:mo> ∂</米米l:mo> <mml:mi> T</米米l:mi> </mml:mrow> </mml:mrow> </mml:mrow> </mml:mfenced> </mml:mrow> <mml:mrow> <mml:mo> - - - - - -</米米l:mo> <mml:mn> 1</米米l:mn> </mml:mrow> </mml:msup> </mml:mrow> </mml:math> </inline-formula>细胞,需要更少的迭代找到根。</pgydF4y2Ba> <p>在毛细管放电的实际模拟,将是乏味的和不必要的<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M159"> <mml:mrow> <mml:mi> δ</米米l:mi> </mml:mrow> </mml:math> </inline-formula>每个细胞在不同的时间点。一个更可行的方法是设置<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M160"> <mml:mrow> <mml:mi> δ</米米l:mi> </mml:mrow> </mml:math> </inline-formula>是一个统一的价值对整个模拟只要可以满足收敛条件。基于上述分析,我们可以发现,最小运行时,这个统一<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M161"> <mml:mrow> <mml:mi> δ</米米l:mi> </mml:mrow> </mml:math> </inline-formula>应选用非常接近的平均值<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M162"> <mml:mrow> <mml:msup> <mml:mrow> <mml:mfenced separators="|"> <mml:mrow> <mml:mrow> <mml:mrow> <mml:mo> ∂</米米l:mo> <mml:mi> E</米米l:mi> </mml:mrow> <mml:mo> /</米米l:mo> <mml:mrow> <mml:mo> ∂</米米l:mo> <mml:mi> T</米米l:mi> </mml:mrow> </mml:mrow> </mml:mrow> </mml:mfenced> </mml:mrow> <mml:mrow> <mml:mo> - - - - - -</米米l:mo> <mml:mn> 1</米米l:mn> </mml:mrow> </mml:msup> </mml:mrow> </mml:math> </inline-formula>在所有的细胞在每一个时间点。因为一个人不能获得<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M163"> <mml:mrow> <mml:msup> <mml:mrow> <mml:mfenced separators="|"> <mml:mrow> <mml:mrow> <mml:mrow> <mml:mo> ∂</米米l:mo> <mml:mi> E</米米l:mi> </mml:mrow> <mml:mo> /</米米l:mo> <mml:mrow> <mml:mo> ∂</米米l:mo> <mml:mi> T</米米l:mi> </mml:mrow> </mml:mrow> </mml:mrow> </mml:mfenced> </mml:mrow> <mml:mrow> <mml:mo> - - - - - -</米米l:mo> <mml:mn> 1</米米l:mn> </mml:mrow> </mml:msup> </mml:mrow> </mml:math> </inline-formula>事先也不平均,<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M164"> <mml:mrow> <mml:mi> δ</米米l:mi> </mml:mrow> </mml:math> </inline-formula>可以将一个相对较小的值,以确保收敛一看到从不等式(<gydF4y2Baxref ref-type="disp-formula" rid="EEq27"> 26</gydF4y2Baxref>)第一次模拟;然后,一个人可以获得的平均值<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M165"> <mml:mrow> <mml:msup> <mml:mrow> <mml:mfenced separators="|"> <mml:mrow> <mml:mrow> <mml:mrow> <mml:mo> ∂</米米l:mo> <mml:mi> E</米米l:mi> </mml:mrow> <mml:mo> /</米米l:mo> <mml:mrow> <mml:mo> ∂</米米l:mo> <mml:mi> T</米米l:mi> </mml:mrow> </mml:mrow> </mml:mrow> </mml:mfenced> </mml:mrow> <mml:mrow> <mml:mo> - - - - - -</米米l:mo> <mml:mn> 1</米米l:mn> </mml:mrow> </mml:msup> </mml:mrow> </mml:math> </inline-formula>从仿真结果。自<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M166"> <mml:mrow> <mml:mi> δ</米米l:mi> </mml:mrow> </mml:math> </inline-formula>设置为一个统一的价值,确切的价值吗<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M167"> <mml:mrow> <mml:mi> δ</米米l:mi> </mml:mrow> </mml:math> </inline-formula>最低运行时也会受到其他因素的影响,如设置迭代初始值的方法。一个可能需要调整<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M168"> <mml:mrow> <mml:mi> δ</米米l:mi> </mml:mrow> </mml:math> </inline-formula>在此基础上的平均值<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M169"> <mml:mrow> <mml:msup> <mml:mrow> <mml:mfenced separators="|"> <mml:mrow> <mml:mrow> <mml:mrow> <mml:mo> ∂</米米l:mo> <mml:mi> E</米米l:mi> </mml:mrow> <mml:mo> /</米米l:mo> <mml:mrow> <mml:mo> ∂</米米l:mo> <mml:mi> T</米米l:mi> </mml:mrow> </mml:mrow> </mml:mrow> </mml:mfenced> </mml:mrow> <mml:mrow> <mml:mo> - - - - - -</米米l:mo> <mml:mn> 1</米米l:mn> </mml:mrow> </mml:msup> </mml:mrow> </mml:math> </inline-formula>额外的模拟测试,直到最低模拟运行时实现。</pgydF4y2Ba> </sec> <sec id="sec5"> <title>5。设置迭代初始值</t我tle> <p>设置迭代初始值的基本原理是把它尽可能接近真实的根,所以需要更少的迭代来找到它。此外,自<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M170"> <mml:mi> h</米米l:mi> <mml:mo stretchy="false"> (</米米l:mo> <mml:mi> T</米米l:mi> <mml:mo stretchy="false"> )</米米l:mo> </mml:math> </inline-formula>是非线性函数,上面的方法设置<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M171"> <mml:mrow> <mml:mi> δ</米米l:mi> </mml:mrow> </mml:math> </inline-formula>也表明,如果迭代初始值太远离真正的根,迭代方程不能满足收敛条件。时间模拟的毛细管放电,等离子体流属性,包括<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M172"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> e</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>,<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M173"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 一个</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>,<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M174"> <mml:mrow> <mml:mi> T</米米l:mi> </mml:mrow> </mml:math> </inline-formula>,是动态的时空数据,所以可能还需要设定的迭代初始值<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M175"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> e</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>,<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M176"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 一个</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>,<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M177"> <mml:mrow> <mml:mi> T</米米l:mi> </mml:mrow> </mml:math> </inline-formula>动态。否则,输入电流的变化,<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M178"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> e</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>,<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M179"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 一个</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>,或<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M180"> <mml:mrow> <mml:mi> T</米米l:mi> </mml:mrow> </mml:math> </inline-formula>目前在一个细胞时间点可能会改变在很大程度上从这些属性这个细胞在以前的时间点的。这也可能导致大量的迭代找到解决方案,甚至使迭代发散的。</pgydF4y2Ba> <p>在毛细管放电仿真,以达到一个高质量的解决方案或满足Courant-Friedrichs-Lewy (CFL)条件下保持稳定的模拟,空间时间步步骤通常是将相对较小的值,所以等离子体特性的两个空间或时间相邻细胞可能非常接近对方。这实际上提供了两种简单的方法来设置动态迭代初始值。如图<gydF4y2Baxref ref-type="fig" rid="fig1"> 1</gydF4y2Baxref>为解决萨哈方程或(<gydF4y2Baxref ref-type="disp-formula" rid="EEq25"> 24</gydF4y2Baxref>)“现在”细胞,可以设定的迭代初始值<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M181"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> e</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>,<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M182"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 一个</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>,或<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M183"> <mml:mrow> <mml:mi> T</米米l:mi> </mml:mrow> </mml:math> </inline-formula>的解决方案(<gydF4y2Baxref ref-type="disp-formula" rid="EEq25"> 24</gydF4y2Baxref>)相应的“最后一次一步”或“最后的细胞。”在这篇文章中,为了方便,设置迭代初始值的解决方案“上次一步”被命名为“Initial-T”和设置迭代初始值的解决方案“最后细胞”命名为“Initial-S”;由于没有“最后细胞”的细胞封闭式的毛细管,“Initial-T”仍然是用于此“Initial-S细胞的方法之一。“另外,比较研究,设置迭代初始值固定命名为“Initial-F。“在这篇文章中,固定的迭代初始值<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M184"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> e</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>,<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M185"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 一个</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>,<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M186"> <mml:mrow> <mml:mi> T</米米l:mi> </mml:mrow> </mml:math> </inline-formula>将他们的平均超过所有的细胞在每一个时间点从先前的模拟结果。</pgydF4y2Ba> <fig id="fig1"> <label>图1</l一个bel> <p>原理图的设置迭代初始值的方法。</pgydF4y2Ba> <graphic xlink:href="//www.newsama.com/downloads/journals/amp/2017/6473234.fig.001"></graphic> </fig> </sec> <sec id="sec6"> <title>6。结果与讨论</t我tle> <p>一些数值试验进行了研究提出了性能的方法。本方法的准确性,进行了两个试验。第一个是程序(一)——(d)。在这个测试中,萨哈方程或(<gydF4y2Baxref ref-type="disp-formula" rid="EEq19"> 18</gydF4y2Baxref>)聚乙烯等离子体单独解决采用程序(一)——(d)与不同的压力和温度。摘要终端标准解决萨哈方程构造的近似相对误差百分比<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M187"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ε</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 一个</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>(<gydF4y2Baxref ref-type="bibr" rid="B20"> 20.</gydF4y2Baxref>]。为<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M188"> <mml:mrow> <mml:mi> T</米米l:mi> </mml:mrow> </mml:math> </inline-formula>、终端标准可以表示为<d我年代p- - - - - -for米ul一个> <mml:math display="block" xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M189"> <mml:mtable> <mml:mlabeledtr id="EEq29"> <mml:mtd> <mml:mtext> (28)</米米l:mtext> </mml:mtd> <mml:mtd> <mml:msub> <mml:mrow> <mml:mi> ε</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 一个</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> =</米米l:mo> <mml:mfenced open="|" close="|" separators="|"> <mml:mrow> <mml:mfrac> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> T</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi mathvariant="normal"> c</米米l:mi> <mml:mi mathvariant="normal"> u</米米l:mi> <mml:mi mathvariant="normal"> r</米米l:mi> <mml:mi mathvariant="normal"> r</米米l:mi> <mml:mi mathvariant="normal"> e</米米l:mi> <mml:mi mathvariant="normal"> n</米米l:mi> <mml:mi mathvariant="normal"> t</米米l:mi> </mml:mrow> </mml:msub> <mml:mo> - - - - - -</米米l:mo> <mml:msub> <mml:mrow> <mml:mi> T</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi mathvariant="normal"> p</米米l:mi> <mml:mi mathvariant="normal"> r</米米l:mi> <mml:mi mathvariant="normal"> e</米米l:mi> <mml:mi mathvariant="normal"> v</米米l:mi> <mml:mi mathvariant="normal"> 我</米米l:mi> <mml:mi mathvariant="normal"> o</米米l:mi> <mml:mi mathvariant="normal"> u</米米l:mi> <mml:mi mathvariant="normal"> 年代</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> T</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi mathvariant="normal"> c</米米l:mi> <mml:mi mathvariant="normal"> u</米米l:mi> <mml:mi mathvariant="normal"> r</米米l:mi> <mml:mi mathvariant="normal"> r</米米l:mi> <mml:mi mathvariant="normal"> e</米米l:mi> <mml:mi mathvariant="normal"> n</米米l:mi> <mml:mi mathvariant="normal"> t</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:mfrac> </mml:mrow> </mml:mfenced> <mml:mo> <</米米l:mo> <mml:mn mathvariant="normal"> 5</米米l:mn> <mml:mo> ×</米米l:mo> <mml:mn mathvariant="normal"> 1</米米l:mn> <mml:msup> <mml:mrow> <mml:mn mathvariant="normal"> 0</米米l:mn> </mml:mrow> <mml:mrow> <mml:mo> - - - - - -</米米l:mo> <mml:mn mathvariant="normal"> 6</米米l:mn> </mml:mrow> </mml:msup> <mml:mo> 。</米米l:mo> </mml:mtd> </mml:mlabeledtr> </mml:mtable> </mml:math> </disp-formula></p> <p>通过这个终端标准,5位有效数字可以实现,这是为目前的研究发现是足够精确的。然后,解决方案,电导率计算并与金正日的结果(<gydF4y2Baxref ref-type="bibr" rid="B11"> 11</gydF4y2Baxref>]。在图<gydF4y2Baxref ref-type="fig" rid="fig2"> 2</gydF4y2Baxref>,聚乙烯导电率计算等离子体温度和压力的函数从现在的结果,<gydF4y2Baxref ref-type="bibr" rid="B11"> 11</gydF4y2Baxref>绘制。可以观察到,现在的结果是发现金完全同意的结果(<gydF4y2Baxref ref-type="bibr" rid="B11"> 11</gydF4y2Baxref>)这可能证明程序的准确性(一)——(d)以及编码工作的正确性。</pgydF4y2Ba> <fig id="fig2"> <label>图2</l一个bel> <p>计算为聚乙烯等离子体电导率温度和压力的函数。实线代表目前的计算结果和符号代表的结果(<gydF4y2Baxref ref-type="bibr" rid="B11"> 11</gydF4y2Baxref>]。</pgydF4y2Ba> <graphic xlink:href="//www.newsama.com/downloads/journals/amp/2017/6473234.fig.002"></graphic> </fig> <p>第二个测试是程序(一)- (D)。在该测试中,使用的毛细管放电等或等发射系统数值模拟与实验数据从投8中系列实验鲍威尔和Zielinski [<gydF4y2Baxref ref-type="bibr" rid="B4"> 4</gydF4y2Baxref>]和程序(一)——(D)是用来解决萨哈方程在每个单元模拟。同时,模拟等离子体流属性提出的模拟结果相比,在鲍威尔和Zielinski [<gydF4y2Baxref ref-type="bibr" rid="B4"> 4</gydF4y2Baxref>]。在图<gydF4y2Baxref ref-type="fig" rid="fig3"> 3</gydF4y2Baxref>,输入电流和使用的毛细管的几何模拟绘制。输入当前数据图<gydF4y2Baxref ref-type="fig" rid="fig3a"> 3(一个)</gydF4y2Baxref>收集从[<gydF4y2Baxref ref-type="bibr" rid="B11"> 11</gydF4y2Baxref>)与图像识别软件,然后安装到ninth-order多项式。从200年开始进行仿真<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M190"> <mml:mrow> <mml:mi> μ</米米l:mi> </mml:mrow> </mml:math> </inline-formula>年代到700<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M191"> <mml:mrow> <mml:mi> μ</米米l:mi> </mml:mrow> </mml:math> </inline-formula>年代避免初期由电气爆炸用来初始化弧(<gydF4y2Baxref ref-type="bibr" rid="B4"> 4</gydF4y2Baxref>,<gydF4y2Baxref ref-type="bibr" rid="B11"> 11</gydF4y2Baxref>]。用于仿真的数学模型(<gydF4y2Baxref ref-type="disp-formula" rid="EEq20"> 19</gydF4y2Baxref>)- (<gydF4y2Baxref ref-type="disp-formula" rid="EEq24"> 23</gydF4y2Baxref>)。宽松方案求解这个模型。网格点和CFL数是固定的100和0.5,分别。设置迭代初始值的方法是“Initial-T。“模拟结果与“Initial-S”或“Initial-F”几乎可以互相区别。在图<gydF4y2Baxref ref-type="fig" rid="fig4"> 4</gydF4y2Baxref>沿轴向方向,模拟等离子体流属性为500<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M192"> <mml:mrow> <mml:mi> μ</米米l:mi> </mml:mrow> </mml:math> </inline-formula>s是策划。可以观察到的,从目前的仿真模拟等离子体流属性略高于鲍威尔和Zielinski [<gydF4y2Baxref ref-type="bibr" rid="B4"> 4</gydF4y2Baxref>]。可能的原因之一可能是,当前的输入电流收集与图像识别软件不是一模一样,用于模拟鲍威尔和Zielinski。其他可能的原因这些细微的差别可以设置不同的网格点或时间步骤,用于模拟不同方案,等等。更多详细信息的模拟和讨论模拟等离子体流的属性,可以发现他们在<gydF4y2Baxref ref-type="bibr" rid="B4"> 4</gydF4y2Baxref>,<gydF4y2Baxref ref-type="bibr" rid="B11"> 11</gydF4y2Baxref>,<gydF4y2Baxref ref-type="bibr" rid="B21"> 21</gydF4y2Baxref>]。</pgydF4y2Ba> <fig-group id="fig3"> <label>图3</l一个bel> <p>输入电流分布和几何的毛细管放电<gydF4y2Baxref ref-type="bibr" rid="B4"> 4</gydF4y2Baxref>]。</pgydF4y2Ba> <fig id="fig3a"> <label>(一)</l一个bel> <graphic xlink:href="//www.newsama.com/downloads/journals/amp/2017/6473234.fig.003a"></graphic> </fig> <fig id="fig3b"> <label>(b)</l一个bel> <graphic xlink:href="//www.newsama.com/downloads/journals/amp/2017/6473234.fig.003b"></graphic> </fig> </fig-group> <fig id="fig4"> <label>图4</l一个bel> <p>模拟等离子体流沿轴向方向属性为500<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M193"> <mml:mrow> <mml:mi> µ</米米l:mi> </mml:mrow> </mml:math> </inline-formula>年代。<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M194"> <mml:mrow> <mml:mi> l</米米l:mi> </mml:mrow> </mml:math> </inline-formula>计算区域的总长度(<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M195"> <mml:mi> l</米米l:mi> <mml:mo> =</米米l:mo> <mml:mn> 60.9</米米l:mn> </mml:math> </inline-formula>毫米)。实线代表目前的计算结果和符号代表的结果(<gydF4y2Baxref ref-type="bibr" rid="B4"> 4</gydF4y2Baxref>]。</pgydF4y2Ba> <graphic xlink:href="//www.newsama.com/downloads/journals/amp/2017/6473234.fig.004"></graphic> </fig> <p>基于上面的模拟,目前方法的效率进行了研究。在图<gydF4y2Baxref ref-type="fig" rid="fig5"> 5</gydF4y2Baxref>,模拟运行时<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M196"> <mml:mrow> <mml:mi> δ</米米l:mi> </mml:mrow> </mml:math> </inline-formula>变化从0.4到0.8,不同设置迭代初始值的方法绘制以选择<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M197"> <mml:mrow> <mml:mi> δ</米米l:mi> </mml:mrow> </mml:math> </inline-formula>最低的运行时。用c++编写模拟程序和Windows 7的CPU上运行英特尔酷睿i7 - 6700总部。大部分的文件编写函数关闭模拟程序的效率。每个最终记录运行时重复模拟运行时平均超过50倍。可以观察到,有一个显著的优势利用“Initial-S”或在“Initial-T Initial-F。“运行时”Initial-T”和“Initial-S”大约是安排从9到13秒20至29秒,分别而运行时以“Initial-F”大约是安排从45到69秒。在图<gydF4y2Baxref ref-type="fig" rid="fig5"> 5</gydF4y2Baxref>还发现,即使<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M198"> <mml:mrow> <mml:mi> δ</米米l:mi> </mml:mrow> </mml:math> </inline-formula>将<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M199"> <mml:mn> 0.7</米米l:mn> <mml:mo> ×</米米l:mo> <mml:msup> <mml:mrow> <mml:mn> 10</米米l:mn> </mml:mrow> <mml:mrow> <mml:mo> - - - - - -</米米l:mo> <mml:mn> 4</米米l:mn> </mml:mrow> </mml:msup> </mml:math> </inline-formula>“Initial-T”总运行时仍然是大约50%不到的最小运行时用“Initial-S仿真。“方法”Initial-S”和“Initial-T,”<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M200"> <mml:mfenced open="|" close="|" separators="|"> <mml:mrow> <mml:mfenced separators="|"> <mml:mrow> <mml:mo> ∂</米米l:mo> <mml:mi> ζ</米米l:mi> <mml:mo> /</米米l:mo> <mml:mo> ∂</米米l:mo> <mml:mi> z</米米l:mi> </mml:mrow> </mml:mfenced> </mml:mrow> </mml:mfenced> <mml:mo> ×</米米l:mo> <mml:mi> Δ</米米l:mi> <mml:mi> z</米米l:mi> </mml:math> </inline-formula>和<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M201"> <mml:mfenced open="|" close="|" separators="|"> <mml:mrow> <mml:mfenced separators="|"> <mml:mrow> <mml:mo> ∂</米米l:mo> <mml:mi> ζ</米米l:mi> <mml:mo> /</米米l:mo> <mml:mo> ∂</米米l:mo> <mml:mi> t</米米l:mi> </mml:mrow> </mml:mfenced> </mml:mrow> </mml:mfenced> <mml:mo> ×</米米l:mo> <mml:mi> Δ</米米l:mi> <mml:mi> t</米米l:mi> </mml:math> </inline-formula>,在那里<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M202"> <mml:mrow> <mml:mi> ζ</米米l:mi> </mml:mrow> </mml:math> </inline-formula>可能是<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M203"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> e</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>,<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M204"> <mml:mrow> <mml:msub> <mml:mrow> <mml:mi> ψ</米米l:mi> </mml:mrow> <mml:mrow> <mml:mi> 一个</米米l:mi> </mml:mrow> </mml:msub> </mml:mrow> </mml:math> </inline-formula>,<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M205"> <mml:mrow> <mml:mi> T</米米l:mi> </mml:mrow> </mml:math> </inline-formula>或其他属性,实际上意味着迭代初始值的距离真正的根,所以这两种方法连接<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M206"> <mml:mfenced open="|" close="|" separators="|"> <mml:mrow> <mml:mfenced separators="|"> <mml:mrow> <mml:mo> ∂</米米l:mo> <mml:mi> ζ</米米l:mi> <mml:mo> /</米米l:mo> <mml:mo> ∂</米米l:mo> <mml:mi> z</米米l:mi> </mml:mrow> </mml:mfenced> </mml:mrow> </mml:mfenced> <mml:mo> ×</米米l:mo> <mml:mi> Δ</米米l:mi> <mml:mi> z</米米l:mi> </mml:math> </inline-formula>和<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M207"> <mml:mfenced open="|" close="|" separators="|"> <mml:mrow> <mml:mfenced separators="|"> <mml:mrow> <mml:mo> ∂</米米l:mo> <mml:mi> ζ</米米l:mi> <mml:mo> /</米米l:mo> <mml:mo> ∂</米米l:mo> <mml:mi> t</米米l:mi> </mml:mrow> </mml:mfenced> </mml:mrow> </mml:mfenced> <mml:mo> ×</米米l:mo> <mml:mi> Δ</米米l:mi> <mml:mi> t</米米l:mi> </mml:math> </inline-formula>计算效率。它可能是<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M208"> <mml:mfenced open="|" close="|" separators="|"> <mml:mrow> <mml:mfenced separators="|"> <mml:mrow> <mml:mo> ∂</米米l:mo> <mml:mi> ζ</米米l:mi> <mml:mo> /</米米l:mo> <mml:mo> ∂</米米l:mo> <mml:mi> z</米米l:mi> </mml:mrow> </mml:mfenced> </mml:mrow> </mml:mfenced> <mml:mo> ×</米米l:mo> <mml:mi> Δ</米米l:mi> <mml:mi> z</米米l:mi> </mml:math> </inline-formula>通常是大于<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M209"> <mml:mfenced open="|" close="|" separators="|"> <mml:mrow> <mml:mfenced separators="|"> <mml:mrow> <mml:mo> ∂</米米l:mo> <mml:mi> ζ</米米l:mi> <mml:mo> /</米米l:mo> <mml:mo> ∂</米米l:mo> <mml:mi> t</米米l:mi> </mml:mrow> </mml:mfenced> </mml:mrow> </mml:mfenced> <mml:mo> ×</米米l:mo> <mml:mi> Δ</米米l:mi> <mml:mi> t</米米l:mi> </mml:math> </inline-formula>。从仿真结果,发现的大小顺序<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M210"> <mml:mrow> <mml:mfenced open="|" close="|" separators="|"> <mml:mrow> <mml:mo> ∂</米米l:mo> <mml:mi> T</米米l:mi> <mml:mo> /</米米l:mo> <mml:mo> ∂</米米l:mo> <mml:mi> t</米米l:mi> </mml:mrow> </mml:mfenced> </mml:mrow> </mml:math> </inline-formula>是<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M211"> <mml:mn> 1</米米l:mn> <mml:mo> ×</米米l:mo> <mml:msup> <mml:mrow> <mml:mn> 10</米米l:mn> </mml:mrow> <mml:mrow> <mml:mn> 7</米米l:mn> </mml:mrow> </mml:msup> </mml:math> </inline-formula>和的大小顺序<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M212"> <mml:mi> Δ</米米l:mi> <mml:mi> t</米米l:mi> </mml:math> </inline-formula>是<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M213"> <mml:mn> 1</米米l:mn> <mml:mo> ×</米米l:mo> <mml:msup> <mml:mrow> <mml:mn> 10</米米l:mn> </mml:mrow> <mml:mrow> <mml:mo> - - - - - -</米米l:mo> <mml:mn> 8</米米l:mn> </mml:mrow> </mml:msup> </mml:math> </inline-formula>,这样的大小顺序<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M214"> <mml:mfenced open="|" close="|" separators="|"> <mml:mrow> <mml:mo> ∂</米米l:mo> <mml:mi> T</米米l:mi> <mml:mo> /</米米l:mo> <mml:mo> ∂</米米l:mo> <mml:mi> t</米米l:mi> </mml:mrow> </mml:mfenced> <mml:mo> ×</米米l:mo> <mml:mi> Δ</米米l:mi> <mml:mi> t</米米l:mi> </mml:math> </inline-formula>大约是1 K。的大小顺序<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M215"> <mml:mrow> <mml:mfenced open="|" close="|" separators="|"> <mml:mrow> <mml:mo> ∂</米米l:mo> <mml:mi> T</米米l:mi> <mml:mo> /</米米l:mo> <mml:mo> ∂</米米l:mo> <mml:mi> z</米米l:mi> </mml:mrow> </mml:mfenced> </mml:mrow> </mml:math> </inline-formula>,它是不同的从0封闭式的毛细管<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M216"> <mml:mn> 3</米米l:mn> <mml:mo> ×</米米l:mo> <mml:msup> <mml:mrow> <mml:mn> 10</米米l:mn> </mml:mrow> <mml:mrow> <mml:mn> 5</米米l:mn> </mml:mrow> </mml:msup> </mml:math> </inline-formula>在出口处。自<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M217"> <mml:mi> Δ</米米l:mi> <mml:mi> z</米米l:mi> </mml:math> </inline-formula>是<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M218"> <mml:mn> 60.9</米米l:mn> <mml:mo> ×</米米l:mo> <mml:msup> <mml:mrow> <mml:mn> 10</米米l:mn> </mml:mrow> <mml:mrow> <mml:mo> - - - - - -</米米l:mo> <mml:mn> 5</米米l:mn> </mml:mrow> </mml:msup> </mml:math> </inline-formula>,<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M219"> <mml:mfenced open="|" close="|" separators="|"> <mml:mrow> <mml:mo> ∂</米米l:mo> <mml:mi> T</米米l:mi> <mml:mo> /</米米l:mo> <mml:mo> ∂</米米l:mo> <mml:mi> z</米米l:mi> </mml:mrow> </mml:mfenced> <mml:mo> ×</米米l:mo> <mml:mi> Δ</米米l:mi> <mml:mi> z</米米l:mi> </mml:math> </inline-formula>最多可能达到约180 K。</pgydF4y2Ba> <fig id="fig5"> <label>图5</l一个bel> <p>模拟运行时不同的<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M220"> <mml:mrow> <mml:mi> δ</米米l:mi> </mml:mrow> </mml:math> </inline-formula>和不同的迭代初始值的设置方法。</pgydF4y2Ba> <graphic xlink:href="//www.newsama.com/downloads/journals/amp/2017/6473234.fig.005"></graphic> </fig> <p>在图<gydF4y2Baxref ref-type="fig" rid="fig6"> 6</gydF4y2Baxref>、时空分布的循环计数程序(i) - (vi)从模拟“Initial-S”和“Initial-T”策划,分别。因为程序的结构(一)——(E)实际上是一个三层的循环和过程(i) - (vi)是最内层循环,循环计数的程序(i) - (vi)可以在一个细胞显示总迭代。的价值<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M221"> <mml:mrow> <mml:mi> δ</米米l:mi> </mml:mrow> </mml:math> </inline-formula>用于模拟图<gydF4y2Baxref ref-type="fig" rid="fig5"> 5</gydF4y2Baxref>最低运行时。在图<gydF4y2Baxref ref-type="fig" rid="fig6a"> 6(一)</gydF4y2Baxref>,因为<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M222"> <mml:mrow> <mml:mfenced open="|" close="|" separators="|"> <mml:mrow> <mml:mo> ∂</米米l:mo> <mml:mi> ζ</米米l:mi> <mml:mo> /</米米l:mo> <mml:mo> ∂</米米l:mo> <mml:mi> z</米米l:mi> </mml:mrow> </mml:mfenced> </mml:mrow> </mml:math> </inline-formula>沿轴向方向总是增加毛细管放电,<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M223"> <mml:mfenced open="|" close="|" separators="|"> <mml:mrow> <mml:mfenced separators="|"> <mml:mrow> <mml:mo> ∂</米米l:mo> <mml:mi> ζ</米米l:mi> <mml:mo> /</米米l:mo> <mml:mo> ∂</米米l:mo> <mml:mi> z</米米l:mi> </mml:mrow> </mml:mfenced> </mml:mrow> </mml:mfenced> <mml:mo> ×</米米l:mo> <mml:mi> Δ</米米l:mi> <mml:mi> z</米米l:mi> </mml:math> </inline-formula>将是更大的在该地区附近的毛细管的退出。这导致迭代初始值从“细胞”成为进一步获得真正的根沿着轴向方向和循环计数也会增加并达到最高水平毛细管的出口附近的地区。在图<gydF4y2Baxref ref-type="fig" rid="fig6b"> 6 (b)</gydF4y2Baxref>分布通常是一致的<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M224"> <mml:mo> ∂</米米l:mo> <mml:mi> 我</米米l:mi> <mml:mo> /</米米l:mo> <mml:mo> ∂</米米l:mo> <mml:mi> t</米米l:mi> </mml:math> </inline-formula>在这<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M225"> <mml:mrow> <mml:mi> 我</米米l:mi> </mml:mrow> </mml:math> </inline-formula>是输入电流。早期和晚期阶段的模拟,可以观察到当前的迅速增加和减少从图<gydF4y2Baxref ref-type="fig" rid="fig3a"> 3(一个)</gydF4y2Baxref>;也就是说,<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M226"> <mml:mrow> <mml:mfenced open="|" close="|" separators="|"> <mml:mrow> <mml:mo> ∂</米米l:mo> <mml:mi> 我</米米l:mi> <mml:mo> /</米米l:mo> <mml:mo> ∂</米米l:mo> <mml:mi> t</米米l:mi> </mml:mrow> </mml:mfenced> </mml:mrow> </mml:math> </inline-formula>是在这两个时期在一个相对较高的水平。这使得等离子体流属性也迅速改变,这样的迭代初始值获得“上次一步”可能相对远离真正的根在早期和后期的仿真。在430年的中期<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M227"> <mml:mrow> <mml:mi> μ</米米l:mi> </mml:mrow> </mml:math> </inline-formula>年代到530<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M228"> <mml:mrow> <mml:mi> μ</米米l:mi> </mml:mrow> </mml:math> </inline-formula>当电流峰值附近,<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M229"> <mml:mrow> <mml:mfenced open="|" close="|" separators="|"> <mml:mrow> <mml:mo> ∂</米米l:mo> <mml:mi> 我</米米l:mi> <mml:mo> /</米米l:mo> <mml:mo> ∂</米米l:mo> <mml:mi> t</米米l:mi> </mml:mrow> </mml:mfenced> </mml:mrow> </mml:math> </inline-formula>大约是0,所以等离子体流属性将在这段时间内保持相对稳定。这使得迭代初始值,获得一份礼物细胞从“上次一步”变得非常接近真正的根,所以中间循环计数可能达到的最低水平。</pgydF4y2Ba> <fig-group id="fig6"> <label>图6</l一个bel> <p>时空分布的循环计数程序(i) - (vi) (a)和“Initial-S”<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M230"> <mml:mi> δ</米米l:mi> <mml:mo> =</米米l:mo> <mml:mn> 0.7</米米l:mn> <mml:mo> ×</米米l:mo> <mml:msup> <mml:mrow> <mml:mn> 10</米米l:mn> </mml:mrow> <mml:mrow> <mml:mo> - - - - - -</米米l:mo> <mml:mn> 4</米米l:mn> </mml:mrow> </mml:msup> </mml:math> </inline-formula>。(b)和“Initial-T”<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M231"> <mml:mi> δ</米米l:mi> <mml:mo> =</米米l:mo> <mml:mn> 0.54</米米l:mn> <mml:mo> ×</米米l:mo> <mml:msup> <mml:mrow> <mml:mn> 10</米米l:mn> </mml:mrow> <mml:mrow> <mml:mo> - - - - - -</米米l:mo> <mml:mn> 4</米米l:mn> </mml:mrow> </mml:msup> </mml:math> </inline-formula>。</pgydF4y2Ba> <fig id="fig6a"> <label>(一)</l一个bel> <graphic xlink:href="//www.newsama.com/downloads/journals/amp/2017/6473234.fig.006a"></graphic> </fig> <fig id="fig6b"> <label>(b)</l一个bel> <graphic xlink:href="//www.newsama.com/downloads/journals/amp/2017/6473234.fig.006b"></graphic> </fig> </fig-group> <p>相似的数据<gydF4y2Baxref ref-type="fig" rid="fig6a"> 6(一)</gydF4y2Baxref>和<gydF4y2Baxref ref-type="fig" rid="fig6b"> 6 (b)</gydF4y2Baxref>是所有循环计数可以达到自己的巅峰毛细管的出口附近的开始模拟。如前所述,一个细胞,解决(<gydF4y2Baxref ref-type="disp-formula" rid="EEq25"> 24</gydF4y2Baxref>)能找到至少迭代通过设置<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M232"> <mml:mrow> <mml:mi> δ</米米l:mi> </mml:mrow> </mml:math> </inline-formula>尽可能的真实<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M233"> <mml:mrow> <mml:msup> <mml:mrow> <mml:mfenced separators="|"> <mml:mrow> <mml:mo> ∂</米米l:mo> <mml:mi> E</米米l:mi> <mml:mo> /</米米l:mo> <mml:mo> ∂</米米l:mo> <mml:mi> T</米米l:mi> </mml:mrow> </mml:mfenced> </mml:mrow> <mml:mrow> <mml:mo> - - - - - -</米米l:mo> <mml:mn> 1</米米l:mn> </mml:mrow> </mml:msup> </mml:mrow> </mml:math> </inline-formula>这个细胞。在这篇文章中,<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M234"> <mml:mrow> <mml:mi> δ</米米l:mi> </mml:mrow> </mml:math> </inline-formula>将一个统一的价值对整个模拟,它可能发生,统一<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M235"> <mml:mrow> <mml:mi> δ</米米l:mi> </mml:mrow> </mml:math> </inline-formula>相对远离他们是真的吗<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M236"> <mml:mrow> <mml:msup> <mml:mrow> <mml:mfenced separators="|"> <mml:mrow> <mml:mo> ∂</米米l:mo> <mml:mi> E</米米l:mi> <mml:mo> /</米米l:mo> <mml:mo> ∂</米米l:mo> <mml:mi> T</米米l:mi> </mml:mrow> </mml:mfenced> </mml:mrow> <mml:mrow> <mml:mo> - - - - - -</米米l:mo> <mml:mn> 1</米米l:mn> </mml:mrow> </mml:msup> </mml:mrow> </mml:math> </inline-formula>在一些时空区域。这将导致收敛速度可能会保持在一个相对低的水平在迭代过程中导致相对高水平的循环计数在这些时空区域。在图<gydF4y2Baxref ref-type="fig" rid="fig7"> 7</gydF4y2Baxref>的时空分布<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M237"> <mml:mo> ∂</米米l:mo> <mml:mi> E</米米l:mi> <mml:mo> /</米米l:mo> <mml:mo> ∂</米米l:mo> <mml:mi> T</米米l:mi> </mml:math> </inline-formula>是绘制。它是发现,<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M238"> <mml:mo> ∂</米米l:mo> <mml:mi> E</米米l:mi> <mml:mo> /</米米l:mo> <mml:mo> ∂</米米l:mo> <mml:mi> T</米米l:mi> </mml:math> </inline-formula>可能达到峰值约<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M239"> <mml:mn> 2.14</米米l:mn> <mml:mo> ×</米米l:mo> <mml:msup> <mml:mrow> <mml:mn> 10</米米l:mn> </mml:mrow> <mml:mrow> <mml:mn> 4</米米l:mn> </mml:mrow> </mml:msup> </mml:math> </inline-formula>在出口处的毛细管的开始模拟和大约达到最低水平<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M240"> <mml:mn> 1.3</米米l:mn> <mml:mo> ×</米米l:mo> <mml:msup> <mml:mrow> <mml:mn> 10</米米l:mn> </mml:mrow> <mml:mrow> <mml:mn> 4</米米l:mn> </mml:mrow> </mml:msup> </mml:math> </inline-formula>封闭式的毛细管430点附近<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M241"> <mml:mrow> <mml:mi> μ</米米l:mi> </mml:mrow> </mml:math> </inline-formula>年代到530<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M242"> <mml:mrow> <mml:mi> μ</米米l:mi> </mml:mrow> </mml:math> </inline-formula>当电流也达到高峰。从仿真结果,发现的平均值<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M243"> <mml:mrow> <mml:msup> <mml:mrow> <mml:mfenced separators="|"> <mml:mrow> <mml:mo> ∂</米米l:mo> <mml:mi> E</米米l:mi> <mml:mo> /</米米l:mo> <mml:mo> ∂</米米l:mo> <mml:mi> T</米米l:mi> </mml:mrow> </mml:mfenced> </mml:mrow> <mml:mrow> <mml:mo> - - - - - -</米米l:mo> <mml:mn> 1</米米l:mn> </mml:mrow> </mml:msup> </mml:mrow> </mml:math> </inline-formula>在所有的细胞在每一个时间点<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M244"> <mml:mn> 0.69921</米米l:mn> <mml:mo> ×</米米l:mo> <mml:msup> <mml:mrow> <mml:mn> 10</米米l:mn> </mml:mrow> <mml:mrow> <mml:mo> - - - - - -</米米l:mo> <mml:mn> 4</米米l:mn> </mml:mrow> </mml:msup> </mml:math> </inline-formula>这是接近的吗<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M245"> <mml:msup> <mml:mrow> <mml:mfenced separators="|"> <mml:mrow> <mml:mn> 1.3</米米l:mn> <mml:mo> ×</米米l:mo> <mml:msup> <mml:mrow> <mml:mn> 10</米米l:mn> </mml:mrow> <mml:mrow> <mml:mn> 4</米米l:mn> </mml:mrow> </mml:msup> </mml:mrow> </mml:mfenced> </mml:mrow> <mml:mrow> <mml:mo> - - - - - -</米米l:mo> <mml:mn> 1</米米l:mn> </mml:mrow> </mml:msup> </mml:math> </inline-formula>,然后一边<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M246"> <mml:msup> <mml:mrow> <mml:mfenced separators="|"> <mml:mrow> <mml:mn> 2.14</米米l:mn> <mml:mo> ×</米米l:mo> <mml:msup> <mml:mrow> <mml:mn> 10</米米l:mn> </mml:mrow> <mml:mrow> <mml:mn> 4</米米l:mn> </mml:mrow> </mml:msup> </mml:mrow> </mml:mfenced> </mml:mrow> <mml:mrow> <mml:mo> - - - - - -</米米l:mo> <mml:mn> 1</米米l:mn> </mml:mrow> </mml:msup> </mml:math> </inline-formula>最远的点的平均值吗<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M247"> <mml:mrow> <mml:msup> <mml:mrow> <mml:mfenced separators="|"> <mml:mrow> <mml:mo> ∂</米米l:mo> <mml:mi> E</米米l:mi> <mml:mo> /</米米l:mo> <mml:mo> ∂</米米l:mo> <mml:mi> T</米米l:mi> </mml:mrow> </mml:mfenced> </mml:mrow> <mml:mrow> <mml:mo> - - - - - -</米米l:mo> <mml:mn> 1</米米l:mn> </mml:mrow> </mml:msup> </mml:mrow> </mml:math> </inline-formula>。的原因,<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M248"> <mml:mo> ∂</米米l:mo> <mml:mi> E</米米l:mi> <mml:mo> /</米米l:mo> <mml:mo> ∂</米米l:mo> <mml:mi> T</米米l:mi> </mml:math> </inline-formula>可以在相对较高的水平靠近出口的启动或后期毛细管放电,可以发现它从您的定义或通过分析您的尺寸。</pgydF4y2Ba> <fig id="fig7"> <label>图7</l一个bel> <p>时空分布<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M249"> <mml:mo> ∂</米米l:mo> <mml:mi> E</米米l:mi> <mml:mo> /</米米l:mo> <mml:mo> ∂</米米l:mo> <mml:mi> T</米米l:mi> </mml:math> </inline-formula>。</pgydF4y2Ba> <graphic xlink:href="//www.newsama.com/downloads/journals/amp/2017/6473234.fig.007"></graphic> </fig> <p>正如上面提到的,准确的值<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M250"> <mml:mrow> <mml:mi> δ</米米l:mi> </mml:mrow> </mml:math> </inline-formula>当最低运行时实现也可能受到设置迭代初始值的方法。在图<gydF4y2Baxref ref-type="fig" rid="fig5"> 5</gydF4y2Baxref>最低运行时实现时,<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M251"> <mml:mrow> <mml:mi> δ</米米l:mi> </mml:mrow> </mml:math> </inline-formula>是<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M252"> <mml:mn> 0.54</米米l:mn> <mml:mo> ×</米米l:mo> <mml:msup> <mml:mrow> <mml:mn> 10</米米l:mn> </mml:mrow> <mml:mrow> <mml:mo> - - - - - -</米米l:mo> <mml:mn> 4</米米l:mn> </mml:mrow> </mml:msup> </mml:math> </inline-formula>而“Initial-T”<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M253"> <mml:mrow> <mml:mi> δ</米米l:mi> </mml:mrow> </mml:math> </inline-formula>是<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M254"> <mml:mn> 0.7</米米l:mn> <mml:mo> ×</米米l:mo> <mml:msup> <mml:mrow> <mml:mn> 10</米米l:mn> </mml:mrow> <mml:mrow> <mml:mo> - - - - - -</米米l:mo> <mml:mn> 4</米米l:mn> </mml:mrow> </mml:msup> </mml:math> </inline-formula>“Initial-S。“前一个不是像“Initial-S”的平均值<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M255"> <mml:mrow> <mml:msup> <mml:mrow> <mml:mfenced separators="|"> <mml:mrow> <mml:mo> ∂</米米l:mo> <mml:mi> E</米米l:mi> <mml:mo> /</米米l:mo> <mml:mo> ∂</米米l:mo> <mml:mi> T</米米l:mi> </mml:mrow> </mml:mfenced> </mml:mrow> <mml:mrow> <mml:mo> - - - - - -</米米l:mo> <mml:mn> 1</米米l:mn> </mml:mrow> </mml:msup> </mml:mrow> </mml:math> </inline-formula>这是<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M256"> <mml:mn> 0.69921</米米l:mn> <mml:mo> ×</米米l:mo> <mml:msup> <mml:mrow> <mml:mn> 10</米米l:mn> </mml:mrow> <mml:mrow> <mml:mo> - - - - - -</米米l:mo> <mml:mn> 4</米米l:mn> </mml:mrow> </mml:msup> </mml:math> </inline-formula>正如上面提到的。一个近似的解释如下:如图<gydF4y2Baxref ref-type="fig" rid="fig6b"> 6 (b)</gydF4y2Baxref>,循环计数在更高水平和在早期和晚期阶段模拟由于这个概要文件的输入电流。与此同时,在这两个时期,<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M257"> <mml:mo> ∂</米米l:mo> <mml:mi> E</米米l:mi> <mml:mo> /</米米l:mo> <mml:mo> ∂</米米l:mo> <mml:mi> T</米米l:mi> </mml:math> </inline-formula>也在相对较高的水平,如图<gydF4y2Baxref ref-type="fig" rid="fig7"> 7</gydF4y2Baxref>,所以<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M258"> <mml:mrow> <mml:mi> δ</米米l:mi> </mml:mrow> </mml:math> </inline-formula>可以调整到低一点的<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M259"> <mml:mrow> <mml:msup> <mml:mrow> <mml:mfenced separators="|"> <mml:mrow> <mml:mo> ∂</米米l:mo> <mml:mi> E</米米l:mi> <mml:mo> /</米米l:mo> <mml:mo> ∂</米米l:mo> <mml:mi> T</米米l:mi> </mml:mrow> </mml:mfenced> </mml:mrow> <mml:mrow> <mml:mo> - - - - - -</米米l:mo> <mml:mn> 1</米米l:mn> </mml:mrow> </mml:msup> </mml:mrow> </mml:math> </inline-formula>,这是<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M260"> <mml:msup> <mml:mrow> <mml:mfenced separators="|"> <mml:mrow> <mml:mn> 2.17</米米l:mn> <mml:mo> ×</米米l:mo> <mml:msup> <mml:mrow> <mml:mn> 10</米米l:mn> </mml:mrow> <mml:mrow> <mml:mn> 4</米米l:mn> </mml:mrow> </mml:msup> </mml:mrow> </mml:mfenced> </mml:mrow> <mml:mrow> <mml:mo> - - - - - -</米米l:mo> <mml:mn> 1</米米l:mn> </mml:mrow> </mml:msup> </mml:math> </inline-formula>,减少循环计数早期和晚期阶段的模拟。这种调整也是原因<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M261"> <mml:mrow> <mml:mi> δ</米米l:mi> </mml:mrow> </mml:math> </inline-formula>成为进一步的一侧<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M262"> <mml:msup> <mml:mrow> <mml:mfenced separators="|"> <mml:mrow> <mml:mn> 1.3</米米l:mn> <mml:mo> ×</米米l:mo> <mml:msup> <mml:mrow> <mml:mn> 10</米米l:mn> </mml:mrow> <mml:mrow> <mml:mn> 4</米米l:mn> </mml:mrow> </mml:msup> </mml:mrow> </mml:mfenced> </mml:mrow> <mml:mrow> <mml:mo> - - - - - -</米米l:mo> <mml:mn> 1</米米l:mn> </mml:mrow> </mml:msup> </mml:math> </inline-formula>在中间的模拟。然而,这一时期的循环计数可能仍然很低,因为目前也达到高峰和等离子体流属性不变的大约如前所述。模拟的“Initial-S”,总不能减少了通过调整运行时<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M263"> <mml:mrow> <mml:mi> δ</米米l:mi> </mml:mrow> </mml:math> </inline-formula>任何一边。循环计数的地区相对较高的出口附近的毛细管,如图<gydF4y2Baxref ref-type="fig" rid="fig6a"> 6(一)</gydF4y2Baxref>。在这个地区,<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M264"> <mml:mo> ∂</米米l:mo> <mml:mi> E</米米l:mi> <mml:mo> /</米米l:mo> <mml:mo> ∂</米米l:mo> <mml:mi> T</米米l:mi> </mml:math> </inline-formula>增加和减少沿时间方向和不高或低,如图<gydF4y2Baxref ref-type="fig" rid="fig7"> 7</gydF4y2Baxref>。有可能是平均的<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M265"> <mml:mo> ∂</米米l:mo> <mml:mi> E</米米l:mi> <mml:mo> /</米米l:mo> <mml:mo> ∂</米米l:mo> <mml:mi> T</米米l:mi> </mml:math> </inline-formula>出口附近的毛细管甚至可能发生<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M266"> <mml:mn> 0.7</米米l:mn> <mml:mo> ×</米米l:mo> <mml:msup> <mml:mrow> <mml:mn> 10</米米l:mn> </mml:mrow> <mml:mrow> <mml:mo> - - - - - -</米米l:mo> <mml:mn> 4</米米l:mn> </mml:mrow> </mml:msup> </mml:math> </inline-formula>,所以没有这样的一面<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M267"> <mml:mrow> <mml:mi> δ</米米l:mi> </mml:mrow> </mml:math> </inline-formula>可以调整为更少的总运行时间与“Initial-S模拟。“为了进一步减少运行时,一个可能的方法可能是整个时空区域划分成3 - 4小区域和设置<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M268"> <mml:mrow> <mml:mi> δ</米米l:mi> </mml:mrow> </mml:math> </inline-formula>这些地区分别,所以<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M269"> <mml:mrow> <mml:mi> δ</米米l:mi> </mml:mrow> </mml:math> </inline-formula>可以接近平均<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M270"> <mml:mrow> <mml:msup> <mml:mrow> <mml:mfenced separators="|"> <mml:mrow> <mml:mo> ∂</米米l:mo> <mml:mi> E</米米l:mi> <mml:mo> /</米米l:mo> <mml:mo> ∂</米米l:mo> <mml:mi> T</米米l:mi> </mml:mrow> </mml:mfenced> </mml:mrow> <mml:mrow> <mml:mo> - - - - - -</米米l:mo> <mml:mn> 1</米米l:mn> </mml:mrow> </mml:msup> </mml:mrow> </mml:math> </inline-formula>相应的时空区域。这可能有助于避免统一的情况<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M271"> <mml:mrow> <mml:mi> δ</米米l:mi> </mml:mrow> </mml:math> </inline-formula>远离是真的吗<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M272"> <mml:mrow> <mml:msup> <mml:mrow> <mml:mfenced separators="|"> <mml:mrow> <mml:mo> ∂</米米l:mo> <mml:mi> E</米米l:mi> <mml:mo> /</米米l:mo> <mml:mo> ∂</米米l:mo> <mml:mi> T</米米l:mi> </mml:mrow> </mml:mfenced> </mml:mrow> <mml:mrow> <mml:mo> - - - - - -</米米l:mo> <mml:mn> 1</米米l:mn> </mml:mrow> </mml:msup> </mml:mrow> </mml:math> </inline-formula>在出口处等这些地区的毛细管的开始模拟。</pgydF4y2Ba> </sec> <sec id="sec7"> <title>7所示。结论</t我tle> <p>方法基于的工作Trayner Głowacki和扎格卢勒放电介绍解决萨哈方程的模拟中使用的毛细管放电等或等发射系统。通过修改特定的能量方程,解决萨哈方程的多维问题降低到一维迭代方程的温度。一个系数<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M273"> <mml:mrow> <mml:mi> δ</米米l:mi> </mml:mrow> </mml:math> </inline-formula>引入迭代方程满足迭代法的收敛条件。为了提高计算效率,包括优化两种方法<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M274"> <mml:mrow> <mml:mi> δ</米米l:mi> </mml:mrow> </mml:math> </inline-formula>介绍了动态和设置迭代初始值。几个数值模拟中使用的毛细管放电等或等发射系统进行测试的详细性能提出了方法。程序的循环计数(i) - (vi)和模拟运行时间记录研究方法的效率。结果表明,仿真与优化运行时可以显著降低<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M275"> <mml:mrow> <mml:mi> δ</米米l:mi> </mml:mrow> </mml:math> </inline-formula>和“Initial-T。Initial-T“首先,使迭代初始值可能更接近真正的根细胞;其次,“Initial-T”可以利用脉冲的输入电流使<我nl在e- - - - - -for米ul一个> <mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" id="M276"> <mml:mrow> <mml:mi> δ</米米l:mi> </mml:mrow> </mml:math> </inline-formula>可以进一步优化的最小运行时间与“Initial-S。”</pgydF4y2Ba> </sec> <back> <sec> <title>的利益冲突</t我tle> <p>作者宣称没有利益冲突有关的出版。</pgydF4y2Ba> </sec> <ref-list> <ref id="B1" content-type="article"> <label>1</l一个bel> <element-citation publication-type="journal"> <person-group person-group-type="author"> <name> <surname> Salge</年代urn一个米e> <given-names> j·g . H。</g我ven- - - - - -n一个米e年代> </name> <name> <surname> 怀斯</年代urn一个米e> <given-names> t·h·G·G。</g我ven- - - - - -n一个米e年代> </name> <name> <surname> Braunsberger</年代urn一个米e> <given-names> U。</g我ven- - - - - -n一个米e年代> </name> <name> <surname> Fien</年代urn一个米e> <given-names> H。</g我ven- - - - - -n一个米e年代> </name> <name> <surname> Loffler</年代urn一个米e> <given-names> M。</g我ven- - - - - -n一个米e年代> </name> <name> <surname> 威特</年代urn一个米e> <given-names> W。</g我ven- - - - - -n一个米e年代> </name> <name> <surname> Zwingel</年代urn一个米e> <given-names> D。</g我ven- - - - - -n一个米e年代> </name> <name> <surname> Zocha</年代urn一个米e> <given-names> K。</g我ven- - - - - -n一个米e年代> </name> </person-group> <article-title> 质量乘以加速度脉冲等离子体</一个rt我cle- - - - - -t我tle> <source> <italic> IEEE磁学</我t一个l我c> <year> 1989年</gydF4y2Bayear> <volume> 25</volu米e> <issue> 1</我年代年代ue> <fpage> 495年</fp一个ge> <lpage> 499年</lp一个ge> <pub-id pub-id-type="doi"> 10.1109/20.22588</pugydF4y2Bab-id> </element-citation> </ref> <ref id="B2" content-type="article"> <label>2</l一个bel> <element-citation publication-type="journal"> <person-group person-group-type="author"> <name> <surname> 勒布</年代urn一个米e> <given-names> 一个。</g我ven- - - - - -n一个米e年代> </name> <name> <surname> 卡普兰</年代urn一个米e> <given-names> Z。</g我ven- - - - - -n一个米e年代> </name> </person-group> <article-title> 物理过程的理论模型在高压放电的电热发射器</一个rt我cle- - - - - -t我tle> <source> <italic> IEEE磁学</我t一个l我c> <year> 1989年</gydF4y2Bayear> <volume> 25</volu米e> <issue> 1</我年代年代ue> <fpage> 342年</fp一个ge> <lpage> 346年</lp一个ge> <pub-id pub-id-type="other"> 2 - s2.0 - 0024479226</pugydF4y2Bab-id> <pub-id pub-id-type="doi"> 10.1109/20.22561</pugydF4y2Bab-id> </element-citation> </ref> <ref id="B3" content-type="article"> <label>3</l一个bel> <element-citation publication-type="journal"> <person-group person-group-type="author"> <name> <surname> Kovitya</年代urn一个米e> <given-names> P。</g我ven- - - - - -n一个米e年代> </name> <name> <surname> Lowke</年代urn一个米e> <given-names> J·J。</g我ven- - - - - -n一个米e年代> </name> </person-group> <article-title> 理论预测ablation-stabilised弧在圆柱管</一个rt我cle- - - - - -t我tle> <source> <italic> 物理学学报D辑:应用物理</我t一个l我c> <year> 1984年</gydF4y2Bayear> <volume> 17</volu米e> <issue> 6</我年代年代ue> <fpage> 1197年</fp一个ge> <lpage> 1212年</lp一个ge> <pub-id pub-id-type="doi"> 10.1088 / 0022 - 3727/17/6/016</pugydF4y2Bab-id> <pub-id pub-id-type="other"> 2 - s2.0 - 0021439802</pugydF4y2Bab-id> </element-citation> </ref> <ref id="B4" content-type="techreport"> <label>4</l一个bel> <element-citation publication-type="gov"> <person-group person-group-type="author"> <name> <surname> 鲍威尔</年代urn一个米e> <given-names> j . D。</g我ven- - - - - -n一个米e年代> </name> <name> <surname> Zielinski</年代urn一个米e> <given-names> 答:E。</g我ven- - - - - -n一个米e年代> </name> </person-group> <article-title> 理论和实验的消融毛细管放电和应用程序electrothermal-chemical枪支</一个rt我cle- - - - - -t我tle> <source> <italic> 弹道研究实验室的报告</我t一个l我c> <year> 1992年</gydF4y2Bayear> <issue> brl - tr - 3355</我年代年代ue> </element-citation> </ref> <ref id="B5" content-type="article"> <label>5</l一个bel> <element-citation publication-type="journal"> <person-group person-group-type="author"> <name> <surname> Zoler</年代urn一个米e> <given-names> D。</g我ven- - - - - -n一个米e年代> </name> <name> <surname> Cuperman</年代urn一个米e> <given-names> 年代。</g我ven- - - - - -n一个米e年代> </name> <name> <surname> Ashkenazy</年代urn一个米e> <given-names> J。</g我ven- - - - - -n一个米e年代> </name> <name> <surname> 癌症</年代urn一个米e> <given-names> M。</g我ven- - - - - -n一个米e年代> </name> <name> <surname> 卡普兰</年代urn一个米e> <given-names> Z。</g我ven- - - - - -n一个米e年代> </name> </person-group> <article-title> 高压放电的时间模型在狭窄的烧蚀毛细血管</一个rt我cle- - - - - -t我tle> <source> <italic> 等离子体物理杂志</我t一个l我c> <year> 1993年</gydF4y2Bayear> <volume> 50</volu米e> <issue> 1</我年代年代ue> <fpage> 51</fp一个ge> <lpage> 70年</lp一个ge> <pub-id pub-id-type="other"> 2 - s2.0 - 84972273810</pugydF4y2Bab-id> <pub-id pub-id-type="doi"> 10.1017 / S0022377800026908</pugydF4y2Bab-id> </element-citation> </ref> <ref id="B6" content-type="phdthesis"> <label>6</l一个bel> <element-citation publication-type="thesis"> <person-group person-group-type="author"> <name> <surname> 赫尔利</年代urn一个米e> <given-names> j . D。</g我ven- - - - - -n一个米e年代> </name> </person-group> <source> <italic> In-bore电热等离子体的诊断和建模发射器(博士学位。论文)</我t一个l我c> <year> 1993年</gydF4y2Bayear> <publisher-loc> 美国罗利数控</pugydF4y2Bablisher-loc> <publisher-name> 北卡州立大学</pugydF4y2Bablisher-name> </element-citation> </ref> <ref id="B7" content-type="phdthesis"> <label>7</l一个bel> <element-citation publication-type="thesis"> <person-group person-group-type="author"> <name> <surname> 拉贾</年代urn一个米e> <given-names> L . L。</g我ven- - - - - -n一个米e年代> </name> </person-group> <source> <italic> 金属蒸汽等离子体放电的理论模型的电热炮[博士。论文)</我t一个l我c> <year> 1996年</gydF4y2Bayear> <publisher-loc> 美国德克萨斯州奥斯汀市</pugydF4y2Bablisher-loc> <publisher-name> 德克萨斯大学的</pugydF4y2Bablisher-name> </element-citation> </ref> <ref id="B8" content-type="phdthesis"> <label>8</l一个bel> <element-citation publication-type="thesis"> <person-group person-group-type="author"> <name> <surname> 温弗瑞</年代urn一个米e> <given-names> a . L。</g我ven- - - - - -n一个米e年代> </name> </person-group> <source> <italic> 数值研究的理想行为、参数和电热的新颖应用等离子体源(博士学位。论文)</我t一个l我c> <year> 2010年</gydF4y2Bayear> <publisher-loc> 美国罗利数控</pugydF4y2Bablisher-loc> <publisher-name> 北卡州立大学</pugydF4y2Bablisher-name> </element-citation> </ref> <ref id="B9" content-type="article"> <label>9</l一个bel> <element-citation publication-type="journal"> <person-group person-group-type="author"> <name> <surname> 金</年代urn一个米e> <given-names> K。</g我ven- - - - - -n一个米e年代> </name> </person-group> <article-title> 时间一维模型的脉冲等离子体放电毛细管等离子体装置</一个rt我cle- - - - - -t我tle> <source> <italic> IEEE等离子体科学</我t一个l我c> <year> 2003年</gydF4y2Bayear> <volume> 31日</volu米e> <issue> 4</我年代年代ue> <fpage> 729年</fp一个ge> <lpage> 735年</lp一个ge> <pub-id pub-id-type="doi"> 10.1109 / TPS.2003.815472</pugydF4y2Bab-id> <pub-id pub-id-type="other"> 2 - s2.0 - 0041927741</pugydF4y2Bab-id> </element-citation> </ref> <ref id="B10" content-type="article"> <label>10</l一个bel> <element-citation publication-type="journal"> <person-group person-group-type="author"> <name> <surname> 温弗瑞</年代urn一个米e> <given-names> a . L。</g我ven- - - - - -n一个米e年代> </name> <name> <surname> 米塔尔</年代urn一个米e> <given-names> 年代。</g我ven- - - - - -n一个米e年代> </name> <name> <surname> Abd Al-Halim</年代urn一个米e> <given-names> m·A。</g我ven- - - - - -n一个米e年代> </name> </person-group> <article-title> 超级喷气式focused-injection锥形毛细管等离子体源的概念</一个rt我cle- - - - - -t我tle> <source> <italic> IEEE等离子体科学</我t一个l我c> <year> 2015年</gydF4y2Bayear> <volume> 43</volu米e> <issue> 6</我年代年代ue> <fpage> 2104年</fp一个ge> <lpage> 2109年</lp一个ge> <pub-id pub-id-type="other"> 2 - s2.0 - 84932650378</pugydF4y2Bab-id> <pub-id pub-id-type="doi"> 10.1109 / TPS.2015.2417878</pugydF4y2Bab-id> </element-citation> </ref> <ref id="B11" content-type="article"> <label>11</l一个bel> <element-citation publication-type="journal"> <person-group person-group-type="author"> <name> <surname> 金</年代urn一个米e> <given-names> K。</g我ven- - - - - -n一个米e年代> </name> </person-group> <article-title> 毛细管等离子体流的数值模拟所产生的大电流脉冲的力量</一个rt我cle- - - - - -t我tle> <source> <italic> 国际热科学杂志》上</我t一个l我c> <year> 2005年</gydF4y2Bayear> <volume> 44</volu米e> <issue> 11</我年代年代ue> <fpage> 1039年</fp一个ge> <lpage> 1046年</lp一个ge> <pub-id pub-id-type="other"> 2 - s2.0 - 25444525847</pugydF4y2Bab-id> <pub-id pub-id-type="doi"> 10.1016 / j.ijthermalsci.2005.04.007</pugydF4y2Bab-id> </element-citation> </ref> <ref id="B12" content-type="article"> <label>12</l一个bel> <element-citation publication-type="journal"> <person-group person-group-type="author"> <name> <surname> 扎格卢勒</年代urn一个米e> <given-names> m·R。</g我ven- - - - - -n一个米e年代> </name> </person-group> <article-title> 制定和有效的算法减少测定平衡组成和配分函数的理想和非理想的复杂的等离子体混合物</一个rt我cle- - - - - -t我tle> <source> <italic> 物理评论E:统计、非线性和软物质物理学</我t一个l我c> <year> 2004年</gydF4y2Bayear> <volume> 69年</volu米e> <issue> 2</我年代年代ue> <fpage> 1</fp一个ge> <lpage> 7</lp一个ge> <pub-id pub-id-type="publisher-id"> 026702年</pugydF4y2Bab-id> <pub-id pub-id-type="doi"> 10.1103 / PhysRevE.69.026702</pugydF4y2Bab-id> <pub-id pub-id-type="other"> 2 - s2.0 - 42749100123</pugydF4y2Bab-id> </element-citation> </ref> <ref id="B13" content-type="article"> <label>13</l一个bel> <element-citation publication-type="journal"> <person-group person-group-type="author"> <name> <surname> Trayner</年代urn一个米e> <given-names> C。</g我ven- - - - - -n一个米e年代> </name> <name> <surname> Głowacki</年代urn一个米e> <given-names> m . H。</g我ven- - - - - -n一个米e年代> </name> </person-group> <article-title> 萨哈方程的新的技术解决方案</一个rt我cle- - - - - -t我tle> <source> <italic> 杂志的科学计算</我t一个l我c> <year> 1995年</gydF4y2Bayear> <volume> 10</volu米e> <issue> 1</我年代年代ue> <fpage> 139年</fp一个ge> <lpage> 149年</lp一个ge> <pub-id pub-id-type="doi"> 10.1007 / BF02087963</pugydF4y2Bab-id> <pub-id pub-id-type="other"> 2 - s2.0 - 0029277714</pugydF4y2Bab-id> </element-citation> </ref> <ref id="B14" content-type="article"> <label>14</l一个bel> <element-citation publication-type="journal"> <person-group person-group-type="author"> <name> <surname> Zoler</年代urn一个米e> <given-names> D。</g我ven- - - - - -n一个米e年代> </name> <name> <surname> Saphier</年代urn一个米e> <given-names> D。</g我ven- - - - - -n一个米e年代> </name> <name> <surname> Alimi</年代urn一个米e> <given-names> R。</g我ven- - - - - -n一个米e年代> </name> </person-group> <article-title> 数值研究进化的烧蚀毛细管放电等离子体参数的双脉冲形式的输入能量</一个rt我cle- - - - - -t我tle> <source> <italic> 物理学学报D辑:应用物理</我t一个l我c> <year> 1994年</gydF4y2Bayear> <volume> 27</volu米e> <issue> 7</我年代年代ue> <fpage> 1423年</fp一个ge> <lpage> 1432年</lp一个ge> <pub-id pub-id-type="other"> 2 - s2.0 - 0028464164</pugydF4y2Bab-id> <pub-id pub-id-type="doi"> 10.1088 / 0022 - 3727/27/7/013</pugydF4y2Bab-id> </element-citation> </ref> <ref id="B15" content-type="article"> <label>15</l一个bel> <element-citation publication-type="journal"> <person-group person-group-type="author"> <name> <surname> Batteh</年代urn一个米e> <given-names> J。</g我ven- - - - - -n一个米e年代> </name> <name> <surname> 鲍威尔</年代urn一个米e> <given-names> J。</g我ven- - - - - -n一个米e年代> </name> <name> <surname> 水槽</年代urn一个米e> <given-names> D。</g我ven- - - - - -n一个米e年代> </name> <name> <surname> Thornhill</年代urn一个米e> <given-names> l</g我ven- - - - - -n一个米e年代> </name> </person-group> <article-title> 的方法计算等离子体混合物的热力学和输运性质等喷油器</一个rt我cle- - - - - -t我tle> <source> <italic> IEEE磁学</我t一个l我c> <year> 1995年</gydF4y2Bayear> <volume> 31日</volu米e> <issue> 1</我年代年代ue> <fpage> 388年</fp一个ge> <lpage> 393年</lp一个ge> <pub-id pub-id-type="doi"> 10.1109/20.364656</pugydF4y2Bab-id> <pub-id pub-id-type="other"> 2 - s2.0 - 0029230306</pugydF4y2Bab-id> </element-citation> </ref> <ref id="B16" content-type="article"> <label>16</l一个bel> <element-citation publication-type="journal"> <person-group person-group-type="author"> <name> <surname> Cuperman</年代urn一个米e> <given-names> 年代。</g我ven- - - - - -n一个米e年代> </name> <name> <surname> Zoler</年代urn一个米e> <given-names> D。</g我ven- - - - - -n一个米e年代> </name> <name> <surname> Ashkenazy</年代urn一个米e> <given-names> J。</g我ven- - - - - -n一个米e年代> </name> </person-group> <article-title> 分析结合放电毛细管烧蚀等离子体临界流的管道系统</一个rt我cle- - - - - -t我tle> <source> <italic> 等离子体科学和技术来源</我t一个l我c> <year> 1994年</gydF4y2Bayear> <volume> 3</volu米e> <issue> 4</我年代年代ue> <fpage> 593年</fp一个ge> <lpage> 601年</lp一个ge> <pub-id pub-id-type="doi"> 10.1088 / 0963 - 0252/3/4/017</pugydF4y2Bab-id> <pub-id pub-id-type="other"> 2 - s2.0 - 5844261494</pugydF4y2Bab-id> </element-citation> </ref> <ref id="B17" content-type="article"> <label>17</l一个bel> <element-citation publication-type="journal"> <person-group person-group-type="author"> <name> <surname> 维斯</年代urn一个米e> <given-names> E。</g我ven- - - - - -n一个米e年代> </name> <name> <surname> Zoler</年代urn一个米e> <given-names> D。</g我ven- - - - - -n一个米e年代> </name> <name> <surname> 瓦尔德</年代urn一个米e> <given-names> 年代。</g我ven- - - - - -n一个米e年代> </name> <name> <surname> 伊莱亚斯</年代urn一个米e> <given-names> E。</g我ven- - - - - -n一个米e年代> </name> </person-group> <article-title> 建模的电气confined-capillary-discharge放电区延长额外的管道</一个rt我cle- - - - - -t我tle> <source> <italic> 物理信</我t一个l我c> <year> 2009年</gydF4y2Bayear> <volume> 373年</volu米e> <issue> 10</我年代年代ue> <fpage> 972年</fp一个ge> <lpage> 975年</lp一个ge> <pub-id pub-id-type="other"> 2 - s2.0 - 59849085305</pugydF4y2Bab-id> <pub-id pub-id-type="doi"> 10.1016 / j.physleta.2009.01.032</pugydF4y2Bab-id> </element-citation> </ref> <ref id="B18" content-type="article"> <label>18</l一个bel> <element-citation publication-type="journal"> <person-group person-group-type="author"> <name> <surname> 扎格卢勒</年代urn一个米e> <given-names> m·R。</g我ven- - - - - -n一个米e年代> </name> <name> <surname> Bourham</年代urn一个米e> <given-names> m·A。</g我ven- - - - - -n一个米e年代> </name> <name> <surname> 难道</年代urn一个米e> <given-names> j . M。</g我ven- - - - - -n一个米e年代> </name> </person-group> <article-title> 半解析建模和模拟进化和Ohmically-heated理想的等离子体流的电热枪支</一个rt我cle- - - - - -t我tle> <source> <italic> 物理学学报D辑:应用物理</我t一个l我c> <year> 2001年</gydF4y2Bayear> <volume> 34</volu米e> <issue> 5</我年代年代ue> <fpage> 772年</fp一个ge> <lpage> 786年</lp一个ge> <pub-id pub-id-type="other"> 2 - s2.0 - 0035819844</pugydF4y2Bab-id> <pub-id pub-id-type="doi"> 10.1088 / 0022 - 3727/34/5/316</pugydF4y2Bab-id> </element-citation> </ref> <ref id="B19" content-type="misc"> <label>19</l一个bel> <element-citation publication-type="other"> <person-group person-group-type="author"> <name> <surname> Kramida</年代urn一个米e> <given-names> 一个。</g我ven- - - - - -n一个米e年代> </name> <name> <surname> Ralchenko</年代urn一个米e> <given-names> Y。</g我ven- - - - - -n一个米e年代> </name> <name> <surname> 读者</年代urn一个米e> <given-names> J。</g我ven- - - - - -n一个米e年代> </name> <etal></etal> </person-group> <article-title> 原子光谱数据库</一个rt我cle- - - - - -t我tle> <comment> 2015年版本5日,<egydF4y2Baxt-link ext-link-type="url" xlink:href="https://www.nist.gov/pml/atomic-spectra-database"> https://www.nist.gov/pml/atomic-spectra-database</egydF4y2Baxt-link> </comment> </element-citation> </ref> <ref id="B20" content-type="incollection"> <label>20.</l一个bel> <element-citation publication-type="book"> <person-group person-group-type="author"> <name> <surname> Chapra</年代urn一个米e> <given-names> s . C。</g我ven- - - - - -n一个米e年代> </name> <name> <surname> Canale</年代urn一个米e> <given-names> r P。</g我ven- - - - - -n一个米e年代> </name> </person-group> <article-title> 开放的方法</一个rt我cle- - - - - -t我tle> <source> <italic> 数值方法为工程师</我t一个l我c> <year> 2010年</gydF4y2Bayear> <edition> 6日</ed我t我on> <publisher-loc> 纽约,纽约,美国</pugydF4y2Bablisher-loc> <publisher-name> 麦格劳-希尔</pugydF4y2Bablisher-name> <fpage> 142年</fp一个ge> <lpage> 173年</lp一个ge> </element-citation> </ref> <ref id="B21" content-type="article"> <label>21</l一个bel> <element-citation publication-type="journal"> <person-group person-group-type="author"> <name> <surname> 鲍威尔</年代urn一个米e> <given-names> j . D。</g我ven- - - - - -n一个米e年代> </name> <name> <surname> Zielinski</年代urn一个米e> <given-names> 答:E。</g我ven- - - - - -n一个米e年代> </name> </person-group> <article-title> 电热的毛细管放电枪</一个rt我cle- - - - - -t我tle> <source> <italic> IEEE磁学</我t一个l我c> <year> 1993年</gydF4y2Bayear> <volume> 29日</volu米e> <issue> 1</我年代年代ue> <fpage> 591年</fp一个ge> <lpage> 596年</lp一个ge> <pub-id pub-id-type="other"> 2 - s2.0 - 0001012234</pugydF4y2Bab-id> <pub-id pub-id-type="doi"> 10.1109/20.195642</pugydF4y2Bab-id> </element-citation> </ref> <ref id="B22" content-type="article"> <label>22</l一个bel> <element-citation publication-type="journal"> <person-group person-group-type="author"> <name> <surname> Zollweg</年代urn一个米e> <given-names> r . J。</g我ven- - - - - -n一个米e年代> </name> <name> <surname> 利伯曼</年代urn一个米e> <given-names> r·W。</g我ven- - - - - -n一个米e年代> </name> </person-group> <article-title> 不理想的等离子体的导电性</一个rt我cle- - - - - -t我tle> <source> <italic> 应用物理杂志</我t一个l我c> <year> 1987年</gydF4y2Bayear> <volume> 62年</volu米e> <issue> 9</我年代年代ue> <fpage> 3621年</fp一个ge> <lpage> 3627年</lp一个ge> <pub-id pub-id-type="other"> 2 - s2.0 - 0003180629</pugydF4y2Bab-id> <pub-id pub-id-type="doi"> 10.1063/1.339265</pugydF4y2Bab-id> </element-citation> </ref> </ref-list> </back> </article> </body> </html>