Interactions of three-dimensional vicous axisymmetric vortex rings with a free surface

The unsteady nonlinear interaction of three-dimensional vortices with a free surface is a great challenge in fluid mechanics, which has deep theoretical significance and important practical background. Applying the three-dimensional VOF method, the interactions of three-dimensional axisymmetric vortex rings with a free surface in an incompressible viscous fluid are numerically simulated. The influence of the Froude number and the surface tension are studied and the evolution of the vorticity, the trajectories of the vortex rings and the baroclinic vorticity on the surface are obtained. The results agreed well with the experiments reported in the literature. The project supported by the National Natural Science Foundation of China     

Study of water waves with submerged obstacles using a vortex method with Helmholtz decomposition

The present study develops a 2‐D numerical scheme that combines the vortex method and the boundary integral method by a Helmholtz decomposition to investigate the interaction of water waves with submerged obstacles. Viscous effects and generation of vorticity on the free surface are neglected. The second kind of Fredholm integral equations that govern the strengths of vortex sheets along boundaries are solved iteratively. Vorticity is convected and diffused in the fluid via a Lagrangian vortex (blob) method with varying cores, using the particle strength exchange method for diffusion, with particle redistribution. A grid‐convergence study of the numerical method is reported. The inviscid part of the method and the simulation of the free‐surface motion are tested using two calculations: solitary wave propagation in a uniform channel and a moving line vortex in the fluid. Finally, the full model is verified by simulating periodic waves travelling over a submerged rectangular obstacle using nonuniform vortex blobs with a mapping of the redistribution lattice. Overall, the numerical model predicts the vortices' evolution and the free‐surface motion reasonably well. Copyright © 2008 John Wiley & Sons, Ltd.     

Coaxial interactions of two vortex rings or of a ring with a body

Inviscid coaxial interactions of two vortex rings, including head-on collisions and leapfrogging motions, are considered using a contour dynamics technique. Interactions of vortex rings with solid bodies are also investigated by combining the contour dynamics technique with a boundary integral equation method. Numerical results show that a clean, successful passage motion is possible for two vortex rings with not too thick cores. In both cases of head-on collisions and leapfrogging motions, very large core deformations are observed when two vortex rings get close to each other. A head-tail structure is formed in the later stage of a head-on collision of two fat vortices. Numerical results also show that a vortex ring will stretch and slow down when it moves toward a solid boundary, will shrink and speed up when it moves away from a solid boundary, and will either translate steadily or approach an oscillating asymptotic state when it is far away from any boundaries. The project supported by The National Education Commission of China and NASA under cooperative grant agreement #NCC5-34.     

Vortex ring generation due to the coalescence of a water drop at a free surface

 Results are presented of an experimental investigation of vortex ring formation by a fluid drop contacting a free surface with negligible velocity. The pool fluid is mixed with fluorescein dye, and a laser sheet is used to illuminate a plane of the flow. A series of representative images is recorded by a CCD camera and speculation is made regarding specific sources of vorticity flux through the free surface. Two scaling analyses previously presented by other investigators are demonstrated to be equivalent under the assumptions of this experiment, and they provide the motivation for a series of test runs in which the duration of the coalescence process, τ_*, is related to variations in drop diameter L and fluid surface tension σ. Experimental results are in agreement with the analyses, showing τ_*∼σ^-1/2 and τ_*∼L ^3/2. Received: 22 December 1995 / Accepted: 15 October 1996     

On variational features of vortex flows

Ideal incompressible fluid is a Hamiltonian system which possesses an infinite number of integrals, the circulations of velocity over closed fluid contours. This allows one to split all the degrees of freedom into the driving ones and the “slave” ones, the latter to be determined by the integrals of motions. The “slave” degrees of freedom correspond to “potential part” of motion, which is driven by vorticity. Elimination of the “slave” degrees of freedom from equations of ideal incompressible fluid yields a closed system of equations for dynamics of vortex lines. This system is also Hamiltonian. The variational principle for this system was found recently (Berdichevsky in Thermodynamics of chaos and order, Addison-Wesly-Longman, Reading, 1997; Kuznetsov and Ruban in JETP Lett 67, 1076–1081, 1998). It looks striking, however. In particular, the fluid motion is set to be compressible, while in the least action principle of fluid mechanics the incompressibility of motion is a built-in property. This striking feature is explained in the paper, and a link between the variational principle of vortex line dynamics and the least action principle is established. Other points made in this paper are concerned with steady motions. Two new variational principles are proposed for steady vortex flows. Their relation to Arnold’s variational principle of steady vortex motion is discussed.      

