双曲型线性方程三阶和四阶TVD格式的新构造

利用Taylor级数理论和总变差减小(TVD)格式的充分条件构造了时间二阶、空间五点三阶和四阶新TVD格式．给出了新TVD格式与传统TVD格式及近期建立的二阶新TVD格式用于线性双曲型方程的计算结果，表明本文新格式特别是四阶TVD格式具有比二阶新TVD格式和传统TVD格式峰值衰减更慢、间断更陡，而计算工作量具有与传统二阶TVD格式相当的良好数值性能．     

一维WENO格式及其在标量守恒方程数值计算中的应用

研究高阶精度加权本质无振荡(WENO)格式及其在标量守恒律方程中的应用.应用五阶WENO空间离散格式和三阶TVD Runger-kutta时间离散格式对一维标量守恒方程以及二维标量守恒方程进行了数值模拟.数值结果表明该格式具有高精度性和本质无振荡性.     

一种求解时间关联型缓坡方程的数值方法

以含底摩阻能量耗散项的时间关联型缓坡方程为控制方程,建立了一种具有二阶精度的数值离散格式.该格式对时间导数使用Euler预测-校正格式离散;对空间导数使用中心差分格式离散.基于统一边界条件表达式,对边界条件进行处理.数值解与物理模型实验值吻合较好,表明数值模拟模型可以有效地模拟波面随时间和空间的变化以及波高的空间分布.     

二维无结构网格浅水波方程的高阶非振荡有限体积法

基于二维浅水波方程,对无结构网格给出了一种三阶精度非振荡有限体积方法,方法的主要思想是先对每一个三角形单元采用最小二乘的思想构造一个二次插值多项式,而在计算交界面的流通量时采用了两点高斯公式以保证格式的整体精度,时间离散采用三阶TVD Runge-Kutta方法。最后用该格式对圆柱形溃坝问题和倾斜水跃问题进行了数值试验,得到了满意的结果。     

和谐WAF格式在带干河床浅水波方程中的应用

采用二阶精度和谐的加权平均通量（WAF）格式对浅水波方程中的间断解和干湿边界问题进行了研究．用基于HLL求解器的WAF格式近似数值通量,中心差分格式离散地形源项,在干湿边界处重新定义数值河床,保证格式是和谐和扩展和谐的．最后模拟了带激波和干湿边界的浅水问题,精确地捕捉到了激波和干湿边界面,验证了该格式具有高分辨率、无振荡以及正确模拟干湿边界的能力．     

求解扩散方程的Pade′逼近方法

对扩散方程提出了精度为O(t3+h2)的差分格式,首先对空间变量中心差分格式离散,所得到常微分方程组利用指数函数的Pade[2/1]逼近,得到空间二阶时间三阶精度的两层绝对稳定的隐式差分格式,并讨论了稳定性.数值结果与Crank-Nicholson格式进行比较,数值结果表明,该格式不但有效地解决初始边界条件间断问题,而且适合于大时间步长问题.     

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

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

相场晶体方程的一个高精度能量稳定数值格式

针对具有周期边界条件的相场晶体方程,本文提出了一个具有能量稳定性的高精度数值格式.该格式基于方程的能量泛函结构,在空间上采用Fourier拟谱逼近,在时间上进行三阶精度的向后差分离散,并在格式中增加Douglas-Dupont正则项,以保证格式的能量稳定性.本文证明了数值解的存在唯一性及数值格式的能量稳定性.数值算例验证了算法的高精度和稳定性.     

一维磁流体力学方程的三阶非交错中心格式

基于交错网格上的三阶无振荡重构和将交错网格单元平均值转化为非交错网格单元平均值,构造了一维理想磁流体力学方程的一类三阶非交错无振荡中心差分格式.给出两个典型的数值算例,验证了格式具有高分辨率和无振荡特性等优点.     

移动网格上解二维Euler方程的WENO型有限体积法

对二维Euler方程在基于弹簧技术的移动非结构三角形网格上给出了一种WENO型有限体积格式,方法的主要思想是先对每一个三角形单元构造插值多项式,而在计算交界面的流通量采用两点高斯积分公式以保证格式的整体精度,时间离散采用三阶TVD Runge-Kutta方法.最后用该格式对一些典型算例进行了数值模拟,结果表明该方法计算速度快,对间断有很好的分辨能力.     

线接触弹流润滑问题的贴体坐标变换和交错网格解法

通过贴体坐标变换把带这界的线接触弹流润滑问题化为定边界问题。在分析了目前常用解法的二阶格式失隐的原因后，提出了一种交错网格二阶格式。算例表明，本二阶格式不令比一阶格式的精度有显著提高，而且具有良好的数值稳定性。     

燃气射流CFD计算软件初版本的研究

