边界条件对对流扩散方程数值稳定性的影响

本文利用数值计算方法对采用均分网格的一维线性无源的对流－扩散方程在各种边界条件下的稳定性进行了分析,燕求出了不同边界条件下一维问题的中心差分和ＱＵＩＣＫ格式的临界网格Ｐｅｃｌｅｔ数。指出按现有方法得出的临界网格Ｐｅｃｌｅｔ数是判别差分格式对流数值稳定性的最苛刻的要求。对中心差分和ＱＵＩＣＫ格式,除两点边值问题以外的其它边界条件下的稳定性范围均不小于或远远大于两点边值问题的稳定性范围。通过计算还得出了格式的数值稳定性主要取决于计算区域下游侧的边界条件类型而与计算区域上游侧的边界条件类型无关的结论。     

低速多组份混合工质掺混流场数值模拟的加速收敛

本文对基于时间推进的高分辨率迎风数值格式采用一种新的加速收敛技巧,对低速掺混流动问题进行数值模拟, 可以有效地加快收敛速度,而且每步计算时间并未增加许多,因此能够有效地减少计算费用。同时改进后的数值格式能够 在处理高速流动问题的时候恢复成为原来的格式,从而使得改进后的数值格式对高速和低速流动问题均能保持良好的计 算效率。     

高分辨率激波捕捉格式及其CFD应用

在曲线坐标下将高分辨率AUSMPW激波捕捉格式扩展到三维问题。为了提高计算精度,采用三阶MUSCL格式,将AUSMPW格式与隐械时间推进法LU－SGS方法结合,应用于可压缩无粘和有粘流动的数散模拟,计算结果与文献中计算结果和实验数据相符,并将该格式与CD、TVD、FVS和FD格式进行了比较。     

抗爆容器内爆炸流场数值模拟

 采用计算流体动力学中的二阶精度TVD差分格式和特殊算子分裂法,按轴对称问题,对半球顶圆柱筒密闭式抗爆容器内部爆炸流场进行了数值模拟。计算得到的容器壁面载荷分布与实验结果基本一致。不同时刻爆炸流场压力分布图像清晰地描述了容器壁面的冲击波加载过程。     

一种结构动力响应分析的快速高阶算法

为了提高基于高阶格式的结构动力响应微分求积分析方法的计算效率,发展了一种求解动力方程的快速算法．利用微分求积原理将结构动力方程转化为标准Sylvester方程的形式,通过对系数矩阵进行矩阵分解,进而将动力响应Sylvester方程化为一系列标准线性方程组,采用相关成熟算法求解这些线性方程组后即可获得结构动力时程响应的全部解答．结构动力响应微分求积分析方法为高阶数值方法,一步计算可以获得多个时点处的动力响应．基于本文快速算法,不必直接对矩阵方程进行求解．数值算例表明,本文快速算法能够准确地计算出结构动力响应,具有数值精度高、收敛性好的优点．     

基于子结构解耦的连接特性识别新方法

连接部件动态特性的准确辨识对预测装配式机械结构的动力学行为有重要意义. 针对传统基于子结构解耦的连接结构动力学特性识别方法难以直接利用可测量数据进行辨识及辨识结果受噪声影响显著等问题, 本文提出了一种新方法. 首先, 提取子结构解耦基本方程在测试自由度上的分量, 并经矩阵变换得到显含连接动刚度矩阵的形式, 而后由真实连接动刚度矩阵分解为已知的初始矩阵与待求的增量矩阵, 推导了具有收敛性质的增量型方程以增强界面自由度较多时辨识的数值稳定性, 并采用多项式拟合动刚度将其转化为了拟合系数的频域估计方程, 按给定准则选取合适的频率点联立方程后, 得到了只需装配体测试自由度上的频响函数来辨识连接特性的迭代公式. 最后, 以若干算例说明了算法的具体流程. 对10自由度弹簧?质量块系统进行了数值仿真, 验证了所提方法的正确性及抗噪性; 对包含一处胶接连接的T形梁结构和包含两处螺栓连接的L形梁结构进行了试验, 所辨识连接结构与残余结构重组的装配体有限元模型计算的频响函数与测量值在较宽频带内吻合较好, 表明了该方法能有效识别实际装配体结构中的连接特性.      

一种基于块雅可比迭代的高阶FR格式隐式方法

最近, 基于非结构网格的高阶通量重构格式(flux reconstruction, FR)因其构造简单且通用性强而受到越来越多人的关注. 但将FR格式应用于大规模复杂流动的模拟时仍面临计算开销大、求解时间长等问题. 因此, 亟需发展与之相适应的高效隐式求解方法和并行计算技术. 本文提出了一种基于块Jacobi迭代的高阶FR格式求解定常二维欧拉方程的单GPU隐式时间推进方法. 由于直接求解FR格式空间和隐式时间离散后的全局线性方程组效率低下并且内存占用很大. 而通过块雅可比迭代的方式, 能够改变全局线性方程组左端矩阵的特征, 克服影响求解并行性的相邻单元依赖问题, 使得只需要存储和计算对角块矩阵. 最终将求解全局线性方程组转化为求解一系列局部单元线性方程组, 进而又可利用LU分解法在GPU上并行求解这些小型局部线性方程组. 通过二维无黏Bump流动和NACA0012无黏绕流两个数值实验表明, 该隐式方法计算收敛所用的迭代步数和计算时间均远小于使用多重网格加速的显式Runge-Kutta格式, 且在计算效率方面至少有一个量级的提升.      

基于五阶WENO格式的燃气在药床中流动过程二维两相流研究

为研究内弹道初始阶段中心点火管燃气在膛内药床中的流动特性和传播规律,设计了可视化点传火实验平台,并进行了膛内假药床的点传火实验。基于加权本质无震荡(weighted essentially non-oscillatory, WENO)格式,构造了膛内轴对称二维内弹道两相流模型,对膛内燃气在假药床中的流动过程进行数值模拟。计算结果与可视化实验结果符合较好,全局压力平均误差为5.35%。表明数值计算准确地描述了燃气流动特性,完整地呈现了点火管燃气在假药床中的发展过程。在点火初始阶段,膛内压力径向效应明显,气相沿径向传播较快,药床药粒基本不会发生运动;随着燃气逐渐在膛内传播,膛内压力呈现径向一致、轴向梯度分布的特征,在压力梯度作用下,气相轴向速度开始占据主导,径向速度在膛底和中部区域减小为零,而固相速度随气相速度变化而变化;气相在到达弹底前,由于固相颗粒的壅塞,会提前出现速度反向波动现象。     

一个改进的三阶WENO-Z型格式

在WENO-Z型格式框架下,基于高阶全局光滑因子,在非线性权建立过程中引入参数,通过收敛性分析确定参数取值范围,兼顾精确性与不振荡性,得到参数最佳取值.最终得到一个低耗散、高分辨率的三阶WENO差分格式,该格式在函数一阶极值点处仍保持预期三阶精度.最后通过精确解算例验证了格式在各种类型极值点处精度恢复情况,并通过一、二维Euler方程组经典算例测试了格式的低耗散、高分辨特性.结果表明,该文格式是一个性能优良的激波捕捉格式.     

一类奇异摄动对流扩散方程组的自适应网格方法的收敛性分析

基于经典的迎风有限差分方法,本文讨论一类奇异摄动对流扩散方程组的自适应网格算法.首先,利用Taylor级数展开,给出离散格式的局部截断误差估计.然后,利用等分布原理和极大值原理,证明基于弧长控制函数的自适应网格算法是一阶收敛的.最后的数值实验验证了本文的理论结果.