充水压力管道受弹体冲击破坏的实验研究

报导了厚壁管 (外径D =4 5mm ,壁厚H =3mm)和薄壁管 (外径D =4 4mm ,壁厚H =1mm)两组充水压力管道的撞击穿透破坏实验结果 .借助落锤装置实现对三跨连续管梁撞击加载 ,冲头形状有平头、半球头和 90°圆锥头三种 .实验为内空管和管内充水加压 0MPa ,5MPa ,10MPa ,15MPa和 2 0MPa六种工况 .实验得到了不同工况下的破坏模式 ,同时记录了撞击过程中管内压力变化时程曲线及冲击力时程曲线 .实验结果表明 ,内充介质显著降低了薄壁管的穿透能量 ,而对厚壁管的降低程度较小 .随内充介质的压力增高 ,穿透能量也有不同程度的降低 .     

方管内气固两相流动速度的LDV测量

介绍了气固两相流动的ＬＤＶ测试方法和有关技术，并给出了方管内网栅后气固两相流动的部分测试结果与分析。对于气相，给出了纯净气流（α＝０）、平均粒径ｄｐ＝３５０μｍ的颗粒（α＝０．１２％、０．２１％、０．３３５％）存在时的气流在管内不同断面上时均速度分布情况及不同颗粒浓度时气流速度结构的比较；对于固相，给出了三种不同固相浓度时颗粒在管内不同断面上时均速度的分布情况，并对三种固相浓度时颗粒与气相的速度分布进行了比较。     

剪切流场中含内流立管横向涡激振动特性

立管是海洋工程中输送油气或其他矿产资源的必备结构, 外部洋流引起的立管涡激振动影响着立管的疲劳寿命, 危害深海资源开发. 本文基于欧拉?伯努利梁方程, 结合半经验时域水动力模型, 建立剪切流与内流耦合作用下海洋立管涡激振动预报模型, 运用有限元方法和Newmark-β逐步积分法求解方程, 首先将数值模拟结果与实验数据进行对比, 验证模型正确性. 然后, 运用此模型, 对剪切流作用下含内流的顶张立管在不同内流速度和密度下的横向涡激振动响应特性进行研究, 主要分析了立管的横向振动模态、振动频率以及均方根位移等涡激振动参数随内流速度和密度等参数的变化规律. 结果表明, 在剪切流场中, 含内流海洋立管在横向上表现出多模态多频率的涡激振动;立管横向振动的最大均方根位移随内流速度和密度的增大而增大, 特别是当内流速度较大时, 横向最大均方根位移增大明显;立管横向振动的主导频率随内流速度和密度的增大而减小, 并且内流密度的增大同样会引起模态转换和频率转换.      

自由面对潜艇尾流场流动特性影响研究

潜艇周围绕流场流动特性会影响潜艇的机动性能, 特别是近水面航行时, 自由面的存在会增大潜艇尾流场的复杂程度.为探究潜艇在近水面航行时自由液面对潜艇尾流场流动特性的影响机理,借助大型水下三维粒子图像测速技术开展潜艇尾流场流动特性研究.首先通过美国泰勒水池标准模型实验结果对试验方法准确性进行验证; 随后,用验证后的模型试验方法对潜艇尾流场进行测量,得到不同潜深工况、不同速度下的桨盘面轴向速度以及脉动速度,同时辅以数值模拟对试验无法测得的兴波波系及中纵剖面速度场加以补充,从兴波角度阐述了自由面对潜艇尾流场流动特性影响机理. 研究结果表明:潜艇在近水面航行时, 随着Fr增大,桨盘面处轴向速度云图中上方等值线整体趋于扁平化, 较4D潜深工况,1.5D潜深工况出现局部脉动速度极大值, 且脉动速度结构整体下移; 自由面存在时,艇体与自由面间流场速度明显增大, 特别在桨盘面区域, 流场速度明显提升.随着Fr增大, 桨盘面处的自由液面高度逐渐降低,这就导致了桨盘面位置出现更大的流体速度, 即造成了桨盘面伴流场挤压现象.      

