三维两相渗流驱动问题迎风区域分裂显隐差分法

对三维两相渗流驱动问题提出了两种迎风区域分裂显隐差分格式．压力方程采用了七点差分格式,为了能达到实际并行计算的要求,对饱和度方程采用了迎风区域分裂差分法,内边界处和各子区域分别对应显隐格式．得到了离散l2模收敛性分析,最后给出数值试验,支撑了理论分析结果．     

一种健壮的低耗散通量分裂格式

随着计算流体力学的快速发展,设计精确、高效并且健壮的数值格式变得尤为重要.通过对3种流行的通量分裂方法(AUSM、Zha-Bilgen和Toro-Vázquez)的对流通量和压力通量进行特征分析,构造了一种简单、低耗散并且健壮的通量分裂格式(命名为R-ZB格式).采用Zha-Bilgen分裂方法将Euler方程的通量分裂成对流通量和压力通量,其中对流通量采用迎风方法来计算,压力通量采用低耗散的HLL格式来计算,从而克服了原始的HLL格式不能精确分辨接触间断的缺点.数值实验表明,该文给出的R-ZB格式不仅保留了原始Zha-Bilgen格式简单高效、能够精确分辨接触间断等优点,而且具有更好的健壮性,在计算二维问题时不会出现数值激波不稳定现象.     

一类奇异摄动边值问题的一致收敛迎风差分方法

对于一类奇异摄动边值问题,基于等分布弧长控制函数构建网格,提出了一种迎风差分方法.利用先验截断误差估计,基于离散比较原理和障碍函数技巧,证明了该方法得到的逼近解在最大模下是不依赖于摄动参数且一阶一致收敛的.收敛性分析是在整个区域上进行的,不需要对区域进行子区域的划分.为了验证理论分析,给出了数值实验结果.     

椭圆型方程的重叠型区域分裂混合元方法

本文研究椭圆型方程的重叠型区域分解混合元方法，对第一边值和第二边值问题，分别给出了离散形式的区域分解混合元格式；证明了区域分裂格式解的存在唯一性和算法的收敛性，并给出数值算例．     

区域分裂法求非线性抛物方程的扩张混合元解

针对非线性抛物方程,给出了全离散的扩张混合元格式,利用一个建立在非重叠型区域分裂技巧上的并行迭代法求解了最后的非线性代数方程组,证明了迭代法的收敛性并给出了最优阶的误差估计.     

关于二阶椭圆方程区域分裂法的最优预处理

§0.引言 区域分裂是与微分方程数值解的并行计算的数学基础密切相关的,预处理共轭梯度法是区域分裂的一个主要途径,寻找好的预处理子是关键问题,本文给出一个较一般性的方法,预处理过程包括一个整体小规模问题和若干个独立的局部子问题,整体问题和局部问题的选取均有极大的任意性,预处理条件数的估计是由整体问题和局部问题的一些特     

非饱和水流问题的迎风差分法及其数值模拟

本文考察了非饱和水流问题模型方程的守恒型迎风差分法.我们基于有限体积方法建立的非饱和流动的守恒形式,分别提出了一阶和二阶迎风差分格式,并对差分格式进行了误差估计,给出了收敛性定理.最后,数值模拟验证了计算格式的有效性.     

沉浮,俯仰及控制面振动引起的跨音速非定常气动载荷计算

本文将文献[9]提出改进的通量分裂方法,应用于随时间变化的贴体网格中,建立了可用于求解非定常Euler方程的通量分裂方法.该方法是以连续的特征值分离为基础,它具有方法简单,便于推广使用的特点.同时克服了Steger-Warming通量分裂方法存在的问题.对通量分裂后的Euler方程.利用MUSCL型迎风差分建立了具有二阶精度的有限体积方程.文中以NACA64A—10翼型为例,对其在跨音速流场中进行沉浮、俯仰及带有振动控制面引起的非定常气动载荷进行了计算.部分计算结果与相应的实验结果进行了比较,吻合良好     

矩阵多分裂迭代法的收敛条件

本文给出了求解非奇异线性方程组的矩阵多分裂并行迭代法的一些新的收敛结果.当系数矩阵单调和多分裂序列为弱正则分裂时,得到了几个与已有的收敛准则等价的条件,并且证明了异步迭代法在较弱条件下的收敛性.对于同步迭代,给出了与异步迭代不同且较为宽松的收敛条件.     

求解粘性流体和热迁移联立方程的迎风局部微分求积法

微分求积方法(DQM)已成功地应用于数值求解流体力学中的许多问题．但是已有的工作大多限于正规区域的流动问题,同时缺少用迎风机制来描述流体流动的对流特性．该文对一个不规则区域中的不可压缩层流和热迁移的耦合问题给出了一种具有迎风机制的局部微分求积方法,对通过边界和坐标不平行的收缩管道中的流体,只用少数网格点得到了比较好的数值解．和有限差分方法（FDM）相比较,这一方法具有计算工作量少、存储量小和收敛性好等优点．     

3阶非单谷Feigenbaum映射的拟极限集

一个从闭区间到自身的连续映射被称为3阶非单谷Feigenbaum映射,如果它是函数方程f~3(λx)=λf(x)的解．本文讨论了3阶非单谷Feigenbaum映射的拟极限集及其Hausdorff维数．3阶非单谷Feigenbaum映射必然产生混沌,混沌的产生使得拟极限集的存在性问题复杂化．文中采用分形几何中的知识方法证明了此类映射的拟极限集的存在性,并相应的对其Hausdorff维数作出了估计．最后给了一个具体的例子,说明确实存在这样的3阶非单谷Feigenbaum映射．     

