【免费下载】计算水动力学2009重点
计算流体动力学分析-CFD软件原理与应用_王福军--阅读笔记
计算流体动力学(简称CFD)是建立在经典流体动力学与数值计算方法基础之上的一门新型独立学科，通过计算机数值计算和图像显示的方法，在时间和空间上定量描述流场的数值解，从而达到对物理问题研究的目的。它兼有理论性和实践性的双重特点。 第一章节 流体流动现象大量存在于自然界及多种工程领域中，所有这些过程都受质量守恒、动量守恒和能量守恒等基本物理定律的支配。本章向读者介绍这些守恒定律的数学表达式，在此基础上提出数值求解这些基本方程的思想，阐述计算流体力学的任务及相关基础知识，最后简要介绍目前常用的计算流体动力学商用软件。 计算流体动力学((Computational Fluid Dynamics简称CFD)是通过计算机数值计算和图像显示，对包含有流体流动和热传导等相关物理现象的系统所做的分析。CFD的基本思想可以归结为:把原来在时间域及空间域上连续的物理量的场，如速度场和压力场，用一系列有限个离散点上的变量值的集合来代替，通过一定的原则和方式建立起关于这些离散点上场变量之间关系的代数方程组，然后求解代数方程组获得场变量的近似值。 CFD可以看做是在流动基本方程(质量守恒方程、动量守恒方程、能量守恒方程)控制卜对流动的数值模拟。通过这种数值模拟，我们可以得到极其复杂问题的流场内各个位置上的基本物理量(如速度、压力、温度、浓度等)的分布，以及这些物理量随时间的变化情况，确定旋涡分布特性、空化特性及脱流区等。还可据此算出相关的其他物理量，如旋转式流体机械的转矩、水力损失和效率等。此外，与CAD联合，还可进行结构优化设计等。 1.1.2计算流体动力学的工作步骤 采用CFD的方法对流体流动进行数值模拟，通常包括如下步骤: (1)建立反映工程问题或物理问题本质的数学模型。具体地说就是要建立反映问题各个量之间关系的微分方程及相应的定解条件，这是数值模拟的出发点。没有正确完善的数 学模型，数值模拟就毫无意义。流体的基本控制方程通常包括质量守恒方程、动量守恒方程、能量守恒方程，以及这些方程相应的定解条件。 (2}}寻求高效率、高准确度的计算方法，即建立针对控制方程的数值离散化方法，如有限差分法、有限元法、有限体积法等。这里的计算方法不仅包括微分方程的离散化方法及求解方法，还包括贴体坐标的建立，边界条件的处理等。这些内容，可以说是c}}的核心。 (3})编制程序和进行计算。这部分工作包括计算网格划分、初始条件和边界条件的输入、控制参数的设定等。这是整个工作中花时间最多的部分。由于求解的问题比较复杂，比如Na}ier-Stakes方程就是一个讨，分复杂的非线性方程，数值求解方法在理论上不是绝对完善的，所以需要通过实验加以验证。正是从这个意义上讲.数值模拟又叫数值试验。应该指出，这部分工作不是轻而易举就可以完成的。 4})显示计算结果。计算结果一般通过图表等方式显示，这对检查和判断分析质量和结果有重要参考意义。 以上这些步骤构成了CFD数值模拟的全过程。其中数学模型的建立是理论
流体力学复习要点(计算公式)
D D y S x e P gh2 gh1 h2 h1 b L y C C D D y x P hc 第一章 绪论 单位质量力： m F f B m = 密度值： 3 m kg 1000=水ρ， 3 m kg 13600=水银ρ， 3 m kg 29.1=空气ρ 牛顿内摩擦定律：剪切力： dy du μ τ=， 内摩擦力：dy du A T μ= 动力粘度： ρυ μ= 完全气体状态方程：RT P =ρ 压缩系数： dp d 1dp dV 1ρρκ= -=V （N m 2 ） 膨胀系数：T T V V V d d 1d d 1ρρα - == (1/C ?或1/K) 第二章 流体静力学+ 流体平衡微分方程： 01;01;01=??-=??-=??- z p z y p Y x p X ρρρ 液体平衡全微分方程：)(zdz ydy xdx dp ++=ρ 液体静力学基本方程：C =+ +=g p z gh p p 0ρρ或 绝对压强、相对压强与真空度：a abs P P P +=；v a abs P P P P -=-= 压强单位换算：水银柱水柱mm 73610/9800012 ===m m N at 2/101325 1m N atm = 注： h g P P →→ρ ； P N at →→2m /98000乘以 2/98000m N P a = 平面上的静水总压力：（1）图算法 Sb P = 作用点e h y D +=α sin 1 ) () 2(32121h h h h L e ++= ρ 若01 =h ，则压强为三角形分布，3 2L e y D == ρ 注：①图算法适合于矩形平面；②计算静水压力首先绘制压强分布图， α 且用相对压强绘制。 （2）解析法 A gh A p P c c ρ== 作用点A y I y y C xc C D + = 矩形12 3 bL I xc = 圆形 64 4 d I xc π= 曲面上的静水总压力： x c x c x A gh A p P ρ==；gV P z ρ= 总压力z x P P P += 与水平面的夹角 x z P P arct an =θ 潜体和浮体的总压力： 0=x P 排浮gV F P z ρ== 第三章 流体动力学基础 质点加速度的表达式??? ? ? ? ??? ??+??+??+??=??+??+??+??=??+??+??+??=z u u y u u x u u t u a z u u y u u x u u t u a z u u y u u x u u t u a z z z y z x z z y z y y y x y y x z x y x x x x A Q V Q Q Q Q Q G A = === ? 断面平均流速重量流量质量流量体积流量g udA m ρρ 流体的运动微分方程： t z t y t x d du z p z d du y p Y d du x p X = ??-=??-=??- ρρρ1;1;1 不可压缩流体的连续性微分方程 ： 0z u y u x u z y x =??+??+?? 恒定元流的连续性方程： dQ A A ==2211d u d u 恒定总流的连续性方程：Q A A ==2211νν 无粘性流体元流伯努利方程：g 2u g p z g 2u g p z 2 2 222 111++=++ρρ 粘性流体元流伯努利方程： w 2 2222111'h g 2u g p z g 2u g p z +++=++ρρ
高考物理二轮复习计算题题型1运动学、动力学类问题练习
计算题题型1 运动学、动力学类问题 角度1:直线运动规律及牛顿运动定律的综合应用 1.(2017·江西吉安一诊)如图所示,在赛车训练场相邻两车道上有黑白两辆车,黑色车辆停在A线位置,某时刻白色车速度以v1=40 m/s通过A线后立即以大小a1=4 m/s2的加速度开始制动减速,黑车4 s后开始以a2=4 m/s2的加速度开始向同一方向匀加速运动,经过一定时间,两车同时在B线位置.两车看成质点.从白色车通过A线位置开始计时,求经过多长时间两车同时在B线位置及在B线位置时黑色车的速度大小. 2.质量M=10 kg的木板A沿水平面向右运动,与水平面之间的动摩擦因数μ1=0.1,当A的速度v0=5 m/s时,在A的左端施加一个恒力F=35 N,如图所示,同时在木板上表面无初速度地放上一个质量m=5 kg的滑块B.已知滑块B右端的木板上表面粗糙,长度为12.5 m,与滑块之间的动摩擦因数μ2=0.1,滑块左端的木板上表面包括滑块所放的位置均光滑,长度为 2.5 m,g 取10 m/s2. (1)至少经过多长时间滑块与木板的速度相等? (2)共经过多长时间滑块与木板分开? 3.(2017·辽宁鞍山一模)如图所示为在某工厂的厂房内用水平传送带将工件的半成品运送到下一工序的示意图.传送带在电动机的带动下保持v=2 m/s的速度匀速向右运动,现将质量
为m=20 kg的半成品轻放在传送带的左端A处,半成品工件与传送带间的动摩擦因数为μ=0.5,设传送带足够长,重力加速度g=10 m/s2.试求: (1)半成品工件与传送带相对滑动所经历的时间; (2)半成品工件与传送带间发生的相对位移大小; (3)若每分钟运送的半成品工件为30个,则电动机对传送带做功的功率因运送工件而增加多少? 角度2:带电粒子(带电体)在电场与磁场中的平衡与运动 1.(2017·黑龙江双鸭山一模)如图所示,一带电荷量为+q、质量为m的小物块处于一倾角为37°的光滑斜面上,当整个装置被置于一水平向右的匀强电场中,小物块恰好静止.重力加速度取g,sin 37°= 0.6,cos 37°=0.8.求: (1)水平向右电场的电场强度; (2)若将电场强度减小为原来的,物块的加速度是多大? (3)电场强度变化后物块下滑距离L时的动能.
《计算流体动力学分析》学习报告
《计算流体动力学分析》学习报告 计算流体力学基础： 本章主要讲解流体动力学的核心思想以及流体动力学的控制方程。 1、计算流体动力学（Computational Fluid Dynamic ）基本思想：把原来在时间和空间上的连续的物理量，用一系列离散点上的变量值来代替，通过一定的原则和方式建立变量之间的代数方程式，求解之后获得变量的近似值。 2、CFD 控制方程： 质量守恒方程 0)·=?+??u t ρρ（ 动量守恒方程（Navier-Stokes 方程） Fz z y x z u w div t w F z y x y u v div t v F z y x x u u div t u zz zx zx y zy yy xy x zx yx xx +??+??+??+??-=+??+??+??+??+??-=+??+??+??+??+??-=+??τττρρρτττρρρτττρρρ)()()()()()( 能量守恒方程 T p S gradT c k div T u div t +=+??)()(T ( ρρ） S T 为粘性耗散项。 方程含有u ，v ，w ，p ，T 和ρ六个未知量，所以还需要一个方程组，才能使其封闭，而这个方程组就是联系P 和ρ的状态方程组：P=（ρ，T ）。 组分质量守恒方程（在一个系统中，可能存在质的交换，或者存在化学组分时使用。） ()s s s s S c grad D div c u div t +=+??)()(c (s ρρρ ） 为便于对控制方程进行计算和分析，对CFD 控制方程写成通用格式： ()S z z y y x x z w y v x u t S grad div u div t +??Γ??+??Γ??+??Γ??=??+??+??+??+Γ=+??)()()()()()())()(φφφφρφρφρρφφφρρφ 依次为瞬态项，对流项，扩散项和源项。 3、湍流控制方程 三维的N-S 方程无论对于层流还是湍流都是是使用的，但由于直接求解三维瞬态的控制方程，对计算机的内存和速度要求很高，因此在工程上广为采用的方法是对瞬态的N-S 方程进行实践平均处理，同时补充反应湍流特性的其他方程，例如湍动能方程以及湍流耗散率方程
地下水动力学其它计算题
第二章计算题 1. 在厚度不等的承压含水层中，沿地下水流方向打四个钻孔（孔1、孔2、 孔3、孔4），如图2—1所示，各孔所见含水层厚度分别为：M 1=14.5，M 2 =M 3 =10m， M 4 =7m，已知孔1—孔2、孔2—孔3、孔3—孔4的间距分别为210m、125m、180m。 2. 图2—2所示，作侧河水已受污染，其水位用H 1 表示，没有受污染的右侧 河水位用H 2 表示。（1）已知河渠间含水层为均质、各向同性，渗透系数未知，
入渗强度为0.01m/d。当含水层中水位至少下降2m时，两侧排水渠水位都为H=6m。试求：（1）排水渠的间距L；（2）排水渠一侧单位长度上的流量Q。 两河间距l=500m，含水层的稳定单宽流量为1.2m2/d。在无入渗补给量的条件下，试求含水层的渗透系数。 5. 水文地质条件如图2—4所示。已知h 1=10m，H 2 =10m，下部含水层的平均 厚度M=20m，钻孔到河边距离l=2000m，上层的渗透系数K 1 =2m/d，下层的渗透系 A B 水层分为上下两层，上层为细砂，A、B两处的含水层厚度分别为h A =5.19m、
h =2.19m，渗透系数为3.6m/d。下层为粗砂，平均厚度M=6.4m，渗透系数为30m/d。B 试求含水层的单宽流量。 7. 图2—5所示，某河旁水源地为中粗砂潜水含水层，其渗透系数为100m/d。 含水层平均厚度为20m，给水度为0.002。以井距30m的井排进行取水，井排与 8. 某水库蓄水后，使岸边潜水产生回水现象，如图2—6所示。设计水库蓄 水后最高水位标高H=28m。在距水库l=5km处有一工厂，其地面标高为25m，已 =8m（以含水层底版算起），渗透系数为10m/d，给水度为0.04。设计灌渠水h 位瞬时抬高1.5m后，使地下水位在一天内最小抬高0.3m。试求灌渠的合理间距。
第07章基元反应动力学习题及答案
第07章基元反应动力学习题及答案
2 第七章 基元化学反应动力学 习题及答案 1. N 2O 5在25℃时分解反应的半衰期为5.7h, 且与N 2O 5的初始压力无关。 试求此反应在25℃条件下完成90%所需时间。 解：由题意知此反应为一级反应 1 1 1216.07 .56932 .06932.02 1-=== h t k t k y 1)11 ln( =- h k y t 9.181216.0/)% 9011ln(/)11ln( 1=-=-= 即完成90％所需时间为18.9h 。 2．异丙烯醚气相异构化成丙烯酮的反应是一级反应，其反应速率系(常)数与温度的关系为：k /s －1 =5.4×1011exp(－122 474 J ·mol －1/RT )，150℃下，反应开始时只有异丙烯醚，其压力为101 325 Pa ，问多长时间后，丙烯酮的分压可达54 kPa ？ 解：k /S －1＝5.4×1011 exp[－122474/8.314×（150＋273）] ＝4.055×10－4 据题意：kt p p t =0 ln
3 t 410005.454000 101325101325 ln -?=- t ＝1877S 3. 双分子反应2A(g)?→?k B(g) + D(g),在623K 、初始浓度为0.400mol dm －3时, 半衰期为105s,请求出 (1) 反应速率系数k (2) A(g)反应掉90%所需时间为多少? (3) 若反应的活化能为140 kJ mol -1, 573K 时的最大反应速率为多少? 解：(1) r = k [A]2 , t 0.5= 1/(2 k [A]0) , k = 0.012dm 3mol -1s -1 (2) 1/[A] – 1/[A]0 =2 k t , t = 945 s (3) ln(k/k ’)=(E a /R )(1/T ’-1/T ) , 573K 时k = 0.00223dm 3mol -1s -1, 最大反应速率r max = k [A]02=3.6×10－4 moldm -3s -1. 4． 450℃时实验测定气相反应3A + B →2C 的速率数据如下； 实验 初压 / Pa 初速率－dp B / dt / (Pa/h) P A,0 P B,0 1． 100 1.00 0.0100 2． 200 1.00 0.0400 3． 400 0.50 0.0800 (1)若反应的速率方程为r = kP A x P B y ，求x 、y 及k 。 (2)求P A =150Pa P B = 0.75Pa 时反应的速率。 解：（1）把1、2两组数据分别代入速率方程r ＝kp A X ·p B Y 并相比，得： r 1/r 2＝(100/200)X 即 0.0100/0.0400＝（100/200）X
药代动力学代表计算题
计算题（Calculation questions ） 1．某患者单次静脉注射某单室模型药物2g ，测得不同时间的血药浓度结果如下： 时间(h) 1.0 2.0 3.0 4.0 5.0 6.0 8.0 10.0 血药浓度(mg/ml) 0.28 0.24 0.21 0.18 0.16 0.14 0.1 0.08 求k ，Cl ，T 1/2，C 0，V ，AUC 和14h 的血药浓度。 【解】对于单室模型药物静脉注射 kt 0e C C -=，t 303 .2k C log C log 0- = log C 对t 作直线回归（注：以下各题直线回归均使用计算器或计算机处理），得： a = 0.4954, b = -0.0610,|r | = 0.999（说明相关性很好） 将a 、b 代入公式0C log 303 .2kt C log +-= 得回归方程： 4954.0t 061.0C log --= ① 1h 1405.0)061.0(303.2b 303.2k -=-?-=?-= ② h 9323.41405 .0693 .0k 693.0T 2/1=== ③ mg/ml 3196.0)4954.0(log C 10=-=- ④ 6.258L ml)(62583196 .02000C X V 00==== ⑤ L/h 8792.0258.61405.0kV Cl =?== ⑥ )(mg/ml h 2747.21405 .03196.0k C AUC 00 ?=== ∞ ⑦ 3495.14954.014061.0C log -=-?-= g/ml 44.7mg/ml)(0477.0C μ== 即14h 的血药浓度为g/ml 44.7μ。 2．某患者单次静脉注射某药1000mg ，定期测得尿药量如下： 时间(h) 1 2 3 6 12 24 36 48 60 72 每次尿药量(mg) 4.02 3.75 3.49 9.15 13.47 14.75 6.42 2.79 1.22 0.52 设此药属一室模型，表观分布容积30L ，用速度法求k ，T 1/2，k e ，Cl r ，并求出80h 的累积药量。 【解】单室模型静脉注射尿药数据符合方程0e c u X k log 303 .2kt t X log +-=??， t X log u ??对c t 作图应为一直线。根据所给数据列表如下： t (h) 1 2 3 6 12 t ? 1 1 1 3 6
CFD—计算流体动力学软件介绍
CFD 流体动力学软件介绍 CFD—计算流体动力学，因历史原因，国一直称之为计算流体力学。其结构为： 提出问题—流动性质（流、外流；层流、湍流；单相流、多相流；可压、不可压等等），流体属性（牛顿流体：液体、单组分气体、多组分气体、化学反应气体；非牛顿流体） 分析问题—建模—N-S方程（连续性假设），Boltzmann方程（稀薄气体流动），各类本构方程与封闭模型。 解决问题—差分格式的构造/选择，程序的具体编写/软件的选用，后处理的完成。 成果说明—形成文字，提交报告，赚取应得的回报。 CFD实现过程： 1.建模——物理空间到计算空间的映射。 主要软件： 二维： AutoCAD： 大家不要小看它，非常有用。一般的网格生成软件建模都是它这个思路，很少有参数化建模的。相比之下AutoCAD的优点在于精度高，草图处理灵活。可以这样说，任何一个网格生成软件自带的建模工具都是非参数化的，而对于非参数化建模来说，AutoCAD应该说是最好的，毕竟它发展了很多很多年！ 三维： CATIA：航空航天界CAD的老大，法国人的东西，NB，实体建模厉害，曲面建模独步武林。本身可以生成有限元网格，前几天又发布了支持ICEM-CFD的插件ICEM-CFD CAA V5。有了它和ICEM-CFD，可以做任何建模与网格划分！ UG：总觉得EDS脑袋进水了，收了I-deas这么久了，也才发布个几百M的UG NX 2.0，还被大家争论来争论去说它如何的不好用！其实，软件本身不错，大公司用得也多，可是就这么打市场，早晚是走下坡路。按CAD建模的功能来说它排不上第一，也不能屈居第二，尤其是加上了I－DEAS更是如虎添翼。现
流体力学计算公式
1、单位质量力：m F f B B = 2、流体的运动粘度：ρ μ=v （μ[动力]粘度，ρ密度） 3、压缩系数：dp d dp dV V ρρκ?=?-=11（κ的单位是N m 2）体积模量为压缩系数的倒数 4、体积膨胀系数：dT d dT dV V v ρρα?-=?=11（v α的单位是C K ?1,1） 5、牛顿内摩擦定律：为液体厚）为运动速度，以应力表示为y u dy du dy du A T (,μτμ== 6、静止液体某点压强：为该点到液面的距离）h gh p z z g p p ()(000ρρ+=-+= 7、静水总压力： )h (为受压面积，为受压面形心淹没深度为静水总压力，A p ghA A p p c ρ== 8、元流伯努利方程;'2221112w h g p z g u g p z ++=++ρρ('w h 为粘性流体元流单位重量流体由过流断面1-1运动至过流断面2-2的机械能损失，z 为某点的位置高度或位置水头，g p ρ为测压管高度或压强水头，g u ρ2是单位流体具有的动能，u gh g p p g u 22'=-=ρ，u gh C g p p g C u 22'=-=ρC 是修正系数，数值接近于1） 9、总流伯努利方程:w h g v g p z g v g p z +++=++222 221221111αραρ(α为修正系数通常取1） 10、文丘里流量计测管道流量：)21)(41()()(42 122211g d d d k h k g p z g p z k Q -=?=+-+=πμρρμ 11、沿程水头损失一般表达式：g v d l h f 22 λ=（l 为管长，d 为管径，v 为断面平均流速，g
化学动力学计算题
第十一章 化学动力学计算题 1．放射性同位素的蜕变速率符合一级反应的规律(蜕变速率与放射性同位素的数量成正比)。210Po 经α蜕变生成稳定的206Pb ：210Po → 206Pb + 4He 。实验测得14 d 后放射性降低了6.85%，试求210Po 的蜕变速率常数和半衰期，并计算它蜕变掉90%时所需要的时间。 解: 2．双分子反应2A(g)?→?k B(g) + D(g), 在623 K 、初始浓度为0.400 mol·dm -3时,半衰期为105 s, 求出: (1) 反应速率常数k ; (2) A(g)反应掉90%所需时间为多少? (3) 若反应的活化能为140 kJ/mol, 573 K 时的最大反应速率为多少? 解：(1) 由基元反应得: v = k [A]2 那么: t 1/2= 1/(2k [A]0), k = 1/(2t 1/2[A]0) =1/(2×105×0.400) dm 3·mol -1·s -1 =0.0120 dm 3·mol -1·s -1 (2) 由1/[A] – 1/[A]0 =2kt 得: 1/0.1[A]0 – 1/[A]0 =2kt 1/0.0400 – 1/0.400 = 2×0.0120×t 所以: t = 945 s (3) 由ln(k/k ')=(E a /R )×(1/T '-1/T ) 得: ln(0.012/k ')=(140000/8.314)×(1/573-1/623) k ' = 0.00223 dm 3·mol -1·s -1, 因此, 最大反应速率v max = k '[A]02= 0.00223×(0.400)2 mol·dm -3·s -1 =3.6×10-4 mol·dm -3·s -1 1 3d 1007.50685.011ln d 14111ln 1--?=-=-=A A x t k d 137d 1007.52ln 2ln 1 32/1=?==--A k t d 45490.011ln d 1007.5111ln 113=-?=-=--A A x k t
第二章计算流体力学的基本知识
第二章计算流体力学的基本知识 流体流动现象大量存在于自然界及多种工程领域中，所有这些工程都受质量守恒、动量守恒和能量守恒等基本物理定律的支配。这章将首先介绍流体动力学的发展和流体力学中几个重要守恒定律及其数学表达式，最后介绍几种常用的商业软件。 2.1计算流体力学简介 2.1.1计算流体力学的发展 流体力学的基本方程组非常复杂，在考虑粘性作用时更是如此，如果不靠计算机，就只能对比较简单的情形或简化后的欧拉方程或N-S方程进行计算。20 世纪30～40 年代，对于复杂而又特别重要的流体力学问题，曾组织过人力用几个月甚至几年的时间做数值计算，比如圆锥做超声速飞行时周围的无粘流场就从1943 年一直算到1947 年。 数学的发展，计算机的不断进步，以及流体力学各种计算方法的发明，使许多原来无法用理论分析求解的复杂流体力学问题有了求得数值解的可能性，这又促进了流体力学计算方法的发展，并形成了"计算流体力学" 。 从20 世纪60 年代起，在飞行器和其他涉及流体运动的课题中，经常采用电子计算机做数值模拟，这可以和物理实验相辅相成。数值模拟和实验模拟相互配合，使科学技术的研究和工程设计的速度加快，并节省开支。数值计算方法最近发展很快，其重要性与日俱增。 自然界存在着大量复杂的流动现象，随着人类认识的深入，人们开始利用流动规律来改造自然界。最典型的例子是人类利用空气对运动中的机翼产生升力的机理发明了飞机。航空技术的发展强烈推动了流体力学的迅速发展。 流体运动的规律由一组控制方程描述。计算机没有发明前，流体力学家们在对方程经过大量简化后能够得到一些线形问题解读解。但实际的流动问题大都是复杂的强非线形问题，无法求得精确的解读解。计算机的出现以及计算技术的迅速发展使人们直接求解控制方程组的梦想逐步得到实现，从而催生了计算流体力
计算流体力学课后题作业
课后习题 第一章 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基于有限体积法，适用于不可压流体和可压流的计算、热力学的计算及非牛顿流的计算。它具有前处理器、求解器、后处理器三大模块，以良好的可视化用户界面把建模、求解及后处理与全部的物理模型和算法结合在一个软件包中。
上计算水力学课的心得
上计算水力学课的心得 水利水电学院水力学及河流动力学 胥慧1030201016 摘要：首先通过计算水力学这门课程的学习，联想到不规则的平面图形面积的求解；还简要说明了从中学到的内容，着重说明了离散的有关问题；最后阐述了自己对这门课程的几点意见。 关键词：面积，区域离散，控制方程离散，意见 1、不规则图形面积求解 上计算水力学这门课程时，我突然想起小时候学过对于一个边界形状不规则的平面图形面积问题的求解方法。当时是先把那个不规则的平面图形誊画在一个透明的玻璃板上，再把一张事先做好的1cm×1cm方格纸铺在玻璃板下边，先记录一下不规则图形里显示完整的小方格数目，对于不完整的小方格，正好满半个格算的两个算一个格，大于半个格计一个格，不满半个格的舍去，这样相加在一起就是这个不规则的几何图形的近似面积。同样的办法，再分别用0.5cm×0.5cm 的方格纸和0.1cm×0.1cm的方格纸对不规则图形面积进行计算。结果不言而喻，必然是用0.1cm×0.1cm的方格纸得到的近似解更接近真实解。通过缩短方格纸的边长，来实现接近真实解的方法。用类比的方法学习了计算水力学这门课。2、学到的内容 在以前的学习中我了解到，描述流体流动及传热等物理问题的基本方程为偏微分方程，想要得它们的解析解或者近似解析解，在绝大多数情况下都是非常困难的，甚至是不可能的，就拿我们熟知的Navier-Stokes方程来说，现在能得到的解析的特解也就70个左右。通过学习计算水力学这么课程，我知道对这些问题进行研究，可以借助于现在已经相当成熟的代数方程组求解方法，对于这种方法简单来说就是将连续的偏微分方程组及其定解条件按照某种方法遵循特定的规则在计算区域的离散网格上转化为代数方程组，以得到连续系统的离散数值逼
CFD计算流体动力学入门教程选择
非流体、热动专业CFD新手入门 首先掌握流体力学基本原理，丁祖荣主编的流体力学这本教材，仔细看两天，这样就会知道gambit中为什么会有边界层设置，边界层厚度如何设置；雷诺系数如何确定来判断层流与湍流；马赫数如何确定来判断流体是可压还是不可压，这样就能解决Fluent，是基于压力还是基于密度求解。能够对实际中一些看似简单的流体现象有深刻的认识，能够准确判断是定常流还是非定常流。 CFD网格划分 网格划分对于初学者所接触案例，其实非常简单。但实际工程中，大项目，特别涉及到整套工程，如环保，飞机，网格质量与数量都要求非常高，往往服务器类的PC才能解决问题，所谓的内存128G，CPU四核主频3.0以上。初学者，简单的管道，一般的机器还是没问题。有机械三维软件基础的，对于gambit建模就非常容易了。往往大项目，复杂的结构gambit 建模显得力不从心，所以对于流体工作者来说，学习三维软件对于建模有莫大的帮助，如Proe。 1.1Gambit介绍 网格的划分使用Gambit软件，首先要启动Gambit，在Dos下输入Gambit
药物动力学计算题
1.计算题：一个病人用一种新药，以2mg/h的速度滴注，6小时即终止滴注，问终止后2小时体血药浓度是多少？（已知k=0.01h－1，V＝10L） 2.计算题：已知某单室模型药物，单次口服剂量0.25g，F=1，K=0.07h－1，AUC=700μg/ml·h，求表观分布容积、清除率、生物半衰期（假定以一级过程消除）。 3.某药静注剂量0.5g，4小时测得血药浓度为 4.532μg/ml，12小时测得血药浓度为2.266μg/ml，求表观分布容积Vd为多少？ 4.某人静注某药，静注2h、6h血药浓度分别为1.2μg/ml和0.3μg/ml（一级动力学），求该药消除速度常数？如果该药最小有效剂量为0.2μg/ml，问第二次静注时间最好不迟于第一次给药后几小时？ 5.病人静注复方银花注射剂2m/ml后，立即测定血药浓度为1.2μg/ml，3h为0.3μg/ml，该药在体呈单室一级速度模型，试求t1/2。 6.某病人一次用四环素100mg，血药初浓度为10μg/ml，4h后为 7.5μg/ml，试求t1/2。 7.静脉快速注射某药100mg，其血药浓度－时间曲线方程为：C=7.14e-0.173t，其中浓度C的单位是mg/L，时间t的单位是h。请计算：（1）分布容积；（2）消除半衰期；（3）AUC。
8.计算题：某药物具有单室模型特征，体药物按一级速度过程清除。其生物半衰期为2h，表观分布容积为20L。现以静脉注射给药，每4小时一次，每次剂量为500mg。 求：该药的蓄积因子 第2次静脉注射后第3小时时的血药浓度 稳态最大血药浓度 稳态最小血药浓度 9.给病人一次快速静注四环素100mg，立即测得血清药物浓度为10μg/ml，4小时后血清浓度为7.5μg/ml。求四环素的表观分布体积以及这个病人的四环素半衰期（假定以一级速度过程消除）。 10.计算题：病人体重60kg，静脉注射某抗菌素剂量600mg，血药浓度－时间曲线方程为：C=61.82e-0.5262t，其中的浓度单位是μg/ml，t的单位是h，试求病人体的初始血药浓度、表观分布容积、生物半衰期和血药浓度－时间曲线下面积。 11.计算题：已知某药物具有单室模型特征，体药物按一级速度方程清除，其t1/2=3h，V=40L，若每6h静脉注射1次，每次剂量为200mg，达稳态血药浓度。求：该药的（1）ss C max （2）ss C m in （3）ss C （4）第2次给药后第1小时的血药浓度
流体力学的计算公式
流体力学的计算公式 众所周知，锅炉之类的热力设备可以采用水封管作为限压装置。设备工作时，水封管内的水将设备与大气隔绝，使设备保持一定的压力。水封管内水的多少，即水封管内的实际水位的高低与设备汽压大小有关。汽压大，水位就高；汽压小，水位就低。当设备压力升至额定蒸汽压力时，水封管内的水位应该达到最高水位；而在设备压力稍微大于额定蒸汽压力时，水封管内的水应能立即被冲除掉，使设备内的蒸汽能够迅速得到排放。但是，水封管的内径必须足够大；以保证它的排汽能力大于设备的最大蒸发量，从而防止设备发生超压事故。《蒸汽锅炉安全技术监察规程》第132条规定：“对于额定蒸汽压力小于等于0.1MPa的锅炉可以采用静重式安全阀或水封式安全装置。水封装置的水封管内径不应小于25mm”。《小型和常压热水锅炉安全监察规定》第二十八条要求：“水封管的内径应当根据锅炉的额定容量和压力确定，且内径不得小于25毫米”。 但是，水封管的内径究竟应该根据什么样的公式去计算？这里没有具体加以规定。水封管的内径为何不能小于25mm？这里没有加以说明。 1 确定水封管内径的一种错误观点 有人说，可以按安全阀的排放量计算公式去间接确定水封管的内径。笔者认为这是不对的。理由如下： 其一，按流体力学，安全阀的局部阻碍因为其进出口呈直角型式，边壁是突变的，主流与边壁之间形成大尺度旋涡，蒸汽排放时能量损失很大；而水封管的边壁是渐变的，又不出现减速增压现象的部位，故蒸汽排放时的能量损失很小。其二，安全阀的出口之外肯定还有排汽管，它的排放量计算公式应该考虑这一额外的因素；而水封管却只相当于它的排汽管。其三，安全阀在排汽时，汽流还自始至终受到弹簧或者重锤施加的反方向作用力，能量损失非常之大；而水封管在排汽时却不会受到这样的反作用力，无此项能量损失。显然，安全阀排汽与水封管排汽区别太大，不能按安全阀排放量计算公式去确定水封管的内径。 2 用流体力学理论建立水封管内径计算方法 为了正确确定水封管的应有内径，我们有必要按流体力学的理论来分析一下水封管的流动阻力和能量损失。蒸汽流过水封管时的阻力损失ΔP等于水封管进口压力P1与出口压力P2之差。管道的阻力损失 ΔP＝Pa（1） 式中λ—管道摩擦系数； d —管道内径，m； ρ—流体的密度，k g／m3； u —流体在管内的流速，m／s； L—平直管段的总长度，m； Σξ—管道中各种局部阻力之和。 管道的摩擦系数λ值取决于流体流动的雷诺数Re 和管壁的相对粗糙度。管壁的相对粗糙度等于管壁的绝对粗糙度（即管壁内凸起高度）K与管内径d之比。雷诺数Re =(2) 式中ν－流体运动粘滞系数，m2／s； 其他符号说明同上。 由于水封管进出口压差通常大于40000Pa，排汽时的流速都大于100m／s ，雷诺数Re 都大于80000000，所以蒸汽的流动状态为紊流流动，而且位于阻力平方区。在此区域内，摩擦阻力系数λ值仅取决于管壁的相对粗糙度，可以用尼古拉兹公式进行计算：λ＝（3）
全国高考动力学计算专题
2010―― 2015年物理高考动力学计算题汇编 1. ( 2010全国大纲I 卷)汽车由静止开始在平直的公路上行驶, 间变化的图线如右图所示。 (1) 画出汽车在 0~60s 内的v-t 图线； (2) 求在这60s 内汽车行驶的路程。 1?解(I)设t=10, 40, 60 s 时刻的速度分别为 V i , V 2 , V 3。 由图知0~10 s 内汽车以加速度 2m|_s 工匀加速行驶，由运动学公式得 V 1 =2x10=20 m / s ① 由图知10~40 S 内汽车匀速行驶.冈此 v 2 =20 m/s ② I ~I _2 由图知40~60 s 内汽车以加速度1m 生 匀减速行 驶.由运动学公式得 根据①②③式，可画出汽车在 0~60 s 内的V -1图线，如 右图所示。 ⑵由右图可知，在这 60 s 内汽车行驶的路程为 30 60 s 20=900m ④ 2 2. ( 2010课标1卷) 短跑名将博尔特在北京奥运会上创造了 100m 和200m 短跑项目的新世 界纪录，他的成绩分别是 9. 69 s 和19 . 30 s 。假定他在100 m 比赛时从发令到起跑的反 应时间是0. 15 S ，起跑后做匀加速运动，达到最大速率后做匀速运动。 200 m 比赛时，反 应时间及起跑后加速阶段的加速度和加速时间与 100 m 比赛时相同，但由于弯道和体力等因 素的影响，以后的平均速率只有跑 100 m 时最大速率的96%。求：(结果保留两位小数) (1) 加速所用时间和达到的最大速率： (2) 起跑后做匀加速运动的加速度。 2.解：(1)设加速所用时间为t (以s 为单位)，迅速运动的速度为 v (以m/s 为单位)， 0 ~60s 内汽车的加速度随时 v 3 = 20 -1 20 =0 ③
流体力学计算题..
水 水银 题1图 1 2 3 题型一：曲面上静水总压力的计算问题（注：千万注意方向，绘出压力体） 1、AB 曲面为一圆柱形的四分之一，半径R=0.2m ，宽度（垂直纸面）B=0.8m ，水深H=1.2m ，液体密度3 /850m kg =ρ，AB 曲面左侧受到液体压力。求作用在AB 曲面上的水平分力和铅直分力。（10分） 解：（1）水平分力： RB R H g A h P z c x ?-==)2 (ργ…….(3分) N 1.14668.02.0)2 2 .02.1(8.9850=??- ??=，方向向右（2分） 。 （2）铅直分力：绘如图所示的压力体，则 B R R R H g V P z ??? ? ????+-==4)(2πργ……….（3分） 1.1542 8.042.014.32.0)2.02.1(8.98502=???? ? ?????+?-??=，方向向下（2分）。 2．有一圆滚门，长度l=10m ，直径D=4.2m ，上游水深H1=4.2m ，下游水深H2=2.1m ，求作用于圆滚门上的水平和铅直分压力。
解题思路：（1）水平分力： l H H p p p x )(2 1 222121-= -=γ 方向水平向右。 （2）作压力体，如图，则 l D Al V p z 4 432 πγγγ? === 方向垂直向上。 3.如图示，一半球形闸门，已知球门的半径m R 1= ，上下游水位差m H 1= ，试求闸门受到的水平分力和竖直分力的 大小和方向。 解： (1)水平分力： ()2R R H A h P c πγγ?+===左,2R R A h P c πγγ?=' =右 右左P P P x -= kN R H 79.30114.31807.92=???=?=πγ， 方向水平向右。 （２）垂直分力： V P z γ=，由于左、右两侧液体对曲面所形成的压力体均为半球面，且两侧方向相反，因而垂直方向总的压力为0。 4、密闭盛水容器，已知h 1=60cm,h 2=100cm ，水银测压计读值cm h 25=?。试求半径R=0.5m 的半球盖AB 所受总压力的水平分力和铅垂分力。
计算流体力学大作业报告(翼型空气动力分析)
课程综合作业课程名称：计算流体力学 专业班级：研究方向： 学生姓名：学号： 完成日期：
计算流体力学课程综合报告 1.简介 计算流体动力学（Computational Fluid Dynamics，简称CFD）是通过计算机数值计算和图像显示，对包含有流体流动和热传导等相关物理现象的系统所做的分析。其基本思想为:把原来在时间域及空间域上连续的物理量的场，如速度场和压力场，用一系列有限个离散点上的变量值的集合来代替，通过一定的原则和方式建立起关于这些离散点上场变量之间关系的代数方程组，然后求解代数方程组获得场变量的近似值。 CFD可以看作是在流动基本方程（质量守恒方程、动量守恒方程、能量守恒方程）控制下对流动的数值模拟。通过这种数值模拟，我们可以得到极其复杂问题的流场内各个位置上的基本物理量（速度、压力、温度、浓度等）的分布，以及这些物理量随时间的变化情况，确定旋涡分布特性、空化特性及脱流区等。还可据此算出相关的其他物理星，如旋转式流体机械的转矩、水力损失和效率等。此外，与CAD联合，还可进行结构优化设计等。 2.计算流体动学的特点: ①流动问题的控制方程一般是非线性的，自变量多，计算域的几何形状和边界条件复杂，很难求得解析解，而用CFD方法则有可能找出满足工程需要的数值解。 ②可利用计算机进行各种数值试验，例如，选择不同流动参数进行物理方程中各项有效性和敏感性试验，从而进行方案比较。 ③它不受物理模型和实验模型的限制，省钱省时，有较多的灵活性，能给出详细和完整的资料，很容易模拟特殊尺寸、高温、有毒、易燃等真实条件和实验中只能接近而无法达到的理想条件。 ④数值解法是一种离散近似的计算方法，依赖于物理上合理、数学上适用、适合于在计算机上进行计算的离散的有限数学模型，且最终结果不能提供任何形式的解析表达式，只是有限个离散点上的数值解，并有一定的计算误差。 ⑤它不像物理模型实验一开始就能给出流动现象并定性地描述，往往需要由原体观测或物理模型试验提供某些流动参数，并需要对建立的数学模型进行验证。