基于DES方法平底水力旋流器内部流场特征的数值分析

采用分离涡模拟方法(Detached Eddy Simulation,DES)针对平底水力旋流器的特殊结构进行了清水流场三维湍流流动的数值模拟.与RNG k-ε模型、雷诺应力模型的计算结果比较后发现,DES方法所预测出的切向速度结果与实验结果吻合最好,从而验证了DES方法在水力旋流器数值模拟方面的可行性.进一步分析了在一定入口流量下平底水力旋流器内部流场的切向速度、轴向速度、径向速度及静压力的分布特征以及局部流动现象,获得了有关平底水力旋流器内流场流动结构的丰富的信息,讨论了圆柱结构内的流动特征对固液分离可能存在的影响.分析表明平底水力旋流器的内流场呈现出兰金涡的特点,并且为非对称分布,其内部的能量损耗主要集中在受强制涡影响的中心轴线附近区域内.本文结果为进一步研究平底水力旋流器的流场特征提供了参考.     

自由面对潜艇尾流场流动特性影响研究

潜艇周围绕流场流动特性会影响潜艇的机动性能,特别是近水面航行时,自由面的存在会增大潜艇尾流场的复杂程度.为探究潜艇在近水面航行时自由液面对潜艇尾流场流动特性的影响机理,借助大型水下三维粒子图像测速技术开展潜艇尾流场流动特性研究.首先通过美国泰勒水池标准模型实验结果对试验方法准确性进行验证;随后,用验证后的模型试验方法对潜艇尾流场进行测量,得到不同潜深工况、不同速度下的桨盘面轴向速度以及脉动速度,同时辅以数值模拟对试验无法测得的兴波波系及中纵剖面速度场加以补充,从兴波角度阐述了自由面对潜艇尾流场流动特性影响机理.研究结果表明:潜艇在近水面航行时,随着Fr增大,桨盘面处轴向速度云图中上方等值线整体趋于扁平化,较4D潜深工况, 1.5D潜深工况出现局部脉动速度极大值,且脉动速度结构整体下移;自由面存在时,艇体与自由面间流场速度明显增大,特别在桨盘面区域,流场速度明显提升.随着Fr增大,桨盘面处的自由液面高度逐渐降低,这就导致了桨盘面位置出现更大的流体速度,即造成了桨盘面伴流场挤压现象.     

膨胀管残余应力测试

膨胀管残余应力是影响膨胀管抗外挤、内压强度的重要因素之一。确定膨胀管的残余应力,测试方法至关重要。在常用测试方法中,机械测量法中的盲孔法体现出测试的优越性。应用盲孔法测试J55、K55两种规格膨胀管的残余应力,并根据实际测试的应力水平和实验误差,采用应变释放系数分级法对试验结果进行修正,以提高测试结果的精度。实验结果表明,膨胀管膨胀后存在环向、轴向的残余应力,且外表面总是压残余应力。将实验结果与仿真结果对比,误差满足工程要求,说明盲孔法测试膨胀管残余应力方便、快捷、准确、有效。     

倾斜非均匀磁场下导电方管磁流体管流的数值模拟研究

热核聚变反应堆液态金属包层应用中的一个重要问题是液态金属在导电管中流动和强磁场相互作用产生的额外的磁流体动力学压降.这种磁流体动力学压降远远大于普通水力学压降.美国阿贡国家实验室ALEX研究小组,对非均匀磁场下导电管中液态金属磁流体动力学效应进行了实验研究,其实验结果成为液态金属包层数值验证的标准模型之一.液态金属包层在应用中会受到不同方向的磁场作用,本文以ALEX的非均匀磁场下导电方管中液态金属管流实验中的一组参数为基础,保持哈特曼数、雷诺数和壁面电导率不变,采用三维直接数值模拟的方法,研究了外加磁场与侧壁之间的倾角对导电方管内液态金属流动的速度、电流和压降分布的影响.研究结果表明:沿流向相同横截面上的速度、电流以及压力分布均随磁场的倾斜而同向旋转.倾斜磁场均匀段,横截面上的高速区位于平行磁场方向的哈特曼层和平行层交叉位置,压力梯度随磁场倾角的增大先增大后减小.倾斜磁场递减段,在三维磁流体动力学效应作用下,横截面上的高速射流位置向垂直磁场方向偏移.磁场递减段的三维磁流体动力学压降随磁场倾角的增大而增大.随磁场倾斜,截面上的射流峰值逐渐减小,二次流增强,引发层流向湍流的转捩.      

