机器推导在空气动力学计算中的应用.doc

1、1空气动力方法在发展相对论中的应用杨新铁 (西北工业大学 503,翼型中心,西安 710072)赵双任(德国 慕尼黑 comroad 卫星定位所 )摘要 首先针对 Navier-Stokes 方程和 Maxwell 方程的相似关系一直缺少一个对应方程的问题分析了原因,借助于引入了非牛顿粘性流体的松弛效应,最后补充上了缺少的这个对应关系,使得不可压的 Navier-stokes(后面简称 NS)方程和 Maxwell 方程完全对应了起来.从而引起了利用粘性可压缩流体方程来进行对Maxwell 方程进一步的非线性化的探讨.于是先从速势流动开始,研究空气动力学的可压缩波动方程的种种变换,找出一种拟洛

2、伦兹时空变换使它变为不可压流的波动方程.说明了在声学波动方程的数学描述上,可压缩流和不可压缩流加上相对论时空变换只不过是相同客体的不同数学表达.后者是空间上二级精度,时间上一级精度.可压缩因子（ ）将根据表达形式不同出现在可压缩流动方程的系数中或者不可压缩加时空变量变换的数学描述的时空伸缩延迟系数中因为 Maxwell 方程的波动方程和不可压流是完全一样的，所以也期望上面得到的结论可以引伸到 Maxwell 方程的讨论中，为了探讨电磁场和引力场可能存在的介质规律.对于无粘可压缩流动,借助于卡门 - 钱学森在空气动力学中应用的切线虚拟气体法,得出了与质能关系类似的规律.对于可压缩粘性流体和电磁场

3、方程的相同数学结构,作者还缺少明了的结果,但是美国宇航工程师 Paul 在 AIAA Paper99-2606,0562 上所发表的平行的表达可以预先把可压缩因子提取到通量项的外面,给出了一种表象的处理方法，本文对 Paul 在对 NS 方程提取了可压缩因子以后，又重复引入洛伦兹变换是否必要的问题进行了讨论.尽管如此 Paul 给出的矢量通量漩涡表达方式是一种新描述，希望这种重新表达的力和漩涡的关系能够对引力场的研究以及对 Maxwell 方程组的强非线性化带来生机.关键词：NS 方程,相对论,超光速,质能关系现在关心相对论的发展的人越来越多了.当我们站在伟大的相对论奠基者建立的丰富和优美的理

4、论基础上的时候不能不回顾历史,发掘他们的足迹.首先让人想起的是麦克斯韦尔在结合法拉第的电力线理论和流体力学的场理论以及亥姆赫兹漩涡理论中作的优美描述.原来这些描述都是在无粘流动中建立的.多年来,流体方程组和 Maxwell 组的相似关系一直缺少一个对应方程组,就是引入了粘性也是如此.借助于引入了非牛顿粘性流体的松弛效应,可以补充上这个对应关系,从而使得不可压的 Ns 方程和 Maxwell 方程对应了起来.这就引起了利用粘性可压缩流体方程会对 Maxwell 方程更进一步的非线性化的兴趣.可以先从速势流动开始,研究空气动力学的可压缩波动方程的种种变换,找出一种拟洛伦兹时空变换使它变为不可压流的

5、波动方程.这就说明了从声学波动方程的数学描述来看, 相对论变换加不可压缩流等于可压缩流.协变不变性以及所谓四维空间度规不变只不过是可压缩流体的不同数学表达.为了探讨电磁场和引力场可能存在的介质规律.对于无粘可压缩流动,我们还可以借助卡门 - 钱学森在空气动力学中应用的切线虚拟气体法,得出了与质能关系类似的规律.虽然对于可压缩粘性流体和电磁场方程的相同数学结构,还缺少明了的结果,但是美国宇航工程师 Paul 在 AIAA Paper 上所发表的平行的的表达,却给与了希望和引起讨论.即这种重新表达的力和漩涡的关系有可能对引力场的研究以及对 Maxwell 方程组的强非线性化带来生机.下面分六个部分

