计算流体力学的大作业
计算流体力学大作业
（2007.05）
一、一维标量问题 标量方程0=??+??x
u a t u ， =0.2。 a 初值为：
???????>?≤
.0 0.17
.00.5 )5.0(5.70.50.2 )10)2.0sin((2.0 0.0)(x x x x x x x u π精确解为。计算时间取为t ＝1s ，计算区域为，网格为均匀网格，网格数100～200。
),(),(0t at x u t x u ?=]1 ,0[（1） 用一阶迎风格式进行差分离散；
（2） 采用Lax －Wendroff 格式进行差分离散；
（3） 采用Harten TVD 格式进行空间离散，用具有TVD 性质二阶
Runge-Kutta 方法进行时间离散。
（4） 比较、分析以上几种差分格式离散数值解的耗散与频散特性。
00.250.50.751
x -1-0.5
0.5
1
1.5
U
A
c c u r a c y E
N O -2E
N O -3M U S C L
图1 t=1s 函数图象
二、一维Lax 问题
考察经典的Riemann 问题，即初始间断的分解问题。
一维守恒律的Euler 方程是0=??+??x
F t U ，方程中相对应的列矢量为： ??????????=E u U ρρρ ??
????????++=u p E p u u F )(2ρρρ()212u p E +?=γρ Lax 问题的初始间断分布
U 0L ＝(ρ=0.445, u=0.698, p=3.528)，
U 0R ＝(ρ=0.5, u=0.0, p=0.571)，
计算条件：给出t=0.15s 流场的密度、速度、压强分布。
（1） 空间差分采用一阶矢通量分裂格式进行离散，时间采用一阶Euler 时
间前差。
（2） 采用MacCormack 两步格式进行离散；
（3） 采用NND 格式进行空间离散，用具有TVD 性质二阶Runge-Kutta 方
法进行时间离散。
（4） 采用Harten TVD 格式，时间方向采用一阶Euler 时间前差。 （5） 比较、分析以上几种差分格式离散数值解的耗散与频散特性；
（6） 比较采用50，100，200个均匀网格计算得到的数值解的分辨率、收
敛性等特性。
本题说明：
(1)解题分组：
? 学号为单数的同学采用(1),(4)方法解题;
? 学号为双数的同学采用(2),(3)方法解题
（2）变换矩阵：
1?Λ=R R A
c u c u u +=?==321,,λλλ
c
为音速。ργp c =2，4.1=γ
ρP E H +=
????????
??????????+????+???????=c uc H c uc H u c c c u c c u u c c c c R 222
122121211
2222
γγγ ()()???????
?????????????????+????+??????=?c u c
u c u c c u u c u u c u R 111211112111222221γγγγγγγ
(3) 精确解
大作业说明：
（1） 两题均要求完成。
（2） 编程语言可以采用C，Fortran，matlab。（3） 可通过电子邮件或纸质方式上交。（4） 学期结束前上交，否则考试扣30分。
2020-2021年中国科学院大学(中科院)计算数学考研招生情况、分数线、参考书目、经验指导
一、中国科学院数学与系统科学研究院简介 中国科学院数学与系统科学研究院由中科院数学研究所、应用数学研究所、系统科学研究所及计算数学与科学工程计算研究所四个研究所整合而成，此外还拥有科学与工程计算国家重点实验室、中科院管理决策与信息系统重点实验室、中科院系统控制重点实验室、中科院数学机械化重点实验室、华罗庚数学重点实验室、随机复杂结构与数据科学重点实验室，以及中科院晨兴数学中心和中科院预测科学研究中心等。2010年11月成立国家数学与交叉科学中心，旨在从国家层面搭建一个数学与其它学科交叉合作的高水平研究平台。数学与系统科学研究院拥有完整的学科布局，研究领域涵盖了数学与系统科学的主要研究方向。共有16个硕士点和13个博士点（二级学科），分布在经济学、数学、系统科学、统计学、计算机科学与技术、管理科学与工程六个一级学科中，可以在此范围内招收和培养硕士与博士研究生。在2006年全国学科评估中，我院数学学科的整体评估得分为本学科的最高分数。数学与系统科学研究院硕士招生类别为硕士研究生、硕博连读生和专业学位硕士研究生。2019年共计划招收122名。 二、中国科学院大学计算数学专业招生情况、考试科目
三、中国科学院大学计算数学专业分数线 2018年硕士研究生招生复试分数线 2017年硕士研究生招生复试分数线 四、中国科学院大学计算数学专业考研参考书目 616数学分析 现行（公开发行）综合性大学（师范大学）数学系用数学分析教程。 801高等代数 [1] 北京大学编《高等代数》，高等教育出版社，1978年3月第1版，2003年7月第3版，2003年9月第2次印刷. [2] 复旦大学蒋尔雄等编《线性代数》，人民教育出版社，1988. [3] 张禾瑞，郝鈵新，《高等代数》，高等教育出版社， 1997. 五、中国科学院大学计算数学专业复试原则 在中国科学院数学与系统科学研究院招生工作小组领导下，按研究所成立招收硕士研究生复试小组，设组长1人、秘书1人。 复试总成绩按百分制计算，其中专业知识成绩占60％，英语听力及口语测试成绩占20％，综合素质成绩占20%。 在面试环节，每位考生有5分钟自述，考查内容主要包括专业知识、外语（口语）水平
流体力学大作业
《计算流体力学》课程大作业 作业内容：3-4人为小组完成数值模拟，在第8次课上每组进行成果展示，并在课程结束后每组上交一份纸质版报告。 数值模拟实现形式：自编程或者使用任意的开源、商业模型。 成果展示要求：口头讲述和幻灯片结合的方式，每组限时10分钟（8分钟讲述，2分钟提问和讨论）。 报告要求：按照期刊论文的思路和格式进行撰写（包括但不限于如下内容：摘要、绪论\引言、数值模型简介、数值结果分析\讨论、结论、参考文献）。 （以下题目二选一） 题目一：固定单方柱扰流问题 根据文章《Interactions of tandem square cylinders at low Reynolds numbers》中的实验进行数值模拟，完成但不局限于如下工作： （1）根据Fig. 2 中的雷诺数和方柱排列形式，进行相同雷诺数不同间距比情况下的方柱绕流数值模拟，并做出流线图和Fig.2中的结果对比。 （2）根据Fig. 3 中的雷诺数和方柱排列形式，进行相同雷诺数后柱不同转角情况下的方柱绕流数值模拟，并做出流线图和Fig.3中的结果对比。 （3）根据Fig. 12, 13 中的雷诺数和方柱间距比的设置进行数值模拟，作出频率、斯特劳哈尔数、阻力系数随雷诺数变化的折线并与图中对应的折线画在同一坐标系下比较。 （中共有4条折线，对应4种不同的方柱排列形式下的物理参数随雷诺数变化的规律，仅需选取单柱模型和其中一种双柱模型进行数值模拟，共计16个工况）。 题目二：溃坝问题 根据文章《Experimental investigation of dynamic pressure loads during dam break》中的实验进行数值模拟，完成但不局限于如下工作： （1）分别完成二维、三维的溃坝的数值建模，讨论二维、三维模型的区别。 （2）分别将二维、三维溃坝的数值模拟结果和Fig. 7,10中各时刻的自由面形态进行对比，并分别观测溃坝前端水舌的位置随时间的变化，其结果和Fig. 12 种的各试验结果放在同一坐标系下进行对比。 （3）根据实验设置数值观测点，分别观测与实验测点相对应的数值观测点上的水体高度、压力随时间的变化曲线，并和Fig.16, 18,21,30,31,32,33,35中的实验结果进行对比。
计算流体力学课后题作业
课后习题 第一章 1.计算流体动力学的基本任务是什么 计算流体动力学是通过计算机数值计算和图像显示，对包含有流体流动和热传导等相关物理现象的系统所做的分析。 2.什么叫控制方程？常用的控制方程有哪几个？各用在什么场合？ 流体流动要受物理守恒定律的支配，基本的守恒定律包括：质量守恒定律、动量守恒定律、能量守恒定律。如果流动包含有不同组分的混合或相互作用，系统还要遵守组分守恒定律。如果流动处于湍流状态，系统还要遵守附加的湍流输运方程。控制方程是这些守恒定律的数学描述。 常用的控制方程有质量守恒方程、动量守恒方程、能量守恒方程、组分质量守恒方程。质量守恒方程和动量守恒方程任何流动问题都必须满足，能量守恒定律是包含有热交换的流动系统必须满足的基本定律。组分质量守恒方程，在一个特定的系统中，可能存在质的交换，或者存在多种化学组分，每种组分都需要遵守组分质量守恒定律。 4.研究控制方程通用形式的意义何在？请分析控制方程通用形式中各项的意义。 建立控制方程通用形式是为了便于对各控制方程进行分析，并用同一程序对各控制方程进行求解。
各项依次为瞬态项、对流项、扩散项、源项。 6.CFD商用软件与用户自行设计的CFD程序相比，各有何优势？常用的商用CFD软件有哪些？特点如何？ 由于CFD的复杂性及计算机软硬件条件的多样性，用户各自的应用程序往往缺乏通用性。 CFD商用软件的特点是 功能比较全面、适用性强。 具有比较易用的前后处理系统和其他CAD及CFD软件的接口能力，便于用户快速完成造型、网格划分等工作。 具有比较完备的容错机制和操作界面，稳定性高。 可在多种计算机、多种操作系统，包括并行环境下运行。 常用的商用CFD软件有PHOENICS、CFX、SRAR-CD、FIDAP、FLUENT。PHOENICS除了通用CFD软件应该拥有的功能外，PHOENICS软件有自己独特的功能：开放性、CAD接口、运动物体功能、多种模型选择、双重算法选择、多模块选择。 CFX除了可以使用有限体积法外，还采用基于有限元的有限体积法。用于模拟流体流动、传热、多相流、化学反应、燃烧问题。其优势在于处理流动物理现象简单而几何形状复杂的问题。 SRAR-CD基于有限体积法，适用于不可压流体和可压流的计算、热力学的计算及非牛顿流的计算。它具有前处理器、求解器、后处理器三大模块，以良好的可视化用户界面把建模、求解及后处理与全部的物理模型和算法结合在一个软件包中。
流体力学实验报告
流体力学 实验指导书与报告 静力学实验 雷诺实验 中国矿业大学能源与动力实验中心
学生实验守则 一、学生进入实验室必须遵守实验室规章制度，遵守课堂纪律，衣着整洁，保持安静，不得迟到早退，严禁喧哗、吸烟、吃零食和随地吐痰。如有违犯，指导教师有权停止基实验。 二、实验课前，要认真阅读教材，作好实验预习，根据不同科目要求写出预习报告，明确实验目的、要求和注意事项。 三、实验课上必须专心听讲，服从指导教师的安排和指导，遵守操作规程，认真操作，正确读数，不得草率敷衍，拼凑数据。 四、预习报告和实验报告必须独自完成，不得互相抄袭。 五、因故缺课的学生，可向指导教师申请一次补做机会，不补做的，该试验以零分计算，作为总成绩的一部分，累计三次者，该课实验以不及格论处，不能参加该门课程的考试。 六、在使用大型精密仪器设备前，必须接受技术培训，经考核合格后方可使用，使用中要严格遵守操作规程，并详细填写使用记录。 七、爱护仪器设备，不准动用与本实验无关的仪器设备。要节约水、电、试剂药品、元器件、材料等。如发生仪器、设备损坏要及时向指导教师报告，属责任事故的，应按有关文件规定赔偿。 八、注意实验安全，遵守安全规定，防止人身和仪器设备事故发生。一旦发生事故，要立即向指导教师报告，采取正确的应急措施，防止事故扩大，保护人身安全和财产安全。重大事故要同时保护好现场，迅速向有关部门报告，事故后尽快写出书面报告交上级有关部门，不得隐瞒事实真相。 九、试验完毕要做好整理工作，将试剂、药品、工具、材料及公用仪器等放回原处。洗刷器皿，清扫试验场地，切断电源、气源、水源，经指导教师检查合格后方可离开。 十、各类实验室可根据自身特点，制定出切实可行的实验守则，报经系(院)主管领导同意后执行，并送实验室管理科备案。 1984年5月制定 2014年4月再修订 中国矿业大学能源与动力实验中心
流体力学 大作业
一.选择题 1.牛顿内摩擦定律适用于（）。 A．任何流体B．牛顿流体C．非牛顿流体 2.液体不具有的性质是（）。 A．易流动性B．压缩性C．抗拉性D．粘滞性 3连续介质假定认为流体（）连续。 A．在宏观上B．在微观上C．分子间D．原子间 4.在国际单位制中流体力学基本量纲不包括（）。 A．时间B．质量C．长度D．力． 5.在静水中取一六面体，作用在该六面体上的力有（） A．切向力、正压力B．正压力C．正压力、重力D．正压力、切向力、重力 6．下述哪些力属于质量力( ) A．惯性力B．粘性力C．弹性力D．表面张力E．重力 7．某点存在真空时，（）（） A．该点的绝对压强为正值B．该点的相对压强为正值c.该点的绝对压强为负值D．该点的相对压强为负值 8.流体静压强的（）。 A．方向与受压面有关B．大小与受压面积有关B．大小与受压面方位无关 9.流体静压强的全微分式为（）。 A．B．C． 10.压强单位为时，采用了哪种表示法（）。 A．应力单位B．大气压倍数C．液柱高度 11.密封容器内液面压强小于大气压强，其任一点的测压管液面（）。A．高于容器内液面B．低于容器内液面C．等于容器内液面 12.流体运动的连续性方程是根据（）原理导出的。 A.动量守恒 B. 质量守恒 C.能量守恒 D. 力的平衡 13. 流线和迹线重合的条件为（）。
A.恒定流 B.非恒定流 C.非恒定均匀流 14.总流伯努利方程适用于（）。 A.恒定流 B.非恒定流 C.可压缩流体 15. 总水头线与测压管水头线的基本规律是：（）、（） A.总水头线总是沿程下降的。 B.总水头线总是在测压管水头线的上方。 C.测压管水头线沿程可升可降。 D.测压管水头线总是沿程下降的。 16 管道中液体的雷诺数与（）无关。 A. 温度 B. 管径 C. 流速 D. 管长 17.. 某圆管直径d=30mm，其中液体平均流速为20cm/s。液体粘滞系数为0.0114cm3/s，则此管中液体流态为（）。 A. 层流 B. 层流向紊流过渡 C.紊流 18.等直径圆管中紊流的过流断面流速分布是（）A呈抛物线分布B. 呈对数线分布 C.呈椭圆曲线分布 D. 呈双曲线分布 19.等直径圆管中的层流，其过流断面平均流速是圆管中最大流速的（） A 1.0倍B.1/3倍C. 1/4倍D. 1/2倍 20.圆管中的层流的沿程损失与管中平均流速的（）成正比. A. 一次方 B. 二次方 C. 三次方 D. 四次方 21..圆管的水力半径是( ) A. d/2 B. d/3 C. d/4 D. d/5. 22谢才公式中谢才系数的单位是（）A. 无量纲B. C. D. . 23. 判断层流和紊流的临界雷诺数是（） A.上临界雷诺数 B.下临界雷诺数 C.上下临界雷诺数代数平均 D.上下临界雷诺数几何平均 24.. 对于管道无压流，当充满度分别为（）时，其流量和速度分别达到最大。 A. 0.5, 0.5 B. 0.95, 0.81 C. 0.81, 081 D. 1.0, 1.0 25.对于a, b, c三种水面线，下列哪些说法是错误（）（） A．所有a、c型曲线都是壅水曲线，即，水深沿程增大。B．所有b型曲线都是壅水曲线，即，水深沿程增大。C．所有a、c型曲线都是降水曲线，即，水深沿程减小。C．所有b型曲线都是降水曲线，即，水深沿程减
计算流体力学作业习题
2014级西安理工大学计算流体力学作业 1.写出通用方程，并说明其如何代表各类守恒定律。 由守恒型对流-扩散方程： ()()() div U div T grad S t φφρφρφφ?+=+? 其中φ为通用变量；T φ为广义扩散系数；S φ为广义原项。 若令1;1;0T S φφφ===时，则得到质量守恒方程（mass conservation equation ） ()()()() 0u v w t x y z ρρρρ????+++=???? 若令;i u φ=时，则得动量守恒方程（momentum conservation equation ） 以x 方向为例分析，设;u P u S S x φφ?==- ?，通用方程可化为： ()()()()(2)u uu vu wu P u divU t x y z x x x ρρρρλη???????+++=-++??????? z v u u w F y x y z z x ηηρ???????????? ??+++++?? ? ????????????????? 同理可证明y 、z 方向的动量守恒方程式 若令;;T p T T S S C φφλ φ===时，则得到能量守恒方程(energy conservation equation) ()()() ()h h div Uh div U div gradT S t ρρρλφ?+=-+++? ()()()T p h div Uh div gradT S t C ρλ ρ?+=+? 证毕 2.用控制体积法离散 0)(=+++s dx dT k dx d dx dT u dt dT ，要求对S 线性化，据你的理解，谈谈网格如何划分？交界面传热系数何如何计算？边界条件如何处理？ 根据守恒型对流-扩散方程： ()()()u T S t x x x ρφρ?φ ????' +=+????，对一维模型 进行分析，则有： 0)(=+++s dx dT k dx d dx dT u dt dT
计算流体力学大作业
1 提出问题 [问题描述] Sod 激波管问题是典型的一类Riemann 问题。如图所示，一管道左侧为高温高压气体，右侧为低温低压气体，中间用薄膜隔开。t=0 时刻，突然撤去薄膜，试分析其他的运动。 Sod 模型问题：在一维激波管的左侧初始分布为：0 ,1 ,1111===u p ρ，右侧分布为：0 ,1.0 ,125.0222===u p ρ，两种状态之间有一隔膜位于5.0=x 处。隔膜突然去掉，试给出在14.0=t 时刻Euler 方程的准确解，并给出在区间10≤≤x 这一时刻u p , ,ρ的分布图。 2 一维Euler 方程组 分析可知，一维激波管流体流动符合一维Euler 方程，具体方程如下： 矢量方程： 0U f t x ??+=?? (0.1) 分量方程： 连续性方程、动量方程和能量方程分别是： 2 22,,p u ρ
() ()()()2 000u t x u u p t x x u E p E t x ρρρρ???+ =?????????++=? ??????+?????+ =????? (0.2) 其中 22v u E c T ρ?? =+ ?? ? 对于完全气体，在量纲为一的形式下，状态方程为： ()2 p T Ma ργ∞ = (0.3) 在量纲为一的定义下，定容热容v c 为： () 21 1v c Ma γγ∞= - (0.4) 联立(1.2)，(1.3)，(1.4)消去温度T 和定容比热v c ，得到气体压力公式为： ()2112p E u γρ??=-- ??? (0.5) 上式中γ为气体常数，对于理想气体4.1=γ。 3 Euler 方程组的离散 3.1 Jacibian 矩阵特征值的分裂 Jacibian 矩阵A 的三个特征值分别是123;;u u c u c λλλ==+=-，依据如下算法将其分裂成正负特征值： () 12 222 k k k λλελ±±+= (0.6) 3.2 流通矢量的分裂 这里对流通矢量的分裂选用Steger-Warming 分裂法，分裂后的流通矢量为 ()()()()()()()12312322232121212122f u u c u c u u c u c w γλλλργλλλγλλγλ?? ? -++ ?=-+-++ ? ? ? -+-+++ ??? ＋＋＋＋＋＋＋ ＋＋ ＋＋ (0.7)
科学与工程计算国家重点实验室(中科院数学与系统科学研究所)
科学与工程计算国家重点实验室 简介 中国科学院科学与工程计算国家重点实验室（简称LSEC）是在已故著名数学家、中国计算数学的奠基人和开拓者冯康院士的倡导、并亲自筹备和组织下，由原中科院计算中心从事计算数学研究的部分课题组成的。实验室筹建于1990年，1993年10月经中科院验收后正式投入运行，1994年向国内外开放，1995年9月和 2005年3月两次通过国家验收。 实验室主要开展科学与工程计算中具有重要意义的基础理论研究，解决科学与工程领域中的重大计算问题，着重研究计算方法的构造、理论分析及实现。研究内容包括：动力系统与数值方法，研究各类保结构算法的理论、算法的构造和数值试验；有限元边界元方法，针对具有应用背景的椭圆边值问题及其它相关问题，提出适合于这些问题的有限元边界元新型高性能计算方法；非线性最优化，主要研究求解非线性规划的新算法以及算法的收敛性；计算流体力学,研究非定常不可压N-S方程和可压缩流的计算方法；并行计算方法和科学计算可视化；非均匀多孔介质中渗流问题的多尺度计算方法。 实验室主任是陈志明研究员。实验室学术委员会主任是中国工程院院士崔俊芝。 实验室建设以来在动力系统几何算法，非线性优化，有限元边界元，数理方程反问题，计算流体力学，并行算法，科学计算可视化等方面取得了大量的研究成果，十分突出的是关于哈密尔顿系统的辛几何算法的研究。其成果荣获“国家自然科学一等奖”。实验室在设备研制方面也取得了显著的成绩。 实验室现有科研人员19人，中科院院士2人（石钟慈、林群）,中国工程院院士1人（崔俊芝），其中研究员16人，此外，实验室还获得多项其它重要奖项，其中石钟慈院士在 2000年获“何梁何利科学与技术进奖”,林群院士获2001年获捷克科学院“数学科学成就荣誉奖”、2004年获“何梁何利科学与技术进奖”。实验室十分重视队伍建设和人才培养工作，尤其注重青年学术骨干的培养和引进。目前通过中科院“百人计划”已引进３位年轻的学科带头人，其中实验室主任陈志明研究员被国家科技部任命为９７３计划项目“高性能科学计算研究”首席科学家，一批优秀青年学术骨干脱颖而出，他们在各自的研究领域取得了可喜的成果，并因此获得了荣誉。例如，袁亚湘研究员曾获1995年首届“冯康科学计算奖”、1996年度“中国青年科学家奖”、“国家杰出青年科学基金”、1998年度“全国十大杰出青年”称号；2005年度“北京市科学技术一等奖”；张林波研究员曾获1995年度“中科院青年科学家二等奖”、1997年度“中科院优秀青年”奖、2000年度“国家科技进步奖二等奖”；白中治研究员获得1998年度“中科院自然科学三等奖”、1999年度“中科院青年科学家二等奖”、“中科院优秀青年”称号、2005年度“国家杰出青年科学基金”；许学军研究员获2000年度“钟家庆数学奖”；陈志明研究员获2000年度“国家杰出青年科学基金”、2001年度“第四届冯康科学计算奖”、2003年度“第七届中科院杰出青年”称号、2004年度“新世纪百千万人才工程国家级人选”、2005年度“海外青年学者合作研究基金”；周爱辉研究员获2004年度“国家杰出青年科学基金”。
计算流体力学课程大作业
《计算流体力学》课程大作业 ——基于涡量-流函数法的不可压缩方腔驱动流问题数值模拟 张伊哲 航博101 1、 引言和综述 2、 问题的提出，怎样使用涡量-流函数方法建立差分格式 3、 程序说明 4、 计算结果和讨论 5、 结论 1引言 虽然不可压缩流动的控制方程从形式上看更为简单，但实际上，目前不可压缩流动的数值方法远远不如可压缩流动的数值方法成熟。 考虑不可压缩流动的N-S 方程： 01()P t νρ??=? ? ??+??=-?+???? U U UU f U (1.1) 其中ν是运动粘性系数，认为是常数。将方程组写成无量纲的形式： 01()Re P t ??=?? ??+??=-?+????U U UU f U (1.2) 其中Re 是雷诺数。 从数学角度看，不可压缩流动的控制方程中不含有密度对时间的偏导数项，方程表现出椭圆-抛物组合型的特点；从物理意义上看，在不可压缩流动中，压力这一物理量的波动具有无穷大的传播速度，它瞬间传遍全场，以使不可压缩条件在任何时间、任何位置满足，这就是椭圆型方程的物理意义。这就造成不可压缩的N-S 方程不能使用比较成熟的发展型．．．偏微分方程的数值求解理论和方法。 如果将动量方程和连续性方程完全耦合求解，即使使用显示的离散格式，也将会得到一个刚性很强的、庞大的稀疏线性方程组，计算量巨大，更重要的问题是不易收敛。因此，实际应用中，通常都必须将连续方程和动量方程在一定程度上解耦。 目前，求解不可压缩流动的方法主要有涡量-流函数法，SIMPLE 法及其衍生的改进方法，有限元法，谱方法等，这些方法各有优缺点。其中涡量-流函数法是解决二维不可压缩流动的有效方法。作者本学期学习了研究生计算流体课程，为了熟悉计算流体的基本方法，选择使用涡量-流函数法计算不可压缩方腔驱动流问题，并且对于不同雷诺数下的解进行比较和分析，得出一些结论。 本文接下来的内容安排为：第2节提出不可压缩方腔驱动流问题，并分析该问题怎样使用涡量-流函数方法建立差分格式、选择边界条件。第3节介绍程序的结构。第4节对于不同雷诺数下的计算结果进行分析，并且与U.GHIA 等人【1】的经典结论进行对比，评述本
流体力学大作业
流体力学-大作业
————————————————————————————————作者: ————————————————————————————————日期：
一．选择题 1．牛顿内摩擦定律适用于()。 A.任何流体B．牛顿流体 C.非牛顿流体 2．液体不具有的性质是()。 A．易流动性Ｂ.压缩性Ｃ.抗拉性 D.粘滞性 3连续介质假定认为流体（)连续。 A．在宏观上 B.在微观上 C.分子间Ｄ．原子间 4.在国际单位制中流体力学基本量纲不包括（）。 A．时间 B.质量 C.长度D．力． 5.在静水中取一六面体，作用在该六面体上的力有（) Ａ．切向力、正压力Ｂ．正压力C．正压力、重力D．正压力、切向力、重力 ６. 下述哪些力属于质量力（) A．惯性力Ｂ.粘性力 C.弹性力Ｄ.表面张力Ｅ．重力 7.某点存在真空时,( ）(） A.该点的绝对压强为正值 B.该点的相对压强为正值c.该点的绝对压强为负值 D.该点的相对压强为负值 ８．流体静压强的( )。 A.方向与受压面有关Ｂ.大小与受压面积有关B．大小与受压面方位无关 9．流体静压强的全微分式为（)。 A. B. C． １0．压强单位为时，采用了哪种表示法（)。 A．应力单位Ｂ.大气压倍数C．液柱高度 11.密封容器内液面压强小于大气压强，其任一点的测压管液面( ）。 A.高于容器内液面B．低于容器内液面 C.等于容器内液面 12．流体运动的连续性方程是根据( ）原理导出的。 A.动量守恒 B. 质量守恒 C.能量守恒 D. 力的平衡 13．流线和迹线重合的条件为（)。
A．恒定流B.非恒定流C.非恒定均匀流 １４．总流伯努利方程适用于（）。 A.恒定流 B.非恒定流Ｃ.可压缩流体 15. 总水头线与测压管水头线的基本规律是：( ）、( ） A.总水头线总是沿程下降的。B．总水头线总是在测压管水头线的上方。 C.测压管水头线沿程可升可降。 D.测压管水头线总是沿程下降的。 １6 管道中液体的雷诺数与（)无关。 A.温度B．管径Ｃ. 流速D．管长 １7．. 某圆管直径d=3０mｍ，其中液体平均流速为20ｃｍ/s。液体粘滞系数为0.0１14cm3/s，则此管中液体流态为( )。 A. 层流 B. 层流向紊流过渡C．紊流 1８．等直径圆管中紊流的过流断面流速分布是（)A呈抛物线分布B．呈对数线分布 C.呈椭圆曲线分布D．呈双曲线分布19．等直径圆管中的层流,其过流断面平均流速是圆管中最大流速的（） A １．０倍 B.1/３倍Ｃ.1/４倍D. 1/2倍 20.圆管中的层流的沿程损失与管中平均流速的（）成正比. A. 一次方 B.二次方 C. 三次方Ｄ. 四次方 ２1..圆管的水力半径是（) A. d/2Ｂ.d/3 C. d/４D. d／5. 2２谢才公式中谢才系数的单位是（）A.无量纲B．Ｃ.D．. 23．判断层流和紊流的临界雷诺数是（） Ａ.上临界雷诺数 B.下临界雷诺数 C.上下临界雷诺数代数平均 D.上下临界雷诺数几何平均 24..对于管道无压流,当充满度分别为( )时，其流量和速度分别达到最大。Ａ．0.５,0.５B．0.95,0．８1 C.0.８1, ０8１ D. １.0，1.０ 25.对于ａ, ｂ，c三种水面线，下列哪些说法是错误( ）（) A．所有a、ｃ型曲线都是壅水曲线，即，水深沿程增大。B．所有
《计算流体力学》结课作业要点.doc
2012～2013学年第1学期 12级研究生《计算流体力学》结课作业 适用专业：供热供燃气通风及空调工程 一、结合某一具体学科，阐述纯理论方法、实验方法及数值方法在科学研究中的各自优缺点，在此基础上论述数值模拟方法的发展前景。（不少于4千字）。 流体力学是力学的一个重要分支, 是研究流体(液体和气体)的力学运动规律及其应用的学科, 主要研究在各种力的作用下，流体本身的静止状态和运动状态特征，以及流体和相邻固体界面有相对运动时的相互作用和流动规律。在人们的生活和生产活动中随时随地都可遇到流体，流体力学与人类的日常生活和生产事业密切相关。按其研究内容的侧重点不同，分为理论流体力学和工程流体力学。其中理论流体力学主要采用严密的数学推理方法，力求准确性和严密性，工程流体力学侧重于解决工程实际中出现的问题，而不追求数学上的严密性。当然由于流体力学研究的复杂性，在一定程度上，两种方法都必须借助于实验研究，得出经验或半经验的公式。 在实际工程的诸多领域流体力学都起着十分重要的作用。如气象、水利的研究，船舶、飞行器、叶轮机械和核电站的设计及其运行，可燃气体或炸药的爆炸，都广泛地用到流体力学知识。许多现代科学技术所关心的问题既受流体力学的指导，同时也促进了流体力学自身的不断发展。1950年后，计算机的发展给予流体力学以极大的推动作用。 目前，解决流体力学问题的方法主要有实验方法、理论分析方法和数值方法三种。 实验方法 同物理学、化学等学科一样，流体力学的研究离不开实验，尤其是对新的流体运动现象的研究。实验能显示运动特点及其主要趋势，有助于形成概念，检验理论的正确性。二百年来流体力学发展史中每一项重大进展都离不开实验。流体力学实验研究方法有实物实验、比拟研究和模型研究三类：实物实验是用仪器实测原型系统的流动参数，适用于较小的原型；比拟实验是利用电场和磁场来模拟流场，实施起来限制条件较多；模型研究是实验流体力学最常用的研究方法。 实验研究的一般过程是：在相似理论的指导下建立实验模型，用流体测量技术测量流动参数，处理和分析实验数据。建立实验模型要求模型与原型满足相似理论，即满足两个流场
计算流体力学大作业报告(翼型空气动力分析)
课程综合作业课程名称：计算流体力学 专业班级：研究方向： 学生姓名：学号： 完成日期：
计算流体力学课程综合报告 1.简介 计算流体动力学（Computational Fluid Dynamics，简称CFD）是通过计算机数值计算和图像显示，对包含有流体流动和热传导等相关物理现象的系统所做的分析。其基本思想为:把原来在时间域及空间域上连续的物理量的场，如速度场和压力场，用一系列有限个离散点上的变量值的集合来代替，通过一定的原则和方式建立起关于这些离散点上场变量之间关系的代数方程组，然后求解代数方程组获得场变量的近似值。 CFD可以看作是在流动基本方程（质量守恒方程、动量守恒方程、能量守恒方程）控制下对流动的数值模拟。通过这种数值模拟，我们可以得到极其复杂问题的流场内各个位置上的基本物理量（速度、压力、温度、浓度等）的分布，以及这些物理量随时间的变化情况，确定旋涡分布特性、空化特性及脱流区等。还可据此算出相关的其他物理星，如旋转式流体机械的转矩、水力损失和效率等。此外，与CAD联合，还可进行结构优化设计等。 2.计算流体动学的特点: ①流动问题的控制方程一般是非线性的，自变量多，计算域的几何形状和边界条件复杂，很难求得解析解，而用CFD方法则有可能找出满足工程需要的数值解。 ②可利用计算机进行各种数值试验，例如，选择不同流动参数进行物理方程中各项有效性和敏感性试验，从而进行方案比较。 ③它不受物理模型和实验模型的限制，省钱省时，有较多的灵活性，能给出详细和完整的资料，很容易模拟特殊尺寸、高温、有毒、易燃等真实条件和实验中只能接近而无法达到的理想条件。 ④数值解法是一种离散近似的计算方法，依赖于物理上合理、数学上适用、适合于在计算机上进行计算的离散的有限数学模型，且最终结果不能提供任何形式的解析表达式，只是有限个离散点上的数值解，并有一定的计算误差。 ⑤它不像物理模型实验一开始就能给出流动现象并定性地描述，往往需要由原体观测或物理模型试验提供某些流动参数，并需要对建立的数学模型进行验证。
计算流体力学作业
计算流体力学课程作业 任课教师：魏文礼 姓名： 学号： 指导老师：
目录 1.写出通用方程，并说明如何代表各类守恒方程。 (1) 2.推导流体运动的质量、动量守恒方程。 (2) 3.简述源项线性化、网格划分问题。 (5) 4.用ddxKeTex+S=0，谈谈边界条件如何处理。 (8) 5.用有限体积法离散ρceTet=eexKeTex，并推广到二维、三维问题，写出过程。 (9) 6.从不同角度对流体运动分类。 (12) 7.谈谈物理模型试验与计算流体力学方法的关系。 (12) 8.讨论离散对流项时离散格式的进化过程。 (13) 9.利用幂函数格式离散二维、三维通用方程的离散方程。 (15) 10.解释交错网格的概念。 (15) 11.简述压力校正法解N-S方程的过程。 (16) 12.思考anbvnb′为什么可以省去。 (17)
1.写出通用方程，并说明如何代表各类守恒方程。 答：(1)写出通用方程。 在Cartesian坐标系下单位体积黏性流动N-S方程组微分形式如下： { eρ et +??(ρV)=0 (1) eρu +??(ρuV)=??(μ?u)+ 1 μ[ e (??V)]? ep +F x+S mx(2a)eρv et +??(ρvV)=??(μ?v)+ 1 3 μ[ e ey (??V)]? ep ey +F y+S my(2b) eρw et +??(ρwV)=??(μ?w)+ 1 3 μ[ e ez (??V)]? ep ez +F z+S mz(2c)eρe et +??(ρeV)=??(k?T)?p(??V)+Φ+Q (3) 上述微分形式黏性流动N-S方程组中，式(1)为连续性方程，式(2a)、 (2b)、(2c)分别为x、y、z方向上的动量方程，式(3)为能量方程。 上述方程组中各个方程具有不同变量，代表不同的守恒定律，但他们的形式都十分相似。若引入一个通用的特征变量?，在不同的方程中?代表不同的变量，就可以把它写为通用变量形式。 非定常通用变量N-S方程为： e(ρ?) et +??(ρ?V)=??(Γ???)+S? 若流场中速度等物理量不随时间变化，则e(ρ?) et =0，可得定常通用变量N-S方程为： ??(ρ?V)=??(Γ???)+S? 其中，?为通用变量，可代表u、v、w、T等求解变量；Γ?为扩散和热传导系数，S?为方程组源项。 (2)用通用方程代表各类守恒方程 用通用方程代表各类守恒方程是，通用变量在各守恒方程中的取值如表1所示。
流体力学-中国科学院海洋研究所研究生部
中科院海洋研究所硕士研究生入学考试 《流体力学》考试大纲 本流体力学考试大纲适用于中国科学院研究生院力学专业的硕士研究生入学考试。流体力学是现代力学的重要分支，是许多学科专业的基础理论课程，本科目的考试内容主要包括流体的物理性质、流体运动学、动力学和静力学，无粘不可压缩、可压缩流动，粘性不可压缩流动及湍流、流体波动和漩涡理论等方面。要求考生对其基本概念有较深入的了解，能够熟练地掌握基本方程的推导，并具有综合运用所学知识分析问题和解决问题的能力。 一、考试内容： （一）流体的物理性质 固液气体的宏观性质与微观结构，连续介质假设及其适用条件，流体的物理性质(粘性、可压缩性与热膨胀性、输运性质、表面张力与毛细现象) ，质量力与表面力。 （二）流体运动学 流体运动的描述(拉格朗日描述与欧拉描述及其间的联系、物质导数与随体导数、迹线、流线及脉线)，流场中的速度分解，涡量，涡量场，涡线、涡管、涡通量，涡管强度及守恒定理。 （三）流体动力学 连续性方程(雷诺输运定理)，动量方程(流体的受力、应力张量)，能量方程(热力学定律)，本构关系，状态方程，流体力学方程组及定解条件，正交曲线坐标系，量纲分析与流动相似理论，流体力学中的无量纲量及其物理意义、相似原理的应用。 （四）流体静力学 控制方程，液体静力学规律，自由面的形状，非惯性坐标系中的静止液体。 （五）无粘流动的一般理论 无粘流动的控制方程，Bernoulli方程，Bernoulli方程和动量定理的应用。 （六）无粘不可压缩流体的无旋流动 控制方程及定解条件，势函数及无旋流动的性质，平面定常无旋流动(流函数、源汇、点涡、偶极子、镜像法、保角变换)，无旋轴对称流动，非定常无旋流动。 （七）液体表面波 控制方程(小振幅水波) 及定解条件，平面单色波，水波的色散和群速度，水波的能量及其传输，速度与压力场特性，表面张力波及分层流体的重力内波，非线性水波理论。 （八）旋涡运动 涡量动力学方程和涡量的产生，涡量场(空间特性、时间特性)，典型的涡模型。 （九）粘性不可压缩流动 控制方程及定解条件，定常的平行剪切流动(Couette流动、Poiseuille流动等)，非定常的平行剪切流动(Stokes第一和第二问题、管道流动的起动问题)，圆对称的平面粘性流动(圆柱Couette流及其起动过程)，小雷诺数粘性流动。 （十）层流边界层和湍流 边界层的概念，层流边界层方程(Blasius平板边界层)，边界层的分离，湍流的发生，层流到湍流的转捩，雷诺方程和雷诺应力。 （十一）无粘可压缩流动 声速和马赫数，膨胀波、弱压缩波的形成及其特点，一维等熵流(定常和非定常)，激波(正激波和斜激波)，拉瓦尔喷管流动的特征。 二、考试要求：
流体力学试验
流體力學實驗 老師：Ａ1班→李宗翰老師；Ｂ1班→楊龍杰老師 Ａ2班→蔡欣正老師；Ｂ2班→邵德文老師 時間：Ａ1班→星期五12、13、14節；Ｂ1班→星期三11、12、13節Ａ2班→星期二11、12、13節；Ｂ2班→星期四12、13、14節 上課進度：
成績計算： 1.作業（30﹪）：上課後10分鐘未交報告者扣總分3分！ 當日無故未交者扣總分10分！ 1.課堂（20﹪）：分組合作精神，數據結果及隨堂口試小考。 2.口試（25﹪）：於考前一週公告口試方式。 3.筆試（25﹪）：於考前一週公告考場。 4.上課遲到10分鐘內扣總分3分！ 無故缺課扣總分10分！缺課3次下學期再見！ ※實驗前每組須備有空白數據表格一份，以方便記錄實驗數據※ 規定事項： 一、預習報告：（限用Ａ4大小的紙書寫，不可用打字） 1.封面：包含實驗名稱、組別、班級、姓名、學號、座號。 2.內容：包含實驗目的、實驗原理、實驗步驟及空白數據表格。 3.每人一份，於實驗前由組長收齊交給助教簽章，並於批閱後取回。 二、結論報告：（限用Ａ4大小的紙書寫） 1.個人結報：每人一份，含實驗心得和討論（心得須300字以上）。上課前將 個人結報及前一次實驗領回的預報合訂在一起，交給組長。 2.整組結報：每組一份，含數據、回歸分析結果，回歸分析圖表。 3.回歸分析須有電腦分析報表結果和座標曲線圖，圖可用手畫或電腦處理，若 用手畫請用方格紙，不可用工學院作業紙的背面。 4.未能及時繳交之作業，也一定要儘快繳交，不可缺交。 三、上課期間： 1.在實驗室內不可抽煙、進食及喝水，並注意安全。 2.不可無故離開實驗室，如有需要請先報備，助教會不定時的抽點。 3.組長負責整組的實驗操作、秩序及做完實驗後的清潔。 4.實驗後的數據表格，須在下課前交給助教檢查才算完成，嚴禁抄襲。
计算流体力学与传热学大作业
########学院 计算流体力学与传热学 学号： 专业： 学生姓名： 任课教师：教授 2013年12月
目录 第一章验证显式格式的稳定性 (4) 1.1 概述 (4) 1.2 数学推导 (4) 1.3 问题描述 (4) 1.4 数值模拟 (4) 1.5 结果及分析 (5) 第二章判断肋片可以按一维问题处理的主要依据 (6) 2.1 概述 (6) 2.2 问题描述及算法 (6) 2.3 数值模拟 (7) 2.4 结果及分析 (8) 第三章三层墙导热 (9) 3.1 概述 (9) 3.2 问题描述 (9) 3.3 TDMA算法 (9) 3.4 结果 (10) 第四章一维无源稳态对流扩散问题 (11) 4.1 公式及初值 (11) 4.2 情况一 (11) 4.3 情况二 (12) 4.4 情况三 (13)
第五章用ADI算法计算长方肋内的温度分布 (14) 5.1 问题描述 (14) 5.2 初始参数 (14) 5.3 情况一，一列列扫 (14) 5.4 情况二，一行行扫 (14) 5.5 情况三，采用ADI算法 (15) 5.6 结果分析 (15) 参考文献 (16)
第一章 验证显式格式的稳定性 1.1 概述 将一维非稳态热传导方程用显式格式差分化为代数方程，在求解的迭代过程中必须满足一定的条件，才能使方程收敛且结果正确。此处即验证β≤?。 1.2 数学推导 方程： 22T t T x α??=?? （1） 显式离散格式： 此处时间向前差分，空间中心差分 111 22n n n n n i i i i i T T T T T t x α+-+--+=?? 1112(2)n n n n n i i i i i t T T T T T x α +-+?-=-+? 令β=2 t x α ??则： 111(2)n n n n n i i i i i T T T T T β+-+-=-+ （2） 误差也应该满足上式，故： ()()1()()()2()()i i i i i Ikx Ikx Ik x x Ikx Ik x x n n n n n T e T e T e T e T e ψψβψψψ----?--+?+??-=-+?? ()()()1()12()()()i i i i Ikx Ikx Ik x x Ik x x n n n n T e T e T e T e ψβψβψψ----?-+?+??=-++?? ()()1()12()()i i i Ikx Ikx Ikx n n Ik x Ik x n T e T e e e T e ψβψβψ---+-??=-++ ()()1() 121() n Ik x Ik x n T e e T ψββψ+-??=-++≤ 因此 β≤?。即当β≤? 时方程（2）才会有收敛的解。 1.3 问题描述 在验证过程中同时可模拟一个实际问题，即冬季里墙壁中的温度分布。此时室内壁温设为Tl=30.0℃，室外壁温Tr=-25.0℃，墙壁以11号楼为例，L=1m ，热扩散系数ɑ=alfa=1.33e-6m 2/s 然后分别取β=0.4，n=10和β=0.6，n=10两种情况，看最后的结果是否收敛和正确。 1.4 数值模拟
计算数学研究方向
计算数学研究方向 网上摘抄：计算数学研究方向及网上资料 计算数学目的为物理学和工程学作计算。主要研究方向包括： 数值泛函分析；连续计算复杂性理论；数值偏微与有限元；非线性数值代数及复动力系统； 非线性方程组的数值解法；数值逼近论；计算机模拟与信息处理等；工程问题数学建模与计算等等。 目前发展最好的方向已经与应用数学的CAGD 方向合二为一。现在最热的方向应该是微分方程的数值求解、数值代数和流形学习，数值计算名校：西安交通大学、北京大学、大连理工大学 从计算数学的字面来看，应该与计算机有密切的联系，也强调了实践对于计算数学的重要性。 也许Parlett 教授的一段话能最好地说明这个问题： How could someone as brilliant as von Neumann think hard about a subject as mundane as triangular factoriz-ation of an invertible matrix and not perceive that, with suitable pivoting, the results are impressively
good Partial answers can be suggested-lack of hands-on experience, concentration on the inverse rather than on the solution of Ax = b -but I do not find them adequate. Why did Wilkinson keep the QR algorithm as a backup to a Laguerre-based method for the unsymmetric eigenproblem for at least two years after the appearance of QR Why did more than 20 years pass before the properties of the Lanczos algorithm were understood I believe that the explanation must involve the impediments to comprehension of the effects of finite-precision arithmetic. ( 引自既然是计算数学专业的学生，就不能对自己领域内的专家不有所了解。早些年华人在计算数学领域里面占有一席之地是因为冯康院士独立于西方，创立了有限元方法，而后又提出辛算法。这里只是列出几位比较年轻的华人计算数学专家，因为他们代表了当前计算数学的研究热点，也反映华人对计算数学的发展的贡献。 侯一钊（加州理工） 研究方向：计算流体力学、多尺度计算与模拟、多相流 鄂维南（Princeton 大学） 北京大学长江学者，研究方向：多尺度计算与模拟 包刚（Michigan 州立大学） 吉林大学长江学者，研究方向：光学与电磁场中的计算等 金石(Wisconsin 大学)