管道流的间断Galerkin方法研究

重力场中的欧拉方程在满足等温条件下保持定常状态。通过对控制方程的源项进行改写,并采用间断Galerkin方法中数值流通量的做法来数值逼近源项积分,同时修改传统的Lax-Friedrichs数值流通量,使其与源项相互平衡,从而使DG格式保持well-balanced的性质。数值实验表明,具有well-balanced性质的Runge-Kutta Discontinuous Galerkin方法可以保持well-balanced性质,可更好的"捕捉"小扰动。     

一维对流扩散方程的四阶精度有限差分法

文章基于Ｈｅｒｍｉｔｅ插值多项式的构造思路，推导出了一维含源扩散方程、扩散反应方程和无源项对流扩散方程的高精度紧致差分格式，并在导出三类特殊方程差分格式的基础上，建立了统一的含源项一维对流扩散方程的差分格式。本文方法是精确的，稳定的，且易应用到其它问题中。数值算例证明了所建立差分格式的有效性。     

非定常对流扩散方程的高阶差分格式

对流扩散方程主要包含对流和扩散两项。在数值计算中,方程中的扩散项一般采用具有优良物理特性和计算精度的中心差分离散格式,而对对流项的处理就稍显困难,处理不当便会产生数值震荡或数值弥散,给数值计算带来困难。针对对流扩散方程,通过引入指数变换将对流扩散方程变为扩散方程,避免对对流项的直接处理。利用四阶紧致差分格式,首先建立三类特殊方程的高精度差分格式,在此基础上建立一维非定常含源对流扩散方程的高阶格式,并进行稳定性分析,所得格式精度高且绝对稳定。数值算例表明了该格式的有效性。     

数值求解一类空间分数阶扩散方程源项系数反问题

数值求解一类空间分数阶扩散方程源项系数反问题.利用函数变换,将源项系数反问题转为对应的定解问题,利用隐式差分格式,求解对应定解问题,然后利用数值积分,求得待定系数函数的数值解,并且证明了隐式差分格式的绝对稳定性.通过数值算例表明,该数值方法具有较高的计算精度.     

应用高精度差分格式求解涡轮叶片热传导方程

为了提高涡轮气热弹耦合计算时叶片热传导的计算精度,给出热传导方程在任意曲线坐标系中的表达形式,并将所得一次项和交叉项作为方程源项,推导出了数值求解三维复杂几何体热传导方程的精度为O(Δτ+Δh4)的ADI紧致格式,研究了源项求解精度为O(Δh4)的紧致格式,并采用Fourier法分析了格式的稳定性。通过数值实验结果验证所得方法的精确性和可靠性,该方法适用于涡轮叶片的热传导的计算,并分析了某涡轮叶片热传导。     

一维对流扩散方程的四阶精度有限差方法

文章基于Ｈｅｒｍｉｔｅ插值多项式的构造思路，推导出了一维含源扩散方程、扩散反应方程和无湖泊项对流扩散方程的高精度紧致差分格式，并在导出三类特殊主程差分格式的基础上，建立了统一的含源项一维对流扩散方程的差分格式。本文方法是精确的。稳定的，且易应用到其它问题中，数值算例证明了所建立差分格式的有效性。     

含源汇项定常对流扩散方程的有限解析差分格式

对于含有源汇项的定常对流扩散方程，本文导出了一种点点精确有限解所差分格式的一般公式，并讨论了源汇项的几种离散形式及对应差分格式的精度,最后利用数值算例验证了本文差分格式的性能．     

稳态不可压流计算方法在稳恒电磁场中的应用

从流体力学和电磁学的相似性出发,揭示了稳态不可压缩流体流动方程、稳恒电场方程和稳恒磁场方程在数学表达上具有相似性,均可用对流项,扩散项和源项来表达.而且在数值求解中也存在类似规律,可利用计算流体力学中对对流项和扩散项的处理方法来处理电磁场方程中的对流项和扩散项.通过分析上(下)风格式与向前(后)差分格式的数学意义和物理意义,提出采用上风格式和下风格式能求解电荷运动问题,而上风格式能求解流体力学问题.     

求解具有复杂地形二维浅水方程的修正HLL格式

采用非结构网格的有限体积方法,对具有复杂地形的二维浅水方程进行数值模拟.采用HLL近似Riemann解计算界面数值通量,基于三角形网格,底坡源项采用简单的斜底模型离散,保证了地形的离散精度,摩阻源项采用全隐方式求解以保证格式的稳定性.采用多维重构及多维限制器的方法获得空间二阶精度的格式,时间离散采用三阶Runge-Kutta法以获得高阶时间精度.为保证格式的和谐性,对经典的HLL格式计算的数值通量中的静水压力项进行了修正.数值计算的结果验证此格式具有良好的高精度捕捉间断的能力,可以应用到地形复杂的二维浅水问题计算中去.     

二维稳态Navier-Stokes方程的有限元方法

Navier-Stokes方程是流体力学中一类重要的数学物理方程,其相关控制方程是非线性的.设计二维Navier-Stokes方程的有限元格式,并实现该算法.对于非线性项采用Newton迭代格式.数值结果表明,该方法不仅具有稳定性,而且具有较好的收敛性.     

QUICK和乘方格式在顶盖驱动方腔流动数值计算中的比较

推导了非均分网络系统中QUICK格式附加源项的表达式，并用QUICK和乘方两种差分格式，分别对不同网格数下的二维、三维顶盖驱动方腔流动进行了数值计算。通过与文献中基准解的比较，考察了两种对流差分格式的数值预测性能。     

