Development of the Raleigh-Tailor instability in inhomogeneous magnetic gas-dynamic flows

Results of the numerical simulation of the development of the hydrodynamic Raleigh-Tailor instability, which occurs when a high-current discharge is moving in a rail accelerator and in MHD flows with a current layer, are presented. It is established that the development of the instability under the conditions of the flow of cold gas around the discharge region leads to the formation of an irregular vortex structure of the magnetic gas-dynamic flow in the channel.     

On the Generation of Mean Flows by the Interaction of Görtler Vortices and Tollmien-Schlichting Waves in Curved Channel Flows

There are many fluid flows where the onset of transition can be caused by different instability mechanisms which compete in the nonlinear regime. Here the interaction of a centrifugal instability mechanism with the viscous mechanism which causes Tollmien-Schlichting waves is discussed. The interaction between these modes can be strong enough to drive the mean state; here the interaction is investigated in the context of curved channel flows so as to avoid difficulties associated with boundary layer growth. Essentially it is found that the mean state adjusts itself so that any modes present are neutrally stable even at finite amplitude. In the first instance the mean state driven by a vortex of short wavelength in the absence of a Tollmien-Schlichting wave is considered. It is shown that for a given channel curvature and vortex wavelength there is an upper limit to the mass flow rate which the channel can support as the pressure gradient is increased. When Tollmien-Schlichting waves are present then the nonlinear differential equation to determine the mean state is modified. At sufficiently high Tollmien-Schlichting amplitudes it is found that the vortex flows are destroyed, but there is a range of amplitudes where a fully nonlinear mixed vortex-wave state exists and indeed drives a mean state having little similarity with the flow which occurs without the instability modes. The vortex and Tollmien-Schlichting wave structure in the nonlinear regime has viscous wall layers and internal shear layers; the thickness of the internal layers is found to be a function of the Tollmien-Schlichting wave amplitude.     

纤维悬浮槽流空间模式稳定性分析

采用扰动的空间发展模式而非通常的时间发展模式,对含有悬浮纤维的槽流进行了线性稳定性分析。建立了适用于纤维悬浮流的稳定性方程并针对较大范围的流动Re数及扰动波角频率进行了数值求解。计算结果表明,纤维轴向抗拉伸力与流体惯性力之比H可以反映纤维对流动稳定性的影响。H增大使临界Re数升高,对应的扰动波数减小,扰动空间衰减率增加,扰动速度幅值的峰值降低,不稳定扰动区域缩小,长波扰动所受影响相对较大。纤维的存在抑制了流场的失稳。     

Weak Nonlinear Long Waves in Channel Flow with Internal Dissipation

This article concerns the evolution of long waves ( O (ε^−1/2) wavelength) of small [ O (ε)] amplitude in channel flow with internal dissipation. We use multiple scale expansions to derive a generalized Kuramoto–Sivashinsky (GKS) equation that governs the dominant asymptotic solution in the limit of small disturbances and marginal linear instability. We compare this solution with numerical integrations of the full quasilinear system, and show that the error is consistent with an asymptotic solution to ε^3/2 over a time interval of order ε^−3/2.     

明渠层流失稳与沙纹成因机理研究

动床水流中,泥沙起动之后,往往要出现沙纹,沙纹成因各家的解释不一。白玉川,罗纪生的观点是:沙纹的尺度较小,主要是由于明渠层流不稳定性波或床面近壁流层中小尺度拟序结构发展演化所致。当床面边界附近扰动波或拟序结构以及水流自身所产生的床面底部切应力大于Shields切应力后,床面即产生响应,形成沙纹;如果扰动所产生的扰动切应力频率接近床面泥沙固有频率,则产生与泥沙颗粒的共振,这种现象也称之为“泥沙的检波性质”,此时床面发生最大响应,沙纹发展速度也最快。     

Asymptotic calculation of inviscidly absolutely unstable modes of the compressible boundary layer on a rotating disk

In this work a long-wavelength asymptotic approach is used to analyze the region of absolute instability in the compressible rotating disk boundary layer flow. Theoretically determined values of branch points for the occurrence of absolute instability in the compressible flow are shown to match onto the ones which are obtained via a numerical solution of the linear inviscid compressible Rayleigh equations.     

