高雷诺数下非混相Rayleigh-Taylor不稳定性的格子Boltzmann方法模拟

本文采用相场格子Boltzmann方法研究了竖直微通道内中等Atwoods数流体的单模Rayleigh-Taylor不稳定性问题,系统分析了雷诺数对相界面动力学行为以及扰动在各发展阶段演化规律的影响.数值结果表明高雷诺数条件下,不稳定性界面扰动的增长经历了四个不同的发展阶段,包括线性增长阶段、饱和速度阶段、重加速阶段及混沌混合阶段.在线性增长阶段,我们计算获得的气泡与尖钉振幅符合线性稳定性理论,并且线性增长率随着雷诺数的增加而增大.在第二个阶段,我们观察到气泡与尖钉将以恒定的速度增长,获得的尖钉饱和速度略高于Goncharov经典势能模型的解析解[Phys.Rev.Lett.200288134502],这归因于系统中产生了多个尺度的旋涡,而涡之间的相互作用促进了尖钉的增长.随着横向速度和纵向速度的差异扩大,气泡和尖钉界面演化诱导产生的Kelvin–Helmholtz不稳定性逐渐增强,从而流体混合区域出现许多不同层次的涡结构,加速了气泡与尖钉振幅的演化速度,并在演化后期阶段,导致界面发生多层次卷起、剧烈变形、混沌破裂等行为,最终形成了非常复杂的拓扑结构.此外,我们还统计了演化后期气泡与尖钉的无量纲加速度,发现气泡和尖钉的振幅在后期呈现二次增长规律,其增长率系数分别为0.045与0.233.而在低雷诺条件下,重流体在不稳定性后期以尖钉的形式向下运动而轻流体以气泡的形式向上升起.在整个演化过程中,界面变得足够光滑,气泡与尖钉在后期的演化速度接近于常数,未观察到后期的重加速与混沌混合阶段.     

冲击加载下“V”形界面的失稳与湍流混合

采用多组分混合物质量分数模型和最小色散可控耗散格式的高分辨率有限体积方法,数值模拟了弱激波冲击不同角度的"V"形空气/SF_6界面的Richtmyer-Meshkov不稳定性问题.激波冲击界面后,在界面附近沉积涡量,形成沿界面规则排列的旋涡结构,同时界面扰动发展形成气泡和尖钉结构.本文统计了界面左端移动速度和界面混合宽度增长率等特征量的演化规律,并与已有的实验结果进行了对比,两者符合较好.讨论了物质界面处的流体向湍流混合发展的过程,随着界面旋涡结构的演化,涡结构之间开始发生相互诱导、并对等现象,并逐渐聚集在几个区域,而多尺度结构也因旋涡的诱导作用在这些区域中产生.通过对由雷诺数定义的惯性尺度进行分析,发现了具有上下边界的惯性尺度区域的形成,对动能能谱的分析发现了-5/3对数率的出现,这同样说明了惯性尺度区域的形成.由于湍流混合转捩与惯性尺度区域的形成是一致的,界面附近流场将发展为湍流.     

理想磁流体中激波与矩形密度界面相互作用的数值研究

发展一套采用三阶WENO格式和混合GLM方法的理想磁流体数值方法,并对激波与矩形密度界面相互作用进行数值研究.通过圆极化阿尔芬波和旋转激波管问题对数值方法的稳定性和可靠性进行验证.在入射激波马赫数为10,界面内外气体密度比为10的情况,对比不同磁场中矩形密度界面的演变过程.结果表明,磁场能够减少界面上涡量的生成从而抑制界面不稳定性,并且磁场对界面的加速过程以及界面内外气体混合率有影响;而界面的存在将会使波后部分区域磁场增强;由于尖角的存在,矩形界面的发展与圆形界面不同.     

激波汇聚效应对球形气泡演化影响的数值研究

利用三维程序,对比研究了汇聚激波及平面激波冲击下SF6球形气泡演化规律的异同,以期发现激波的汇聚效应对界面演化的影响.三维程序采用多组分可压缩欧拉方程,基于有限体积法,利用MUSCL-Hancock格式进行数值求解,可以达到时间和空间的二阶精度.相比平面激波,汇聚激波由于存在曲率,且激波强度以及壁面效应在汇聚激波运行的过程中逐渐增强,使得激波冲击后的流场演化有较大的不同.计算结果表明:汇聚激波作用下,气泡界面的涡结构更加尖锐;气泡内部的透射激波聚焦程度更强,在界面下游附近形成的最高压力大于平面激波算例,由此产生的射流运动速度更快;由于汇聚激波曲率及激波强度的变化,导致界面上涡量的分布规律以及涡量幅值产生较大变化.通过界面上产生的环量以及界面内外气体混合速度的对比表明,汇聚激波更有助于涡量的产生以及气体的混合.因此激波的汇聚效应对气泡界面演化具有重要影响.     

