一种改进的分层剪切湍流亚格子尺度模型及其应用

提出了一种新型的分层剪切湍流应力和湍流热通量的动力学亚格子尺度(SGS)模型,并就热分层槽道湍流问题,计算验证了该模型的正确性.采用大涡模拟方法和该动力学亚格子尺度模型,数值研究了在浮力和剪切双重作用下的稳定分层和不稳定分层槽道热剪切湍流,探讨了某些湍流特性随Richardson数的变化规律,并分析了相关的流动物理机制.     

分层剪切湍流大涡模拟的一种动力学亚格子尺度模型

基于Yoshizawa涡粘性模型提出了一种适用于热分层剪切湍流大涡模拟的动力学亚格子尺度 (SGS)模型 ,包括SGS湍流应力和湍流热通量模型 ,计算验证了该模型的正确性 .进一步采用大涡模拟方法和该动力学SGS模型研究了稳定分层和不稳定分层槽道湍流的湍流特性 ,以及湍流主要变量关联量的变化规律 .同时 ,根据计算结果 ,预测了热分层槽道湍流的临界Richardson数 ,与理论分析结果基本相符 .     

湍流不同尺度间的近程和共振作用规律及应用

从Navier-Stokes方程出发,研究了湍流不同尺度间的相互作用规律,给出相近尺度间近程粘性应力的积分和微分表达式．引入极相近尺度之间共振相互作用的概念,得到共振粘性应力的微分表达式．利用共振粘性应力张量获得不含经验关系和常数、近似封闭的大涡模拟（LES）方程组．利用近程和共振粘性应力张量获得不含经验关系和常数、近似封闭的湍流多尺度方程组．讨论了湍流多尺度方程的性质及用于湍流计算的优点,尺度间相互作用的近程特性说明:多尺度模拟是湍流计算很有价值的方法,并列举了算例．     

双尺度二阶矩颗粒相湍流模型经验系数的影响

基于将颗粒脉动分成湍流引起的大尺度脉动和颗粒间碰撞产生的小尺度脉动的概念,建立了双尺度二阶矩两相湍流模型．用该模型对下行床内两相流动进行了数值模拟,颗粒体积浓度、平均速度的计算结果和实验数据吻合较好．分析了双尺度二阶矩两相湍流模型经验系数变化对预报结果的影响:在经验系数的一定变化范围内,预报结果并无明显的影响,但是变化范围增大,预报结果会产生较大变化．     

基于时空最优低维动力系统的多尺度可压缩湍流数值模拟方法（英文）

按照周培源教授关于研究湍流数值模拟建模时必须分析和求解脉动速度场的思想,该研究基于第一性原理,系统地建立了基于时空低维最优动力系统的多尺度可压缩湍流数值模拟方法(LMS方法),并将其应用于多次冲击Richtmyer-Meshkov问题的数值模拟中,首次得到了可压缩湍流的中尺度流场和不同于DNS近似解的湍流近似解.数值结果表明,LMS方法可以用较少的网格获得更精确的湍流近似解.首先解决了研究中遇到的几个问题,为LMS方法的构建铺平了道路.这些问题是:基于湍流的物理特性,提出了湍流大、中、小尺度分解的新概念;找到了box滤波空间相关性的计算方法;指出了湍流建模理论中长期存在的逻辑错误,提出了多尺度湍流模型的概念;讨论了湍流封闭问题的本质和关键,给出了克服湍流封闭问题的数值方法.采用box滤波方法/空间网格平均方法且在大尺度网格的意义下,LMS方法的本质是一种将RANS、LES、DES和DNS等湍流数值模拟方法统一的全新湍流数值模拟方法.需要指出的是,LMS方法也可以作为湍流模型研究的辅助工具,以检验SGS尺度方程/脉动方程中各项所对应的湍流模型是否正确.     

宾汉流体稠密两相湍流流动中的二阶矩模型

建立的Bingham流体稠密两相流动的二阶矩-颗粒动力论湍流模型（USM-theta模型）既体现了两相的作用,又体现了屈服应力所引起的附加项,并提出了USM-theta模型下考虑浓度修正值影响的两相湍流流动的算法．利用该模型对圆管内Bingham流体的单相湍流流动、稠密液固两相的湍流流动进行了计算,并和五方程湍流模型进行了比较,结果表明该模型的预测效果更好．利用USM-theta模型对含颗粒的Bingham流体的两相湍流流动进行了模拟,随着屈服应力的增加,Bingham流体相与颗粒相在管道中心附近的主流速度减小．液固两相湍流和Bingham流体两相湍流的计算结果表明屈服应力引起的附加项对流动有很重要的影响．     

突然扩张方管中三维湍流流动的数值模拟

