Solitary waves on vortex lines in Ginzburg-Landau models for the example of Bose-Einstein condensates

Axisymmetric disturbances that preserve their form as they move along the vortex lines in uniform Bose-Einstein condensates are obtained numerically by the solution of the Gross-Pitaevskii equation. A continuous family of such solitary waves is shown in the momentum (p)-substitution energy (Epsilon) plane with p-->0.09 rho kappa(3)/c(2), Epsilon-->0.091 rho kappa(3)/c as U-->c, where rho is the density, c is the speed of sound, kappa is the quantum of circulation, and U is the solitary wave velocity. It is shown that collapse of a bubble captured by a vortex line leads to the generation of such solitary waves in condensates. The various stages of collapse are elucidated. In particular, it is shown that during collapse the vortex core becomes significantly compressed, and after collapse two solitary wave trains moving in opposite directions are formed on the vortex line.     

油滴撞击油膜层内气泡的变形与破裂过程的数值模拟

旋转工作的机械零部件和机械设备的润滑系统工作过程中普遍存在着油滴和油膜的碰撞行为,这一行为易引起气泡夹带现象.气泡将对油滴撞击油膜时的运动过程和附壁油膜层的形成质量造成不可忽视的影响.基于耦合的水平集-体积分数方法,对油滴撞击含气泡油膜的行为进行数值模拟研究,考察油膜层内气泡的变形运动过程,分析气泡大小和位置等因素对撞击过程中气泡变形特征参数的影响规律,并探讨气泡破裂的动力学机制.研究表明,随着气泡直径的增大,油滴撞击含气泡油膜后气泡会依次出现自由表面破裂、稳定变形以及油膜内部破裂等现象;直径d=20μm的气泡能较稳定地存在于油膜层内,同时该值也是气泡发生自由表面破裂和油膜内部破裂的临界值.此外,气泡所在位置同样对气泡变形历程有一定影响,气泡越接近油膜表面,其变形量越大;位于油膜底层的气泡会附着在壁面上.在自由表面破裂和油膜内部破裂过程中,气泡破裂是由气-液界面不稳定引起的,表面张力对这两种现象起重要作用;而黏性剪切力对油膜内部破裂现象也有着不可忽视的影响.     

Lagrangian coherent structures analysis of gas-liquid flow in a bubble column

The objective of this paper is to apply a new identifying method to investigating the gas-liquid two-phase flow behaviors in a bubble column with air injected into water.In the numerical simulations,the standard k-?turbulence model is employed to describe the turbulence phenomenon occurring in the continuous fluid.The Finite-Time Lyapunov Exponent(FTLE)and Lagrangian Coherent Structures(LCS)are applied to analyze the vortex structures in multiphase flow.Reasonable agreements are obtained between the numerical and experimental data.The numerical results show that the evolution of gas-liquid in the column includes initial and periodical developing stages.During the initial stage,the bubble hose is forming and extending along the vertical direction with the vortex structures formed symmetrically.During the periodical developing stage,the bubble hose starts to oscillate periodically,and the vortexes move along the bubble hose to the bottom of column alternately.Compared to the Euler-system-based identification criterion of a vortex,the FTLE field presents the boundary of a vortex without any threshold defined and the LCS represents the divergence extent of infinite neighboring particles.During the initial stage,the interfaces between the forward and backward flows are highlighted by the LCS.As for the periodical developing stage,the LCS curls near the vortex centers,providing a method of analyzing a flow field from a dynamical system perspective.     

过热液体中蒸汽泡上升过程的格子Boltzmann三维数值模拟

应用格子Boltzmann相变模型,在三维空间研究蒸汽泡在过热液体中生长、上升和变形等动力学行为.为研究传热传质对蒸汽泡运动的影响,对比模拟相同条件下气泡在等温环境中上升的物理过程.结果表明：蒸汽泡在过热液体中上升发生的变形程度较小,意味着相变对蒸汽泡的影响和表面张力一样使汽泡保持初始的形状.蒸汽泡在过热液体中的上升速度较小,说明随着汽泡生长拖拽力的影响比浮力大.蒸汽泡生长率在初始阶段达到最大值,随后会趋于一个恒定的值.随着汽泡体积增大和上升速度的增加,其对流场的扰动也越来越剧烈.蒸汽泡生长和上升引起的对流运动对温度场的演化造成很大的影响.     

