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

 在条纹沟槽表面减阻的研究发展及其优点的基础上 对条纹沟槽表面减阻进行水洞实验,研究在航行器外表面加工条纹沟槽所具有的阻 力特性. 通过对光体、0.1mmV型条纹沟槽表面和0.2mmV型条纹沟槽表面3 种不同模型在零攻角、不同水流速度下进行阻力测量, 测得条纹沟槽表面减阻特 性曲线. 从特性曲线中可以看出减阻效果约为8.3%. 实验的结果表明, 条纹沟槽 表面能明显降低水下航行器的阻力.     

水下条纹沟槽表面近壁区流场的数值分析

近年来对于水下条纹沟槽表面减阻的研究多集中在试验方面,理论分析和数值计算进行的却不多。而由于受实验条件的限制,条纹沟槽表面减阻的微观机理也一直没有搞清楚。本文正是通过理论方面的分析和计算来进一步揭示条纹沟槽表面的微观流场。首先基于“第二涡群”论对条纹沟槽表面流场的减阻特性进行了理论分析,阐述了其减阻机理;在此基础上采用B．L两层代数模型,建立了条纹沟槽表面水下流场的简化数学模型;并利用有限差分法对近壁区流场进行了数值计算．首次通过理论计算的方法获得了条纹沟槽表面水下微观流场的速度分布;为以后进一步从理论和实验上研究条纹沟槽减阻特性奠定了重要的基础。     

脊状表面航行器模型减阻特性的风洞实验研究

针对航行器提高航程和航速的需要,开展脊状表面湍流边界层减阻的实验和数值仿真研究。在航行器模型的外表面加工具有特定形状、尺寸的脊状结构,导致湍流边界层的流动稳定性增强,壁面摩擦阻力降低。在风洞中对具有光滑表面和脊状表面的航行器模型在不同风速和攻角下进行阻力测试,得到其减阻特性曲线。实验结果表明,具有横向脊状表面的航行器模型在一定来流速度范围内具有很好的减阻效果,实验获得的最大减阻量为23.5%。数值仿真结果则发现,在脊状结构内形成了稳定的"二次涡",边界层内湍动能和湍流猝发强度降低,很好地揭示了减阻机理。     

超疏水沟槽表面通气减阻实验研究

减阻是解决航行体提速和增程的主要技术途径之一, 对缓解日益严峻的能源危机极为重要. 在重力式管道实验系统中, 测试给出了湍流状态下不同通气速率时减阻率随雷诺数及沟槽无量纲间距的变化规律和气膜铺展状态, 对比分析了单纯超疏水表面与超疏水沟槽表面上通气时减阻效果的差异.实验板材质为无色亚克力, 沟槽结构采用机械方法加工, 并在表面喷涂超疏水涂层. 结果表明, 持续通气能解决超疏水沟槽表面气膜层流失问题, 实现气膜层长时间稳定维持; 恒定雷诺数下, 随通气速率增大, 超疏水沟槽表面气膜铺展更趋均匀, 减阻率上升; 由于通气速率影响气膜横向扩展能力, 致使恒定通气速率下, 减阻率随雷诺数的变化呈现两种模式; 在固定雷诺数及通气速率时, 减阻率随沟槽尺寸的扩大先增后减, $S^{+}\approx 76$时减阻率最大. 分析其原因在于, 沟槽结构增大沾湿面积的同时, 显著提升了通气状态下超疏水表面气膜层的稳定性, 因而展示出与超疏水表面和沟槽表面均不相同的减阻规律, 且效果更佳.      

超疏水沟槽表面通气减阻实验研究

减阻是解决航行体提速和增程的主要技术途径之一,对缓解日益严峻的能源危机极为重要.在重力式管道实验系统中,测试给出了湍流状态下不同通气速率时减阻率随雷诺数及沟槽无量纲间距的变化规律和气膜铺展状态,对比分析了单纯超疏水表面与超疏水沟槽表面上通气时减阻效果的差异.实验板材质为无色亚克力,沟槽结构采用机械方法加工,并在表面喷涂超疏水涂层.结果表明,持续通气能解决超疏水沟槽表面气膜层流失问题,实现气膜层长时间稳定维持;恒定雷诺数下,随通气速率增大,超疏水沟槽表面气膜铺展更趋均匀,减阻率上升;由于通气速率影响气膜横向扩展能力,致使恒定通气速率下,减阻率随雷诺数的变化呈现两种模式;在固定雷诺数及通气速率时,减阻率随沟槽尺寸的扩大先增后减, S~+≈76时减阻率最大.分析其原因在于,沟槽结构增大沾湿面积的同时,显著提升了通气状态下超疏水表面气膜层的稳定性,因而展示出与超疏水表面和沟槽表面均不相同的减阻规律,且效果更佳.     