6、来叙述这个想法.首先我们回顾以下这些熟知的事实. 一牛顿流体框架内已经有三个 Maxwell 方程和流体力学方程相似,2电动力学基本方程组是： (E) = (E)/t = H + E ( H)/t = - E ( H) = 0 其中 是电荷密度,E 是电场强度，H 是磁场强度，是磁导率， 是电介常数，是电导率。下面让我们在连续介质力学领域里寻找和它相类似的方程.由于电动力学方程组创始人 Maxwell 在建立他之初,就利用了流体力学的亥姆霍兹定律和法拉第的电力线理论,所以上面的方程 1,2,4 自然在流体力学方程组中有它的对应表达形式:对 Maxwell 方程组的 1）式来说,在连续介质力学的

7、类似表达,可以用涡的散度为零来对应: = 0 其中是介质流场速度的旋度:= V对 Maxwell 方程组的 4 式可以类比重力场的守恒法则得到F = 其中 F 是加速度， 是和万有引力有关的系数， 是质量密度。对 Maxwell 方程组的 2 式对应表达,，需要一个涡的变化也能产生力的旋度的表达形式，亥姆霍兹定律正是这样一个关系，对黏性流体可以通过考虑柯罗柯-兰姆形式的动量方程把亥姆霍兹定律写成更一般的形式，一般的牛顿流体动量方程为:V/t +VV =F-1/ P+ 1/( V)+ 1/( (2) 其中 V 是速度矢量, F 是引力,P 是压力, 是密度(在电场里面已经用来表示电荷密度), 和

8、 是粘性系数,和体粘性系数, 表示应变率张量。| 11 12 13 | = | 21 22 23 | | 31 32 33 |利用 VV=(VV/2)+ V 把上方程可以改写成柯罗柯兰姆形式 (Crocco-Lamb)如下:V/t +(VV/2)+ V =F-1/P+ 1/(V)+ 1/(2) 对上式两边取旋度得：/t + (V)= F-(1/P)+ 1/(V)+ 1/2亦即： /t = F +-( V)-1/P+1/(V)+1/2 由于方程最右面的第一层括号里面实际是压力,离心力和粘性力之和，把它合起来称为 F2。于是有：/t = F + F2 其中 F2=-( V)-1/P+1/(V)+1

9、/2 对不可压流体,中间的这一项 1/(V)为零.至此电动力学一共四个方程,我们已经找到了其中的三个在粘性流体力学中的对应形式.但是寻找最后一个(方程 3)的对应形式却遇到了在牛顿流体框架内不可克服的困难.因而我们不得不考虑如何突破它 3二,不得不放弃牛顿流体同构关系寻找方程 3 的对应方程遇到困难的原因是牛顿流体的线性假设所导致的.所以为了得到完全和电动力学一致的方程组而又不想引入刚体介质这样的约束,就应当考虑放开牛顿流体对应力应变率关系假设的限制，牛顿流体假设的要害是应力和应变率成线性关系。即主应力方向上表达式为:= - P+ 11+ 22+ 33 其中 a, b 为常数， ii为应变率主

10、应力方向张量元素。在这个假设中，缺乏我们所需要的应力和应变的非线性关系，显然这是牛顿流体假设所限定的.近年来非牛顿流体得到了长足的发展,使得我们可以考虑,如果我们利用考虑非牛顿流体的办法来放宽这个限定.现在我们可以选用最简单的一种非牛顿流体应力应变关系,这种非牛顿流体在无限小的时间间隔内应力和应变成线性相关的，这样无限小的时间间隔内主应力方向上表达式为:= - P +a * 11 + b* 22 + b * 33 其中 *， * 亦为常数， ii为应变张量元素，它代表流体内部名义上的位移应变的积累效应.所谓积累效应是指,在流体内部,应变并不可能永远保持下去, 由于分子的自由运动和湍流涡交换的作

11、用,在流动的每一个微运动距离上产生的微应变都是要经过衰减的,然后叠加上新产生的应变.最后总的效应应当是这些不同位置上产生的应变又经过松驰衰减后的总和. 衰减模式可以采用最简单的一阶衰减,也就是指数衰减模式,对于稳定流动只有经过一段相对长时间 的对这种衰减效应的积累以后 ,才能够得到真正的名义上的应变,从而利用它和应力相关的假设得出应力.由于这种非牛顿流体的模型是在牛顿流体概念上的延拓,所以对稳定的平直流动,从长时间平均的角度上来看,应当又回到牛顿流体的同构关系。我们可以用平行层流剪切运动来说明这个问题，在每一个距离微元 ds 上,在时间 间隔内产生的应变是 Vi/Xj (其中 V 是 X1,X

12、2,X3 和哑时间参数 的函数).用随时间在介质中任何瞬时产生的应变衰减的最简单的模型-指数衰减模型来处理，于是应变的实际值就应当是在整个路程上产生的应变的积累：其中 k 为表示衰减快慢的常数.比如空气 k 就应当取比较大的值,说明它衰减的很快,而流体 玻璃取比较小的值,衰减的很慢,k=0 不衰减.引入非牛顿流体应力关系式以后就可以从类似于方程的非牛顿黏性流体动量方程的柯罗柯-兰姆形式得到:V/t+(VV/2)+ V = F-1/P+ 1/2 其中 是基于应变的粘性系数,上式也可以写成F+ F3 = -1/2 其中: F3 = -(V/t +1/P +(VV/2)+ V) = -(V/t +1

13、/P + VV) F3 表示出了流体中除了引力以外的其他的力:惯性力,压力,迁移加速度引起的惯性力,以及部分粘性力.在上面的 15 式等号两边同时取对时间的偏微分得：F/+ F3/t =- /2/ 而方程左边的每一个元素对时间的偏导数均如下: t)-(k)( tdexVxVijjitij tdtdt4实际的微分结果说明,各量对时间的变化，在一个很小时间间隔内，松驰效应还显示不出来，也就是: ij)/ = (V i/xj+Vj/xi);故有 2)/ = V + (V) T ,此处上标 T代表转置 其中 在后面我们让它吸收 1/ 以后仍然把它记做 .也就是说：新的方程中间的应变张量那部分表达式对时

