波浪作用下直立结构物附近强湍动掺气流体运动的数值模拟

波浪破碎卷入气体易对建筑物受力产生压力振荡,了解波浪作用下建筑物附近掺气水流的运动特性是精确计算建筑物受力的前提.基于OpenFOAM开源程序包和修正速度入口造波方法建立三维数值波浪水槽,模型采用S-A IDDES湍流模型进行湍流封闭,并采用修正的VOF方法捕捉自由液面,数值模拟了规则波在1:10的光滑斜坡上与直立结构物的相互作用过程,重点分析了结构物附近的水动力和掺气水流运动特性.结果表明,建立的数值模型能精确地捕捉波浪作用下直立结构物附近的自由液面的变化以及气泡输运过程,较好地描述气体卷入所形成的气腔形态以及多气腔之间的融合、分裂等过程;波浪与直立结构物相互作用产生强湍动掺气水流,其运动过程十分复杂;掺气流体输运过程中水气界面周围一直伴随着涡的存在,其中,气泡的分裂与周围正负涡量剪切作用密切相关,且其输运轨迹主要受周围流场的影响;研究揭示了结构物附近湍动能与掺气特性的关系,发现波浪作用下直立结构物附近湍动能的分布与掺气水流特征参数(气泡数量、空隙率)整体呈现一定的线性关系.     

水翼端部间隙泄漏流的压降及黏性损失机理

叶顶泄漏流产生的局部压降及黏性损失是导致轴流式水力机械效率下降和轮缘间隙空化的主要原因.为探明间隙泄漏流的黏性损失特性和低压形成机制,以NACA0009水翼为对象,采用超大涡模拟方法(VLES)对翼端间隙流动进行数值模拟,基于平均流动动能转换与输运分析,提出了间隙区黏性损失定量计算模型,研究了翼端间隙区湍动能生成、黏性损失和压降的产生机理及主要影响因素.结果表明,间隙区存在间隙分离涡(TSV)、间隙泄漏涡(TLV)和诱导涡(IV)等流动结构;湍动能生成是导致TSV内压降的主导因素,TLV内压降则主要受湍动能生成和平均动能的对流和扩散效应影响;湍动能耗散导致的翼端区域黏性损失占间隙区黏性损失总量的91.2%.间隙区不同流动结构对湍动能生成的影响存在明显差异,水翼吸力面的强剪切效应主要生成湍动能的■分量,而TLV, TSV和IV等间隙涡结构则主要生成湍动能的■和■分量;湍动能产生机制分析表明,湍动能生成项分量P_vw是TLV和TSV中湍动能生成的主导因素,减小TSV和TLV内的速度梯度■,可有效降低湍动能生成,进而减少翼端区域因湍流耗散导致的黏性损失.研究结果可为间隙流...     

考虑当地涡旋运动修正的湍流模型在非定常空化湍流流场计算中的应用

为了更精确地模拟当地涡旋运动对复杂非定常空化湍流流动的影响,考虑当地涡旋运动特性对标准k-ε湍流模型进行了修正;并采用该修正模型对二维Clark-Y水翼周围非定常空化湍流场进行了计算;通过与实验结果进行对比,分析并验证了该修正模型计算结果的有效性。研究表明:该修正湍流模型根据当地涡旋效应对湍动能输运方程中的生成项进行修正,可以有效解决标准k-ε模型对湍动能的过度预测问题,得到的空化核心区域时均湍动能和湍流黏性的预测结果降低了约30%,与实际值更为接近;标准k-ε模型计算所得的时均空泡长度不大于0.6倍弦长,而修正模型计算所得的时均空穴可以覆盖水翼吸力面;修正模型可以准确预测水翼吸力面表面空泡的非定常脱落过程,包括附着空穴的拉伸和断裂,以及脱落空泡的逐渐消失。     

高超声速激波湍流边界层干扰直接数值模拟研究

