可压缩流体的Rayleigh-Taylor和Kelvin-Helmholtz不稳定性

在状态方程为压力是密度的任意单值函数形式情况下，运用小扰动分析和奇异摄动法，给出了流体微扰方程渐近解和界面不稳定性的色散关系. 分析表明：对Rayleigh-Taylor不稳定性，在重力场作用下流体可压缩性形成的密度分布是致稳因素；而扰动流体的膨胀收缩效应助长不稳定性的发展；上层重流体的可压缩性是稳定因素，下层轻流体可压缩性是失稳因素. 而对Kelvin-Helmholtz不稳定性，流体可压缩性助长扰动的发展，是不稳定因素.     

基于FTM方法的二维Kelvin-Helmholtz不稳定性数值模拟

使用界面跟踪法FTM（Front Tracking Method）对二维不混溶、不可压缩流体的K-H（Kelvin-Helmholtz）不稳定性进行数值模拟。研究表明,速度梯度层越厚,界面在水平分量中移动越快,卷起越少;初始水平速度差越大,界面卷起越多,内扰动增长速度越快,K-H不稳定性的特征形式更加明显;此外,在Neumann边界条件（即无滑移边界条件）下界面的扰动发展得比Dirichlet边界条件（即对称边界条件）下的扰动快。由于Dirichlet边界中的边界层,在开始时刻涡量扩展到两侧,影响了K-H不稳定性的生长速率;而在Neumann边界条件下涡量由于初始水平速度差,在界面中心聚集。最后,研究了不同边界条件下各种理查德森数对K-H不稳定性的影响。     

Flutter instability of a fluid-conveying fluid-immersed pipe affixed to a rigid body

A set of simplified boundary conditions for a flexible beam connected to a rigid body at one end and free at the other end is presented and applied to the case of a fluid-conveying, fluid-immersed pipe. These boundary conditions represent an analytically tractable approximation to those of a submersible which uses a combination of jet action and flutter instability induced tail motion to produce thrust. The boundary conditions are made non-dimensional, and the effect of the non-dimensional mass of the rigid body on system stability is assessed. The neutral stability of this system is determined within a two-parameter space consisting of the velocity of the fluid within the tail, and the forward speed of the submersible. Equations in the literature, derived using slender-body theory, were used to compute the sign of the thrust produced by the tail and the tail's Froude efficiency for the neutrally stable waveforms of the beam.     

GENERAL SECOND ORDER FLUID FLOW IN A PIPE

ＧＥＮＥＲＡＬＳＥＣＯＮＤＯＲＤＥＲＦＬＵＩＤＦＬＯＷＩＮＡＰＩＰＥＨｅＧｕａｎｇｙｕ（何光渝）（ＤｅｐａｒｔｍｅｎｔｏｆＰｅｔｒｏｌｅｕｍＥｎｇｉｎｅｅｒｉｎｇ，Ｘｉ＇ａｎＰｅｔｒｏｌｅｕｍＩｎｓｔｉｔｕｔｅ，Ｘｉ＇ａｎ７１００６１，Ｐ．Ｒ．Ｃｈｉ...     

Impracticality of MHD Convection in a Porous Medium

The subject of the title is discussed.     

Numerical solutions of pulsating flow and heat transfer characteristics in a pipe

Numerical studies are made of flow and heat transfer characteristics of a pulsating flow in a pipe. Complete time-dependent laminar boundary-layer equations are solved numerically over broad ranges of the parameter spaces, i.e., the frequency parameter β and the amplitude of oscillation A. Recently developed numerical solution procedures for unsteady boundary-layer equations are utilized. The capabilities of the present numerical model are satisfactorily tested by comparing the instantaenous axial velocities with the existing data in various parameters. The time-mean axial velocity profiles are substantially unaffected by the changes in β and A. For high frequencies, the prominent effect of pulsations is felt principally in a thin layer near the solid wall. Skin friction is generally greateer than that of a steady flow. The influence of oscillation on skin friction is appreciable both in terms of magnitude and phase relation. Numerical results for temperature are analyzed to reveal significant heat transfer characteristics. In the downstream fully established region, the Nusselt number either increases or decreases over the steady-flow value, depending on the frequency parameter, although the deviations from the steady values are rather small in magnitude for the parameter ranges computed. The Nusselt number trend is amplified as A increases and when the Prandtl number is low below unity. These heat transfer characteristics are qualitatively consistent with previous theoretical predictions.     

圆管湍流减阻电磁力控制的直接数值模拟

利用直接数值模拟研究了圆管湍流减阻的电磁力控制问题. 在圆管表面布置一层交替排 列的条形电极和磁极,施加电压后在导电的流体介质中产生周向随时间变化的电磁力,通过 改变电磁力的大小和变化周期等参数可以使阻力大幅减少. 研究减阻率和控制参数的 关系,并研究了电磁力减阻的机理. 电磁力控制存在最优周期,阻力最大降 低24.2{\%}; 电磁力控制抑制了流向和法向速度的脉动,而周向脉动速度随振荡周 期的增大而增大.     

Kelvin-Helmholtz instability of a compressible plasma in a magnetic field

