Pr数对湍流被动标量输运的影响

采用直接数值模拟的方法,研究分子Pг数对湍流被动标量输运的影响,并提供充分的证据证明,湍流Pг数明显依赖于分子Pг数.在算例中,湍流雷诺平均PгT数与分子Pг数的倒数呈线性关系;湍流亚格子Pгt数与分子Pг数的关系较为复杂,在分子Pг数为1附近时,湍流亚格子Pгt数出现极小值.     

一种用于分析封闭腔内湍流自然对流换热的k-ε新模型

由于目前用于求解湍流自然对流流动与传热的k-ε模型在应用过程中存在不足,结合高雷诺数k-ε模型需要借助壁面函数法来确定壁面上相关参数值和低雷诺数k-ε模型在近壁区布置更多节点以便获得粘性底层详细信息的特点,重新定义了湍流普朗特数σt的计算式,提出了一种修正的k-ε新模型;利用该模型对封闭方腔内的湍流自然对流流动与传热进行了数值分析。结果表明:与文献中数值模拟结果相比,当108≤Ra≤1014时本文模型所得壁面平均努塞尔特数更接近文献中的实验值,与实验值之间的相对误差在8%以内;壁面的局部努塞尔特数与文献中的实验值吻合得较好。这说明本文模型用于求解封闭腔内湍流流动与传热问题是合适的,比其它湍流模型更能准确地描述封闭腔内湍流自然对流换热中边界层发展与壁面传热特性之间的内在联系。     

矩形液池热毛细对流转捩途径研究

主要研究矩形液池热毛细对流的分岔转捩. 通过测量流体内部温度振荡情况, 详细研究了热毛细对流的转捩过程和转捩途径. 实验发现, 矩形液池热毛细对流的转捩过程依次经历了定常、规则振荡、不规则振荡的阶段. 对于不同普朗特数的硅油在不同长高比情况下, 通向混沌的途径不同. 在转捩过程中, 随着温差的增加, 普朗特数在16 (1cSt) 以下和普朗特数为25 (1.5cSt)、长高比为26 的硅油热毛细对流主要以准周期分岔的转捩方式为主;而普朗特数为25 以上的则以倍周期分岔的转捩方式为主;两种分岔有时还会伴随有切分岔形式的出现.实验中还观察到了表面波动和对流涡胞振荡等现象.      

可压缩自由剪切流混合转捩大涡模拟

针对湍流气动光学效应与冲压发动机气体混合机理问题,开展了可压缩混合层流动空间模式大涡模拟和时间模式直接数值模拟研究.通过对流场(包含亚/亚混合、超/亚混合两种情况)失稳、转捩直至完全湍流的空间发展过程的研究表明,对流Mach数0.4状态下流场失稳以二维最不稳定扰动为主;非线性发展中,基频涡对并及展向涡撕裂主控流动转捩,流场发生混合转捩;转捩后脉动流场基本达到局部各向同性,此时,湍流Mach数低于0.3,流动压缩性可近似忽略.     

侧加热腔内的自然对流

开展侧加热腔内自然对流的研究具有重大的环境及工业应用背景. 总结侧加热腔内水平温差驱动的自然对流的最新研究进展, 并概述相应的流动性质、动力机制和传热特性以及对不同无量纲控制参数的依赖也有重要的科学价值. 已取得的研究结果显示突然侧加热的腔内自然对流的发展可包括初始阶段、过渡阶段和定常或准定常阶段. 不同发展阶段的流动依赖于瑞利数、普朗特数及腔体的高宽比, 且定常或准定常阶段的流态可以是定常层流流动、非定常周期性流动或者湍流流动. 此外, 回顾了对流流动失稳机制的研究成果以及湍流自然对流方面的新进展. 最后, 展望了侧加热腔内的自然对流研究的前景.      

高径比对GaAs熔体液桥热毛细对流失稳的影响

采用基于谱元法线性稳定性分析方法,研究了高径比对GaAs熔体(Pr=0.068)液桥热毛细对流失稳的影响,同时结合能量分析揭示了热毛细对流的失稳机制.研究结果表明:与典型低普朗特数(例如Pr=0.011)熔体静态失稳模式和典型高普朗特数(例如Pr>1)熔体振荡失稳模式不同,GaAs熔体热毛细对流失稳模式依赖于液桥高径比...     

亚网格尺度稳定化有限元求解不可压黏性流动

从亚网格尺度稳定化方法的基本原理出发, 提出了适合时间推进求解非定常Navier-Stokes方程获得定常解的SGS稳定化方法. 基于一定程度的近似和简化, 获得了与时间步长相关的稳定化参数, 从而排除了传统SGS稳定化方法在求解高Re数、小时间步长问题时所引发的数值不稳定性. 把SGS稳定化方法应用于求解不可压湍流, 结合标准k-\varepsilon湍流模型和壁面函数法估计湍流黏性系数, 详细讨论了壁面函数法的实施、湍流输运方程的求解和保证湍流变量非负性的限制策略, 发展了时间推进求解不可压湍流的分离式算法. 二维外掠后台阶层流和湍流计算结果表明,该方法求解不可压黏性流动是可行的, 并且具有稳定性好、计算精度高的特点.      

