汶川地震滑坡与地震参数及地质地貌因素之间的相关关系

采用大涡模拟方法,模拟了槽道湍流,得到了不同雷诺数下槽道湍流的结果. 在此基 础上,研究了平均速度、雷诺应力、脉动动能和脉动速度均方根的分布;讨论了平均速度的 壁面律问题;给出了雷诺应力、脉动动能和脉动速度均方根随雷诺数的变化规律,其中雷诺 应力、脉动动能给出了定量公式.     

壁湍流结构的二维子波分析

对槽道湍流直接数值模拟的脉动速度场进行了二维子波分析。研究了近壁区低速条带结构的多尺度特性。利用定义的局部雷诺应力测度结合平坦因子的方法对湍流脉动速度分解,得到脉动速度的间歇成份和高斯成份。用局部雷诺应力测度正的峰值为特征,通过条件统计平均,得到近壁区典型的拟序结构。     

雷诺应力各向异性涡黏模型的层析TRPIV测量

利用层析TRPIV测量水洞中平板湍流边界层3D-3C速度场的高分辨率时间序列数据库. 提出了空间局部平均多尺度速度结构函数的新概念, 描述湍流多尺度涡结构的空间拉伸、压缩、剪切变形和旋转. 用空间局部平均多尺度速度结构函数对湍流脉动速度进行了空间多尺度分解. 用空间流向局部平均多尺度速度结构函数, 根据湍流多尺度涡结构在流向的拉伸和压缩物理特征, 提出了新的湍流相干结构条件采样方法, 检测并提取了层析TRPIV数据中相干结构“喷射”和“扫掠”事件中的脉动速度、平均速度变形率、雷诺应力等物理量的空间拓扑形态. 通过研究平均速度变形率各分量与雷诺应力各分量之间的空间相位差异,肯定了壁湍流相干结构雷诺应力各向异性复涡黏模型的合理性.      

湍流边界层相干结构剩余脉动雷诺应力项的实验测量

用IFA300恒温热线风速仪和X形二分量热线探针以采样间隔小于湍流耗散时间尺度的分辨率精细测量了风洞中平板湍流边界层不同法向位置的瞬时流向、法向速度分量的时间序列信号。用条件采样和相位平均技术提取了相干结构猝发过程中相干结构剩余脉动雷诺应力和随机脉动对相干结构贡献的雷诺应力的条件相位平均波形。基于理论上对湍流相干结构动量方程中随机脉动对相干结构贡献的雷诺应力和相干结构剩余脉动雷诺应力项的分析,对两种雷诺应力项进行了对比研究。研究发现,相干结构剩余脉动雷诺应力项在数值上具有和随机脉动对相干结构贡献的雷诺应力相同的数量级,表明在相干结构动力学模型方程中,相干结构剩余脉动雷诺应力项并不像以前估计的那样可以忽略不计。     

各向同性湍流通过正激波的演化特征研究

激波与湍流相互作用(shock-turbulence interaction,STI)是空气动力学研究中的一个基础问题.基于格心有限差分法(cell-centered finite difference method,CCFDM)求解器Helios,采用五阶加权紧致非线性格式(weighted compact nonlinear scheme,WCNS)对各向同性湍流通过正激波的情形进行直接数值模拟(direct numerical simulation,DNS).对湍流相关物理量进行统计,分析结果表明,在湍流中波后的密度、温度和压力较无湍流情形下略小,而速度则略大,均在波后呈现短暂过冲然后缓慢向理论值逼近的变化趋势;波后流向雷诺应力突降随之快速增长又衰减,呈现非单调变化趋势,线性相互作用分析(linear interaction analysis,LIA)将其归结为波后能量从声模式转移为涡模式方式,与流向不同,横向雷诺应力突增后单调衰减,波后雷诺应力各向异性明显且随下游距离逐渐增强;波后湍动能突增后呈现非单调变化趋势;泰勒微尺度和Kolmogorov尺度过激波后均明显减小,说明波后湍流长度尺度变小,从而对波后网格的分辨率提出了更高的要求;密度、温度和压力过激波后脉动均方根均增加,密度和压力脉动强度减小,温度脉动强度增大.      

各向同性湍流通过正激波的演化特征研究

