多孔介质壁面剪切湍流速度时空关联的研究

作为一个基础统计量,时空关联函数在湍流问题的研究中有着广泛的应用,是研究湍流噪声、湍流中物质扩散和大涡模拟亚格子模型等问题的重要参考.本文通过建立三维多孔结构壁面剪切湍流模型,采用含Darcy-Brinkman-Forchheimer作用力项的格子Boltzmann方程对无穷大多孔介质平行板之间壁湍流进行了数值模拟,进而研究其速度脉动时空关联函数的统计特性.一方面,根据计算得到的流场数据,对比分析了常规槽道湍流与多孔介质壁面槽道湍流的时间关联函数.另一方面,计算并讨论了不同孔隙率和渗透率的多孔介质壁面对速度脉动时空关联性的影响.通过研究表明:多孔结构壁面剪切湍流的时空关联函数等值线与椭圆理论相符;在研究参数范围内,多孔介质壁面的速度时空关联系数随着孔隙率增大而增大,随着渗透率增大而减小.同时发现在槽道壁面的近壁区、过渡区、对数律区和中心区等不同位置处,速度时空关联呈现较大差异性:越远离壁面位置(对数律区和中心区),其时空关联函数所呈现的关联等值线椭圆越细长,高值相关等值线越集中.多孔介质主要改变速度时空关联椭圆图像的椭圆率,说明多孔介质壁面主要影响湍流横扫速度.     

大涡模拟研究展向旋转槽道湍流流动

本文用大涡模拟预测了以不同转速做展向旋转的槽道湍流流动,统计平均的流向速度型在壁面附近与已有实验数据符合很好,在通道中部的预测差异也能给出合理解释,对比不同转速的计算结果,表明展向旋转通道的湍流应力和壁面摩擦力在压力面附近提高、在吸力面附近降低,这些高阶湍流统计量的变化规律可以结合湍流应力输运方程加以解释,漩涡识别技术显示了近壁条带结构,其形态和猝发率受旋转附加力的影响发生改变,进而影响壁面摩擦速度的数值和分布,进一步考察垂直流动方向的截面内速度分布,发现旋转引起了垂直壁面方向的流动,形成正负相间排列的流向涡对,并随着转速的增加向压力面靠近。      

出口条件对圆形湍流射流远场区自保持性的影响

报道出口条件对圆形湍流射流自保持性影响的实验研究结果. 对来自渐缩和长管两种不同结构喷嘴的射流,在相同雷诺数条件下,沿轴线进行了速度测量; 研究的统计量包括平均速度、湍流强度、高阶矩、能谱和积分尺度. 实验结果表明,渐缩喷嘴射流比长管射流发展得更快、更容易达到自保持状态. 通过对比发现,在两射流的速度(温度)场中,平均速度(温度)、湍流强度、偏斜因子和平坦度因子都存在明显的异同. 同时发现两射流的积分尺度随轴向距离的增加都成线性增长,且在渐缩喷嘴射流中增长得更快. 通过对比两射流的边界层厚度、径向与轴向湍流强度的比值、湍动能能谱图并结合前人的研究结果,对两射流湍流场所表现出的不同的统计学行为给出了合理的解释.      

TRPIV MEASUREMENT OF DRAG-REDUCTION IN THE TURBULENT BOUNDARY LAYER OVER RIBLETS PLATE

采用高时间分辨率粒子图像测速技术对沟槽壁面平板湍流边界层速度矢量场的时间序列及其统计量进行了实验测量,讨论了在同一来流速度下沟槽壁面对平均速度剖面﹑雷诺切应力及湍流强度的影响. 用流向速度分量的多尺度空间局部平均结构函数辨识壁湍流多尺度相干结构,用条件采样和相位平均技术提取壁湍流多尺度相干结构喷射和扫掠事件的脉动速度、展向涡量的二维空间拓扑形态. 结果表明,与同材料光滑壁面对比,沟槽壁面实现了10.73%的摩阻减小量;沟槽壁面湍流边界层湍流强度及雷诺切应力皆比光滑平板湍流边界层对应统计量小,说明沟槽壁面有效降低了湍流边界层内流体的脉动. 通过比较壁湍流相干结构猝发事件各脉动速度分量与展向涡量的空间分布特征,肯定了沟槽壁面的减阻效果,发现沟槽壁面通过抑制相干结构猝发事件实现减阻.     

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

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

可压缩各向同性衰减湍流直接数值模拟研究

采用五阶有限差分WENO格式直接模拟了高初始湍流Mach数的可压缩均匀各向同性湍流,主要分析了湍流的统计特性 和压缩性的影响,包括能谱特征、激波串、耗散率、标度律等. 研究表明,湍动能主要来自于速度场螺旋分量的贡献;各向同性湍流的小尺度脉动对压缩性更为敏感,并且压缩性的增强加快了湍流大 尺度脉动向小尺度脉动的湍动能输运;随着湍流Mach数的升高,胀量(压缩)耗散率所占比率也显著增长. 标度律分析表明,强可压缩湍流的横向速度结构函数仍然具有扩展自相似性;当阶数较高(p ≥ 5)时,纵向速度结构函数的扩展自相似性则不再成立. 对于压缩性较弱的湍流,与不可压缩湍流一致,横向湍流脉动的间歇性要强于纵向湍流脉动;而对于强可压缩湍流,纵向湍流脉动的 间歇性要强于横向湍流脉动.     

