磁场对多孔介质方腔内空气热磁对流的影响

数值分析了微重力下圆形载流线圈倾斜时多孔介质方腔内空气热磁对流. 磁场计算采用毕奥--萨伐定律求解; 动量方程与能量方程分别采用达西模型与局部热非平衡模型求解. 计算结果表明随着磁场力数\gamma 数和Da数的增加, 方腔内对流变得越来越强. 线圈倾斜角x_{\rm euler}从0^\circ到90^\circ变化时, 对流结果关于x_{\rm euler}=45^\circ呈现对称分布. Nu_{\rm m}数随线圈倾斜角的改变而变化且每个工况下局部最大Nu_{\rm m}数出现在x_{\rm euler}=45^\circ. 局部最小Nu_{\rm m}数出现在x_{\rm euler}=0^\circ, 90^\circ.      

细长椭球体在水中自由下落的运动特性实验研究

固体颗粒在液体中的运动现象在日常生活和工业应用领域广泛存在,其中因蕴含着丰富的流体力学现象而受到学者们的广泛关注.本文通过实验研究了细长椭球体在水中受浮力影响的下落特性.实验中采用带有两台相互垂直的高速摄像机和光源组成的运动跟踪平台并结合荧光染色技术对细长椭球体下落过程中的运动轨迹和尾涡结构进行研究.文中选用的细长椭球体与环境流体的密度比为1.2,其长短轴比范围为2～10,相应的阿基米德数范围为400～1400,对应实现的终态雷诺数范围为120～1350.实验过程中我们观察到细长椭球体在水中下落过程中产生的5种典型路径,分别为:小振幅不规则运动、小振幅高频振荡运动、大振幅低频振荡运动、高度非线性运动以及直线运动,并得到了对应的速度振荡以及倾斜角的演化规律.进一步地,分析了细长椭球体运动过程中受到的阻力系数与雷诺数之间的关系.随后采用荧光可视化技术清晰获得了颗粒下落过程中的尾涡结构特性,并结合颗粒的运动状态详细分析了涡脱落过程对颗粒运动状态转捩的影响.最后,通过对比前人关于圆柱体下落的运动特性的相关结果,获得了细长椭球体和细长圆柱体运动特性之间的异同点以及其潜在的物理机理.     

砂-膨润土混合屏障材料渗透性影响因素研究

建立了一个新的结构-尾流振子耦合模型. 流场近尾迹动力学特征被模化为非线性阻尼 振子,采用van der Pol方程描述. 以控制体中结构与近尾迹流体间受力互为反作 用关系来实现流固耦合. 采用该模型进行了二维结构涡激振动计算,得到了合理的 振幅随来流流速的变化规律和共振幅值,并正确地预计了共振振幅值$A_{\max}^\ast$ 随着质量阻尼参数$\left( {m^\ast + C_A } \right)\zeta $的变化规律,给出了预测$A_{\max }^\ast $值的拟合公式. 采用该模型计算了三维柔性结构在均匀来流和简谐波形来流作用下的VIV 响应. 结构在均匀来流作用下振动呈现由驻波向行波的变化过程, 并最后稳定为行波振动形态. 在简谐波形来流作用下,结构呈现混合振动形态,幅值随时间呈周期变化.     

圆柱涡激振动的结构-尾流振子耦合模型研究

建立了一个新的结构-尾流振子耦合模型. 流场近尾迹动力学特征被模化为非线性阻尼振子,采用van der Pol方程描述. 以控制体中结构与近尾迹流体间受力互为反作用关系来实现流固耦合. 采用该模型进行了二维结构涡激振动计算,得到了合理的振幅随来流流速的变化规律和共振幅值,并正确地预计了共振振幅值$A_{\max}^\ast$随着质量阻尼参数$\left( {m^\ast + C_A } \right)\zeta$的变化规律,给出了预测$A_{\max }^\ast$值的拟合公式. 采用该模型计算了三维柔性结构在均匀来流和简谐波形来流作用下的VIV响应. 结构在均匀来流作用下振动呈现由驻波向行波的变化过程, 并最后稳定为行波振动形态.在简谐波形来流作用下,结构呈现混合振动形态,幅值随时间呈周期变化.      