水平管气液两相间歇流含气率研究

对大直径水平管内气液二相流动进行了实验，实验结果表明大直径水平管与小直径水平管具有不同的气液二相流动特征。分析和提出了适用于大直径水平管内间歇流液弹含气率模型，其计算值与实验值非常接近。     

倾斜非均匀磁场下导电方管磁流体管流的数值模拟研究

热核聚变反应堆液态金属包层应用中的一个重要问题是液态金属在导电管中流动和强磁场相互作用产生的额外的磁流体动力学压降.这种磁流体动力学压降远远大于普通水力学压降.美国阿贡国家实验室ALEX研究小组,对非均匀磁场下导电管中液态金属磁流体动力学效应进行了实验研究,其实验结果成为液态金属包层数值验证的标准模型之一.液态金属包层在应用中会受到不同方向的磁场作用,本文以ALEX的非均匀磁场下导电方管中液态金属管流实验中的一组参数为基础,保持哈特曼数、雷诺数和壁面电导率不变,采用三维直接数值模拟的方法,研究了外加磁场与侧壁之间的倾角对导电方管内液态金属流动的速度、电流和压降分布的影响.研究结果表明:沿流向相同横截面上的速度、电流以及压力分布均随磁场的倾斜而同向旋转.倾斜磁场均匀段,横截面上的高速区位于平行磁场方向的哈特曼层和平行层交叉位置,压力梯度随磁场倾角的增大先增大后减小.倾斜磁场递减段,在三维磁流体动力学效应作用下,横截面上的高速射流位置向垂直磁场方向偏移.磁场递减段的三维磁流体动力学压降随磁场倾角的增大而增大.随磁场倾斜,截面上的射流峰值逐渐减小,二次流增强,引发层流向湍流的转捩.     

Numerical study on mitigating severe slugging in pipeline/riser system with wavy pipe

The gas/liquid two-phase flow in pipeline/wavy-pipe/riser systems was investigated numerically with CFD. A CFD model of the pipeline/wavy-pipe/riser system was obtained by adding a wavy pipe to the model of the pipeline/riser system verified by the experimental data previously. The effects of the geometrical parameters and location of the wavy pipe on its performance of slug mitigation and flow characteristics in pipeline/wavy-pipe/riser systems were examined through the CFD models. With the increase of the amplitude or length of the wavy pipe, the slug in the pipeline/riser system becomes shorter. The optimum location of the wavy pipe in the pipeline exists for a pipeline/riser system and a wavy pipe at given operating conditions. The CFD modelling provides a feasible and flexible way to investigate the effectiveness of the wavy pipes on mitigating severe slugging in pipeline/riser systems.     

含颗粒污染物的油液均质流的数值分析

含颗粒污染物的油液是浓度很稀、粒径极小的伪均质流,为了掌握颗粒污染物在输送过程中的浓度分布,利用一维扩散方程构建了污染油液的数学模型;通过特征线法数值求解,获得了污染油液中各相的动态特征。结果表明,油液压力和速度沿管长呈脉动规律运动,且随着时间的延长逐渐衰减;颗粒污染物对油液速度具有极好的跟随特性;颗粒污染物的浓度分布也随着油液流速的变化而呈现规律性的变化;在不同运行时间内油液压力沿管长的衰减趋势不同,油液速度沿管长的变化趋势与压力的趋势相反;颗粒污染物速度和浓度分布沿管长与油液速度具有紧密的联系。     