）IRECT NUMERICAL SIMULATION OF DECAYING COMPRESSIBLE ISOTROPIC TURBULENCE

采用五阶有限差分WENO格式直接模拟了高初始湍流Mach数的可压缩均匀各向同性湍流,主要分析了湍流的统计特性和压缩性的影响,包括能谱特征、激波串、耗散率、标度律等.研究表明,湍动能主要来自于速度场螺旋分量的贡献;各向同性湍流的小尺度脉动对压缩性更为敏感,并且压缩性的增强加快了湍流大尺度脉动向小尺度脉动的湍动能输运;随着湍流Mach数的升高,胀量（压缩）耗散率所占比率也显著增长.标度律分析表明,强可压缩湍流的横向速度结构函数仍然具有扩展自相似性;当阶数较高（p≥5）时,纵向速度结构函数的扩展自相似性则不再成立.对于压缩性较弱的湍流,与不可压缩湍流一致,横向湍流脉动的间歇性要强于纵向湍流脉动;而对于强可压缩湍流,纵向湍流脉动的间歇性要强于横向湍流脉动.     

入口扰动频率对矩形射流影响的大涡模拟

为研究入口扰动对矩形射流的影响,对不同扰动频率下的三维矩形射流进行了大涡模拟,用分布投影法求解动量方程,对压力泊松方程采用FACR直接解法.计算结果预报的速度在横向的单峰分布和速度半宽沿流向的近似线性增长与实验结果吻合良好.计算结果表明扰动频率为射流的本征频率fe=0.22时射流的卷吸作用最强,而高频扰动则抑制卷吸作用;激励频率接近本征频率时射流提前出现速度在展向的马鞍型双峰分布,更早趋于圆射流形态.进一步的流场二维和三维瞬态涡结构图谱表明:上述现象与f=0.1和f=0.22时射流近场更频繁卷起、更丰富的大涡结构存在密切关系,正是这些大涡运动不断地在横向卷吸入外界流体,实现射流在横向的迅速扩展.     

利用O型环抑制圆柱绕法流涡脱落的数值研究

圆柱绕流涡脱落诱发较大的振动和声,如何有效地抑制值得关注.利用大涡模拟技术求解了Navier-Stokes方程,得到了涡脱落频率,升力脉动幅值及平均阻力系数.计算表明二维模拟不能体现流动基本特征,三维计算与实验吻合较好.为了抑制涡脱落,在直径为D的圆柱表面装入间距为1D,直径为0.0167D的O型环.通过升力、速度谱分析以及柱向横截面流场分析可知,在光滑圆柱外表面加入O型环能诱发流体边界层分离,有效地抑制涡脱落现象,升力脉动和观测点速度脉动幅值几乎完全消失,阻力系数也略微降低,适合在实际工程中采用.     

基于Stereo-PIV技术的三维发卡涡结构定量测量研究

发卡涡是湍流相干结构研究中最为关注的内容,实现发卡涡三维结构的定量测量并进行流体动力学分析,对深入研究湍流相干结构、实现湍流精准控制等具有重要意义.本研究通过对合成射流装置进行合理控制,使得层流边界层中产生了规则的人造发卡涡结构,进而用体视图像粒子测速仪(Stereo-PIV)锁相实验技术对发卡涡结构所在的三维空间流场进行了定量测量,并得到了一个完整周期内形成的发卡涡三维结构的空间流场. 结果发现,重构所得的三维发卡涡结构质量较高, 实验技术和方案具有可行性.发卡涡结构所在空间流场情况,符合目前人们对于发卡涡、高低速条带、喷射和扫掠事件的常规认识. 此外,对近壁二次流向涡、展向涡量集中区域的展向涡头和强剪切区域、与低速喷射流体相关的汇聚流动和发散流动等有了更细致的认识.同时, 也探讨了"基于二维脉动流场的相关特征去重构发卡涡三维流场"的可行性.为进一步定量探究发卡涡结构的形成演化、不同涡结构的融合及二次诱导等壁湍流相干结构问题提供思路.      

一种基于直接计算高阶奇异积分的断裂力学双边界积分方程分析法

相对于有限元法,边界单元法在求解断裂问题上有着独特的优势,现有的边界单元法中主要有子区域法和双边界积分方程法.采用一种改进的双边界积分方程法求解二维、三维断裂问题的应力强度因子,对非裂纹边界采用传统的位移边界积分方程,只需对裂纹面中的一面采用面力边界积分方程,并以裂纹间断位移为未知量直接用于计算应力强度因子.采用一种高阶奇异积分的直接法计算面力边界积分方程中的超强奇异积分;对于裂纹尖端单元,提供了三种不同形式的间断位移插值函数,采用两点公式计算应力强度因子.给出了多个具体的算例,与现存的精确解或参考解对比,可得到高精度的计算结果.      