入射和反射激波诱导重气泡变形和失稳的三维数值研究

 采用大涡模拟方法,对入射激波及其反射激波诱导球形重气泡的变形失稳过程进行了三维数值模拟,利用已有实验验证了计算模型的可靠性,重点考察了反射激波与已经失稳的气泡界面的再次作用,讨论了涡环的形成及其三维失稳的过程。研究结果显示：入射和反射激波与球形重气泡作用产生斜压效应,会在流场中产生旋转方向截然相反的多个涡环;反射激波诱导的涡环具有较小的强度,故更加容易失稳,甚至能完全形成具有流向涡量的复杂小尺度涡结构。     

基于光学涡旋相移技术的离面位移测量

设计了基于光学涡旋相移技术的离面位移测量实验方案,实现了电子散斑干涉中相移的数字控制.该方法利用输入液晶空间光调制器中的叉形光栅产生涡旋光束,通过涡旋光束绕轴的旋转产生相移;同时,产生的涡旋光束又作为参考光与物光干涉.实验中,在物体发生离面位移前后依次输入四幅叉形光栅,产生相移步长为π/2的涡旋光束,利用CCD获得涡旋光与物光的干涉光场,从而获得离面位移场的包裹相位;再通过解包裹,获得物体离面变形的相位变化.光学涡旋相移法可应用于离面位移测量.     

Numerical Simulations of Two Coaxial Vortex Rings Head-on Collision

Vortex rings have been a subject of interest in vortex dynamics due to a plethora of physical phenomena revealed by their motions and interactions within a boundary. The present paper is devoted to physics of a head-on collision of two vortex rings in three dimensional space, simulated with a second order finite volume scheme and compressible. The scheme combines non-iterative approximate Riemann-solver and piecewise-parabolic reconstruction used in inviscid flux evaluation procedure. The computational results of vortex ring collisions capture several distinctive phenomena. In the early stages of the simulation, the rings propagate under their own self-induced motion. As the rings approach each other, their radii increase, followed by stretching and merging during the collision. Later, the two rings have merged into a single doughnut-shaped structure. This structure continues to extend in the radial direction, leaving a web of particles around the centers. At a later time, the formation of ringlets propagate radially away from the center of collision, and then the effects of instability involved leads to a reconnection in which small-scale ringlets are generated. In addition, it is shown that the scheme captures several experimentally observed features of the ring collisions, including a turbulent breakdown into small-scale structures and the generation of small-scale radially propagating vortex rings, due to the modification of the vorticity distribution, as a result of the entrainment of background vorticity and helicity by the vortex core, and their subsequent interaction.     

Numerical prediction of tip vortex cavitation behavior and noise considering nuclei size and distribution

The tip vortex cavitation behavior and sound generation were numerically analyzed. A numerical scheme combining Eulerian flow field computation and Lagrangian particle trace approach was applied to simulate tip vortex cavitation. Flow field was computed by using hybrid method which combines Reynolds-Averaged Navier-Stokes solver with Dissipation Vortex Model. The trajectory and behavior of each cavitation bubble were computed by Newton’s second law and Rayleigh-Plesset equation, respectively. According to nuclei population data, the cavitation nuclei were distributed and convected into the tip vortex flow. Calculated volume of the cavitation bubble and the trajectory were used as the input of cavitation bubble noise analysis. The relationship of cavitation inception, sound pressure level, and cavitation nuclei size was studied at several cavitation numbers. It was found that cavitation inception of smaller nuclei is more sensitive to the change of cavitation number and cavitation noise due to the cavitated smallest nuclei has the most influence on overall tip vortex cavitation noise.     

Shape stability of a gas bubble in a soft solid

Predicting the onset of non-spherical oscillations of bubbles in soft matter is a fundamental cavitation problem with implications to sonoprocessing, polymeric materials synthesis, and biomedical ultrasound applications. The shape stability of a bubble in a Kelvin-Voigt viscoelastic medium with nonlinear elasticity, the simplest constitutive model for soft solids, is analytically investigated and compared to experiments. Using perturbation methods, we develop a model reducing the equations of motion to two sets of evolution equations: a Rayleigh-Plesset-type equation for the mean (volume-equivalent) bubble radius and an equation for the non-spherical mode amplitudes. Parametric instability is predicted by examining the natural frequency and the Mathieu equation for the non-spherical modes, which are obtained from our model. Our theoretical results show good agreement with published experiments of the shape oscillations of a bubble in a gelatin gel. We further examine the impact of viscoelasticity on the time evolution of non-spherical mode amplitudes. In particular, we find that viscosity increases the damping rate, thus suppressing the shape instability, while shear modulus increases the natural frequency, which changes the unstable mode. We also explain the contributions of rotational and irrotational fields to the viscoelastic stresses in the surroundings and at the bubble surface, as these contributions affect the damping rate and the unstable mode. Our analysis on the role of viscoelasticity is potentially useful to measure viscoelastic properties of soft materials by experimentally observing the shape oscillations of a bubble.     

