高超声速激波湍流边界层干扰直接数值模拟研究

高超声速激波与湍流边界层干扰会导致飞行器表面出现局部热流峰值,严重影响飞行器气动性能和飞行安全. 针对高马赫数激波干扰问题,以往数值研究多采用雷诺平均方法,而在直接数值模拟方面的相关工作较为少见. 开展高超声速激波与湍流边界层干扰的直接数值模拟研究,有助于进一步提升对其复杂流动机理认识和理解,同时也将为现有湍流模型和亚格子应力模型的改进提供理论依据. 采用直接数值模拟方法对来流马赫数6.0,34°压缩拐角内激波与湍流边界层的干扰问题进行了研究. 基于雷诺应力各向异性张量,分析了高超声速湍流边界层在压缩拐角内的演化特性. 通过对湍动能输运方程的逐项分析,系统地研究了可压缩效应对湍动能及其输运的影响机制. 采用动态模态分解方法,探讨了干扰流场的非定常运动历程. 研究结果表明,随着湍流边界层往下游发展,近壁湍流的雷诺应力状态由两组元轴对称状态逐渐演化为两组元状态,外层区域则由轴对称膨胀趋近于各向同性. 干扰流场内存在强内在压缩性效应(声效应),其对湍动能输运的影响主要体现在压力--膨胀项,而对膨胀--耗散项影响较小. 高超声速下压缩拐角内的非定常运动仍存在以分离泡膨胀/收缩为特征的低频振荡特性,其物理机制与分离泡剪切层密切相关.      

高超声速激波湍流边界层干扰直接数值模拟研究

高超声速激波与湍流边界层干扰会导致飞行器表面出现局部热流峰值,严重影响飞行器气动性能和飞行安全.针对高马赫数激波干扰问题,以往数值研究多采用雷诺平均方法,而在直接数值模拟方面的相关工作较为少见.开展高超声速激波与湍流边界层干扰的直接数值模拟研究,有助于进一步提升对其复杂流动机理认识和理解,同时也将为现有湍流模型和亚格子应力模型的改进提供理论依据.采用直接数值模拟方法对来流马赫数6.0,34?压缩拐角内激波与湍流边界层的干扰问题进行了研究.基于雷诺应力各向异性张量,分析了高超声速湍流边界层在压缩拐角内的演化特性.通过对湍动能输运方程的逐项分析,系统地研究了可压缩效应对湍动能及其输运的影响机制.采用动态模态分解方法,探讨了干扰流场的非定常运动历程.研究结果表明,随着湍流边界层往下游发展,近壁湍流的雷诺应力状态由两组元轴对称状态逐渐演化为两组元状态,外层区域则由轴对称膨胀趋近于各向同性.干扰流场内存在强内在压缩性效应(声效应),其对湍动能输运的影响主要体现在压力-膨胀项,而对膨胀-耗散项影响较小.高超声速下压缩拐角内的非定常运动仍存在以分离泡膨胀/收缩为特征的低频振荡特性,其物理机制与分离泡剪切层密切相关.     

并列三柱体风致干扰DDES模拟分析

采用基于SSTk-ω湍流模型的延迟分离涡模拟方法DDES(Delayed detached eddy simulation),对三并列线性布置柱体的风荷载和流场特性进行了数值模拟.为验证方法的有效性,对单个三维方柱绕流进行了计算分析,比较计算值与试验结果的风荷载特征,两者符合较好.然后采用DDES研究了三并列方柱对称布置时,间距比对中间受扰建筑风致干扰的问题.数值结果表明,间距比较小时对受扰建筑风荷载分布有较大影响.总的来看,当间距比约为3.0时,受扰建筑的影响达到最大.同时,数值分析了三方柱周围瞬态流场特征,探讨了临近建筑之间的干扰机理.结果 表明,SST-DDES方法能够较好地预测绕方柱流动分离问题,可为群体高层建筑气动相互干扰研究提供参考.     

下击暴流作用下群体高层建筑风压干扰特性的数值分析

针对布局变化对高层建筑下击暴流场气动干扰特性研究的缺乏,采用RNG k-ε湍流模型封闭求解N-S方程的RANS方法模拟下击暴流场,在与已有单体高层建筑下击暴流场试验数据进行对比验证的基础上,模拟研究两幢高层建筑布局的下击暴流风效应场。考虑建筑布局距风暴中心距离变化以及布局纵横向间距变化时,施扰建筑对受扰建筑立面风压的影响,获取立面风压干扰规律与特性。结果表明:受扰建筑迎风面干扰因子随风暴中心距离的增大而增大,突出表现在较大距离区间(r2D),且干扰作用随布局纵横向间距的变化存在较明显差异;当距风暴中心较近时(r2D),布局纵横方向间距变化对受扰建筑干扰因子的影响较小;在径向最大风速位置处(r/D=1.0),临近施扰建筑的受扰建筑侧风面风压相对单体情况有所减小,其干扰因子随布局纵向间距的增大而减小,随横向间距的增大而增大。     

