Vortex motion on a sphere: barrier with two gaps

Techniques based on the conformal mapping and the numerical method of contour dynamics are presented for computing the motion of a finite area patch of constant vorticity on a sphere in the presence of a thin barrier with two gaps. Finite area patch motion is compared with exact point vortex trajectories and good agreement is found between the point vortex trajectories and the centroid motion of finite area patches when the patch remains close to circular. Patch centroids are, in general, closely constrained to follow point vortex trajectories. However, Kelvin’s theorem constrains the circulation about the barrier to be a constant of the motion, thus, forcing a time-dependent volume flux through the gaps. More exotic motion is observed when the through-gap flow forces the vortex patch close to an edge of a barrier, resulting in the vortex splitting with only part of the patch passing through the gap. As the gap width is decreased this effect becomes more dramatic.     

Conformal mapping and efficient boundary element method without boundary elements for fast vortex particle simulations

In this paper, a revitalization of conformal mapping methods applied to fluid flows in two dimensions is proposed. The present work addresses several important issues concerning their application for vortex particle flow solvers. Difficulties of past conformal based method are reviewed. One difficulty concerns the ability of a mapping procedure to represent complicated shapes. The present paper improves past algorithms to be able to map new shapes, including multiply connected domains. A new fast procedure allows transferring a set of points in the mapped simplified plane to the complicated domain and vice versa. After a mapping construction, it is demonstrated how basic exact solutions to potential flow problems with vortices can be put in a new form which provides a faster and more accurate computation than with distributed singularity methods.     

环空流道液流失速水击特性及其成因的研究

环空流道与圆管的结构差异,使得其失速水击特性及其成因亦有不同。为了对其进行分析研究,利用PIV对套管环空和圆管流水击流场进行拍摄,并通过Tecplot显示流场,提取轴向速度和径向速度、等速度线等参数加以分析;同时采用高精度智能动态压力传感器采集了套管环空内、外壁及内管内壁压力,对环空中水击特性进行了进一步研究。结论是:水击发生时环空断面外壁面水击压力大于内壁面水击压力;同初速情况下环空液流水击压力大于圆管流水击压力;环空中水击压力衰减速度快于圆管中;涡流是水击压力衰减以及速度变化的主要因素;水击压力振荡变化主要是压力涡流引起的断面能量的相互转化形成的。     

Computation of Hele-Shaw free boundary problems near obstacles

The time-dependent evolution of source-driven Hele-Shaw free boundary flows in the presence of an obstacle is computed numerically. The Baiocchi transformation is used to convert the Hele-Shaw Laplacian growth problem into a free boundary problem for a streamfunction-like variable u(x, y, t) governed by Poisson’s equation (with constant right-hand side) with the source becoming a point vortex of strength linearly dependent on time. On the free boundary, both u and its normal derivative vanish, and on the obstacle, the normal derivative of u vanishes. Interpreting u as a streamfunction, at a given time, the problem becomes that of finding a steady patch of uniform vorticity enclosing a point vortex of given strength such that the velocity vanishes on the free boundary and the tangential velocity vanishes on the obstacle. A combination of contour dynamics and Newton’s method is used to compute such equilibria. By varying the strength of the point vortex, these equilibria represent a sequence of source-driven growing blobs of fluid in a Hele-Shaw cell. The practicality and accuracy of the method is demonstrated by computing the evolution of Hele-Shaw flow driven by a source near a plane wall; a case for which there is a known exact solution. Other obstacles for which there are no known exact solutions are also considered, including a source both inside and outside a circular boundary, a source near a finite-length plate and the interaction of an infinite free boundary impinging on a circular disc.     

The evolution of an initially circular vortex near an escarpment. Part I: analytical results