Schardin问题的数值研究

激波绕过三角楔(Schardin问题)时会产生激波马赫反射与绕射、 三角楔尾涡与涡串等复杂物理现象. 本文利用三阶精度加权基本无振荡(WENO)格式、 结构化矩形网格的自适应加密方法与沉浸边界法对Schardin问题进行了数值模拟. 数值结果清晰地显示了激波与三角楔相互作用, 在楔面发生马赫反射以及在楔角绕射诱导主涡的过程, 并与Schardin等的实验结果及相关数值结果完全符合. 另外, 数值结果还详细反映了先前实验与数值结果没有详细讨论的主涡滑移层上的涡串生成机理, 以及激波与涡串相互作用和产生声波的     

Schardin问题的数值研究

激波绕过三角楔(Schardin问题)时会产生激波马赫反射与绕射、 三角楔尾涡与涡串等复杂物理现象. 本文利用三阶精度加权基本无振荡(WENO)格式、 结构化矩形网格的自适应加密方法与沉浸边界法对Schardin问题进行了数值模拟. 数值结果清晰地显示了激波与三角楔相互作用, 在楔面发生马赫反射以及在楔角绕射诱导主涡的过程, 并与Schardin等的实验结果及相关数值结果完全符合. 另外, 数值结果还详细反映了先前实验与数值结果没有详细讨论的主涡滑移层上的涡串生成机理, 以及激波与涡串相互作用和产生声波的过程.     

磁控条件下激波冲击三角形气柱过程的数值研究

本文基于磁流体动力学方程组,在保证磁场散度为零的条件下,采用CTU+CT（corner transport upwind+constrained transport）算法,对有无磁场控制下激波与重质或轻质三角形气柱相互作用过程进行数值研究.结果表明：无论有无磁场,两气柱在激波冲击下均具有完全不同的波系结构和射流现象.其中,入射激波与重气柱发生常规折射,形成介质射流,而与轻气柱作用则发生非常规折射,形成反相空气射流.无磁场时,气柱在激波冲击下,产生Richtmyer-Meshkov和Kelvin-Helmholtz不稳定性,界面出现次级涡序列,重气柱上下角卷起形成主涡对,轻气柱空气射流穿过下游界面后形成偶极子涡.施加横向磁场后,次级涡序列、主涡对以及偶极子涡均消失.进一步研究表明,在磁场作用下,洛伦兹力将不稳定性诱导产生的涡量向界面两侧的Alfvén波上输运,减少界面涡量沉积,抑制界面卷起失稳.最终,涡量沿界面两侧形成相互远离的涡层,界面不稳定性得到控制.此外,定量分析表明磁场能加快两气柱上游界面的运动,抑制下游界面的运动,且对轻气柱的控制效果更好.     

激波诱导的液柱气液界面RM不稳定性的研究

为了进-步研究气/液界面上的Richtmyer-Meshkov(RM)不稳定性现象,在水平激波管中进行了激波与单、双排液柱相互作用的实验.通过使用驱动气体为氮气、被驱动气体为空气,获得了马赫数为1.10与1.25的激波;通过更换不同尺寸的试件获得了初始直径不同的液柱。使用高速摄影仪记录了激波作用之后液柱失稳变形的详细过程,测量分析了不稳定性发展过程中,尖钉高度、气泡深度及混合区域宽度随时间的发展规律,并与理论分析中的公式进行了比较。     

平面弱激波加载下球形气泡演化的实验研究

 利用高速摄影机直接观测了激波加载下球形气泡的演化过程。采用白色的烟雾颗粒对气泡内透明的测试气体进行染色,直接观察到了气泡界面。通过对平面弱激波加载下轻气体气泡和重气体气泡进行研究,成功验证了轻气体气泡和重气体气泡演化过程中出现的典型界面结构,获得了弱激波加载轻气体气泡中背风涡环的环状几何结构。运用相关方法分析了轻气体气泡实验初期的流场分布,得到的结果与理论分析结果吻合得很好,对于此类实验的数据处理具有一定的启发意义。同时,实验中使用的烟雾颗粒法为今后在球形气泡实验中引入更为精确的平面激光诱导荧光技术（PLIF）、激光粒子图像测速技术（PIV）等实验测试系统提供了布撒示踪粒子的范例。     

