计算流体力学第一章绪论

《计算流体力学第一章绪论》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、第一章绪论§1.3模型方程以及在计算流体力学中的应用§1.4计算流体力学的控制方程§1.2流体力学微分方程的数学性质§1.1计算流体力学简介9/21/20211§1.1计算流体力学简介（一）什么是计算流体力学1.以电子计算机作为模拟手段，运用一定的计算技术寻求流体力学各种复杂问题的离散化数值解。★数值解而不是解析解★计算理论和计算技术起关键作用★与计算机的发展紧密相关9/21/202122．计算流体力学、理论流体力学、实验流体力学是流体力学研究工作的三种主要手段――既互相独立又相辅相成▲理论分析具有普遍性――各种影响因素清晰可见、为实验和计算研究提供依据▲实验研究

2、仍是研究工作的基石，数值研究的许多方面都密切依赖于实验研究：实验提供数据；计算结果需由实验验证；观察实验现象分析实验数据以建立计算模型等等▲数值模拟是特殊意义下的实验，也称数值实验9/21/20213(二）计算流体力学研究工作的方向◆与计算理论、计算技术的发展相关联的研究方向：（基本方程与计算物理，计算力学发展、图形学、网格技术等）★基本方程性质、简化、改型等★座标变换理论▲网格技术▲分区处理技术、并行处理与计算技术9/21/20214◆与离散数学的理论研究相关连的研究方向；▲离散化与精度分析理论▲稳定性分析理论，特别是边界条件数值处理的稳定性分析理论▲守恒性及其

3、数值分析中的守恒性理论▲数值分析的收敛性理论9/21/20215◆解决众多相关学科的的科研工作和工程实际提出的与流体力学问题有关的各类复杂的问题◆在一些相关学科的边缘上寻求新的发展点◇机械、航天航空、气象、海洋、船舶、车辆、石油、化工、环境（包括气动噪音控制）、建筑、土木、水利、核能等等▲与计算机技术的发展相关的一些领域▲交通流、信息流、生物流等的数值模拟9/21/20216WindaroundCoolingTowers9/21/20217（三）计算流体力学研究工作的优势、存在的问题和困难1．优势：☆“数值实验”比“物理实验”具有更大的自由度和灵活性，例如“自由”

4、地选取各种参数等☆“数值实验”可以进行“物理实验”不可能或很难进行的实验；例如：天体内部地温度场数值模拟，可控热核反应地数值模拟☆“数值实验”的经济效益极为显著，而且将越来越显著9/21/202189/21/202199/21/2021102．问题与不足△流动机理不明的问题，数值工作无法进行△数值工作自身仍然有许多理论问题有待解决△离散化不仅引起定量的误差，同时也会引起定性的误差，所以数值工作仍然离不开实验的验证9/21/202111§1.2流体力学微分方程的数学性质△当微分方程转化为差分方程并用数值方法求解时，不同类型的微分方程，其数值处理方法各异，其中包括提法

5、的适定性、物理解的性质、差分格式的适用性等△在一些特殊的问题中，甚至通过差分格式的特技巧来改变方程的数学性质9/21/202112（一）一阶拟线性微分方程组的分类对于一阶拟线性微分方程组的向量形式：其中：U为n阶向量，A为n阶矩阵若：A的特征值为:9/21/202113⑴.当n个特征值全部为复数时，称方程在(t,xi)平面上为纯椭圆型；⑵.n个特征值全部为互不相等的实数时，称方程在(t,xi),平面上为纯双曲型；而当n个特征值全部为实数，但有部分为相等的实数时，称方程(t,xi)在平面上为双曲型；⑶.当n个特征值全部为零时，称方程在(t,xi)平面上为纯抛物型；⑷

6、．当n个特征值部分为复数、部分为实数时，称方程在(t,xi)平面上为双曲椭圆型；二阶拟线性方程组，可以通过降阶法进行类似的分析。(例见后)9/21/2021141.二维定常理想流体流动的Euler方程写成向量形式：（一）流体力学控制方程数学分类的举例：9/21/2021159/21/202116求矩阵C的特征值得：如果：9/21/2021172。二维非定常理想流体流动的Euler方程求C的特征值，结论与定常相同：得到在X-Y平面的方程性质；9/21/202118求D的特征值,得:为四个实根，即方程在x-t平面为双曲型；所以Euler方程可以在时间座标方向推进，而在

7、定常问题中能否推进计算，必须根据流动是否为超音速（M与1的关系）来定。9/21/2021193.定常不可压缩Navier–Stokes方程的数学分类降阶法,令：9/21/2021209/21/202121接下来的分析与一阶拟线性方程组的讨论相似，结论为定常N-S方程为椭圆型。9/21/2021224．非定常不可压缩Navier–Stokes方程5．定常可压缩Navier–Stokes方程6．非定常可压缩Navier–Stokes方程7．抛物化N–S方程利用边界层流动的概念，设x方向为主流方向，即考虑有：把流动方向的二阶偏导数略去。（注意与边界层方程的区别）结论是定

8、常N-S方