Higher-order Accurate and High-resolution Implicit Upwind Finite Volume Scheme for Solving Euler/Reynolds-averaged Navier-Stokes Equations

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｅｍｅｒｇｅｎｃｅｏｆＣＦＤａｂｏｕｔ15ｙｅａｒｓａｇｏｐｒｏｖｉｄｅｄａｍａｊｏｒｉｍｐｅｔｕｓｔｏｓｏｌｖｅｔｈｅＥｕｌｅｒａｎｄＮａｖｉｅｒＳｔｏｋｅｓｅｑｕａｔｉｏｎｓｇｏｖｅｒｎｉｎｇｔｈｅｆｌｏｗｆｉｅｌｄｉｎｅｘｔｅｒｎａｌａｎｄｉｎｔｅｒｎａｌｆｌｏｗｓ.ＭａｊｏｒｐｒｏｇｒｅｓｓｈａｓｂｅｅｎｍａｄｅｉｎｔｈｅｆｉｅｌｄｏｆＣＦＤａｌｇｏｒｉｔｈｍｄｅｖｅｌｏｐｍｅｎｔ,ｓｕｃｈａｓｍｕｌｔｉｇｒｉｄａｃｃｅｌｅｒａｔｉｏｎｔｅｃｈｎｉｑｕｅｓ,ｄｉ…     

基于能量不变二次化法的Cahn-Hilliard方程的数值误差分析

基于能量不变二次化方法,构造了一个求解Cahn-Hilliard方程的线性数值格式,该线性数值格式对非线性项半显式处理,每步迭代相应的半离散化方程只需要求解一个线性方程;证明了该线性数值格式是无条件能量稳定的,而且是唯一可解的;讨论了该线性数值格式在时间方向的误差估计.数值例子表明:该线性数值格式的数值解在时间方向上基本达到二阶精度, 能够有效模拟相位变化过程.     

复杂地形上坝堤溃决洪水演进的非结构有限体积模型

针对复杂地形上坝堤溃决洪水演进的二维数值模拟,综合考虑底坡项、摩阻项和干湿界面对水流的影响,基于浅水二维水动力学理论建立了非结构有限体积模型.围绕传统底坡项处理技术应用于复杂地形时的不足,通过计算底坡项对物质通量和动量通量的影响,并进行相应的分解,提出了集成底坡项的近似Riemann求解器.分析并采用改进的隐式格式解决了摩阻项可能引起的刚性问题.通过求解干湿单元界面处"镜像湿单元"Riemann问题,结合变量重构技术,有效处理了动边界问题.算例研究结果表明,模型稳定健全,和谐性好,具有较好的推广应用价值.     

Rosenau-RLW方程的拟紧致C-N守恒差分格式

本文对Rosenau-RLW方程的初边值问题进行了数值研究,提出了一个隐式拟紧致C-N差分格式,该格式很好地模拟了问题的守恒性质.通过Brouwer不动点定理,本文得到了差分解的存在性,给出了解的先验估计和误差估计,并通过离散能量法分析了该格式的稳定性、二阶收敛性和解的唯一性.数值算例表明,该格式是可行的,且相对于一般的二层格式具有更好的计算精度.     

窃听双向协作网络机会中继选择方案的保密性能分析

对于多中继、窃听双向协作网络,研究使用机会中继选择和模拟网络编码技术来增强窃听系统物理层安全性.具体地,窃听双向协作网络的通信过程共3个时隙.第1和第2时隙,两信源分别发送消息到中继节点;第3时隙,根据最大最小中继选择方案,选取一个成功解码两信源消息的中继节点.使用模拟网络编码技术将两信源消息混合,把混合消息广播给两信...     

数值模拟入射斜激波／平板湍流边界层干扰流动

应用GAO-YONG可压缩湍流模型数值计算了入射斜激波／平板湍流边界层相互干扰现象。计算程序中的对流项、扩散项分别采用AUSM格式和中心差分格式离散,并用多步Runge-Kutta显式时间推进法求解空间离散后的控制方程。计算中包含了无分离流动、初始分离流动以及较大分离流动等多个情况,比较了平板壁面压力、法向平均速度剖面、壁面摩阻系数Cf以及壁面斯坦顿数St等的计算结果与实验值。结果发现：GAO-YONG可压缩湍流模型能够很好地预测入射斜激波／平板湍流边界层相互干扰下的无分离以及小分离流动,对高马赫数下的大分离流动也能得到较合理的结果,但再附点之后的壁面摩阻系数以及斯坦顿数的计算值不够理想。     

一个多维混沌系统的最终有界集和正向不变集及其在同步之中的应用

研究了一个新的不同于Lorenz系统、Chen系统和Lu系统、T系统、Qi系统的五维自治混沌系统,该系统有1个参数,6个非线性乘积项;得到了该系统的一个五维椭球最终有界集和正向不变集界估计;最后,结果应用到混沌同步之中去,数值仿真显示了方案的可行性．     

热源周期振荡条件下一维融化问题的数值研究

借助空间坐标变换,把移动区域模型转化为固定区域模型,通过构造显式和隐式两种有限差分格式求解热源周期振荡条件下的一维融化问题.对所构造的两种差分格式分别研究它们的数值稳定性,比较它们的计算量和计算效率;应用这两种差分格式分别数值模拟融化过程中移动边界的运动及液态介质内温度场的分布.数值实验结果表明,这两种差分格式的数值结果吻合得非常好,而隐式差分格式的计算效率要明显优于显式差分格式.