用守恒元和解元法数值模拟二维溃坝洪水波

本文采用守恒元和解元方法(简称CE/SE法)来求解二维浅水方程,该方法从守恒积分型方程出发,严格遵守物理意义上的守恒律,通过设立守恒元和解元,从而提出一种新的数值模拟方法。用该方法计算得到的全溃的溃坝洪水波的数值解和理论解与文献中的实验值符合的较好。在此基础上,进一步对二维局部溃坝洪水波的演进过程进行模拟,结果表明,CE/SE法精度较高,稳定性较好,其溃坝洪水波的间断具有很高的分辨率,能很好的模拟溃坝洪水波的运动过程。

二维翼型绕流的CE/SE方法

采用时空守恒元和解元CE/SE(space-time Conservation Element and Solution Element method)法,求解二维Euler方程,开展了翼型绕流的无粘数值模拟研究.用非敏感克朗数计算格式消除克朗数过小引起的数值耗散对解的污染,结合当地时间步长法,解决网格不均匀引起的当地克朗数变化跨度大的问题.对NACA0012翼型的无激波流场进行了二维数值模拟,并与AGARD算例做了对比.结果表明:CE/SE方法的计算结果与AGARD结果吻合得很好,为该数值计算方法对翼型绕流数值模拟的进一步应用奠定了基础.

CE/SE法模拟爆震波点火及传播过程

运用时空守恒元和解元(CE/SE)法模拟一维爆震波的形成及传播过程,分析产生爆震波的不同机理和DDT(爆燃-爆震转捩)过程的影响因素.采取激波点火和热点点火两种点火方式,分别对应SDT(激波-爆震转捩)和DDT过程,热点点火初始场采用线性温度分布和阶梯性温度分布,整个管道均为常压.用隐式梯形方法实现刚性源项积分.结果表明,SDT和DDT对应于不同的爆震波机理,并且DDT受点火区不同参数的影响,为研究爆震波点火和DDT过程提供了理论依据.

用CE/SE法对弯曲与分叉河道的溃坝洪水波的数值研究

对作者的二维溃坝洪水波的数学模型进一步推广,得到了一般形式的基于任意四边形网格的时空守恒元和解元方法(简称CE/SE法)的新的格式。CE/SE法从守恒积分型浅水方程出发,设立守恒元和解元,严格保证其物理意义上的守恒律,并且构造思想简单,格式通用性好。首先采用CE/SE法计算等宽矩形河道的溃坝洪水波,并与Stoker解析解进行比较,在此基础上,数值模拟了180度强弯曲河道、45度三支分叉河道的二维溃坝洪水波的演进过程,揭示了溃坝洪水波在弯曲河道中内外两岸速度与水位的变化,在分叉河道中自动进行流量与动量的再分配,在分叉点处形成旋涡,水位变化剧烈等复杂的运动特征,算例结果表明基于任意四边形网格的CE/SE法精度高,稳定性好,该格式对各种不规则几何区域内的溃坝问题具有较强的适应性,对溃坝洪水波的间断具有较高的分辨率。

有、无障碍物爆震管起爆过程的CE/SE模拟

用时空守恒元和解元(CE/SE)法模拟有无障碍物情况下爆震管二维爆震波的形成和传播过程,分析爆震波的不同机理,研究在较低能量下,点火区温度、爆震管管径以及障碍物对爆燃-爆震转捩过程(DDT)的影响.控制方程采取无量纲形式,用拓展的隐式梯形法实现刚性源项积分.研究表明,点火能量与点火压力相关,在较低的点火能量下,可以通过提高点火温度、减小爆震管径及加入障碍物加强爆震过程,减小爆燃-爆震转捩过程DDT时间和距离.为研究爆震波点火和DDT过程提供了理论依据.