单液滴撞击倾斜液膜飞溅过程的耦合Level Set-VOF模拟

采用耦合水平集--体积分数法(CLSVOF)对液滴撞击倾斜表面液膜后液膜的形态演化及飞溅过程进行数值模拟, 并对液滴撞击液膜过程中形成的空气卷吸现象进行研究并探讨了撞击角对此的影响, 分析了液滴撞击后液体内部的压力和速度分布, 对液滴撞击倾斜表面液膜的飞溅过程进行讨论, 并与实验结果进行了对比, 验证了CLSVOF方法研究液滴撞击倾斜液膜的可行性. 结果表明, 液滴撞击倾斜液膜时前后两部分飞溅现象产生的机理不同, 前半部分飞溅是由于压差引起的颈部射流, 而后半部分则是由液膜径向流动产生的飞溅现象. 随着撞击角的增大, 空气卷吸气泡数量减少.     

Flow visualization of vortex structure in a pulsed rectangular jet

Pulsating jet is visualized using hydrogen bubble method to clarify the vortex nature in the near field of the jet. This study focused on the development in space and time of vortex structures evolution in low aspect-ratio rectangular jet with pulsation. Pulsation means large-amplitude, low-frequency excitation which is expected to increase the mixing and spreading of the jet and to accelerate its transition from a rectangular form to an axisymmetric form. It was deemed appropriate to investigate whether jet characteristics of a pulsating, submerged jet flow can be altered by including pulsations. The difference of the vortex deformation process is discussed in relation to pulsating conditions. Consequently, the pulsation leads to the formation of vortices at regular intervals, which are larger than those occurring in a steady jet. The results show that the streamwise interaction, between leading vortex and trailing vortex rolled up at nozzle lips, strengthens with increasing pulsating frequency. The spanwise drift of the vortex becomes strongly apparent at large amplitude and high frequency conditions. The drifting start position does not change regardless of pulsating condition. The convection velocity of vortex increases at lower frequency and larger amplitude.     

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

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

Numerical simulation of bubble dynamics in a Phan-Thien–Tanner liquid: Non-linear shape and size oscillatory response under periodic pressure

We study non-linear bubble oscillations driven by an acoustic pressure with the bubble being immersed in a viscoelastic, Phan-Thien–Tanner liquid. Solution is provided numerically through a method which is based on a finite element discretization of the Navier–Stokes flow equations. The proposed computational approach does not rely on the solution of the simplified Rayleigh–Plesset equation, is not limited in studying only spherically symmetric bubbles and provides coupled solutions for the velocity, stress fields and bubble interface. We present solutions for non-spherical bubbles, with asphericity being addressed by means of Legendre polynomials or associated Legendre functions. A parametric investigation of the bubble dynamical oscillatory response as a function of the fluid rheological properties shows that the amplitude of bubble oscillations drastically increases as liquid elasticity (quantified by the Deborah number) increases or as liquid viscosity decreases (quantified by the Reynolds number). Extensive numerical calculations demonstrate that increasing elasticity and/or viscosity of the surrounding liquid tend to stabilize the shape anisotropy of an initially non-spherical bubble. Results are shown for pressure amplitudes 0.2–2 MPa and Deborah, Reynolds numbers in the intervals of 1–8 and 0.094–1.256, respectively.     

Large eddy simulation of bluff-body stabilized swirling non-premixed flames

Large eddy simulations (LES) using a subgrid mixing and combustion model are carried out to study two bluff-body stabilized swirling non-premixed flames (SM1 and SMA2). The similarities and differences between the two flames are highlighted and discussed. Flow features, such as, the recirculation zone (RZ) size and the flame structure are captured accurately in both cases. The SM1 flame shows a toroidal RZ just behind the bluff body and a vortex breakdown bubble (VBB) downstream. In addition, a highly rotational non-recirculating region in-between the RZ and VBB is observed as well. On the other hand, the SMA2 shows a single elongated recirculation zone downstream the bluff body. Flame necking is observed downstream the bluff body for the SM1 flame but not for the SMA2 flame. The time-averaged velocity and temperature comparison also shows reasonable agreement. The study shows that the sensitivity of the flame structure to inflow conditions can be captured in the present LES without requiring any model changes.     

