Non-linear analysis of creeping flow on the inclined permeable substrate plane subjected to an electric field

The effect of an externally applied electric field on the stability of a thin fluid film over an inclined porous plane is analyzed using linear and non-linear stability analysis in the long wave limit. The principle aim of this study is to illustrate the influence of electric field on the non-linear stability of a thin liquid layer flow down incline substrate when the plane is porous. The driving force for the instability under an electric field is an electrostatic force exerted on the free charges accumulated at the dividing interface. The coupled non-linear evolution equations for the local film thickness and the interfacial charge for two-dimensional disturbances are derived to analyze the effect of long-wave instabilities. The method of multiple scales is applied to obtain approximate solutions and analyze the stability criteria. Numerical simulations of this system of non-linear evolution equations are performed. It is found that the permeability parameter as well as the inclination of the plane plays a destabilizing role in the stability criteria, while the damping influence is observed for increasing of the electrical conductivity in both linear and non-linear behavior.     

The electrohydrodynamics of superimposed fluids subjected to a nonuniform transverse electric field

This study aims to investigate electrohydrodynamics of two superimposed fluids that are confined between a pair of two-dimensional flat plates and are exposed to a sinusoidal electric field in zero gravity. The goal is to identify the parameters that affect the flow structure and interface deformation using a simple closed form solution. The governing electrohydrodynamic equations are solved analytically for Newtonian and immiscible fluids in the framework of leaky-dielectric theory and in the limit of small electric field and fluid inertia. A detailed analysis of the electric and flow fields is presented and it is shown that the electric field induces sinusoidal electrical stresses at the interface, which lead to periodic convection cells. The parameters affecting the sense of flow circulation and strength are investigated and it is shown that the former depends on the relative magnitude of the electric permittivity and conductivity ratios while the latter is controlled by the relative thicknesses of the fluid layers and the ratio of the electric conductivities and viscosities of the fluids. The maximum flow strength is achieved at a relative thickness that is set by the competition between the electric and hydrodynamic effects. For small deformation, the distortion of the interface is examined using a normal stress balance at the interface, and it is shown that the degree of interface deformation scales with the square of the amplitude of the electric potential nonuniformity, while its wavenumber is twice that of the imposed potential nonuniformity. Furthermore, a zero-deformation curve is found, which delineates the region in the permittivity-conductivity space according to the sense of interface deformation. The results show that for certain ranges of fluid layer thicknesses and permittivity ratios, the interface will remain flat, despite the action of the nonuiform field.     

Modulation instability of electrohydrodynamic waves on liquid jet

An analytical study of slow modulation has been made of cylindrical interface between two inviscid streaming fluids, in the presence of a relaxation of electrical charges at the interface, and stressed by an axial electric field. A new technique based on the perturbation theory, to derive the non-linear evolution equations has been introduced. These equations are combined to yield a non-linear Ginzburg–Landau equation and a non-linear modified Schrödinger equation describing the evolution of wave packets. The linear analysis showed that the streaming has a destabilizing effect and the electric field has stabilizing influence associated with parameters condition involving the electric conductivity and permittivity of the fluids. While the non-linear approach indicated that the streaming may become unstable for sufficiently high velocities, with a new condition on the material properties, involving weak electric relaxation times in both fluids.     

Onset of shear-layer instability at the interface of parallel Couette flows

A non-planar or a bilateral mixing-layer is studied by means of a series of direct numerical simulations (DNSs). This mixing-layer forms at the interface of two co-current plane Couette flows of different Reynolds numbers. The current DNS study determined the conditions for the onset of shear-layer instability at the interface. The influence of different Reynolds number (of the co-current plane Couette flows) and their Reynolds number ratio on the mixing-layer is studied. A critical Reynolds number of about 500 (or more particularly one of the co-current plane Couette flows must be turbulent) and a Reynolds number ratio greater than 2 is required for the genesis of this bilateral shear-layer instability. Independent of the Reynolds number and the Reynolds number ratio, the temporal evolution of the shear-layer instability followed the same pattern. In addition, the oscillation frequency of the instability was found to increase with increasing Reynolds number and increasing Reynolds number ratio. Further, influence of instability on the local skin friction and the two-point correlation is elaborated on.     

