Linear and nonlinear approaches and digital signal processing in aeroacoustics

The linear and nonlinear approaches to the calculation of small acoustic disturbance propagation and evolution in nonuniform flows are compared. In the conventional linear approach it is the linearized equations of time-dependent, ideal (inviscid and non-heat-conducting) or viscous gas flow that are integrated. In the nonlinear approach the original nonlinear equations governing the same time-dependent flow (Euler equations for an ideal gas) are integrated; these are the same equations that, together with time relaxation procedure, are used in the linear approach for calculating the stationary background. It is shown that the application of digital signal processing, widely used in acoustic experiments, makes it possible to isolate the harmonic acoustic waves from the results of integration of the nonlinear equations, though their intensity is smaller than that of the noise due to computational errors, including inadequate attainment of the stationary background.     

A Navier-Stokes code for S-shaped diffusers—a review

Three-dimensional, compressible, internal flow solutions obtained using a thin-layer Navier-Stokes code are presented. The code, formulated by P.D. Thomas, is based on the Beam-Warming implicit factorization scheme; the boundary conditions also are formulated implicitly. Turbulent flow is treated through the use of the Baldwin-Lomax two-layer, algebraic eddy viscosity model. Steady-state solutions are obtained by solving numerically the time-dependent equations from given initial conditions until the time-dependent terms become negligible. The configuration considered is a rectangular cross-section, S-shaped centreline diffuser duct with an exit/inlet area ratio of 2.25. The Mach number at the duct entrance is 0.9, with a Reynolds number of 5.82 × 10⁵. Convergence to the final results required about 2700 time steps or 11 hours of CPU time on our CRAY-1M computer. The averaged residuals were reduced by about two orders of magnitude during the computations. Several regions of separated flow exist within the diffuser. The separated flow region on the upper wall, downstream of the second bend, is by far the largest and extends to the exit plane.     

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.     

Two-equation (k, ϵ) turbulence computations by the use of a finite element model

A finite element technique is presented and applied to some one- and two-dimensional turbulent flow problems. The basic equations are the Reynolds averaged momentum equations in conjunction with a two-equation (k, ?) turbulence model. The equations are written in time-dependent form and stationary problems are solved by a time iteration procedure. The advection parts of the equations are treated by the use of a method of characteristics, while the continuity requirement is satisfied by a penalty function approach. The general numerical formulation is based on Galerkin's method. Computational results are presented for one-dimensional steady-state and oscillatory channel flow problems and for steady-state flow over a two-dimensional backward-facing step.     

Liquid flow in a gap between a cylinder and a wall in motion

The results of an experimental investigation and calculations of the location of the minimum pressure point are presented for the case, when a cylindrical body moves along a wall in the presence of a small gap. The pressure on the cylindrical body surface is measured in the confusor and diffuser regions. It is shown that with decrease in the gap the minimum pressure point is displaced toward the minimum gap line, with increase in the pressure drop. An increase in the velocity of the motion at a constant gap leads only to a pressure increase in the diffuser region, while the location of the minimum pressure point remains the same. It is established that an increase in the inner cylinder radius moves the minimum pressure location away from the minimum gap line. The formation of two return flow regions in the confusor and diffuser regions near the cylindrical surface is detected. It is shown that the return flow in the pressure drop region reaches the stage of incipient cavitation bubbles. The results obtained can be useful in lubrication theory, as well as in medicine and biology.     

Implicit application of non‐reflective boundary conditions for Navier–Stokes equations in generalized coordinates

The non‐reflective boundary conditions (NRBC) for Navier–Stokes equations originally suggested by Poinsot and Lele (J. Comput. Phys. 1992; 101 :104–129) in Cartesian coordinates are extended to generalized coordinates. The characteristic form Navier–Stokes equations in conservative variables are given. In this characteristic‐based method, the NRBC is implicitly coupled with the Navier–Stokes flow solver and are solved simultaneously with the flow solver. The calculations are conducted for a subsonic vortex propagating flow and the steady and unsteady transonic inlet‐diffuser flows. The results indicate that the present method is accurate and robust, and the NRBC are essential for unsteady flow calculations. Copyright © 2005 John Wiley & Sons, Ltd.     

On Rayleight instability in decaying plane Couette flow