脊状表面减阻机理研究

针对脊状表面流场的特点,通过实验测量和数值模拟的方法对脊状表面微观流场进行了深入研究,获得了脊状表面湍流边界层的时均速度分布曲线、湍流度分布曲线和微观流场结构.为了得到脊状结构对壁面物性的影响,对脊状表面进行了疏水性测试,获得了液滴在脊状表面上的表观接触角,并通过水洞试验验证了脊状表面的减阻效果.研究表明,与光滑表面相比,脊状表面微观流场结构中存在"二次涡",近壁区的黏性底层厚度比平板的要厚得多,湍流度显著降低,且脊状表面表现出明显的疏水性.由此提出了基于壁面隔离效应、增大湍流阻尼效应和改变壁面物性效应的减阻机理.     

高层建筑标准模型风效应的阻塞效应试验研究

对一组阻塞比为1.7%～10.9%的高层建筑标准模型在均匀低湍流强度风场和均匀高湍流强度风场中进行了多点同步测压试验,研究了不同风场下阻塞效应对高层建筑风致响应和等效静力风荷载的影响。结果表明:当阻塞比不大于4.5%时,可以忽略高层建筑风致响应和等效静力风荷载的阻塞效应;当阻塞比大于4.5%时,在均匀低湍流强度风场中,高层建筑的风致响应和等效静力风荷载均随阻塞比的增大而增大;在均匀高湍流强度风场中,湍流强度的增大不同程度地削弱了阻塞效应的影响。基于试验结果,提出了高层建筑顶部平均位移、位移均方根和加速度均方根的阻塞效应修正公式,可应用于修正等效静力风荷载,修正效果较好。     

??????????潨????绷?????????????

采用数值模拟方法,结合风洞模型试验对群体低层四坡屋面房屋周围的 风场及表面风压进行了计算和分析.数值模拟基于Reynolds时均方程和可实现 $k$-$\varepsilon$湍流模型,采用了具有良好拓扑性能的非结构四面体网 格,运用流体软件FLUENT实现流场的求解．在单体计算结果和风洞试验结果有较好吻合的前 提下,获得了有相邻建筑干扰的情况下,低层四坡屋面房屋的表面风压的变化规律,结论可 直接供工程设计参考．     

湍流横掠固定颗粒的全尺度直接数值模拟

气固两相流模拟中,当固相尺度接近或大于Kolmogorov尺度时,普通的点源模型将不再适用,固体相的体积效应和表面效应将对流体相产生显著的影响。通过采用直接数值模拟方法,结合内嵌边界方法对湍流中不同湍流强度流体横掠大于Kolmogorov尺度的固相颗粒进行了全尺度模拟,讨论分析了在两种湍流度下方形颗粒对湍流的调制影响以及颗粒的受力情况。     

两个相邻建筑物周围风环境的数值模拟

高层建筑物及其群体周围所引起的风环境问题正日益受到人们的重视。本文用数值模拟方法对两幢串列布置高层建筑周围的流场进行了模拟，讨论了建筑物高度、间距比等因素对风环境特性的影响。为了与两个相邻建筑物的流场进行比较，本文同时计算了单个建筑物周围流动的分布情况。数值模拟结果对揭示高层建筑群区域内再生风载荷的形成机理提供了一定的分析方法，为城市规划和高层建筑住宅小区设计提供了一定的理论参考和依据。     

低密度流体界面不稳定性大涡模拟

采用拉伸涡亚格子尺度应力模型对湍流输运中的亚格子作用项进行模式化处理,发展了适用于可压多介质黏性流动和湍流的大涡模拟方法和代码MVFT（multi-viscous flow and turbulence）。利用MVFT代码对低密度流体界面不稳定性及其诱发的湍流混合问题进行了数值模拟。详细分析了扰动界面的发展,流场中冲击波的传播、相互作用、湍流混合区边界的演化规律,以及流场瞬时密度和湍动能的分布和发展。数值模拟获得的界面演化图像和流场中波系结构与实验结果吻合较好。三维和二维模拟结果的比较显示,两者得到的扰动界面位置、波系及湍流混合区边界基本一致,只是后期的界面构型有所不同,这也正说明湍流具有强三维效应。     

An experimental study on wind loads acting on a high-rise building model induced by microburst-like winds

