激波与不同形状界面作用的数值模拟

该文数值研究了激波与不同形状界面相互作用下的Richtmyer-Meshkov不稳定性发展．数值方法采用自适应的非结构四边形网格来精细刻画界面的演变．在激波Ma=1．2,界面内为sR气体的条件下,通过数值模拟激波与矩形、三角形以及椭圆界面相互作用的过程,分析并对比了这3种形状界面的波系、涡量和界面演变．结果表明,在激波作用下,不同形状流体界面的波系结构是不同的,主要表现在波系的形状以及发展两方面．此外,各种形状界面涡量的产生、分布和大小也有所不同,进而导致各种形状界面演变过程也有一定差异,因此不同形状的流体界面引起的Richtmyer-Meshkov不稳定性有不同的机理．     

三维飞船绕流问题的高精度迎风紧致差分格式

从迎风紧致逼近出发,给出了一个求解三维可压Navier-Stokes方程的一种高精度的数值方法,利用Steger-Warming的通量分裂技术将守恒型方程中的流通向量分裂成2部分,在此基础上据风向构造逼近于无粘项的三阶迎风紧致有限差分格式,对方程中的粘性部分采用通常的二阶差分逼近,所建立的差分格式被用来数值求解了三维飞船绕流问题。     

求解一维扩散反应方程的隐式高精度紧致差分格式

提出了一维扩散反应方程的一种隐式高精度紧致差分格式,空间二阶导数采用四阶紧致差分格式进行离散,时间导数采用四阶向后欧拉公式进行离散,格式截断误差为O（τ⁴+h⁴）,即时间和空间都可以达到四阶精度,最后通过数值实验验证了本文方法的精确性和可靠性.     

控制体有限元方法中质量加权迎风格式的迁延修正

控制体有限元方法对于处理流动问题提供了一种较好的数值模拟方法，但对于大Ｒｅｙｎｏｌｄｓ数的回流流动问题，该方法在采用各种迎风格式时仍存在着严重的假扩散．为此，对二维质量加权迎风格式，给出了一种由插值单元及其一个与插值点最靠近的相邻单元所确定的多项式插值格式迁延修正的形式．通过对Ｒｅ＝４００、１０００、５０００三种情形的二维空穴流问题的计算，表明该格式可有效地减少控制体有限元方法对回流问题数值模拟的假扩散．     

二维粘弹介质五点八阶超紧致有限差分声波方程数值模拟

本文首次将五点八阶超紧致有限差分格式(CCD8)用于粘弹介质声波方程的数值模拟中。此文根据泰勒级数展开粘滞声波方程,建立了位移场时间二阶离散格式,将CCD8用于对位移场空间导数的求取。然后,对CCD8格式进行稳定性研究,频散的压制效果对比分析,及截断误差对比。将CCD8格式运用于均匀介质模型以及水平层状介质模型中进行数值模拟,最后运用于对Marmousi模型的数值模拟中。研究结果表明：(1)CCD8与普通八阶紧致差分相比,具有更小的截断误差、更高的模拟精度以及低数值频散的优势;(2)CCD8具有较高稳定性;(3)采用完全匹配层(PML)对人工边界进行处理后对均匀介质、含有层状介质以及Marmousi模型进行粘弹声波方程的数值模拟,发现模拟效果不错,从而验证了方法的实用性和有效。     

Numerical Simulation of Planar 4:1 Contraction Flow of a Viscoelastic Fluid Using a Higher-order Upwind Finite Volume Method

New form of the constitutive equations for the Oldroyd-B model, which have physical meaning, is developed to facilitate theoretical analysis. The new equations are used to simulate planar 4∶1 contraction flow of a Maxwell fluid using a third-order upwind finite volume method. The numerical results compare well with the theoretical solutions and the results of other references to show the effectiveness of the numerical method. Numerical experiments suggest that the present method not only converges fairly rapidly, but can also generate a highly resolved approximation to an Oldroyd-B fluid flow at a high Weissenberg number.     

Euler方程有限差分方法数值模拟

研究了定义在二维水槽上带非线性自由面边界条件的Euler方程的数值解.通过合适的σ坐标变换对不规则的水槽液体区域变换为一个规则的正方形区域,建立流场变量的差分耦合迭代的算法,运用交错网格求解了无粘不可压缩的Euler方程的数值解,数值结果表明,与之前结果和解析解比较数值解较好,对水平激励和垂直激励下非线性的效果和波拍的现象非常明显.     

三维水动力学模型高精度差分格式和解法研究

采用高精度差分格式对三维水动力学问题的时均化Navier-Stokes方程进行了数值模拟，进而采用广义共轭剩余法（GCR方法）求解压力泊松方程，并采用显式三级二阶Runge-Kutts格式模拟了时间步进过程。傅立叶分析表明，文中所采用的三阶迎风紧致差分格式具有较高的精度，数值实验进一步验证了上述数值模型的准确性和有效性。     

三层隐式差分格式在一维常系数扩散方程中的应用

讨论了具有二阶精度且无条件稳定的三层隐式差分格式在一维常系数扩散方程中的应用,通过对算例使用三种不同的方法分别计算第2时间层的近似值.对近似值及误差进行比较,得出向后差分格式是计算第2时间层较好的方法.     

用速度-涡量方法解具有表面抽吸的圆柱绕流问题

用有限差分法计算了不可压缩粘性流体绕具有表面抽吸圆柱的流动，控制方程采用涡量－速度形式的N－S方程，并引入对数有坐标变换，以使近壁处的网格加密，其中涡量输运方程采用二步R－K方法求解，速度泊松方程采用LSOR方法求解，并对速度进行了无散修正，计算表明合理选择抽吸的位置与强度，可有效减少阻力和升力，得到较好的流动控制效果。