The influence of electric fields and surface tension on Kelvin–Helmholtz instability in two-dimensional jets

We consider nonlinear aspects of the flow of an inviscid two-dimensional jet into a second immiscible fluid of different density and unbounded extent. Velocity jumps are supported at the interface, and the flow is susceptible to the Kelvin–Helmholtz instability. We investigate theoretically the effects of horizontal electric fields and surface tension on the nonlinear evolution of the jet. This is accomplished by developing a long-wave matched asymptotic analysis that incorporates the influence of the outer regions on the dynamics of the jet. The result is a coupled system of long-wave nonlinear, nonlocal evolution equations governing the interfacial amplitude and corresponding horizontal velocity, for symmetric interfacial deformations. The theory allows for amplitudes that scale with the undisturbed jet thickness and is therefore capable of predicting singular events such as jet pinching. In the absence of surface tension, a sufficiently strong electric field completely stabilizes (linearly) the Kelvin–Helmholtz instability at all wavelengths by the introduction of a dispersive regularization of a nonlocal origin. The dispersion relation has the same functional form as the destabilizing Kelvin–Helmholtz terms, but is of a different sign. If the electric field is weak or absent, then surface tension is included to regularize Kelvin–Helmholtz instability and to provide a well-posed nonlinear problem. We address the nonlinear problems numerically using spectral methods and establish two distinct dynamical behaviors. In cases where the linear theory predicts dispersive regularization (whether surface tension is present or not), then relatively large initial conditions induce a nonlinear flow that is oscillatory in time (in fact quasi-periodic) with a basic oscillation predicted well by linear theory and a second nonlinearly induced lower frequency that is responsible for quasi-periodic modulations of the spatio-temporal dynamics. If the parameters are chosen so that the linear theory predicts a band of long unstable waves (surface tension now ensures that short waves are dispersively regularized), then the flow generically evolves to a finite-time rupture singularity. This has been established numerically for rather general initial conditions.     

The influence of electric fields and surface tension on Kelvin–Helmholtz instability in two-dimensional jets

We consider nonlinear aspects of the flow of an inviscid two-dimensional jet into a second immiscible fluid of different density and unbounded extent. Velocity jumps are supported at the interface, and the flow is susceptible to the Kelvin–Helmholtz instability. We investigate theoretically the effects of horizontal electric fields and surface tension on the nonlinear evolution of the jet. This is accomplished by developing a long-wave matched asymptotic analysis that incorporates the influence of the outer regions on the dynamics of the jet. The result is a coupled system of long-wave nonlinear, nonlocal evolution equations governing the interfacial amplitude and corresponding horizontal velocity, for symmetric interfacial deformations. The theory allows for amplitudes that scale with the undisturbed jet thickness and is therefore capable of predicting singular events such as jet pinching. In the absence of surface tension, a sufficiently strong electric field completely stabilizes (linearly) the Kelvin–Helmholtz instability at all wavelengths by the introduction of a dispersive regularization of a nonlocal origin. The dispersion relation has the same functional form as the destabilizing Kelvin–Helmholtz terms, but is of a different sign. If the electric field is weak or absent, then surface tension is included to regularize Kelvin–Helmholtz instability and to provide a well-posed nonlinear problem. We address the nonlinear problems numerically using spectral methods and establish two distinct dynamical behaviors. In cases where the linear theory predicts dispersive regularization (whether surface tension is present or not), then relatively large initial conditions induce a nonlinear flow that is oscillatory in time (in fact quasi-periodic) with a basic oscillation predicted well by linear theory and a second nonlinearly induced lower frequency that is responsible for quasi-periodic modulations of the spatio-temporal dynamics. If the parameters are chosen so that the linear theory predicts a band of long unstable waves (surface tension now ensures that short waves are dispersively regularized), then the flow generically evolves to a finite-time rupture singularity. This has been established numerically for rather general initial conditions.     

An Analytic Approach for Calculating Absolutely Unstable Inviscid Modes of the Boundary Layer on a Rotating Disk

