纳威―斯托克斯方程推导专题教学改革探索.docVIP

纳威―斯托克斯方程推导专题教学改革探索 　　【摘要】水力学是一门技术基础科学，它是力学的一个分支。N-S（Navier-Stokes）方程是水力学课程的重要内容，是联系本科生水力学与研究生高等流体力学课程的关键环节。但现实授课过程中经常听到学生反映很难理解N-S方程，并且觉得没有应用价值，原因是本科教学中没有用其做过习题。因此，为了更好地让学生理解N-S方程的原理，掌握方程的推导过程，我们对此开展专题性教学改革探索，进而推进水力学重、难点的教学改革与创新，使学生更愿意学，更容易懂，更能真正理解和应用。 　　【关键词】水力学 N-S方程 研究型本科 专题教学 教学探索 　　【中图分类号】G642.0 【文献标识码】A 【文章编号】2095-3089（2014）11 -0060-02 　　水力学是高等工科院校水利、土木、环境等诸多专业都需要掌握的一门重要技术基础课，是学生专业能力形成与未来职业发展必不可少的依托[1]。水力学的任务主要是研究液体（主要是水）的平衡和机械运动规律及其工程应用[2]，是一门既有较强理论性又有较强工程实际意义的课程，具有理论推导复杂，概念多且易混淆特点[3]。另外，水力学具有对学生的高等数学、大学物理以及综合分析问题的能力要求较高的特点[4]。作为专业技术基础，它不仅为后续专业课学习提供理论基础，而且注重理论联系实际[5]。比如用于描述粘性不可压缩流体动量守恒的纳维-斯托克斯方程（N-S方程），是不可压缩流体动力学最有用的一组方程之一，也是学习高等流体力学、计算流体力学等课程的重要基础[6]。N-S方程可用于描述与分析大量对学术和经济有用现象的物理过程，包括模拟管道中的水流、翼型周围的气流、污染效应的分析、电站的设计、天气、洋流、星系中恒星的运动、飞行器和车辆的设计、血液循环的研究等[7]。 　　因此，考虑到N-S方程在上述专业领域本科教学过程中的重要地位，本文从专题教学的角度积极探索本科生容易接受的N-S方程推导方法，把握水力学研究型本科教学内涵，着力提升学生的理论基础和实践能力，为学生创新精神和研究能力的培养提供可靠支撑。 　　一、液体运动微分方程的推导过程 　　液体运动微分方程（又称N-S方程）是水动力学的基础。目前，在《水力学》教材中，推导N-S方程的方法很多，本文按照由平衡到运动、由理想到实际的专题教学思路推导N-S方程，帮助学生牢固掌握方程的意义。 　　（一）液体平衡微分方程式 　　在静止液体内，任取一点M，该点压强为p（x，y，z），以M为中心做微元直角六面体，正交的三个边分别与坐标轴平行，长度为dx，dy，dz。微元六面体静止，各方向的作用力平衡（图1），以y方向为例，对其进行受力分析，对表面力取泰勒（Taylor）级数展开式的前两项，列y方向表面力与质量力的平衡方程 ， 　　上式即为液体平衡微分方程，是瑞士数学家和力学家欧拉于1755年导出的，又称为欧拉平衡微分方程。方程表明，在静止液体中各点单位质量流体所受表面力和质量力相平衡。 　　图1 静止液体平衡微元六面体y方向的应力图示 　　（二）理想液体运动微分方程 　　在运动的理想流体中，同样取微小平行六面体如图1所示，分析该微小六面体y方向的受力和运动情况。根据牛顿第二定律，建立y方向的动力学方程。 　　上式即理想流体运动微分方程式，又称欧拉运动微分方程式。其物理意义：作用在单位质量流体上的质量力与表面力之代数和等于其加速度。适用于无黏性流体，可压缩流体及不可压缩流体，恒定流及非恒定流。 　　（三）实际液体运动微分方程 　　在有粘性的实际液流中，不仅有压应力，还有切应力存在，其表面力的方向不与作用面相垂直，可分解成相互正交的一个法向应力和两个切向应力（图2）。以y方向为例，进行受力分析，质量力可表示为，表面力中压应力和切应力大小如图2所示。 　　图2 实际液流微元六面体各表面的应力图示 　　根据牛顿第二定律，同理可得x和z方向的分量式，将其整理后一并写为： 　　式（6）即为以应力形式表示的运动微分方程，简称应力微分方程。 　　根据推导的粘性液体切应力与压应力特性关系式，将应力与变形率的关系式及不可压缩液体连续性方程代入方程式（6），并将加速度项展开，整理得沿x，y，z方向的方程： 　　式（7-9）即为不可压缩液体运动微分方程。它是由纳维（Navier）于1821年首先提出，后由斯托克斯（Stokes）于1845年完善而成，故称为纳维-斯托克斯方程，简称N-S方程。 　　二、分析与讨论 　　在推导N-S方程的过程中，离不开与液体平衡微分方程和理想液体运动微分方程（又称欧拉方程）进行对比分析。在静止液体中，由于液体质点间没有相对运动，液体的粘滞性不起作用，质点只承受压应力即静水压强；在运动的理想