细长椭球体在水中自由下落的运动特性实验研究

固体颗粒在液体中的运动现象在日常生活和工业应用领域广泛存在, 其中因蕴含着丰富的流体力学现象而受到学者们的广泛关注. 本文通过实验研究了细长椭球体在水中受浮力影响的下落特性. 实验中采用带有两台相互垂直的高速摄像机和光源组成的运动跟踪平台并结合荧光染色技术对细长椭球体下落过程中的运动轨迹和尾涡结构进行研究. 文中选用的细长椭球体与环境流体的密度比为1.2, 其长短轴比范围为2$\sim $10, 相应的阿基米德数范围为400$\sim$1400, 对应实现的终态雷诺数范围为120$\sim$1350. 实验过程中我们观察到细长椭球体在水中下落过程中产生的5种典型路径, 分别为: 小振幅不规则运动、小振幅高频振荡运动、大振幅低频振荡运动、高度非线性运动以及直线运动, 并得到了对应的速度振荡以及倾斜角的演化规律. 进一步地, 分析了细长椭球体运动过程中受到的阻力系数与雷诺数之间的关系. 随后采用荧光可视化技术清晰获得了颗粒下落过程中的尾涡结构特性, 并结合颗粒的运动状态详细分析了涡脱落过程对颗粒运动状态转捩的影响. 最后, 通过对比前人关于圆柱体下落的运动特性的相关结果, 获得了细长椭球体和细长圆柱体运动特性之间的异同点以及其潜在的物理机理.      

基于椭球曲面拟合的三维磁罗盘误差补偿算法

针对现有三维磁罗盘误差补偿速度慢、需专用转台、野外使用不便的问题,提出了一种建立总磁场的不等边椭球曲面数学模型进行误差补偿的方法.首先根据Poisson方程描述的总磁场测量模型,得出一般椭球曲面模型,然后,采用牛顿迭代法对椭球曲面方程线性化,分别使用最小二乘法和总体最小二乘法计算线性化方程,推导椭球曲面拟合的数学公式,计算椭球曲面参数,最后,利用最小二乘法提出并推导参数初值的计算公式,给出了选取数据点的方法.实验结果表明,仅手持磁罗盘在各象限旋转即可实现误差补偿,航向测量精度可达0.8°,补偿精度与传统的转台补偿基本相同,而补偿方式更为灵活,适用于无转台设备的野外环境测量使用.     

超空泡航行体前部线形对空泡生成临界通气量影响的实验研究

基于航行体超空泡理论和格兰威尔线形设计方法,设计了三种具有不同前部线形的航行体模型.并针对所设计的三种模型和具有锥形前部外型的航行体模型在水洞中进行了生成超空泡临界通气量的实验研究.结果表明:航行体前部线型对超空泡临界通气量的影响程度与航行体的轴向斜率分布有关,模型表面斜率轴向分布曲线越平坦或变化率越小,越有利于减少超空泡生成所需的临界通气量.实验数据显示设计的三种格兰威尔前部线型航行体与锥形前部外型航行体相比,生成超空泡临界通气量都有明显减小,最大减小量达24%,为降低超高速航行体超空泡生成所需的临界通气量提供了一条技术途径和研究方法.     

砂土串囊式充气锚杆承载力的研究

在砂土地层中,串囊式充气锚杆的研究还比较少,其承载特性及受力机理尚不明确。本文基于莫尔-库仑模型和Vesic圆孔扩张理论法,分别对圆柱体、球体、组合体、椭球体假设下的串囊式充气锚杆的扩大段进行计算分析。并将计算结果与试验得到的实测值进行对比。结果表明:四种形状假设中椭球体的形状假设理论值与实测值的误差最小,仅为8.35%。通过拟合试验数据,并引入与端阻力和侧摩阻力有关的两个系数对承载力公式进行修正,得到了抗拔承载力的经验公式。     

