一种用边界元方法求解二维不可压缩粘性流动的方法

本文用边界元方法求解了二维不可压缩粘性流动的涡量——速度方程,利用求解区域边界上的速度法向导数和速度值直接得到了涡量的边界条件,克服了利用涡量方程求解二维不可压缩粘性流动时涡量边界条件(主要是壁面边界条件)难以给定的困难,算例表明:这种方法比较简单且计算结果基本上是令人满意的。     

离心压缩机蜗壳型线的一种二维逆命题计算方法

本文提出一个蜗壳型线的二维逆命题计算方法。该方法考虑了蜗壳进口流动参数沿圆周方向分布不均匀和蜗壳周向平面内损失分布的影响；改进蜗壳与叶轮或扩压器之间的匹配关系，改善离心压缩机和风机的性能，提高蜗壳的设计质量。     

任意回转面上,可压缩流动轴流式叶栅最佳壁面速度分布的计算方法

本文在边界层理论分析和叶栅损失计算的基础上,通过近代应用数学“最优控制论”的工具,首先提供了一个二元、不可压、稳定流动叶栅壁面最佳速度分布的理论计算方法。又通过一些转换的方法把它推广到了可压缩流动和任意回转面上。同时证明了二元不可压稳定流动叶栅壁面最佳速度分布经推广到可压缩流动任意回转面上去后,仍然是后者叶栅壁面的最佳速度分布。并且可以推广到二元扩压器,喷管以及轴对称细长旋成体表面最佳速度分布的计算。      

采用不同黏性处理方法的宽无叶扩压器不稳定流动研究

径向无叶扩压器的全局稳定性可能受到核心主流失稳,出口回流与壁面边界层分离等因素影响,对于宽无叶扩压器,无黏核心主流与壁面边界层流动对不稳定扰动诱发的作用机理是当前研究的重点.本文首先通过数值计算获得了大宽度比孤立无叶扩压器平均流动,然后基于小扰动理论和周向均质假设,分别对欧拉方程与 Navier-Stokes 方程进行线性化,建立了基于无黏核心流动的稳定性分析方法,以及基于涡黏性与分子黏性的混合稳定性分析方法;通过与实验结果的对比,验证了混合稳定性分析方法预测所得流动失稳频率和全局直接模态的准确性;最后基于伴随方法获得了特征值的结构敏感性,揭示了不同黏性处理条件下宽无叶扩压器内全局不失稳扰动的源发区域.在只考虑核心主流的无黏条件下,宽无叶扩压器内流动不稳定扰动来源于流场中部,为二维的离心失稳;在同时考虑核心主流与边界层的作用时,宽无叶扩压器不稳定扰动不仅来源于扩压器流场中部的核心主流,壁面回流对于不稳定扰动的产生了重要影响.      

采用不同黏性处理方法的宽无叶扩压器不稳定流动研究

径向无叶扩压器的全局稳定性可能受到核心主流失稳,出口回流与壁面边界层分离等因素影响,对于宽无叶扩压器,无黏核心主流与壁面边界层流动对不稳定扰动诱发的作用机理是当前研究的重点.本文首先通过数值计算获得了大宽度比孤立无叶扩压器平均流动,然后基于小扰动理论和周向均质假设,分别对欧拉方程与Navier-Stokes方程进行线性化,建立了基于无黏核心流动的稳定性分析方法,以及基于涡黏性与分子黏性的混合稳定性分析方法;通过与实验结果的对比,验证了混合稳定性分析方法预测所得流动失稳频率和全局直接模态的准确性;最后基于伴随方法获得了特征值的结构敏感性,揭示了不同黏性处理条件下宽无叶扩压器内全局不失稳扰动的源发区域.在只考虑核心主流的无黏条件下,宽无叶扩压器内流动不稳定扰动来源于流场中部,为二维的离心失稳;在同时考虑核心主流与边界层的作用时,宽无叶扩压器不稳定扰动不仅来源于扩压器流场中部的核心主流,壁面回流对于不稳定扰动的产生了重要影响.     

受限亚音速气流中倒置悬臂壁板静气弹稳定性的理论及实验研究