高超声速激波与湍流边界层干扰会导致飞行器表面出现局部热流峰值,严重影响飞行器气动性能和飞行安全. 针对高马赫数激波干扰问题,以往数值研究多采用雷诺平均方法,而在直接数值模拟方面的相关工作较为少见. 开展高超声速激波与湍流边界层干扰的直接数值模拟研究,有助于进一步提升对其复杂流动机理认识和理解,同时也将为现有湍流模型和亚格子应力模型的改进提供理论依据. 采用直接数值模拟方法对来流马赫数6.0,34°压缩拐角内激波与湍流边界层的干扰问题进行了研究. 基于雷诺应力各向异性张量,分析了高超声速湍流边界层在压缩拐角内的演化特性. 通过对湍动能输运方程的逐项分析,系统地研究了可压缩效应对湍动能及其输运的影响机制. 采用动态模态分解方法,探讨了干扰流场的非定常运动历程. 研究结果表明,随着湍流边界层往下游发展,近壁湍流的雷诺应力状态由两组元轴对称状态逐渐演化为两组元状态,外层区域则由轴对称膨胀趋近于各向同性. 干扰流场内存在强内在压缩性效应(声效应),其对湍动能输运的影响主要体现在压力--膨胀项,而对膨胀--耗散项影响较小. 高超声速下压缩拐角内的非定常运动仍存在以分离泡膨胀/收缩为特征的低频振荡特性,其物理机制与分离泡剪切层密切相关.      

FST方法分析二维Rayleigh-Bénard湍流热对流系统的能量输运

本文应用空间滤波方法：FST(Filter-space technique)方法，研究二维Rayleigh-Bénard(RB)湍流热对流系统中湍动能、热能和拟涡能的能量输运．研究中Rayleigh数(Ra)选取为1x10^8、1x10^9和1x10^10，Prandtl数(Pr)固定为4.38．我们展示了的结果表明，在二维RB系统中，三个Ra数下全场的平均湍动能和平均拟涡能在不同滤波尺度下的能量输运与Kraichnan在1967年预测的二维湍流中的级串理论有所偏差，而中心区域的能量都是向小尺度输运的．结果还揭示了瞬时能量输运的一些局部特性，包括它们在小尺度上不对称的分布．     

高超声速激波湍流边界层干扰直接数值模拟研究

高超声速激波与湍流边界层干扰会导致飞行器表面出现局部热流峰值,严重影响飞行器气动性能和飞行安全.针对高马赫数激波干扰问题,以往数值研究多采用雷诺平均方法,而在直接数值模拟方面的相关工作较为少见.开展高超声速激波与湍流边界层干扰的直接数值模拟研究,有助于进一步提升对其复杂流动机理认识和理解,同时也将为现有湍流模型和亚格子应力模型的改进提供理论依据.采用直接数值模拟方法对来流马赫数6.0,34?压缩拐角内激波与湍流边界层的干扰问题进行了研究.基于雷诺应力各向异性张量,分析了高超声速湍流边界层在压缩拐角内的演化特性.通过对湍动能输运方程的逐项分析,系统地研究了可压缩效应对湍动能及其输运的影响机制.采用动态模态分解方法,探讨了干扰流场的非定常运动历程.研究结果表明,随着湍流边界层往下游发展,近壁湍流的雷诺应力状态由两组元轴对称状态逐渐演化为两组元状态,外层区域则由轴对称膨胀趋近于各向同性.干扰流场内存在强内在压缩性效应(声效应),其对湍动能输运的影响主要体现在压力-膨胀项,而对膨胀-耗散项影响较小.高超声速下压缩拐角内的非定常运动仍存在以分离泡膨胀/收缩为特征的低频振荡特性,其物理机制与分离泡剪切层密切相关.     

STUDY OF NON-SUBMERGED GROIN TURBULENCE FLOW IN A SHALLOW OPEN CHANNEL BYLES

采用动态亚格子模式和浸没边界法,对宽浅槽道中的丁坝群绕流的水动力学特性进行了三维大涡模拟研究. 利用丁坝绕流,试验中采用粒子图像测速仪(particle image velocimetry, PIV)测量的试验中自由水面处的时间平均流速和湍动强度数据对模型进行率定,结果表明计算结果与试验数据吻合良好. 丁坝长度与丁坝之间距离的比值L/D对丁坝周围的水流流动形式、湍流强度、涡量分布有显著影响. 在L保持不变并且L/D较大时,丁坝之间的距离D较小,这限制了混合层的发展,因此混合层中的湍动强度和涡量都较小;同时丁坝之间的回流区的流线形式也发生明显变化. 此外,还给出了涡体在丁坝坝头附近产生,发展并向下游输运的动态过程.     