基于深度神经网络的横流转捩预测

基于深度神经网络DNN构建了从层流流场无量纲速度梯度、流向涡强度等物理量到横流转捩模态下间歇因子间的映射关系,获得一种新的数据驱动转捩模型.通过将数据驱动转捩模型与SST k-ω湍流模型耦合,有效简化了转捩模型输运方程求解,实现高效、准确的亚音速三维边界层横流转捩流场计算. DNN训练数据来自变雷诺数的NLF(2)-0415无限展长后掠翼计算结果,并以两种工况进行测试,数据驱动转捩模型预测精度与γ-Reθ转捩模型近似.将数据驱动转捩模型用于其他典型横流转捩算例的计算,以验证其泛化能力.对于变后掠角的NLF(2)-0415后掠翼,数据驱动转捩模型与γ-Reθt-CF模型预测的转捩位置几乎一致,并且能够预测出后掠角从45°增长到65°的过程中,转捩位置先向前再向后移动的现象;对于标准椭球体,使用低分辨率网格进行计算,数据驱动转捩模型依然能够实现转捩位置预测,对椭球体表面Cf的计算结果与多个平台的横流转捩模型、实验结果基本一致.研究表明,以横流转捩相关物理量作为输入对DNN进行训练,并将获得的数据驱动转捩模型与SST k-ω湍流模型耦合,可以实现对横流转捩的有效预测,且具有较强的泛化能力.数...     

一种电力架空线垂度的简算方法

以输配电工程中带有绝缘子的电力架空线为研究对象, 简化绝缘子的力学模型, 采用纵横弯 曲理论推导了钢芯铝铰线的挠度方程, 计算了在自重和温度载荷共同作用下的架空线垂度和 水平张力. 计算结果表明, 文中算法简单有效, 架空线跨距和温度等是架空线垂度的主要 影响因素. 架空线垂度对温度的平均变率为2.55mm/$^{\circ}$C; 当跨距减 小1%时, 垂度增加近60%.     

The influence of gas phase velocity fluctuations on primary atomization and droplet deformation

The effects of grid-generated velocity fluctuations on the primary atomization and subsequent droplet deformation of a range of laminar liquid jets are examined using microscopic high-speed backlit imaging of the break-up zone and laser Doppler anemometry of the gas phase separately. This is done for fixed gas mean flow conditions in a miniature wind tunnel experiment utilizing a selection of fuels, turbulence-generating grids and two syringe sizes. The constant mean flow allows for an isolated study of velocity fluctuation effects on primary atomization in a close approximation to homogeneous decaying turbulence. The qualitative morphology of the primary break-up region is examined over a range of turbulence intensities, and spectral analysis is performed in order to ascertain the break-up frequency which, for a case of no grid, compares well with the existing literature. The addition of velocity fluctuations tends to randomize the break-up process. Slightly downstream of the break-up region, image processing is conducted in order to extract a number of metrics, which do not depend on droplet sphericity, and these include droplet aspect ratio and orientation, the latter quantity being somewhat unconventional in spray characterization. A turbulent Weber number $We^{\prime}$ which takes into account gas phase fluctuations is utilized to characterize the resulting droplet shapes, in addition to a mean Weber number <We _d>. Above a $We^{\prime}>0.05$ a clear positive relationship exists between the mean aspect ratio of droplets and the turbulent Weber number where $We^{\prime}$ is varied by altering all relevant variables including the velocity root mean square, the initial droplet diameter, the surface tension and the density.     

旋转弹丸入水侵彻规律

