低温气液两相流数值计算

低温工质两相流动在很多场合都存在,研究深冷剂在两相流动时参数的分布有着很重要的作用.文中采用一维均相流模型对液氮循环流动过程进行稳态数值模拟,得到了循环流动过程中液氮含汽率、压力和温度沿管程方向的分布情况,并对由于不同回路条件引起的不同结果进行了分析.     

汶川震区北川9.24暴雨泥石流特征研究

2008年9月24日汶川震区的北川县暴雨导致区域性泥石流发生,这次9.24暴雨泥石流灾害导致了42人死亡,对公路和其他基础设施造成严重损毁。本研究采用地面调查和遥感解译方法分析地震与暴雨共同作用下的泥石流特征,获取的气象数据用于分析泥石流起动的临界雨量条件。本文探讨了研究区泥石流起动和输移过程,并根据野外调查,分析了泥石流形成的降雨、岩石和断层作用,特别是强降雨过程与物源区对泥石流发生的作用。根据应急调查发现北川县境内暴雨诱发的泥石流72处,其分布受岩石类型、发震断层和河流等因素控制。根据对研究区震前和震后泥石流发生的临界雨量和雨强的初步分析,汶川地震后,该区域泥石流起动的前期累积雨量降低了14.8%～22.1%,小时雨强降低25.4 %～31.6%。震区泥石流起动方式主要有二种,一是由于暴雨过程形成的斜坡表层径流导致悬挂于斜坡上的滑坡体表面和前缘松散物质向下输移,进入沟道后转为泥石流过程;二是"消防水管效应"使沟道水流快速集中,并强烈冲刷沟床中松散固体物质,导致沟床物质起动并形成泥石流过程。调查和分析发现沟内堆积的滑坡坝对泥石流的阻塞明显,溃决后可导致瞬时洪峰流量特别大。研究结果表明了汶川震区已进入一个新的活跃期。因此,应该开展对汶川地震区的泥石流风险评估和监测、早期预警,采取有效的工程措施控制泥石流的发生和危害。     

多孔介质——自由流界面应力跳跃条件下流动特性解析解

对非对称多孔介质--自由流复合通道内多孔介质内部及多孔介质与自由流体界面处复杂质量、动量输运特性进行研究. 在多孔介质区采用Brinkman-extended Darcy模型并结合速度连续,剪切应力跳跃的界面条件对此复合通道内流体的传递现象进行求解,提出了考虑界面应力跳跃时非对称复合通道各区域流体运动速度及摩擦系数的解析式,分析了界面应力跳跃系数,达西数及无量纲多孔层偏心厚度对流体速度及摩擦系数的影响. 结果表明:改变界面性质可在一定条件下明显控制各区域流体速度分布;在达西数、多孔层偏心厚度一定情况下,界面应力系数的增大会使界面流速减小,而使流体摩擦系数增大,特别是界面应力系数小于0的情况下变化更明显,此时若不考虑界面应力系数则会造成较大误差. 当界面应力系数及多孔层偏心厚度均为较小负数值时,改变多孔层偏心厚度对界面速度的影响要大于改变界面应力系数的情况;而当界面应力系数及多孔层偏心厚度为较大正数值时,情况则相反. 较大达西数下,界面应力系数及多孔层偏心厚度对流体摩擦系数的影响均较大,继续减小达西数至一定程度时,界面应力系数对流体摩擦系数的影响可忽略不计而认为只与多孔层偏心厚度相关,且对较大多孔层偏心厚度更敏感.      

Open and traction boundary conditions for the incompressible Navier–Stokes equations

We present numerical schemes for the incompressible Navier–Stokes equations (NSE) with open and traction boundary conditions. We use pressure Poisson equation (PPE) formulation and propose new boundary conditions for the pressure on the open or traction boundaries. After replacing the divergence free constraint by this pressure Poisson equation, we obtain an unconstrained NSE. For Stokes equation with open boundary condition on a simple domain, we prove unconditional stability of a first order semi-implicit scheme where the pressure is treated explicitly and hence is decoupled from the computation of velocity. Using either boundary condition, the schemes for the full NSE that treat both convection and pressure terms explicitly work well with various spatial discretizations including spectral collocation and C⁰ finite elements. Moreover, when Reynolds number is of O(1) and when the first order semi-implicit time stepping is used, time step size of O(1) is allowed in benchmark computations for the full NSE. Besides standard stability and accuracy check, various numerical results including flow over a backward facing step, flow past a cylinder and flow in a bifurcated tube are reported. Numerically we have observed that using PPE formulation enables us to use the velocity/pressure pairs that do not satisfy the standard inf–sup compatibility condition. Our results extend that of Johnston and Liu [H. Johnston, J.-G. Liu, Accurate, stable and efficient Navier–Stokes solvers based on explicit treatment of the pressure term. J. Comp. Phys. 199 (1) (2004) 221–259] which deals with no-slip boundary conditions only.     

超燃冲压发动机等截面隔离段中流动的LBM模拟