板壳结构在航空航天、高速列车、能量采集等诸多工程领域已经得到了广泛应用.将悬臂壁板倒置于轴向气流中并在壁板周围流场中设置刚性壁面可有效地调控壁板的失稳速度,是俘能器优化设计的重要措施之一.但针对刚性壁面作用下亚音速气流中倒置悬臂壁板的失稳机制仍需要开展深入研究.本文以受限亚音速气流中倒置的二维悬臂壁板为对象,以理论分析及风洞实验为手段,研究了单侧刚性壁面效应对倒置悬臂壁板静态失稳特性的影响规律.在理论分析中,首先应用镜像函数法来处理壁面约束条件,基于算子理论研究获得了以Possio积分方程为表征的壁板气动力,壁面效应实际表征为一包含移位Tricomi算子的复合算子;然后将壁板失稳方程的求解问题转化为定区间上的函数逼近问题;最后,依据Wererstrass定理并利用最小二乘法求解该最优函数,以获得系统的失稳临界参数.在试验研究中依据压杆稳定原理设计了壁板静态失稳的测试方法并完成了风洞实验.理论分析结果表明,壁板会发生发散(静气动弹性)失稳,临界动压随壁板与壁面间距的增加而增大并最终趋于稳定(无壁面情况);通过理论与风洞实验结果的对比分析,验证了本文气动力及理论分析的适用性及准确性.针对倒置悬臂壁板结构的气动弹性失稳问题,本文提出的方法不涉及系统方程的离散及特征值求解问题,而是将其转化为了定区间上的函数逼近问题进行求解,这为弹性结构静气动弹性失稳问题的研究提供了一个可行的新思路.     

受限亚音速气流中倒置悬臂壁板静气弹稳定性的理论及实验研究

板壳结构在航空航天、高速列车、能量采集等诸多工程领域已经得到了广泛应用. 将悬臂壁板倒置于轴向气流中并在壁板周围流场中设置刚性壁面可有效地调控壁板的失稳速度, 是俘能器优化设计的重要措施之一. 但针对刚性壁面作用下亚音速气流中倒置悬臂壁板的失稳机制仍需要开展深入研究. 本文以受限亚音速气流中倒置的二维悬臂壁板为对象, 以理论分析及风洞实验为手段, 研究了单侧刚性壁面效应对倒置悬臂壁板静态失稳特性的影响规律. 在理论分析中, 首先应用镜像函数法来处理壁面约束条件, 基于算子理论研究获得了以Possio积分方程为表征的壁板气动力, 壁面效应实际表征为一包含移位Tricomi算子的复合算子; 然后将壁板失稳方程的求解问题转化为定区间上的函数逼近问题; 最后, 依据Wererstrass定理并利用最小二乘法求解该最优函数, 以获得系统的失稳临界参数. 在试验研究中依据压杆稳定原理设计了壁板静态失稳的测试方法并完成了风洞实验. 理论分析结果表明, 壁板会发生发散(静气动弹性)失稳, 临界动压随壁板与壁面间距的增加而增大并最终趋于稳定(无壁面情况); 通过理论与风洞实验结果的对比分析, 验证了本文气动力及理论分析的适用性及准确性. 针对倒置悬臂壁板结构的气动弹性失稳问题, 本文提出的方法不涉及系统方程的离散及特征值求解问题, 而是将其转化为了定区间上的函数逼近问题进行求解, 这为弹性结构静气动弹性失稳问题的研究提供了一个可行的新思路.      

???????ζ?ALE??????????????

利用自由液面曲率参量，在传统求解方法基础上对相应结点网格速度引入修正，以提高计算精度.域内网格速度计算，由传统线性关系改用分段函数表达，以保证网格剖分的质量. 将二维矩形贮箱内流体晃动问题作为算例，求解激励响应结果并与已有物理实验进行对比,二者符合较好.     

绕壁面小障碍物的湍流边界层研究

采用RNGK－ε模型，对统平板壁面上的条柱、方柱等多种不同形状二维小障碍物的固壁剪切湍流进行数值研究，获得了小障碍物下游回流区及其发展区的时均速度场、回流分离点长度和湍流边界层厚度，计算结果与实验比较吻合良好．结果表明，小障碍物在壁面上的设置，明显地改变了平板壁面边界层的湍流形态和结构．     

亚临界雷诺数二维圆柱绕流边界层流场动态特性的实验研究

本文给出了关于亚临界雷诺数二维圆柱绕流动态特性的实验结果。应用热膜、热线和压力传感器测量了壁面剪切应力脉动、壁面压力脉动和流场的速度脉动,给出了壁面剪切应力脉动频率在驻点附近和分离前后变化的特征,计算了这些脉动量在圆柱面上任意两点间的相关特性。实验结果表明,在亚临界雷诺数二维圆柱绕流的边界层流场中存在着一个频率与涡脱落频率相同的整体同步脉动。     