回转体高速倾斜入水的流场特性及结构响应

回转体高速入水过程涉及液体和固体的耦合作用，是一个复杂的非线性、非定常过程。为研究回转体高速入水的结构动响应及流场演变规律，本文中基于STAR-CCM+和ABAQUS平台，建立了回转体高速入水的双向流固耦合数值模型，开展了不同入水速度的回转体高速倾斜入水流固耦合数值计算。结果表明:数值计算的入水速度、位移曲线和空泡形态与实验结果良好吻合，验证了流固耦合方法的有效性；回转体倾斜高速入水的载荷先集中在触水部分边缘处，后集中于回转体底部中心处；流固耦合方法的入水冲击载荷峰值小于刚体的，弹性回转体的载荷曲线产生明显波动；撞水阶段，回转体空泡呈现不对称形态，随着入水加深，空泡不对称性变弱；入水速度60 m/s下，空泡发生表面闭合，回转体入水初速度越快，空泡表面闭合越晚；冲击载荷与入水速度有关，入水速度越大，峰值出现越早，震荡越明显，速度超过100 m/s时，回转体产生塑性形变。     

湍流流动中鲨鱼皮表面流体减阻研究进展

快速游动的鲨鱼, 其皮肤表面沿流动方向有序地排列着沟槽状结构, 人们认为这种结构能在湍流流动 中减小表面摩擦阻力. 人们仿真这种生物结构来进行实验研究和应用, 通过复制和改善鲨鱼皮肤表面沟槽状 结构, 使得摩擦阻力最大减小了近10%. 在实验和模拟仿真中, 人们不断讨论和研究湍流流动阻力的形成机制 和沟槽减阻的理论特性. 本文综述了沟槽减阻理论特性的一些研究方法, 并且归纳定义了沟槽减阻最优几何 形状及其尺寸; 详细考虑流体流动的特点, 给出了一种用来选择最优沟槽形状及其尺寸的方法; 综述了目前的 沟槽加工制造技术. 由于鲨鱼皮肤表面存在少量黏液, 从仿生学的角度, 文章最后综述并展望了通过局部应用 疏水性材料来改变沟槽附近流场属性, 从而达到更大程度上减小阻力的目标.      

脊状表面减阻特性的风洞试验研究

利用热线风速仪,对光滑表面和多个脊状表面在低速风洞中进行了表面流场测试。基于测得的边界层速度分布数据,利用对数律区速度分布公式,编程分别计算出光滑表面和脊状表面的壁面摩擦速度和虚拟原点。研究发现,脊状表面最大减阻量达13．5％;有减阻效果的脊状表面使边界层速度曲线上移、湍流强度下降;与光滑表面相比,脊状表面的位移厚度和动量损失厚度明显减小,也表明脊状表面具有减阻效果;位移厚度和动量损失厚度减少量随槽间距s^＋的增加呈现先变大后变小的趋势,在S^＋=12时达到最大。     

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

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

翼型表面脊状结构减阻特性的数值仿真研究

对翼型表面脊状结构进行了数值仿真研究,初步发现了翼型表面脊状结构的减阻特性,并进一步研究了间隔对于翼型表面脊状结构减阻效果的影响规律。针对翼型表面脊状结构的特点,在数值计算过程中对其计算域、计算网格及流动参数进行了合理的处理。结果表明:脊状结构在翼型表面仍然具有很好的减阻效果,其表面边界层内湍流强度、湍动能明显低于光滑表面;有间隔的脊状结构减阻效果优于无间隔的脊状结构,间隔大小与脊状结构尺寸相当时,减阻效果最佳。     

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

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

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


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).