N-S方程推导.docx

模拟flunte软件模拟T型血管内的血流动力学 模型的建立 2.1 N-S方程 N-S方程的最初形式 作用在单元体上的力 作用力有两类，即质量力和表面力。 1.1 质量力 质量力是作用在每一个流体质点上，大小与流体的质量成正比。工程流体力学中，会遇到两种质量力：重力和惯性力。惯性力是一个很特殊的称谓，原来中学教程中认为惯性力并不是力，但是实际上，在出现加速度的时候，惯性力的作用同普通力完全一样的，只是惯性力会随着加速度的消失而消失。如果认为惯性力是一种力，那么牛顿第二定律（1）也可以认为是力的平衡。式的右端就是惯性力，左端就是其他的常规力。其实观察一下重力，，同惯性力的本质上是一致的，本身就是重力加速度。 假设单位质量流体上的质量力在各个坐标轴的分量分别为X，Y，Z。图1流体单元体的质量为：。则作用在流体单元体上的质量力在坐标轴的分量分别为：、、。 1.2 表面力 作用在隔离流体（也就是所取的研究流体单元）的表面，和作用的面积成正比的力。分为垂直于作用面的压力和眼沿作用面方向的切力。表面力可以使作用于流体界面的压力、切力。单位作用面的压应力、切应力即为图1中的、（第一个下标表示作用面的法线方向，第二个下标表示力的方向）。 以X方向为例，流体单元受到的力： (2) 即: (3) y，z方向同理可获得。 （4） （5） 单元体的加速度和重量 加速度和质量的乘积（（1）式右侧）在三个方向上的分量分别为： （6） （7） （8） 将（3） （6）式代入（1）式，X方向有： （9） 即： （10） 同样： (11) (12) 应力形式化简 1、切应力与应变的关系 （13） （14) (15) 法向应力与应变的关系 （16） （17） （18） 将（13）、（16）代入（10）， （19） （20） （21） 即： （22） 同理： (23) （24） (25) 矢量形式： (26) 不可压缩流体的N-S方程 连续性方程的基本推导原理就是，单元体流出、流入质量差等于该时间段内质量的变化。原理是很简单的。没有流入流出质量就不会变化，流入流出有了差值，说明单元体的质量变化了。 仍以x方向为例。 左侧质量流速（一般的流速是体积流速，m/s,为了推到质量的变化需要引入质量流速，质量流速的定义就是单位时间内通过单位横截面的流体质量）为，质量流速是位置的函数，因此在右侧面流出的质量流速为。时间段内流出、流入单元体的质量差为： （27） 同理，改时间段dt内y方向，z方向的流出流入质量差为： （28） （29） 因此，时间段内单元体六面流出、流入的质量差为： （30） 该时间段内单元体质量的变化体现在密度随时间的变化上，开始时间密度为，时间末密度为。质量的变化为： （31） 根据质量守恒，（30）式等于（31）式，即 （32） 化简: (33) (34) 如果不可压缩流体，密度=，密度项可以提取出来，散度为： （34） 将（34)式代入（26）式，不可压缩流体的N-S方程为： (35) 加速度项的处理 流动中，不仅不同位置的点具有不同的速度，就是在同一点，不同时刻速度也可能不同。速度即时位置的函数也是时间的函数。因此，加速度有两部分组成，迁移加速度和当地加速度。 以x方向为例，加速度的表达式为： (36) 式中，单位