Heat Transfer Influenced by Turbulent Airflow Inside an Axially Rotating Diffuser

The paper presents a study of heat transfer between the turbulent airflow and the inner wall surface of an axial diffuser rotating around its longitudinal axis. Heat transfer was assessed through the measurement of a time-dependent temperature field of the diffuser inner wall surface. Measurements of the instantaneous flow velocity components were performed by a laser–Doppler anemometry system, which delivered information on mean velocity components as well as on the turbulence intensity. A significant increase of all three mean velocity components was observed near the rotating diffuser wall in comparison with a non-rotating diffuser. Temperature field measurements were carried out by means of infrared thermography. The experiment showed a significant dependence of the temperature field on the turbulent flowfield induced by diffuser rotation. A strong influence of the flow separation and reattachment on the temperature distribution was observed, while rotation was found to suppress the occurrence of flow separation from the diffuser wall. Properties of the velocity field such as turbulent kinetic energy were directly coupled with the temperature distribution in order to gain the information on how to enhance or reduce heat transfer by changing the integral parameters of the diffuser (e.g. rotation frequency or amount of flow).     

Effect of humidity on the flow characteristics of a two-phase medium in a diffuser

The effect of the initial humidity of a two-phase flow on the aerodynamic characteristics of a diffuser is considered. A series of calculations is carried out which permit the pressure distribution in the diffuser to be investigated in the case of change of the initial flow velocity, coefficient of slip of the phases, and the initial humidity of the flow. It is shown that, with increase of the initial humidity in the diffuser, there appears a sharply expressed region which has a negative pressure gradient, where the steam velocity can reach the critical value. The numerical solution is compared with the experimental results.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 6, pp. 44–49, November–December, 1973.     

The use of perforated plates to control the flow emerging from a wide-angle diffuser,with application to electrostatic precipitator design

The flow through a wide-angle, pyramidal diffuser of area ratio 6.8, in which two perforated plates are used to control the emergent velocity distribution was investigated. (Wide-angle diffusers combined with perforated plates find application in electrostatic precipitator flow systems. The efficacy of these systems depends critically on the uniformity of the gas flow within the collection chamber downstream of the diffuser outlet plane.) The main results of the study are (i) the establishment of the main factor influencing the flow at the diffuser exit plane, (ii) the determination of plate characteristics which produce a uniform velocity profile in the collection chamber, and (iii) the establishment of the corresponding pressure drop characteristics of the plate—diffuser combinations. The results also extend the range of diffuser geometries for which two perforated plates provide uniform velocity profiles at exit.     

Heterodyne Doppler global velocimetry

Doppler Global Velocimetry (DGV) is an imaging flow measurement technique which allows the measurement of the velocity distribution in a plane. In DGV the frequency shift of scattered light from moving particles within the flow is used to determine the local flow velocity. Heterodyne Doppler Global Velocimetry (HDGV) is a new approach which combines the imaging and geometrical characteristics of DGV with the measurement principles of reference beam laser Doppler anemometry. The frequency shifted scattered light from the flow tracers is heterodyned with a reference beam from the same light source. Due to interference the result of this superposition is a harmonic intensity modulated signal. This signal is detected using a smart pixel detector array to obtain the velocity distribution. Two different experiments are presented. The first experiment compares the measured velocity distribution of a rotating disk with its actual velocity. The second experiment demonstrates the capability of the technique to measure a real flow.     

Turbulent flow in a conical diffuser: A review

This is a review of experimental studies of turbulent flow in a conical diffuser by eight Ph.D. students, eleven M.Sc. students, one M.Eng. student, and myself in the past 29 years. During this time, two conical diffusers were constructed: the first was of cast aluminum construction, and the second was of plastic fabrication. These two diffusers were basically the same in geometry except that the pipe section was constructed as an integral part of the plastic diffuser to avoid the lip at the junction of the inlet pipe and the diffuser. The conical diffuser had a total divergence angle of 8°, an area ratio of 4:1, and an inlet diameter of 0.1016 m (4 in.).

The flow at the inlet of the diffuser was usually fully developed pipe flow, but sometimes it was boundary layer grown on the pipe wall. Hot-wire and pulse-wire anemometry together with computer facilities were used to obtain the results of complex flow present in the conical diffuser. Mean velocity profiles were obtained throughout the diffuser, which in turn were used to obtain strain rates and their principal direction. Turbulence moments up to fourth order were measured. The results were used to assess momentum, turbulent kinetic energy, and shear stress equations. Other features such as instantaneous flow reversals in the wall region, relative strength of large eddies, extra strain rate, and the production of kinetic energy also were investigated to find the dynamical picture in the diffuser flow.     

