一类半导体模型的几种数值方法比较

来源 :山东大学 | 被引量 : 0次 | 上传用户:liusheng123321
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
间断有限元方法集合了高分辨率有限差分方法和有限体积方法的优点,它是采用完全分片的多项式空间对近似解和试验函数进行空间离散,使用Runge-kutta方法进行时间分解的有限元方法。自八十年代末以来引起了一些数学家的注意,也因此得到了很好的发展.本文主要应用这种方法来分析一类半导体模型,构造相应的LDG格式,给出了稳定性分析,最后并给出了数值实验。  本文共分为五章。  第一章为引言,简单介绍了间断有限元方法的由来和发展,并对这种数值方法的有点做出总结.  第二章用局部间断有限元方法讨论如下半导体模型:  nt+(μnE)x=((τ)θnx)x,[x,t]∈[a,b]×[0,T],(0.0.1)  Φxx=e/ε(n-nd).(0.0.2)  -Φx=E.(0.0.3)。  边界条件由电势Φ给出,Dirichlet边界条件:  Φ(0,t)=Φ0(t),Φ(0.6,t)=Φ(0,t)+vbias,(0.0.4).  vbias表示电压差。n满足周期性边界条件,且满足初始条件:  n(x,0)=n0(x),(0.0.5)。  其中,n表示空穴浓度,nd表示掺杂浓度,即杂质造成的间断点,μ表示电子的迁移率,(τ)=mμ/e表示松弛参数,m表示电子有效质量,e表示电子电量,ε介电常数,θ=kh/mT,kb表示Boltzmann常数,T为初始温度,Φ表示电势,E为电场强度。  我们对方程(0.0.1)构造出半离散格式:  引入辅助变量q=√(τ)0nx,找到nh,qh,Φh和Eh∈Vhk,使得对(∨)v(x),w(x),z(x),r(x)∈Vkh,满足:  ∫Ijnhtvdx-∫Ij(μnhEh-√(τ)θqh)vxdx+(μ(n)(E)-(q))j+1/2v-j+1/2-(μ(n)(E)-(q))j-1/2v+j-1/2=0,(0.0.6)  ∫Ij√(τ)θqhwdx+√(τ)θnhwxdx-(√(τ)θ(n))j+1/2w-j+1/2+(√(τ)θ(n))j-1/2w+j-1/2=0,(0.0.7)  ∫IjEhzxdx-(E)j+1/1z-j+1/2+(E)j-1/2z+j-1/2=∫Ije/ε(nh-nd)zdx,(0.0.8)  ∫Ij(-Φhx)rxdx-((Φ)x)j+1/1r-j+1/2+((Φ)x)j-1/2r+j-1/2=∫IjEhrd,(0.0.9).  其中nh,qh,Φh和Eh是在空间上的某种近似,V-j+1/2表示v(x-j+1/2),V+j-1/2表示v(x+j-1/2).在上式中的(n),(q),(Φ),(E)项被称为数值流通量,是由于数值解在剖分边界的不连续造成的。求解此方程的关键是寻找合适的数值流通量。文献介绍了一种选择交替数值流通量的方法。在本文中,给出了适合本模型方程的数值流通量选取。  随后,作者对耦合方程组的半离散格式进行了稳定性分析和收敛性分析,证明了半离散格式的稳定性,得到了K+1/2阶的收敛阶。本章的最后,作者对半离散格式进行时间分解,给出了全离散格式。本文使用的时间分解主要采用了Shu提出的Runge-Kutta分解方法,该方法在保持了高精度的同时一定程度上提高了稳定性.通过选择合适的数值流通量我们可以证明该格式的稳定性。  第三章给出了方程的有限差分格式:  φkj+1-2φkj+φkj-1/h2=e/εnkj-(nd)j,(0.0.10).  Ekj+1/2=-φkj+1-φkj/h,(0.0.11).  nkj+1-nkj/r+μnk+1j+1/2Ekj+1/2-nk+1j-1/2Ekj-1/2/h=(τ)θnk+1j+1-nk+1j-nk+1j-nk+1j-1/h.(0.0.12).  第四章我们针对该模型应用间断有限元方法和有限差分方法给出了具体的数值试验,列出数值结果,给出了近似解的图形,从图形中我们通过比较可以看出有限差分法和LDG方法求近似解的优劣。  第五章是全文的总结,提出了一些希望。
其他文献
本文主要是对次线性期望框架下的G-正态分布及G-布朗运动进行数值模拟并对所用方法进行数值误差分析。  在金融中的风险度量以及波动率不确定性的研究中,次线性期望的概念是
库存控制自从被提出以来就受到了广大学者的关注。确定性库存下经典的EOQ模型已经不再适用于复杂多变的当今生产消费市场。社会生产力的高速发展、居民的购买力不断上涨带动
摘 要:通过对庄36井区长8储层的岩石学特征、物性特征和孔喉特征等进行了深入分析,表明:庄36井区长8储层属低孔—特低孔、特低渗—超低渗储层,岩性主要为细粒岩屑长石砂岩,发育粒间孔和溶蚀孔,孔喉结构类型属小孔隙、微细喉道型,其储层物性主要受沉积作用和成岩作用的共同影响,造成水下分流河道和河口砂坝物性明显好于河道侧翼。  关键词:鄂尔多斯盆地 庄36井区 储层  一、沉积背景  鄂尔多斯盆地是一个整
数学物理逆问题是现代数学中的一个热点研究领域,研究它的难点在于它的不适定性.在本文中,我们考虑一类经典的逆问题,即热传导问题的侧边值问题(SHE),具体的我们考虑:{uxx=ut,x≥0
全息数据存储系统很久已有大存储容量、短存取时间和高数据传输率的希望。诸如铌酸锂(LiNbO3)之类光折变材料可在适中激光功率下用于全息图记录(图1)。图1 数字全息记录光装
基因芯片可以快速检测成千上万条基因,从而对生物细胞基因水平做出一个全面的认识。近年来,越来越多学者将基因芯片技术应用到肿瘤与癌症的研究中。针对维数高、样本数量少的微
对于整数r,如下定义pr(n):  ∞∑n=0pr(n)qn=∞Πn=1(1-qn)r.研究pr(n)是有意义的.例如,当r=-1时,我们得到经典的拆分恒等式  ∞∑n=0p-1(n)qn=∞Πn1/1-qn.当r=1时,我们得到Euler
自Stefan在十九世纪末对移动边界问题开始研究以来,经过了一百多年各国学者的共同努力,已经取得了许多成果,但是它的应用潜力还是无限的,例如对新型材料的开发和研究.在本篇论文
自从1832年法拉第第一次提出了磁流体动力学(Magnetohydro Dynamic,MHD)问题以来,磁流体动力学的理论研究被国内外诸多学者所研究[1-49]。作为流体力学一个重要的分支,磁流体动
系统发生学是一门研究生物进化规律和物种间遗传关系的学科,利用系统演化树来描述自然界中物种之间遗传关系,得到了较为广泛的关注和研究.用图论的方法研究系统发生学的问题,一