首页/文章/ 详情

Navier-Stokes方程都看不懂，还怎么做模拟？

1年前浏览2423

导读：本文转自热流工程坊原创文章，已授权，旨在传播知识，希望对学习者有所帮助。以下是正文：

曾几何时，我们简要分析了数学中的“▽”到底是干嘛用的，并最后通过表格对梯度、散度和旋度进行了对比。

没关系，我直接在下面贴出来供大家复习。

有了上面的基础，今天就以不可压缩流体为例来看看，流体力学里面的连续性方程和Navier-Stokes方程到底是啥意思？

连续性方程

Navier-Stokes方程

一、连续性方程

简单来说，连续性方程就是质量守恒定律在流体力学中的具体表述形式。

如上图所示，假设流体流动速度为v且与平面AB垂直，考虑图中的平面AB，则在dt的时间内，通过平面AB的质量流量dΦ为

其中的ρ表示流体密度。

进一步，如果速度v和平面AB的方向不垂直，而是存在一个夹角θ，则在dt的时间内，通过这个平面的质量流量dΦ就成了

对于一个闭合曲面而言，我们只需对每一个平面微元的质量流量进行求和，即可得到在dt时间内流体通过该闭合曲面的质量流量Φ为

也就是说：

因为流体流动，流出闭合曲面的流体质量为Φ。反之，流入闭合曲面的流体质量为-Φ。

此外，闭合曲面内流体的质量为

在dt的时间内，闭合曲面内流体的质量增加量dm为

由于流体质量守恒，流体的质量增加量dm就应当等于流入闭合曲面的流体质量-Φ。

变形整理一下就是

根据散度的定义，有

因此

因为流体密度ρ除了是时间t的函数以外，还是位置x,y,z的函数，所以这里用偏导数代替导数，得到

也就是说，对于任意区域来说，流体质量的增加量就等于流入该区域的流体质量。

换算到单位体积（其实就是去掉积分符号），自然就得到了

这就是我们平时见到的连续性方程。

进一步，如果流体不可压缩，那流体密度ρ就成了常数且不随时间发生变化，因此连续性方程就简化为了

二、Navier-Stokes方程

类比连续性方程，我们自然也可以知道，对于任意区域来说，流体动量的增加量就等于流入该区域的流体动量。

质量表达式为

动量表达式为

也就是说，动量与质量相比，也就是多乘了个v而已。

所以动量守恒方程就变成了

这比Navier-Stokes方程方程简单太多了，你这是错的吧！

的确，以上情况只在流体完全不受力的情况下成立。

牛顿爵士告诉我们，动量为p的质点，在外力F的作用下，其动量随时间的变化率等于该质点所受的外力。

牛顿

也就是说，除了流入该区域的流体动量，流体区域内所受的力也会造成流体动量的变化。

其中一定存在的力有两个，一个是压力，而另一个就是粘性力。压力主要用于平衡流体之间的受力，而粘性力则是由于流体的流动而引起的。对于不可压缩流体而言，它们的大小就等于

由于推导实在是太过复杂，这里就不详细推导了。我们直接把前人的结果拿来用就好了。想具体了解的同学可以参考周光坰流体力学上册的1.6和4.5节。

至于其他力，比如说重力，电磁力等等，既然没有表达式，不妨直接把它们写到一起，记做

力的存在效果和流入该区域的动量是等价的，因此N-S方程最终就成了

也就是说，对于任意区域来说，流体动量的增加量就等于流入该区域的流体动量加上该区域所受的力。

换算到单位体积（其实就是去掉积分符号），自然就得到了

这就是我们平时见到的Navier-Stokes方程。

三、万物流量皆相等（我瞎编的）

事实上，以上推导逻辑可以应用到任何关于流动的场景。

对于任意变量Φ，最终的方程基本上都可以写成：

物理量的增加量=物理量的流入量+引起物理量增加的其他原因。

变成公式（其实就是把连续性方程的密度ρ换成变量Φ）就是

里面的U表示unkown，意思是我也不知道它是个啥。

是不是还挺有意思的，大家快都来学习流体力学吧！比如说，从关注下面的二维码开始！

流体力学基础21讲

（完）

来源：仿真秀App

电磁力物流曲面

著作权归作者所有，欢迎分享，未经许可，不得转载

首次发布时间：2023-05-09

最近编辑：1年前

技术圈粉知识付费学习强国

获赞 10051粉丝 21523文章 3526课程 218

点赞

还没有评论

相关推荐

Creo/proe-beats耳机-高级曲面造型-一加一教育 SolidWorks结构分析入门基础课—7天快速掌握结构分析基础知识、结构分析一般过程 Maxwell高级进阶案例7-直线感应电机设计和电磁仿真 2、Creo曲面建模教程第二例 CST Studio Suite建模和仿真20讲：获得滤波器、连接器、全腔仿真和双工器等仿真关键技能

最新文章

螺旋折流板式换热器

FLAC3D6.0 | 如何导入壳型结构单元

中国版的TechSoft3D，国产CAD/CAE开发组件提供商有哪些？

OptFuture | CFD气动分布载荷在结构仿真中的加载方式介绍

不同骨强度下全膝置换过程中发生股骨前皮质切迹的三维有限元分析

热门文章

仿真笔记——ANSYS APDL命令汇总（收藏备用）

STAR CCM 案例｜电池包散热

CloudCompare---一款开源点云处理软件

仿真天地 • 汽车｜AVL的动力电池探索创新之路

基于WEB的仿真平台建设方案-上海幂知科技有限公司

其他人都在看

盘点·近十年来国外各公司推出的碳纤维产品

几种常见的热仿真软件

Abaqus分析常见问题及解决方法（2）：零主元和过约束

ABAQUS中Cohesive粘聚力模型的2种定义方式（附案例操作步骤）

仿真工作者必须知道的15款开源软件！

VIP会员学习福利任务兑换礼品帮助与反馈

返回顶部