Numerical analysis of broadband noise reduction with wavy leading edge

Large Eddy Simulation(LES) is performed to investigate the airfoil broadband noise reduction with wavy leading edge under anisotropic incoming turbulence. The anisotropic incoming turbulence is generated by a rod with a diameter of 10 mm. The incoming flow velocity is 40 m/s and the corresponding Reynolds numbers based on airfoil chord and rod diameter are about 397000 and 26000, respectively. The far-field acoustic field is predicted using an acoustic analogy method which has been validated by the experiment. A straight leading edge airfoil and a wavy leading edge airfoil are simulated. The results show that wavy leading edge increases the airfoil lift and drag whereas the lift and drag fluctuations are substantially reduced. In addition, wavy leading edge can significantly change the flow pattern around the leading edge and a pair of counter-rotating streamwise vortices stemming from each wavy leading edge peak are observed.An averaged noise reduction of 9.5 dB is observed with the wavy leading edge at the azimuthal angle of 90°. Moreover, the wavy leading edge can mitigate noise radiation at all the azimuthal angles without significantly changing the noise directivity. The underlying noise reduction mechanisms are then analyzed in detail.     

增减材复合制造内部缺陷的涡流检测

涡流检测（ECT）技术具有非接触、无需耦合剂、检测灵敏等特点，适用于加工环境较为特殊的增减材复合制造（ASHM）中。本文建立了无缺陷半无限大试样内部涡流分布的解析模型，开展了预制人工缺陷的钛合金增材试样检测实验，研究了ECT深度与激励频率、提离量之间的关系。理论分析与实验结果均表明，内部缺陷较深时，低激励频率条件下缺陷产生的电抗增量信号较大，不同提离量下的电抗增量信号相差不大，因此检测位置较深的内部缺陷时可采用较低的激励频率并适当提高提离量。在本文实验条件下，ECT最佳激励频率为90 kHz；提离量增加到0.97 mm时，有效检测深度略有减小。这一结论可为ECT技术与ASHM的集成提供理论依据。     

多段翼型缝翼前缘结冰大迎角分离流动数值模拟

应用基于SST(Shear-Stress-Transport)湍流模型的IDDES(Improved Delayed Detached Eddy Simulation)方法,对大迎角状态下多段翼型缝翼前缘典型角状冰引起的复杂分离流动进行了数值模拟研究。采用后台阶流动标准算例和干净无冰多段翼型分离流动算例对数值方法的可靠性和适用性进行了验证。缝翼结冰状态下的数值模拟结果表明:来流迎角较大时,前缘角状冰将会导致结构相对稳定的流动分离泡产生,分离泡的非定常尾迹会对主翼前缘附近流场产生较为强烈的干扰,抑制了缝道流动的加速效应,使得缝翼增升效率降低。在失速点附近,由于分离泡回流强度随来流迎角而增长,同时脱落旋涡的输运方向逐渐向远离壁面方向偏移,使得尾迹影响区域范围和强度均有所增加。     

压气机叶栅端壁叶尖涡系结构非定常特性研究

为了研究压气机带间隙平面叶栅近失速工况下叶尖涡系结构特点以及其非定常流动特性,本文采用大涡模拟(Large eddy simulation,LES)方法对典型叶栅进行数值计算并结合Q准则分析叶尖涡系结构特点,探索流动规律及叶尖涡系耦合过程。研究表明:与额定工况相比,近失速工况叶尖流场更为复杂,并通过大涡模拟观察到了次泄漏涡的存在;额定工况下泄漏涡不发生破碎,主次泄漏涡在近失速条件下均发生破碎,破碎后形成的低能流体与尾缘分离涡是造成叶尖堵塞及损失的主要原因;次泄漏涡在不同时刻生成点位置及与弦长夹角周期性变化,次泄漏涡的摆动与叶尖角区分离涡团的周期性脱落是叶尖非定常性的主要原因。     

竖直平板Blasius流层流边界层流动与传热耦合解

对考虑辐射传热、流体黏度随温度变化、有和无滑移边界条件下的可渗透竖直平板Blasius流层流边界层的无量纲速度场与温度场进行了深入研究.经相似变换将描述速度场与温度场耦合的偏微分方程组转换成非线性常微分方程组,用Runge-Kutta法对常微分方程组进行了数值求解.研究了无量纲参数对无量纲速度场及温度场的影响,着重分析了滑移边界条件下速度和温度随无量纲参数的变化规律.结果表明:吸入时边界层变薄,喷注时边界层加厚;对比于无滑移边界条件,滑移边界条件下速度、温度边界层变薄;随着变黏度参数a或喷注与吸入参数的增大,壁面摩擦因数、局部努塞尔数Nu增大,速度和温度边界层变薄;随着普朗特数Pr、热辐射参数R的增加,或毕渥数Bi,布林克曼数Br的减小, 温度边界层变薄.      

