Steady vortices in plasmas and geophysical flows

A number of two-dimensional fluid models in geophysical fluid dynamics and plasma physics are examined to find out whether they have steady and localized monopole vortex solutions. A simple and general method that consists of two steps is used. First the dispersion relation is calculated, to find all possible values of the phase velocity of the linear waves. Then an integral relation that determines the center-of-mass velocity of localized structures must be found. The existence condition is that this velocity should be outside the region of linear phase velocities. After a presentation of the method, previous work on the plasma drift wave model and the shallow-water equations is reviewed. In both cases it is found that the center-of-mass velocity is larger than the maximum phase velocity of the linear waves if the amplitude is large enough, and steady localized vortices can therefore exist. New results are then obtained for a number of two-field models. For the coupled ion acoustic-drift modes in plasmas, it is found that the center-of-mass velocity depends on the ratio between the parallel ion velocity component and the electrostatic potential in the vortex. If this ratio is large enough, the vortex can be steady. For the drift-Alfven mode the "center-of-charge" velocity is proportional to the ratio between the parallel current and the total charge in the vortex. It can therefore be steady if this ratio satisfies the appropriate conditions. For the quasigeostrophic two-layer equations, describing stratified flow on a rotating planet, it is found that the center-of-mass velocity is determined by the ratio between the baroclinic and the barotropic components in the vortex. If a baroclinic component with an appropriate sign is added to a barotropic vortex, it propagates faster than the barotropic Rossby waves, and can be steady. Finally, the existence conditions for a vortex in an external zonal flow are examined. It is found that the center-of-mass velocity acquires an additional westward contribution in an anticyclonic shear zone in the framework of the shallow-water equations, and also that an easterly jet south of this shear zone partly shields a vortex situated in the shear zone from the dispersive influence of the fast Rossby waves on the equatorward side.     

On the damping of unsteady flow by compliant boundaries

This paper concerns unsteady motion of an incompressible inviscid fluid near a flexible surface which, in responding to the surface pressure field, absorbs energy. The modification of the flow consequent on energy removal at the boundaries is examined. Energy absorption always occurs when the mechanical surface properties include an element of dissipation. But surface dissipation is not essential; surface waves have a similar property. Unsteady fluid induced forces excite surface waves which carry with them energy that must have originated in the flow. The question of how flow characteristics change as energy is gradually given up to the boundary is examined through a particular model problem from which it becomes evident that surface motion draws vorticity towards the surface. The model chosen is that of a rectilinear vortex adjacent to a weakly responding boundary. Surface motion induces a velocity perturbation which is shown to move the vortex towards the surface whenever the fluid gives energy to that surface.     

斜压二层模式框架下的准地转理论及其应用

准地转理论是短期天气预报的理论基础.针对中纬度大尺度运动得到的准地转运动方程组是准地转理论的核心.在斜压二层模式框架下,由准地转运动方程组得出的准地转位势倾向方程和准地转垂直运动方程,可以更为明确地解释高空和地面天气系统的移动和发展.结果显示:500 hPa高空槽的移动取决于本层的涡度平流,500 hPa高空槽的发展取决于本层的涡度平流和250与750h Pa层的微差温度平流;地面气旋系统的移动和发展取决于500hPa的正温度平流,250与750 hPa层的微差涡度平流.由实际天气过程对上述应用进行的分析表明,在斜压二层模式框架下得到的准地转位势倾向方程及准地转垂直运动方程,可以更好地确定中纬度天气尺度斜压发展系统的特征,并有助于对天气系统发展客观规律的了解和对数值预报产品内含的物理依据的认识.     

Relativistic vortex dynamics in axisymmetric stationary perfect fluid configuration

Relativistic formulation of Helmholtz’s vorticity transport equation is presented on the basis of Maxwell-like version of Euler’s equation of motion. Entangled characteristics associated with vorticity flux conservation in a vortex tube and in a stream tube are displayed on basis of Greenberg’s theory of spacelike congruence of vortex lines and \(1+1+(2)\) decomposition of the gradient of fluid’s 4-velocity. Vorticity flux surfaces are surfaces of revolution about the rotation axis and are rotating with fluid’s angular velocity due to gravitational isorotation in a stationary axisymmetric perfect fluid configuration. Fluid’s angular velocity, angular momentum per baryon, injection energy, and invariant rotational potential are constant on such vorticity flux surfaces. Gravitation causes distortion of coaxial cylindrical vorticity flux surfaces in the limit of post-Newtonian approximation. The rotation of the fluid with angular velocity relative to vorticity flux surfaces generates swirl which causes the stretching of material vortex lines being wrapped on vorticity flux surfaces. Fluid helicity which is conserved in the fluid’s rest frame does not remain conserved in a locally nonrotating frame because of the existence of swirl. Vortex lines are twist free in the absence of meridional circulations, but the twisting of spacetime due to dragging effect leads to the increase in vorticity flux in a vortex tube.     

