弯道水流数值模拟研究

采用边界拟合坐标技术,生成边界处正交的曲线网格,在曲线坐标系中进行弯道水流数值计算。借鉴现有关于弯道水流流速分布的研究成果,对曲线坐标系平面二维浅水方程作了修正,修正后的方程组能够计入弯道环流引起的横向的动量交换,即考虑了二次流对流线弯曲的复杂水力特性的影响。通过连续弯道实例计算,验证了该模型的可靠性和实用性。     

弯道水流二维浅水模型的修正研究

借鉴有关弯道水流流速分布的研究成果,计入深度平均流速与真实流流速分布差值引起的扩散效应,在正交曲线坐标系下建立了平面二维浅水模型.采用以标准κ-ε模型为基础的曲率效应修正紊流模型模拟紊动应力项,在一定程度上考虑了流线弯曲水流紊动应力的各向异性.应用控制体积法和交错网格法离散方程,并用SIMPLEC算法求解离散方程;同时采用修正后的模型对90°弯道水流进行了数值模拟,并与原模型的计算结果及实测资料进行了比较,结果表明该模型能够有效地模拟流线弯曲水流的水力特性.     

含污染物的弱弯曲明渠弯道湍流数值模拟

建立了可模拟弯道中含污染物湍流的三维部分抛物型代数应力模型。针对左右两岸分别泄放污染物的环流非充分发展弯道流进行了计算。分析了水流结构及污染物浓度分布的特点。     

正交数值网格的生成及平面二维流场的数值模拟

提出了一种边界正交曲线网格的生成方法。在内域,求解拉普拉斯方程组生成二维正交曲线网格;在计算域边界,提出了正交曲线网格边界正交化处理的分段拉格朗日插值边界滑动法,并对控制方程N-S方程进行曲线坐标变换,使用SIMPLEC算法求解曲线坐标下的k-ε双方程湍流模型。文中以长江重庆九龙坡弯道河段正态模型为资料进行验证,计算了流场,数值模拟结果与实测结果吻合较好。     

明渠层流近床面颗粒拖曳力研究

泥沙颗粒受到的拖曳力是泥沙运动的主要驱动力,而当前应用于计算流体力学-离散颗粒法(CFD-DPM)耦合模型进行水沙运动模拟的泥沙颗粒拖曳力公式均没有考虑明渠流底床边壁作用的影响。求解不可压缩Navier-Stokes方程,对明渠层流不同雷诺数条件下床面附近不同高度处颗粒所受拖曳力进行了模拟,根据模拟结果变化规律,提出了综合考虑床面和水流惯性对标准拖曳力影响的修正拖曳力计算公式。与常用的单颗粒标准拖曳力公式和考虑遮蔽效应的多颗粒拖曳力公式相比,采用本文修正公式得到的水沙作用力更接近高精度数值解,应用于CFD-DPM输沙模拟获得的输沙结果与输沙率公式结果一致,应用分析表明输沙模拟应当采用粗糙底床边界。     

90°弯道明渠横向环流的大涡模拟研究

采用气液两相流大涡模型模拟90°弯道明渠水流的水力特性。使用有限体积法对控制方程进行离散;使用压力隐式算子分裂PISO(Pressure-Implicit with Splitting of Operators)算法求解速度与压力的耦合;采用了VOF(Volume of Fluid)法模拟自由水面。通过计算得到弯道处流速、压强系数、横向环流强度等水力参数的分布规律。分别比较弯道纵向流速与横向流速的模拟值与实验值,二者吻合较好。通过分析弯道环流的流态特征得出:因受重力作用与水流在流经弯道时发生的离心现象的影响,凸岸附近较凹岸存在较大的二次流强度的结果;凸岸与凹岸压强系数均沿程增大,当达到弯顶时,水面横比降达到最大值;弯道横向环流强度随弯道角度增大而逐渐增大,其最大值出现在明渠底部。     

曲线坐标下连通域河道平面二维水流的数值模拟

在边界拟合曲线坐标系下,运用B型交错网格模式和动边界扫描技术建立了基于连通域的二维水流数学模型,并提出了模型中有关参数的处理方法.采用贴体坐标变换将复杂的物理域变换成规则的计算域,在计算域上采用控制容积法离散方程,应用SIMPLEC算法计算速度-压力耦合.研究结果表明:采用控制容积法和SIMPLEC算法离散求解方程,具有良好的守恒性和稳定性; 该模型能够较准确地模拟连通域河段的流场变化、水位变化等过程,可供实际工程应用.     