Electrohydrodynamic instability of a dielectric compressible liquid sheet streaming into an ambient stationary compressible gas

The effect of compressibility of fluids on the linear electrohydrodynamic instability of a dielectric liquid sheet issued from a nozzle into an ambient dielectric stationary gas in the presence of a horizontal electric field is investigated. It is found that increasing the Mach number from subsonic to transonic causes the maximum growth rate and the dominant wavenumber of the disturbances to increase, and the increase is higher in the presence of the electric field. Liquid compressibility has been found to have a minimal effect on instability. At constant wavenumber and electric field values, the growth rate of disturbances increases as the gas Mach number tends to 1, and then begins to decrease with further increase in the gas Mach number. At small values of wavenumber, antisymmetrical disturbances grow faster than symmetrical ones, while the growth rate of both types of disturbances approach each other at large wavenumbers, which increases by increasing the electric field values. At small Weber numbers, antisymmetrical disturbances exhibit a higher maximum growth rate and a lower dominant wavenumber than symmetrical disturbances. However, the maximum growth rate and dominant wavenumber of the two types of disturbances are almost identical when both Weber number and electric field values become large. An increase in the gas to liquid density ratio enhances the instability, and this effect is enhanced for high electric field values. Surface tension and electric fields always oppose and increase the development of instability, respectively; and they have opposite effects for long wavelengths and high Weber numbers.     

横向交流电场下液膜参数不稳定性分析

当将运动的平面液膜置于横向的交流电场之间时会产生参数振荡现象.为了得到交流电场下平面液膜的色散关系并为液膜的破碎行为分析提供理论基础,本文基于漏电介质模型对液体的电学特性进行假设,对平面液膜在直流和交流电场下的参数不稳定性进行了分析.由于主流是基于时间的流动, 在稳定性分析中引入了Floquet理论. 在文中,将电场定义为部分交流电场和部分直流电场共同耦合后的混合型电场. 最后,波数和不稳定增长率之间的无量纲色散方程可由矩阵的形式表示.本文考虑了多种参数对不稳定的影响, 包括气液密度比、韦伯数、雷诺数、欧拉数、 松弛时间以及衡量交流电场占比的参数及频率参数,并得知欧拉数同时影响毛细不稳定性及参数不稳定性,交流电场占比对不稳定性的影响体现在恒定电场力上, 而交流电场频率主要影响参数不稳定性. 为了在实验中更容易地寻求参数振荡现象,增大电欧拉数及减小交流电场频率是有效的方法.      

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     

Boundary effects on the Bénard-Marangoni instability under an electric field

This article theoretically studies the Bénard-Marangoni instability problem for a liquid layer with a free upper surface, which is heated from below by a heating coil through a solid plate in ana.c. electric field. The boundary effects of the solid plate, which include its thermal conductivity, electric conductivity and thickness, have great influences on the onset of convective instability in the liquid layer. The stability analysis in this study is based on the linear stability theory. The eigenvalue equations obtained from the analysis are solved by using the fourth order Runge-Kutta-Gill's method with the shooting technique. The results indicate that the solid plate with a higher thermal or electric conductivity and a bigger thickness tends to stabilize the system. It is also found that the critical Rayleigh numberR _c, the critical Marangoni numberM _c, and the criticala.c. Rayleigh numberE _ac become smaller as the intensity of thea.c. electric field increases.     

双流体雾化液滴荷电特性分析