Simple and efficient numerical methods for vortex sheet motion with surface tension

We present two simple and efficient explicit methods for the vortex sheet with surface tension. The first one is the standard point vortex method, which has been known to be unstable in the presence of surface tension, due to spurious growth of waves of high modes. We show, for the first time, that the standard point vortex method is able to calculate the vortex sheet motion with surface tension by employing a Fourier filtering. The second method is a modification of the Pullin method using central differences for numerical differentiations. This method is more convenient to implement than other spectral methods and is free from the aliasing instability. We give a linear stability analysis for the numerical methods and show results for the long‐time evolution of the vortex sheet. We also propose a new redistribution procedure to control point clustering, by setting limits of minimum and maximum distances between neighboring points. This procedure is found to be very efficient for long‐time computations of the explicit methods. Copyright © 2013 John Wiley & Sons, Ltd.     

Normal and oblique collisions of a vortex ring with a wall

The oblique collision of a vortex ring with a solid wall, atRe=/=1389, has been analysed by the direct simulation of the Navier-Stokes equations in Cartesian coordinates. In accordance with a previous experimental study [1], the secondary vorticity produced at the wall is organized into a loop-like vortex in the region of the ring furthest away from the wall. As the ring approaches the wall, the region closest is subjected to a high rate of stretching which increases the vorticity in the core. The vorticity gradients along the core generate bi-helical vortex lines continually displaced towards the region of the ring furthest away from the wall. The analysis of the vorticity and straining fields revealed that the pressure gradient along the core is responsible for the convective motion that displaces these vortex lines and accumulates secondary vorticity in the region far from the wall. This vorticity rolls up and forms a secondary structure which by self-induction moves away from the wall.The fundamental role of the differential stretching has been demonstrated by comparing the case of oblique collision with that of normal collision and with the collision of a two-dimensional vortex pair with an oblique wall.

---

Sommario L'interazione di un vortice ad anello con una parete obliqua, aRe=1389, è stata analizzata mediante la simulazione diretta delle equazioni di Navier-Stokes in coordinate cartesiane. In accordo con un precedente esperimento [1] è stato evidenziato che la vorticità secondaria, prodotta alla parete, si organizza in una strutura vorticosa a loop nella regione dell'anello più lontana dalla parete. Quando il vortice si avvicina alla parete, la parte più vicina è soggetta ad un'elevata deformazione che aumenta il valore della vorticità nel core. La distribuzione non uniforme di vorticità lungo il core del vortice genera delle linee di vorticità elicoidali che vengono transportate verso la regione dell'anello più lontana dalla parete. L'analisi dei campi di vorticità e di deformazione ha rivelato che il gradiente di pressione, dovuto al campo di deformazione non uniforme lungo il core del vortice, è responsabile di un moto convettivo che trasporta le linee di vorticità ed accumula la vorticità secondaria nella regione del vortice più lontana dalla parete, dove la struttura secondaria viene generata.Il ruolo fondamentale della deformazione non uniforme è stato evidenziato mediante il confronto della collisione obliqua coni casi di collisione normale e di collisione di una coppia di vortici bidimensionali con una parete obliqua.

     

Dynamics of ellipsoidal cavities in fluid

A mathematical model of the hydrodynamics of free closed surfaces in a fluid is expounded. It is used for studying the dynamics of ellipsoidal cavities during their development. The model is based on a system of differential equations that accounts for the influence exerted on the dynamics of cavities by various perturbations such as gravity, surface tension, viscosity, and geometrical features of the cavity. Solving this system makes it possible to determine the hydrodynamic characteristics of the flow around the cavity and to plot cavity shapes depending on time and flow regimes. Characteristic features of the development of such cavities under gravity and surface tension are established __________ Translated from Prikladnaya Mekhanika, Vol. 42, No. 2, pp. 24–31, February 2006.     