The motion of a quasigeostrophic, equivalent-barotropic, initially circular vortex patch near an infinitely long topographic escarpment is studied using f-plane dynamics. There are two time scales in the problem: the advective time scale associated with the vortex, and the time scale for topographic vortex stretching. Analytical progress is possible when these two time scales are well-separated and results are presented here.If topographic vortex stretching dominates advection by the vortex the vortex is said to be ‘weak’. The vortex patch remains circular on the topographic time scale, and dispersive topographic waves rapidly propagate the initial disturbance away from the vicinity of the vortex. Subsequently cross-isobath motion is inhibited, and the vortex moves as though the escarpment were a plane wall. The same behaviour was observed for the motion of a weak singular vortex near an escarpment by Dunn, McDonald and Johnson [7], who named the phenomenon the ‘pseudoimage’ of the vortex.If advection dominates over topographic effects, the vortex is said to be ‘intense’. The vortex also remains circular to leading order, but the relative vorticity produced by the swirl of the vortex is less able to escape the vicinity of the vortex. The vortex follows a similar curved trajectory to those observed for intense vortices on the β-plane. The dipolar mechanism for this behaviour is described. Large time solutions are inhibited by the form of the escarpment topography, but examination of the equations leads to the conclusion that the leading order solution may be predict the motion for times beyond its formal range of validity.     

Anisotropic thermopiezoelectric solids with an elliptic inclusion or a hole under uniform heat flow

The two-dimensional problem of a thermopiezoelectric material containing an elliptic inclusion or a hole subjected to a remote uniform heat flow is studied. Based on the extended Lekhnitskii formulation for thermopiezoelectricity, conformal mapping and Laurent series expansion, the explicit and closed-form solutions are obtained both inside and outside the inclusion (or hole). For a hole problem, the exact electric boundary conditions on the hole surface are used. The results show that the electroelastic fields inside the inclusion or the electric field inside the hole are linear functions of the coordinates. When the elliptic hole degenerates into a slit crack, the electroelastic fields and the intensity factors are obtained. The effect of the heat flow direction and the dielectric constant of air inside the crack on the thermal electroelastic fields are discussed. Comparison is made with two special cases of which the closed solutions exist and it is shown that our results are valid.     

Pseudo-canonical formulation of 3-dimensional vortex motion and vorton model analysis

From the author's pseudo-canonical formulation of 3-dimensional vortex motion, the basic equations for the vorton model are derived. The division of vorton and reduction of meshes of time integration are assumed, in numerical calculations. This method is applied to numerical studies on evolutions of a connected 3-ring vortex and a vortex ring inside or outside a sphere. These seem to show cutting of the vortex filament and the following reconnection or attachment to the wall. The dissipation effect due to the approximation scheme and the singularity of vorton are discussed.     

Numerical study of vortex interactions behind two circular cylinders in bistable flow regime

A two-dimensional discrete-vortex model was used to investigate vortex interaction inside the near wakes of two circular cylinders in side-by-side arrangement within bistable flow regime. Two phases of vortex evolution are mainly identified in the near wakes: a symmetric shedding phase, characterized by two antiphase vortex streets, and a flip-flopping phase, characterized by biased gap flow, switching at irregular intervals. For the flip-flopping phase, vortex amalgamation, vortex pairing and dipole are found. Vortex dynamics of the flow is presented and its possible effects on the flow parameters are discussed. The initiation and transition from the symmetric to flip-flopping phase are caused by the asymmetry of one of the gap vortices. Flow visualization and quantitative results obtained seem to support the findings from the model.     

Review: Formulating efficient finite-element codes for flows in regular domains

Many fluid flow problems of current interest occur in domains that are mappable to a rectangle or a box; conformal mappings are particularly useful in this regard. We are concerned here with the efficient solution of such problems using finite elements. The central issue is the element choice, and this issue is addressed in terms of operation counts, computer memory and I/0 requirements, and the extent to which code vectorization is possible. It is concluded that rectangular (box) elements generally lead to more efficient algorithms that triangular (tetrahedral) elements. A synthesis of algorithms, based on bilinear (trilinear) elements, is presented. The algorithms have the attributes of simplicity, accuracy, stability and straightforward incorporation of boundary conditions. For bilinear and trilinear elements, it is found that product and first-derivative terms are well-handled by the Galerkin FE method, but that it is advantageous to go outside of the Galerkin framework when treating second-derivative terms. It is particularly important to consider the form of the governing equations, vis-à-vis the choice of staggered, non-staggered and/or mixed-order elements, and to choose an appropriate time scheme. The described techniques have been successfully applied to a variety of problems in regular domains, including the solution of the three-dimensional time-dependent hydrostatic primitive equations; these are stiff and include first and second derivative terms, non-linearities and variable coefficients due to a conformal mapping.     