Turbulent separated reattached flow in a two-dimensional curved-wall diffuser

A turbulent separation-reattachment flow in a two-dimensional asymmetrical curved-wall diffuser is studied by a two-dimensional laser doppler velocimeter. The turbulent boundary layer separates on the lower curved wall under strong pressure gradient and then reattaches on a parallel channel. At the inlet of the diffuser, Reynolds number based on the diffuser height is 1.2×10⁵ and the velocity is 25.2m/s. The results of experiments are presented and analyzed in new defined streamline-aligned coordinates. The experiment shows that after Transitory Detachment Reynolds shear stress is negative in the near-wall backflow region. Their characteristics are approximately the same as in simple turbulent shear layers near the maximum Reynolds shear stress. A scale is formed using the maximum Reynolds shear stresses. It is found that a Reynolds shear stress similarity exists from separation to reattachment and the Schofield-Perry velocity law exists in the forward shear flow. Both profiles are used in the experimental work that leads to the design of a new eddy-viscosity model. The length scale is taken from that developed by Schofield and Perry. The composite velocity scale is formed by the maximum Reynolds shear stress and the Schofield-Perry velocity scale as well as the edge velocity of the boundary layer. The results of these experiments are presented in this paper.     

Minimization of the Wing Airfoil Drag Coefficient Using the Optimal Control Method

The problem of designing a wing airfoil starting from the surface velocity (or pressure) distribution, given in multi-parameter form, is considered. In the diffuser region, the velocity decrease law is determined from the conditions of minimum drag and separationless flow for a given Reynolds number. The case of boundary layer suction for the purpose of improving the aerodynamic characteristics of the airfoil is studied. A solution is obtained using optimal control theory.     

Sensitivity of an asymmetric 3D diffuser to plasma-actuator induced inlet condition perturbations

Experiments were conducted for the flow in a straight-walled 3D diffuser fed by a fully developed turbulent duct flow. Previous work found that this diffuser has a stable 3D separation bubble whose configuration is affected by the secondary flows in the upstream duct. Dielectric barrier discharge plasma actuators were used to produce low-momentum wall jets to determine if the separation behavior could be modified by weak forcing. Actuators producing a streamwise force along the wall where separation occurred in the baseline flow had a relatively small effect. However, spanwise acting plasma actuators that produced a pair of streamwise vortices in the inlet section of the diffuser had a strong effect on the diffuser pressure recovery. The diffuser performance could be either improved or degraded depending on the actuation parameters, including the actuator modulation frequency, duty cycle, and drive voltage. Velocity profile measurements in the diffuser inlet showed that the streamwise vortices affect the uniformity of the streamwise mean velocity accounting for some of the performance changes. However, phase-locked hotwire measurements at the diffuser exit indicate that the periodic nature of the forcing also plays an important role for cases with enhanced pressure recovery.     

Numerical prediction of turbulent flow in a conical diffuser usingk-ε model

Fully developed incompressible turbulent flow in a conical diffuser having a total divergence angle of 8° and an area ratio of 4∶1 has been simulated by ak-ε turbulence model with high Reynolds number and adverse pressure gradient. The research has been done for pipe entry Reynolds numbers of 1.16×10⁵ and 2.93×10⁵. The mean flow velocity and turbulence energy are predicted successfully and the advantage of Boundary Fit Coordinates approach is discussed. Furthermore, thek-ε turbulence model is applied to a flow in a conical diffuser having a total divergence angle of 30° with a perforated screen. A simplified mathematical model, where only the pressure drop is considered, has been used for describing the effect of the perforated screen. The optimum combination of the resistance coefficient and the location of the perforated screen is predicted for high diffuser efficiency or the uniform velocity distribution.     

Separation control in a conical diffuser with an annular inlet: center body wake separation

In many practical applications of conical diffusers, the flow is fed by an annular flow passage formed by a center body. Flow separation, which occurs if the center body ends abruptly, is undesirable because it degrades the diffuser performance. The present experiment utilizes magnetic resonance velocimetry to acquire three-component mean velocity measurements for a set of conical diffusers with an annular inlet. The results show strong coupling between the diffuser wall boundary layer development and the wake of the center body. Coanda blowing is used to mitigate the center body wake separation. The diffuser wall boundary layer is thick in the absence of the central separation bubble and separates when Coanda blowing is too strong.