Free-surface deformation due to spiral flow owing to a source/sink and a vortex in Stokes flow

The free-surface shape and cusp formation are analyzed by considering a viscous flow arising from the superposition of a source/sink and vortex below the free surface where the strength of the source and vortex are arbitrary. In the analysis, Stokes’ approximation is used and surface tension effects are included, but gravity is neglected. The solution is obtained analytically by using conformal mapping and complex function theory. From the solution, shapes of the free surface are obtained, and the formation of a cusp on the free surface is discussed. Above some critical capillary number with a sink, the free-surface shape becomes singular and an apparent cusp should form on the free surface below a real fluid. On the other hand, no cusp would occur for sources of zero or positive strength. Typical streamline patterns are also shown for some capillary numbers. As the capillary number vanishes, the solution is reduced to a linearized potential flow solution.     

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.      

Steady Free-surface Flow Past an Uneven Channel Bottom

New free-surface flows past a semi-infinite ‘step’ in the bottom of a channel are considered. Surface tension is neglected but gravity is included in the dynamic boundary condition. Fully nonlinear solutions are computed by boundary integral equation methods. Additional weakly nonlinear solutions are derived analytically. A thorough analysis of the weakly nonlinear problem provides a systematic approach to identify all the possible types of solutions and the number of independent parameters.     

波涡相互作用研究的某些进展(Ⅱ)

<正> 5 波涡共振 从第3节的感受性问题再前进一步,自然要问在什么条件下入射波激发起涡中之波的最大响应.这就导致了入射波与层状或轴状涡中受激波之间共振的概念,简称波涡共振.一般说来,在流体内部若有两个或多个波相会,它们将互相穿透而沿原来的方向离去.但若它们的波矢量和频率满足一定的关系(参见Craik 1985),就会在相会点产生新的波.2阶扰动的振幅可达到1阶扰动振幅的量级,而且流场中会出现一些重要的独特性质.这就是流体内部波共振,波涡共振是其一类情形.      

Flow of an Inclined Surface-Charged Layer in a Longitudinal Electric Field

The motion of an inclined layer under the influence of an electric force acting in the opposite direction to the gravity force is investigated theoretically and experimentally. The motion of the layer is driven by the tension tangential to the free boundary associated with the interaction between the surface charge and the longitudinal electric field. The condition of ascent of the fluid along an inclined dielectric plane is found. The domains of stable layer motion are calculated.     

A model of stratified entrainment using vortex persistence

A new model is proposed for the entrainment rate by vortices across stratified interfaces. In the model, different entrainment regimes are distinguished by the conventional parameters Richardson, Reynolds, and Schmidt number as well as a new parameter, the “vortex persistence”. Vortex persistence is defined as the number of rotations a vortex makes during the time it moves its own diameter with respect to the interface. It is further proposed that the concept of vortex persistence is important whenever a vortex is near any kind of surface, either stratified or solid. The model is in accord with most field and laboratory observations in a variety of stratified and bounded flows, including measurements of wall heat transfer and vortex formation in starting jets.     

Convective mass transfer across fluid interfaces in straight angular pores

Steady convective mass transfer to or from fluid interfaces in pores of angular cross-section is theoretically investigated. This situation is relevant to a variety of mass transport process in porous media, including the fate of residual non-aqueous phase liquid ganglia and gas bubbles. The model incorporates the essential physics of capillarity and solute mass transfer by convection and diffusion in corner fluid filaments. The geometry of the corner filaments, characterized by the fluid–fluid contact angle, the corner half-angle and the interface meniscus curvature, is accounted for. Boundary conditions of zero surface shear (‘perfect-slip’) and infinite surface shear (‘no-slip’) at the fluid–fluid interface are considered. The governing equations for laminar flow within the corner filament and convective diffusion to or from the fluid–fluid interface are solved using finite-element methods. Flow computations are verified by comparing the dimensionless resistance factor and hydraulic conductance of corner filaments against recent numerical solutions by Patzek and Kristensen (J. Colloid Interface Sci 236, 305–317 2001). Novel results are obtained for the average effluent concentration as a function of flow geometry and pore-scale Peclet number. These results are correlated to a characteristic corner length and local pore-scale Peclet number using empirical equations appropriate for implementation in pore network models. Finally, a previously published “2D-slit” approximation to the problem at hand is checked and found to be in considerable error.     