目的 发展一个燃气射流数值计算软件。方法 用时间相关法求解全ＮＳ方程及Ｒｅｙｎｏｌｄｓ平均应力方程,方程的时间项采用了ＬＵ分解及４步Ｒｕｎｇｅ－Ｋｕｔｔａ方法,对流项采用了Ｒｏｅ差分分裂基础上的３阶ＴＶＤ级２阶、３阶ＥＮＯ格式。结果与结论软件可应用于轴对称,三维燃气射流计算,格式具有高分辩和健壮稳定性,能够捕获复杂波系,计算出射流中波系与剪切层的相互作用和湍流特征。     

基于非结构化网格有限体积的LBM

提出了一种基于非结构化网格有限体积的LBM.采用Cell-vertex有限体积法离散控制方程.该方法在时间上采用伪、实二时间步,其中伪时间步采用向前差分,实时间步采用二阶向后差分方法;空间上采用edge-based通量计算方法,采用高阶TVD格式计算控制体边界通量.离散后的控制方程采用隐式迭代,控制变量采用五层二阶Runge-Kutta方法求解.二维同心圆环内圆柱间Couette流与顶盖驱动方腔流的数值结果显示该方法为一种有效求解不规则边界流体动力特性的实用工具.     

多物质流体动力学界面追踪的Particle Level Set方法

为验证Particl Level Set方法界面处理的高精度,将Particle Level Set方法与流体动力学控制方程进行了耦合,用5阶WENO格式进行空间离散,采用3阶TVD-Ringe-Kutta方法离散时间导数,并用C++语言编制了程序. 运用该程序对激波与气泡相互作用问题进行了数值模拟. 与不同网格尺寸下的Level Set方法对比,结果表明Particle Level Set方法比Level Set方法界面处理精度更高.     

求解Hamilton-Jacobi方程的一类高精度差分格式

构造了一、二维非线性Hamilton-Jacobi方程的一类新的高精度高分辨率差分格式.首先将计算区域划分为互不重叠的子单元,再根据格式的精度要求分割子单元为细小于单元,其次通过子单元上各个细小子单元节点的函数值构造空间导数的高阶插值逼近,为避免由此产生的数值振荡,对空间导数在各节点左右侧的值进行TVD/TVB校正,利用高阶Runge-Kutta TVD时间离散方法得到一维Hamilton-Jacobi方程的高阶全离散格式并推广到二维情况,最后给出了几个典型的数值算例,验证了格式具有计算简单、高分辨间断导数、无振荡等特性.     

气-液两相流界面迁移现象的数值模拟研究

借助于LevelSet函数,建立了气-液两相的统一控制方程组,并在交错网格中进行离散.用两种格式,即Superbee TVD格式和5阶WENO格式求解LevelSet函数的输运方程,用SIMPLER算法的思想对主流场控制方程的求解方法进行改进.数值实验结果表明,在求解LevelSet的控制方程时,5阶WENO方法比Superbee TVD格式的结果更准确;用改进的数值算法可成功实现对密度比大于1000/1的气-液两相流界面迁移问题的数值模拟.对几种典型大密度比气-液两相流问题的计算结果与实际问题的物理规律完全一致,验证了该方法的有效性和可靠性.     

三维非定常不可压N-S方程一个高精度求解格式

对三维、非定常、不可压Navier-Stokes方程提出了一种新的数值计算方法.在时间离散上应用3阶精度混合显-隐相结合的分裂格式,空间离散在x及y轴向采用非等间距网格的紧致有限差分格式与z向应用Fourier谱展开相杂交的数值方法逼近.经平板边界层流的验证表明,该算法具有计算精度高、稳定性好、收敛速度快等特点.同时也研究了三维、非定常流体运动下游边界问题,提出了无反射出流边界条件,以减少在有限计算区域内人工出流边界反射引起的数值误差,保证直接数值模拟的精度和准确性.该算法的提出对于求解边界层、射流及混合层等流动中的转捩与演化问题具有重要的理论意义.     

低耗散TVD格式的构造

目的　为进一步提高双曲型方程数值解的精度,减少数值耗散,同时保证稳定性,处理有强间断的流动问题; 方法　从总变差减小（ＴＶＤ）定义出发,分析保证数值解ＴＶＤ的条件,找出满足研究目的的可能性和途径; 结果　首次提出和证明双曲型方程数值解ＴＶＤ的充分必要条件,并据此简化导出适用于时间显示格式的ＴＶＤ充分条件,并据此构造了一种低耗散ＴＶＤ 格式; 结论　进行的数值实验证明,高分辨格式可以进一步减小数值耗散,提高数值解精度;     

广义对称正则长波方程的一个拟紧致守恒差分算法

对一类广义对称正则长波(GSRLW)方程的初边值问题进行了数值研究,提出了一个三层拟紧致平均隐式差分格式,格式模拟了初边值问题的守恒性质,得到了差分解的存在唯一性,并利用离散泛函分析方法分析了该格式的二阶收敛性与稳定性.数值结果表明,该格式的精度明显好于一般的二阶格式.