本文通过理论分析,提出了双流体荷电雾化过程中,增大电场强度或减小液滴粒径均能提高荷电效果的论点。针对不同的电极布置情况建立了数值模拟模型,计算结果表明,双流体喷嘴接地时,喷嘴与环状电极之间的电场强度比喷嘴不接地时大很多。通过实验,验证了电场强度、液滴粒径等参数对荷电效果的影响规律,在环状电极荷电情况下,双流体喷嘴接地时,液体荷质比高于双流体喷嘴不接地时的荷质比;电场强度增大,荷质比增加;液滴粒径越小,荷质比越大;在一定的电导率范围内,电导率越大,荷质比越小。实际喷雾中,运动液滴发生二次破碎的临界荷电量小于Rayleigh极限。     

电迁移诱发表面扩散下沿晶微裂纹的演化

随着微电子技术的快速发展,铜内连导线的失效问题日益受到关注。基于表面扩散和蒸发—凝结的经典理论及其弱解描述,建立描述电迁移下微结构演化的有限单元法,对铜材料内沿晶微裂纹在电迁移诱发表面扩散下的不稳定外形演化进行了有限元模拟。详细讨论了电场、形态比和晶界能与表面能比值对沿晶微裂纹演化的影响。结果表明：沿晶微裂纹在沿着晶界和电场方向发生漂移的过程中存在分节与不分节两种形貌演化分叉趋势,且裂腔分节存在临界形态比 和临界电场值 。当 或 时,沿晶微裂纹会沿着晶界分节成一大一小两个小沿晶微裂纹。沿晶微裂纹分节时间随形态比和电场的增大而减小,即形态比和电场的增大都使沿晶微裂纹加速分节。而临界电场值 随着形态比的增大而减小;临界形态比 随电场的增大而减小。也就是说形态比和电场的增大还将有助于沿晶微裂纹分节。此外,沿晶微裂纹分节时间要比晶内微裂纹的小,即晶界的存在有助于加速裂腔分节。     

Direct numerical simulations of two-layer viscosity-stratified flow

Two-dimensional simulations of flow instability at the interface of a two-layer, density-matched, viscosity-stratified Poiseuille flow are performed using a front-tracking/finite difference method. We present results for the small-amplitude (linear) growth rate of the instability at small to medium Reynolds number for varying thickness ratio n, viscosity ratio m, and wavenumber. We also present results for large-amplitude non-linear evolution of the interface for varying viscosity ratio and interfacial tension. For the linear case, the interfacial mode is neutrally stable for as predicted by analysis. The growth rate is proportional to Reynolds number for small Re, and increases with viscosity ratio. The growth rate also increases when the thickness of the more viscous layer is reduced. Strong non-linear behavior is observed for relatively large initial perturbation amplitude. The higher viscosity fluid is drawn out as a finger that penetrates into the lower viscosity layer. The simulated interface shape compares well with previously reported experiments. Increasing interfacial tension retards the growth rate of the interface as expected, whereas increasing the viscosity ratio enhances it. Drop formation at the small Reynolds number considered in this study is precluded by the two-dimensional nature of the calculations.     

Effect of streamwise and spanwise electric fields on transient growth in a two-fluid channel flow

A linear model of a two-fluid channel flow under streamwise/spanwise electric field is built. Both the fluids are assumed to be incompressible, viscous and perfectly dielectric. The effect of the streamwise and spanwise electric fields on transient behavior of small three-dimensional disturbances is studied. The numerical result shows that the streamwise electric field suppresses transient growth of the disturbance with spanwise uniform wave number. The spanwise electric field diminishes transient growth of the disturbance with streamwise uniform wave number. Two peaks of optimal growth are detected in the wave number plane. The peak at relatively large spanwise wave number is dominated by the lift-up mechanism and little influenced by electric field. Differently, the peak at relatively small wave number is associated with the characteristic of the interface and possibly influenced by electric field. The effect of the Weber number, the Reynolds number and the relative electrical permittivity on optimal growth is studied as well. A scaling law is obtained for relatively small Weber numbers and relatively large Reynolds numbers.     

激波与椭圆形重气柱相互作用的PLIF实验