超声速平板边界层斜波失稳转捩过程研究

以5阶迎风和6阶对称紧致格式混合差分求解三维可压缩滤波Navier-Stokes方程，对Mach 数为4.5, Reynolds数为10000的空间发展平板边界层湍流进行了大涡模拟. 时间推进采用 紧致存储3阶Runge-Kutta方法，亚格子尺度模型为修正Smagorinsky涡黏性模型. 通过在 入口边界叠加一对线性最不稳定第一模态斜波扰动，数值模拟得到了平板层流边界层失稳转 捩直至湍流的演化过程. 对流场转捩过程中瞬时量及统计平均量的分析表明，数值模拟结果 与理论吻合，得到的Y型剪切层、交替\Lambda涡结构以及转捩后期的发卡涡结构的发展 变化与相关文献结果一致，湍流流谱定性合理.     

介绍《普朗特流体力学基础》

德国力学家普朗特(Ludwig Prandl,1875年2月4日生于弗赖辛,1953年8月15日卒于哥廷根),是近代力学的奠基人之一,他在流体力学方面作出过许多开创性贡献,如边界层理论、风洞实验技术、机翼的举力线和举力面理论、湍流混合长度理论等,都对流体力学的发展产生了巨大的影响.     

????????????????????????????

介绍了求解复杂截面闭口薄壁杆件扭转问题的网络理论解法,以普朗特应力函 数解法为基础,通过与电路问题比拟,借鉴了电路理论中的网络理论解法. 本文方法对解决复 杂截面闭口薄壁杆件扭转的工程问题有一定的应用价值. 在教学上,可以起到拓宽学生思路 的作用.     

Steady turbulent flow and heat transfer downstream of a sudden enlargement in a pipe of circular cross-section

Stream-line and temperature contours, and the corresponding fluxes at the walls, are computed by numerical solution of the elliptic transport equations of vorticity, enthalpy and turbulence energy, together with auxiliary relations comprising a turbulence model similar to that of Prandtl [19]. The length-scale distribution is determined empirically in order to ensure that the recirculation region has the right length, and the maximum of the wall heat flux occurs at the right place, but the other empirical inputs have values which are determined from quite different experiments. — Agreement between predictions and experimental data of Krall and Sparrow [13] is good. In particular, the correct exponent is predicted for the Stanton number ? Reynolds number law. This exponent is uninfluenced by the length-scale distribution. — For practical use, it is argued, the Prandtl turbulence model needs to be replaced by one embodying two differential equations for turbulence quantities.     

Effect of particles on turbulent thermal field of channel flow with different Prandtl numbers

The direct numerical simulation(DNS) of heat transfer in a fully developed non-isothermal particle-laden turbulent channel flow is performed.The focus of this paper is on the modulation of the particles on turbulent thermal statistics in the particle-laden flow with three Prandtl numbers(P r = 0.71,1.5,and 3.0) and a shear Reynolds number(Reτ = 180).Some typical thermal statistics,including normalized mean temperature and their fluctuations,turbulent heat fluxes,Nusselt number and so on,are analyzed.The results show that the particles have less effects on turbulent thermal fields with the increase of Prandtl number.Two reasons can explain this.First,the correlation between fluid thermal field and velocity field decreases as the Prandtl number increases,and the modulation of turbulent velocity field induced by the particles has less influence on the turbulent thermal field.Second,the heat exchange between turbulence and particles decreases for the particle-laden flow with the larger Prandtl number,and the thermal feedback of the particles to turbulence becomes weak.     

Scalar turbulence model investigation with variable turbulent Prandtl number in heated jets and diffusion flames

Traditional turbulence models using constant turbulent Prandtl number fail to predict the experimentally observed anisotropies in the thermal eddy diffusivity and thermal turbulent intensity fields. Accurate predictions depend strongly on the turbulence model employed. Consequently, the objective of this paper is to assess the performance of turbulence model with variable turbulent Prandtl number in predicting of thermal and scalar fields quantities. The model is applied to axisymmetric turbulent round jet with variable density and in turbulent hydrogen diffusion flames using the flamelet concept. The k − ɛ turbulence model is used in conjunction with thermal field; the model involves solving supplemental scalar equations for the temperature variance and its dissipation rate. The model predictions are compared with available experimental data for the purpose of validating model. In reacting cases, velocity and scalar (including temperature and mass fractions) predictions agree relatively well in the near field of the investigated diluted hydrogen flames.     

DNS of turbulent heat transfer in a channel flow with a high spatial resolution