An analytical treatment of inviscidly absolutely unstable modes is pursued using the long-wavelength asymptotic approach. It is shown using the inviscid Rayleigh scalings in conjunction with the linear critical layer theory that the rotating-disk boundary layer flow undergoes a region of absolute instability for some small azimuthal wave numbers. The analytically calculated branch points for the absolute instability are found to be in good agreement with those obtained via a numerical solution of the inviscid Rayleigh equation.     

Nonlinear stability characterization of thin Newtonian film flows traveling down on a vertical moving plate

The paper presents both the linear and nonlinear stability theories for the characterization of thin Newtonian film flows traveling down along a vertical moving plate. The linear model is first developed to characterize the flow behavior. After showing the inadequacy of the linear model in representing certain flow characteristics, the nonlinear kinematics model is then developed to represent the system. The long-wave perturbation method is employed to derive the generalized kinematic equations with free film surface condition. The linear model is solved by using the normal mode method for three different, namely, the quiescent, up-moving and down-moving, moving conditions. Subsequently, the elaborated nonlinear film flow model is solved by the method of multiple scales. The modeling results clearly indicate that both subcritical instability and supercritical stability conditions are possible to occur in the film flow system. The effect of the down-moving motion of the vertical plate tends to enhance the stability of the film flow.     

Effects of roughness on nonlinear stationary vortices in rotating disk flows

Primary instability of rotating disk boundary layer flow over a rough surface for stationary modes was investigated by using the weakly nonlinear theory where the Reynolds number R is close to its critical value R_c as determined by linear theory. Both the single mode case, where the wave vector K equals its critical K_c at the onset of stationary primary instability, and the bimodal case, where the wave vectors K_n (n = 1, 2) are close to K_c for the primary instability of the flow, are considered. The analysis leads to stable solutions for particular roughness forms and magnitude, and particular wave vectors ˜K_n (n = 1, 2) of the surface roughness.     

The stability of Couette flow in a toroidal magnetic field

The stability of the hydromagnetic Couette flow is investigated when a constant current is applied along the axis of the cylinders. It is shown that if the resulting toroidal magnetic field depends only on this current, no linear instability to axisymmetric disturbances is possible.     

对冲圆射流生成径向扩张液膜研究

采用两股互相冲击的圆射流可以形成环形的液体薄膜,液膜在径向扩展到一定的临界半径距离会破碎．数值模拟了液膜在周围气体中形成和破碎的非定常过程．考虑了液体和气体都是不可压缩Newton流体的轴对称问题．液体和气体的界面采用Level set函数来跟踪,Navier-Stokes 控制方程和物理边界条件采用有限差分格式离散求解．计算结果给出了环形液体薄膜形成并在其环形边缘处破碎,并缓慢运动的过程．液膜的厚度随着液膜在轴向的扩展会逐渐变薄,因此定义的局部Weber数会在径向逐渐减小,这里的局部Weber数定义为ρu2h/σ,其中ρ和σ分别为液体的密度和界面的张力,u和h分别为在径向某个位置的液膜的平均径向速度和半液膜厚度．数值结果表明就像实验中所观察到的那样,液膜径向扩展的过程的确会在局部Weber数趋向于1的时候终结而停止扩张．根据空间-时间线性稳定性理论,液膜的破碎最初是由正弦模式在临界局部Weber数Wec=1引起的,在临界局部Weber数小于1时会发生绝对不稳定性．在线性理论中另一个独立的模式,所谓的余弦模式,则增长比正弦模式要慢,从而会推测到正弦模式主导破碎的结论．然而,这里的数值结果却表明,余弦模式在界面波的非线性发展阶段实质的超越了正弦模式的增长,并对液膜的最终阶段的破碎起主导作用．这验证了线性理论只能够对触发时扰动波的性质进行预测,而对失稳后情况和结果的预测则不一定正确．     

Absolute and Convective Instabilities of Semi-bounded Spatially Developing Flows

We analyse the absolute and convective instabilities of, and spatially amplifying waves in, semi-bounded spatially developing flows and media by applying the Laplace transform in time to the corresponding initial-value linear stability problem and treating the resulting boundary-value problem on ?⁺ for a vector equation as a dynamical system. The analysis is an extension of our recently developed linear stability theory for spatially developing open flows and media with algebraically decaying tails and for fronts to flows in a semi-infinite domain. We derive the global normal-mode dispersion relations for different domains of frequency and treat absolute instability, convectively unstable wave packets and signalling. It is shown that when the limit state at infinity, i.e. the associated uniform state, is stable, the inhomogeneous flow is either stable or absolutely unstable. The inhomogeneous flow is absolutely stable but convectively unstable if and only if the flow is globally stable and the associated uniform state is convectively unstable. In such a case signalling in the inhomogeneous flow is identical with signalling in the associated uniform state.     