本文运用SIMPLEC算法计算了突然扩张方管中的三维湍流流动,湍流模型采用k-ε模型。计算结果详细反映了突然扩张方管中三维湍流流场。从本文结果可以看出,由于突然扩张方管几何形状非轴对称,且尺寸有限,边壁对流场的作用是不可忽略的。以往文献中常见的二维突然扩张湍流的数值模拟结果与三维情况有较大差别,在靠近边壁的区域差别很大,因此对于突然扩张方管中湍流流动的数值模拟应用三维模拟。本文计算所得突然扩张截面后主回流区长度与实验结果接近。本文方法可为数值模拟突然扩张方管中湍流流场及各物理参数的分布提供有效工具。     

二维剪切流的粘性-无粘湍流干扰理论

对二维不可压缩近壁剪切湍流,本文提出一个粘性-无粘湍流干扰理论.主要内容有:从分子粘性考虑出发确定干扰湍流的流动结构及其物理尺度,导出空间为小尺度的局部流动结构随顺流距离的演变规律,导出支配干扰湍流流动的简化Reyno-lds(SR)方程和扩散抛物化K-ε方程.该SR方程是作者早先提出的简化Navier-Stokes(SNS)方程的湍流形式,它的重要性质是“简化运算”和时间Reynolds平均运算的顺序可以交换.关于最大湍流剪应力、本理论计算值与实验测量值很好相符.经典湍流边界层理论、Clauser平衡湍流边界层以及湍流分离Triple-deck理论均是本文理论的特例.证实了顺流方向长度尺度随干扰增强而显著减小的实验结论.     

近壁区Reynolds应力及湍能耗散率输运特性的理论模型 *

提出在湍流边界层的近壁区采用三维不稳定波来描述湍流相干结构 ,然后根据理论模型对Reynolds应力及湍能耗散率的输运特性进行系统的计算和定量的分析 ,并计算了平均速度分布 .所得理论结果与直接数值模拟结果符合较好 ,表明本文方法正确地体现了湍流边界层近壁区的物理实质 .这不仅有益于对湍流机制的了解 ,而且可能为湍流的近壁模型化开辟一条新的途径 .     

湍流边界层近壁区相干结构起因的研究

利用流动稳定性理论中的一般共振三波的概念,提出了一种湍流边界层近壁区相干结构产生机理的理论模型。由此所得的相干结构的空间形态,展向尺度,传播速度等都与数值模拟所得结果相近。特别是计算了多数是不对称的各种流向涡对环量差的概率密度分布曲线,得到了与数值模拟结果相比基本满意的结果。     

Assessment of an improved turbulence model in simulating the unsteady flows around a D-shaped cylinder and an open cavity

Most engineering flows are still predicted by the conventional Reynolds-averaged Navier-Stokes method because of the low requirements of the computational quantities. However, the resolution capability of Reynolds-averaged Navier-Stokes models is still open to deliberation, especially in the recirculation and wake regions, where the vortical flows dominate. In the present work, an improved turbulence model derived from the original shear stress transport k-ω model is proposed and its superiority is assessed by our modeling the unsteady flows around a D-shaped cylinder and an open cavity, corresponding to two different Reynolds numbers. The results are compared with results from experiments and other turbulence models in terms of the flow morphology and mean velocity profiles. This shows that the predictive accuracy of the modified turbulence model is increased significantly in the bluff body wake flows and in the shear layer and separation flows of the cavity. Some special vortex structures can be captured in the open cavity, in which the secondary vortex emerging from the shear layer and the separation vortex near the trailing edge can induce large flow instability, and this phenomenon should be eliminated in engineering applications. It is believed that this improved turbulence model can be used for the more complex turbomachinery flows with better prediction of the hydrodynamic/aerodynamic performance and the unsteady vortical flows, which can provide some guidelines to design or optimize rotating machines.     

Smoke flow phenomena and turbulence characteristics of tunnel fires

Gravity currents are similar in behavior with smoke flows. This work aims to provide evidence justifying the use of gravity current approach to model smoke flows downstream of the fire source. The turbulence solver available in almost all commercial CFD codes solves RANS for the flow field. To find out how well the nature of smoke flow be accurately modeled using RANS that is widely used for incompressible flows. The feasibility of using both Reynolds- and Favre-averaging schemes was numerically compared and examined in this paper. In this work, numerical simulations of a fire occurred in a 400-m longitudinally ventilated tunnel have been successfully performed using FDS version 4. Large eddy simulation is employed in this study. Although the ranges of fire size and ventilation velocity vary respectively from 0 MW to 100 MW and 0 m/s to 10 m/s, this paper focuses on the general flow and temperature fields and the turbulence characteristics. Furthermore, the turbulence kinetic energy levels of the flow in the tunnel at several locations were investigated. Since the flow field is generally induced by mechanical ventilation and combustion, the main contribution to the turbulence kinetic energy comes from its longitudinal, vertical, or their combination.     

反映强流动曲率效应的非线性湍流模型