基于超声速边界层转捩的WENO格式对比分析

基于高精度weighted essentially non-oscillation(WENO)格式对马赫数为4.5,雷诺数为10000的超声速来流条件下平板边界层转捩过程进行了大涡模拟.无黏通量分别用5阶、7阶、9阶WENO格式进行离散,黏性通量离散采用4阶中心差分格式,时间推进采用具有总变差递减(TVD)性质的3阶精度Runge-Kutta方法.通过在入口边界叠加一对等幅最不稳定第一模态谐波扰动,分别采用3种WENO格式计算得到了平板层流边界层失稳转捩的演化过程.结果表明:5阶WENO格式的数值耗散明显高于7阶和9阶WENO格式,在捕捉湍流涡和流场脉动特性上存在明显不足.9阶WENO格式的耗散小,能够捕捉到流场更小尺度的涡和高频脉动.研究具有高频脉动特性的问题或者欲捕捉精细涡结构时,建议采用7阶以上的高精度格式.      

汽车发电机冷却风扇旋转噪声预测方法

以汽车发电机冷却风扇为对象,针对其修改叶片分布角度前后的旋转噪声预测问题,提出一套较精确且节约计算机时的预测方法.该方法结合了声类比法和矢量合成方法,首先,用大涡模拟(LES)和Ffowcs-Williams和Hawkings(FW-H)方程相结合的声类比方法对原风扇总噪声和主要阶次旋转噪声幅值进行预测;然后,针对只改变叶片分布角度情况下,提出一种矢量合成方法,用于对修改叶片分布前后主要阶次旋转噪声变化量的预测;最后,得到修改后主要阶次旋转噪声的幅值.计算和实验结果表明,原风扇总噪声最大预测误差4.3dB,第12阶和第18阶主要阶次旋转噪声幅值预测误差为1.24dB和4.26dB;修改后风扇第12阶和第18阶主要阶次旋转噪声分别变小了9.3dB和10.5dB,其变化量预测误差分别为0.36dB和0.43dB.结果表明,这一整套对修改前后风扇旋转噪声进行预测的方法是可行的,且大大节省计算时间,为风扇叶片周向分布角度设计提供了很好的依据.      

整体次加筋壁板屈曲载荷近似计算方法

为了在初步设计阶段能够快速计算整体次加筋板的失稳载荷,在一些合理假设的基础上,提出了一种简单的近似计算方法.以无缺陷的四边简支的矩形次加筋板为研究对象,针对该结构的3种失稳形式,利用传统加筋板理论分别计算相应的屈曲载荷,并以3种失稳形式中最小的临界屈曲载荷作为整体次加筋板的近似屈曲载荷.应用ABAQUS软件的屈曲线性摄动步方法分别计算了两组有限元模型:一组用来验证3种失效形式理论公式计算的准确度;另一组是整体次加筋板有限元模型,用以验证所提出的次加筋板屈曲载荷计算方法的适用性.以上研究均考虑了纵向压缩载荷和压剪组合载荷两种工况.计算结果表明,理论近似计算方法能够准确地计算次加筋板的失稳载荷,有一定的工程应用价值.      

Construction of Crack Perturbation Model and Forward Semi-analytical Model of Attached Eddy Current Sensor

Fatigue damage monitoring is critical metallic structure health monitoring of aircraft.The sensor should be high sensitive,easy to be integrated into structure and well adaptable for poor working conditions.Therefore,an attached eddy current sensor with flexible plane is put forward and its characteristics are analyzed.By extracting material′s conductivity as the crack features,forward semi-analytical model is established and parameter optimizations are carried out.Crack perturbation model of attached eddy current sensor is constructed,and perturbation voltages of sensing channels under three-dimension structural crack are obtained.To verify the sensor′s performance,monitoring experiment on crack extension is conducted under condition of 3 MHz frequency.The validation experimental results show that perturbation model of 2A12-T4 aluminum alloy agrees well with experiment results,and perturbation model errors of four sensing channels are within 25%.The attached eddy current sensor is capable of testing the crack nondestructively and measuring the crack extension quantitatively with the accuracy of 1mm.     

一种确定求积公式余项的新方法

目前，各类文献中尚未见有确定一般形式数值求积公式余项的一般方法，正是基于此，本文提出了确定二类求积公式余项的一种新的通用的计算方法，得到了几个引理和定理，并且通过几个实例进行了说明。结果表明，该方法仅仅依赖于求积公式的代数精度，且计算简单。