汽车涡轮增压器同步谐波噪声仿真与优化

针对增压发动机急加速急减速时产生的增压器同步谐波噪声问题,该文通过噪声仿真技术进行分析与优化.首先,分析该噪声的特征与传播路径;其次,建立流场仿真模型.利用剪切应力输运湍流模型与分离涡流模拟湍流模型对增压器进行稳态与非稳态瞬态流场分析,提取非稳态流场的叶轮与压气机流道表面的偶极子声源;最后,建立噪声传播模型,计算该增压...     

Large-eddy simulations of temporally accelerating turbulent channel flow

The unsteady turbulent channel flow subject to the temporal acceleration is considered in this study. Large-eddy simulations were performed to study the response of the turbulent flow to the temporal acceleration. The simulations were started with the fully developed turbulent channel flow at an initial Reynolds number of Re₀ = 3500 (based on the channel half-height and the bulk-mean velocity), and then a constant temporal acceleration was applied. During the acceleration, the Reynolds number of the channel flow increased linearly from the initial Reynolds number to the final Reynolds number of Re₁ = 22,600. The effect of grid resolution, domain size, time step size on the simulation results was assessed in a preliminary study using simulations of the accelerating turbulent flow as well as simulations of the steady turbulent channel flow at various Reynolds numbers. Simulation parameters were carefully chosen from the preliminary study to ascertain the accuracy of the simulation. From the accelerating turbulent flow simulations, the delays in the response of various flow properties to the temporal acceleration were measured. The distinctive features of the delays responsible for turbulence production, energy redistribution, and radial propagation were identified. Detailed turbulence statistics including the wall shear stress response during the acceleration were examined. The results reveal the changes in the near-wall structures during the acceleration. A self-sustaining mechanism of turbulence is proposed to explain the response of the turbulent flow to the temporal acceleration. Although the overall flow characteristics are similar between the channel and pipe flows, some differences were observed between the two flows.     

湍流扩散火焰的非稳态火焰面模拟

采用稳态的和非稳态的火焰面模型同时对一个湍流甲烷射流扩散火焰进行了数值模拟,比较了两者对湍流平均火焰结构、活性自由基和污染物(氮氧化物)排放的模拟效果。速度场采用κ-ε模型计算,守恒标量混合物分数的分布通过其概率密度函数(PDF)输运方程的求解得到。稳态的火焰面结构由查询火焰面数据库得到,而非稳态的火焰面结构由火焰面方程和流场方程耦合求解来计算。采用详细的GRI—Mech 3.0机理描述甲烷的氧化和氮氧化物的形成。数值模拟结果和实验数据作了广泛的对比,验证了火焰面模型对湍流扩散燃烧的定量模拟能力。     

Numerical simulation of turbulence transition and sound radiation for flow through a rigid glottal model

Large eddy simulation (LES)-based computational aeroacoustics techniques were applied to a static model of the human glottis, idealized here as a planar channel with an orifice, to study flow-acoustic interactions related to speech. Rigid models of both converging and diverging glottal passages, each featuring a 20 deg included angle and a minimal glottal diameter of 0.04 cm, with an imposed transglottal pressure of 15 cm H2O, were studied. The Favre-filtered compressible Navier-Stokes equations were integrated for this low-Mach-number flow using an additive semi-implicit Runge-Kutta method and a high-order compact finite-difference scheme with characteristic-based nonreflecting boundary conditions and a multiblock approach. Flow asymmetries related to the Coanda effect and transition to turbulence, as well as the far-field sound, were captured. Acoustic-analogy-based far-field sound predictions were compared with direct simulations and showed that dipole sources, arising from unsteady flow forces exerted on the glottal walls, are primarily responsible for the tonal sound observed in the divergent glottis case.     

对旋轴流水轮机的非定常数值模拟研究

为了研究采用对旋转轮的新型轴流式水轮机的非定常性能,应用三维非定常数值计算方法,对该新型水轮机进行了数值模拟研究.数值计算中,采用了全三维全流道的数值计算方法,基于DES(Detached Eddy Simulation)方法获得了对旋轴流水轮机的非定常流场和能量性能.计算结果显示,该新型水轮机在一个转动周期的不同时刻,其内部流场存在显著的非定常特性.     