The inflexion point criterion of Rayleigh is one of the most well-known results in hydrodynamic stability theory but cannot easily be demonstrated experimentally in wall bounded flows. For plane Couette flow, where both walls move with equal speed in opposite directions, it is possible to establish a (time-dependent) inflectional velocity profile if both walls are brought momentarily to rest. If the Reynolds number is high enough a growing stationary instability develops. This situation is ideally suited for flow visualization of the instability. In this paper we show flow visualization experiments and stability calculations of the developing transverse roll cell instability in such a flow at low Reynolds numbers. Although the stability calculations are based on a quasi-stationary velocity profile, the measured and most amplified wave length obtained from the calculations are in excellent agreement.     

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

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

Turbulent flow in a circular separationless diffuser at Reynolds numbers smaller than 2000

The results of experimental and numerical investigation of flow in a circular conical diffuser with a small conicity angle ensuring separationless flow are presented. The measurements are carried out in an air flow with the Reynolds number Re₂ in the diffuser exit section ranging from 600 to 3000. A considerable effect of the channel expansion on the flow pattern is found to exist. It is shown that, as distinct from a tube, in which only laminar flow can be realized as steady for Re < 2000, in the exit section of a diffuser with the generator slope of 0.3° and a length equal to 70 entry diameters a developed turbulent flow is formed for Re₂ > 1000. For Re₂ > 1300 this flow is steady, that is, almost independent of the turbulence level at the entry, and is determined by the Reynolds number Re₂ in the exit section. For Re₂ ≈ 1000 the turbulent flow continuously goes over into a laminar flow. The flow parameters measured at the diffuser exit correspond to calculations in accordance with the threeequation turbulence model.     

无壁面参数低雷诺数非线性涡黏性模式研究

建立了一个低雷诺数的非线性涡黏性湍流模式,该模式的一个显著特性是它不包含壁面参数（如y^ ,n等）,因而特别适用于复杂几何流场的计算,本模式在几种包括回流、分离、激波等典型流动中进行了验证,结果令人满意。     

LES study on the influence of the diffuser inlet angle of a centrifugal pump on pressure fluctuations

Large-Eddy Simulations are conducted on a centrifugal pump at design and reduced flow-rates for three diffuser geometries, to investigate the effect of changing the diffuser inlet angle on the overall performance and the pressure fields. In particular, pressure fluctuations are investigated, which affect the unsteady loads acting on the pump, as well as vibrations, noise and cavitation phenomena. The considered modification of the diffuser geometry is targeted at decreasing the incidence angle at the off-design flow-rate by rotating the stationary blades of the pump around their leading edge. Results are compared against those of an earlier study, where the same modification of the diffuser inlet angle was achieved by increasing also the radial gap between impeller and diffuser, whose blades were rotated relative to their mid camber location. The comparisons across cases demonstrate that the radial gap between the trailing edge of the impeller blades and the leading edge of the diffuser blades has a more profound influence on pressure fluctuations, compared to the angle of incidence on the diffuser blades of the flow coming from the impeller.     

微喷管流动粒子仿真与Navier-Stokes数值模拟比较研究

运用DSMC（Direct Simulation Monte—Carlo）方法从分子运动论层次对大膨胀比、喉部转角为尖角的微喷管流动现象进行模拟,考察来流总压对喷管性能的影响,并与Navier—Stokes方程运算结果、实验结果进行比较。研究表明：在模拟微型喷管的流动现象时,DSMC方法比N—S方程更加适用。     

Experimental and computational study of the flow induced by a plasma actuator

