The shock–vortex interaction patterns affected by vortex flow regime and vortex models

We have used a third-order essentially non-oscillatory method to obtain numerical shadowgraphs for investigation of shock–vortex interaction patterns. To search different interaction patterns, we have tested two vortex models (the composite vortex model and the Taylor vortex model) and as many as 47 parametric data sets. By shock–vortex interaction, the impinging shock is deformed to a S-shape with leading and lagging parts of the shock. The vortex flow is locally accelerated by the leading shock and locally decelerated by the lagging shock, having a severely elongated vortex core with two vertices. When the leading shock escapes the vortex, implosion effect creates a high pressure in the vertex area where the flow had been most expanded. This compressed region spreads in time with two frontal waves, an induced expansion wave and an induced compression wave. They are subsonic waves when the shock–vortex interaction is weak but become supersonic waves for strong interactions. Under a intermediate interaction, however, an induced shock wave is first developed where flow speed is supersonic but is dissipated where the incoming flow is subsonic. We have identified three different interaction patterns that depend on the vortex flow regime characterized by the shock–vortex interaction.      

Experimental investigations of a swirling jet in both stationary and rotating surroundings

The ‘plug’ flow emerging from a long rotating tube into a large stationary reservoir was used in the experimental investigation of swirling jets with Reynolds numbers, Re = 600, 1,000 and 2,000, and swirl numbers, S = ΩR/U, in the range 0–1.1, to cover flow regimes from the non-rotating jet to vortex breakdown. Here Ω is the nozzle rotation rate, R is the radius of the nozzle exit, and U is the mean mass axial velocity. The jet was more turbulent and eddies shed faster at larger Re. However the flow criticality and shear layer morphology remained unchanged with Re. After the introduction of sufficient rotation, co-rotating and counter-winding helical waves replaced vortex rings to become the dominant vortex structure. The winding direction of the vortex lines suggests that Kelvin–Helmholtz and generalized centrifugal instability dominated the shear layer. A quantitative visualization study has been carried out for cases where the reservoir was rotating independently with S _a  = Ω _a R/U = ±0.35, ±0.51 and ±0.70 at Re = 1,000 and 2000, where Ω _a is the rotation rate of the reservoir. The criterion for breakdown was found to be mainly dependent on the absolute swirl number of the jet, S. This critical swirl number was slightly different in stationary and counter-swirl surroundings but obviously smaller when the reservoir co-rotated, i.e. S _c  = 0.88, 0.85 and 0.70, respectively. These results suggest that the flow criticality depends mainly on the velocity distributions of the vortex core, while instabilities resulting from the swirl difference between the jet and its ambient seem to have only a secondary effect.     

Laminar co-rotating Batchelor vortex merging

The dynamics of laminar co-rotating vortex pairs without axial flow have been recently thoroughly studied through theoretical, experimental and numerical studies, which revealed different instabilities contributing to the decay of the vortices. In this paper, the objective is to extend the analysis to the case of co-rotating vortices with axial flow at low Reynolds numbers. A high-order incompressible Navier–Stokes flow solver is used. The momentum equations are spatially discretized on a staggered mesh by finite differences and all derivatives are evaluated with 10th order compact finite difference schemes with RK-4 temporal discretization. The initial condition is a linear superposition of two co-rotating circular Batchelor vortices with q = 1. It is found that there is an initial evolution that resembles the evolution that single q = 1 vortices go through. Azimuthal disturbances grow and result in the appearance of large-scale helical sheets of vorticity. With the development of these instability waves, the axial velocity deficit is weakened. The redistribution of both angular and axial momentum between the core and the surroundings drives the vortex core to a more stable configuration, with a higher q value. After these processes, the evolution is somewhat similar to a pair of co-rotating Lamb–Oseen vortices. A three-dimensional instability develops, with a large band of unstable modes, with the most amplified mode corresponding scaling with the vortex initial separation distance. P. J. S. A. Ferreira de Sousa wishes to acknowledge the support of FCT—SFRH/BD/1129/2000 and SFRH/BPD/21778/2005.     

Fine structure of a stratified flow near a flat-plate surface

The pattern of disturbances arising during the motion of a strip along a horizontal surface in a continuously stratified fluid with identified upstream and attached internal waves, boundary layers, and edge singularities is calculated in the liner approximation. The flow pattern behind a flat plate moving with a constant velocity in a continuously stratified fluid is studied with the use of the optical schlieren technique; transformation of waves and finely structured elements of the flow with increasing plate velocity is analyzed. The calculated and experimentally observed patterns of internal waves at low velocities are demonstrated to be in good agreement. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 48, No. 6, pp. 77–91, November–December, 2007.     

