剪切流下发生涡激振动的柔性立管阻力特性研究

利用缩尺模型试验的方法研究了线性剪切流下涡激振动发生时柔性立管的阻力特性.文中基于光纤光栅应变传感器测得的模型应变信息,采用梁复杂弯曲理论计算了立管的平均阻力,继而分析了阻力系数沿管长方向和雷诺数的分布特性以及涡激振动对阻力系数的放大效应,并提出了用于估算柔性立管发生涡激振动时阻力系数的经验公式.结果表明:涡激振动对阻力系数有放大效应,使得立管局部阻力系数高达3.2;平均阻力系数在1.0×104到1.2×105的雷诺数区间内的值为1.3~2.0,并随雷诺数的增大而减小.本文提出的经验公式可准确估算高雷诺数下涡激振动发生时柔性立管的阻力系数,此经验公式考虑了流速、涡激振动主导模态以及主导频率对阻力系数的影响.      

中等Reynolds数平板绕流的动理论分析

郭永怀先生1953年给出的中等Reynolds数下、不可压缩流体有限长平板绕流的解析解是边界层理论中的经典工作.许多研究者对平板绕流阻力系数的郭永怀公式以及后续工作进行了 评估,评估的依据是Janour与Schaaf和Sherman的实验数据.本文的动理论分析和计算表 明: Schaaf和Sherman在低亚声速条件下(郭永怀先生1953年给出的中等Reynolds数下、不可压缩流体有限长平板绕流的解析解是边 界层理论中的经典工作. 许多研究者对平板绕流阻力系数的郭永怀公式以及后续工作进行了 评估, 评估的依据是Janour与Schaaf和Sherman的实验数据. 本文的动理论分析和计算表 明: Schaaf和Sherman在低亚声速条件下(郭永怀先生1953年给出的中等Reynolds数下、不可压缩流体有限长平板绕流的解析解是边 界层理论中的经典工作. 许多研究者对平板绕流阻力系数的郭永怀公式以及后续工作进行了 评估, 评估的依据是Janour与Schaaf和Sherman的实验数据. 本文的动理论分析和计算表 明: Schaaf和Sherman在低亚声速条件下(郭永怀先生1953年给出的中等Reynolds数下、不可压缩流体有限长平板绕流的解析解是边 界层理论中的经典工作. 许多研究者对平板绕流阻力系数的郭永怀公式以及后续工作进行了 评估, 评估的依据是Janour与Schaaf和Sherman的实验数据. 本文的动理论分析和计算表 明: Schaaf和Sherman在低亚声速条件下(郭永怀先生1953年给出的中等Reynolds数下、不可压缩流体有限长平板绕流的解析解是边界层理论中的经典工作.许多研究者对平板绕流阻力系数的郭水怀公式以及后续工作进行了评估,评估的依据是Janour与Schaaf和Sherman的实验数据.本文的动理论分析和计算表明:Schaaf和Sherman在低亚声速条件下(0.16相似文献   

Richtmyer-Meshkov不稳定性流体混合区发展的实验研究

使用矩形激波管，在马赫数分别为$M=1.5$和1.7的条件下实验研究了气/液界面 上(即Atwood数$A$接近1时)由Richtmyer-Meshkov不稳定性引起的流体混合现象. 得到 了气/液界面上Richtmyer-Meshkov不稳定性后期流体混合区域宽度随时间的发展呈现出线 性关系的结果，即$h \propto t$. 比较了不同马赫数和初始扰动下的发展情况，发现当 马赫数增加时，同一时间混合区域 宽度随之增加，混合区域宽度增长变快；而相比于波长差别不大的弱多模态初始扰动(无人 为干扰界面), 当界面初始扰动获得随机外界干扰时，界面混合区域具有较大的宽度以及增 长速度. 并且增加激波马赫数和初始扰动多模态性，流体混合程度更为剧烈.     

Air／SF6界面RM不稳定性诱导的物质混合诊断

不同流体伴随激波的作用,在不同密度介质的界面处往往存在激波诱导的界面不稳定性,即RM(Richtmyer-Meshkov)不稳定性,由此不稳定诱导了物质间的相互混合。文中研究了低马赫(1.23)激波作用Air/SF6界面RM不稳定性问题。Air/SF6初始界面由厚度为1～2μm的硝化纤维薄膜相隔得到,利用阴影测试法研究了Air/SF6界面在1.23马赫数激波冲击下,界面混合宽度随时间的发展过程,以及反射激波作用后混合宽度的再增长。实验结果表明混合宽度线性发展前期与经验公式吻合较好。     

弹体侵彻混凝土靶体的尺寸效应分析

由于混凝土靶体抗刚性弹侵彻实验大多基于缩比弹体展开，侵彻深度相似律是否成立显得尤为重要。在侵彻相似模型基础上，综合分析已有侵彻实验数据及经验公式，发现侵彻深度在通常情况下存在尺寸效应，且无量纲侵彻深度随弹体尺寸变大而增大。但如果模型以及原型实验中弹体与混凝土靶体（包括粗骨料）严格等比例设计，侵彻深度相似律是成立的。不变的骨料特征（粗骨料未随弹体尺寸缩放）是引起侵彻实验以及侵彻经验公式中尺寸效应的主要原因。为研究由粗骨料引起的侵彻尺寸效应，开发了混凝土二维细观有限元建模程序，细观数值实验成功地反映出了侵彻尺寸效应，考虑模拟尺寸效应后的侵彻公式能较好地预测不同尺寸侵彻实验。     

上覆岩土条件下孤石爆破预处理炸药单耗实验

盾构机在风化花岗岩地层掘进会遇到高强度孤石,为降低施工风险、避免刀盘及刀具严重损坏,常采用地面钻孔爆破对前方孤石进行预处理,但其炸药单耗一般由经验公式确定,导致爆破效果不理想。结合实际工程对爆破块度小于30 cm的特殊要求,采用爆破模型实验对炸药单耗进行研究。实验结果表明:满足特殊块度要求的陆地钻孔爆破炸药单耗是常规钻孔爆破的5.4~6.5倍;上覆地层厚度与炸药单耗呈线性递增关系;实验得出的炸药单耗是瑞典经验公式计算值的3.4~4.9倍。由此提出满足工程需要的炸药单耗修正公式,应用于实际工程并通过验证取得良好效果,为类似工程孤石爆破预处理炸药单耗的计算提供依据和参考。     

考虑惯性效应的渗流方程及相关阻力标度律

本文基于体积平均法推导得到了多孔介质中考虑惯性效应时的局部线化宏观流动方程,由此可以递推得到较大雷诺数Re 时Navier-Stokes 方程的解,从而避免了直接求解Navier-Stokes 方程所带来的计算成本高和计算稳定性差的问题．针对正方形周期排列模型的算例表明,平均速度方向与宏观压力梯度方向并不总是一致,一般情况下,随着Re 增大,二者差异也会增大．当固定平均速度方向v ? 时,压差阻力在较大的Re 范围内存在标度律,标度指数约为3．该标度指数与弱惯性区域标度指数一致,但弱惯性区域Re 范围仅为0相似文献   

高应变率下断裂韧性实验的数值模拟

采用有限元软件ANSYS/LS-DYNA程序对静态和冲击荷载作用下的含裂纹半圆弯曲(SCB)实验进行了数值模拟。根据静态实验的模拟结果,提出了适合复合型加载的Ⅰ型应力强度因子拟合公式,采用该公式计算应力强度因子的最大误差不超过10%。动态实验的模拟结果表明:对于纯Ⅰ型加载的SCB实验,动态应力强度因子随着试样半径、支座间距以及相对裂纹长度的变化呈现规律性变化;当试样半径小于60mm、相对支座间距为1.2、相对裂纹长度在0.1~0.4范围内时,惯性效应的影响较小,采用静态拟合公式计算裂尖的动态应力强度因子的误差约10%;对于复合型加载的SCB实验,当相对裂纹长度为0.2~0.4、裂纹倾角在10°~40°范围内时,采用静态拟合公式计算裂尖的动态应力强度因子的误差小于10%。     

关于湍流标度律的争鸣

研究湍流结构函数的标度律。实验或数值模拟得到的湍流结构函数的标度指数是奇异的。很多学者认为:这一实验事实否定Kolmogorov1941年(K41)提出的正常标度律,各向同性湍流惯性区的标度律是奇异的。近年来作者发表一系列文章,提出不同的观点:由于有限雷诺数效应,有限雷诺数湍流的标度指数不等于真正的惯性区标度指数,湍流结构函数的标度指数的实验数据并不否定K41正常标度律,各向同性湍流惯性区的标度律可能是正常的。惯性区奇异标度律和正常标度律对应的湍流物理本质是完全不同的,因而研究解决这个争论具有重要的意义。      

无伞末敏弹气动特性风洞实验研究

为实现末敏弹的无伞稳态扫描运动,设计了一种短圆柱轴向非对称尾翼子弹气动外形。设计了模型和实验装置,进行了大迎角低速风洞实验。获得了模型在有/无减旋翼、固定和自由旋转条件下的气动力数据,测量了模型在气动力作用下的转速。实验结果表明:设计的尾翼能够为模型提供较大的阻力,旋转时模型的阻力系数有所减小;在有减旋翼的情况下,模型能够保持较低的稳定转速;模型的静稳定性较差;模型旋转时稳定性有较大的改善。需要进一步改进尾翼,提高静稳定性,为扫描角的稳定提供保证。实验结果可为无伞末敏弹的设计和改进提供参考。     

乘波体前体/进气道一体化设计研究

为了研究乘波体几何外形参数和飞行参数对前体/进气道一体化设计的影响,采用理论分析和数值模拟相结合的方法,以马赫数Ma=6和攻角α=0为设计状态、进气道总压恢复系数和前体阻力系数为目标函数,对乘波体前体/进气道进行了优化设计,并在此基础上研究了攻角、马赫数、前缘半径、前体宽度对气动参数的影响。结果表明:该乘波体前体/进气道构型具有良好的攻角特性,总压恢复系数比基准构型提高17.79%,阻力系数比基准构型降低78.5%,符合高超声速飞行器高升力、低阻力的要求,且非常适合小攻角高超声速巡航飞行;为了得到较高升阻比的前体,在前缘半径R≤2mm的范围内进行流场反设计时,可以将设计马赫数的取值比预期低一些。     

Experimental study of the Richtmyer-Meshkov instability induced by a Mach 3 shock wave

An experimental investigation of a shock-induced interfacial instability (Richtmyer-Meshkov instability) is undertaken in an effort to study temporal evolution of interfacial perturbations in the late stages of development. The experiments are performed in a vertical shock tube with a square cross-section. A membraneless interface is prepared by retracting a sinusoidally shaped metal plate initially separating carbon dioxide from air, with both gases initially at atmospheric pressure. With carbon dioxide above the plate, the Rayleigh-Taylor instability commences as the plate is retracted and the amplitude of the initial sinusoidal perturbation imposed on the interface begins to grow. The interface is accelerated by a strong shock wave (M = 3.08) while its shape is still sinusoidal and before the Kelvin-Helmholtz instability distorts it into the well known mushroom-like structures; its initial amplitude to wavelength ratio is large enough that the interface evolution enters its nonlinear stage very shortly after shock acceleration. The pre-shock evolution of the interface due to the Rayleigh-Taylor instability and the post-shock evolution of the interface due to the Richtmyer-Meshkov instability are visualized using planar Mie scattering. The pre-shock evolution of the interface is carried out in an independent set of experiments. The initial conditions for the Richtmyer-Meshkov experiment are determined from the pre-shock Rayleigh-Taylor growth. One image of the post-shock interface is obtained per experiment and image sequences, showing the post-shock evolution of the interface, are constructed from several experiments. The growth rate of the perturbation amplitude is measured and compared with two recent analytical models of the Richtmyer-Meshkov instability.PACS: 52.35.Py, 52.35.Tc     

A parameteric study of the Busemann biplane phenomena

This paper presents a numerical investigation of the flow field around a Busemann biplane. Several parameters were investigated numerically in order to find their effect on the drag coefficient of the Busemann biplane at various flow conditions. It was found that under most flow conditions the Busemann biplane has a lower drag coefficient than that of a similar single wing.      

Development of the Richtmyer-Meshkov instability on an interface with continuous density variation

The development of the Richtmyer-Meshkov instability driven by a shock wave reflected from a rigid wall is investigated. It is shown that the perturbation amplitude growth rate depends on the nature of the shock wave refraction on the interface between gases of different density. In the case of regular-soft refraction the reduction in the growth rate of the Richtmyer-Meshkov instability associated with the continuous change in density on the interface is measured.Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 5, pp. 52–58, September–October, 1996.     

小扰动界面形态对RM不稳定性影响的数值分析

采用Navier-Stokes 方程对入射激波及其反射激波连续诱导小振幅扰动界面的Richtmyer-Meshkov 不稳定性增长过程进行了二维数值模拟,分析了单模和随机多模两种不同初始形态的界面上钉结构和泡结构在反射激波作用前后的发展特性. 研究结果发现:单模扰动的初始界面形态对反射激波前、后界面的扰动增长都有影响,反射激波前的界面形态信息可以通过钉和泡结构之间的反转传递到反射激波过后. 扰动界面上钉结构的发展速度控制了界面混合区总体的发展速度,反射激波前界面上发展成具有完整冠部形态的钉,在反射激波后会反转成复杂的泡结构,此泡结构对反射激波后钉的发展不利. 随机多模界面显示了与单模界面类似的发展规律,但随机多模界面上的复杂泡结构分布的不对称性使得其对钉结构增长的拖曳效应相对要弱,这导致了相似扰动波长下多模随机界面的扰动发展相对单模界面扰动发展要快.      

NUMERICAL STUDY ON THE EFFECTS OF SMALL-AMPLITUDE INITIAL PERTURBATIONS ON RM INSTABILITY

采用Navier-Stokes 方程对入射激波及其反射激波连续诱导小振幅扰动界面的Richtmyer-Meshkov 不稳定性增长过程进行了二维数值模拟,分析了单模和随机多模两种不同初始形态的界面上钉结构和泡结构在反射激波作用前后的发展特性. 研究结果发现：单模扰动的初始界面形态对反射激波前、后界面的扰动增长都有影响,反射激波前的界面形态信息可以通过钉和泡结构之间的反转传递到反射激波过后. 扰动界面上钉结构的发展速度控制了界面混合区总体的发展速度,反射激波前界面上发展成具有完整冠部形态的钉,在反射激波后会反转成复杂的泡结构,此泡结构对反射激波后钉的发展不利. 随机多模界面显示了与单模界面类似的发展规律,但随机多模界面上的复杂泡结构分布的不对称性使得其对钉结构增长的拖曳效应相对要弱,这导致了相似扰动波长下多模随机界面的扰动发展相对单模界面扰动发展要快.     

A drag coefficient model for Lagrangian particle dynamics relevant to high-speed flows

A blended drag coefficient model is constructed using a series of empirical relations based on Reynolds number, Mach number, and Knudsen number. When validated against experiments, the drag coefficient model produces matching values with a standard deviation error of 2.84% and a maximum error of 11.87%. The model is used in a Lagrangian code which is coupled to a hypersonic aerothermodynamic CFD code, and the particle velocity and trajectory are validated against experimental results. The comparative results agree well and show that the new blended drag coefficient model is capable of predicting the particle motion accurately over a range of Reynolds number, Mach number, and Knudsen number.     

Porous Medium Modeling of Air-Cooled Condensers

This article presents a porous media transport approach to model the performance of an air-cooled condenser. The finned tube bundles in the condenser are represented by a porous matrix, which is defined by its porosity, permeability, and the form drag coefficient. The porosity is equal to the tube bundle volumetric void fraction and the permeability is calculated by using the Karman–Cozney correlation. The drag coefficient is found to be a function of the porosity, with little sensitivity to the way this porosity is achieved, i.e., with different fin size or spacing. The functional form was established by analyzing a relatively wide range of tube bundle size and topologies. For each individual tube bundle configuration, the drag coefficient was selected by trial and error so as to make the pressure drop from the porous medium approach match the pressure drop calculated by the heat exchanger design software ASPEN B-JAC. The latter is a well-established commercial heat exchanger design program that calculates the pressure drop by using empirical formulae based on the tube bundle properties. A close correlation is found between the form drag coefficient and the porosity with the drag coefficient decreasing with increasing porosity. A second order polynomial is found to be adequate to represent this relationship. Heat transfer and second law (of thermodynamics) performance of the system has also been investigated. The volume-averaged thermal energy equation is able to accurately predict the hot spots. It has also been observed that the average dimensionless wall temperature is a parabolic function of the form drag coefficient. The results are found to be in good agreement with those available in the open literature.     

A unified view of energetic efficiency in active drag reduction,thrust generation and self-propulsion through a loss coefficient with some applications

An analysis of the energy budget for the general case of a body translating in a stationary fluid under the action of an external force is used to define a power loss coefficient. This universal definition of power loss coefficient gives a measure of the energy lost in the wake of the translating body and, in general, is applicable to a variety of flow configurations including active drag reduction, self-propulsion and thrust generation. The utility of the power loss coefficient is demonstrated on a model bluff body flow problem concerning a two-dimensional elliptical cylinder in a uniform cross-flow. The upper and lower boundaries of the elliptic cylinder undergo continuous motion due to a prescribed reflectionally symmetric constant tangential surface velocity. It is shown that a decrease in drag resulting from an increase in the strength of tangential surface velocity leads to an initial reduction and eventual rise in the power loss coefficient. A maximum in energetic efficiency is attained for a drag reducing tangential surface velocity which minimizes the power loss coefficient. The effect of the tangential surface velocity on drag reduction and self-propulsion of both bluff and streamlined bodies is explored through a variation in the thickness ratio (ratio of the minor and major axes) of the elliptical cylinders.