φ型风力机空气动力学特性研究

本文采用均匀多流管模型和多流管模型,建立了基于Troposkien曲线的Φ型垂直轴风力机气动性能的计算方法和程序。在多流管模型中加入了附加阻力系数修正和动态失速修正,提高了计算精度。分析了流管数目对流管模型计算结果,以及实度和高径比对Φ型风力机空气动力学特性的影响。     

流管实验装置中声传播计算的模态方法

流管实验装置是测量有流动情况下航空发动机消声短舱内声衬声阻抗的主要装置。本文发展了一种解析的模态匹配方法进行在平均流有声衬条件下矩形流管中声传播的计算。用同伦方法求解特征值问题,并与用环绕积分求解的结果进行比较。声场通过轴向阻抗间断面的声压和声质点速度积分相等计算。第一个算例是无流动、硬壁、有限长、考虑端口反射的情况,并与北航流管实验台测量数据进行了对比;第二个算例为有流动情况下有限长声衬管道不考虑端口反射的声场计算,它与文献中NASA流管实验结果和CAA计算结果符合得很好。     

矩形表面波探头声场的高斯声束叠加法

利用矩形压电晶片和有机玻璃楔块折射可激励出超声表面波,广泛用于固体近表面缺陷检测和材料特性测量.由于描述表面波三维声场的理论方法还鲜有报道,因而主要采用简化的表面波二维声场模型来定量分析这类问题.高斯声束模型近些年被广泛应用于解决超声体波传播的各种复杂问题,然而,目前还没有将其扩展应用到超声表面波的声场的计算中.通过结合表面波格林方程和矩形换能器的高斯声束模型,推导出基于高斯声束叠加的表面波三维声场解析解.进一步,将该方法与点源叠加的数值解进行了分析比较,计算结果表明表面波声场的高斯声束叠加方法在具有较好计算精度的同时,还具有更快的计算效率.     

螺旋内肋管内对流换热的三维数值模拟

针对螺旋内肋管内壁面结构复杂的特点,发展了一种新的网格划分方法,采用结构化的六面体网格,提高计算精度的同时又可节省计算量。应用Fluent软件对螺旋内肋管内的湍流流动和换热进行了三维数值模拟,数值模拟结果与Jensen等人的实验数据吻合良好。在其他参数相同的条件下对矩形、三角形、半圆形三种顶端外形轮廓的肋片性能进行了数值分析比较。     

矩形微通道中环状冷凝的数值模拟

建立了一套恒热流边界条件下矩形微通道内环状冷凝过程的一维稳态模型,进行了数值计算,并将数值模拟结果与正三角形中的冷凝过程进行了比较.研究显示,在不同截面的微通道中,液相毛细半径和流速的沿程变化趋势都是相似的.同条件下,矩形通道的冷凝段长度大于正三角形通道.在矩形微通道中,通道入口蒸气压力和水力直径越大或者接触角越小,则冷凝段长度越长.     

矩形介质光波导的条形传递函数方法

将条形传递函数方法引入矩形介质光波导传播特性的分析计算,从标量波动方程的变分表达式出发,推导和建立了条形传递函数方法的基本方程和两种应用模型,作为算例,给出了阶跃型折射率分布矩形介质光波导的数值计算结果,并与其它方法进行了比较,表明了此方法的精度及对光波导问题的适用性。     

新型自由矩形螺旋线慢波结构高频特性

 在Sheath模型下,采用严格的场匹配法,结合积分形式的边界条件,推导了自由矩形螺旋线的色散方程和耦合阻抗表达式,并与近似理论进行对比。结果表明:与近似理论相比, 严禁场匹配法具有更高的准确性,且采用场匹配法的数值计算结果与3维商业电磁仿真软件结果吻合得很好。从而证明了所采用理论方法的有效性。同时分析了矩形螺旋线横截面尺寸、螺距、螺旋角、纵横比对色散特性和耦合阻抗的影响,结果表明：只有当矩形螺旋线横截面纵横比大于4时,才可忽略横截面的宽度对高频特性的影响,通过调节结构的参数可以改善色散和提高耦合阻抗。     

扁平多分支管内两相流流量分配的数值模拟