CFD modelling of the cross-flow through normal triangular tube arrays with one tube undergoing forced vibrations or fluidelastic instability

A CFD methodology involving structure motion and dynamic re-meshing has been optimized and applied to simulate the unsteady flow through normal triangular cylinder arrays with one single tube undergoing either forced oscillations or self-excited oscillations due to damping-controlled fluidelastic instability. The procedure is based on 2D URANS computations with a commercial CFD code, complemented with user defined functions to incorporate the motion of the vibrating tube. The simulation procedure was applied to several configurations with experimental data available in the literature in order to contrast predictions at different calculation levels. This included static conditions (pressure distribution), forced vibrations (lift delay relative to tube motion) and self-excited vibrations (critical velocity for fluidelastic instability). Besides, the simulation methodology was used to analyze the propagation of perturbations along the cross-flow and, finally, to explore the effect on the critical velocity of the Reynolds number, the pitch-to-diameter ratio and the degrees of freedom of the vibrating cylinder.     

Computational Fluid Dynamics (CFD) modelling of transfer chutes: Assessment of viscosity,drag and turbulence models

Computational Fluid Dynamics (CFD) has been successfully applied to evaluate potential dust emissions from bulk material transfer chutes. The implementation of appropriate models and modelling parameters is shown to be critical to the overall accuracy of the simulation results. This paper presents the influence of different models on the CFD simulation of transfer chutes and follows from an earlier study that details the influence of model parameters. The aim of this paper is to offer guidance to select models, and provide a better understanding of their influence in order to evaluate the most appropriate viscosity model, drag model and turbulence model for this application. A two-phase three-dimensional Euler–Euler model in commercial CFD software ANSYS FLUENT has been selected to model the granular and air flow in the transfer chute. The simulated air velocity profiles are discussed by comparing with each other and against experimental data obtained from Particle Image Velocimetry (PIV) results. The simulated particle velocity distributions were compared with results obtained using a well-established continuum method which was developed by Roberts. Furthermore, the air mass flow rates were analysed to evaluate the influence of different models. The results show that the granular viscosity model had a strong influence on the predictions of both air and particle flow. It was observed that both the drag model and turbulence model had limited influence on the outlet air velocities. The results also indicate that the Di-Felice drag model and SST k–ω turbulence model provide solutions closer to the experimental values than the other models investigated.     

Characteristics of developing vertical bubbly flow under normal and microgravity conditions

The axial development of the void fraction profile, interfacial area concentration and Sauter mean bubble diameter of adiabatic nitrogen-water bubbly flows in a 9 mm-diameter pipe were measured using stereo image processing under normal and microgravity conditions. The flow measurements were performed at four axial locations (axial distance from the inlet, z normalized by the pipe diameter, D, z/D = 5, 20, 40 and 60) and with various flows: superficial gas velocity of 0.00840-0.0298 m/s, and superficial liquid velocity of 0.138-0.914 m/s. The effect of gravity on radial distribution of bubbles and the axial development of two-phase flow parameters is discussed in detail based on the obtained database and visual observation. Following Serizawa-Kataoka’s phase distribution pattern criteria under normal gravity conditions, the phase distribution pattern map was developed. Similar to normal gravity two-phase flows, wall, core and intermediate void peak patterns are observed under microgravity conditions but a transition void distribution pattern is not observed in the current experimental conditions. The data obtained in the current experiment are expected to contribute to the benchmarking of CFD simulation of phase distribution pattern and interfacial area concentration in forced convective pipe flow under microgravity conditions.     

Performance analysis of wick-assisted heat pipe solar collector and comparison with experimental results

The performance of heat pipe solar collector is investigated theoretically and experimentally. The system employs wick-assisted heat pipe for the heat transfer from the absorber (evaporator) to a heat exchanger (condenser). The heat pipe is made with a copper tube and the evaporator section is finned with aluminium plate. Theoretical model predicts the outlet water from heat exchanger, heat pipe temperature and also the thermal efficiency of solar collector. The results are compared with experimental data.     