In the present study, an experimental investigation is conducted to quantify the characteristics of the microburst-induced wind loads (i.e., both static and dynamic wind loads) acting on a high-rise building model, compared to those with the test model placed in conventional atmospheric boundary layer (ABL) winds. The experimental study is performed by using an impinging-jet-based microburst simulator available at Iowa State University. In additional to conducting flow field measurements to quantify the flow characteristics of the microburst-like wind, both mean and dynamic wind loads acting on the test model induced by the microburst-like wind are assessed in detail based on the quantitative measurements of the surface pressure distributions around the test model and the resultant aerodynamic forces. It is found that the microburst-induced wind loads acting on high-rise buildings would be significantly different from their counterparts in conventional ABL winds. Both the static and dynamic wind loads acting on the high-rise building model were found to change significantly depending on the radial locations and the orientation angles of the test model in respect to the oncoming microburst-like wind. The dynamic wind loads acting on the test model were found to be mainly influenced by the periodical shedding of the primary vortices and the high turbulence levels in the microburst-like wind. The findings derived from the present study are believed to be useful to gain further insight into the underlying physics of the flow–structure interactions of high-rise buildings in violent microburst winds for a better understanding of the damage potential of microburst winds to high-rise buildings.     

Turbulence Intensity and Length Scale Effects on Premixed Turbulent Flame Propagation

The effects of varying turbulence intensity and turbulence length scale on premixed turbulent flame propagation are investigated using Direct Numerical Simulation (DNS). The DNS dataset contains the results of a set of turbulent flame simulations based on separate and systematic changes in either turbulence intensity or turbulence integral length scale while keeping all other parameters constant. All flames considered are in the thin reaction zones regime. Several aspects of flame behaviour are analysed and compared, either by varying the turbulence intensity at constant integral length scale, or by varying the integral length scale at constant turbulence intensity. The turbulent flame speed is found to increase with increasing turbulence intensity and also with increasing integral length scale. Changes in the turbulent flame speed are generally accounted for by changes in the flame surface area, but some deviation is observed at high values of turbulence intensity. The probability density functions (pdfs) of tangential strain rate and mean flame curvature are found to broaden with increasing turbulence intensity and also with decreasing integral length scale. The response of the correlation between tangential strain rate and mean flame curvature is also investigated. The statistics of displacement speed and its components are analysed, and the findings indicate that changes in response to decreasing integral length scale are broadly similar to those observed for increasing turbulence intensity, although there are some interesting differences. These findings serve to improve current understanding of the role of turbulence length scales in flame propagation.

     

Y型柱体绕流特性及风荷载

通过风洞实验研究了二种典型尺寸的二维Y型柱体在不同风向角下的绕流特性及风荷载。实验结果表明:当风向顺Y型柱体某一肢时主要是来流脉动引起柱体较弱横向振动;当风向顺Y型柱体两肢分角线时则由涡脱落而引起柱体强烈的横向振动。来流湍流度的增加使脉动升力和阻力都大幅度增加,然而却使涡脱落引起的振动相对减弱。     

Inlet grid-generated turbulence for large-eddy simulations

A new technique of generating turbulence in large-eddy simulations (LES) has been investigated and results compared with previous studies for validation. The proposed gridInlet technique uses a grid pattern on the inlet boundary patch to produce grid-generated turbulence as used in wind tunnel experiments. This allows the turbulence integral length scale to be controlled by changing the grid size, while the turbulence intensity is controlled by changing the inlet distance. The objective of this paper is to investigate domain and mesh requirements to implement the gridInlet technique. This technique is most suited to studies on the influence of high-intensity isotropic turbulence on objects, particularly if comparisons are to be made to experimental data obtained with grid-generated turbulence.     

Behavior of homogeneous turbulence in channels with variable cross-sectional area

On the basis of the equation describing the behavior of the spectral tensor of the energy of the pulsation velocity for one-dimensional flow in a channel with variable cross-sectional area, we obtained a system of equations for the meansquare components of the pulsation velocity vector and integral scales of the turbulence length in various directions, enabling us to use these parameters in the initial section of the channel for determining their behavior along the channel. We made use of some ideas of A. N. Kolmogorov and J. Rotta concerning the possibility of describing viscous and nonlinear terms in the equations for the components of the tensor of Reynolds stresses in terms of the energy of pulsation motion and the integral scale of the turbulence length. The resulting system, in the special cases of very low intensity of turbulence, leads to the results of the linear theory; for constant cross-sectional area (or for a very high intensity of turbulence, when the damping affects the turbulence much more strongly than does the deformation effect) it describes the known empirical laws of the degeneration of turbulence beyond grids. We made a comparison with the data available in the literature on the behavior of the characteristics of turbulence in channels with variable cross-sectional area.     

公寓铝幕墙风荷载体型系数精细化分析