A complementary experimental and computational study of the flow field evoked by a plasma actuator mounted on a flat plate was in focus of the present work. The main objective of the experimental investigation was the determination of the vector force imparted by the plasma actuator to the fluid flow. The force distribution was presently extracted from the Navier–Stokes equations directly by feeding them with the velocity field measured by a PIV technique. Assuming a steady-in-mean, two-dimensional flow with zero-pressure gradient, the imbalance between the convective term and the momentum equation’s right-hand-side terms reveals the desired resulting force. This force-distribution database was used afterwards as the source term in the momentum equation. Furthermore, an empirical model formulation for the volume-force determination parameterized by the underlying PIV-based model is derived. The model is tested within the RANS framework in order to predict a wall jet-like flow induced by a plasma actuator. The Reynolds equations are closed by a near-wall second-moment closure model based on the homogeneous dissipation rate of the kinetic energy of turbulence. The computationally obtained velocity field is analysed along with the experimental data focussing on the wall jet flow region in proximity of the plasma actuator. For comparison purposes, different existing phenomenological models were applied to evaluate the new model’s accuracy. The comparative analysis of all applied models demonstrates the strength of the new empirical model, particularly within the plasma domain. In addition, the presently formulated empirical model was applied to the flow in a three-dimensional diffuser whose inflow was modulated by a pair of streamwise vortices generated by the present plasma actuator. The direct comparison with existing experimental data of Grundmann et al. (2011) demonstrated that the specific decrease of the diffuser pressure corresponding to the continuous forcing was predicted correctly.     

Numerical simulation of swirling flow in diffusers

A reduced form of Navier–Stokes equations is developed which does not have the usual minimum axial step size restriction. The equations are able to predict accurately turbulent swirling flow in diffusers. An efficient single sweep implicit scheme is developed in conjunction with a variable grid size domain-conforming co-ordinate system. The present scheme indicates good agreement with experimental results for (1) turbulent pipe flow, (2) turbulent diffuser flow, (3) turbulent swirling diffuser flow. The strong coupling between the swirl and the axial velocity profiles outside of the boundary layer region is demonstrated.     

The structure of instantaneous reversals in highly turbulent flows

Structure of instantaneous flow reversals has been measured in a highly turbulent axisymmetric diffuser flow using pulsed-wire anemometry. In this 8° nominal included angle conical diffuser, the adverse pressure gradient (APG) is strong enough to cause appreciable instantaneous flow reversals (instantaneous backflow up to 30% of the time), but the time-averaged flow is non-separated. The results are compared with the other severe APG separating flows reported in literature. An increase in entry Reynolds number indicated a decrease in the size of near-wall instantaneous reversals region as well as a decrease in the magnitude of instantaneous backflow. Also, the region of instantaneous reversals moves slightly downstream at appreciably higher Reynolds numbers. The initiation and growth of instantaneous reversals in a conical diffuser was found to strongly influence the wall-layer and the central region. Present results also suggest that the instantaneous backflow should be considered for modelling of instantaneously-separating diffuser flows. In the final stages of a conical diffuser, the magnitudes of cross-stream pressure gradient were found to be appreciably larger than that of the longitudinal pressure gradient, indicating that accurate representation of a conical diffuser flow can not be achieved without considering V-momentum equation. A comparison of various separating flows revealed remarkable similarity of instantaneous reversals regions and distributions even in different flow configurations.     

NUMERICAL STUDY OF PERIODICAL FLOWS OF PIEZOELECTRIC VALVELESS MICROPUMP FOR BIOCHIPS

IntroductionMicro-pump is one of the most important fluid-driving elements in biochips. Workingprinciple of micropumps is quite different from traditional macro-scale pumps. Currently,valveless piezoelectric micropump is one of the most popular micropumps…     

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).     

Simplified system of equations for calculating the stability of the boundary layer gas flow on a swept wing

A composite system of equations for calculating the stability of a compressible flow in a three-dimensional boundary layer is developed. A comparison with the calculations in accordance with the complete system of equations shows that the system proposed almost accurately predicts the growth rates of both the Tollmien-Schlichting mode and the transverse flow instability mode over a wide range of frequencies and transverse wavenumbers at moderate Mach numbers. At the same time, the computation time reduces by an order.     

Flow characteristics in the downstream region of a conical diffuser

Using the Navier-Stokes equations in conjunction with the k-? model of turbulence, the characteristics of flow in the region downstream of a conical diffuser with 5° angle of inclination are calculated. Two representative stations 1D₂ and 10D₂ after the diffuser exit are selected for comparison against experimental results. The calculations indicate an underestimation of mean velocity and turbulence kinetic energy at the first station, while satisfactory agreement is obtained for the mean velocity at the second station. The use of a modified k-? model sensitive to adverse pressure conditions improves the predictions considerably. The effect of inlet properties and Reynolds number on the flow characteristics at the above stations is studied using various inlet profiles and a range of Reynolds numbers based on the inlet diameter from 50 000 to 280 000.