在水平激波管中,采用平面激光诱发荧光(planar laser-induced fluorescence, PLIF)方法对椭圆形重气柱界面的Richtmyer-Meshkov不稳定性进行实验。气柱由SF₆混入一定比例的丙酮蒸气构成,环境气体为空气。通过改变椭圆形气柱的长短轴比值,得到了激波马赫数为1.25时,3种初始界面的演化形态。通过相对体积分数标定,得到了界面失稳演化过程中的相对体积分数分布,观察到了激波作用后界面气体聚集、转移、消散等现象。实验结果发现,对于流向轴长与展向轴长之比较大的气柱界面,初始界面产生的涡量更大且分布更广,其界面不稳定性发展得越迅速和剧烈。失稳发展迅速的界面甚至出现涡对碰撞并产生尾部射流结构的现象。初始界面直接决定了失稳发展初期形成的涡对强度和间距,并对后期演化有重要影响。     

Axisymmetric and non-axisymmetric instability of an electrically charged viscoelastic liquid jet

A linear analysis is performed to investigate the competition between axisymmetric and non-axisymmetric instability of an electrically charged viscoelastic liquid jet. The liquid is assumed to be a dilute polymer solution modeled by the Oldroyd-B constitutive equation. As to its electric properties, the liquid is assumed to be of finite electrical conductivity and is described by the Taylor–Melcher leaky dielectric theory. An analytical dispersion relation is derived and the temporal growth rate is solved numerically. Two viscoelastic liquids, i.e. a PEO aqueous solution and a PIB Boger fluid, are taken as examples to study the effects of electric field and electrical conductivity on jet instability. The result shows that electric field basically destabilizes both the axisymmetric and the non-axisymmetric mode. On the other hand, the effect of electrical conductivity on the modes is quite limited. An energy analysis shows that elasticity enhances both axisymmetric and non-axisymmetric jet instability and its destabilization effect on the axisymmetric mode is more profound. For viscoelastic jets of high Deborah numbers the combined effect of viscosity and elasticity is possibly characterized by an equivalent Reynolds number related only to the viscosity of solvent.     

Thin liquid film morphology driven by electro-static field

The development of stationary patterns on a thin polymer surface subject to an electric field is studied by means of the hexagonal-planform weakly nonlinear stability analysis and numerical simulations.The time evolution of the interface between the air and the polymer film on the unbounded spatial domain is described by a thin film equation,incorporating the electric driving force and the surface diffusion.The nonlinear interfacial growth includes the amplitude equations and superposition of one-dimensional structures at regular orientations.The pattern selection is driven by the subcritical instability mechanism in which the relative thickness of the polymer film plays a critical role.     

Stability of conducting liquid flowing down an inclined plane at moderate Reynolds number in the presence of constant electromagnetic field

The stability of a conducting viscous film flowing down an inclined plane at moderate Reynolds number in the presence of electromagnetic field is investigated under induction-free approximation. Using momentum integral method a non-linear evolution equation for the development of the free surface is derived. The linear stability analysis of the evolution equation shows that the magnetic field stabilizes the flow whereas the electric field stabilizes or destabilizes the flow depending on its orientation with the flow. The weakly non-linear study reveals that both the supercritical stability and subcritical instability are possible for this type of thin film flow. The influence of magnetic field on the different zones is very significant, while the impact of electric field is very feeble in comparison.     

扩散型气柱界面R-M失稳中混合率的实验研究