The Kelvin-Helmholtz instability of a compressible plasma in a uniform magnetic field with respect to disturbances propagating along the flow is considered. First, the case with the magnetic field parallel to the direction of streaming is considered. The result given by Sen [4] that the compressibility effects destabilize an otherwise neutrally stable system even in the hydrodynamic limit is apparently erroneous. Re-examination of the dispersion relation in the limit of small compressibility effects shows that the latter reduce the growth rate of an otherwise unstable disturbance. Attention is also drawn to errors in the calculations of Fejer [2] in the limit of small compressibility effects. Next, the case with the magnetic field transverse to the direction of streaming is considered. It is found that the transverse magnetic field does influence the stability of the system when the compressibility effects are present, contrary to the result given by Chandrasekhar [1] for the case of an incompressible plasma. However, interestingly enough, the compressibility effects are effectively reduced if a transverse magnetic field is present! It is further shown that the transverse magnetic field reduces the stability of the system.     

Generation of compressible modes in MHD turbulence

Astrophysical turbulence is magnetohydrodynamic (MHD) in nature. We discuss fundamental properties of MHD turbulence and in particular the generation of compressible MHD waves by Alfvénic turbulence and show that this process is inefficient. This allows us to study the evolution of different types of MHD perturbations separately. We describe how to separate MHD fluctuations into three distinct families: Alfvén, slow, and fast modes. We find that the degree of suppression of slow and fast modes production by Alfvénic turbulence depends on the strength of the mean field. We review the scaling relations of the modes in strong MHD turbulence. We show that Alfvén modes in compressible regime exhibit scalings and anisotropy similar to those in incompressible regime. Slow modes passively mimic Alfvén modes. However, fast modes exhibit isotropy and a scaling similar to that of acoustic turbulence both in high and low plasmas. We show that our findings entail important consequences for star formation theories, cosmic ray propagation, dust dynamics, and gamma ray bursts. We anticipate many more applications of the new insight to MHD turbulence and expect more revisions of the existing paradigms of astrophysical processes as the field matures. PACS 47.65.+a; 52.30.Cv; 52.35.Ra; 95.30.Qd     

进动圆筒内液体流动不稳定的非线性特征

在建立进动充液圆筒内液体偏差流动方程的基础上，结合液体惯性波和轴向二次流动线性解，通过对定常二次流动的线性稳定性分析，提出了函数空间表达的流动不稳定性非线性分岔分析方程. 对非惯性坐标系下液体流动的Navier-Stokes方程进行了数值求解，并对惯性波发生破裂(实验提供的3种主模态下得出的共振破裂现象)时的压力时间序列进行分析，得出了液体流动不稳定的基本非线性特征.     

Transition from Planar to Whirling Oscillations in a Certain Nonlinear System

A single-mass two-degrees-of-freedom system is considered, witha radially oriented nonlinear restoring force. The latter is smooth andbecomes infinite at a certain value of a radial displacement. Stabilityanalysis is made for planar oscillation, or motion along a givendirection. As long as this motion is periodic, the nonlinearity in therestoring force provides a periodic parametric excitation in thetransverse direction. The linearized stability analysis is reduced tostudy of the Mathieu equation for the (infinitesimal) motions in thetransverse direction. For the case of free oscillations in the givendirection an exact solution is obtained, since a specific analyticalform is used for the (strongly nonlinear) restoring force, which permitsexplicit integration of the equation of motion. Stability of the planarmotion in this case is shown to be very sensitive to even slightdeviations from polar symmetry in the restoring force (as well as to theamplitude of oscillations in the given direction). Numerical integrationof the original equations of motion shows the resulting motion to be awhirling type indeed in case of the transversal instability. For thecase of a sinusoidal forcing in the given direction solution for the(periodic) response is obtained by Krylov–Bogoliubov averaging. Thisresults in the transmitted Ince–Strutt chart – namely, stabilitychart for transverse direction on the amplitude-frequency plane of theexcitation in the original direction.     

Boundary element method solution of MHD flow in a rectangular duct

The magnetohydrodynamic (MHD) flow of an incompressible, viscous, electrically conducting fluid in a rectangular duct with an external magnetic field applied transverse to the flow has been investigated. The walls parallel to the applied magnetic field are conducting while the other two walls which are perpendicular to the field are insulators. The boundary element method (BEM) with constant elements has been used to cast the problem into the form of an integral equation over the boundary and to obtain a system of algebraic equations for the boundary unknown values only. The solution of this integral equation presents no problem as encountered in the solution of the singular integral equations for interior methods. Computations have been carried out for several values of the Hartmann number (1 ? M ? 10). It is found that as M increases, boundary layers are formed close to the insulated boundaries for both the velocity and the induced magnetic field and in the central part their behaviours are uniform. Selected graphs are given showing the behaviours of the velocity and the induced magnetic field.     

On the Electric Potential Induced by a Homogeneous Magnetic Field in a Turbulent Pipe Flow