Stability of a viscous subsonic swirl flow

The results of calculating the stability of a three-dimensional swirl flow of a viscous heat-conducting gas are presented. The stability characteristics are determined using the linear time-dependent theory of plane-parallel flow stability. The main undisturbed axisymmetric vortex flow was determined numerically using a quasi-cylindrical approximation for the complete set of Navier-Stokes equations. The circulation of the peripheral velocity in the cocurrent flow surrounding the viscous vortex core was assumed to be constant. In analyzing the stability, nonaxisymmetric perturbations in the shape of waves traveling along the vortex axis with both positive and negative wavenumbers were considered; in these two cases the perturbation rotation is either the same or opposite in sense to the rotation in the vortex core. Neutral stability curves are determined for various values of the swirling parameter and the cocurrent flow Mach number. Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 3, pp. 50–59, May–June, 1998.     

Isolated Singularities of the 1D Complex Viscous Burgers Equation

The Cauchy problem for the 1D real-valued viscous Burgers equation u _t +uu _x  = u _xx is globally well posed (Hopf in Commun Pure Appl Math 3:201–230, 1950). For complex-valued solutions finite time blow-up is possible from smooth compactly supported initial data, see Poláčik and Šverák (J Reine Angew Math 616:205–217, 2008). It is also proved in Poláčik and Šverák (J Reine Angew Math 616:205–217, 2008) that the singularities for the complex-valued solutions are isolated if they are not present in the initial data. In this paper we study the singularities in more detail. In particular, we classify the possible blow-up rates and blow-up profiles. It turns out that all singularities are of type II and that the blow-up profiles are regular steady state solutions of the equation.     

Flow Bifurcations in a Thin Gap Between Two Rotating Spheres

This paper presents the use of a parameter continuation method and a test function to solve the steady, axisymmetric incompressible Navier–Stokes equations for spherical Couette flow in a thin gap between two concentric, differentially rotating spheres. The study focuses principally on the prediction of multiple steady flow patterns and the construction of bifurcation diagrams. Linear stability analysis is conducted to determine whether or not the computed steady flow solutions are stable. In the case of a rotating inner sphere and a stationary outer sphere, a new unstable solution branch with two asymmetric vortex pairs is identified near the point of a symmetry-breaking pitchfork bifurcation which occurs at a Reynolds number equal to 789. This solution transforms smoothly into an unstable asymmetric 1-vortex solution as the Reynolds number increases. Another new pair of unstable 2-vortex flow modes whose solution branches are unconnected to previously known branches is calculated by the present two-parameter continuation method. In the case of two rotating spheres, the range of existence in the (Re ₁ , Re ₂ ) plane of the one and two vortex states, the vortex sizes as a function of both Reynolds numbers are identified. Bifurcation theory is used to discuss the origin of the calculated flow modes. Parameter continuation indicates that the stable states are accompanied by certain unstable states. Received 26 November 2001 and accepted 10 May 2002 Published online 30 October 2002 Communicated by M.Y. Hussaini     

On the formation of Taylor–G?rtler vortices in RMF-driven spin-up flows

This paper is concerned with a liquid metal flow driven by a rotating magnetic field inside a stationary cylinder. We consider especially the secondary meridional flow during the time when the fluid spins up from rest. The developing flow is investigated experimentally and by direct numerical simulations. The vertical profiles of the axial velocity are measured by means of the ultrasound Doppler velocimetry. Evolving instabilities in the form of Taylor–G?rtler vortices have been observed just above the instability threshold (Ta ≥ 1.5· Ta ^cr). The rotational symmetry may survive over a distinct time even if a first Taylor–G?rtler vortex pair has been formed as closed rings along the cylinder perimeter. The transition to a three-dimensional flow in the side layers results from the advection or a precession and splitting of the Taylor–G?rtler vortex rings. The predictable behaviour of the Taylor–G?rtler vortices disappears with increasing magnetic field strength. The numerical simulations agree very well with the flow measurements.     

Application of digital particle image velocimetry to insect aerodynamics: measurement of the leading-edge vortex and near wake of a Hawkmoth

Some insects use leading-edge vortices to generate high lift forces, as has been inferred from qualitative smoke visualisations of the flow around their wings. Here we present the first Digital Particle Image Velocimetry (DPIV) data and quantitative analysis of an insect’s leading-edge vortex and near wake at two flight speeds. This allows us to describe objectively 2D slices through the flow field of a tethered Tobacco Hawkmoth (Manduca sexta). The near-field vortex wake appears to braodly resemble elliptical vortex loops. The presence of a leading-edge vortex towards the end of the downstroke is found to coincide with peak upward force production measured by a six-component force–moment balance. The topology of Manduca’s leading-edge vortex differs from that previously described because late in the downstroke, the structure extends continuously from wingtip across the thorax to the other wingtip.     