为了建立高速弹体入水弹道的模型,利用数字式高速录像机实验研究了球形与普通手枪两种弹丸在三个入水角、六种发射速度下斜入水的水中弹道轨迹与空泡。实验结果表明,弹丸形状对入水弹道的稳定性有重要影响。球形弹丸在斜入水时弹道稳定性较好,而普通制式弹丸的弹道不稳定。在一定的入水角与速度范围内,球形弹丸入水初期弹道的空泡特性、弹道轨迹以及速度衰减规律具有一定的相似性。而对于水中速度随时间的衰减规律,则两种弹丸都具有一定的相似性,且表现出极强的速度衰减特性。给出了弹丸水中速度衰减规律的数学预报模型,并与实验结果进行了比较,理论结果与实验结果吻合较好。     

Visualization of high-speed gas jets and their airblast sprays of cross-injected liquid

A planar and instantaneous visualization study of high-speed gas jets and their airblast sprays was performed to qualitatively examine the different atomization performances of different gas nozzles. For the visualization of high-speed gas jets (with no liquid injected), Nd:YAG pulsed laser sheets imaged the clustered vapor molecules in the Rayleigh range (d?λ), condensed from the natural humidity during the isentropic gas expansion through a nozzle. This method visualized both underexpanded sonic gas jets from a converging nozzle (SN-Type) and overexpanded supersonic gas jets from a converging-diverging nozzle (CD-Type). When liquid is cross-injected, the same laser sheet images the spray droplets of relatively large sizes (d?λ). The present visualization results show that the SN-Type nozzle develops a wider spray than the CD-Type nozzle, quite probably because the SN-Type nozzle has a wider gas jet (in the absence of liquid) than the CD-Type. Also, the wider spray of the SN-Type nozzle lowers the probability of droplet coalescence and generates finer sprays compared to the CD-Type nozzle. These visualization results qualitatively agree with the previous quantitative finding of the different atomization characteristics of the two types of nozzles (Park et al. 1996).     

Adaptive collision meshing and satellite droplet formation in spray simulations

Improved numerical methods and physical models have been applied to droplet collision modeling. Numerically, an adaptive collision mesh method is developed to calculate collision rate. This method produces a collision mesh that is independent of the gas phase mesh and adaptively refined according to local parcel number density. An existing model describing the satellite droplet formation during the collision process is improved to reflect the experimental findings that the satellite droplets are much smaller than the parent droplets. The adaptive collision mesh and the improved satellite model have been used to simulate three impinging spray experiments. The model was able to qualitatively predict the occurrence of small satellite drops and bi-modal post-collision drop size distributions. The effect of the collision mesh and the satellite droplet model on a high-speed non-evaporating diesel spray is also assessed.     

Experimental estimation of evaporation rates of water droplets in high-temperature gases

Evaporation rates of water droplets in high-temperature gases were experimentally determined using high-speed video recording cameras and low-inertia thermocouples (for heated air flow as an example). The experiments were carried out for droplets of initial size (radius) 1–3 mm at an air temperature of 500–1000 K. Dependences of the evaporation rate of water droplets on time and gas temperature were obtained for various initial droplet sizes.     

Numerical modeling and experimental measurements of a high speed solid-cone water spray for use in fire suppression applications

Experimental measurements and numerical simulations of a high-speed water spray are presented. The numerical model is based on a stochastic separated flow technique that includes submodels for droplet dynamics, heat and mass transfer, and droplet–droplet collisions. Because the spray characteristics near the nozzle are difficult to ascertain, a new method for initialization of particle diameter size is developed that assumes a Rosin–Rammler distribution for droplet size, which correctly reproduces experimentally measured Sauter and arithmetic mean diameters. By relating the particle initialization to lower moments of the droplet statistics, it is possible to take advantage of measurements without substantial penalties associated with the greater experimental uncertainty of individual droplet measurements. Overall, very good agreement is observed in the comparisons of experimental measurements to computational predictions for the streamwise development of mean drop size and velocity. In addition, the importance of modeling droplet–droplet collisions is highlighted with comparison of selected droplet–droplet collision models.     

Multi-frame visualization for detonation wave diffraction