14、间的偏导数刚好和原来牛顿流体动量方程中的应变率的表达形式是一样的。所以对不可压流就有:/2 / 2 = 2V = 2V 通过同样的推导过程，我们可以得到和新的粘弹性连续介质力学理论中的类似于电动力学方程 3 的表达形式如下:F/t + F3/t = rot 这个公式的物理意义是引力和流体内各种内力的时间变化率和涡强度的旋度成正比三,两个不同领域物理现象的相似对比和物理意义 现在我们终于可以把这两种不同的方程组按对应项写在下方了:电动力学方程组 连续介质力学方程组(E) = F = (E)/t = H + E F+F3)/t = ( H)/t = - E / t = (F+F2)( H) = 0

15、 = 0上述的两个方程尽管有不同的地方,比如连续介质描述中在力的表示中增加了压力梯度,粘性和牵连惯性力的影响,但是它和电磁场的数学描述还是在形式上相似的.右面这个方程组的物理意义是，力的随时间变化产生介质的涡的环量，涡的脉动又可以生成力的环量。也就是说,电磁场中的磁场 H 和不可压缩流体力学中的涡是相似的，而电磁场中的电场强度和不可压缩连续介质场中的力地位是相同的。因此我们就会联想Maxwell 方程是否可以有 Navier-Stokes 方程类似的性质.由于不可压 NS 方程还可以延拓到可压缩流动的 NS方程,而可压缩的影响自然引入了(1- 2),那么我们是否可以考虑麦氏方程的类似可压缩性的