Experimental investigation of primary and secondary features in high-mach-number shock-bubble interaction

Experiments to study the compression and unstable evolution of an isolated soap-film bubble containing helium, subjected to a strong planar shock wave (M=2.95) in ambient nitrogen, have been performed in a vertical shock tube of square internal cross section using planar laser diagnostics. The early phase of the interaction process is dominated by the formation of a primary vortex ring due to the baroclinic source of vorticity deposited during the shock-bubble interaction, and the mass transfer from the body of the bubble to the vortex ring. The late time (long after shock interaction) study reveals the presence of a secondary baroclinic source of vorticity at high Mach number which is responsible for the formation of counterrotating secondary and tertiary vortex rings and the subsequent larger rate of elongation of the bubble.     

区域谱元正压大气模式台风移动数值模拟试验

针对地图投影坐标系下的正压原始方程组,将计算区域按三角形元进行分解,在三角形元内用三角形截断的勒让德多项式的积为插值函数对变量进行谱分解,发展出区域正压谱元大气模式.采用固定边界条件,以2006年5月15日08时500 hPa位势高度和风场为初值,在勒让德多项式最高阶数为3和7这两种情形下开展0601号台风"珍珠"移动的数值模拟试验.结果表明,数值模式模拟的风压场关系合理,数值模式的实现是成功的.     

Symmetry analysis and explicit solutions of the (3+1)-dimensional baroclinic potential vorticity equation

<正>This paper investigates an important high-dimensional model in the atmospheric and oceanic dynamics-(3+l)- dimensional nonlinear baroclinic potential vorticity equation by the classical Lie group method.Its symmetry algebra, symmetry group and group-invariant solutions are analysed.Otherwise,some exact explicit solutions are obtained from the corresponding(2+1)-dimensional equation,the inviscid barotropic nondivergent vorticy equation.To show the properties and characters of these solutions,some plots as well as their possible physical meanings of the atmospheric circulation are given.     

Vortex dynamics in different combustion regions of a cavity-based scramjet

A hybrid RANS/LES study of a cavity-based scramjet was performed and reasonable agreements were found between simulation results and experimental measurements. In the current case, the flame was stabilized by the subsonic cavity shear layer and propagated downstream into the supersonic flow. The vortex dynamic in the flow, mixing, and combustion regions was comparatively investigated. The averaged vorticity in the combustion regions was lower by 55% compared to the mixing region, primarily due to dilatation as a result of the heat release. Furthermore, the combustion zone was decomposed into four regions based on premixed/diffusion flame and subsonic/supersonic combustion. Then the vorticity and its transport in the four regions were compared. The averaged vorticity in the premixed combustion regions was only slightly larger than that in the diffusion combustion regions. However, the averaged heat release rate was nearly 3 times larger in the premixed regions, leading to higher contributions of dilatation and baroclinic torque in the premixed regions, with an overall weak positive impact on the vorticity generation. In the subsonic combustion regions, the vorticity was three times larger than that in the supersonic combustion regions, despite similar heat release rates on average. It could be explained by the relatively large magnitude of dilatation and baroclinic torque in the supersonic flow. Vortex stretching and dilatation were comparable in the supersonic flame but the former became two times larger than the latter in the subsonic flame. Moreover, the baroclinic torque had larger contributions than diffusion in the supersonic flame whereas the opposite trend was found in the subsonic flame. The results highlight that the subsonic combustion regions in the cavity shear layer and near the walls significantly contribute to the vortex dynamics and mixing process, in addition to flame stabilization.     

Type-N,shear-free,perfect-fluid spacetimes with a barotropic equation of state

We present the class of Petrov-typeN, shear-free, perfect-fluid solutions of Einstein's field equations in which the fluid satisfies a barotropic equation of statep=p(w) andw+p0. All solutions are stationary and possess a three- parameter, abelian group of local isometries which act simply transitively on timelike hypersurfaces. Furthermore, there exists one Killing vector parallel to the vorticity vector and another parallel to the four-velocity. This class of solutions is identified as part of a larger class present in the literature.     

