建筑物周围流场格子Boltzmann仿真

基于格子BGK方程的三维15点格子Boltz-mann模型,对建筑物周围流场进行了数值仿真,对其流场的划分进行了分析,并对不同雷诺数下建筑物周围的流场进行了数值模拟计算,得到了比较满意的结果,为环境质量的评估提供了新的途径.     

并列三柱体风致干扰DDES模拟分析

采用基于SSTk-ω湍流模型的延迟分离涡模拟方法DDES(Delayed detached eddy simulation),对三并列线性布置柱体的风荷载和流场特性进行了数值模拟.为验证方法的有效性,对单个三维方柱绕流进行了计算分析,比较计算值与试验结果的风荷载特征,两者符合较好.然后采用DDES研究了三并列方柱对称布置时,间距比对中间受扰建筑风致干扰的问题.数值结果表明,间距比较小时对受扰建筑风荷载分布有较大影响.总的来看,当间距比约为3.0时,受扰建筑的影响达到最大.同时,数值分析了三方柱周围瞬态流场特征,探讨了临近建筑之间的干扰机理.结果 表明,SST-DDES方法能够较好地预测绕方柱流动分离问题,可为群体高层建筑气动相互干扰研究提供参考.     

建筑群环境风场的特性及模拟--风环境模拟研究之一

建筑风环境是研究建筑物与城市风场相互作用的一门空气动力学与大气科学的交叉学科、将空气动力学惯用的风洞实验和流场显示研究手段，直接用于城市大气边界层和建筑物周边流场变化规律的研究．充分利用数值计算的广度和深度与风洞实验相结合，相互验证、相互补充、弥补了目前实验手段尚无法实现的研究．本以北京商务中心区(CBD)为研究背景，论述了建筑风环境中大气边界层的模拟和具体建筑物风场的显示方法，并与数值计算进行对比，取得了非常一致的结果．为该地区的城市规划和环境保护，提出了科学依据．     

相邻高层建筑表面风压的数值模拟

采用数值模拟方法对一幢由主楼和配楼组成的复杂高层建筑的表面风压进行了模拟计算，考虑了相邻一幢高层建筑对风压分布的影响。计算得到的风压值还与风洞试验值作了比较。结果表明，数值模拟较好地反映了复杂高层建筑的周围流线和表面风压的分布情况；同时，由于相邻建筑物的存在，使得原有建筑的时均流线分布尤其是涡流的数目和尺度以及建筑表面的时均风压分布均发生了较明显的改变。     

高湍流度风场中两个高层建筑间的风干扰效应

来流湍流度对群体高层建筑结构间的干扰效应有很大的影响。但由于流场模拟的复杂性，以往多数干扰效应研究所用的流场均不同程度地存在湍流度低于荷载规范建议值的问题，而且多数文献主要致力于研究在郊区地貌下的建筑群体间的干扰效应，甚少详细研究高湍流度的城市中心地貌的干扰效应。本文首先用被动模拟方法较好地模拟出乎均风剖面指数为0．3的城市中心地貌而湍流度不同的4种流场，所模拟出的流场的湍流度在模型高度处分别为10％，15％，20％和26％．在此基础上用高频动态力天平技术对两个高层建筑间的静态和动态干扰效应进行了风洞试验研究。文末对城市中心地貌下两个高层建筑间的动静态干扰以及风振系数的取值提出了一些建议。     

大涡模拟建筑物对近源大气污染物扩散的影响

采用大涡模拟(LES)方法对一个由上风低源排放的烟流在孤立建筑物周围扩散的现象进行了数值模拟。模拟中,离散网格的尺度控制在100倍Kolmogorov湍流尺度的量级以内,亚格子应力由Smagorinsky模型加以模拟,并采用了随机数方法生成符合指定湍流特征的入流边界条件。数值模拟结果表明:建筑物尾流区平均浓度场的分布和大小不仅取决于局地流场的影响,在很大程度上也受到因建筑物阻挡出现的烟流抬升和马蹄涡对烟流的捕获的影响。通过对比风洞试验数据显示:大涡模拟方法能较好地再现气流绕建筑物的流动形态和建筑物周围的浓度分布,因而具有广泛的应用前景。     

群体建筑风环境的数值模拟及分析

群体建筑的布局优化及结构抗风设计都需要准确预测建筑群的风环境．针对风环境研究的必要性, 将气体流场数值模拟(CFD)方法引入建筑群风场效应的预测和研究中,侧重模拟分析了建筑群的布局改变而引起的风速及风压场的变化,获得了特定情形下较合理的气动布局;同时表明数值模拟方法可以作为风洞试验的前期预测手段,并为群体建筑的布局优化及结构抗风设计提供科学依据．     

某商业街区建筑风压及风环境数值模拟

针对上海市某商业街区建立了简化物理模型,采用离散化的数值模拟方法对该街区风场进行三维模拟研究.首先采用CAARC标准模型进行可行性验证,探讨了流体力学数值模拟的适用性和精确度.然后对该街区风场及建筑物的风场进行了数值模拟,分析了由于狭长T字型建筑布局而引起的巷道风效应.最后,讨论了街区群楼间间距及建筑布局等因素的改变对风场和风荷载的影响.研究结果表明,商业街区规划时应避免巷道风效应,并建议将行人高度风环境纳入考虑的范畴.     

带制退器膛口流场的三维数值模拟

利用三维任意坐标 系下Euler方程组和LU分解的数值方法模拟了不带膛口装置的某步枪膛口流场，计算结果与 实验结果比较吻合. 对某口径火炮带有制退器的膛口流场进行了数值模拟，其结果反映了制 退器内外流场的主要特征和冲击波波系之间的相互作用，在数值模拟的基础上分析 了3排侧 孔对制退器性能的影响. 本文提供了有关膛口流场的详细数据，将对制退器总体性能研究有 重要的参考意义.     

带制退器的膛口射流噪声数值模拟与实验研究

为了研究膛口装置对膛口噪声气动特性的影响，对带膛口制退器的某小口径武器的膛口射流噪声进行了数值模拟和实验研究。采用计算流体力学CFD (computational fluid dynamics)-计算气动声学CAA (computational aeroacoustics)耦合算法对膛口噪声进行数值模拟，即对膛口流场进行瞬态CFD模拟，获取流场数据，然后利用所得到的结果采用声学方程模拟声源信息求解声场。基于数值模拟结果，分析了膛口流场变化及噪声的指向性分布，并与实验结果进行了对比。研究表明:膛口制退器的安装改变了膛口流场结构，影响了膛口射流噪声的指向性分布。计算结果与实验结果的误差小于9%，验证了该计算方法的可行性。研究结果可为膛口射流噪声的预测及膛口制退器的设计提供一定的参考。     

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

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

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

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