In this paper we study a turbulent pipe flow of a weakly electrical conducting fluid subjected to a homogeneous magnetic field which is applied perpendicular to the flow. This configuration forms the basis of a so-called electromagnetic induction flow meter. When the Hartmann number is small so that modification of flow by the Lorenz force can be neglected, the influence of the magnetic field results only in a spatially and temporally varying electric potential. The magnitude of the potential difference across the pipe is then proportional to the flow rate and this constitutes the principle of the flow meter. In this study the flow and electric potential are computed with help of a numerical flow simulation called Large-Eddy Simulation (LES) to which we have added an equation for the electrical potential. The results of the LES have been compared with experiments in which the electric potential is measured as a function of time at several positions on the circumference of the pipe. Both the experimental and numerical results for the mean potential at the pipe wall agree very well with an exact solution that can be obtained in this particular case of a homogeneous magnetic field. Furthermore, it is found that fluctuations in the electric potential due to the turbulence, are small compared to the velocity fluctuations. Based on the results we conclude that electrical-magnetic effects in pipe flow can be accurately computed with LES.     

Instability of the interface in two-layer flows with large viscosity contrast at small Reynolds numbers

The Kelvin–Helmholtz instability is believed to be the dominant instability mechanism for free shear flows at large Reynolds numbers. At small Reynolds numbers, a new instability mode is identified when the temporal instability of parallel viscous two fluid mixing layers is extended to current-fluid mud systems by considering a composite error function velocity profile. The new mode is caused by the large viscosity difference between the two fluids. This interfacial mode exists when the fluid mud boundary layer is sufficiently thin. Its performance is different from that of the Kelvin–Helmholtz mode. This mode has not yet been reported for interface instability problems with large viscosity contrasts.These results are essential for further stability analysis of flows relevant to the breaking up of this type of interface.     

NONLINEAR SATURATION OF BAROCLINIC INSTABILITY IN THE GENERALIZED PHILLIPS MODEL (Ⅰ) -THE UPPER BOUND ON THE EVOLUTION OF DISTURBANCE TO THE NONLINEARLY UNSTABLE BASIC FLOW

On the basis of the nonlinear stability theorem in the context of Arnol’s second theorem for the generalized Phillips model, nonlinear saturation of baroclinic instability in the generalized Phillips model is investigated. By choosing appropriate artificial stable basic flows, the upper bounds on the disturbance energy and potential enstrophy to the nonlinearly unstable basic flow in the generalized Phillips model are obtained, which are analytic completely and without the limitation of infinitesimal initial disturbance.     

Nonlinear instability of wavetrain under influences of shear current with varying vorticity and air pressure

Under the assumption of weak shear current with varying vorticity in water and weak air pressure the Zakharov theory is extended to include the effects of vorticity and air pressure on the modulation of water waves. This new equation is used to examine the influence of current and wind on the Benjamin-Feir sideband instability and long-time evolution of wavetrain. As strength of the current increases the bandwidth is found broadened, and the maximum growth rate of sidebands decreased. Periodic solution of sidebands in the presence of current is indicated, which means that shear current does not affect the downshift of wave spectrum peak. Energy input by imposing the air pressure leads to the enhancement of the lower sideband, which is in agreement with the finding of Hara and Mei (1991). The project supported by the National Natural Science Foundation of China     

Three-phase upward flow in a vertical pipe

In this paper we develop an approach to design a three-phase, gas–solid–liquid flow system that transports pneumatically scarified solid particles, including sticky ones, through a vertical pipe. The proposed system permits the introduction and maintenance of a liquid film that coats the pipe’s inner wall and acts as a lubricant that ensures sticky particles continue to move upward without permanently adhering to the pipe wall. The system’s operating conditions fall within the boundaries of the annular dispersed region on a typical flow pattern map of vertical flow of a gas–liquid mixture. High gas superficial velocities combined with low liquid superficial velocities characterize such a region. A combination of a modified one-dimensional, two-fluid annular dispersed flow model and a one-dimensional pneumatic conveying model is shown to describe this transport process satisfactorily. Solution of the combined models produces all the necessary design parameters including power requirements and superficial velocities of the two-fluid media needed to transport a given amount of solid particles. Results of model calculations are compared with rare three-phase flow data obtained prior to the development of the present model, by an independent experimental team that used the physical conditions of the present approach. Reasonable agreement justifies the use of the combined model for engineering design purposes.     

Nonlinear waves in a liquid film on a near-horizontal surface

Nonlinear waves in a liquid film on a slightly inclined rigid plane are studied. A mathematical model is reduced to a system of two evolutionary equations for the layer thickness and the local fluid mass flow. In addition to viscous forces, gravity, and surface tension, the pressure difference over the layer thickness, induced by the gravity force projection on the normal to the underlying surface, is also taken into account. Spatially periodic solutions developing with time from small initial disturbances into regular nonlinear waves are considered. A spectral representation of the solution, the Galerkin method with respect to the uniform coordinate, and subsequent numerical calculation of the corresponding dynamic system on large time intervals are employed. Different variants in the space of the three governing parameters are calculated and some basic mechanisms of nonlinear dynamics of the two-dimensional waves are detected. The calculation results are compared with the existing experimental data. It is shown that the theoretical conclusions can be used to interpret and predict experiments.