采用耦合的双分布函数格子Boltzmann方法模拟超燃冲压发动机中隔离段的流动.计算结果反映了隔离段内流场的基本结构,表明LBM在隔离段的模拟中是-种具有潜力的数值计算方法.     

Effect of variable duty cycle flow pulsations on heat transfer enhancement for an impinging air jet

A series of experiments has been conducted in which a pulsed air jet is impinged upon a heated surface for the purpose of enhancing heat transfer relative to the corresponding steady air jet. Traditional variables such as jet to plate spacing, Reynolds number, and pulse frequency have been investigated. One additional flow variable – the duty cycle – representing the ratio of pulse cycle on-time to total cycle time is introduced and shown to be significant in determining the level of heat transfer enhancement. Specifically, heat transfer enhancement exceeding 50% is shown for a variety of operating conditions. In each case, the duty cycle producing the best heat transfer is shown to depend upon each of the other flow parameters. Recommendations are made for further experimentation into optimizing the duty cycle parameter for any particular application.     

叶栅可压缩流场的数值分析

将一维ＡＵＳＭ＾＋格式与三阶精度ＭＵＳＣＬ格式融合,给出了其在任意风线坐标下的二维形式,采用有限体积法,将ＡＵＳＭ＾＋格式与Ｒｕｎｇｅ－Ｋｕｔｔａ格式结合,对叶栅可压缩流场进行了数值模拟,典型算例的计算结果与有关献结果相符很好,表明ＲＫ－ＡＵＳＭ＾＋混合格式对叶栅可压缩流场进行数值模拟是可行的。     

Numerical methods for conjunctive two‐dimensional surface and three‐dimensional sub‐surface flows

Sophisticated catchment runoff problems necessitate conjunctive modeling of overland flow and sub‐surface flow. In this paper, finite difference numerical methods are studied for simulation of catchment runoff of two‐dimensional surface flow interacting with three‐dimensional unsaturated and saturated sub‐surface flows. The equations representing the flows are mathematically classified as a type of heat diffusion equation. Therefore, two‐ and three‐dimensional numerical methods for heat diffusion equations were investigated for applications to the surface and sub‐surface flow sub‐models in terms of accuracy, stability, and calculation time. The methods are the purely explicit method, Saul'yev's methods, the alternating direction explicit (ADE) methods, and the alternating direction implicit (ADI) methods. The methods are first examined on surface and sub‐surface flows separately; subsequently, 12 selected combinations of methods were investigated for modeling the conjunctive flows. Saul'yev's downstream (S‐d) method was found to be the preferred method for two‐dimensional surface flow modeling, whereas the ADE method of Barakat and Clark is a less accurate, stable alternative. For the three‐dimensional sub‐surface flow model, the ADE method of Larkin (ADE‐L) and Brian's ADI method are unconditionally stable and more accurate than the other methods. The calculations of the conjunctive models utilizing the S‐d surface flow sub‐model give excellent results and confirm the expectation that the errors of the surface and sub‐surface sub‐models interact. The surface sub‐model dominates the accuracy and stability of the conjunctive model, whereas the sub‐surface sub‐model dominates the calculation time, suggesting the desirability of using a smaller time increment for the surface sub‐model. Copyright © 2000 John Wiley & Sons, Ltd.     

纳米尺度圆柱绕流尾迹区流动形式模拟研究

采用非平衡分子动力学模拟方法,对微尺度低{Re}数下的圆柱绕流问题进行了研究,模拟结果表明:当{Re}<12时,圆柱下游形成对称、无分离的定常流;当{Re}>20时,圆柱下游形成周期性交替出现的对称涡;当12相似文献   

Unified semi‐analytical wall boundary conditions for inviscid,laminar or turbulent flows in the meshless SPH method

Wall boundary conditions in smoothed particle hydrodynamics (SPH) is a key issue to perform accurate simulations. We propose here a new approach based on a renormalising factor for writing all boundary terms. This factor depends on the local shape of a wall and on the position of a particle relative to the wall, which is described by segments (in two‐dimensions), instead of the cumbersome fictitious or ghost particles used in most existing SPH models. By solving a dynamic equation for the renormalising factor, we significantly improve traditional wall treatment in SPH, for pressure forces, wall friction and turbulent conditions. The new model is demonstrated for cases including hydrostatic conditions for still water in a tank of complex geometry and a dam break over triangular bed profile with sharp angle where significant improved behaviour is obtained in comparison with the conventional boundary techniques. The latter case is also compared with a finite volume and volume‐of‐fluid scheme. The performance of the model for a two‐dimensional laminar flow in a channel is demonstrated where the profiles of velocity are in agreement with the theoretical ones, demonstrating that the derived wall shear stress balances the pressure gradient. Finally, the performance of the model is demonstrated for flow in a schematic fish pass where both the velocity field and turbulent viscosity fields are satisfactorily reproduced compared with mesh‐based codes. Copyright © 2012 John Wiley & Sons, Ltd.