Modeling fluid-particle interaction in dilute-phase turbulent liquid-particle flow simulation

The present work examines the predictive capability of a two-fluid CFD model that is based on the kinetic theory of granular flow in simulating dilute-phase turbulent liquid-particle pipe flows in which the inter-stitial fluid effect on the particle fluctuating motion is significant.The impacts of employing different drag correlations and turbulence closure models to describe the fluid-particle interactions(i.e.drag force and long-range interaction)are examined at both the mean and fluctuating velocity levels.The model pre-dictions are validated using experimental data of turbulent liquid-particle flows in a vertical pipe at different particle Reynolds numbers(ReP > 400 and ReP < 400),which characterize the importance of the vortex shedding phenomenon in the fluid-phase turbulence modulation.The results indicate that(1)the fluctuating velocity level predictions at different ReP are highly sensitive to the drag correlation selec-tion and(2)different turbulence closure models must be employed to accurately describe the long-range fluid-particle interaction in each phase.In general,good agreement is found between the model predic-tions and the experimental data at both the mean and fluctuating velocity levels provided that appropriate combinations of the drag correlation and the turbulence closure model are selected depending on Rep.     

Asymptotic solutions of the flow of a Johnson-Segalman fluid through a slowly varying pipe using double perturbation strategy

A double perturbation strategy is presented to solve the asymptotic solutions of a Johnson-Segalman (J-S) fluid through a slowly varying pipe. First, a small parameter of the slowly varying angle is taken as the small perturbation parameter, and then the second-order asymptotic solution of the flow of a Newtonian fluid through a slowly varying pipe is obtained in the first perturbation strategy. Second, the viscoelastic parameter is selected as the small perturbation parameter in the second perturbation strategy to solve the asymptotic solution of the flow of a J-S fluid through a slowly varying pipe. Finally, the parameter effects, including the axial distance, the slowly varying angle, and the Reynolds number, on the velocity distributions are analyzed. The results show that the increases in both the axial distance and the slowly varying angle make the axial velocity slow down. However, the radial velocity increases with the slowly varying angle, and decreases with the axial distance. There are two special positions in the distribution curves of the axial velocity and the radial velocity with different Reynolds numbers, and there are different trends on both sides of the special positions. The double perturbation strategy is applicable to such problems with the flow of a non-Newtonian fluid through a slowly varying pipe.     

适用于浸入边界法的大涡模拟紊流壁面模型

基于近壁定常剪切应力假设,提出了一种新的适用于浸入边界法的大涡模拟紊流壁面模型。通过引入壁面滑移速度,修正了线性速度剖面计算得到的壁面剪切应力,使之满足Werner-Wengle模型。将其应用于平板紊流和高Re数圆管紊流的数值模拟,对比采用和不采用壁面模型的结果得知,采用此模型的速度剖面与实验值吻合良好,验证了此模型的有效性。研究了不同欧拉/拉格朗日网格相对位置对结果的影响,证明了此模型具有较好的鲁棒性,以及可根据局部流动状态和网格精度自动开闭的特点。     

掺氢比对甲烷-氧气爆轰特性的影响

含氢多组分燃料由于其优良的燃烧特性逐渐成为研究关注的重点。为了对掺氢燃料的爆轰特性作进一步的研究，设计了长3 000 mm、管径30 mm的圆柱形半封闭燃烧室，对不同初压下的CH₄-2O₂、6CH₄-H₂-12.5O₂、3CH₄-H₂-6.5O₂（掺氢比分别为0%、5.1%、9.5%）3种预混合气的爆轰特性进行了实验研究，并采用烟熏膜、离子探针和压力传感器分别探测胞格结构、火焰位置和内部压力。结果表明，甲烷/氧气掺氢后可以有效提高爆轰波的传播速度，且掺氢浓度越高，传播速度越快；同时，氢气的掺入可减少管道出口处的速度亏损并在初始压力较低时加速火焰和激波的耦合，降低胞格尺寸，提高爆轰敏感性。