汶川地震诱发文家沟巨型滑坡 碎屑流基本特征及成因机制初步分析


2008年5月12日,汶川M80地震在四川省绵竹市清平乡文家沟内诱发一巨型滑坡。通过现场调查得知,滑坡前后缘高差455m,厚度20～30m,滑面为基岩层面,初始方量2750×107m3。滑体在运动中转化为碎屑流。滑坡-碎屑流总的水平运动距离为4022m,垂直运动距离为1443m,遗留的堆积物体积达5×107m3。滑坡距映秀—北川断裂仅36km,位于其下盘,地震烈度达XI度。滑坡导致文家沟中48人遇害,并形成一条完整的地震次生地质灾害链。初步分析表明滑坡启动速度快,滑坡向碎屑流转化过程明显、地点明确。碎屑流运动过程复杂,伴有强烈的“气垫效应”和“前缘气浪冲击效应”。作者认为,文家沟滑坡的高启动速度是长持时强烈地震动作用的结果,与山体的猛烈碰撞是导致滑体解体并转化为碎屑流的原因。     

Preface: Orbits around asteroid

One of the core issues in modern celestial mechanics is the orbital dynamics in the near-regime gravitational field of as- teroids, which provides deep insights into the mathematical nature of a class of nonlinear systems, and plays as a critical basis for in situ explorations of different science goals. Lots of efforts have been made to reveal the characteristics of orbital motion in the vicinity of asteroids, and to improve the skills of asteroid research in methodology.     

Effet de l'anisotropie élastique cristalline sur la distribution des facteurs de Schmid à la surface des polycristauxInfluence of crystalline elasticity anisotropy on Schmid factor distribution at the free surface of polycrystals

Experimental studies of the plasticity mechanisms of polycrystals are usually based on the Schmid factor distribution supposing crystalline elasticity isotropy. A numerical evaluation of the effect of crystalline elasticity anisotropy on the apparent Schmid factor distribution at the free surface of polycrystals is presented. Cubic elasticity is considered. Order II stresses (averaged on all grains with the same crystallographic orientation) as well as variations between averages computed on grains with the same crystallographic orientation but with different neighbour grains are computed. The Finite Element Method is used. Commonly studied metals presenting an increasing anisotropy degree are considered (aluminium, nickel, austenite, copper). Concerning order II stresses in strongly anisotropic metals, the apparent Schmid factor distribution is drifted towards small Schmid factor values (the maximum Schmid factor is equal to 0.43 instead of 0.5) and the slip activation order between characteristic orientations of the crystallographic standard triangle is modified. The computed square deviations of the stresses averaged on grains with the same crystallographic orientation but with different neighbour grains are a bit higher than the second order ones (inter-orientation scatter). Our numerical evaluations agree quantitatively with several observations and measures of the literature concerning stress and strain distribution in copper and austenite polycrystals submitted to low amplitude loadings. Hopefully, the given apparent Schmid factor distributions could help to better understand the observations of the plasticity mechanisms taking place at the free surface of polycrystals. To cite this article: M. Sauzay, C. R. Mecanique 334 (2006).     

Nonlinear behavior of matrix-inclusion composites under high confining pressure: application to concrete and mortar

This paper is devoted to a micromechanics-based simulation of the response of concrete to hydrostatic and oedometric compressions. Concrete is described as a composite made up of a cement matrix in which rigid inclusions are embedded. The focus is put on the role of the interface between matrix and inclusion which represent the interfacial transition zone (ITZ). A plastic behavior is considered for both the matrix and the interfaces. The effective response of the composite is derived from the modified secant method adapted to the situation of imperfect interfaces. To cite this article: T.H. Le et al., C. R. Mecanique 336 (2008).     

Dynamics and Breakup of Viscoelastic Liquids (A Review)

The phenomena of hydrodynamic breakup of liquid jets, drops, films, bridges, and filaments are reviewed for liquids with viscoelastic properties. The reasons for breakup are capillary instabilities, collisions with rigid obstacles, and other forms of dynamic action. The relationship between the properties of the liquids and the features of the breakup process is discussed.     

Flow stability analysis and excitation using pulsating jets

Classical flow stability applied to transition from laminar to turbulent flow may also describe the behavior of vorticity fluctuations created by a pulsating jet placed along a solid boundary. A numerical laminar flow experiment involving a pulsating jet placed along the surface of a duct with flow separation downstream, resulted in eliminating most part of the separated flow region. Applying the same approach to a turbulent flow, it was possible to develop a turbulent stability flow formulation and apply successfully turbulent pulsating jet flow separation control. To cite this article: D. Skamnakis, K. Papailiou, C. R. Mecanique 333 (2005).