Vorticity dynamics and sound generation in two-dimensional fluid flow

An approximate solution to the two-dimensional incompressible fluid equations is constructed by expanding the vorticity field in a series of derivatives of a Gaussian vortex. The expansion is used to analyze the motion of a corotating Gaussian vortex pair, and the spatial rotation frequency of the vortex pair is derived directly from the fluid vorticity equation. The resulting rotation frequency includes the effects of finite vortex core size and viscosity and reduces, in the appropriate limit, to the rotation frequency of the Kirchhoff point vortex theory. The expansion is then used in the low Mach number Lighthill equation to derive the far-field acoustic pressure generated by the Gaussian vortex pair. This pressure amplitude is compared with that of a previous fully numerical simulation in which the Reynolds number is large and the vortex core size is significant compared to the vortex separation. The present analytic result for the far-field acoustic pressure is shown to be substantially more accurate than previous theoretical predictions. The given example suggests that the vorticity expansion is a useful tool for the prediction of sound generated by a general distributed vorticity field.     

Vibration of a liquid with a free surface in a spinning spherical tank

The problem of forced fluid vibrations in a partially filled spinning spherical tank is solved numerically by using the finite element method. The governing equations include Coriolis acceleration and spatially homogeneous vorticity. An exponential instability is detected in the present simulation for fill ratios below 0·5 and centrifugal acceleration to thrust ratios less than 1·7. This fictitious instability appears in the model as a result of the homogeneous vortex assumption since the free slosh equations are neutrally stable in the Liapunov sense.     

Observation of vector solitons with hidden vorticity

This letter reports the first experimental observation, to our knowledge, of optical vector solitons composed of two incoherently coupled vortex components. We employ nematic liquid crystal to generate stable vector solitons with counterrotating vortices and hidden vorticity. In contrast, the solitons with explicit vorticity and corotating vortex components show azimuthal splitting.     

A high resolution spatially adaptive vortex method for separating flows. Part I: Two-dimensional domains

A grid-free high-resolution spatially-adaptive vortex method for two-dimensional incompressible flow in bounded domains is presented. The computational algorithm is based on operator splitting in which convection and diffusion are handled separately every time step. In the convection step, computational elements are convected with velocities obtained by fast approximations of the Biot–Savart superposition with second-order Runge–Kutta time integration scheme. Diffusion is performed using the smooth redistribution method that employs a Gaussian basis function for vorticity in the interior. Near solid walls, the core functions are modified to conserve circulation. The no-slip boundary condition is enforced by creating of a vortex sheet that is redistributed to neighboring elements using the redistribution method. The proposed method enables accurate and smooth recovery of the vorticity and does not require explicit use of vortex images or occasional re-meshing. Algorithms for reduction in computational cost by accurately removing elements in overcrowded regions and for spatial adaptivity that allows for variable core sizes and variable element spacing are presented. Computations of flow around an impulsively started cylinder for Reynolds number values of 1000, 3000, and 9500 are preformed to investigate various aspects of the proposed method.     

湿斜压热动力耦合强迫激发辐合增长和暴雨维持的一种机制

打破将散度演化分析局限于水平动力学框架内的传统认识思维,针对大气的斜压性与层结特征将对大气运动散度场演化产生显著影响的事实,从湿斜压原始方程出发,将大气的斜压热力作用通过热力学方程引入散度方程,在替换平衡近似（alternative balance）下,导出了显式包含湿非地转Q矢量与垂直风切变耦合强迫作用的新型散度方程.阐明了潮湿大气中非地转斜压强迫激发散度场演化的物理机制,为利用常规观测资料诊断散度场演化提供了一种有效的数学工具.通过对流层中低层气柱内总体辐合增长影响因子的分 关键词： 暴雨天气 散度演变 正压非平衡强迫 湿斜压热动力耦合强迫     

The GHP II Package with Applications

We present an advanced version of the Maple package GHP called GHPII. In it we provide a number of additional sophisticated tools to assist with problems formulated in the Geroch-Held-Penrose (ghp) formalism. The first part of this article discusses these new tools while in the second part we shall apply the ghp formalism, using the GHPII routines, to vacuum Petrov type D spacetimes and shear-free perfect fluids. We prove that for all shear-free perfect fluids with a barotropic equation of state, where two of the principal null directions are coplanar with the fluid four-velocity and vorticity then either the expansion or vorticity of the fluid must be zero.     