一种基于湍动能方程的转捩判定方法简

转捩现象是阻碍阻力高精度求解的主要问题之一. Menter 和Langtry 所提出的γ-θ转捩模型通过引入涡量雷诺数和间歇因子输运方程来驱动转捩,但是其中很多经验公式的理论立足点有待商榷. 驱使层流转变到湍流依赖的仍然是平均速度的一阶和二阶相关量,它们组合构成了湍动能方程的耗散尺度. 在湍动能方程中做合适的耗散平衡后,仅仅依靠湍动能方程可以有效地捕捉转捩现象. 采用自然转捩和旁路转捩测试算例进行了验证,结果证明该方法与试验值匹配较好,具有一定的工程实用价值.     

一种基于湍动能方程的转捩判定方法

转捩现象是阻碍阻力高精度求解的主要问题之一. Menter 和Langtry 所提出的γ-θ转捩模型通过引入涡量雷诺数和间歇因子输运方程来驱动转捩,但是其中很多经验公式的理论立足点有待商榷. 驱使层流转变到湍流依赖的仍然是平均速度的一阶和二阶相关量,它们组合构成了湍动能方程的耗散尺度. 在湍动能方程中做合适的耗散平衡后,仅仅依靠湍动能方程可以有效地捕捉转捩现象. 采用自然转捩和旁路转捩测试算例进行了验证,结果证明该方法与试验值匹配较好,具有一定的工程实用价值.      

槽道湍流近壁区多尺度输运特性研究

利用槽道湍流直接数值模拟的数据库和离散正交子波，对近壁湍流的多尺度输运特性进 行了研究. 通过在流向和展向分别进行子波多尺度分解，得到了近壁区湍动能在流向和展向 多尺度传输的不同性质，发现流向传输以能量的反传为主，而在展向能量存在明显的正传， 并且当过滤尺度较大时以正传为主. 近壁湍流能量传输的各向异性为进一步构造各向异 性大涡模拟亚格子模式提供了必要的参考.     

Specific features of the modeling of boiling-fluid flows

Experiments show that in low-and high-velocity flows the boiling process is fundamentally different: in the former, the fluid boils on the walls, and in the latter in the volume. In high-velocity flows, the boiling intensity is orders of magnitude greater. In modeling fast and slow flows, the number of bubbles, which is a free parameter of the model and must be specified, differs by orders of magnitude. When high-speed flows of different kinds are modeled (vessel depressurization, nozzle flows) the number of bubbles specified also differs by orders of magnitude. In this study, we formulate the hypothesis that in both kinds of flows the process of boiling starts similarly, namely, on the walls. However, in high-speed flows the number of bubbles increases by orders of magnitude due to bubble fragmentation. As a result of intense fragmentation, the system “forgets” the initial number of bubbles and the process becomes volume boiling. This approach makes it possible to construct a universal model of boiling. To test this hypothesis, we constructed a mathematical model which takes into account the possibility of bubble fragmentation due to the instability developing under the action of centrifugal accelerations of the bubble surface. This model was used to calculate the process of depressurization of a high-pressure vessel. The calculations demonstrated that, for any initial number of bubbles, 1 ms after depressurization the bubble number attains the same level. Bubble fragmentation takes place in “self-sustained detonation waves”. The stationary structure of detonation waves in a boiling fluid is investigated. A scheme of the wave structure according to which the wave consists of a shock wave and a relaxation zone is proposed. Calculations of a boiling-fluid flow through a Laval nozzle reveal the periodic appearance of detonation waves. Accordingly, nozzle flows should be accompanied by significant oscillations of the parameters.     

An experimental study on turbulent coherent structures near a sheared air-water interface

