《中科院计算流体力学最新讲义CFD117讲有限体积法1.ppt》由会员分享，可在线阅读，更多相关《中科院计算流体力学最新讲义CFD117讲有限体积法1.ppt（25页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、中科院计算流体力学最新讲义CFD117讲有限体积法1 Still waters run deep.流静水深流静水深,人静心深人静心深 Where there is life,there is hope。有生命必有希望。有生命必有希望Copyright by Li Xinliang2知识回顾：知识回顾：Roe 格式格式平均斜率平均斜率线性化，以平均增长率代替瞬时增长率j,j+1区间内连续，且连续，且可通过相似变换对角化可通过相似变换对角化 应当具有的性质应当具有的性质常系数方程的Riemann解Copyright by Li Xinliang3知识回顾知识回顾2：LU-SGSStep 1:求解S
2、tep 2:求解Copyright by Li Xinliang4 7.1 结构网格有限体积法结构网格有限体积法有限体积法主要优势：有限体积法主要优势：处理复杂网格处理复杂网格差分法处理复杂外形差分法处理复杂外形 坐标变换坐标变换坐标变换函数坐标变换函数必须足够光滑必须足够光滑 否则损失精度否则损失精度实际问题：实际问题：外形复杂，外形复杂，光滑的结构网格生成困难光滑的结构网格生成困难差分法差分法有限体积法有限体积法优点优点简单、计算量小、易简单、计算量小、易于提高精度于提高精度本身包含几何信息，本身包含几何信息，易处理复杂网格易处理复杂网格不足不足差分离散与几何解耦，差分离散与几何解耦，难以
3、处理复杂网格难以处理复杂网格复杂、不易提高精度复杂、不易提高精度Copyright by Li Xinliang51.基本概念基本概念1）控制体控制体节点（中心）型控制体与网格型控制体节点（中心）型控制体与网格型控制体Copyright by Li Xinliang62）积分型控制方程积分型控制方程在控制体上积分物理含义：物理含义：控制体内总质量控制体内总质量/动量动量/能量的增加能量的增加=穿过控制体边界流入穿过控制体边界流入的净质量的净质量/动量动量/能量能量无粘通量粘性通量Copyright by Li Xinliang7物理含义：物理含义：单位时间内，无粘流动流过单位时间内，无粘流动流
4、过垂直于垂直于n方向的单位面积的质量、动量和方向的单位面积的质量、动量和能量能量法向速度Copyright by Li Xinliang83）有限体积法中物理量的含义）有限体积法中物理量的含义4）残差残差残差残差=净通量净通量=右端项右端项含义：含义：控制体内的平均量控制体内的平均量（平均质量密度、平均动量密度、（平均质量密度、平均动量密度、平均能量密度）平均能量密度）控制体控制体几何中心处几何中心处的当地密度、的当地密度、动量密度、能量密度动量密度、能量密度二阶精度二阶精度近似近似Copyright by Li Xinliang92.无粘通量的计算无粘通量的计算常用方法常用方法（流过（流过A
5、B边的通量）：边的通量）：a.利用周围点的值，计算出利用周围点的值，计算出(I+1/2,J)点处的物理量；点处的物理量；b.利用该处的物理量，计算出流过利用该处的物理量，计算出流过AB边的流通量边的流通量方法方法1：中心型有限体积法中心型有限体积法人工粘性项人工粘性项Copyright by Li Xinliang10方法方法2：迎风型有限体积法迎风型有限体积法Step 1:利用（偏）左侧点及（偏）右侧点利用（偏）左侧点及（偏）右侧点的值，计算出的值，计算出I+1/2,J点的值点的值计算方法：计算方法：与差分法完全相同与差分法完全相同各种差分格式，均可直接使用各种差分格式，均可直接使用也称为也
6、称为“差分格式差分格式”该过程称为该过程称为“重构重构”（很多文献中称为（很多文献中称为“插值插值”）有限体积与有限差有限体积与有限差分共通之处分共通之处,可直可直接使用差分格式接使用差分格式Copyright by Li Xinliang11常见的差分格式：常见的差分格式：2阶NND格式minmod(a,b):a,b符号相反时取0，符号相同时取绝对值小的3阶迎风3阶MUSCL格式TVD,WENO,GVC,保单调格式保单调格式Copyright by Li Xinliang12重构方式：重构方式：原始变量、守恒变量及特征变量原始变量、守恒变量及特征变量以NND格式为例：守恒变量重构原始变量重构
7、特征变量重构先算出先算出UI+1/2(可用可用UI和和UI+1的算术平均或的算术平均或Roe平均），再平均），再利用该值算出利用该值算出SI+1/2 （称为数值流通量）（称为数值流通量）的含义的含义Copyright by Li Xinliang13重要概念澄清：重构与插值A.有限差分法：有限差分法：j+1/2切线切线j-1/2jj-1 注意：注意：与与 f 在在xj+1/2点的值含义不同！点的值含义不同！用周围几个点的值用周围几个点的值 计算计算 的过程称为的过程称为“重构重构”，不能理，不能理解为用解为用 来来插值插值记号记号 确实容易混淆，让人容易联想起确实容易混淆，让人容易联想起 。记
8、为。记为 更好些更好些否则，最高只能否则，最高只能达到达到2阶精度了！阶精度了！是控制体内的平均值是控制体内的平均值 （称为数值流通量）（称为数值流通量）的含义的含义Copyright by Li Xinliang14重要概念澄清：重构与插值B.有限体积法：有限体积法：j+1/2j-1/2确实为确实为f在在xj+1/2点的值点的值！通常做法：通常做法：1）用用 计算出计算出 2）u在xj+1/2点的值！关键：关键：是用是用 计算计算 （称为（称为重构重构），而不是用，而不是用 计算计算 （是标准的（是标准的插值插值）；否则最高也只能达到）；否则最高也只能达到2阶精度。阶精度。Copyright
9、 by Li Xinliang15Step 2:利用利用 ，计算出通量计算出通量 方法方法1：流通矢量分裂：流通矢量分裂（FVS）方法方法2：通量差分分裂：通量差分分裂(FDS)通过通过Riemann解，获得通量解，获得通量可利用近似可利用近似Riemann阶（阶（Roe,HLL,HLLC)Roe格式：格式：xyxy方法方法3：AUSM方法方法压力项单独处理，其余项类似压力项单独处理，其余项类似Van Leer 分裂分裂Copyright by Li Xinliang小技巧：小技巧：利用局部坐标系，计算通量时简化利用局部坐标系，计算通量时简化为为1维问题维问题AByxy“扩展的”一维问题x,y
10、坐标系下的Riemann问题（二维）x,y 坐标系下的Riemann问题（一维）x,y 坐标系下，切向速度坐标系下，切向速度v表现相当于表现相当于被动标量被动标量Copyright by Li Xinliang173.粘性通量的计算粘性通量的计算关键问题：关键问题：计算计算 对对 的导数同样计算的导数同样计算Copyright by Li Xinliang18方法方法1.利用利用Jacobian变换变换方法方法2.利用利用Green积分公式计算积分公式计算最终公式与方法最终公式与方法1相同相同Copyright by Li Xinliang19 7.2 边界条件及网格分块边界条件及网格分块1.
11、边界条件的处理方法边界条件的处理方法常用方法：常用方法：虚网格虚网格设置设置1或或2层虚网格（如果采用高层虚网格（如果采用高精度格式，需设置更多重虚网格）精度格式，需设置更多重虚网格）虚网格点上的值根据边界条件给定虚网格点上的值根据边界条件给定 绝热固壁边界条件绝热固壁边界条件超声速入口边界条件超声速入口边界条件超声速出口边界条件超声速出口边界条件Copyright by Li Xinliang20亚声速入口边界条件亚声速入口边界条件外部计算域内部边界特征分析：内点提供一个边界条件，外部提供其余边界条件调节作用，内点速度降低时，增加入口压力；内点速度增加时，降低入口压力；Copyright b
12、y Li Xinliang21亚声速出口边界条件亚声速出口边界条件外部计算域内部边界特征分析：特征分析：外部提供一个边界条件；外部提供一个边界条件；通常情况下，指定背压通常情况下，指定背压如果边界处压力比内点高，则边界如果边界处压力比内点高，则边界处密度增加，速度降低；处密度增加，速度降低；如果边界压力比内点低，则边界处如果边界压力比内点低，则边界处密度降低，速度升高；密度降低，速度升高；Copyright by Li Xinliang222.多块网格多块网格复杂外形情况下，单块结构网格很难实现，需要由多块网格实现复杂外形情况下，单块结构网格很难实现，需要由多块网格实现Copyright by
13、 Li Xinliang231）对接网格对接网格对接网格示意图对接网格示意图相邻的网格块共享交界网格点相邻的网格块共享交界网格点处理方法：处理方法：利用虚网格传递信息；利用虚网格传递信息；利用连接信息，将内点的物理量利用连接信息，将内点的物理量复制给对应的虚网格点；复制给对应的虚网格点；Block2 的值的值 Block 1虚网格的值虚网格的值Copyright by Li Xinliang242）重叠网格重叠网格重叠网格示意图重叠网格示意图处理方法：处理方法：使用虚网格交换信息；使用虚网格交换信息；利用插值，用网格内点利用插值，用网格内点的信息计算出另一套网格的信息计算出另一套网格虚网格上的
14、信息。虚网格上的信息。Copyright by Li Xinliang25作业作业 7.1 编制有限体积程序编制有限体积程序,并计算图示并计算图示钝楔绕流问题钝楔绕流问题r=1几何及流动参数：几何及流动参数：钝楔的半楔角为钝楔的半楔角为5。来。来流流Mach数数6，Reynold数数10000（以头半（以头半径及来流参数度量）。径及来流参数度量）。来流攻角为来流攻角为0，壁面，壁面温度为来流温度的温度为来流温度的4.2倍。倍。计算中，粘性系数可用Sutherland公式计算：本计算中，来流温度设定为要求：要求：编写程序，给出壁面压力、温度的二维分布云图，并给编写程序，给出壁面压力、温度的二维分布云图，并给出壁面压力及热流。出壁面压力及热流。网格可从“流体中文网”下载（tecplot格式）