第9讲-有限体积法1.pptVIP

Copyright by Li Xinliang * 知识回顾： Riemann问题精确解 Riemann问题 问题描述： 初始时刻，物理量分布存在单个间断；间断两侧物理量为常数。 求解思路： 采用积分方程 单个间断，且间断两侧物理量为常数情况下： 积分方程转化为代数方程 代数方程： 质量、动量、能量守恒 计算出 ， 将 与这三个值进行比较，判断会产生的情况。具体见下图： Copyright by Li Xinliang * Riemann问题的具体计算步骤 （全流场） 1. 判断可能会出现的情况（五种情形之一） a. 定义函数 b. 进行判断 情况5 情况4 情况3 情况1 情况5 情况4 情况2 情况1 单调增函数，性质很好 Copyright by Li Xinliang * 2. 求解中心区的压力和速度 单未知数的代数方程，迭代求解（例如Newton法，F(p)性质好，求解不困难） 3. 确定中心区接触间断两侧的密度 以及左、右波传播的速度 a. 左波为激波的情况（情况1,3） b. 左波为稀疏波的情况 （情况2，4，5） 中心区 接触间断左侧的物理量 膨胀波的波头及波尾速度 激波的传播速度 对于情况（5），波尾速度为： 中心区为真空，音速 无定义，改由该式计算 Copyright by Li Xinliang * c. 右波为激波的情况（情况1，2） 中心区 接触间断右侧的物理量 b. 右波为稀疏波的情况 （情况2，4，5） 4. 计算稀疏波区域的值（如果有稀疏波的话） a. 左稀疏波 b. 右稀疏波 情况2，4 情况5： 注意： 教科书32页c的公式有误！ Copyright by Li Xinliang * 有限体积法 “扩展的”Riemann问题的计算方法 （中心线x=0处） 迎风型有限体积法，需要求解“扩展的”一维Riemann问题 x y 物理问题分析： 所有物理量均沿x方向一维分布，沿y方向均匀分布。 仅需计算t时刻x=0处各物理量的值 v和w跟随流体运动，相当于“被动标量” 穿过激波及稀疏波，切向速度不变 Copyright by Li Xinliang * 求解t时刻x=0处物理量的具体步骤 Step 1: 求解 得到中心区压力 Step 2: 计算中心区的速度 Step 3: 根据 及 判断会出现哪种情况（五种情况之一） Step 4: 根据具体情况（左、右波是激波还是膨胀波）计算出中心区接触间断两侧的密度 及 Step 5: 如果中心区x方向速度0, 则中心线（x=0）处的密度及切向速度为接触间断右侧的值，否则为接触间断左侧的值。 具体公式见本PPT33-34页， 本页中上标“L”和“R”分别对应原先的下标“1”和“2” x t v和w没有给计算过程带来任何麻烦 先无视v和w的存在，求解标准1维Riemann问题。再根据u的符号确定中心线的v和w Copyright by Li Xinliang * 作业9.1 求解 “坐标旋转的Sod激波管问题” x y 物理问题描述：如右图示，有一条与x轴夹角120°的直线，左右两侧充满同种介质的无粘完全气体。初始时刻，左右两侧的气体状态为： 左侧： ，右侧 试计算t=0.14时刻的流动分布。 计算要求： 1） 计算域 网格 2）初值设置如图所示； 3）空间离散采用有限体积法计算。采用线性重构（见上页）。时间推进采用3阶Runge-Kutta方法 4) 分别采用FVS方法及Roe-FDS （又称Roe- Riemann近似解，见本PPT 22页）两种不同的方法计算通量。 5） 绘制出t=0.14时刻密度、速度u,v及压力的二维分布。 6） 绘制出t=0.14时刻x’轴（垂直于初始间断面，见右图）上的密度、沿x’方向的速度及压力的一维分布，并同精确解进行比较。 x y * * Copyright by Li Xinliang Copyright by Li Xinliang 计算流体力学讲义 第九讲 有限体积法（1） 李新亮 lixl@ ；力学所主楼219； 知识点： * 讲义、课件上传至 (