16、延拓?在这一方面,尤其是在无粘性的欧拉方程以及速势方程和电磁场方程的相似性方面已经有很多人从不同角度做了工作.为了避免可压缩粘性流动的复杂性,我们可以用简单的用速度势表达的可压缩波动方程来进行分析,由于不可压缩流体介质的波动方程和电场用的波动方程形式都是一样的,所以从速度势表示的波动方程出发来进行分析的结果也可望用在相同数学结构的电场波动方程上.下面让我们来寻找可压缩性在方程结构方面的影响的表现的一些规律性.在这方面我们发现夏皮洛曾经早就作过工作,他在可压缩动力学与热力学一书中作过这样的论述,即有一种空气动力学里面不常用到的变换-洛伦兹变换,用它可以把线化的可压流波动方程变到不可压缩流动.也就

17、是从声学上来看,线化的可压缩流的波动方程和不可压流波动方程加上相对论时空变换在数学上来说描述的是一个客体.但是经过严格的推演发现这个变换和严格的洛伦兹变换还在时间项上面有一点出入,后来夏皮洛在此书以后的版本中又删去了这一段.今天我们可以很容易用数学软件(maple)推出这个变换来,从形式上来看,它和洛伦兹变换相去不远的,在空间描述上有和洛伦兹变换相同的精度,时间描述上和洛伦兹变换相比是一阶精度,差别仅仅发生在时间项上面,好在时间方面还没有有利的试验能说明,这种区别,所以我们下面来详细讨论它,希望把它应用到全部方程组里面去.t)-(k)(/)( tdexVxVttij ijjit 5四可压缩流动

18、中的波动方程等效于声速不变原理加上拟洛伦兹时空关系为简单起见,讨论采用无量纲形式,这样线化亚音速可压缩流的波动方程为： 而不可压流以及引力场的波动方程为： 其中 =v/c=马赫数,c 式波动的速度,首先让我们利用待定系数的方法寻找一个时间空间的线性变换, 能够把把右边的可压缩流动方程变换到左边的不可压缩形式方程.另外还要补充两个条件:第一,这个变换对待不可压缩流(静止系统)要有罗仑兹变换尺缩的性质.第二它还要对不可压缩流动的时间膨胀的性质.按照这些给定的约束条件建立的方程,我们可以很容易用机器推理(maple6)求出这个变换的待求系数,于是得到其中合乎意义的变换如下:显然这个变换不可能是协变不

19、变的,因为我们要把可压缩流动变换到不可压缩,方程形式必须改变,那么它就不可能是协变的,而协变不变性只能用在不可压波动方程之间的变换,亦即从达朗贝尔方程到达朗贝尔方程,这个把可压变到不可压,也就是把声速从可以超越变到不可超越的变换的逆变换为:我们不妨把这种变换称为拟洛伦兹变换.从空气动力学上看,把波动方程的时空变量采用拟洛伦兹变换,就可以把它变成一个在洛伦兹时空里面的不可压缩流的波动方程_达朗贝尔方程, 在这个相似变换里面,不可压波动方程的系数(1- 2)可以通过时间空间的(旋转)洛伦兹量换把它变换走,变换到时空关系里面去了,从而变出来一个标准的达朗贝尔方程(不可压流).这是一个非常有意思的事情

20、, 他说明对可压缩流动的一个客体的数学描述,其(1- 2)的因子既可以在方程的导数项里面出现,也可以在时空的相对变换项里面出现,只不过是要看我们采用什么不同的数学系统来描述它.用可压缩方程描述客观的时候,(1- 2)的因子就出现方程的导数项里面,而用不可压方程来表示可压缩流动的时候,它就出现在相对时空变换里面,不论在何处,都可以把他看成可压缩性的影响因子,前者就形成了协变不变的一系列数学描述的框架,后者就是真实的迦里洛时空中的可压缩流动. 甚至同构的数学表述还不止这两种,比如我们还可以引入普朗特变换等等,这又是一种新的时空构架.这也就是说,我们常识所依赖的普通空间的可压缩流动的波动方程,和类似

21、相对论空间的不可压波动方程,实际上是同样一回事情,起码他们的数学表述是一样的. 于是使人自然产生这样的想法,可压缩流动与不可压缩流动的相似变换关系,和电磁场波动方程以及它所遵循的时空变换关系,虽然看起来是相差很远两个领域的效22222 .)1( YXtYXt 222211tzyxtzyxtzyxtzyx111226VW1vd-dw 3W22W /*c /dwd /*d /w /dvdWWV 应，但是有一样的数学描述.这里面有没有物理上的内在联系?既然电磁场波动方程和不可压流加上拟洛伦兹变换雷同,那么它是不是也隐含着一种带有可压缩性并且确立于实在的迦里洛空间的数学描述?更进一步,整个洛伦兹时空中

22、的 Maxwell 方程组是否也和实际上迦里洛空间的可压缩粘性流体方程组有类似关系?能否借助空气动力学方法来探索 Maxwell 方程的这种强非线性化的表达形式?五,亚光速下的质能关系既然希望利用空气动力学的一些方法,那么首先应当把质量和能量联系起来，首先用最粗旷的假设,把压力看成能量的一种量度,那么质量和能量之间的关系犹如理想气体的质量和压力一样的关系一样有:dW/d = C 2 其中 W 代表能量, 代表质量然而质量和能量之间的关系我们事先并不清楚,但是我们总可以假设 W和 或者说 1/ 之间有一个单值关系,如果把这个关系用条曲线 W = f(1/)来代表.然后我们利用卡门_钱学森的切线近

23、似假设,即在 =v/c 把 21 式代入 22 式得到,下式中下标 0 代表滞止状态,下标 表示没有扰动时的状态,实际对应于流场无穷远处,:dw/d(1/) = 2 C 2 = 2 C 2 = 常数 = 0 2 C 0 2 下面再考虑流动介质的欧拉方程: dw = -VdV 那么对方程 24 两边除以 ,然后对两边从 W 到 W 的积分就可以写成：而又因为所以就有了考虑第六式中假设的关系 c 2* 2 是个常数 可以提到积分号以外，于是上方程右边积分外边的部分可以写7 322W /*cd/ d成 c 2* 2，于是右边就继续等于：积分后两边分别得:v 2 - v 2 = -C 2 * 2*(1

24、/ 2 -1/ 2)= C 2 - C 2进一步,由于滞止时的速度 V 0 = 0, 而我们假想这种静止滞止状态下光速对应值称之为 C0,同样道理,就有: v 2 = C 2 C0 2 把此式两边同时除以 C 2,得到: V 2/C 2= 1-(C0/C) 2 = 1 - (/ 0) 2改写一下, 于是我们从另一条路,得到了和相对论相同的质能关系. 系静止质量, 0 系表征能量的总质量其关系为:值得指出的是, 得到的质能关系是沿着扩大了的连续介质方程的可压缩性(扩大了的麦克思维尔方程)走下来的.它的结果在 v/c1 时和相对论的结果相容.六.美国宇航专家也用粘流方程延拓相对论美国 AIAA P

25、eper 99-0562 发表宇航工程师,AIAA 会员 Paul A.Murad 也发表文章提出用粘性流体方程消除相对论比表达是里面的奇点,并把质能关系延拖到超光速领域,NS 方程的一般表的形式为: f/ t + Fx/ x + Gy/ y = Hx/ x + Hy/ y .其中符号 ( ) 表示偏微分.里面含有 (1- 2) 刚好是爱因斯坦的质能关系和洛伦兹变换因子,假设真空还是有物质存在,并且近似满足纳维尔斯托克斯方程,如果能够从他里面把(1- 2)提取出来,就造成了新的相对论方程,但是可以包括超光速.他假设提取因子後的势矢量为, 于是从上面的 NS 方程可以定义:(1- 2)/ * /

26、 x = Fx-Hx; 1/ * / y = Gy-Hy;这里 =v/C, 而 是个积分因子.把的定义带入 NS 方程得: f/ t + (1- 2) * / x/ x + / y/ y进一步写成* f/ t + (1- 2) * 2/ x 2 + 2/ y 2=1/*(1- 2) / x * / x + / y* / y对于定常流:(1- 2) * 2/ x 2 + 2/ y 2=1/*(1- 2) / x * / x + / y* / y把这种关系利用到涡动力学上,以便和电磁场的涡方程相联系,找出相同规律,从 collocos 涡理论,基于不可压流涡动理论的基础上延拓出来.新的因子提出了建

27、立在连续,动量,能量原始定义基础上的广义可压缩旋涡矢量的定义.这样定义的广义涡采用 定义 ,写出来就是:cvand/.;1/20 8* = (Fx-Hx)/ y -(1- 2) (Gy-Hy)/ x= / x*(Fx-Hx)- / y *(1- 2)* (Gy-Hy)/ x特殊情况下无旋流动方程简化为: (Fx-Hx)/ y -(1- 2) (Gy-Hy)/ x=0当涡在流场改变时其表达式变得复杂:= 1, 2, 3, 4; 且 x = x/(1- 2)(1/2) 其中: 1= u/ y - v/ x + u / y-v / x 是密度,当然可以验证,不可压流,涡强就是速度的旋度.后面以此用在

28、电磁场理论中.可以看出 Paul A.Murad 的理论没有光速奇点,可以发展到超光速,显然 paul 已经注意到了相对论的质能关系转换因子和可压缩因子都是一回事情,已经开始从粘性流体方程着手,去提取这个因子,但是由于 NS 方程的复杂性,他也没有昨晚这个工作,只好把(1- 2)这个因子先除在矢量通量的外面,并且定义了可压缩流的广义漩涡,这是非常创新和独特的.的推导过程，已经隐含着使得方程得到了 1 2 的因子，但是，为了不突破相对论的框架，在最后还是使用了罗伦兹变换，实际上如前波动方程的分析中所述，可压缩波动方程用拟洛伦兹时空的不可压缩体系来描述的时候，自动就会引入尺缩和时间膨胀效应，对于欧

29、拉方程和方程尽管有(1 2)这个因子存在,奇点效应会在变换中抵消掉.所以按作者的观点没有必要再引入洛伦兹变换.可惜 NS 方程的线化如同 Paul 所说,会得到一个非常复杂的表达,作者一直没有能够把里面的可压缩因子分离出来,这里还恳请对 Maxwell 方程进一步非线性化的关心和指导.参考文献1 连续介质力学， 朗道 19622 No Nuwtonian Fluids w.l.Wilkinson 19603。 VORTEX FLOW IN NATURE AND TECHNOLOGY，HANS J。 LUGT，19824。 No Nuwtonian Fluid mechnics G. Boehm

30、e 19865。 狭义相对论的实验基础。张元仲 19796。 近代物理问题，相对论质疑，郑铨 学术出版社 JSBN7-80045-038-227不可压缩流体运动规范场论的若干理论进展及其应用， 冉政 力学 20008Diskusstion about the similarity of Maxwell equations and mechanical equations of continuous mass, X.Yang; PIERS 1999.9. 转捩的新判定准则, 杨新铁, 力学 2000,10麦克斯韦尔方程和连续介质方程的相似性 杨新铁 力学2000,11.粒子超(广义)非线性速度,

31、时间的相对不对称效应和随体动力学公式 数学力学学报 杨文熊12.高速粒子的质量守恒性 数学力学学报 杨文熊13. Paul A.Murad,AIAA Peper 99-056214. Paul A.Murad,AIAA Peper 99-26069Subject:The application of aerodynamic method in development of spatial relative theory.Abstract:From the similarity between Navier-Stokes equation and Maxwell equation, which

32、lake one corresponds equation from whole 4 equations, an analysis is deduced, depend using the concept of relax effect of non-Newtonian viscose Fluid, the last missing equation is constructed. Through deduce give a concept to consider the Maxwell equation and NS equation as similar object but in dif

33、ferent viewpoint. On this background discussed the equivalence between compressible wave equation and incompressible wave equation with Lorentz. By using the Kaman-Tsian virtual gas method, the relation of mass and energy of relativity theory is given. In the last a parallel viewpoint of Mr. Paul in AIAA Pepper 99-2606 and its defect is discussed.Due to the mass-energy relation can be used to gravity, so the theory will satisfy all the experiment till now have made to validate relativity theory. Keyword: NS equation, Relativity, Super light speed, relation of mass and energy,