在2014年,我与好几个学友谈论过流体力学中的NS方程。眼下有空,也就把自身近几年的思考简述如下。

就NS方程的推导及其所反映的客观现象而言,NS方程是对流体微元在瞬时意义上变形运动的描述。在流体力学本构方程中的压力是天外来客,在力学本质上,压力的空间梯度是微元体惯性力的表征。然而,在本构方程中,它是被置于变形应力的地位。从而,在哲学上,惯性力(牛顿意义上的)是没有独立地位的。这是因为,如果研究稳态流动,则压力是产生变形力(流速空间梯度)的唯一原因(反之也然)。

微元体本身的绝对惯性速度(整体平移)是对流体应变没有贡献的,这个假疵也就由压力来作为补赏项以各向同性应力的形式而出现了。就最原始的流体力学而言,压力是由流体的物性及整体平移速度所确定的。而对于漩涡的研究也表明,漩涡的转速(或是旋度)也是由局部整体的转动速度所决定的。而这类局部整体转动在局部效果上等价于局部整体平移,从而也反映为压力的局部特征。由于在这个局部意义上的整体是有强烈的局部尺度依赖性质的,所以对流体动力学学而言,尺度性是非常明显的。

在大尺度上,如果整体的惯性平移被压力所包容(速度梯度已消除的空间不变速度部分),那么在此尺度下的局部整体转动(小涡)就被忽略了。而在方程中的忽略只是表明方程本身的局限性,并不能被看成是正确的(物理上真实的)。因而,总可以通过把尺度变小来把适当尺度的小涡考察进来。得到的好处是应力(或压力)更接近于局部实测量,但是,由此所付出的代价是在NS方程中不断的补充新的应力(或压力)项。方程趋向于无限的复杂化。最后干脆就是把流体动力学运动看成是一个统计学现象,从而,在哲学上是对其原出发点NS方程的否定。

对于固体变形,整体的平移、转动对应变没有贡献,从而也对应力没有贡献。这基本上是把运动中的固体等价于静止的固体,哲学上认为:整体的平移、转动对局部变形没有贡献。这基本上是可以接受的(近似意义上)。但是,对流体这是明显的无法接受的。然而,流体力学以引入静水压力的方式,在速度梯度为零时保持了应力的存在。这是区别于固体变形力学的。因而,我们可以认为,静水压力项在动力学方程上的意义是把整体的平移、转动的物理效应以间接的方式给强行纳入进来。

事实上,理性力学家在1965年前后对于这个问题的看法是:在整个连续介质力学中,以本构方程来区别流体物质和固体物质是形式上协调的,但是在物质定义上是不协调的。间接的认为,流体的本构方程是非理性的。其中,Truesdell则明确的在其专著中表明静水压力项的存在是力学理论无法讲清的,只能作为一个经验项被动的接受下来。

在回头看NS方程的推导,压力项是以隐含的方式被应力项所隐盖的,从而其有效性是在回避压力的独立物理意义而取得的,在实质上把流体微元等价于高速运动中的可变形固体微元。在这样一个哲学意义下也就明确了:连续性方程是一个事后的补充办法。

因而,我们可以看到这类做法是比较有成绩的:把连续性方程看成是一个独立方程(再本质上是补充一个关于局部整体惯性速度的约束方程),而NS作为变形力学方程。

事实上,如果客观的看待一百年来的流体力学研究,那么Truesdell对涡量的研究是要把局部的整体转动从NS方程中剥离出来,认为它满足其它的动力学方程(类似于角动量守恒类),而欧拉的流体力学方程基本上就是立足于角动量守恒的考虑。当然,NS方程也可以接近圆满的表达涡动,但是付出的代价是变尺度的把涡搞成尺度上的分形属性(嵌套结构)。在多个尺度上重复使用NS方程在力学理论上看是一种变通办法,也是理论局限性的特征体现。

另一类比较有效的办法是先引入一个流矢量函数,既由它来定义流动的整体属性,也由它来定义流体微元的变形属性,这样,与NS方程不同的是,对于流函数,有独立的物理方程作为控制。目前,此类理论还在发展之中。可以把它概括为:修订NS方程,至少引入一组关于流动整体属性的动力学方程。

虽然有关的研究工作没有明确的说明:光是NS方程不足于描述流体动力学方程(从而NS方程无唯一解,或无解),但是其隐含的理念是如此的。我记得多数的国际流体力学权威期刊 (Fluiddynamics) 接受文章的前提是:无条件接受NS方程,不接受反对NS方程的论文。不知现在是否还是如此。这是造成NS当前绝对主导地位的基本原因,而非是NS本身的威力。

因此,我认为(这也是我在研究中的论题),流体力学未来的进展在于:

  1. 在给定尺度下,整体的平移、转动虽然对应变没有贡献,但是对于流动是有决定性作用的,需要有专门的运动方程来表达;
  2. 压力的变化应与变形应力区别开来,从而,压力不应是作为静水压力来表征,而是作为以整体流速、涡度为自变量的二阶张量来表征;
  3. 固体的连续性概念不能简单的移值到流体中,需要有新的数学表达形式,而这种表达方式更类似于以物理属性为标准的办法,而不是以实质性的微元流体来表征。

总而言之,流体运动的复杂性在于我们无法把固体力学理论(变形力学)以一一对应的方式移值过去,也无法把质点系力学理论(拉氏力学,哈氏力学)以一一对应的方式移值过去,而必须是二者的某种组合。而这种组合又必须是物理学或力学原理所决定的,而不是我们简单的分解为独立的两种运动的叠加(加法的或是乘法的)结构。