Numerical simulation of a turbulent magnetic dynamo

We present numerical simulations of a turbulent magnetic dynamo mimicking closely the Riga-dynamo experiment at Re approximately 3.5x10(6) and 15< or =Rem< or =20. The Reynolds-averaged Navier-Stokes equations for the fluid flow and turbulence field are solved simultaneously with the direct numerical solution of the magnetic field equations. The fully integrated two-way-coupled simulations reproduced all features of the magnetic self-excitation detected by the Riga experiment, with frequencies and amplitudes of the self-generated magnetic field in good agreement with the experimental records, and provided full insight into the unsteady magnetic and velocity fields and the mechanisms of the dynamo action.     

Experimental characterization of the convective heat transfer in a vortex-wall interaction

The development of turbulence models and wall laws for the numerical simulation of flows in complex geometries requires a detailed experimental analysis of turbulence and of the phenomena that appear in turbulent boundary layers. There is a strong need to develop new measurement systems allowing the determination of unsteady wall heat transfer coefficients. In order to improve the knowledge of the unsteady phenomena occurring in perturbed boundary layers, a fundamental study is conducted on the interaction of a single vortex with a flat plate. An experimental methodology using a specific thermal sensor whose surface temperature is measured by an infrared thermography system is presented. It allows the characterization of the unsteady convective heat transfer coefficient whose evolution is compared with the fluctuations of the wall friction coefficient, calculated from velocity profiles measured by laser Doppler velocimetry.     

三维边界层内诱导横流失稳模态的感受性机理

边界层感受性问题是层流向湍流转捩的初始阶段,在转捩过程中起关键性作用,尤其是三维边界层流动.因此,研究三维边界层感受性问题对进一步理解层流向湍流转捩机理以及湍流成因具有重要的理论意义.采用数值方法研究自由来流湍流与三维壁面局部粗糙相互作用下三维边界层的感受性问题,确定是否能在三维边界层内寻找一种新的横流失稳模态;确定在何种条件下三维边界层内能诱导出定常、非定常的横流失稳模态;探索自由来流湍流的强度、展向波数和法向波数以及三维壁面局部粗糙的大小和结构类型等因素在自由来流湍流与三维壁面局部粗糙作用下三维边界层内被激发出的感受性过程中有何影响,并确定何种横流失稳模态在三维边界层感受性过程中占据何种地位.对自由来流湍流与三维壁面局部粗糙作用激发三维边界层内感受性问题的深入研究,将有助于完善流动稳定性与湍流理论,为层流向湍流转捩过程的预测与控制提供合理的理论依据.     

Adaptation of mesoscale turbulence parameterisation schemes as RANS closures for ABL simulations

The present study presents different k-ε turbulence closures for atmospheric boundary layer flows using computational fluid dynamics (CFD) simulations that are consistent with inflow conditions from numerical weather prediction (NWP) simulations. Eight different mesoscale turbulence parameterisation schemes of the Weather Research and Forecasting (WRF) model are covered. To ensure consistency between the NWP and CFD simulations, different closure coefficients of the k ? ε turbulence model for each NWP scheme are proposed. This is achieved by combining production–dissipation closure coefficient relationships based on the Monin–Obukhov similarity theory and the formulation based on the Coriolis parameter proposed by Detering and Etling. The proposed methodology has been implemented in the open source CFD toolbox OpenFOAM and is demonstrated at near-neutral stability conditions for the classical Askervein Hill case.     

叶栅内定常及非定常粘性流动数值模拟

本文给出了一个模拟叶栅内准三维定常和非定常粘性流动的数值方法。对于定常流动,采用ＴＶＤ Ｌａｘ－Ｗｅｎｄｒｏｆｆ格式和代数湍流模型求解雷诺平均Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程,使用当地时间步长和多网格技术使计算加速收敛到定常状态;对于非定常流动,使用双时间步长和全隐式离散,采用与求解定常流动相似的多网格方法求解隐式离散方程。文中给出了ＶＫＩ透平叶栅内的定常流结果和１．５级透平叶栅内的非定常数值结果。     