三维势流场的比例边界有限元求解方法

比例边界有限元法(SBFEM)是线性偏微分方程的一种新的数值求解方法。该方法只对计算域边界利用Galerkin方法进行数值离散,相对于有限元方法(FEM)减少了一个空间坐标的维数,而在减少的空间坐标方向利用解析方法进行求解;相对于边界元法(BEM),比例边界有限元方法不需要基本解,避免了奇异积分的计算,所以它结合了有限元和边界元方法的优点。本文建立了利用比例边界有限元法求解三维Laplace方程的数值模型并用于计算三维物体周围的水流场,将计算结果与解析解和边界元方法进行了对比,结果表明此方法可以很好地模拟水流场,且具有较高的计算精度。     

常曲率弯道二维流动稳定性与非线性演化规律分析

国内外对弯曲河流的研究从线性特性转移到不稳定性与非线性特性上来, 弯曲河流作为一个不稳定的系统, 主要是受到水流和边界的不稳定性的影响.针对弯曲河流中存在的水流动力的不稳定性与非线性特征,利用弱非线性理论与摄动法进行展开, 在微弯条件下建立了时间模式下的常曲率 弯道二维水流扰动幅值与扰动幅角的非线性演化方程.研究了在不同弯曲度下的扰动发展特征,探讨了扰动流场在弯曲度影响下的时空分布的规律,分析了扰动流场中出现的扰动漩涡的运动过程, 阐述了弯曲度对弯曲边界内水流不稳定性的影响, 具体表现为在微弯的情况下, 弯道的弯曲度增大会提高河湾内水流的稳定性, 扰动振幅与扰动幅角会随弯曲度的增大其衰减速度更快, 并且扰动流速的对称性在弯曲度较大时减弱, 逐渐向弯道的凸岸偏斜, 在中性状态附近 的弯道水流动力具有对流不稳定性和非线性衰减的特征.本文的研究成果为构建水流动力非线性与床面形态几何非线性相互作用的模型提供了思路.      

基于滑移速度壁模型的复杂边界湍流大涡模拟

采用滑移速度壁模型实现了浸入边界方法与壁模型相结合的大涡模拟.本文首先分别采用平衡层模型和非平衡壁模型对周期山状流进行数值模拟,以考查在壁模型中考虑切向压力梯度的作用.数值结果表明,流场的压力对本文所采用的壁模型形式并不敏感,但是考虑切向压力梯度可以显著改进壁面摩擦力的计算结果,并且能够准确的预测强压力梯度区以及分离区内的流动平均统计特性.不考虑压力梯度效应的平衡层模型显著低估了壁面摩擦力的分布,同时无法准确预测分离区内的平均速度剖面.非平衡模型的修正项正比于切向压力梯度和壁面法向距离,因此在强压力梯度区或者网格较粗时,计算得到的平均压力和摩擦力分布以及流动的低阶统计量均与参考的实验和计算结果吻合.在此基础上,通过回转体绕流的大涡模拟考查了该方法用于模拟高雷诺数壁湍流的适用性,非平衡壁模型可以准确地捕捉流动的物理结构并较准确地预测其水动力学特性.结果表明,将浸入边界方法与非平衡滑移速度壁模型相结合的大涡模拟,有望成为数值模拟复杂边界高雷诺数壁湍流的工具.      

The 7th International Conference on Fluid Mechanics May 24-27,2015,Qingdao,China

正http://www.icfm7.org First Announcement and Call for PapersThe objective of International Conference on Fluid Mechanics(ICFM)is to provide a forum for researchers to exchange new ideas and recent advances in the fields of theoretical,experimental,computational Fluid Mechanics as well as interdisciplinary subjects.It was successfully convened by the Chinese Society of Theoretical and Applied Mechanics(CSTAM)in Beijing(1987,     

INFORMATION FOR CONTRIBUTORS

Contributions： The Journal, Acta Mechanica Solida Sinica, is pleased to receive papers from engineers and scientists working in various aspects of solid mechanics. All contributions are subject to critical review prior to acceptance and publication.     

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

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

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

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