我认为,也只有在有了上述的大致结构后,热力学量(熵,温度)、化学量(化学反应)、及分子动力学效应(多相,多组分)等才能有比较可靠的理论途径来进入对复杂流动问题的研究中。

无论如何看待NS方程,流体力学的高难度决定了它是本世纪的热门研究课题。而把NS作为单一路径的研究方案注定是令人失望的,而几乎所有的研究者对修改NS的各类方案都是持怀疑态度的,更不用说形成前赴后续的连续攻关形势了,而达成共识就是很遥远的未来的事了。因而,多数力量还将是针对NS方程,但是无非是后人把半个世纪前后的前人所做的全部研究工作再来一遍而已,研究结论本身是不会变的,变的只不过是说法。从科学哲学上看,这是一个以NS方程为中心的逻辑怪圈。

来源:肖建华科学网博客

解ns方程_流体动力学NS方程的哲学缺陷相关推荐

  1. python迭代法求解方程_迭代法求解方程(组)的根

    首先,迭代法解方程的实质是按照下列步骤构造一个序列x0,x1,-,xn,来逐步逼近方程f(x)=0的解: 1)选取适当的初值x0: 2)确定迭代格式,即建立迭代关系,需要将方程f(x)=0改写为x=φ ...

  2. Python(分治算法)问题 E: 解方程_求方程f(x)=2^x+3^x-4^x=0在[1,2]内的根。

    问题 E: 解方程 题目描述 求方程f(x)=2^x+3^x-4^x=0在[1,2]内的根. 输入 输入m(0<=m<=8),控制输出精度 输出 三个实根(根与根之间留有空格)输出方程f( ...

  3. matlab 定义一个有自变量的方程_常微分方程:(第四章) 高阶微分方程

    参考<常微分方程>第三版(王高雄) 常微分方程王高雄 第四章 高阶微分方程_哔哩哔哩 (゜-゜)つロ 干杯~-bilibili​www.bilibili.com 对于高阶微分方程,线性部分 ...

  4. 数值方法3:偏微分方程1:使用有限差分法解一维热传导(扩散)方程

    初学练习,看b站课程,教学为matlab代码,自己写的Python代码,后面会放b站课程链接,感兴趣的同学可以学习观看. 说明:Python初学者,代码可能不够漂亮,欢迎大家批评指正.本系列代码用no ...

  5. l20范数最小化求解系数方程_贪婪组稀疏方法(Greedy group sparsity)

    l20范数最小化求解系数方程_贪婪组稀疏方法(Greedy group sparsity) 本文章部分参考Fast group sparse classification l20范数最小化求解系数方程 ...

  6. 参考平面及其高度_柱面及其方程

    柱面及其方程 学习李建平教授空间解析几何 柱面的定义 柱面方程及其特点 准线位于坐标面母线平行于坐标轴的柱面方程及特点 常见柱面 一般柱面方程 柱面的参数方程 例1 例2 参考资料 回顾:圆柱面 圆柱 ...

  7. python迭代法求解方程_第一部分:趣味算法入门;第六题牛顿迭代法求一元三次方程的根...

    100个不同类型的python语言趣味编程题 在求解的过程中培养编程兴趣,拓展编程思维,提高编程能力. 第一部分:趣味算法入门:第六题SRE实战 互联网时代守护先锋,助力企业售后服务体系运筹帷幄!一键 ...

  8. 【转】流体动力学控制方程(详细推导)

    [转]流体动力学控制方程(详细推导) diyhoos 原文链接:https://blog.csdn.net/qq_42020563/article/details/80940387 CFD建立在流体力 ...

  9. matlab中solve解方程,matlab解方程组.我的每一个方程都很长,如果用solve,solve里的参数就会很长.有什么办法呢?...

    问题描述: matlab解方程组.我的每一个方程都很长,如果用solve,solve里的参数就会很长.有什么办法呢? s= (((22*r)/5 - 479/50)*(7*r^2 + (149*r)/ ...

最新文章

  1. JQuery-学习笔记05【高级——JQuery动画和遍历】
  2. 电子商务实战课程-张晨光-专题视频课程
  3. python中list是什么类型_Python中的list数据类型
  4. dhcp服务环境搭建
  5. 【撸码师的备忘录】JedisPool.returnResource()遭弃用
  6. 一个简易的C++性能分析工具
  7. 链表节点的删除(链表data升序有重复)
  8. poj3159 Candies 2012-09-07
  9. Ubuntu 系统安装显卡驱动教程
  10. QT实现低延迟的RTSP、RTMP播放器
  11. 英伟达各版本CUDA地址以及与驱动对应关系
  12. Matlab中批读取DCM文件
  13. 一键自动下载百度美女图片
  14. 文件分配方式-索引分配
  15. /etc/shadow文件是什么?
  16. MFC (如何通过点击botton打开一个文件夹/文件)
  17. 全球喷气发动机收入预计2028年达到1189.6亿美元
  18. 蚂蚁研究员玉伯:我的技术人生答案
  19. cp: omitting directory xxx
  20. 转贴:[铿锵发金石 幽眇感鬼神] 高瞻“家园”诗4首浅议[蛮桂]

热门文章

  1. PHP 报错 Use of undefined constant prop_values - ass...
  2. magento 在产品页添加评论 Add Review Form in Magento Product View Page
  3. 类似 Observer Pattern 的 NSNotificationCenter (实例)
  4. SELinux入门:了解和配置SELinux
  5. Linux——Vi/Vim如何消除搜索后的关键字高亮
  6. 《剑指Offer》 变态跳台阶
  7. 可迭代对象与迭代器 0318 草稿
  8. 看脸 高效学英语 下
  9. linux指令:软连接与历史命令
  10. iOS实现字符串动画