Surface Water Waves and Tsunamis

Because of the enormous earthquake in Sumatra on December 26, 2004, and the devastating tsunami which followed, I have chosen the focus of my mini-course lectures at this year’s PASI to be on two topics which involve the dynamics of surface water waves. These topics are of interest to mathematicians interested in wave propagation, and particularly to Chilean scientists, I believe, because of Chile’s presence on the tectonically active Pacific Rim. My first lecture will describe the equations of fluid dynamics for the free surface above a body of fluid (the ocean surface), and the linearized equations of motion. From this, we can predict the travel time of the recent tsunami from its epicenter off of the north Sumatra coast to the coast of nearby Thailand, the easy coasts of Sri Lanka and south India, and to Africa. In fact, the signal given by ocean waves generated by the Sumatra earthquake was felt globally; within 48 h distinguishable tsunami waves were measured by wave gages in Antarctica, Chile, Rio di Janeiro, the west coast of Mexico, the east coast of the United States, and at Halifax, Nova Scotia. To describe ocean waves, we will formulate the full nonlinear fluid dynamical equations as a Hamiltonian system [19], and we will introduce the Greens function and the Dirichlet-Neumann operator for the fluid domain along with the harmonic analysis of the theory of their regularity. From an asymptotic theory of scaling transformations, we will derive the known Boussinesq-like systems and the KdV and KP equations, which govern the asymptotic behavior of tsunami waves over an idealized flat bottom. When the bottom is no longer assumed to be perfectly flat, a related theory [6, 13] gives a family of model equations taking this into account. My second lecture will describe a series of recent results in PDE, numerical results, and experimental results on the nonlinear interactions of solitary surface water waves. In contrast with the case of the KdV equations (and certain other integrable PDE), the Euler equations for a free surface do not admit clean (‘elastic’) interactions between solitary wave solutions. This has been a classical concern of oceanographers for several decades, but only recently have there been sufficiently accurate and thorough numerical simulations which quantify the degree to which solitary waves lose energy during interactions [3, 4]. It is striking that this degree of ‘inelasticity’ is remarkably small. I will describe this work, as well as recent results on the initial value problem which are very relevant to this phenomenon [14, 18].     

Evolution of the free surface of a plane channel during filling with a viscous fluid

The behavior of the free surface of a viscous incompressible fluid is studied in the process of filling of a plane channel inclined to the gravity direction. The problem is numerically solved within the framework of the Navier-Stokes equations under the assumption that the fluid is Newtonian and capillary effects can be neglected. The numerical technique is based on the joint use of the SIMPLE and invariant methods. Two filling regimes are found to exist, namely, that with a steady shape of the free surface without the formation of internal cavities and that with cavity formation in the initial filling stage characterized by spouting behavior of the free boundary.     

Numerical modeling of vortex/shock wave interaction and its transformation by localized energy deposition

Three-dimensional unsteady Euler simulations are presented for the interaction of a streamwise vortex with an oblique shock of angle β = 23.3° at Mach 3 and 5. The flowfield features are analyzed for weak, moderate and strong interaction regimes. The details of the free recirculation zone at conditions of subsonic and supersonic flow on the vortex axis are considered. The vortex breakdown under conditions of a subsonic vortex core is characterized by a continuous growth and gradual degeneration of the region, unlike the supersonic core condition wherein a steady recirculation zone is achieved. The possibility of using a localized steady and pulsed periodic energy deposition on the vortex axis for stimulating the breakdown is demonstrated for various interaction regimes. It is shown that the formation of a subsonic wake downstream of an energy source lying on the vortex axis contributes to a more significant growth of the dimensions of the recirculation zone compared to the case when the vortex core remains supersonic. The possibility of achieving the effects similar to the steady case is demonstrated by the effect of a pulsed periodic energy source on the flow under consideration for corresponding equivalence parameters.