Geometric Properties of the Maxwell Set and a Vortex Filament Structure for Burgers Equation

The inviscid limit of the stochastic Burgers equation is discussed in terms of the level surfaces of the minimising Hamilton–Jacobi function, the classical mechanical caustic and the Maxwell set and their algebraic pre-images under the classical mechanical flow map. We examine the geometry of the Maxwell set in terms of the behaviour of the pre-Maxwell set, the pre-caustic and the pre-level surfaces. In particular, contrary to the ideas of Helmholtz and Lord Kelvin, we prove that even if initially the fluid flow is irrotational, in the inviscid limit, associated with the advent of the Maxwell set a non-zero vorticity vector forms in the fluid with vortex lines on the Maxwell set. This suggests that in quite general circumstances for small viscosity there is a vortex filament structure near the Maxwell set for both deterministic and stochastic Burgers equations.      

A nonlinear scenario for development of vortex layer instability in gravity field

A Hamiltonian version of contour dynamics is formulated for models of constant-vorticity plane flows with interfaces. The proposed approach is used as a framework for a nonlinear scenario for instability development. Localized vortex blobs are analyzed as structural elements of a strongly perturbed wall layer of a vorticity-carrying fluid with free boundary in gravity field. Gravity and vorticity effects on the geometry and velocity of vortex structures are examined. It is shown that compactly supported nonlinear solutions (compactons) are candidates for the role of particle-like vortex structures in models of flow breakdown. An analysis of the instability mechanism demonstrates the possibility of a self-similar collapse. It is found that the vortex shape stabilizes at the final stage of the collapse, while the vortex sheet strength on its boundary increases as (t ₀ ? t)^?1, where t ₀ is the collapse time.     

Jupiter's Great Red Spot and zonal winds as a self-consistent,one-layer,quasigeostrophic flow

We present the point of view that both the vortices and the east-west zonal winds of Jupiter are confined to the planet's shallow weather layer and that their dynamics is completely described by the weakly dissipated, weakly forced quasigeostrophic (QG) equation. The weather layer is the region just below the tropopause and contains the visible clouds. The forcing mimics the overshoot of fluid from an underlying convection zone. The late-time solutions of the weakly forced and dissipated QG equations appear to be a small subset of the unforced and undissipated equations and are robust attractors. We illustrate QG vortex dynamics and attempt to explain the important features of Jupiter's Great Red Spot and other vortices: their shapes, locations with respect to the extrema of the east-west winds, stagnation points, numbers as a function of latitude, mergers, break-ups, cloud morphologies, internal distributions of vorticity, and signs of rotation with respect to both the planet's rotation and the shear of their surrounding east-west winds. Initial-value calculations in which the weather layer starts at rest produce oscillatory east-west winds. Like the Jovian winds, the winds are east-west asymmetric and have Karman vortex streets located only at the west-going jets. From numerical calculations we present an empirically derived energy criterion that determines whether QG vortices survive in oscillatory zonal flows with nonzero potential vorticity gradients. We show that a recent proof that claims that all QG vortices decay when embedded in oscillatory zonal flows is too restrictive in its assumptions. We show that the asymmetries in the cloud morphologies and numbers of cyclones and anticyclones can be accounted for by a QG model of the Jovian atmosphere, and we compare the QG model with competing models.     

Counter rotating vortex pair structure in a reacting jet in crossflow

This paper analyzes the time averaged flow structure of a reacting jet in cross flow (RJICF), emphasizing the structure of the counter-rotating vortex pair (CVP) by using simultaneous tomographic particle image velocimetry (TPIV) and hydroxyl radical planar laser induced fluorescence (OH-PLIF). It was performed to determine the extent to which heat release, and the associated effects of gas expansion and baroclinic vorticity production, impact the structure of the CVP. These results show the clear presence of a CVP in the time averaged flow field, whose trajectory lies below the jet centerline on either side of the centerplane. Consistent with other measurements of high momentum flux ratio JICF in nonreacting flows, there is significant asymmetry in strength of the two vortex cores. The strength and structure of the CVP was quantified with vorticity and swirling strength (λ_ci), showing that some regions of the flow with high shear are not necessarily accompanied by large scale bulk flow rotation and vice-versa. The OH PLIF measurement allows for correlation of the flame position with the dominant vortical structures, showing that the leeward flame branch lies slightly above, as well as, in the region between the CVP cores.