×

等离子体鞘层的数值稳定性。 (英语。法语摘要) Zbl 1452.82022号

摘要:我们有兴趣开发一种数值方法来捕获等离子体与金属壁接触形成的固定鞘层。本工作基于[第一作者等,Kinet.Relat.Models 9,No.4,621-656(2016;Zbl 1360.35166号)]. 本工作中的主要问题是了解经典含时Vlasov-Ampère数值格式在长时间内如何准确地保持这些具有发射/吸收边界条件的非均匀定态解。在高阶半拉格朗日方法的背景下,由于其模板较大,域边界附近的插值也需要特殊处理。

MSC公司:

82D10号 等离子体统计力学
65平方米 含偏微分方程初值和初边值问题数值方法的稳定性和收敛性
PDF格式BibTeX公司 XML格式引用

参考文献:

[1] B.Afeyan、F.Casas、N.Crouseilles、A.Dodhy、E.Faou、M.Mehrenberger和E.Sonnendr¨ucker,用变速分辨率网格和高阶时间分裂模拟动态静电电子非线性波,《欧洲物理杂志》D,68(2014),第295页。
[2] M.Badsi,双动力学鞘层模型的线性电子稳定性,《数学分析与应用杂志》,453(2017),第954-972页·兹比尔1373.82085
[3] M.Badsi、M.Campos Pinto和B.Despr´es,双动力学鞘的最小化公式,Kinet。相关。模型,9(2016)·Zbl 1360.35166号
[4] N.Besse和M.Mehrenberger,vlasov-poisson系统高阶半拉格朗日格式类的收敛性,计算数学,77(2008),第93-123页·Zbl 1131.65080号
[5] D.Coulette和G.Manfredi,《等离子体壁相互作用的欧拉-弗拉索夫密码》,《物理学杂志:会议系列》,561(2014)。
[6] N.Crouseilles、M.Mehrenberger和F.Vecil,Vlasov-Poisson的间断Galerkin半拉格朗日方法,ESAIM:Proc。,32(2011年),第211-230页·Zbl 1302.76099号
[7] E.Sonnendr¨ucker,弗拉索夫方程的数值方法,准备中。图标题27图标题图1。时间t=0(顶部)和t=4(底部)时的电子分布函数。(29)中给出的参数。图2。t=0(顶部)和t=4(底部)时的离子分布函数。(29)中给出的参数。图3。时间t=4(左)和t=8(右)时,电子分布函数fse(0,x,v)−fse(t,x,v)的误差。(29)中给出的参数。图4。时间t=4(左)和t=8(右)时,离子分布函数fsi(0,x,v)−fsi(t,x,v)的误差。(29)中给出的参数。图5。时间t=8时电子(左)和离子(右)的分布函数。(29)中给出的参数。图6。输入x=0时电流密度J(t,x=0)的时间演变,并放大时间间隔[0,0.1],[4,4.01]和[8,8.01]。(29)中给出的参数。图7。总能量E(t)(左)和总电荷ρtot(t)的时间演化(右)。(29)中给出的参数。图8。电子(左)和离子(右)的总密度和分布函数的L1范数的时间演化。(29)中给出的参数。图9。电子(左)和离子(右)分布函数L2范数的时间演化。(29)中给出的参数。图10。总能量E(t)的时间演化。左图:d=1、d=4和d=5边界处蝶形程序和无中心插值的比较。右:无中心插值的两种实现方式的比较(无中心:使用边界2d+1最近点上的插值多项式计算外推值,然后使用与蝶形过程相同的算法,无中心实现2:直接实现插值)。图11。入口x=0(左)和总能量(右)处电流密度的比较,将∆t=10−5变为∆t=10−4。图12。入口x=0(左)和总能量(右)处电流密度的比较,将∆t=10−5变为∆t=10−4,Nxto 256和Nvto 512。图13.插值次数d=8和d=0时,入口x=0(左)处的电流密度和总能量(右)的比较。28图标题图14。对于d=0、∆t=10−4和不同的空间/速度参数,入口x=0处的电流密度(左:短时间;右:长时间)图15。d=0、∆t=10−4和不同空间/速度参数的总能量。图29图1。时间t=0(顶部)和t=4(底部)时的电子分布函数。(29)中给出的参数。30图2。t=0(顶部)和t=4(底部)时的离子分布函数。(29)中给出的参数。图31图3。时间t=4(左)和t=8(右)时,电子分布函数fse(0,x,v)−fse(t,x,v)的误差。(29)中给出的参数。图4。时间t=4(左)和t=8(右)时,离子分布函数fsi(0,x,v)−fsi(t,x,v)的误差。(29)中给出的参数。图5。时间t=8时电子(左)和离子(右)的分布函数。(29)中给出的参数。32图6。输入x=0时电流密度J(t,x=0)的时间演变,并放大时间间隔[0,0.1],[4,4.01]和[8,8.01]。(29)中给出的参数。图7。总能量E(t)(左)和总电荷ρtot(t)的时间演化(右)。(29)中给出的参数。图33 1电子分布函数的范数图8。总密度和L(左)和离子(右)的时间演化。(29)中给出的参数。2图9。电子(左)和离子(右)分布函数Lnorm的时间演化。(29)中给出的参数。34图10。总能量E(t)的时间演化。左图:d=1、d=4和d=5边界处蝶形程序和无中心插值的比较。右:无中心插值的两种实现方式的比较(无中心:使用边界2d+1最近点上的插值多项式计算外推值,然后使用与蝶形过程相同的算法,无中心实现2:直接实现插值)。图11。入口x=0(左)和总能量(右)处电流密度的比较,将∆t=10−5变为∆t=10−4。图35图12。入口x=0(左)和总能量(右)处电流密度的比较,将∆t=10−5变为∆t=10−4,Nxto 256和Nvto 512。图13。对于插值次数d=8和d=0,比较入口x=0(左)处的电流密度和总能量(右)。36图−4和不同的图14。对于d=0,∆t=10空间/速度参数,入口处的电流密度x=0(左:短时间;右:长时间)图15。d=0、∆t=10−4和不同空间/速度参数的总能量。
此参考列表基于出版商或数字数学图书馆提供的信息。其项与zbMATH标识符进行启发式匹配,可能包含数据转换错误。在某些情况下,zbMATH Open的数据对这些数据进行了补充/增强。这试图尽可能准确地反映原始论文中列出的参考文献,而不要求完整或完全匹配。