利用FLUENT软件,对一水平矩形扁平多分支管内,氮气和水为工质的气液两相流流量分配问题进行了数值模拟,并与实验数据进行了对比.数值模拟计算分析了相间曳力作用模型、气液进口条件、壁面粗糙度、离散相直径、结构参数对两相流流量分配的影响,指出为了提高模拟计算准确度,需要改进和完善现有模型.     

轴流式模型水轮机压力脉动试验与数值计算预测

本文将利用RNG κ-ε湍流模型对轴流式模型水轮机进行了全流道三维非定常湍流计算,预测模型水轮机的压力脉动性能.同时,本文对轴流式模型水轮机进行了压力脉动试验,通过对模型水轮机压力脉动试验结果与数值计算结果的对比,以验证数值计算预测压力脉动性能的准确性和可行性.     

弯曲微小通道流阻特性的数值模拟

用经典的N-S方程对流体在矩形截面、弯曲的微小通道中的流动特性进行了数值研究,发现计算结果与试验结 果存在较大的差异。在对流体流动特性分析的基础上引入了粗糙粘度模型来对经典的N-S方程进行修正,计算结果表明 用粗糙粘度模型计算的结果与试验值吻合较好。     

Side-branch resonators modelling with Green׳s function methods

This paper deals with strategies for computing efficiently the propagation of sound waves in ducts containing passive components. In many cases of practical interest, these components are acoustic cavities which are connected to the duct. Though standard Finite Element software could be used for the numerical prediction of sound transmission through such a system, the method is known to be extremely demanding, both in terms of data preparation and computation, especially in the mid-frequency range. To alleviate this, a numerical technique that exploits the benefit of the FEM and the BEM approach has been devised. First, a set of eigenmodes is computed in the cavity to produce a numerical impedance matrix connecting the pressure and the acoustic velocity on the duct wall interface. Then an integral representation for the acoustic pressure in the main duct is used. By choosing an appropriate Green?s function for the duct, the integration procedure is limited to the duct–cavity interface only. This allows an accurate computation of the scattering matrix of such an acoustic system with a numerical complexity that grows very mildly with the frequency. Typical applications involving Helmholtz and Herschel–Quincke resonators are presented.     

The attenuation of the higher-order cross-section modes in a duct with a thin porous layer

A numerical method for sound propagation of higher-order cross-sectional modes in a duct of arbitrary cross-section and boundary conditions with nonzero, complex acoustic admittance has been considered. This method assumes that the cross-section of the duct is uniform and that the duct is of a considerable length so that the longitudinal modes can be neglected. The problem is reduced to a two-dimensional (2D) finite element (FE) solution, from which a set of cross-sectional eigen-values and eigen-functions are determined. This result is used to obtain the modal frequencies, velocities and the attenuation coefficients. The 2D FE solution is then extended to three-dimensional via the normal mode decomposition technique. The numerical solution is validated against experimental data for sound propagation in a pipe with inner walls partially covered by coarse sand or granulated rubber. The values of the eigen-frequencies calculated from the proposed numerical model are validated against those predicted by the standard analytical solution for both a circular and rectangular pipe with rigid walls. It is shown that the considered numerical method is useful for predicting the sound pressure distribution, attenuation, and eigen-frequencies in a duct with acoustically nonrigid boundary conditions. The purpose of this work is to pave the way for the development of an efficient inverse problem solution for the remote characterization of the acoustic boundary conditions in natural and artificial waveguides.     

消声管道声学性能计算的二维等几何有限元方法及应用

将等几何有限元方法应用于消声管道的声学性能计算,使用二维等几何有限元方法求解管道截面的声学特征值,考虑了存在穿孔边界和吸声材料边界的情况,进而使用特征值计算消声管道的传递损失。对包覆式消声管道进行传递损失的计算,结果与二维有限元方法吻合较好.对圆形截面的特征值计算结果表明,在计算量相同的情况下,等几何有限元方法取得了比传统有限元方法更好的计算精度.在不同结构参数条件下对消声管道的声学性能进行计算,结果与三维数值方法吻合良好。方法能够在宽频范围内较好地预测消声管道的传递损失。      

Sound source in the presence of wall shear layers

Numerical solution techniques for evaluating the acoustic field generated by a single line source located inside or outside a wall shear layer of an infinitely long lined rectangular duct are presented. A formula for calculating wave attenuation due to an acoustic lining is given.     

