两类波方程基于Birkhoff插值的谱配置法

来源 :哈尔滨工业大学 | 被引量 : 0次 | 上传用户:hlxcun3e5
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
波方程被广泛地用来刻画各种波传播现象。本文将研究其中两类典型的方程,即Schr?dinger方程和Korteweg-de Vries方程。这两种方程出现在量子力学、非线性光学、等离子体、超导体、晶体学等领域中。谱配置法因其指数收敛性,是一种被广泛采用的数值方法。然而,传统的基于多项式基函数的谱配置法产生的系数矩阵是稠密和病态的。本文将基于Birkhoff插值问题,构造新的非多项式函数系,并以此为基函数建立谱配置法。该谱配置法对一类边值问题产生的系数矩阵是良态的,并且精确地满足边界条件。另一方面,在很多情况下物理区域是无界的。当考虑无界区域上波方程的数值解法时,由于计算资源的有限性,需要将计算区域限定为有界区域。因此需要在有界区域上施加合适的人工边界条件。如果施加的边界条件不合适,向外传播的波在边界上会被反射回计算区域,这与无界区域上的波方程解的行为不符。因此必须使这种虚假的反射尽可能的小。如果一种人工边界条件使得在边界处没有虚假反射,则称之为透明边界条件。针对有界区间上的非线性Schr?dinger方程,基于二阶广义Birkhoff插值,构造了一组非多项式基函数。为了讨论该非多项式基函数的插值误差,首先给出了Jacobi-Gauss-Lobatto多项式插值算子在非一致加权Sobolev空间中的稳定性估计和收敛性分析。基于此估计,建立了非多项式基函数的Birkhoff插值在加权Sobolev空间中的误差估计。考虑一维无界区间上的线性Schr?dinger方程,在进行时间离散后,给出了时间半离散的透明边界条件,并证明了在此边界条件下,时间半离散格式的稳定性。最后对时间半离散格式应用基于非多项式基函数的谱配置法进行空间离散,得到全离散格式。全离散格式精确地满足边界条件,并且通过适当地选取基函数中的自由参数,所得到的全离散格式是显式的。最后通过数值实验,验证了格式的高精确性和对反射波的透明性。针对一维无界区间上的线性Korteweg-de Vries方程,进行时间离散后,原无界区间上的问题转化为等价的内问题和外问题。外问题可以由Z变换解析地求解,从而得到内问题的透明边界条件。建立了带有透明边界条件的时间半离散格式的稳定性。对Korteweg-de Vries方程,透明边界条件中的逆Z变换没有显式表达式,给出了一种稳定且精确计算逆Z变换的算法。然后考虑Korteweg-de Vries方程时间半离散格式的空间离散,基于三阶广义Birkhoff插值,构造了另一组非多项式基函数。建立了该基函数在非一致加权Sobolev空间中的插值误差估计。然后以此非多项式函数为基底,构造了Korteweg-de Vries方程的谱配置法全离散格式。该格式保持了边界条件对反射波的透明性。同时,在基函数的构造中包含了两个自由参数,通过适当地选取这两个自由参数,最终的全离散格式是显式的。数值实验也验证了所构造格式的有效性和高精确性。
其他文献
含有驱动冗余的过约束并联机构具有承载能力大、各向同性良好等特点,但其同时含有主动过约束和被动过约束,受力分析相当复杂。本文以4-PRRR冗余驱动过约束并联机构为对象对该
目的探讨免疫抑制剂治疗慢性肾病并发重症肺炎的临床特点及治疗措施。方法回顾分析18例免疫抑制剂治疗慢性肾病并发重症肺炎的临床资料,分析其致病原因、病原菌特点及耐药情
为解决泡沫铝填充管短流程近终成形问题,本文结合连铸工艺中间包底吹氩技术与挤压工艺中汇流模具设计理念,提出了一种底吹反挤复合工艺,并且基于相似性原理,自主搭建了成形容
目的总结连续性肾脏替代治疗(CRRT)重症挤压综合征患者的效果。方法对18例重症挤压综合征导致的急性肾功能衰竭患者进行了CRRT治疗及综合治疗。结果 CRRT治疗后患者症状均有不
目的比较单、双侧椎弓根入路经皮椎体后凸成形术(PKP)治疗老年新鲜骨质疏松性椎体压缩骨折(OVCFs)的临床效果。方法对胸腰椎椎体压缩性骨折38例(47椎)应用PKP治疗,包括单球囊
磁控溅射是一种物理气相沉积成膜技术,采用纯铝靶材进行磁控溅射具有成本低、工艺稳定等优点,是薄膜晶体管液晶显示器(TFT-LCD)金属薄膜的主要成膜方式之一。但纯铝的热膨胀系数远大于玻璃,在受热时铝薄膜的体积会发生膨胀,当铝膜的膨胀在玻璃一侧受到限制时,限制膨胀产生的压缩应力会通过原子扩散的方式释放,沿着晶粒边界横向扩散,形成小丘(hilock)。在TFT-LCD生产中,小丘的产生会引起栅极与源极、
近年来,随着CMOS工艺的不断演进,先进的工艺节点推动着半导体电路向更高速、更高能效的方向不断发展。模数转换器(Analog to Digital Converter,ADC)作为连接自然界和人类数