首先定性地分析了流线曲率效应对流场湍流结构的影响,然后以U型槽道流为典型算例,对多种湍流模型进行了评估．评估的模型包括:线性涡粘性模型,二阶和三阶非线性涡粘性模型,二阶显式代数应力模型和Reynolds应力模型．评估结果表明,性能良好的三阶非线性涡粘性模型,如黄于宁等人发展的HM模型以及CLS模型,可以较好地描述流线的曲率效应对湍流结构的影响,如凸曲率作用下内壁附近湍流强度的衰减和凹曲率作用下外壁附近湍流的增强,并且较好地确定了管道下游的分离点位置和分离泡长度,其预测的结果和实验符合较好,与Reynolds力模型的结果十分接近,因此可以较好地应用于具有曲率效应的工程湍流的计算．     

Comparison of turbulence models in simulating swirling pipe flows

Swirling flow is a common phenomenon in engineering applications. A numerical study of the swirling flow inside a straight pipe was carried out in the present work with the aid of the commercial CFD code fluent. Two-dimensional simulations were performed, and two turbulence models were used, namely, the RNG k–ε model and the Reynolds stress model. Results at various swirl numbers were obtained and compared with available experimental data to determine if the numerical method is valid when modeling swirling flows. It has been shown that the RNG k–ε model is in better agreement with experimental velocity profiles for low swirl, while the Reynolds stress model becomes more appropriate as the swirl is increased. However, both turbulence models predict an unrealistic decay of the turbulence quantities for the flows considered here, indicating the inadequacy of such models in simulating developing pipe flows with swirl.     

Non-linear k– turbulence model results for flow over a building at full-scale

Turbulence modelling is a crucial question in the application of CFD to flows over buildings. The impinging flow and anisotropic nature of the turbulence present severe challenges. This paper presents a comparison of CFD against full-scale results. It differs from previous work which has concentrated on the wind-tunnel scale. In order to better account for the production of turbulent kinetic energy and the anisotropic nature of the turbulence a non-linear k– model is implemented. The results are discussed for different turbulence models and for the comparison of computed results with the measurements from full-scale.     

固液两相流中的K-ε双方程湍流模式*

建立了固液两相流中更一般的K-ε双方程湍流模式。模化了固相和液相的连续方程、动量方程及K方程和ε方程。该湍流模型考虑了固液两相间速度的滑移,颗粒间的作用及相间作用。使用本文所建立的湍流模型,数值预测了一管湍流中的沙水混合流动,其预测结果与实验结果比较一致。     

Modified k – ω Model for Simulation of Cavitating Flows

Although cavitating flows are generally turbulent, the interaction of turbulence and cavitation dynamics is still not well understood. In general, two‐equation models are employed, which were originally developed for single‐phase flows. Therefore they fail by handling cavitation based flow phenomena with very high density variations (dependent on operating condition up to 40000:1). This sudden change of the density causes strong pressure gradients, secondary flows and local compressibility. The aim of this study is to enhance the Wilcox's k‐? model with empirical correlations in order to simulate turbulent cavitating flows more precisely and effciently.     

Comparison of Inviscid and Viscous Two‐Phase Calculations of Cavitating Pump Flow

In cavitating flows there is a strong interaction between the fine dispersed vapor bubbles and turbulence. Therefore in two‐phase flow calculations the prediction of turbulence is a matter of great difficulties and uncertainties. To get an idea of the influence of turbulence modelling on the calculated result, viscous and inviscid two‐phase calculations of cavitating pump flow were performed. In the inviscid calculations the influence of cavitation is isolated from turbulence effects. In the viscous calculations the effect of turbulence is modelled with a k – ϵ turbulence model. The results show that the influence of viscous effects on the flow field is weak in comparison to cavitation. However, in contrast to the steady cavitation behaviour predicted by the viscous calculations, the inviscid calculations show unsteady behaviour of the cavitation (as can be seen in the experiment).     

Numerical modelling of subcritical open channel flow using the K-ε turbulence model and the penalty function finite element technique

A numerical model has been developed that employs the penalty function finite element technique to solve the vertically averaged hydrodynamic and turbulence model equations for a water body using isoparametric elements. The full elliptic forms of the equations are solved, thereby allowing recirculating flows to be calculated. Alternative momentum dispersion and turbulence closure models are proposed and evaluated by comparing model predictions with experimental data for strongly curved subcritical open channel flow. The results of these simulations indicate that the depth-averaged two-equation k-ε turbulence model yields excellent agreement with experimental observations. In addition, it appears that neither the streamline curvature modification of the depth-averaged k-ε model, nor the momentum dispersion models based on the assumption of helicoidal flow in a curved channel, yield significant improvement in the present model predictions. Overall model predictions are found to be as good as those of a more complex and restricted three-dimensional model.     

Influence of turbulence model constants and boundary conditions on the computed hydrodynamic variables in gas-stirred reactors

A finite difference procedure has been employed to obtain numerical predictions of recirculating flows in gas-stirred reactors. Sensitivity of alternative boundary conditions at gas-liquid interface as well as empirical constants of the turbulence model on the predicted flow field have been investigated.