应用多介质PPM方法计算斜激波与物质交界面的相互作用

利用多介质PPM方法研究斜激波与物质交界面的相互作用.采用与体积分数耦合的Euler方程组作为计算模型,用双波近似来求解一般刚性气体状态方程Riemann问题.通过体积分数的计算来获得界面的位置,在整个流场采用统一的高阶PPM格式进行计算.文中对斜激波与不同物质界面相互作用进行了数值模拟,并给出了交界面上由于斜压效应产生的涡列的演化过程,特别是强斜激波与不同物质界面的相互作用的情况.     

气固两相圆柱绕流的直接数值模拟

本文运用谱元方法对气固两相圆柱绕流进行了直接数值模拟。在获得高精度计算流场信息基础上,进行颗粒扩散运动的研究。通过研究颗粒大小对颗粒扩散运动的影响来刻画气固两相圆柱绕流的物理特征。研究时发现,由于圆柱尾流独特的涡结构特征,颗粒在流场中心轴线附近区域的扩散机理主要是决定于吸力作用,该作用会驱使小颗粒卷吸进入圆柱尾部近壁区,甚至与圆柱后壁发生碰撞,这种扩散机理与混合层中颗粒扩散机理明显不同。     

EFFECTS OF THE LENGTH SCALE OF FREE-STREAM TURBULENCE AND CYLINDER SIZE ON LOCAL HEAT TRANSFER IN LAMINAR SEPARATED FLOWS

An experimental study has been conducted to investigate the combined influence of length scale of free-stream turbulence, A, and cylinder diameter, D, on heal transfer from slightly heated circular cylinders in cross flow. By using a hot-wire anemometer, spectrum analysis of its signal with a Fourier analyzer is employed to investigate the separated shear layer formed behind the circular cylinder. A resistance thermometer is used to measure the time-averaged temperature behind the cylinder and the power-spectral density of the temperature fluctuation in the thermal boundary layer. It is found that (I) the Nusselt number at the rear stagnation point varies with the relative scale A / D, at the same Reynolds number, (2) the laminar-to-turbulence transition region in the separated shear layer moves in the streamwise direction with the change in the length-scale ratio, and (3) a streamwise shift of the transition region in the separated shear layer causes a substantial increase in the power spectrum of the temperature fluctuation. A correlation equation with A/ Das a parameter is derived to predict the heat transfer rate at the rear stagnation point.     

Structure of detonation waves in pressed 2′,2′,2′-trinitroethyl-4,4,4-trinitrobutyrate

The reaction zones of steady detonation waves in pressed 2′,2′,2′-trinitroethyl-4,4,4-trinitrobutyrate (TNETB) charges with densities within 1.23–1.79 g/cm³ and initial powder dispersities of 5 and 80 μm were studied with the use of a VISAR interferometer. It was demonstrated that the boundaries of the initial density range within which a pressure growth is observed instead of a chemical spike, predicted by the theory, are determined by the HE dispersity. The measured dependence of the detonation velocity on the initial HE density exhibits features suggestive of an underdriven detonation wave without a chemical spike. The peculiarities observed are accounted for by the heterogeneous structure of pressed HEs, the decomposition of which is of hot-spot character, occurring partly in the compression wave front.     

利用低杂波改善约束的实验研究

介绍了在HT-6M托卡马克上,利用双低杂波的组合,成功地实现了准稳态的高约束模式运行. 能量约束时间提高两倍以上,粒子的约束在较高密度下依然改善3倍以上.通过应用波扩散及 电流径向扩散方程计算低杂波电流传播的方法,对一组典型的数据进行数值模拟.计算表明 ,在HT-6M低杂波实验中,由于纵场较低,密度较高,低杂波的能量沉积在离磁轴较远的位 置,使等离子体电流密度分布成为反剪切位形,内部输运垒地形成,大大提高了等离子体的 约束状况.实验数据也给出了反剪切的证据. 关键词： 托卡马克 低杂波 约束改善     

Materials science of the gel to fluid phase transition in a supported phospholipid bilayer

We report the results of in situ AFM measurements examining the phase transition of bilayers formed from the zwitterionic phospholipid, DMPC, 1,2-dimyristoyl-sn-glycero-3-phosphocholine, supported on mica. The images show that the fluid to gel phase transition process features substantial tearing of the bilayer due to the density change between the two phases. The gel to fluid transition is strongly affected by the resultant stress introduced into the gel phase, which changes the degree of cooperativity, the shape of developing fluid phase regions, and the course of the transition.     