前缘曲率对三维边界层内被激发出非定常横流模态的影响研究

三维边界层感受性问题是三维边界层层流向湍流转捩的初始阶段,是实现三维边界层转捩预测与控制的关键环节.在高湍流度的环境下,非定常横流模态的失稳是导致三维边界层流动转捩的主要原因;但是,前缘曲率对三维边界层感受性机制作用的研究也是十分重要的课题之一.因此,本文采用直接数值模拟方法研究在自由来流湍流作用下具有不同椭圆形前缘三维（后掠翼平板）边界层内被激发出非定常横流模态的感受性机制;揭示不同椭圆形前缘曲率对三维边界层内被激发出非定常横流模态的扰动波波包传播速度、传播方向、分布规律、感受性系数以及分别提取获得一组扰动波的幅值、色散关系和增长率等关键因素的影响;建立在不同椭圆形前缘曲率情况下,三维边界层内被激发出非定常横流模态的感受性问题与自由来流湍流的强度和运动方向变化之间的内在联系;详细分析了不同强度各向异性的自由来流湍流在激发三维边界层感受性机制的物理过程中起着何种作用等.通过上述研究将有益于拓展和完善流动稳定性理论,为三维边界层内层流向湍流转捩的预测与控制提供依据.     

U-RANS／PDF方法计算钝体后漩涡脱落

本研究发展了U-RANS/PDF混合算法研究湍流和化学反应相互作用对燃烧稳定性的影响,采用有限体积和Monte Carlo相结合的方法在非结构网格中求解相容的U-RANS方程和脉动速度-湍流频率-标量的联合PDF方程.本文对钝体火焰驻定器后冷态流场进行了计算,结果表明此混合算法能够捕捉流场中非稳态的漩涡脱落现象.着重研究了湍流频率模型系数的改变对漩涡脱落频率以及拟序结构在动量输运中的作用的影响.     

A comparative study of turbulence models for non-premixed swirl-stabilized flames

ABSTRACT

The accuracy of turbulent swirl-stabilized flame simulation strongly depends on the choice of turbulence model. In this study, four 3D unsteady turbulence closures, including large eddy simulation, scale-adaptive simulation, and two detached eddy simulation variants, along with four RANS models, including RNG k??, SST k?ω, transition SST, and RSM, are examined for moderate- and high-swirl case studies. It is observed that the scale-adaptive simulation provides the most accurate results for almost all variables and both swirl conditions in the reactive flow. Only the 3D unsteady models predict the vortex breakdown bubble and flame attachment state correctly. However, based on our error analysis, the flow and composition fields predicted by the RANS models are in acceptable agreement with the experimental fields, especially the ones of transition SST when higher swirl number cases or minor species concentration are of interest. Moreover, it is concluded that the viscosity ratio criterion is a better measure of the local LES quality than the turbulent kinetic energy ratio, and the accuracy of a hybrid simulation may be much more dependent on the ability of the model to operate close to the RANS mode where the grid resolution is not sufficient for a resolving simulation than the fraction of the resolved kinetic energy. Finally, the propriety of the base (RANS) model of a DES for the application of interest is important, such that DES with realizable k?? outperforms the commonly used DES with SST k?ω model.     

Large-eddy simulation of impinging jets with embedded azimuthal vortices

We have carried out large-eddy simulations of an impinging jet with embedded azimuthal vortices, a model of the wake of a helicopter hovering in ground effect. The azimuthal vortices are generated by sinusoidal forcing of the velocity at the jet exit. They strengthen while they are advected towards the ground; when they are close to the solid surface, a layer of opposite-sign vorticity is formed at the wall, and lifted up to form a secondary vortex that interacts with the primary one. Regions of reversed flow are caused by the strong, localised, adverse pressure gradient. After this interaction, the primary vortices begin to decay, mostly due to the Reynolds shear stresses, which contribute to the turbulent diffusion of vorticity term in the budget of the phase-averaged azimuthal vorticity. This mechanism is extremely robust, and plays the most important role in the vortex decay even if no turbulence is initially present in the jet, or if the no-slip condition is removed. A three-dimensional instability also plays a role: removing it leads to slower decay. Our results also point out some challenges for turbulence models for the unsteady Reynolds-averaged Navier–Stokes equations.     