The effect of wall cooling on compressible Görtler vortices

It is known that in adiabatic boundary layer flow over a curved surface the detailed structure of the spanwise periodic Görtler vortex instability varies markedly over the range of spanwise wavelength. At short wavelengths the modes tend to be concentrated in a well-defined thin zone located within the boundary layer. As the vortex wavenumber diminishes so the region of vortex activity is first driven to the bounding wall but subsequently expands to cover the entire boundary layer at which stage the modes take on a principally inviscid form. At yet longer wavelengths the vortices are given by the solution of an interactive multi-deck structure which has some similarities with that for Tollmien–Schlichting waves.In this work we investigate how the application of wall cooling affects the above scenario. It is shown how cooling both restricts the range of mode types and gives rise to two new structures. The first, for moderate cooling and which relates to longer wavelengths, is interactive in nature. Here the viscous–inviscid interaction between an essentially inviscid Görtler problem, albeit for an effective basic flow which in its general form has a non-standard near-wall structure, and a viscous sublayer is provided by novel boundary conditions. Shorter wavelength vortices are largely unaffected by wall cooling unless this is quite severe. However when this degree of cooling is applied, the vortices take on a fully viscous form and are confined to a thin region next to the bounding wall wherein the basic flow assumes an analytic form. Numerical solutions are obtained and we provide evidence as to how the two new structures are related both to each other and to the previously known uncooled results.     

Vortex dynamics in a model left ventricle during filling

Analysis of the properties of the left ventricle flow during filling (diastole) is known to allow early detection of potential malfunctions during the fundamental heart pumping phase (systole). Diagnoses are now usually based on Doppler measurement of the flow velocity for which interpretative schemes and quantitative references are sought. The flow inside an ideal model of the left ventricle is here studied numerically by a finite difference method in prolate spheroid moving coordinates. The axisymmetric assumption is employed in this first study. The vortex dynamics is characterised by the formation of a wake vortex attached at the valvular edge which is shed at the end of the ventricle expansion. Major features are its translation by self-induced velocity and vortex-induced separations from the ventricle internal wall. Results are represented as, and compared with, clinical data showing a good general agreement and allowing an improved physical interpretation of the latter.     

Assessment of ghost-cell based cut-cell method for large-eddy simulations of compressible flows at high Reynolds number

Large-eddy simulations (LES) of high Reynolds number flows are performed using a non-body conformal method in conjunction with a wall model. We use a simple wall function to model the wall-shear stress and the truncation error of the numerical discretization to model the sub-grid scale turbulence (implicit LES), although these can be easily replaced if necessary. The validation cases are: turbulent flow through an inclined channel, turbulent flow over a wavy surface, and supersonic flow over a circular cylinder. Since the near-wall grids are naturally coarse, the key is to use a method that is capable of capturing the flow dynamics accurately in the vicinity of the interface. Towards the purpose, we develop a Cartesian cut-cell method, referred to as the ghost-cell based cut-cell method (GC-CCM), in the context of fully compressible solutions of Navier–Stokes equations. This method employs ghost-cells inside the solid interface such that the local spatial reconstruction remains consistent everywhere including in the vicinity of the boundary. In order to capture the near-wall flow behavior more accurately with coarse grids, this method decomposes cell faces of merged cells and computes fluxes through each decomposed segment separately. The objective of this work is to qualify whether the proposed method can accurately represent the high Reynolds number flows in the vicinity of immersed interfaces. To analyze the performance of the proposed method, we compare the results to the corresponding numerical results from the two other non-body conformal methods, namely the ghost-cell based immersed boundary method (GCIBM) and standard cut-cell method (S-CCM), that are implemented in the same numerical solver. The comparison demonstrates that the proposed method is capable of capturing near-wall flows relatively accurately with coarse grids.     

Exact calculation of the flow of a staggered tandem cascade with moving second blade row