Scalar mixing in the field of a gaseous laminar line vortex

An experimental study of scalar mixing in a laminar vortex is presented for vortices generated between two gas streams flowing parallel to each other in a rectangular flow channel. An isolated line vortex is initiated on demand by momentarily increasing one stream velocity in relation to the other using an electromagnetically actuated piston. The temporal piston motion profile is tailored to generate vortices of different strengths corresponding to vortex Reynolds numbers, Re≡Γ/2πν=130–210. Evolution of mixing is monitored by laser-induced fluorescence of acetone vapor premixed into one of the gas streams as the vortex structure evolves with increasing downstream distance from its point of origin. Vortex is generated by pulsing either of the gas streams (seeded or unseeded stream). Vortex initiation process affects the abundance of the gas in the vortex core from the pulsed stream. Spatial mixing statistics are obtained by determining scalar concentration probability density functions (pdf) and the mean mixed fluid concentrations obtained from these pdfs. It is found that the interfacial area generation as a result of vortex kinematics and molecular diffusion along this interface are principally responsible for mixing. The mean mixed fluid concentration in the vortex interaction region scales with the product of vortex circulation and the elapsed time of interaction. These results are similar to those found in liquid mixing experiments, but the rate of mixing is significantly higher due to higher diffusivity of gases.     

Explosive discharge of a gas-saturated liquid from channels and tanks

A mathematical model for the discharge of a gas-saturated liquid from cylindrical channels is developed. Two limiting cases of linear and quadratic, relations between the flow friction force and the flow velocity are considered. It is established that the process of evacuation, from a semi-infinite channel consists of two stages. In the initial stage, the flow drag can be ignored, and the process of discharge is described by a Riemann wave solution. For the next stage, in which inertia is insignificant, nonlinear equations are obtained and self-similar solutions are constructed for them. The problem of flow through a slot in a tank of finite volume is solved. It is shown that the discharge proceeds either in a gas-dynamic choking regime or in a subsonic regime, depending on the conditions inside the tank and at the outlet. Examples of numerical calculations are given. Institute of Mechanics, Ufa Scientific Center, Russian Academy of Sciences, Ufa 450000 Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 40, No. 1, pp. 64–73, January–February, 1999.     

Flow structures and force characteristics for flat plate in oscillatory flows withKc number from 2 to 40 and in combined flows

The evolution of wake structures and variation of the forces on a flat plate in harmonic oscillatory and in-line combined flows are obtained numerically by improved discrete vortex method. For the oscillatory oncoming flow cases, wyenKc number varies from 2 to 40, the vortex pattern changes from a “harmonic wave” shaped (in a range of smallKc number) to a slight inclined “harmonic wave” shaped (in a range of moderateKc numbers), then to inclined vortex clusters with an angle of 50° to the oncoming flow direction (atKc=20), at last, asKc number becomes large, the vortex pattern is like a normal Karman vortex street. The well predicted drag and inertia force coefficients are obtained, which are more close to the results of Keulegan & Carpenter's experiment as compared with previous vortex simulation by other authors. The existence of minimum point of inertia force coefficientC _m nearKc=20 is also well predicted and this phenomenon can be interpreted according to the vortex structure. For steady-oscillatory in-line combined flow cases, the vortex modes behave like a vortex street, exhibit a “longitudinal wave” structure, and a vortex cluster shape corresponding to the ratios ofU _m toU ₀ which are ofO (10⁻¹)O(1) andO(10), respectively. The effect on the prediction of forces on the flat plate from the disturbance component in a combined flow has been demonstrated qualitatively. In addition to this, the lock in phenomenon of vortex shedding has been checked. The project supported by National Natural Science Foundation of China & LNM, Institute of Mechanics, CAS     

Non-linear dynamics of semi-dilute polydisperse polymer solutions in microfluidics: effects of flow geometry

The non-linear dynamics of a semi-dilute (c/c* = 15) polydisperse polyethylene oxide (PEO) solution in microfluidics are studied experimentally using benchmark contraction–expansion flow geometries with three contraction–expansion ratios (4:1:4, 8:1:8 and 16:1:16) and two narrow channel lengths (L _c/D _h = 53 and 5.3, where L _c is the length of the narrow channel and D _h is its hydraulic diameter). Complex flows over a range of elasticity numbers (El), Weissenberg numbers (Wi) and Reynolds numbers (Re) are characterized using micro-particle image velocimetry ( \upmu\upmu-PIV) and pressure drop measurements. The evolution of vortex formation and dynamics has been visualized through a step-flow-rate experiment. Various flow dynamics regimes have been quantified and are presented in a Wi–Re diagram. The experimental results reveal that the contraction ratio can result in qualitatively different vortex dynamics of semi-dilute polymer solutions in microfluidics, whereas the length of the narrow channel merely affects the dynamics at a quantitative level. A single elasticity number, if defined by the size of the narrow channel, is not sufficient to account for the effects of contraction ratio on the non-linear vortex dynamics.     