张量封闭模型与三维取向分布对纤维悬浮槽流稳定性的影响

应用3种不同的纤维方向张量封闭模型,数值模拟了纤维悬浮槽流的流动稳定性问题,从而研究封闭模型和纤维的三维取向分布对稳定性分析的影响．结果发现,采用3种不同封闭模型所得到的流动稳定特性与纤维参数之间的关系是相同的,但采用三维混合封闭模型时,由于纤维的取向与流向的偏差程度较大,所以纤维对流动的不稳定性具有最强的抑制作用．而采用二维混合封闭模型时,由于纤维在平面取向条件下,其轴线整体上趋于呈流向排列,使得对流体的作用削弱,导致纤维对流动不稳定性抑制的作用最弱．     

On Disturbances to a Stratified Mixing Layer

Using a nonlinear critical layer analysis, we examine the behavior of disturbances to the Holmboe model of a stratified shear layer for Richardson numbers 0相似文献   

Non-linear stability characterization of the thin micropolar liquid film flowing down the inner surface of a rotating vertical cylinder

Linear and non-linear stability analysis for characterization of micropolar film flowing down the inner surface of a rotating infinite vertical cylinder is given. A generalized non-linear kinematic model is derived to represent the physical system and is solved by the long wave perturbation method in the following procedure. First, the normal mode method is used to characterize the linear behaviors. Then, an elaborated non-linear film flow model is solved by using the method of multiple scales to characterize flow behaviors at various states of sub-critical stability, sub-critical instability, supercritical stability, and supercritical explosion. The modeling results indicate that by increasing the rotation speed, Ω, and the radius of cylinder, R, the film flow will generally stabilize the flow system.     

Minimax principle for magnetohydrodynamic channel flow

A variational principle, applicable to a wide range of linear non-self-adjoint problems, is used to derive maximum and minimum principles for a magnetohydrodynamic channel flow. Truncated Fourier series for the velocity and magnetic field are found which give reasonable approximations at small Hartmann numbers.     

Stability of a deflagration front in compressible media

A two-dimensional linear analysis of the planar flame front stability for a compressible fluid is presented. The analysis shows that there are two types of perturbations. The first type, corresponding to waves in incompressible media, has already been studied by Landau. It predicts absolute instability of the flame front. The second type of perturbations is due to fluid compressibility and the dependence on upstream flow parameters of the flame front velocity. Three different regimes for these perturbations are possible: stable, acoustically unstable, and absolutely unstable. The instability results in a pronounced pressure wave generation.A one-dimensional analysis of the interaction of the flame front with flow boundaries is performed. Under some circumstances, this interaction is shown to cause exponential growth of the perturbations.     

Numerical simulation of spinning detonation in circular section channels

Numerical simulation of three-dimensional structures of gas detonation in circular section channels that emerge due to the instability when the one-dimensional flow is initiated by energy supply at the closed end of the channel is performed. It is found that in channels with a large diameter, an irregular three-dimensional cellular detonation structure is formed. Furthermore, it is found that in channels with a small diameter circular section, the initially plane detonation wave is spontaneously transformed into a spinning detonation wave, while passing through four phases. A critical value of the channel diameter that divides the regimes with the three-dimensional cellular detonation and spinning detonation is determined. The stability of the spinning detonation wave under perturbations occurring when the wave passes into a channel with a greater (a smaller) diameter is investigated. It is found that the spin is preserved if the diameter of the next channel (into which the wave passes) is smaller (respectively, greater) than a certain critical value. The computations were performed on the Lomonosov supercomputer using from 0.1 to 10 billions of computational cells. All the computations of the cellular and spinning detonation were performed in the whole long three-dimensional channel (up to 1 m long) rather than only in its part containing the detonation wave; this made it possible to adequately simulate and investigate the features of the transformation of the detonation structure in the process of its propagation.