奇摄动对流-扩散偏微分方程的混合算法

讨论具有一般边界层的奇摄动对流-扩散偏微分方程,这类问题会在边界层附近出现剧烈振荡现象,产生所谓的边界层函数,其解析解无法求出.本文提出混合算法,其主要思想是引入二个过渡点将区域分为粗网格区域、中等网格区域和细网格区域,在这三个网格区域我们采用等步长.在粗网格区域采用Il'in差分格式,在细网格区域采用一般差分格式,在中等网格区域采用渐近解,新方法的总体误差是O(N-1 M-1 ε).混合算法结合了渐近解、数值解和BVT法的优势,是一个实用、有效的算法.     

非线性双曲型方程的混合有限元两层网格算法

针对一类非线性双曲型方程, 利用混合有限元法,构造了1种混合有限元两层网格算法, 给出了两网格方法的误差分析. 结果表明, 当两层网格算法所选取的粗网格和细网格步长满足H=б(h^1/2)时,能获得渐近最优的离散逼近解. 并用数值例子验证了该混合有限元两层网格算法的有效性.     

基于有限元自动生成系统(FEPG)的组合网格算法

针对工程计算中经常出现的局部特性(特别是奇性)问题,和以往解决此类问题的算法的局限性,提出一种基于有限元自动生成系统(FEPG)的组合网格算法.该算法采用两套网格求解,在整个求解区域采用较粗网格,不考虑奇异的影响;而在奇异附近区域采用较细的网格,考虑奇异的影响;整体粗网格求解和局部细网格求解反复迭代,求得最终结果.该算法用于实际工程计算的迭代次数少,与常用的有限元方法所求得的解相符合,为求解大型实际复杂问题提供了一个好的算法和思路.     

渗流方程的渗透率自适应权重网格粗化算法

为了提高渗流方程的计算速度和精度,将渗透率自适应网格技术应用于三维非均匀非稳态渗流方程的网格粗化算法中。对于渗透率或孔隙度变化异常的区域,采用精细网格直接求解其压力分布;而在其他区域,采用不均匀网格粗化的方法计算其压力分布。用自适应权重网格粗化算法计算了三维非均匀非稳态渗流场的压力分布。结果表明,三维非均匀非稳态渗流方程的三维不均匀自适应网格粗化算法的解在渗透率或孔隙度异常区域的压力分布规律非常逼近精细网格算法的解,在其他区域的压力分布规律非常逼近粗化算法的解。与采用精细网格算法相比,其计算速度大大提高。     

渗流方程的渗透率自适应权重网格粗化算法

为了提高渗流方程的计算速度和精度,将渗透率自适应网格技术应用于三维非均匀非稳态渗流方程的网格粗化算法中。对于渗透率或孔隙度变化异常的区域,采用精细网格直接求解其压力分布;而在其他区域,采用不均匀网格粗化的方法计算其压力分布。用自适应权重网格粗化算法计算了三维非均匀非稳态渗流场的压力分布。结果表明,三维非均匀非稳态渗流方程的三维不均匀自适应网格粗化算法的解在渗透率或孔隙度异常区域的压力分布规律非常逼近精细网格算法的解,在其他区域的压力分布规律非常逼近粗化算法的解。与采用精细网格算法相比,其计算速度大大提高。     

Fredholm型泛函积分方程基于■插值的两层网格解法

本文探讨了一类Fredholm型泛函积分方程基于■插值的两层网格解法.利用Banach不动点原理,给出了其解析解存在唯一性的充分条件;在粗网格上采用高效数值积分公式结合配置方法对积分方程进行离散化,并给出了粗网格上的■插值解及收敛性的结果;再次采用不动点迭代的思想,得到了在细网格上以粗网格■插值为初始解的两网格迭代解及收敛性的结果;最后通过数值实验验证了理论分析的有效性和可靠性.     

半线性椭圆方程的一个新的双重网格差分算法

用所提出的双重网格算法研究了半线性椭圆方程,其对粗网格（可以很粗）的非线性解在细网格上进行了几次线性修正．无需求解细网格上的非线性解,且重复算法最后两步,可使解的误差估计达到任意阶精度,并提出了相应的数值算例．     

非线性反应扩散问题两网格混合有限元法的数值分析

对于一类非线性反应扩散问题,给出了一种两网格混合有限元解法.首先在粗网格上求解非线性方程组,然后在细网格上采用了牛顿迭代求解.从数值分析的角度对两网格混合有限元算法进行了研究.数值算例结果表明,与混合有限元方法相比,两网格混合有限元方法在不降低解的精度阶数的条件下,提高了计算速度.     

一类半线性椭圆方程的二重网格差分算法

应用二重网格差分算法处理了一类半线性椭圆问题。无需求细网格上的非线性解,对粗网格(可以很粗)上的数值解在细网格上进行几次线性修正即可,且重复算法的最后一步可以按粗网格步长任意阶地逼近细网格上的非线性解。算法提高了计算效率但不降低精度,有数值算例加以验证。     

改进格式的组合网格法