Direct numerical simulations of turbulent heat transfer in a channel flow are performed to investigate the effects of Reynolds and Prandtl numbers on higher-order turbulence statistics such as a turbulent Prandtl number and the budget for the dissipation rate of the temperature variance. The Reynolds numbers based on the friction velocity and the channel half width are 180 and 395, and the molecular Prandtl numbers Pr’s 0.71–10.0. Careful attention is paid to ensure accuracy of the higher-order statistics through the use of a high spatial resolution comparable to Batchelor length scale. The wall-asymptotic value of the turbulent Prandtl number is mostly independent of Reynolds number for the current range of Pr’s. The budget for the dissipation rate of the temperature variance has been computed, and the negligible effect of a Reynolds number on the sum of all source and sink terms in near-wall region in the current computational range is found. This result is quite similar to the one in the budget for the dissipation rate of turbulent energy. In addition, a priori test for existing models is also performed to assess the Pr dependence on the individual terms and their summations in the budget.     

Direct simulation of turbulent flow and heat transfer in a channel. Part I: Smooth walls

Interest in the use of supercomputers for the direct numerical calculation of turbulence prompts the development of efficient numerical techniques so that calculation at higher Reynolds numbers might be made. This paper presents an efficient pseudo-spectral technique, similar to but different from others that have recently appeared, for the calculation of momentum and heat transfer to a constant-property, turbulent fluid in a two-dimensional channel with walls at different, uniform temperature. The code uses no empiricism, although periodic boundary conditions are used for fluctuating quantities in the streamwise and spanwise directions. Calculations were made for a Prandtl number of 0·72 and Reynolds number based on friction velocity and channel half-height of 180 or 2800 based on channel half-height and average velocity. Calculations of mean velocity profile, turbulence intensities, skewness, flatness, Reynolds stress and eddy diffusivity of heat near a wall compare favourably with experimental results. Representative contour plots of the temperature field near the wall and of the spanwise and streamwise two-point velocity correlations are given. Deficiencies are that the calculation requires many hours on a fast computer with a large high-speed memory and that the grid size in each direction for appropriate resolution is approximately proportional to the square of the Reynolds number and to the Prandtl number raised to some power greater than one.     

Passive heat transfer in a turbulent channel flow simulation using large eddy simulation based on the lattice Boltzmann method framework

In this paper, a large eddy simulation based on the lattice Boltzmann framework is carried out to simulate the heat transfer in a turbulent channel flow, in which the temperature can be regarded as a passive scalar. A double multiple relaxation time (DMRT) thermal lattice Boltzmann model is employed. While applying DMRT, a multiple relaxation time D3Q19 model is used to simulate the flow field, and a multiple relaxation time D3Q7 model is used to simulate the temperature field. The dynamic subgrid stress model, in which the turbulent eddy viscosity and the turbulent Prandtl number are dynamically computed, is integrated to describe the subgrid effect. Not only the strain rate but also the temperature gradient is calculated locally by the non-equilibrium moments. The Reynolds number based on the shear velocity and channel half height is 180. The molecular Prandtl numbers are set to be 0.025 and 0.71. Statistical quantities, such as the average velocity, average temperature, Reynolds stress, root mean square (RMS) velocity fluctuations, RMS temperature and turbulent heat flux are obtained and compared with the available data. The results demonstrate great reliability of DMRT–LES in studying turbulence.     

Passive heat transfer in a turbulent channel flow simulation using large eddy simulation based on the lattice Boltzmann method framework

In this paper, a large eddy simulation based on the lattice Boltzmann framework is carried out to simulate the heat transfer in a turbulent channel flow, in which the temperature can be regarded as a passive scalar. A double multiple relaxation time (DMRT) thermal lattice Boltzmann model is employed. While applying DMRT, a multiple relaxation time D3Q19 model is used to simulate the flow field, and a multiple relaxation time D3Q7 model is used to simulate the temperature field. The dynamic subgrid stress model, in which the turbulent eddy viscosity and the turbulent Prandtl number are dynamically computed, is integrated to describe the subgrid effect. Not only the strain rate but also the temperature gradient is calculated locally by the non-equilibrium moments. The Reynolds number based on the shear velocity and channel half height is 180. The molecular Prandtl numbers are set to be 0.025 and 0.71. Statistical quantities, such as the average velocity, average temperature, Reynolds stress, root mean square (RMS) velocity fluctuations, RMS temperature and turbulent heat flux are obtained and compared with the available data. The results demonstrate great reliability of DMRT–LES in studying turbulence.     

Application of compressible Reynolds-averaged governing equations to turbulent mixed convection in supercritical fluids in heated vertical tubes

Accurate estimation of thermal-hydraulic characteristics of supercritical flows has long been an attractive but elusive subject to many researchers in spite of tremendous effort devoted to the development of suitable turbulence models. One of the key reasons for the difficulty is a lack of measured turbulence data, which might have been used to formulate adequate turbulence models suitable for highly buoyant fluids. Turbulence models are typically based on the log-law, while the velocity profile in buoyant fluids substantially deviates from the log-law because of significant density variation in a turbulent boundary layer. In this paper, axisymmetric compressible Reynolds-Averaged governing equations were employed together with the property-dependent turbulent Prandtl number to reproduce experimental data representing heat transfer deterioration and consequential sudden temperature increase. The additional turbulence terms associated with turbulent mass flux appeared in the governing equations were modeled using the simple gradient diffusion hypothesis (SGDH). The proposed model successfully reproduced the experimental data. The various turbulence properties are presented and discussed.