Summary  The plane flow around a tandem cascade of flat plates is calculated by means of conformal mapping. The blades of the two rows are perpendicular to each other. The first row is stationary, the second row moves with constant velocity. The conformal mapping will be constructed by a “mapping flow”. The blades of one row are stream lines and those of the other row are potential lines of the flow. By conformal mapping, the physical flow around the tandem cascade of the physical ζ-plane is converted into a flow between infinitely long straight walls in the z-plane, each wall corresponding to one of the blades. The conditions far upstream and far downstream of the cascade are represented by source-vortices. In the z-plane, the boundary conditions may be easily fulfilled by reflection and repetition of the source-vortices, and the flow may be calculated by well-known methods. The physical flow searched for is then obtained by inverse mapping. Received 24 July 2000; accepted for publication 6 December 2000     

Dual cylinder inviscid vortex sound of low Mach number

The unsteady motions of an inviscid vortex under the influence of a cylinder pair in the presence of a low Mach number mean flow and the corresponding sound generation are examined in the present study. The two cylinders are in close proximity. A semi-analytical approach using the conformal mapping together with the potential theory is adopted. The results show that the vortex will interact intensively with the cylinders under the right combinations of mean flow direction and initial vortex position. Such interactions result in a high rate of change of vortex propagation velocity, strong fluctuating forces on cylinder and strong sound radiations. However, it is found that much stronger acoustic energy radiation will result when the vortex approaches the cylinder pair from the bottom than from the top, unless the mean flow is nearly perpendicular to the horizontal cylinder pair axis. Stronger sound radiation is also observed for the identical cylinder cases in general, except the flow direction is close to some critical values.     

A new calculus for two-dimensional vortex dynamics

This article provides a user’s guide to a new calculus for finding the instantaneous complex potentials associated with point vortex motion in geometrically complicated planar domains, with multiple boundaries, in the presence of background flows. The key to the generality of the approach is the use of conformal mapping theory together with a special transcendental function called the Schottky–Klein prime function. Illustrative examples are given.     

Finite Difference Method for the Acoustic Radiation of an Elastic Plate Excited by a Turbulent Boundary Layer: A Spectral Domain Solution

A finite difference method is developed to study, on a two-dimensional model, the acoustic pressure radiated when a thin elastic plate, clamped at its boundaries, is excited by a turbulent boundary layer. Consider a homogeneous thin elastic plate clamped at its boundaries and extended to infinity by a plane, perfectly rigid, baffle. This plate closes a rectangular cavity. Both the cavity and the outside domain contain a perfect fluid. The fluid in the cavity is at rest. The fluid in the outside domain moves in the direction parallel to the system plate/baffle with a constant speed. A turbulent boundary layer develops at the interface baffle/plate. The wall pressure fluctuations in this boundary layer generates a vibration of the plate and an acoustic radiation in the two fluid domains. Modeling the wall pressure fluctuations spectrum in a turbulent boundary layer developed over a vibrating surface is a very complex and unresolved task. Ducan and Sirkis [1] proposed a model for the two-way interactions between a membrane and a turbulent flow of fluid. The excitation of the membrane is modeled by a potential flow randomly perturbed. This potential flow is modified by the displacement of the membrane. Howe [2] proposed a model for the turbulent wall pressure fluctuations power spectrum over an elastomeric material. The model presented in this article is based on a hypothesis of one-way interaction between the flow and the structure: the flow generates wall pressure fluctuations which are at the origin of the vibration of the plate, but the vibration of the plate does not modify the characteristics of the flow. A finite difference scheme that incorporates the vibration of the plate and the acoustic pressure inside the fluid cavity has been developed and coupled with a boundary element method that ensures the outside domain coupling. In this paper, we focus on the resolution of the coupled vibration/interior acoustic problem. We compare the results obtained with three numerical methods: (a) a finite difference representation for both the plate displacement and the acoustic pressure inside the cavity; (b) a coupled method involving a finite difference representation for the displacement of the plate and a boundary element method for the interior acoustic pressure; (c) a boundary element method for both the vibration of the plate and the interior acoustic pressure. A comparison of the numerical results obtained with two models of turbulent wall pressure fluctuations spectrums - the Corcos model [3] and the Chase model [4] - is proposed. A difference of 20 dB is found in the vibro-acoustic response of the structure. In [3], this difference is explained by calculating a wavenumber transfer function of the plate. In [6], coupled beam-cavity modes for similar geometry are calculated by the finite difference method. This revised version was published online in July 2006 with corrections to the Cover Date.     