The turbulence structures near a sheared air-water interface were experimentally investigated with the hydrogen bubble visualization technique. Surface shear was imposed by an airflow over the water flow which was kept free from surface waves. Results show that the wind shear has the main influence on coherent structures under air-water interfaces. Low- and high- speed streaks form in the region close to the interface as a result of the imposed shear stress. When a certain airflow velocity is reached, “turbulent spots” appear randomly at low-speed streaks with some characteristics of hairpin vortices. At even higher shear rates, the flow near the interface is dominated primarily by intermittent bursting events. The coherent structures observed near sheared air-water interfaces show qualitative similarities with those occurring in near-wall turbulence. However, a few distinctive phenomena were also observed, including the fluctuating thickness of the instantaneous boundary layer and vertical vortices in bursting processes, which appear to be associated with the characteristics of air-water interfaces. The project supported by the National Natural Science Foundation of China (Grant No.19672070)     

Air transport in chute flows

The air bubble rise velocity in still water depends mainly on the bubble size and is basically influenced by buoyancy, viscosity and surface tension. In high-speed flows the number of forces acting on air bubbles increases with turbulence, non-hydrostatic pressure gradient, shear forces, bubble clouds and free-surface entrainment. Air bubbles in these flows are used for cavitation protection of hydraulic structures such as chutes, spillways and bottom outlets. Here, air is normally added by means of aerators upstream of regions where the cavitation number falls below a critical value mainly to reduce the sonic velocity of the fluid and cushion the cavitation bubble collapse process. The distance between successive aerators depends basically on the bubble rise velocity. Until today, the bubble rise velocity in high-speed flows was not thoroughly investigated because of limited laboratory instrumentation. The present project focused on the streamwise development of air concentrations in high-speed flows along a 14 m long model chute. The bubble rise velocity was indirectly derived from the air detrainment gradient of the air concentration contour lines downstream of an aeration device. It accounts for the main hydraulic parameters chute slope, Froude number and air concentration. It is demonstrated that the bubble rise velocity in high-speed flow and stagnant water differ significantly due to fracturing processes, turbulence, and the ambient air concentration.     

Full-scale two-phase flow measurements on Athena research vessel

Measurements of gas volume fraction, bubble velocity, chord length and bubble size distributions were performed on the research vessel Athena II operating in Saint Andrew Bay in the gulf coast near Panama City, FL. Double tipped sapphire optical local phase-detection probes were used to acquire indicator functions downstream of the breaking bow wave, behind the masker and at the stern. These indicator functions were also taken at different depths, operating speeds and headings respect to the waves. The data processing includes the computation of velocity of individual bubbles and chord lengths, resulting in chord length distributions. These chord length distributions are used to obtain bubble size distributions using a novel procedure described in detail herein. Uncertainty analysis is performed for gas volume fraction, average bubble velocity and chord length. The results indicate that air entrainment increases with ship speed and sailing against the waves at all positions. The bow wave exhibits unsteady breaking that creates bubble clouds, which were characterized and identified by signal processing. At the stern a very strong dependence of bubble size with depth was found, with evidence that small bubbles (smaller than 500 μm) are transported through the bottom of the hull and reach the transom. The roller present at the transom, the associated strong unsteadiness and bubble entrainment are well captured, as indicated by the stern results, showing the frothy nature of the upper layer.     

Air entrainment in two-dimensional turbulent shear flows with partially developed inflow conditions

In plunging jet flows and at hydraulic jumps, large quantities of air are entrained at the intersection of the impinging flow and the receiving body of water. The air bubbles are entrained into a turbulent shear layer and strong interactions take place between the air bubble advection/diffusion process and the momentum shear region. New air-water flow experiments were conducted with two free shear layer flows: a vertical supported jet and a horizontal hydraulic jump. The inflows were partially developed boundary layers, characterized by the presence of a velocity potential core next to the entrapment point. In both cases, the distributions of air concentration exhibit a Gaussian distribution profile with an exponential longitudinal decay of the maximum air content. Interestingly, the location of the maximum air content and the half-value band width are identical for both flow situations, i.e. independent of buoyancy effects.     

Flow structures and dynamics in the wakes of sliding bubbles