等截面消声管道传递损失计算的简化方法

本文提出了等截面消声管道传递损失计算的简化方法，方法利用消声管道截面形式的特点，将三维声学计算问题简化为二维问题，消声管道的传递损失可以表示为与轴向波数有关的表达式，轴向波数可以通过计算截面的特征值得到。对于规则截面结构，使用传递矩阵法计算特征值；复杂非规则截面的特征值使用二维有限元方法得到，进而可以计算消声管道的传递损失。仿真结果与文献中的数值方法及实验值在较宽的频率范围内吻合较好，说明了方法的正确性，此外，该方法可以考虑均匀流对消声管道声学性能的影响。方法的计算效率高，对消声管道的前期优化设计具有实际意义。     

小阻尼封闭空间声压级灵敏度分析与优化设计

本文研究了小阻尼界面封闭空间低频声学有限元分析、灵敏度分析和优化设计问题。分别用模态法和直接法计算了封闭空间内声压级响应，并推导了声压级响应对声空间边界形状控制参数的灵敏度分析公式，在此基础上建立了小阻尼空间声学问题的优化模型，同时给出了优化求解算法，并在JIFEX软件中进行了程序实现。本文提出的灵敏度分析和优化设计方法可以使声场的边界布局更为合理，从而达到改进小阻尼界面封闭空间声学性能的目的。数值算例验证了本文提出的灵敏度分析和优化算法的有效性。     

Identification of a dynamical model for an acoustic enclosure using the wavelet transform

This paper presents results of research work in the identification of a dynamical model for an acoustic enclosure, a duct with rectangular cross-section, closed ends, and side-mounted speaker enclosures. An acoustic enclosure is excited randomly and random decrement functions are built to convert the random responses to free acoustic responses. It is shown that the estimation of resonance frequencies is possible using the wavelet transform of the system’s free response. Using a particular form of the son wavelet function, results are improved compared to those obtained with the traditionally Morlet wavelet function. An optimal value of a parameter of the son wavelet function is obtained by minimization of the wavelet entropy. The accuracy of this new technique is confirmed by applying it to a numerical example, and to an acoustic enclosure. The advantage of using the wavelet transform method over the Fourier-based modal analysis that would normally be used for the enclosure modes problem is established.     

Bulk reaction modeling of ducts with and without mean flow

A general formulation for analysis of sound field in a uniform flow duct lined with bulk-reacting sound-absorbing material is presented here. Presented theoretical model predicts the rate of attenuation for symmetric as well as asymmetric modes in rectangular duct lined with loosely bound (bulk-reacting) sound-absorbing material, which allows acoustic propagation through the lining. The nature of attenuation in rectangular ducts lined on two and four sides with and without mean flow is discussed. Computed results are compared with published theoretical and experimental results. The presented model can be used as guidelines for the acoustic design of silencers, air-conditioning ducts, industrial fans, and other similar applications.     

Finite difference computation of acoustic scattering by small surface inhomogeneities and discontinuities

The use of finite difference schemes to compute the scattering of acoustic waves by surfaces made up of different materials with sharp surface discontinuities at the joints would, invariably, result in the generations of spurious reflected waves of numerical origin. Spurious scattered waves are produced even if a high-order scheme capable of resolving and supporting the propagation of the incident wave is used. This problem is of practical importance in jet engine duct acoustic computation. In this work, the basic reason for the generation of spurious numerical waves is first examined. It is known that when the governing partial differential equations of acoustics are discretized, one should only use the long waves of the computational scheme to represent or simulate the physical waves. The short waves of the computational scheme have entirely different propagation characteristics. They are the spurious numerical waves. A method by which high wave number components (short waves) in the wave scattering process is intentionally removed so as to minimize the scattering of spurious numerical waves is proposed. This method is implemented in several examples from computational aeroacoustics to illustrate its effectiveness, accuracy and efficiency. This method is also employed to compute the scattering of acoustic waves by scatterers, such as rigid wall acoustic liner splices, with width smaller than the computational mesh size. Good results are obtained when comparing with computed results using much smaller mesh size. The method is further extended for applications to computations of acoustic wave reflection and scattering by very small surface inhomogeneities with simple geometries.