Stokes方程的一种预处理方法

来源 :数值计算与计算机应用 | 被引量 : 1次 | 上传用户:luzihao009
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
本文研究了Stokes方程基于混合有限元离散的预处理方法.基于三角形剖分上速度场的连续P_2元以及压力场的P_0元,我们采用局部子区域逆矩阵建立了预处理子.通过分析预处理矩阵的特征值分布,我们证明了该预处理子大大改善原问题系数矩阵的条件数.本文给出了相应的数值结果.
其他文献
本文基于微元电阻和积分中值定理导出了仅有3个待定参数的超深亚微米MOSFET的漏/源电阻模型,而这3个参数可以用多元线性回归方法得到.论文用数值模拟数据,拟合了衬底掺杂是1×10~(15)~1×10~(16)cm~(-3)、沟道长度是45~2000 nm的平面n MOSFET和n LDD-MOSFET的漏/源电阻模型的待定参数,得到了它们的半经验公式,最大误差仅有9.5%.漏/源电阻模型和它的半经
期刊
回旋管基于相对论电子回旋脉塞原理工作,可在毫米波–太赫兹频段产生高功率相干电磁辐射.脉冲磁体体积小巧,易产生强磁场,相比于超导磁体可极大地降低回旋管系统体积并提供电子回旋所需的强磁场,使回旋管在兼具小型化特性的同时其辐射频率有望突破1 THz.传统脉冲磁体回旋管通常采用开放腔式电路高Q值单频点工作,需要近似稳恒的磁场,使得脉冲磁体仅在峰值场强附近可用于支持回旋管工作,导致回旋管单脉冲工作时间极短,
期刊
本文提出了一种基于分层复杂性度量的形状描述方法 (简称HCMD),以应用于一般的形状图像检索任务.该方法属于基于区域的形状描述方法,不需要专门提取形状的边界完成特征抽取任务.本文的主要贡献是提出了一种分层的形状描述框架,基于该框架,形状沿各个方向被迭代地切割成小的区块,每一层级的形状区块被施以各种几何特性的度量,以刻画其形状复杂性.这种分层抽取形状复杂性特征的描述机制使得HCMD具有由粗到细的形状
期刊
本文讨论非PH(Pythagorean hodograph)曲线的一类五次OR曲线(curves with rational offsets)的构造方法.OR曲线是具有有理形式的等距线的一类参数曲线,在CAD中有着广泛的应用.本文采用参数曲线的复数表示形式,根据导矢曲线的因式分解的不同,将非PH曲线的五次OR曲线分为两种类型,并分别给出这两类曲线的构造方法.在给定C1连续的初始条件下,通过指定一个
期刊
~~
期刊
本栏目新闻资讯从2016年10期起以全新模式呈现给读者。以国内、国际两个视点出发,多领域体会行业发展的魅力所在,进而思考行业激烈竞争态势的发展方向。丰富内容敬请关注和浏览,欢迎扫码参与!(温馨提示:此码扫描仅适用于安卓系统手机,苹果系统版即将推出,敬请期待!)
期刊
信息传播算法求解可满足问题时有良好的有效性,使得难解区域变窄.然而,信息传播算法不总有效,常表现为不收敛.对于这种现象,至今缺少系统的理论解释.调查传播(survey propagation,SP)算法是最为有效的信息传播算法,对SP算法的收敛性研究是设计其他信息传播算法的重要基础,并为信息传播算法的广泛应用提供理论依据.为了分析SP算法的收敛性,通过对消息更新方程取双曲正切,将消息取值从[0,1
期刊
本文采用Legendre Galerkin-Chebyshev配置方法求解广义空间分数阶Burgers方程.该方法基于Legendre Galerkin变分形式,但是非线性项与右端源项采用Chebyshev-Gauss插值逼近.首先,通过在空间方向采用Legendre Galerkin-Chebyshev配置方法离散,时间方向采用leap-frog/Crank-Nicolson格式离散,得到了方程
期刊
研究了L(0,1)模态导波在管道弯头处的反射特性,提出了一种新的基于导波模态置信准则计算的L(0,1)模态导波弯头反射分析方法.使用数值模拟方法研究了L(0,1)模态导波在弯头处的模态转换,并分析了导波激励频率、弯曲半径和弯曲角度3种因素对L(0,1)模态导波反射的影响;利用导波模态置信准则表示导波弯头反射程度,反射程度随频率和弯曲半径变化趋势与数值模拟结果进行了对比,最后进行了实验验证.结果表明
期刊
二次正交Arnoldi方法的存储量小并能保持与标准Arnoldi方法类似的数值稳定性和收敛性,因此成为求解二次特征值问题的重要方法.算法的运行过程中计算量和存储量会不断增加,将算法进行重新启动是算法在实际使用中的必然需求,该方法的特殊分解形式对算法的重启提出了新要求.本文分析了该方法所形成子空间的性质和重启时子空间应具有的形式和性质,提出了一种能够保持算法特殊子空间结构且简便易实现的重启方法.在此
期刊