学位论文 > 优秀研究生学位论文题录展示

基于拟压缩方法的二维低速流动数值模拟

作　者: 叶军科
导　师: 李栋
学　校: 西北工业大学
专　业: 流体力学
关键词: 数值模拟低雷诺数非定常流动拟压缩方法有限体积法双时间法
分类号: O302
类　型: 硕士论文
年　份: 2004年
下　载: 217次
引　用: 2次
阅　读: 论文下载

内容摘要

本文运用拟压缩方法，发展了一套适用于二维低速流动的定常和非定常数值计算方法。通过在不可压连续性方程中引入拟压缩项，使控制方程成为一个封闭且可以沿时间方向推进求解的双曲型偏微分方程组。拟压缩项的引入使不可压方程在满足速度散度为零的同时耦合了速度和压强的变化。控制方程采用中心格式有限体积法进行空间离散，对于定常流动，运用Runge-Kutta显式多步法进行时间推进求解，非定常流动采用隐式时间离散的“双时间法”(Dual-Time Stepping Method)进行推进求解。通过对二维圆柱、平板、翼型绕流的数值模拟，验证了本文方法的正确性和有效性。本文的主要工作如下： 1．参照可压缩流动计算的中心格式有限体积法完成了拟压缩方法二维定常不可压Euler方程的推导和流动计算程序的编写。推导了Euler方程通量矢量的Jacobian矩阵在直角坐标系和曲线坐标系下的特征值，并由此导出相应的稳定性条件。根据理想不可压流体满足Bernoulli方程的原理，在程序中添加了压强阻尼项，加快了解的收敛速度。利用Euler程序计算并得到了二维圆柱表面的压力分布，结果与理论解一致。 2．推导了引入拟压缩项的二维定常不可压N-S方程，编写了基于中心格式有限体积法的二维定常N-S方程计算程序，给出了二维平板层流N-S方程计算结果和Blasius解的比较，验证了本文方法的正确性。进行了低雷诺数下绕翼型流动的数值模拟，得到了正确的计算结果。 3．推导了拟压缩方法二维非定常不可压Euler／N-S方程，编写了二维非定常Euler／N-S方程计算程序。该程序基于“双时间法”，适用于翼型作任意运动的非定常低速流动计算。 4．研究了在二维定常不可压流动状态下雷诺数和翼型厚度等参数对翼型气动特性的影响。 5．对于非定常不可压流动，分析了运动参数对翼型气动特性的影响，为进一步研究绕三维扑翼的非定常流动打下基础。本文的计算结果与分析表明，本文所使用的方法是正确的，且定常计算结果可靠。运用本文方法分析了各参数对翼型气动特性的影响，进一步结合叶素理论还可以估算扑翼的气动特性，进而指导微旋翼或扑翼微型飞行器的初期设计。

全文目录

摘要  3-4
Abstract  4-6
目录  6-8
符号表  8-10
前言  10-13
第一章二维不可压Euler方程的拟压缩方法求解  13-29
  1.1 控制方程  13-14
  1.2 方程的无量纲化  14
  1.3 空间离散  14-18
    1.3.1 通量项的计算  14-17
    1.3.2 人工耗散项的计算  17-18
  1.4 时间积分  18
  1.5 稳定性条件  18-19
  1.6 边界条件  19-21
    1.6.1 物面边界条件  20
    1.6.2 割缝条件  20-21
    1.6.3 远场边界条件  21
  1.7 加速收敛技术  21-23
  1.8 算例分析  23
  1.9 结论  23-29
第二章二维定常不可压N-S方程的拟压缩方法求解  29-58
  2.1 控制方程  29-30
  2.2 方程的无量纲化  30
  2.3 空间离散  30-31
  2.4 粘性应力的计算  31-32
  2.5 时间步长  32-33
  2.6 物面边界条件  33
  2.7 算例分析  33-35
  2.8 结论  35-58
第三章二维非定常不可压Euler／N-S方程的拟压缩方法求解  58-71
  3.1 控制方程  58-59
  3.2 空间离散和时间推进  59-61
  3.3 边界条件  61
  3.4 几何守恒定理  61-62
  3.5 算例分析  62-65
  3.6 结论  65-71
第四章总结与展望  71-72
  4.1 研究结论  71
  4.2 下一步工作重点和目标  71-72
附录A 拟压缩方程特征值的推导  72-77
  1 直角坐标系下二维Euler方程的特征值  72-73
  2 曲线坐标系下二维Euler方程的特征值  73-77
附录B 一维流动的拟压缩不变量  77-80
参考文献  80-82
致谢  82-83

基于拟压缩方法的二维低速流动数值模拟

内容摘要

全文目录

相似论文