Turbulent structures in an optimal Taylor–Couette flow between concentric counter-rotating cylinders

The turbulent structures formed in a Taylor–Couette (TC) flow established between two concentric counter-rotating cylinders were explored numerically. The shear Reynolds number was set to Re_shear = 8000 and the radius ratio was set to r_i/r_o = 0.5. An optimal flow corresponding to the maximal angular velocity transport between the cylinders was selected for the TC flow. The mean velocity profile reached its steepest value near the cylinders in the optimal TC flow. The streamwise velocity correlations at the outer cylinder in the gap exceeded those at the inner cylinder. The large convective transport of angular velocity in the gap generated a maximal angular velocity flux to achieve the optimal flow. The angular velocity flux generated by the momentum source exceeded that generated by the momentum sink. The vorticity dispersion was larger near the inner cylinder than near the outer cylinder, but vorticity stretching near the outer cylinder exceeded than that near the inner cylinder. The skin friction coefficient budgets were plotted using the velocity–vorticity correlation. The vortex stretching contributions dominated the skin friction budgets. The area near the inner cylinder was populated by stronger vortices, but their population density was smaller than the population density of the vortices near the outer cylinder. The probability density functions of the wall-normal and streamwise velocity fluctuations delineated the presence of the large wall-normal velocity fluctuations near the outer cylinder.     

Sound scattering by a Karman street consisting of large-scale vortices

The scattering of acoustic waves by a vortex street formed behind a cylinder in an air flow is studied both theoretically and experimentally for the case of the sound wavelength being much less than the vortex size. The theoretical calculations show that, at flow velocities well below the sound velocity, the vortex street can be considered as a moving phase screen. The spectrum of scattered sound in the far zone is shown to consist of harmonics whose frequencies differ by a multiple of the vortex rate. The computational results agree well with the experimental data obtained for the diffraction of ultrasound of the wavelength λ=3 mm by the Karman street formed behind a circular cylinder with an 8 mm diameter at a flow velocity of 7 m/s.     

Three-dimensional numerical simulation of crown spike due to coupling effect between bubbles and free surface

The motion of gas bubbles beneath a free surface will lead to a spike of fluid on the free surface. The distance of the bubbles to the free surface is the key factor to different phenomena. When the inception distance varies in some range, crown phenomenon would happen after the impact of weak buoyancy bubbles, so this kind of spike is defined as crown spike in the present paper. Based on potential flow theory, a three-dimensional numerical model is established to simulate the motion of the free-surface spike generated by one bubble or a horizontal line of two in-phase bubbles. After the downward jet formed near the end of the collapse phase, the simulation of the free surface is performed to study the crown spike without regard to the toroidal bubble＇s effect. Calculations about the interaction between one bubble and free surface agree well with the experimental results conducted with a high-speed camera, and relative error is within 15%. Crown spike in both single- and two-bubble cases are simulated numerically. Different features and laws of the motion of crown spike, depending on the bubble-boundary distances and the inter-bubble distances, have been investigated.     

Fraunhofer Diffraction Pattern by a Slit Aperture between a Knife Edge and a Metal Cylinder

The Fraunhofer diffraction pattern of an aperture formed between a knife edge and a smooth surface of a turned metal cylinder or half cylinder is different from the diffraction pattern of an ideal slit aperture. The diffraction pattern is simulated by adding the light reflected by the cylinder surface to the pattern formed by a slit, and is measured using a cylinder that has been partly cut away so that it can represent either a full or half cylinder. The theoretical and measured patterns agree well with each other.     

Properties of surface waves in an elastic cylinder filled with a liquid

The properties of harmonic surface waves in an elastic cylinder filled with a liquid are studied. The case of elastic material for which the shear wave velocity is higher than the sound velocity in a liquid is considered. The wave motion is described based on the complete system of equations of the dynamic theory of elasticity and the equation of motion of an ideal compressible liquid. The asymptotic analysis of the dispersion equation in the region of large wave numbers and qualitative analysis of the dispersion spectrum showed that in such a waveguiding system there exist two surface waves, the Stoneley and the Rayleigh waves. The lowest normal wave forms the Stoneley wave on the internal surface of the cylinder. In this waveguide phase, velocities of all normal waves, except for the lowest one, have the velocity of sound in the liquid as their limit. Therefore, the Rayleigh wave on the external surface of the cylinder is formed by all normal waves in the range of frequencies and wave numbers in which phase velocities of normal waves of the composite waveguide and the lowest normal wave of the elastic hollow cylinder coincide.