声传播推延时间计算对飞行器部件噪声预测的影响
摘要
飞机的噪声在航空交通中的影响越发重要，想要到达减小和控制噪声的目的首先要学会怎样去准确而有效的预测飞机飞行时而至的噪声。其中1个必须斟酌的关键就是声传播进程中推延时间的计算会对视察者所接遭到的声压带来怎样的影响。本文内容是在描写飞机噪声的产生和传播机制的基础上，通过建立单极子运动声源计算模型，研究如何求解推延时间方程和将其应用到预测实例中。
关键词：气动声学，推延时间方程，噪声预测
Abstract
The noise of the aircraft become more and more important in the air traffic, if we want to reach the purpose to reduce and control noise, we must first learn how to accurately and efficiently predict the noise of the aircraft. Which one must consider the key is the postpone time calculation will bring what kind of impact of the sound pressure received for observer in the process of sound propagation. The content of this article is to study how to solve the postpone time equation through the establishment of the monopole moving sound source on the basis of describing aircraft noise generation and propagation mechanisms and applied it to predict instance.
Key words: Aeroacoustics, The postpone time equation, Noise Prediction
目录
摘要 ………………………- 3 -
引言 ………………………- 3 -
第1章基本概念 ………………………- 3 -
1.1. 声压、声压级 ………………………- 3 -
1.2. 声波传播速度 ………………………- 4 -
1.3. 流体动力声源的分类 ………………………- 5 -
第2章相干理论 ………………………- 8 -
2.1. 脉动球源（单极子声源） ………………………- 8 -
2.1.1. 脉动球源声场 ………………………- 8 -
2.1.2. 声场对脉动球源的反作用 ………………………- 13 -
2.1.3. 单极子声源 ………………………- 15 -
2.2. 运动声源的多普勒效应 ………………………- 20 -
2.3. 格林函数 ………………………毛病!未定义书签。
2.4. 运动声源的声场 ………………………- 22 -
第3章算例 ………………………- 30 -
第4章结论 ………………………毛病!未定义书签。
致谢 ………………………- 34 -
参考文献 ………………………- 34 -
引言
从初期到现在，飞机的噪声已渐渐成为航空交通工具中的不良副产品。在航空技术发展的早期，人们将精力主要放在如何提高飞机的性能上，故而忽视了飞行器产生噪声的问题。但是由于航空技术的日渐成熟，开始投入了更多的资源以减少气动噪声。而且随着航空技术的成熟，公众和监管机构把注意力都集中在了安全，排放和噪音，而不是性能或效力方面。
随着时间的推动和科学的进步，声学作为物理学的1个重要组成部份，人们渐渐发现了它在利用科学、学术研究、国防技术、文化生活和社会利用等诸多方面的重要性，因此使得声学逐步渗透到科学技术的其他研究领域，并使之成为1个许多学科的边沿分支。而且声学不论从物理学的其他分支，还是从若干新兴的学科领域，它的确具有现代各个科学相互交叉的边沿学科的特点，已成为许多学科领域所关注的对象。而到了20世纪80年代后，随着对电脑的开发研究，电子计算能力的突飞猛进，法拉赛特等人终究用解方程的办法将对声场的计算广泛的带到生活中来，为声学在当代人们的生活中利用奠定了基础。
在传统声学观念中，声波的产生界说为物体表层的振动传递至流体中的模样，当弹性体中的振动波传送到其边部时，它相邻的体积元间遭到膨胀或紧缩使其外部元素产生相应的位移。而物体不同边界处的声源强度可以根据其当地表面的振动情况加以肯定。通过将它们理想化，把它们近似看成是由1系列点声源组成，对小振幅线性声场来讲，声场解就是这些点声源单独作用而至声场的迭加。这样既避免了繁琐的数学，又能揭露出声源辐射的基本规律。
对运动声源来讲，空气动力性噪声的计算都将触及到声源的发射和接收问题。当声源尺寸较大，声源为多个点声源的叠加时，由于不同位置的声源，它们在到达视察者位置时所历经的传播延迟时间是不同的。因此，在计算不同位置声源对视察者噪声的响应时，必须首先计算它们各自推延时间，并找到它们在各自的推延时刻发射声源的大小，以便将这些声源对视察者的噪声响应进行正确的叠加。
第1章基本概念
1.1. 声压、声压级
在传统声学观念中，声波的产生界说为物体表层的振动传递至流体中的模样。如图1⑴⑴所示，当弹性体中的振动波传送到其边部时，它相邻的体积元A遭到膨胀或紧缩是其外部元素产生的位移致使的，从而使靠近未扰动体积元B 的流体压力和体积元A内的流体压力间出现压力差，且在它的作用下，使其相邻的体积元B遭到膨胀或紧缩；后来由于体积元B内的流体压力与靠近未扰动体积元C 的流体压力间出现压力差，因而就造成振动波在流体中可以从近到远地传播出去。这类在流体间从近到远传波的疏密波便是声波。因而可知，声波是1种流体体积元振动方向与其传播方向相同的纵向疏密波。
图 1⑴⑴ 空气体积元的振动方向与传播方向相同的纵向疏密波
（1⑴⑴）
在公式（1⑴⑴）中：流体介质在t点的瞬时气压为P,此时在没有遭到声波干扰时的气压为Po。瞬时声压是指在声场中的1声压在某1时间的声压值。而有效声压是指在1特定时段，对瞬时声压的2次方取t的平方根，即
(1⑴⑵)
其实在平时生活中我们所听到的声音很广泛，从2×10⑸到2×102帕（牛顿/米2），而且人对声压大小的感觉和人的听觉能力是成正相干的，因此比较声压的大小1般使用声压级。声压级的定义为人耳可听阀声压的平方（pref =2×10⑸帕）的比值取对数的无量纲值和当地流体媒质中有效声压取平方后乘以10的结果，即
(dB) (1⑴⑶)
式中用有效声压取平方的缘由是声波的能量或功率与声压的均方值关联，故用其来肯定声波的强弱。分贝作为声压级的单位，记作dB。由上式(1⑴⑶)可知，声压每加1倍，声压级就会上升6分贝。对声压级这个粗略而抽象的数量概念，我们在这里可以举1些实际的例子。比如人耳对频率为1kHz声音的可听阀（理论上人耳能听见的最微弱的音量）为pe =2×10⑸帕，0分贝；很安静环境下微风吹动树叶的声音约为pe =2×10⑷帕，20分贝；放着广场舞音乐的操场上的声量约为pe =0.2帕，80分贝；常说的震耳欲聋的声音约为pe =20帕，120分贝；客机发动机附近5米处的音量约为pe =200帕，140分贝。
1.2. 声波传播速度
在探讨声波的传播速度时可以借助这样1个实例，空气和水相比较需要更长时间才能将体积元内的疏密波传传递到下1个相邻体积元，这是由于水媒质的紧缩性比空气小，所以在空气中声波的传播速度要小。因此我们可以推断，声波的传播速度与媒质的紧缩性有关。
在绝热环境小振幅声波作用下的理想气体的物态方程得到声速与当地空气的压强和密度的关系式为
（1⑵⑴）
当理想气体的密度、当地压强和比热比分别为1.27千克/米3、1.013×105牛顿/米2和1.402时，声速为
米/秒 （1⑵⑵）
这是拉普拉斯在1816年用绝热等熵条件得到的结果，他不但就此纠正了1687年牛顿应用波义尔定律，采取等温条件（）导出的声速米/秒的结果，而且还使实验值与理论值高度吻合。
由克拉柏龙公式我们可以推导出声速c0与当地媒质温度T0的关系为
（1⑵⑶）
式中：T、P、V为M千克气体的绝对温度、压强、体积；μ为气体摩尔量, R=8.31焦尔/开尔文·摩尔为气体常数。因此，声速也能够描写为
（1⑵⑷）
式中：t为当地摄氏温度,为当地绝对温度。
1.3. 流体动力声源的分类
根据英国科学家莱特希尔提出的声学类似理论，依照流体动力声源所对应的物理模型和发声机制，可对流体动力声源作以下分类：
第1类，脉动球源。脉动球源是指在经典声学中可以将由物体运动所致使对其边界上的流体产生紧缩膨胀作用的声源。例如，旋转叶片的运动进程，就是使位于桨盘平面上的流体连续地被具有1定厚度的旋转叶片移出和插入，在这1进程中，运动的叶素体积会呈周期性的占有当地流体空间，从而使得当地流体产生周期性地体积存缩而致使密度增大，这类密度增大可视作在当地流体内添加了1个流体质量，大小为被运动物体的体积。同理，质量源的生成也是由于流体活动然后不断生成和幻灭的气泡而至，因此也能够看成是脉动球源。作为第1类流体动力声源的物理模型，脉动球源的本质可视作1个插入流体的脉动质量源。为了使理论研究更简捷，物理上的脉动球源需要用数学思惟中的质量点源做1个理想假定，在这类理想假定下的脉动球源的模型就是1个单极子声源。
第2类，振动球源。振动球源是指在经典声学中由物体运动而至的物体表面升力对其边界上的流体产生推力作用的声源。例如，旋转叶片的运动进程，就是使旋转叶片不停的对位于桨盘平面上的流体作用上升力。此时，当地叶素升力就会对当地流体产生作用致使质点速度的周期性地变化，而当地流体压力的变化是遭到质点速度变化影响的，终究致使声波的产生。这类声源跟经典声学中所提及振动球源发声机理非常类似，但它的发声机制和在流体中的运动物体有关。为了使理论研究更简捷，物理上的振动球源可用数学思惟中的力点源做理想假定。而脉动球源即偶极子球源又被称作第2类流体动力声源。
第3类，4极子声源。由偶极子声源可知，由于对任1流体模块施加1个脉动推力就会构成1个偶极子声源，所以如果把由于物体运动所引发的物体周围紊流流体间的力看成是两个流体间的相互作用，就构成了4极子声源，在经典声学中把它定义为第3类流体动力声源。第3类流体动力声源的物理模型便可理解为流体之间的1对方向相反大小相等的脉动力源，而作用本质是紊流流体内的应力。为了使理论研究更简捷，物理上的振动球源可用数学思惟中的力点源做理想假定。这样4极子声源就能够通过理想环境下的1对脉动点力源的模型加以理解。
莱特希尔的理论下，采取经典声学的观点，在点源的假定下，分别对形如脉动球源那样的质量源、振动球源和湍流中流体检间相互作用应力源给予数学上的描写，终究得到3类流体动力声源的数学模型，这就是之前所说的单极子、偶极子和4极子声源。在下面的讨论中，重点分析两个方面的问题：1是声源产生声波的辐射特性，比如声压随距离的变化和声源的指向性，声压与声源在声场中的关系等；2是由声源引发的声场反作用于声源振动的影响，也能够形容为因辐射声波引发的附加在声源上的辐射阻抗。
1.4. 章节总结
由于专业的小方向是结构强度，所以在专业课方面所学的多以飞行力学，结构力学为主，在声学方面的触及其实不多。当1开始拿到选到的课题时，看到声传播，部件噪声这些词语真的有些蒙。赶快联系导师，导师让我在寒假时先多看1些外文资料，因而我在寒假时看了很多和声学有关文献，特别是噪声有关的文章。其中1篇简述冲击噪声的文章我进行了仔细的理解。
由于最近有人提出了1种新的联合对角化分数低阶时空DOA矩阵法，那就是估计不相干的窄带信号在脉冲噪声中存在的2维DOA。这类新的方法保存了原始ST-DOA矩阵方法的优势，那就是可以既不用搜索光谱灵敏度峰值的方式也不用成对匹配的方式来估计2维DOAs，另外，它可以和1维角度共同处理来源。摹拟结果表明，提出的方法为了相较于ST-DOA矩阵方法，能够更好地表现限制强烈的脉冲噪声，特别是对低信噪比。其实关于DOAs估计的2维方向问题现在引发了很多关注，特别是在雷达，通讯和智能天线等领域，固然也包括我课题中的主角飞行器。由于当下时髦用高分辨率技术来辨别多种密切空间的来源，在和特点结构配置方法相结合。比如说，2维音乐类型方法和2维欧洲信息技术战略研究计划方法。但是，这些方法遭受了很多困难，比如，在2维领域中，多维搜索光谱灵敏度峰值，非线性优化，1维配置和成对匹配，因此它不合适实现实时化。为了克服这些缺点，有人建议利用开发两个平行均匀线性阵列，发现1种可以以高效力计算的普遍的DOA矩阵方法，著名的欧洲信息技术战略研究计划方法已被证明是DOA矩阵方法中1特殊案例。传统的2维DOAs估计算法是以高斯假定为主导，不幸的是，在信号处理中有很多不是由高斯假定主导的现象，比如说大气噪声，脉冲噪声，人为噪声等等。最近，脉冲噪声可以作为1个复杂的对称α的稳定进程的模本（SαS），由于当2α的最小值的分散标准对评估SαS进程是1个很好的选择时，SαS没有有限方差。在1些利用中，可以通过介绍时域情况获得有效信息，在脉冲噪声环境中使用分数低阶时空矩阵的DOA估计方法和1个时空DOA矩阵（ST-DOA）算法，这是以时间领域中的推理为优先的。但是，不单单是2维欧洲信息技术战略研究计划方法，还有DOA矩阵方法和ST-DOA矩阵方法不能估计1维角度或1些曲线平面的信号，这下降那些算法的估计性能。
在那篇文献中，我们提出了1个新的联合对角化分数低阶时空DOA矩阵方法，这类方法既不用光谱灵敏度峰值，也不用成对匹配，就能够直接获得基于联合对角化的2维DOAs。虽然那篇外文文献讲的理论知识有些深，最开始翻译时也遇到了很多的困难，但我还是收获了很多。首先，让我对现今世界对噪声方面的研究有了1个认识，其次，我也更加深信的肯定我课题研究内容的实用性和对我未来工作的帮助。
第2章相干理论
2.1. 脉动球源（单极子声源）
由于脉动球源是1种球面在不停做均匀涨缩振动的声源，所以在球源表面上的各点都可理解是在做沿径向同相位、同振幅的振动。特别是当用1个点球源理想化替换脉动球源时，那末每个表面做任何振动的物体都可以看作为由无数个小脉动球源组成。对运动物体来讲，由于它表面每点的厚度不1样，所以运动时表面每点和所在地流体缩张幅度也存在差异。故而我们可以把运动表面的每一个点看成是彼此不同的脉动球源。由于物体表面每点跟所在地流体缩张运动及程度有关，同时触及到流体力学知识，所以称它为流体动力声源。且物体表面每点跟所在地流体的缩张程度本质是使当地流体密度产生改变，故将它定义为质量源，这也就是的第2类发声机理。在经典声学中的脉动球源是质量园最基础的物理模型，所以本章节将对脉动球源构成的声场进行讨论并解说第1类流体动力声源所构成的声场特点。
2.1.1脉动球源声场
用最原始的书本知识可知声波波动方程为
(第1节1式)
r 接收点处的声压p(r,t)r球心为原点r边界条件ra
图2⑴⑴ 脉动球源而至的声场
为了更好的解释球源表面的1些特殊条件，求出1个r为a的单极子球源的构成的声场，参照下面的公式，应用拉普拉斯算子的情势，得到1个新的方程。:
(第1节2)
然后重新用坐标定义1个新的球心，这个时候上面的方程就能够表达成下面的情势
而我们说的单极子球源呢，由于它构成的声场是特殊的，是1种对称模式，像上面公式中的声压p 单纯只是r和t的函数。所以说，上式就变成了下面的模样
(第1节3式)
由，(2⑴⑶)式也可写成
(第1节4式)
可以得到上式的另外一种通俗的解，为
(第1节5式)
式中：A、B为两个需要单独肯定的随情况变化的常数。将上式左侧r移到右侧，再结合上上个式子我们就又得到1个公式新的情势
(第1节6式)
由于在自由没有拘束的空间中是没有反射波的，所以这个时候的常数B =0。这样上面的式子呢就成了下面的模样
(第1节7式)
当中A1般为复数，A/r的绝对值即为声压的振幅。
将上面的1个式子代进质点速度和声压相干的式子中，可以得到
(第1节8式)
其中的意思就是质点速度的振幅。
对脉动球源所构成的声场的1个很形象的解释就是说，对小波动的能用函数来表示的声场来讲，照旧是有常规性存在的，可以单纯的来研究单极子球源所构成声场的1个单色波的特解。这时候我们来设出1个新的公式：
而球源表面流体媒质的质点速度可以描写为
(第1节9式)
如果将1定的条件式代入质点速度的1般解式，就会得到下面的这个公式：
通过1番整理，就可以够得到1个待定着的数值为
(第1节10式)
其中
，
把得到的1个A值放到本节7式便可就能够得到1个终究的单极子球源辐射声场的声压解，为
(第1节11式)
把A值代入本节的8式当中，便能够得到1个终究的单极子球源来辐射声场的单点V的解，为
(第1节12式)
把本节的7式和本节的8式结合到1块，就可以都得到反应声波传播优点的声场的声阻抗率，为
(第1节13式)
由13式可以清晰的看到，这个时候球面声场的情况很特殊，它和kr都没关系，不是1个固定的数。而且只有当r接近无穷时跟平面波才相像，所以我们又得到下面这个式子
这个结果实际上是完全可以预感到的。由于如果r非常大的话，波阵面的所有位置都能够等效看作是平面。
最后，将本节的7式和8式代入
(第1节14式)
可得到脉动球源辐射声场的声强，为
(第1节15式)
该式表明，球面声场的声强随距离r的平方反比地减小。而平均声功率的计算公式可推导为
(第1节16式)
由该平均声功率的表达式，我们可以得到很多帮助。
2.1.2声场对脉动球源的反作用
前面我们根据单极子球源所构成声场的情况描写了其声传播的特点，现在我们应用辐射阻抗这样1个新的量词来探讨脉动球源的声辐射特性。
首先我们从做功的角度来定义1个新的运动方程，即下式
(第1节18式)
在这当中呢，是单极子球源的1个振动V。事实上，在单极子球源带动它边界上的流体媒质产生疏和密交换并生成声波的同时，单极子球源自己本身也会遭到声场对它的反向影响，所以有了下面的式子
(第1节19式)
当是1个色彩的单极子球源的情况时，将本节7式和本节10式代入本节的19式就能够得到又1新的表达情势
(第1节20式)
将本节20式代入本节的18式当中，由声学边界条件，把球源表面流体媒质的质点速度用球源表面的振动速度代替，也就是，那末本节的18式就可以够整理为
(第1节21式)
如果用力学的视角来看，假定令
(第1节22式)
则本节的21式就可以够表达为
(第1节23式)
式中：、和分别为辐射阻、辐射抗和辐射阻抗。
不能不说，在我们所研究的问题当中，辐射阻抗是1相当重要的参数。
(第1节24式)
所有我们可以得到，如果单极子球源表面的振速始终不变，那末平均辐射声功率只是单单由单极子动球源的辐射阻所决定，和别的条件没有关系，因而我们还可以根据单极子球源辐射阻的式子
中看出，这个公式终究又取决于声传播媒质的特性阻抗、单极子球源的几何大小a和振动频率。
2.1.3单极子声源
当出现球源的半径远超过声波的波长时，我们将这类现象称作大单极子球源情况，有
(第1节25式)
此时本节的10式当中的不肯定定常数就可以写作：
当出现球源的半径远不及声波的波长时，我们将这类现象称作小单极子球源情况，有
这个时候本节的10式当中不肯定的常数就可以写作:
(第1节26式)
这样，应用数学上的1个思想把单极子球源理想成1个质量的点源方法，也能够理解为对单极子球源所构成声场的声压、质点速度和声强的解，应用本节的第10，12，15式，并都将式中的，然后求出质量点源而至声场的声压、质点速度和声强的解为
(第1节27式)
(第1节28式)
(第1节29式)
可以明显看出，由本节的27式和28得出的是1个色彩质量点源所构成声场的声场解，如果将上面的1个色彩的频率因子添加到函数中，而且让，这样就可以够得到不是1个色彩质量点源所构成声场的声场解，为
(第1节30式)
(第1节31式)
将(4⑵⑵6)式所描写的声强在离球心任意距离r的波阵面上，对其所包围球源的封闭面4πr2进行积分，便可得到单极子声源所辐射的平均声功率，为
(第1节32式)
这样，对物体表面不同厚度的点都可以看成是各自的1个单极子声源。在物体于流体中运动时，由于物体表面厚度的变化多是不同的，所以它们对相邻流体产生紧缩膨胀作用是不同的。因此，我们就用他们各自本身的单极子声源来表示它们，并求出它们所构成的声场。
为了便于描写单极子声源在任意位置而至声场的声场解，和便于数学上对多个不同位置处的单极子声源而至声场的声场解进行迭加，我们常需要将于球坐标下单极子声源而至声场声场解的描写情势转化为在直角坐标下的描写情势。
接收点声源
图2⑴⑵ 单极子声源与接收点的位置
将单极子声源声场解(2⑴⑵7)式~(2⑴⑶0)式改写成直角坐标下的情势，只需把接收点的位置和单极子声源的位置分别改写成直角坐标下的描写情势和，把接收点与单极子声源之间距离r改写成便可，参见图2⑴⑵，这样(2⑴⑵7)式~(2⑴⑶0)式在直角坐标下的情势就为
(2⑴⑶3)
(2⑴⑶4)
(2⑴⑶5)
(2⑴⑶6)
(2⑴⑶7)
关于球坐标与直角坐标间相互转换结果的正确性，可以通过把上述方程解和直角坐标下的拉普拉斯算子
代入到声波波动方程便可加以验证。验证时利用以下对球坐标下位置变量微分时的链式微分法则（描写采取求和约定的方法）：
(2⑴⑶8)
2.1.4格林函数
不失1般性，可以从单极子单色声源的声场解
(2⑶⑴)
看出，其求解得到的声场中某接收点处声波的振幅和相角仅与该接收点到声源发射点间的距离有关，并且通过变换，我们可以把声场解描写成相对位置的函数情势，为
(2⑶⑵)
式中：
(2⑶⑶)
为接收点与发射点间相对位置的函数，并称这个函数为格林函数。而接收点与发射点间的相对位置就为格林函数中的变量，常记作。由(2⑶⑵)式可见，声场的格林函数得到了，声场解也就得到了。因此，求解声场也能够说是通过数学方法来获得声场的格林函数。
格林函数除象(2⑶⑵)式那样，经常使用于描写由已肯定位置（或它在空间的运动轨迹）和时间波形的声源所产生的声场。另外，格林函数也可用于描写流体中的声波波动方程，例如，将(2⑶⑵)式代入声源而至声场的声波波动方程，即可得到格林方程为
(2⑶⑷)
式中： 为狄拉克函数，狄拉克函数在声场接收点与声源发射点不堆叠时为零，堆叠时对它的体积分为1。其直角坐标的描写情势为
明显，格林函数满足格林方程。格林方程与波动方程1样，也经常使用于描写声场。不论采取它们中的哪一种方程，终究取得的声场解是相同的。
2.2. 运动声源的多普勒效应
与经典声学不同，流体动力声源常随活动物体的运动而相对接收者作相对运动。当声源或接收者运动时，接收者地方视察到的、来自流体媒质中所传播过来的声波频率将产生变化，人们称这类由于声源与接收者间的相对运动而使声波频率产生变化的物理现象为多普勒效应。
图2⑵⑴ 由运动声源发射的波阵面
接收者所听到的声波的频率fL产生变化是由于，它与声波相对接收者的传播速度c和声波向接收者方向传播的声波波长有关。接收者所听到的声波的频率fL可描写为
(2⑵⑴)
式中的c为声波相对接收者的传播速度，它与接收者相对声源的运动有关。由于接收者的运动会使声波提早或滞后到达接收者处，有
(2⑵⑵)
这样，空气媒质中声波的波长便为
(2⑵⑶)
将(2⑵⑶)式、(2⑵⑵)式代入(2⑵⑴)式，便可得到接收者所接收到的声波的频率为
(2⑵⑷)
若接收者或声源的运动方向不在接收者与声源的连线上，则(2⑵⑷)式中的VL或VS应取它们各自在接收者与声源连线上的投影，即有
(2⑵⑸)
式中：接收者或声源的运动方向与声源到接收点辐射矢径间的夹角；和为接收者和声源的运动马赫数。
由(2⑵⑷)式、(2⑵⑸)式看出，接收者听到声波的频率其实不是声源发射声波的频率，并且这类差异与声源相对接收者的相对运动有关。理论上把这类声波频率随声源相对接收者的运动而变化的现象称为多普勒效应。
另外，还可以斟酌当传播媒质运动时，如空气中存在风速VM时，则(2⑵⑷)式将变成
(2⑵⑹)
式中：V M为风速在接收者与声源连线上的投影，VM为正表示风速向声源方向运动。
固然，多普勒效应也能够通过对接收点所感遭到的声波的相位变化加以描写，由于单位时间内声波相位的变化次数所反应的就是声波的频率。例如，1个频率为的单色单极子声源以马赫数运动，则声波传播到1个固定的接收点时的相位总可以描写为
(2⑵⑺)
如果声源以亚音速运动，则由(4⑺⑴9)式可以得到接收点接收到该声源的频率为
(2⑵⑻)
式中：
(2⑵⑼)
为多普勒效应因子。可以看出如果接收点位于声源之前，为正，接收到的频率增大；反之，频率将变小。如果声源以超音速运动，由于有两个辐射半径，所以接收点上将接收到两个频率
和 (2⑵⑴0)
的声波。这将致使两个频率不同的声波会同时落在该接收点上，从而使得该点的声压出现很大的起伏，并由此产生声的相干现象。
实际上，不但是频率，声波的强度也会由于声源与接收者间的相对运动而产生变化，例如，我们把称为运动声源前飞效应的声波强度因子。它也与多普勒效应因子有关。
2.3. 运动声源的声场
1个静止于空间的强度为的单极子声源而至的声场，其速度势为
式中：为接收点的空间位置；为接收点的接收到声波的时刻；为声源发射声波时的空间位置。根据声波的传播速度，很容易得出式中为声源发射声波的时刻。以后为便于描写运动声源而至的声场，我们把声源发射声波的时刻用表示。明显，声源发射声波的时刻、接收点的接收到声波的时刻和声波在传播路径上的传播延迟时间之间具有以下关系
(2⑷⑴)
明显，这个式子的作用是用来计算运动声源推延时间的，所以这个式子又常常被称作是推延时间的方程。
由狄拉克函数的筛分性质，1个强度为以速度运动的单极子声源可以表示为
(2⑷⑵)
该式由狄拉克函数定义了声源在发射声波时刻它所在的空间位置，即。明显，该式描写的是1个时位于坐标原点的运动声源。依照已前述的方法，该点声源而至声场的声波波动方程可以描写为
(2⑷⑶)
该方程对应的格林函数为
(2⑷⑷)
式中：为声源发射声波时刻的位置到接收点间的距离。方程(2⑷⑶)的解为
(2⑷⑸)
对该式关于空间位置变量的积分，可由狄拉克函数的筛分性质求解完成，即
(2⑷⑹)
注意，此时式中狄拉克函数中的变量为关于发射时刻变量的隐函数的描写情势，该隐函数实际为1个推延时间方程，所以(2⑷⑹)式可以写成
(2⑷⑺)
为了利用狄拉克函数的筛分性质，需要将(2⑷⑺)式中的积分变量变换为狄拉克函数中的变量。它们之间的变换由推延时间方程(2⑷⑴)很容易得到，为
(2⑷⑻)
式中：声源运动方向与声源到接收点辐射矢径间的夹角，如图4⑺⑴所示，为声源运动马赫数在辐射矢径上的投影。把(2⑷⑻)式代入(2⑷⑺)式得
(2⑷⑼)
式中：是推延时间方程的根，且、是其发射时间的函数。下面讨论推延时间方程的求解问题。
图2⑷⑴ 亚音速运动声源辐射声波的几何关系
求取(2⑷⑼)式的解，首先根据是想法求得推延时间方程的根。对1个初始时刻（）位于坐标原点的运动声源来讲，其推延时间方程可以描写为
(2⑷⑴0)
不失1般性，可以假定声源以速度沿轴的正方向作均速直线运动，这样(2⑷⑴0)式即可写成
(2⑷⑴1)
从而求得推延时间方程(2⑷⑴1)的根值为
(2⑷⑴2)
值得注意的是，式中开方可能出现无意义的解。应当指出，只有为实数，且的解才是推延时间方程(2⑷⑴0)的真解。
然后根据(2⑷⑴2)式求得推延时间就能够求解运动声源的辐射距离，为
(2⑷⑴3)
由于正实数的才有物理意义，所以对亚音速运动，(2⑷⑴3)式的根号前应取正号，这样就有
(2⑷⑴4)
式中：。由(2⑷⑴4)式和图2⑷⑴很容易证明R可以写成
(2⑷⑴5)
1旦推延时间和辐射距离得到后，即可根据(2⑷⑻)式求得发射角和声源运动马赫数在辐射矢径上的投影，它们分别为
， (2⑷⑴6)
由图2⑷⑴还可以看出运动声源的接收角为
(2⑷⑴7)
由上述求得的参数即可方便地获得运动声源而至声场的声压解，为
(2⑷⑴8)
式中：利用(2⑷⑴4)式、(2⑷⑴5)式得
(2⑷⑴9)
(2⑷⑵0)
最后将(2⑷⑴9)式、(2⑷⑵0)式代入(2⑷⑴8)式，便可得到亚音速运动声源而至声场的声压解，为
(2⑷⑵1)
式中：为运动声源的强度关于时间的导数。对(2⑷⑵1)式，当时，就变成静止声源而至的声场，为
(2⑷⑵2)
该式与前面推导的单极子声源而至声场的声压解完全相同。将(2⑷⑵1)式与(2⑷⑵2)式比较看出，(2⑷⑵1)式中的对传播距离r作了修正，而对传播距离修正的根本缘由是由于运动声源发射声波时刻的位置是它发射时刻和运动轨迹的函数。通常，可以把(2⑷⑵1)式中的称为运动声源的前飞效应。
2.4. 章节总结
在大4下这个学期中，虽然各模块的学分已基本修读终了，但为了更好的完成毕业设计的全部内容，我选修了我的导师张强老师的《气动声学基础》这门课，在课堂的学习中，我觉得在传统声学观念中，声波的产生界说为物体表层的振动传递至流体中的模样，当弹性体中的振动波传送到其边部时，它相邻的体积元间遭到膨胀或紧缩使其外部元素产生相应的位移。由于物体不同边界处的声源强度可以根据其当地表面的振动情况加以肯定，对小振幅线性声场来讲，声场解就是这些点声源单独作用而至声场的迭加，而如果当声源的尺寸较大时，声源为多个点声源的叠加时，由于不同位置的声源，它们在到达视察者位置时所历经的传播延迟时间是不同的。
湍流是流体的1种不规则运动。流体的这类不规则运动可以通过流体与物体表面的磨擦产生，也能够通过不同速度流体层间的相互磨擦产生。这两种方式产生的湍流为两种不同的湍流。前者为壁湍流，其湍流的产生与发展都受固壁影响，后者为自由湍流，其湍流的产生与发展都受当地环境的影响。但不管怎样，湍流的产生与发展都与黏性有关。1般，在真实的黏性流体中，如果没有连续的外部能量的补充，湍流运动就会在黏性作用下使活动的动能转变成热，终究衰变而消失。因此，湍流和黏性流体的1切活动1样本质上是耗散的。也正是粘性的作用，使湍流更加均匀和使它更少的依赖于方向，在极端情况下，流场各部份的湍流定量地具有相同的结构。
实际上湍流不但为随时间的不规则运动，而且同时也为对空间的不规则。因此，必须用时间和空间两种变量来描写湍流，并且在后面辨识湍流旋涡的尺度时必须用到时间和空间合成的相干函数概念。湍流的不规则性和随机性可以通过视察流场中的1点速度的时间波形加以得到。如果从丈量到的这个时间波形中肯定出具有指定的振幅和振幅几率，那末对各向同性湍流会得到1个高斯散布。实际上，仅从湍流的时间波形就能够视察到湍流速度平均值的存在，由于对湍流区内1给定点来讲，湍流的脉动时间会或多或少规则的重复出现情势相同本质不同的波形。另外，在1给定时刻，湍流的脉动对空间也会或多或少规则的重复出现情势相同本质不同的波形。因此，概括的说，在全部所考察的区域内湍流具有一样的整体结构。
比较不同的湍流运动，会发现它们的时间波形会有差别。这意味着，为了定量地描写湍流运动，必须引进湍流尺度概念，即时间尺度和空间尺度。这些尺度的大小由产生湍流的装置的大小和其中的速度肯定。例如，对管道出口流出的湍流，可以预期，它的时间尺度具有管理直径与平均流速之比的量级，它的空间尺度具有管道直径的量级。
第3章算例及结论
3.1. 算例
对运动声源的模型而言，我们可以把接收点的位置和声源之间的距离分成接收距离S、推延距离Sσ和有效距离Se的3种，以下面图中描写所示。对这3种不同的距离可作以下说明：
图3⑴⑴ 、和间的几何关系（飞机以速度V0前飞）
（1）1般位于飞行器上声源的所在位置1般常表示声波到达视察者看以看到的位置时，它和视察者也就是上图所描写的接收点间的距离。
（2）根据飞行速度和飞行方向，由第2章的公式可以推算出飞机上声源在它发射声波时刻的所在位置与接收点间的距离，为
，
（3）斟酌到声源运动时多普勒效应而至的声压增大作用，这相当于发射点与接收点间的距离的缩短，所以可以把它称为实际有效距离，由(2⑷⑴5)式，它可以描写为
式中为发射角。可以看出，当视察者接收到来自发射的声波时，声源已运动到，如图3⑴⑴所示，在声辐射方向的投影为
它正是多普勒效应所缩短的距离。
对这样的1个案例。首先，需要斟酌声源从发射点发射声波到到达接收点声波在传播路径上所花费的时间，由于两点间的距离为，所以相对接收点接收到声波的时间来讲，于发射点发射声波的时间为；其次，需要斟酌运动声源的多普勒效应，有声源辐射的声压增大。因此，对运动声源来讲，在接收点处于时间所接收到的声压瞬时值应当写成
假定声源是1架声压恒定为pA = 200Pa的飞机，那末视察者听见它声音的均方根值就能够写成
下面以两架声压恒定为pA = 200Pa（140dB）的相距200米以V0 = 260米/秒直线水平飞越视察者上空的飞机噪声为例，来讲明计算推延时间对噪声的预测的影响。
如图3⑴⑴所示，假定视察者为于处，试预测当第1架飞机正好在头顶时接收到的声音是多少分贝。实际上，当第1架飞机正好在头顶时，此时接收到它的声音应当是它位于点发射的，其发射角为；同理，当第1架飞机正好在头顶时，其后续200米处的第2架飞机位于，此时接收到第2架飞机的声音应当是第2架飞机位于点发射的，其发射角为。根据图3⑴⑴的几何关系可以建立以下方程：
图3⑴⑴
，，，
即
，，，
联立得
，
从而解得：米，米，米，米，度，度。由此得到接收者接收到的声压均方值（假定两架飞机发射的声波不相干）为
Pa
其相应的声压级为
dB
3.2. 结论
1．在接收点听到两架相距200米的飞机它们产生的声音其实不是它们在同1时刻产生的，如上例，第1架飞机是在米处发射的，而第2架飞机是在米处发射的，明显离接收点越远的那架飞机，其推延时间越长，即推延距离越长。
2．从上例中看出，第1架飞机在米处发射声波的多普勒效应为，而第2架飞机在米处发射声波的多普勒效应为，二者是不同的。
3．声传播推延时间的计算对正确预测运动声源的噪声具有特别重要的意义。
3.3. 章节总结
2014年5月上旬，环境保护部公布了《北京新机场项目环境影响报告书》，北京新机场将落地永定河北岸。北京市政协召开新机场建设提案办理协商会，北京市发改委相干负责人流露，新机场计划于2018年完成主体工程，2019年正式投入使用。建设北京第2机场应当说很合适京津空中运行特点，为空管运行提供了多种可行方案，最大限度地利用了北京地区紧张的空域资源，减少航空器地面滑行时间。
作为北京人，我深知当下北京生活环境的状态，雾霾的侵袭，春季的沙尘暴，包括汽车过量都为北京市民的生活质量在减分，而噪音污染一样难以忽视。北京首都国际机场建成于1958年，那时的顺义是绝对的郊区，所以机场所酿成的噪声对北京城区的影响寥寥无几。而如今几10年过去了，北京的主城区不知向外拓展了多少圈，曾今的鲜有人烟如今已经是住宅区密集，飞行器气动噪声的影响也在逐步增加。因此，在发展较落后的南城郊建设新机场是非常必要也是非常及时的。北京作为我们国家的首都，只是1个典型案例，现如今很多大城市都面临着城市拓展的问题，在人们愈来愈关注生活质量的今天，我想今天所研究的这个课题在未来的明天1定会有更大的实用价值。
致谢
大学这4年的时间真的过的很快。头几天晚上，我们几个室友熄灯后躺在床上聊天，其中1个说到，真不相信自己还有1个月就要毕业了，感觉从大学入学到现在简直就是1眨眼的工夫。确切，此时此刻我拿着行李箱初到南京从江宁校区大西门进入学校的情形记忆犹心，4年的时间如白驹过隙。在南京航空航天大学的这段岁月里，我非常高兴也非常荣幸的遇到了很多优秀的老师，结识了既是师长又是朋友的导员，更交到了很多可以推心置腹的好友。独在异乡为异客，出门到外求学的经历其实不是人人都有的，但我很高兴当我以后回想起这4年时光时，我会感到很充实，很成心义，很有味道。对我们这些学子来讲，6月绝对是1个收获的季节，而对我个人来讲，这篇论文就是我收获的果实。在我完成毕设，编写论文的进程中，我感觉我恍如把大学的这4年重新回顾了1遍，计算数据让我想起了大1时学高数这些基础课的窘境，绝对是当头1棒；画图时让我想起了学机械制图时画图纸的艰辛，深入体会工科生的不容易；解答案例时，让我想起了后来学专业课的刻苦研究。今天的果实离不开过去4年努力耕耘的结果，在论文完成之际，我在这里衷心的感谢在这几个月中给予过我帮助的学长，跟我1起研究讨论的同学，更要感谢我的导师张强老师，是您让我可以为大学生活和学习画上美满的句号。
如果说4年前初到南京航空航天大学时是对未来毫无掌控的未知，内心乃至有那末1点点不安；如今当我离开时绝对是心存感恩和对母校满怀不舍。不管将来在任什么时候间，任何地点，面对任何人，我想我都会自信的说到我的大学时光在南京航空航天大学度过！最后再次对这4年所有帮助过我的任课老师，同学，导员，学长，毕设导师，班主任说上1声谢谢！
参考文献：
[1] 张强等，气动声学基础，自编讲义，2010
[2] 姚起航，飞机噪声工程，西北工业大学出版社，1998
[3] Farassat F, Brown TJ. A new capability for predicting helicopter rotor and propeller noise including the effect of forward motion. NASA TM74037,June 1977.
[4] Farassat F. Acoustic radiation from rotating blades - The Kirchhoff method in aeroacoustics. Journal of Sound and Vibration Vol.239,No.4,2001
[5] Kenneth S. Brentner; F. Farassat. Modeling aerodynamically generated sound of helicopter rotors. Progress in Aerospace Sciences. vol.39; 2/3. 2003
[6] 姚起航，飞机噪声工程，西北工业大学出版社，1998
[7] Ffowcs Williams, J.E., and Hawkings, D.L., “Sound Generated by Turbulence and Surfaces in Arbitrary Motion,” Philosophical Transactions of the Royal Society, 1969, Vol.A264(1151):321⑶42
[8] Lighthill, M.J., “On Sound Generated Aerodynamically, I: General Theory,” Proceedings of the Royal Society, London, 1952, Vol.A211:564⑸87.