激波与湍流相互作用(shock-turbulence interaction,STI)是空气动力学研究中的一个基础问题.基于格心有限差分法(cell-centered finite difference method,CCFDM)求解器Helios,采用五阶加权紧致非线性格式(weighted compact nonlinear scheme,WCNS)对各向同性湍流通过正激波的情形进行直接数值模拟(direct numerical simulation,DNS).对湍流相关物理量进行统计,分析结果表明,在湍流中波后的密度、温度和压力较无湍流情形下略小,而速度则略大,均在波后呈现短暂过冲然后缓慢向理论值逼近的变化趋势;波后流向雷诺应力突降随之快速增长又衰减,呈现非单调变化趋势,线性相互作用分析(linear interaction analysis,LIA)将其归结为波后能量从声模式转移为涡模式方式,与流向不同,横向雷诺应力突增后单调衰减,波后雷诺应力各向异性明显且随下游距离逐渐增强;波后湍动能突增后呈现非单调变化趋势;泰勒微尺度和Kolmogorov尺度过激波后均明显减小,说明波后湍流长度尺度变小,从而对波后网格的分辨率提出了更高的要求;密度、温度和压力过激波后脉动均方根均增加,密度和压力脉动强度减小,温度脉动强度增大.     

基于人工神经网络的亚格子应力建模

亚格子(SGS)应力建模在湍流大涡模拟(LES)中有着极为重要的作用. 传统亚格子应力模型存在相对误差较大、耗散过强等问题. 近年来, 计算机技术的发展使得人工神经网络(ANN)等机器学习方法逐渐成为亚格子应力建模型的新研究范式. 本文着重考虑滤波宽度及雷诺数影响, 在不可压缩槽道湍流中建立了亚格子应力的ANN模型. 该模型以滤波后的直接数值模拟(fDNS)流场物理量及滤波尺度为输入信息, 相应滤波尺度下的亚格子应力为输出量. 通过对不同滤波尺度及不同雷诺数数据的训练, ANN模型能够给出与直接数值模拟(DNS)高度吻合的亚格子应力. 此外, 模型在亚格子耗散等非ANN建模量上也有着优异的预测性能, 与基于DNS获得的对应物理量的相关系数大都在0.9以上, 较梯度模型及Smagorinsky模型有明显提升. 在后验测试中, ANN模型对流向平均速度剖面的预测同样优于梯度模型、Smagorinsky模型及隐式大涡模拟(ILES)等传统LES模型. 在脉动速度均方根预测方面, 除了某些法向位置外ANN模型的性能整体上相对其他3个模型有所提升. 然而, 随着网格尺度的增大ANN模型预测的结果与fDNS结果的偏差逐渐增大. 总之, ANN方法在发展高精度亚格子应力模型上具有很大的潜力.      

体育场环状悬挑屋盖脉动风压数值模拟

本文提出了预测悬挑屋盖上脉动风压的实用方法,用数值模拟中稳态计算得到的平均压力、平均速度和湍动能的分布计算大跨屋盖的脉动风压。以昆山体育场为例,通过数值模拟的平均压力系数、湍动能分布及脉动风压与风洞实验结果的比较,验证了用数值模拟方法预测环状悬挑屋盖脉动风压的有效性,并建议了需选用的合适湍流模型和参考点的风速。     

壁湍流猝发过程中速度分量的相位差对雷诺应力影响的实验研究

用IFA300恒温热线风速仪和×形二分量热线探针,以采样间隔小于最小湍流时间尺度的分辨率,精细测量了风洞中平板湍流边界层不同法向位置的瞬时流向、展向速度分量的时间序列信号.用子波分析辨识壁湍流相干结构猝发事件的能量最大准则,确定壁湍流相干结构猝发事件的时间尺度;用条件相位平均技术提取了相干结构猝发过程中流向、展向脉动速度分量条件相位平均波形,用互相关方法研究了相干结构猝发过程中流向、展向脉动速度分量条件相位平均波形的相位差关系及其对雷诺应力的影响,发现在缓冲层和对数律区,展向脉动速度与流向脉动速度的条件相位平均波形具有不同的相位;当两者相位基本一致时,雷诺应力达到正的最大值,此时湍流相干结构的产生非常活跃;当两者相位差分别集中在90°和270°附近时,雷诺应力的幅值减小并接近于零,此时湍流相干结构的产生和猝发都得到了抑制.     

矩形槽道内表面活性剂减阻流体流场特性

采用粒子图像测速仪对矩形槽道内表面活性减阻流体在流动方向(x方向)与壁面垂直方向(y方向)所在平面的流场进行了测量,分析了速度、涡量、速度脉动相关量在流场内的瞬态分布,以及对500幅相同工况的流场进行了统计平均. 结果显示: 与牛顿流体相比, 表面活性剂减阻流体接近于层流流动,横向速度脉动被大幅减弱,导致湍流输运减弱,雷诺应力远远小于水. 减阻流体流向速度脉动呈条带特征,沿流动方向发展,反映了减阻流体不同于水的湍流输运特征.      

Research on particle dispersion in a plane mixing layer with coherent structures