Coupling of turbulence on the ignition of multicomponent sprays

This study examines the effect of turbulence on the ignition of multicomponent surrogate fuels and its role in modifying preferential evaporation in multiphase turbulent spray environments. To this end, two zero-dimensional droplet models are considered that are representative of asymptotic conditions of diffusion limit and the distillation limit are considered. The coupling between diffusion, evaporation and combustion is first identified using a scale analysis of 0D homogeneous batch reactor simulations. Subsequently, direct numerical simulations of homogeneously dispersed multicomponent droplets are performed for both droplet models, in decaying isotropic turbulence and at quiescent conditions to examine competing time scale effects arising from evaporation, ignition and turbulence. Results related to intra-droplet transport and effects of turbulence on autoignition and overall combustion are studied using an aviation fuel surrogate. Depending on the characteristic scale, it is shown that turbulence can couple through modulation of evaporation time or defer the ignition phase as a result of droplet cooling or gas-phase homogenization. Both preferential evaporation and turbulence are found to modify the ignition delay time, up to a factor of two. More importantly, identical droplet ignition behavior in homogeneous gas phase can imply fundamentally different combustion modes in heterogeneous environments.     

Dynamics of passive-scalar turbulence

We present the first study of the dynamic scaling or multiscaling of passive-scalar turbulence. For the Kraichnan version of passive-scalar turbulence we show analytically, in both Eulerian and quasi-Lagrangian frameworks, that simple dynamic scaling is obtained but with different dynamic exponents. By developing the multifractal model we show that dynamic multiscaling occurs in passive-scalar turbulence only if the advecting velocity field is itself multifractal. We substantiate our results by detailed numerical simulations in shell models of passive-scalar advection.     

Qualitative comparison between numerical and experimental results of unsteady flow in a radial diffuser pump

Comparison between numerical simulation and experimental results for unsteady flow field in a radial diffuser pump is presented for the design operating point. The numerical result is obtained by solving three-dimensional, unsteady Reynolds-averaged Navier-Stokes equations by the commercial CFD code CFX-10 withk-ω based shear stress transport turbulence model. Two-dimensional PIV measurements are conducted to acquire the experiment result. The phase-averaged velocity and turbulent kinetic energy fields are compared in detail between the results by the two methods in the impeller, diffuser and return channel regions. The qualitative comparison between CFD and PIV results is quite good in the phase-averaged velocity field. Although the turbulence level by PIV is higher than that by CFD generally, the main turbulence features are nearly the same. Furthermore, the blade orientation effect and other associated unsteady phenomena are also examined, in order to enhance the understanding on impeller-diffuser interaction in a radial diffuser pump.     

壁面在展向作周期运动的槽道湍流的大涡模拟

分别采用3种亚格子模式:传统的Smagorinsky模式、动力Smagorinsky模式和Cui(2004)基于Kolmogorov方程所提出的新模式,对壁面在展向作周期运动的槽道湍流进行了大涡模拟,以考察这3种模式对平均运动为三维、非定常的湍流流动的模拟能力.通过对湍流基本统计量的分析,发现动力模式和新模式都可以较好地预测这种三维非定常的湍流流动;对相位平均的湍流统计量,动力模式的结果略优于新模式;传统的Smagorinsky模式对这种流动的预测结果是最差的.     

应用分离涡模型计算斜圆柱孔气膜冷却

分离涡模型(DES,Detached eddy simulation)兼有雷诺时均湍流模型计算量较小和大涡模拟计算精度高的优点。本文利用DES对平板单斜圆柱孔的气膜冷却进行了数值模拟,并将结果同RANS湍流模型的计算结果以及实验数据相比较,表明DES能有效弥补RANS湍流模型在计算三维不均匀非定常湍流场的不足,更接近实际地反映了气膜冷却中的流动和换热的本质规律。