定常不可压Navier—Stokes方程迎风格式

本文利用齐次定解条件对定常不可压Navier—Stokes方程的非线性项进行处理,给出了相应的一种迎风Galerkin有限元算法;针对这种迎风Galerkin有限元算法,在迎风参数满足一定条件下,利用其三项式具有的一些很好性质,更简单地证明了该问题解的存在唯一性。     

用摄动配置方法求解含时薛定谔方程

将摄动配置方法应用到含时薛定谔方程,在计算实现的基础上结合摄动配置的特征提出了一类新的数值积分方法,并给出了一个2级2阶和一个3级4阶的辛摄动配置方法对含时薛定谔方程的数值算例.为了检验新的数值积分方法,我们还给出了与两个辛摄动配置格式在理论上等价的辛龙格-库塔方法以及同阶的非辛方法的数值模拟.展示了一些数值结果,并给出了一些分析.     

16阶二次幻方的膨胀构造法

2017年詹森构造了6个异基因的8阶二次幻方兼完美幻方,根据它们的特殊性质,创立用一个4阶矩阵代替原有元素的膨胀法,构造出16阶二次幻方兼完美幻方;并对另外2个具有相似性质的8阶二次幻方,也通过膨胀法构造出了16阶二次幻方.     

非线性算子方程的泰勒展式算法

本文的目的是给出一种解Ｈｉｌｂｅｒｔ空间中非线性方程的ｋ阶泰勒展式算法（ｋ１）．标准Ｇａｌｅｒｋｉｎ方法可以看作１阶泰勒展式算法，而最优非线性Ｇａｌｅｒｋｉｎ方法可视为２阶泰勒展式算法．我们应用这种算法于定常的Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程的数值逼近．在一定情景下，最优非线性Ｇａｌｅｒｋｉｎ方法提供比标准Ｇａｌｅｒｋｉｎ方法和非线性Ｇａｌｅｒｋｉｎ方法更高阶的收敛速度．     

激波捕捉差分方法研究

在迎风型格式和矢通量分裂技术的基础之上,对捕捉激波方法进行一种新的尝试．该方法首先对原始格式在特征方向上进行投影,然后用限制器对这些特征分量的变化幅值进行限制以抑止非物理波动,最后再把它转换成守恒形式,得到了基本上无振荡的激波捕捉格式．用该方法对两种迎风显示格式(二阶和三阶)和3种迎风紧致格式(三阶、五阶和七阶)进行处理,并在一维和二维的情况下进行了应用测试．通过与高阶WENO、MP、Compact-WENO等格式的比较,表明该方法在光滑捕捉激波的前提下仍有较高精度和分辨率．     

Uniformly convergent finite difference methods for singularly perturbed problems with turning points

In this work one-dimensional singular perturbation problems with turning points are considered. To resolve these problems numerically we consider a family of finite difference schemes, which includes classical methods in literature, such as the upwind method, the Samarskii method and exponential fitting type methods. Once the uniform convergence of the upwind method on irregular meshes has been established, the same property is easily shown on all the elements of the family.Work supported by a grant of the Diputación General de Aragón.     

THE UPWIND FINITE ELEMENT SCHEME AND MAXIMUM PRINCIPLE FOR NONLINEAR CONVECTION-DIFFUSION PROBLEM

In this paper, a kind of partial upwind finite element scheme is studied for twodimensional nonlinear convection-diffusion problem. Nonlinear convection term approximated by partial upwind finite element method considered over a mesh dual to the triangular grid, whereas the nonlinear diffusion term approximated by Galerkin method. A linearized partial upwind finite element scheme and a higher order accuracy scheme are constructed respectively. It is shown that the numerical solutions of these schemes preserve discrete maximum principle. The convergence and error estimate are also given for both schemes under some assumptions. The numerical results show that these partial upwind finite element scheme are feasible and accurate.     

UPWIND SPLITTING SCHEME FOR CONVECTION-DIFFUSION EQUATIONS

1　ＩｎｔｒｏｄｕｃｔｉｏｎＣｏｎｖｅｃｔｉｏｎ ｄｉｆｆｕｓｉｏｎｅｑｕａｔｉｏｎｉｓａｆｕｎｄａｍｅｎｔａｌｅｑｕａｔｉｏｎｄｅｓｃｒｉｂｉｎｇｔｈｅｐｒｏｃｅｓｓｏｆｆｌｕｉｄｔｒａｎｓ ｆｅｒ,ｆｏｒｅｘａｍｐｌｅ ,ｕｎｄｅｒｇｒｏｕｎｄｗａｔｅｒｃｏｎｔａｍｉｎａｔｉｏｎ ,ｄｉｓｐｌａｃｅｍｅｎｔｉｎｐｏｒｏｕｓｍｅｄｉａ[1,2 ] ,ａｎｄｓｏｏｎ .Ｆｏｒｏｖｅｒｃｏｍｉｎｇｔｈｅｎｕｍｅｒｉｃａｌｉｎｓｔａｂｉｌｉｔｙｏｆｆｉｎｉｔｅｄｉｆｆｅｒｅｎｃｅｍｅｔｈｏｄｏｒｆｉｎｉｔｅｅｌｅｍｅ…