Richtmyer-Meshkov(R-M)不稳定性普遍存在于众多工程问题中,激波管实验是研究R-M失稳问题的主要手段.高精度的平面激光诱导荧光(planar laser-induced fluorescence,PLIF)技术具有分子量级的示踪能力,可获得界面气体浓度(摩尔分数)分布,为研究界面失稳混合问题提供了有力工具.在弱激波(Ma=1.25)冲击扩散型气柱界面实验中,采用PLIF技术对R-M失稳引起的SF6-Air界面混合问题进行了研究.通过改变椭圆形初始界面的长短轴比,得到了3种扩散型初始界面失稳演化过程中气体摩尔分数,观察到了斜压机制下界面的简单拉伸、二次不稳定性、挤压射流等现象.利用浓度分布进一步得到了界面的瞬时混合率,通过瞬时混合率、界面整体平均混合率以及混合率的概率密度分布,分析了界面在不同演化阶段的界面混合特征,初步讨论了界面失稳混合的机制.演化初期,界面在斜压涡的作用下发生拉伸卷曲,通过增大浓度梯度来促进界面的混合.当演化进一步发展,二次不稳定性出现后,界面通过小尺度对流的方式达到湍流混合状态,而浓度梯度驱使的分子间混合逐渐减弱.由浓度梯度引起的扩散与由二次不稳定性引起的对流存在着"竞争"关系,二者共同主导了界面的混合.     

圆柱形汇聚激波诱导 Richtmyer-Meshkov不稳定的 SPH 模拟

汇聚激波诱导不同物质界面的Richtmyer-Meshkov(RM)不稳定现象在惯性约束核聚变领域有重要的学术意义和工程背景.基于网格离散的宏观流体力学方法由于数值扩散问题往往需要高阶精度算法才能准确追踪界面演化,且对大变形和破碎合并等复杂界面追踪也极为困难.光滑粒子流体动力学(smoothed particlehydrodynamics,SPH)方法采用纯拉格朗日算法,可以有效克服上述难点.但经典SPH算法需采用人工黏性处理强间断,在激波间断处往往会出现严重的非物理振荡,对于涉及强冲击不稳定性问题,很难达到理想的模拟效果.本文采用基于HLL黎曼求解器的SPH算法,实现了对强激波和大密度比物质界面的有效分辨和追踪.一维数值校核证明了代码的可靠性、健壮性,并进一步模拟了二维圆柱形汇聚冲击波冲击四边形轻/重气界面诱导的RM不稳定性问题,与已有实验结果进行了对比,发现模拟结果与实验结果吻合.通过分析界面演化过程中的密度及压力变化,发现本文所采用的方法可准确地追踪激波与界面作用的复杂界面和波系演化规律.研究结果为进一步理解和解释汇聚冲击条件下的RM不稳定性机理奠定了基础.      

Bifurcation mechanism of interfacial electrohydrodynamic gravity-capillary waves near the minimum phase speed under a horizontal electric field

We investigate the evolution of interfacial gravity-capillary waves propagating along the interface between two dielectric fluids under the action of a horizontal electric field. There is a uniform background flow in each layer, and the relative motion tends to induce Kelvin–Helmholtz(KH) instability. The combined effects of gravity, surface tension and electrically induced forces are all taken into account. Under the short-wave assumption, the expansion and truncation method of Dirichlet-Neumann(DN) operators is applied to derive a reduced dynamical model. When KH instability is suppressed linearly by a considerably large electric field, our numerical results reveal that in certain regions of parameter space, nonlinear symmetric traveling wave solutions can be found near the minimum phase speed. Additionally, the detailed bifurcation structures are presented together with typical wave profiles.     

柱形汇聚几何中内爆驱动金属界面不稳定性

采用自研的高保真度爆轰与冲击动力学程序，对柱形汇聚几何中内爆驱动金属材料界面不稳定性的动力学行为，进行了数值模拟研究。结果表明，首次冲击后至约12 μs，界面发展以RM（Richtmyer-Meshkov）不稳定性为主；12 μs后至冲击波聚心反弹加载前，界面聚心运动处于加速减速状态，界面发展由RT (Rayleigh-Taylor)不稳定性主导；冲击波聚心反弹加载后，界面发展又由RM不稳定性主导。另外，还研究了初始条件（初始振幅、初始波长、钢壳初始厚度和几何构型）对柱形内爆驱动金属材料界面不稳定性的影响。结果显示:初始振幅较大时振幅增长也较大；初始波长较小（模数较大）时振幅增长较小，而且存在一个截止波长；钢壳厚度会抑制扰动增长，也存在一个截止厚度；几何汇聚效应会使扰动增长速度更快。