N-S方程程序.docVIP

N-S方程的物理和几何意义及其应用 物理意义：运动的流体能量守恒的表现，即流体微团的动量变化率等于作用在流体微团上的合力。 油的泄漏流量的计算： 如右图所示 求一下柱塞缝间的泄漏： 用N-S方程分析： 假设活塞是不动的，这是个定常流动问题，所以 如不计质量力作用 根据N-S方程式得到-+v 0 速度仅仅是y的函数，压强仅仅是x的函数，所以 积分的u dp/ dx*y2+c1y+c2 由边界条件u/y o o, u/y h 0 得到C1 - C2 0 所以得到速度分布为：u （h-y）*y 沿流向长度为L的距离上两端压差p p1-p2 所以u （h-y）*y 计算条件：必须知道油缸的相对压力P，以及油的动力黏度u，套筒的长度L,柱塞的直径d,套筒与柱塞的间隙h 计算结果： 由dq udV uπd*dy 两边积分得Q πd/ 12μ h3 如果知道各数据就可以计算出流量 黏性流体动量平衡方程 纳维-斯托克斯方程（Navier-Stokes? equations） 1.动量平衡的定义 流体在流动过程中遵守能量守恒定律，称为能量平衡 根据牛顿第二定律：? 作用力的合力 单位时间内动量的变化量 作用力形式 动量形式 [动量传入量] [动量传出量] +[系统作用力的总和] [动量蓄积量] [动量传入量] [动量传出量] + [系统作用力的总和] 0 稳定流动系统： 不稳定流动系统： 动量收支差量 动量收支差量 ⒉ 动量传递方式 1 黏性动量传输 黏性动量通量 2 对流动量传输 对流动量传输 ⒊ 作用力的形式 体积力 表面力 压力 重力 作用力 ⒋ 动量平衡方程的推导 元体分析法 牛顿第二定律分析法 建立方法 建立依据 Y 在直角坐标系中由于有三个方向的分速度，所以共有九个动量通量。 以vx为准：动量通量 动量通量收支差量 ⑴ 对流动量收支差量 x方向的速度、x方向的动量通量 对流动量收支差量为 同理，以vx为准，y方向、z方向的对流动量收支差量： 以vx为准，元体对流动量收支差量为 同理，以vy、vz为准，元体对流动量收支差量为 vx vy、vz ⑵ 黏性动量收支差量 黏性动量通量同样由九个分量组成 以vx为准，C、D面上的黏性动量通量为 黏性动量通量收支差量 黏性动量收支差量 同理，vx在y、x方向的黏性动量收支差量分别为 以vx为准，元体黏性动量收支差量为 同理，以vy、vz为准的黏性动量收支差量为x y、z 牛顿流体： ⑶ 作用力的总和 x方向：PA x方向合压力为 x方向的总压力为 同理，y、z方向的总压力为x y、z PA A B o x y z 重力 ⑷ 动量蓄积量 z方向 x方向 y方向 单位时间内元体动量的变化量 [动量传入量] [动量传出量] +[系统作用力的总和] [动量蓄积量] ⒌ 动量平衡方程式 将以上式子代入下式，整理得：N-S方程 简化： ，连续性方程 ⑵ ，牛顿黏性定律 ⑴ 动量收支差量 黏性力引起 压力 体积力 积累 动量收支差量 ⒍ 动量平衡方程的讨论 对流动量 动量蓄积量 黏性动量 压力 重力 （1）方程的物理意义：运动的流体能量守恒的表现 全微分 惯性力 黏性力 压力 重力 流体在运动中以作用力及动量形式表现能量平衡 关系是统一的 ⑵ 适用条件 黏性流体、不稳定流动、不可压缩流体 元体范围内 、层流流动 理想流体： 没有黏性的流体 简化： ②　稳定流动， ③　单位质量流体 ①　 时，N-S方程简化为欧拉方程 理想流体、稳定流动、不可压缩流体 元体范围内 适用条件 对非高速流，气液界面上， 液相速度梯度为零。 流体具有粘性，在与壁面 接触处流体速度为零。 液液界面两侧的速度或切应力相等。 液体－液体边界： 液体－气体边界： 固壁－流体边界： 边界条件 根据问题特点对一般形式的运动方程进行简化，获得针对具体问题的微分方程或方程组。 提出相关的初始条件和边界条件。 　　初始条件：非稳态问题 流动微分方程的应用求解步骤 圆管中的层流 两平行板间的层流 设无限大水平平行平板均固定不动，间距为b。y轴垂直板面， 不可压缩牛顿流体在x方向的恒定压强梯度（?常数）作用下作充分发展层流， 可简化为x y平面上的流动。