肖建华的个人博客分享 http://blog.sciencenet.cn/u/肖建华

博文

湍流的解析方程与湍流的普遍性质

已有 946 次阅读 2020-5-28 09:14 |系统分类:科研笔记

 

 

       湍流的解析方程和相应的解析解,被定义为世纪难题。

       由于科学网上有不少热爱科研的网民,其中有关心湍流问题的,因此本博文把我最近几年的研究结果给出一个哲学性论述。希望对于有关研究者开拓思路有所帮助。

 

       从实验上看,湍流现象类的实验很多,经验性的关系总结也很多。从数值模拟上看,大多数实验结果也能大致的由数值解得到。因此,对于实验力学和计算力学来说,湍流的方方面面的性质都已经被实验所揭示。

       从理论上看,则缺乏关于湍流的理论解。已经知道的是:湍流的运动方程是非线性的,所以也基本上不认为有解析解。

      

       2019年,应用陈至达的变形张量S+R分解理论,对于零压力梯度的壁面流动,得到速度剖面Uy)理论方程的形式为U’’=  -1/v(1-cosA)U。它是一个以函数为系数的谐波方程。式中,A依赖于速度的一次偏导数,为描述微元体涡旋性质的局部转动角参数。

       已经证实的是,普让特的对数律速度就是方程的理论解。因此,可以认为:对于理想的壁面流动,理论解与实验解是吻合的。

      

       由这个方程,随着壁面距离的增大,湍流的尺度是从超高波数的微小尺度演化为趋于零波数的超大尺度。而这种演化局部转动角函数的变化。有三个特征是明显的:

       1)大局部转动角区域,cosA几乎为零。充分发育的湍流在本质上就是谐波方程的解。因此,湍流的谱描述是普遍有效的。

       2)微小局部转动角区域,cosA几乎为1。湍流不发育,流动解表现为近似的线性解。这是经典流体力学。

       3)最为复杂的,也就是最为关心的湍流区域是cosA 从不能近似为1演化到接近于0的区域。这个区域的已知实验解是功率谱的-5/3律,把波数作为自变量。

 

       目前,能确切的外推的一般性湍流方程形式为:在涡旋中心系下的普遍方程U’’=       -1/v(1-cosA)U + F(U),,这里,F(U)指其它次要因素项。

       式中,速度为3维矢量形式。A3维的微元体局部转动角。

 

       通过这个壁面流理论方程特例:可以确认的是,用陈至达的变形张量S+R分解理论给出的加速度取代欧拉方程的加速度,就可以得到湍流的普遍有效方程。但是,由于局部转动角函数需要首先确定(假设、或是经验性关系确定),因此,普遍有效的解析解是难于得到的。

       这就需要实验研究的介入。由于局部转动角是涡旋的几何特征,所以会出现如下问题:

       4)对于一个涡旋群,在小尺度微元体度量时,局部转动角很大,而在取大尺度微元体时,局部转动角很小。实验上,局部转动角表现为尺度的周期性函数。

       这个问题被批判为:尺度依赖性。

       它的后果是:湍流解依赖于尺度的选择(所取的局部转动角函数)。

       哲学上看,这就会有争论:对于客观上的同一个湍流运动,在不同尺度选择下,给出不同的理论解。

       是是非非?

       5)数学上,严格的加速度公式是用李导数来证明的。因此,用S+R导出的微元体加速度与李导数虽然在本质上一致,但是在力学(物理)解释上区别很大。

       而普遍接受的是基于李导数的欧拉方程,或是NS方程,因此,对于这里给出的壁面流方程以及湍流的普遍方程,在理论界几乎没有支持性文献。

 

       因此,何去何从取决于个人的哲学判断

 



http://blog.sciencenet.cn/blog-39419-1235294.html

上一篇:数值求解法逼退解析求解法
下一篇:科学发展阻力在个体上的现实表现

3 朱伯靖 宁利中 黄河宁

该博文允许注册用户评论 请点击登录 评论 (2 个评论)

数据加载中...

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2020-7-13 04:32

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部