Convergence of Ginzburg–Landau Functionals in Three-Dimensional Superconductivity

In this paper we consider the asymptotic behavior of the Ginzburg–Landau model for superconductivity in three dimensions, in various energy regimes. Through an analysis via Γ-convergence, we rigorously derive a reduced model for the vortex density and deduce a curvature equation for the vortex lines. In the companion paper (Baldo et al. Commun. Math. Phys. 2012, to appear) we describe further applications to superconductivity and superfluidity, such as general expressions for the first critical magnetic field H _c1, and the critical angular velocity of rotating Bose–Einstein condensates.     

Single cavitation bubble generation and observation of the bubble collapse flow induced by a pressure wave

This study utilizes a U-shape platform device to generate a single cavitation bubble for a detailed analysis of the flow field characteristics and the cause of the counter jet during the process of bubble collapse caused by sending a pressure wave. A high speed camera is used to record the flow field of the bubble collapse at different distances from a solid boundary. It is found that a Kelvin–Helmholtz vortex is formed when a liquid jet penetrates the bubble surface after the bubble is compressed and deformed. If the bubble center to the solid boundary is within one to three times the bubble’s radius, a stagnation ring will form on the boundary when impinged by the liquid jet. The fluid inside the stagnation ring will be squeezed toward the center of the ring to form a counter jet after the bubble collapses. At the critical position, where the bubble center from the solid boundary is about three times the bubble’s radius, the bubble collapse flow will vary. Depending on the strengths of the pressure waves applied, the collapse can produce a Kelvin–Helmholtz vortex, the Richtmyer–Meshkov instability, or the generation of a counter jet flow. If the bubble surface is in contact with the solid boundary, the liquid jet can only move inside-out without producing the stagnation ring and the counter jet; thus, the bubble collapses along the radial direction. The complex phenomenon of cavitation bubble collapse flows is clearly manifested in this study.     

Spatial distribution of the gas phase in an axisymmetric submerged impact jet

Experimental data on the spatial distribution of the gas phase in an axisymmetric impact jet are obtained by the particle image velocimetry/laser-induced fluorescence (PIV/LIF) method. It is shown that the distribution of bubbles in the flow is determined by the dynamics of vortex structures. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 50, No. 4, pp. 33–38, July–August, 2009.     

Steady temperature field associated with the flow of gassy oil in a porous medium

The steady thermal field associated with the flow of gassy oil through a porous medium is investigated with allowance for the Joule-Thomson and degassing effects. A formula is obtained for estimating the temperature anomalies at the well bottom on oil inflow intervals which correspond to a bottom pressure lower than the saturation pressure. Ufa. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No.1, pp. 113–116, January–February, 1994.     

Approximate analytical calculation of the mach configuration of steady shock waves in a plane constricting channel

An approximate analytical model for calculation of the parameters of a steady gas flow inside a plane constricting channel formed by two symmetrically positioned wedges is suggested. A Mach configuration of shock waves (triple point) is formed in the channel when the wedge angles are larger than some critical value. The flow calculation in a constricting channel reduces to the solution of the iterative problem for a system of nonlinear algebraic equations. The configurations of shock waves, the slipstream, and the sonic line are described by the proposed model of a gas flow. A comparison of the results obtained using this model allows a fairly accurate calculation of the Mach stem and the length of the subsonic-flow region. Institute of Theoretical and Applied Mechanics, Siberian Division, Russian Academy of Sciences, Novosibirsk 630090. Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 39, No. 3, pp. 52–58, May–June, 1998.     

Submodel of rotational motions in gas dynamics

An invariant submodel of the equations of gas dynamics constructed on a one-dimensional subalgebra consisting of the sum of operators of rotation and translation in time is studied within the framework of the SUBMODELS program. The system of equations of the submodel is brought to symmetric form. Hyperbolicity conditions for the system are derived. Group analysis is performed and an invariant solution is considered. Isobaric flows are listed. For the simplest of them, characteristics and strong discontinuities are considered. Necessary conditions for existence of solutions without singularities on the axis are derived. Institute of Mechanics, Ural Scientific Center, Russian Academy of Sciences, Ufa 450000. Translated from Prikladnya Mekhanika i Tekhnicheskaya Fizika, Vol. 39, No. 6, pp. 37–45, November–December, 1998.