An experimental investigation is reported for the flow structures in the wake of an air bubble sliding under an inclined surface in quiescent water. Time-resolved particle image velocimetry (PIV) is used to study the wakes of sliding bubbles for a range of measurement planes, bubble diameters and surface inclination angles. Additionally, key aspects of the bubble’s motion are measured simultaneously using a novel method that accounts for the motion of the bubble’s interface. Thus, vortex shedding may be linked to changes in the bubble shape and path.Analysis of the measured velocity and vorticity fields reveals a wake structure consisting of a near wake that moves in close proximity to the bubble, shedding vorticity at the inversion points of the bubble path. Downstream of the bubble in the far wake, these structures evolve into asymmetrical, oppositely-oriented hairpin vortices that are generated in the near wake. These hairpin vortices bear similarities to those observed behind freely rising bubbles and near-wall bluff bodies and are found to cause significant motion of the bulk fluid. This bulk fluid motion has the potential to offer significant convective cooling of adjacent heated surfaces, such as submerged electronics components.     

气泡及其破碎兴波对浮动冲击平台影响探究

基于LS-DYNA软件中的ALE算法，对近水面水下爆炸气泡脉动过程进行数值模拟并与实验结果进行对比，验证了近水面近壁面混合边界有限元模型和参数设置的正确性。设置不同爆炸工况，对气泡及其破碎兴波对浮动冲击平台影响进行探究，结果表明:在水下爆炸过程中，气泡、自由面、浮动冲击平台会发生强烈的耦合作用，在气泡脉动阶段，气泡会诱导出涌流和水冢效应，影响浮动冲击平台的安全性和使用性；冲击波是影响浮动冲击平台冲击环境的主要因素，由于气泡的低频性，气泡脉动及水冢对浮动冲击平台的直接冲击作用，会小幅度增加浮动冲击平台冲击环境的谱速度值、谱位移值，对谱加速度值几乎无影响；水冢抨击水面所形成的波浪和气泡破碎兴波，对浮动冲击平台造成的激励载荷呈周期性，其周期与波浪周期相同。波浪的激励载荷仅通过激励其对应频率的浮动冲击平台共振来改变平台的冲击环境。波浪载荷很小，对浮动冲击平台的冲击环境影响较小。     

Modeling air entrainment and transport in a hydraulic jump using two-fluid RANS and DES turbulence models

Both RaNS (Reynolds-averaged Navier–Stokes) and DES (Detached Eddy Simulation) type turbulence models were used in conjunction with a two-fluid model of bubbly flow and a new subgrid air entrainment model to predict air entrainment and transport in a hydraulic jump. It was found that the void fraction profiles predicted by both methods are in agreement with the experimental data in the lower shear layer region, which contains the air bubbles entrained at the so-called toe of the hydraulic jump. In contrast, in the upper roller region behind the toe, the averaged results of the DES turbulence model gives accurate predictions while a RaNS turbulence model does not. This is because the DES turbulence model successfully captures the strong fluctuations on the free surface which allows it to entrain air near the top of the roller region. In contrast, RaNS type turbulence model results in a steady, smooth interface which fails to capture the wave-induced bubble sources in that region. To our knowledge, this study is the first successful quantitative numerical simulation of the overall void fraction profiles in a hydraulic jump.     

Gas entrainment by a liquid film falling around a stationary Taylor bubble in a vertical tube

Gas entrainment by a liquid film falling around a stationary Taylor bubble in a 0.1 m diameter vertical tube is studied experimentally with the purpose of validating a model formulated in an earlier phase of our research. According to this model for a fixed liquid velocity the gas entrainment should be proportional to the waviness of the film (its intermittency) and the wave height and inversely proportional to the film thickness. For Taylor bubble lengths ranging from 1D to 15D these film parameters have been measured with a Laser Induced Fluorescence technique. The gas entrainment has been determined from the net gas flux into the liquid column underneath the Taylor bubble by using data on gas re-coalescence into the rear of the Taylor bubble. These data are available for lengths ranging from 4.5D to 9D. The model results with the measured film characteristics compare well with the observed gas entrainment. The fact that the net gas flux becomes constant for long Taylor bubbles, whereas the wave height still increases, warrants further study.