The numerical simulation with two-way coupling was performed in a liquid -particle mixing layer and the corresponding experiment study was made. In the process of vortex rolling up and vortices pairing, the particles with different St number have a very different pattern of dispersion. The mean velocity of particle with St = 1 is higher than that of the fluid phase on the low-speed side, and lower than that of the fluid phase on the high-speed side. The RMS of particle approaches that of the fluid phase with decreasing particle St number. The RMS in the transverse direction is smaller than that in the streamwise direction. The velocity fluctuation correlation of particle is smaller than the Reynolds shear stress, the “overshoot“ phenomenon that the velocity fluctuation correlation of particle is larger than the Reynolds shear stress does not appear. The larger the St number of particle is, the wider the range of the particle dispersion will be. The computed results are in agreement with the experimental ones.     

Velocity distribution of flow with submerged flexible vegetations based on mixing-length approach

By choosing a PVC slice to simulate flexible vegetation, we carried out experiments in an open channel with submerged flexible vegetation. A 3D acoustic Doppler velocimeter （micro ADV） was used to measure local flow velocities and Reynolds stress. The results show that hydraulic characteristics in non-vegetation and vegetation layers are totally different. In a region above the vegetation, Reynolds stress distribution is linear, and the measured velocity profile is a classical logarithmic one. Based on the concept of new-riverbed, the river compression parameter representing the impact of vegetation on river is given, and a new assumption of mixing length expression is made. The formula for time-averaged velocity derived from the expression requires less parameters and simple calculation, and is useful in applications.     

STUDY OF AXIAL VELOCITY IN GAS CYCLONES BY 2D-PIV, 3D-PIY, AND SIMULATION

The axial velocity distribution in a gas cyclone has been examined with two-dimensional particle image velocimetry (2D-PIV) and three-dimensional particle image velocimetry (3D-PIV) experiments in this study. Due to the limitation of 2D-PIV configuration, the contamination generated by the strong tangential velocity in the cyclone can be registered in the axial velocity detected by 2D-PIV. Efficient methods are proposed in this work to remove this contamination. The contamination-removed 2D-PIV data agree well with 3D-PIV results. The distributions of the axial velocity are also computed by the Reynolds stress model (RSM) and verified using the PIV experimental results. Reasonable agreements are obtained.     

Comparative study of lattice‐Boltzmann and finite volume methods for the simulation of laminar flow through a 4:1 planar contraction

In the present paper, a comparative study of numerical solutions for Newtonian fluids based on the lattice‐Boltzmann method (LBM) and the classical finite volume method (FVM) is presented for the laminar flow through a 4:1 planar contraction at a Reynolds number of value one, Re=1. In this study, the stress field for LBM is directly obtained from the distribution function. The calculations of the stress based on the FVM‐data use the evaluations of velocity gradients with finite differences. The stress field for both LBM and FVM is expressed in the present study in terms of the shear stress and the first normal stress difference. The lateral and axial profiles of the velocity, the shear stress and the first normal stress difference for both methods are investigated. It is shown that the LBM results for the velocity and the stresses are in excellent agreement with the FVM results. Copyright © 2004 John Wiley & Sons, Ltd.     

圆柱尾流雷诺应力与平均运动变形率的空间弛豫效应

针对圆柱尾流中沿流向存在的Karman涡街周期性涡旋结构,对湍流雷诺应力与平均运动变形率之间的空间弛豫效应进行了实验研究.在回流式水槽中,放入不同直径的圆柱模型,获得不同雷诺数下的圆柱尾流,利用二维高时间分辨率粒子图像测速(TRPIV)技术测量圆柱尾流二维瞬时速度空间分布图像的时间序列.经过数字图像处理,获得二维雷诺应...     

垂直湍流液-固流中大颗粒的相对速度

通过量纲分析和实验测量,对于垂直、局部均匀的湍流稀态液一固流中,大颗粒的相对速度,建立了无量纲参数表达式．用分析和实验相结合的方法,确定了表达式中无量纲参数的幂次及有关系数．实验中用激光多普勒分相测量技术,分别测出流体和颗粒的时均速度结果表明,大颗粒相对速度强烈依赖于流体雷诺数,当流体雷诺数较高时,其幂次渐近于1.5。     

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.     

钝体分离流动的实验研究

本文报告了楔形钝体分离流区域的湍流实验,提供了分离区内的时均速度、压力,湍流度和雷诺应力的分布,并对分离区的湍流特征进行了分析。实验表明,分离区内时均速度具有很大的横向梯度。湍流度和雷诺应力的分布曲线很相纵,它们在回流区变化较为平缓,而在混合区,当它们达到极大值之后,便以指数形式向(?)衰减。压力在回流区内变化也不大,但在混合区却具有明显的横向梯度。