Experimental investigation of the interaction between a plane wall jet and a parallel offset jet

The dual-jet flow generated by a plane wall jet and a parallel offset jet at an offset ratio of d/w = 1.0 has been investigated using Particle Image Velocimetry (PIV). The particle images are captured, processed, and subsequently used to characterize the flow in terms of the 2D velocity and vorticity distributions. Statistical characteristics of the flow are obtained through ensemble averaging of 360 instantaneous velocity fields. Also presented is a time series of instantaneous flow fields to illustrate the dynamic interaction between the two jets. Results reveal that the near field of the flow is characterized by a periodic large-scale Karman-like vortex shedding similar to what would be expected in the wake of a bluff body. The existence of the Karman-like vortices results in periodic interactions between the two jets; in addition, these vortices produce noticeable impact on the jet outer layers, i.e., the free shear layer of the offset jet and the wall boundary layer of the wall jet. A schematic of vortex/shear layer interaction is proposed to illustrate the flow pattern.     

船用燃气轮机进气滤清器惯性级内的流场计算和实验验证

本文提出了一套求解船用燃气轮机进气滤清器流道流场的数值方法,成功地计算了流道内流场的状态,给出了各种不同型号流道的气动特性,对指导滤清器的设计有较大的现实意义, 在这套方法中,我们应用上风差分来逼近二维、非定常、粘性、不可压缩流体非守恒型的N-S方程,提出了一种可计算雷诺数高达上万的粘性流的差分方程,考察了这种差分方程的稳定性,收敛性、精度和人工粘性,本文还提出了处理一些边界拐点处壁涡的计算方法,实际算例表明,使用本文提出的差分方程和壁涡处理方法给出的计算结果和实验吻合良好。     

Vortex dynamics of a sphere wake in proximity to a wall

Toward getting the vortex dynamics characteristics and wake structure of a sphere in proximity to a wall, the effect of a proximal flat plate on the wake of a stationary sphere is investigated via direct numerical simulation. The vortex shedding process and the significant variation of the wake structure are described in detail. The drag coefficient reduces and the wake structure of the sphere becomes complex due to the combined effect of the wake flow and the wall. A jet flow forms between the sphere and the flat plate, which suppresses the vortex separation on the bottom of the sphere. The asymmetric distributions of the coherent structures and the recirculation region behind the sphere are discussed. Besides vortex shedding patterns, the time-averaged velocity distribution, vortex dynamics, distribution regularities of turbulent kinetic energy and enstrophy are investigated.     

CALCULATION OF STEADY AND OSCILLATING FLOWS IN TUBES USING A VORTICITY TRANSPORT ALGORITHM

Steady and oscillating axisymmetric tube flows are modelled using a vorticity transport algorithm. The axisymmetric convective –diffusive Navier–Stokes equations are solved using a splitting technique. Axisymmetric ring vortex filaments are introduced on the walls and subsequently convected and diffused throughout the flow field. An axisymmetric equation similar to the Oseen diffusion equation is used to diffuse the ring vortex filaments. Vorticity is reflected from the tube walls using two techniques. Results are presented for the developing Poiseuille flow and for the developed flow in the form of the entrance length and the axial velocity and vorticity profiles. Good agreement is achieved with a finite difference method in the developing region of Poiseuille flow. The developed flow results are compared with the analytical solutions. The developed profiles of velocity and vorticity have errors of less than 0ċ3 per cent for both methods of dealing with reflection of diffusion at the bounding surfaces and similar accuracy is obtained for the velocity profiles in oscillating flow except at the wall. Oscillating flow is produced with a discretized sinusoidal piston motion. Velocity profiles, boundary layer thickness and entrance length are presented for oscillating flow. Good agreement is achieved for low-Womersley-number non-dimensional frequency. At higher values of this parameter, flows are inaccurately simulated, because the number of piston positions used to discretize the piston motion is inversely proportional to the non-dimensional frequency.