本文以北京冬奥运村公寓建设项目为工程背景,采用计算流体力学,基于大涡模拟对铝幕墙空腔装饰构件的风荷载体型系数进行研究．研究发现,在不同工况下,装饰构件法向风荷载体型系数整体大于切向风荷载体型系数,拐角区域的装饰构件会比平直区域的装饰构件受到更大的风压作用;建筑的风荷载体型系数的变化幅度会随着所在高度的增加而愈发明显．研究结论对该工程的幕墙设计及其连接件设计提供了精细化的荷载数据．     

大气边界层流场的模拟与大气边界层风洞

动态风荷载测试数据的可靠性依赖于大气边界层流场的正确模拟。本文根据作者在ＦＤ－４大气边界层风洞中,进行大气边界层流场调试时,就如何准确模拟自然风的四个参数方面所积累经验的总结。     

横摆角下汽车尾部湍流特征量的PIV试验分析

汽车尾部湍流场是汽车压差阻力的主要来源,在HD-2汽车模型风洞中,首先使用测力天平和测压系统,对横摆角工况下汽车模型的气动六分力和纵对称截面48个测点的表面压力进行了测量,然后利用PIV测量技术对模型在横摆角分别为0°、15°的尾部湍流场进行了测量,获得该模型尾流场的速度场、涡量场和雷诺应力流场信息,通过计算得出尾流场区域空间相关系数和湍流积分尺度。结果表明:在横摆角工况下,汽车模型尾部涡流的结构呈现向上发展的趋势;尾流场拖拽涡的范围和强度的增大导致了模型气动力出现较大的增加;湍流积分尺度的变化表明,尾部涡流区的分离噪声与涡流分离位置有关,在汽车尾部造型设计中,要尽量推迟尾部涡流的分离。     

Evaluation of bubble-induced turbulence using direct numerical simulation

The presented research evaluates the interaction between a single bubble and homogeneous turbulent flow using direct numerical simulation (DNS) approach. The homogeneous single-phase turbulence is numerically generated by passing a uniform flow through grid planes. The turbulence decay rate is compared with experiment-based correlation. The single phase turbulence is then used as an inflow boundary condition for a set of single bubble studies. By estimating the turbulent field around the fully resolved bubble, the effects of bubble deformability, turbulent intensity and relative velocity on the bubble-induced turbulence are investigated. The existence of bubble creates new vortices in the wake region and the enhancement of turbulence is observed in the region behind the bubble. The results show that the magnitude of the turbulence enhancement would increase as the bubble encounters larger liquid turbulent intensity or higher relative velocity. Set of bubble Weber numbers from 0.34 to 3.39 are used to investigate the effect of bubble deformability. The more deformable bubble is the higher the increase in the magnitude of the turbulence enhancement behind the bubble. This research provides systematic insight on the bubble-induced turbulence (BIT) mechanism and is important for multiphase computational fluid dynamics (M-CFD) closure model development.     

Modeling and Simulation of Effects of Turbulence on Vaporization, Mixing and Combustion of Liquid-Fuel Sprays

The objective of this work is twofold. Firstly, the effects of turbulence intensity variations on the turbulent droplet dispersion, vaporization and mixing for non-reacting sprays (with and without swirl) are pointed out. Secondly, the effects of the coupling of the turbulence modulation with external parameters, such as swirl intensity, on turbulent spray combustion are analyzed in configurations of engineering importance. This is achieved by using advanced models for turbulence, evaporation and turbulence modulation implemented into FASTEST-LAG3D-codes: (1) To highlight the influence of turbulence modulation on some spray properties, a thermodynamically consistent modulation model has been considered besides the standard assumption and the well known Crowe's model. For turbulent droplet dispersion, we rely on the Markov-sequence formulation. (2) In order to characterize phase transition processes ongoing on droplets surfaces, a non-equilibrium evaporation model shows better agreement with experiments in comparison with the quasi-equilibrium-based evaporation models often used. (3) The results of turbulence intensity variations reveal the existence of a limited range out of which the increase or decrease of the turbulence intensity affects no more the efficiency of the heat and mass transfer. A derived characteristic number, a vaporization Damkhöler number, possesses a critical value which separates two different behavior regimes with respect to the turbulence/droplet vaporization interactions. (4) Under reacting conditions, it is shown how the evaporation characteristics, mixing rate and combustion process are strongly influenced by swirl intensity and turbulence modulation. In particular, the turbulence modulation modifies the evaporation rate, which in turn influences the mixing and the species concentration distribution. In the case under investigation, it is demonstrated that this effect cannot be neglected for low swirl intensities (Sw.Nu. ≤ 1) in the region far from the nozzle, and close to the nozzle for high swirl number intensities. In providing these particular characteristics, a reliable control of the mixing of gaseous fuel and air in evaporating and reacting sprays, and a possible optimization of the mixing process can tentatively be achieved.