When a detonation wave emerges from a tube into unconfined space filled with a gas mixture, detonation wave diffraction occurs due to abrupt changes in the cross-sectional area. In the present study, we focused on the local explosion in reinitiation and propagation of a transverse detonation wave by performing comprehensive and direct observation with high time resolution visualization in a two-dimensional rectangular channel. Using the visualization methods of shadowgraph and multi-frame, short-time, open-shutter photography, we determined where the wall reflection point is generated, and also determined where the bright point is originated by the local explosion, and investigated the effects of the deviation angle and initial pressure of the gas mixture. We found that the reinitiation of detonation had two modes that were determined by the deviation angle of the channel. If the deviation angle was less than or equal to 30\(^{\circ }\), the local explosion of reinitiation might occur in the vicinity of the channel wall, and if the deviation angle was greater than or equal to 60\(^{\circ }\), the local explosion might originate on the upper side of the tube exit. With a deviation angle greater than 60\(^{\circ }\), the position of the wall reflection point depended on the cell width, so the radial distance of the wall reflection point from the apex of the tube exit was about 12 times the cell width. Similarly, the bright point (local explosion point) was located a distance of about 11 times the cell width from the apex of the tube exit, with a circumferential angle of 48\(^{\circ }\).     

Spark-generated bubble collapse near or inside a circular aperture and the ensuing vortex ring and droplet formation

This paper aims to study the oscillation of a sparkgenerated submerged bubble located near or inside a circular aperture made in a flat plate using high-speed visualization technique. In the case of a bubble oscillating near an aperture the initial free surface of the water was set at the bottom surface of the plate. The effects of aperture size and bubblefree surface distance on the bubble behavior as well as on the ensuing droplet dynamics are investigated. It was found that the direction of the bubble reentrant jet was towards the aperture or away from it respectively when the normalized aperture size was smaller or greater than a certain critical value. In addition, a toroidal vortex ring was observed to form, which rotated inwards as it moved away from the aperture. It was also found that if the bubble was incepted at a distance sufficiently away from a supercritical size aperture a single droplet could be produced. In the case of a bub- ble initiated in the middle of a circular aperture submerged just beneath the water free surface, the bubble was found to take the shape of an ellipsoid during its expansion. Then a reentrant jet was initiated and pierced the bubble from its top side.     

Microbubble formations in MEMS-fabricated rectangular channels: A high-speed observation

A microfluidic device that consists of MEMS-fabricated rectangular channels is developed for stable and sequential production of monodispersed microbubbles. The central inlet for gas phase is located between two inlets for liquid phases, where the device works as a two-fluid atomizer. The interfacial flow mechanism of microbubble formation at the junction of the inlets in the device is investigated using a high-speed visualization technique and digital image processing. The periodic formation process is successfully realized by the consideration of the wettability between the microbubble and the channel wall. The produced microbubbles are relatively uniform in size, and the size is controlled from 113 to 153 μm by changing the flow rates of the liquid and gas phases. Furthermore, a simple theoretical model to predict the equivalent diameter of microbubbles is developed by considering the mass balance of the gas phase in the formation process, where the experimental and theoretical results are in reasonable agreement.     

Extension of k–ε models for a turbulent falling film

An extension applicable to various k–ε models is proposed to account for the damping of turbulence due to the surface tension. It consists of a sink in the equation for $ {{\overline{D} k} \mathord{\left/ {\vphantom {{\overline{D} k} {\overline{D} t}}} \right. \kern-0em} {\overline{D} t}} $ and a source in the equation for $ {{\overline{D} \varepsilon } \mathord{\left/ {\vphantom {{\overline{D} \varepsilon } {\overline{D} t}}} \right. \kern-0em} {\overline{D} t}} $ both derived from dimensional analysis. First numerical experiments are undertaken with a commercial CFD software. Reasons are given why the tuning of the model should be performed with thermal experiments. The proposed model is practical, since it only needs the programming of source terms. On account of it’s mathematical structure it needs a very small closure coefficient. The intention of this article is to stimulate the numerical turbulence research in a way that is applicable to the engineering practice. A comparison to experimental data is not done here.