Phase field simulation of single bubble behavior under an electric field

Based on the Cahn-Hilliard phase field model, a three-dimensional multiple-field coupling model for simulating the motion characteristics of a rising bubble in a liquid is established in a gas-liquid two-phase flow. The gas-liquid interface motion is simulated by using a phase-field method, and the effect of the electric field intensity on bubble dynamics is studied without electric field, or with vertical electric field or horizontal electric field. Through the coupling effect of electric field and flow field, the deformation of a single rising bubble and the formation of wake vortices under the action of gravity and electric field force are studied in detail. The correctness of the results is verified by mass conservation, and the influences of different electric field directions and different voltages on the movement of bubbles in liquid are considered. The results show that the ratio of the length to axis is proportional to the strength of the electric field when the air bubble is stretched into an ellipsoid along the electric field line under the action of electrostatic gravity and surface tension. In addition, the bubble rising speed is affected by the electric field, the vertical electric field accelerates the bubble rise, and the horizontal direction slows it down.     

入射和反射激波诱导轻气泡变形和失稳的计算

采用NS方程对激波诱导球形轻气泡的变形失稳过程进行三维数值模拟,验证本文模型的可靠性,讨论气泡变形、涡环形成及其三维失稳过程.结果发现:球形轻气泡在入射和反射激波先后作用下,可以变形并沿流向形成多个旋转方向不同的涡环,而斜压效应是产生这些涡环的主要原因,涡环不仅在自诱导效应和流场作用下运动,而且自身会受到扰动而发生方位角不稳定.当涡环形成以小尺度流向涡为主的复杂结构时,即有可能诱发流场湍流化.     

非球形水下爆炸气泡坍塌机制

 为研究非球形水下爆炸气泡的坍塌机制和射流特性,采用水下爆炸实验和数值模拟相结合的手段,研究了典型的柱形装药自由场水下爆炸气泡的运动过程。实验结果与数值模拟结果符合较好,通过对比发现,非球形水下爆炸气泡坍塌所形成的射流与球形气泡坍塌射流的行为存在差异。研究结果表明：非球形气泡的坍塌过程与气泡初始形态的局部曲率半径有关,初始气泡表面曲率半径大的区域的初始膨胀速度较曲率半径小的区域的小,坍塌时刻也较早,导致气泡产生非球形坍塌。     

多晶银纳米线拉伸变形的分子动力学模拟研究

纳米尺度金属Ag以其独特的导电和导热性,广泛应用于微电子、光电子学、催化等领域,特别是在纳米微电极和纳米器件方面的应用.本文采用分子动力学方法模拟了不同晶粒尺寸下多晶银纳米线的拉伸变形行为,详细分析了晶粒尺寸对多晶银纳米线弹性模量、屈服强度、塑性变形机理的影响.发现当晶粒尺寸小于13.49 nm时,多晶Ag纳米线呈现软化现象,出现反Hall-Petch关系,此时的塑性变形机理主要以晶界滑移、晶粒转动为主,变形后期形成五重孪晶;当晶粒尺寸大于13.49 nm时,塑性变形以位错滑移为主,变形后期产生大量的孪晶组织.     

A study of cavitation nucleation in pure water using molecular dynamics simulation

To discover the microscopic mechanism responsible for cavitation nucleation in pure water, nucleation processes in pure water are simulated using the molecular dynamics method. Cavitation nucleation is generated by uniformly stretching the system under isothermal conditions, and the formation and development of cavitation nuclei are simulated and discussed at the molecular level. The processes of energy, pressure, and density are analyzed, and the tensile strength of the pure water and the critical volume of the bubble nuclei are investigated. The results show that critical states exist in the process of cavitation nucleation. In the critical state, the energy, density, and pressure of the system change abruptly, and a stable cavitation nucleus is produced if the energy barrier is broken and the critical volume is exceeded. System pressure and water density are the key factors in the generation of cavitation nuclei. When the critical state is surpassed, the liquid is completely ruptured, and the volume of the cavitation nucleus rapidly increases to larger than 100 nm³; at this point, the surface tension of the bubble dominates the cavitation nucleus, instead of intermolecular forces. The negative critical pressure for bubble nucleation is -198.6 MPa, the corresponding critical volume is 13.84 nm³, and the nucleation rate is 2.42×10³² m^-3·^-1 in pure water at 300 K. Temperature has a significant effect on nucleation: as the temperature rises, nucleation thresholds decrease, and cavitation nucleation occurs earlier.