针对组合网格法在求解特殊类型的工程问题时求得数值解精度低的缺陷,提出一个改进计算格式的组合网格法.粗细两套网格都是在各自区域上单独剖分,两套网格互不影响,在粗网格上达到粗网格的精度,在细网格上达到细网格的精度.该方法的网格剖分为粗网格和细网格,能够大大提高数值解的精确度.     

二维磁场积分方程奇异点处理的有效方法

基于积分方程的矩量法求解电磁散射问题需要奇异积分计算,而且奇异阻抗矩阵的计算是影响矩量法计算精度的重要因素之一.论文将基于二维磁场积分方程的脉冲函数作为基函数、δ函数作为检验函数,对奇异矩阵元素的计算进行特别处理,将对角元素的计算分为两个子部分,每个子部分的贡献采用非奇异传统方法计算,于是对角元素的值为奇异值与两个子部分的贡献之和.数值结果表明:该方法用于曲线散射体求解具有有效性和正确性.     

对流扩散问题非均匀网格上的部分半粗化多重网格方法

结合非均匀网格上的 HOC 格式与部分半粗化的多重网格方法对具有边界层的2维对流扩散问题进行了求解,并基于面积率构造了部分半粗化多重网格方法的插值算子和限制算子。数值实验表明:对于只需要在1个方向进行网格加密的边界层问题,基于部分半粗化的网格分布策略及多重网格算法可以大大减少无边界层方向的网格数,从而较完全粗化的网格分布策略及多重网格算法具有更高的计算精度和求解效率。     

Allen-Cahn型相场模型的移动网格方法

针对2种流体的Alleh-Cahn型相场模型提出了一种自适应移动网格方法.移动网格方法包括网格重分布和偏微分方程求解2个相对独立的部分.通过求解一个类似于Poisson方程的偏微分方程组获得网格重分布,大量的网格点聚集在2种流体的界面附近从而提高了分辨率,而在其他的区域则仅有比较稀疏的网格点.Allen-Cahn模型用...     

不确定广义时滞系统的时滞相关鲁棒镇定新判据

针对线性不确定广义时滞系统,给出了新的时滞相关鲁棒稳定性判据和镇定方法.假定参数是可能时变的范数有界不确定,通过选取一个新的Lyapunov-Krasovskii泛函将整个时滞区间分为2段.每段区间上定义了不同的能量函数,给出了状态反馈控制器设计算法,保证相应的闭环系统对所有容许的不确定性参数是正则、因果和鲁棒稳定的....     

一维非线性对流扩散方程特征有限元的两重网格算法

针对一维非线性对流扩散方程，构造了特征有限元两重网格算法。该算法只需要在粗网格上进行非线性迭代运算，而在所需要求解的细网格上进行一次线性运算即可。对于非线性对流占优扩散方程，不仅可以消除因对流占优项引起的数值振荡现象，更重要的是可以加快收敛速度、提高计算效率。误差估计表明，只要选取粗网格步长与细网格步长的平方根同阶，就可以使两重网格解与有限元解保持同样的计算精度。     

基于五次B样条的子区间法

重构细分划了时域,细分划拓展了0次B样条醒的定义,对高次B样条的递推式进行了拓展,获得了细分划拓展的均匀分划的分段式五次B样条函数,因而拓展了展开定理,构造了五次B样条基函数．该基函数与现有的五次B样奈基函数相比,表征的物理概念更清晰和简洁．基于五次B样条基函数,提出和推导了结构动力响应研究中的位移元子区间法和子区间法的嵌套方法递推格式．通过精度比较,得出子区间法的嵌套方法要优越于位移元子区间法．     

自适应网格下齐次奇异摄动问题的一致收敛性分析

利用自适应移动网格方法求解齐次奇异摄动边值问题,通过常数与解的一阶导数幂的线性组合构造控制函数来进行网格自适应,分析齐次奇异摄动边值问题在此自适应非均匀网格上的收敛性.利用极值原理证明离散问题解的存在唯一性.对离散问题的数值解及其分段线性插值进行误差估计,分别得到一个与ε无关的一阶误差界.有效地解决齐次奇异摄动方程难以...     

简化的WENO重构在h自适应RKDG方法中的应用

将简化的WENO重构方法应用到了一个基于坏单元指示子的h自适应RKDG算法,并从数值上证明了这个算法能有效提高间断附近的数值解质量。数值试验还表明这个重构方法对一维h自适应网格也能很好地消除数值振荡,同时还能保持光滑区域的精度。     

含间断系数弹性结构的多重网格方法与数值模拟

提出了求解含间断系数弹性力学问题的界面保持粗化多重网格方法,该粗化方法在选取粗网格节点时保证在每一个网格层上能保持界面的实际形状,同时可以捕获位移解函数沿界面处导函数的不连续行为,这样只需要构造简单的插值算子,并选取点块Gauss-Seidel作磨光迭代,就能达到理想的多重网格收敛效率.数值实验结果表明,这种界面保持粗化多重网格方法的收敛性不依赖于网格规模及间断系数的大小,具有很好的数值稳定性.