Comparison of analytical and numerical solution for optimal vortex diffuser design

The contribution deals with optimization of wing tip devices, so called vortex diffusers. A comparison is given between an analytical approach for obtaining the optimal circulation loading and the results of a numerical investigation using a lifting line method. The purpose of most wing tip devices is to reduce the induced drag of the main wing by converting vortex energy into thrust. In order to achieve an optimal design, a variational formulation originally proposed by Betz and Prandtl for air screws is applied to the circulation distribution of the diffuser blades. In extension to the inviscid formulation, a viscous correction is applied in order to account for frictional forces. In an effort to validate the analytical results, a comparison is given with numerical solutions from a lifting line method. The loading of the diffuser blades is parametrized and optimized with respect to resulting thrust by use of a quasi-Newton gradient method. Comparison shows that, knowing the velocity distribution in the near wake of the wing, considerable decrease of induced drag may be achieved making use of vortex diffusers. Although actual circulation loading may differ between the analytical and numerical estimation, resulting thrust agrees within a few percent. (© 2005 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Third-order accurate conservative method on unstructured meshes for gasdynamic simulations

A third-order accurate finite-volume method on unstructured meshes is proposed for solving viscous gasdynamic problems. The method is described as applied to the advection equation. The accuracy of the method is verified by computing the evolution of a vortex on meshes of various degrees of detail with variously shaped cells. Additionally, unsteady flows around a cylinder and a symmetric airfoil are computed. The numerical results are presented in the form of plots and tables.     

Numerical simulation and visualization of vortical structure transformation in the flow past a sphere at an increasing degree of stratification

Based on numerical simulation and visualization, the vortex structure of the flow past a sphere moving uniformly and horizontally in a linearly (density) stratified viscous fluid with an increasing degree of stratification (with the internal Froude number Fr decreasing from infinity to 0.005) at Re = 100 is analyzed in detail for the first time. The classification of the flow regimes is refined. The direct numerical simulation is based on the method of splitting with respect to physical factors (MERANZH) with an explicit hybrid finite-difference scheme, which is second-order accurate in space, monotone, and has a minimal numerical viscosity and dispersion.     

Vortex pairs and dipoles

Point vortices have been extensively studied in vortex dynamics. The generalization to higher singularities, starting with vortex dipoles, is not so well understood.We obtain a family of equations of motion for inviscid vortex dipoles and discuss limitations of the concept. We then investigate viscous vortex dipoles, using two different formulations to obtain their propagation velocity. We also derive an integro-differential for the motion of a viscous vortex dipole parallel to a straight boundary.     

The Equation of Motion of a Vortex Layer of Small Thickness

The equation governing the evolution of a vortex layer whose thickness is small compared to its radius of curvature and in which the vorticity is everywhere the same is obtained, viscous and compressibility effects being neglected. The method of matched asymptotic expansions is applied and results in an extension, which is unexpectedly simple, of Birkhoff's integro-differential equation for vortex sheets. The equation is applied to long waves on a straight vortex layer of uniform thickness.     

A parallel computational scheme with ninth-order multioperator approximations and its application to direct numerical simulation

The properties of ninth-order multioperator compact schemes based on known third-and fifth-order compact approximations are examined. The domains where the multioperators have fixed signs are determined numerically. The numerical results are compared with the exact solution to the Burgers equation. The multioperator schemes are applied to the problem of vortex sheet roll-up.     

Analysis of quasi-steady acceleration in circular Couette flow

The effect of inner cylinder acceleration in the transition from circular Couette flow to Taylor vortex regime was numerically investigated. The solution of the eigenvalue problem when the Couette flow was perturbed allowed the determination of the relationship between the inner cylinder acceleration rate and the acceleration time duration expressed in terms of viscous diffusion time. The analysis of the numerical results also allowed the concept of quasi-steady acceleration in the Couette system to be qualitatively described.     

Simulation of particle dynamics in a viscous fluid near a plane wall

A method was developed (and implemented as a software package) for the numerical simulation of the dynamics of a large number of particles in a viscous flow near a plane wall. The fluid motion is considered in the Stokes approximation with allowance for the hydrodynamic interaction of all the particles. The velocities and trajectories of the particles are computed, and fluid streamlines are constructed. In the case of even two particles falling onto the wall, it is found that a complex vortex flow is formed between them. This gives rise to the particles having a velocity component along the wall. The dynamics of a group of several dozen particles moving parallel and perpendicular to the wall are simulated.     

A fourth-order approximate projection method for the incompressible Navier–Stokes equations on locally-refined periodic domains

In this follow-up of our previous work [30], the author proposes a high-order semi-implicit method for numerically solving the incompressible Navier–Stokes equations on locally-refined periodic domains. Fourth-order finite-volume stencils are employed for spatially discretizing various operators in the context of structured adaptive mesh refinement (AMR). Time integration adopts a fourth-order, semi-implicit, additive Runge–Kutta method to treat the non-stiff convection term explicitly and the stiff diffusion term implicitly. The divergence-free condition is fulfilled by an approximate projection operator. Altogether, these components yield a simple algorithm for simulating incompressible viscous flows on periodic domains with fourth-order accuracies both in time and in space. Results of numerical tests show that the proposed method is superior to previous second-order methods in terms of accuracy and efficiency. A major contribution of this work is the analysis of a fourth-order approximate projection operator.     

The Velocity of Viscous Vortex Rings

The motion of vortex rings of small cross section is considered. A formula is derived for the velocity of a ring in an ideal fluid with an arbitrary distribution of vorticity in the core. The effects of viscosity are then examined. A definition of the velocity of an unsteady diffusing ring is given and it is shown that the method used to calculate the ring speed in an ideal fluid can be extended to give the velocity of a vortex ring subject to viscous diffusion.     

The limiting vortex in the similarity solution for a swirling flow

It is known that the similarity solution for a viscous swirling flow over a stationary disk does not exist if the driving vortex far away from the disk is a potential vortex, while the solution exists for a rigid body vortex. Previously, the breakdown has been determined to occur if the azimuthal velocity of the driving vortex decreases faster than a certain power of the radial distance from the axis of symmetry. The decay parameter at which the similarity solution ceases to exist is computed here by a more direct method, and the reason for the breakdown becomes apparent. The analysis confirms (and slightly improves) the known value of the parameter. The case where the fluid, now assumed to be conducting, is subject to an axial magnetic field and the asymptotic behavior of the solution far away from the axis are also briefly discussed.     

Variable-Elliptic-Vortex Method for Incompressible Flow Simulation

A variable-elliptic-vortex method, which is a generalization of the elliptic-vortex method proposed by the author in [1], is presented for the numerical simulation of incompressible flows. The most attractive feature of the new method is that the numerical vortex blobs used in this model like actual vortex blobs can be translated, rotated and deformed in elliptic shape. The new method provides a more reasonable and more accurate approach for flow simulation than the fixed-vortex methods. Numerical examples are presented to demonstrate the performance of the new method.     

An application of a boundary element method to natural convection

The direct boundary element method is applied to the numerical modelling of thermal fluid flow in a transient state. The Navier-Stokes equations are considered under the Boussinesq approximation and the viscous thermal flow equations are expressed in terms of stream function, vorticity, and temperature in two dimensions. Boundary integral equations are derived using logarithmic potential and time-dependent heat potential as fundamental solutions. Boundary unknowns are discretized by linear boundary elements and flow domains are divided into a series of triangular cells. Charged points are translated upstream in the numerical evaluation of convective terms. Unknown stream function, vorticity, and temperature are staggered in the computational scheme.

Simple iteration is found to converge to the quasi steady-state flow. Boundary solutions for two-dimensional examples at a Reynolds number 100 and Grashoff number 10⁷ are obtained.     

涡旋演化的小波自适应模拟

该文考察了小波自适应方法用于涡旋运动的演化过程．首先,通过两个初边值问题,说明小波方法具有可精度可控和局部结构自动捕捉的能力．然后,计算了涡旋的合并过程,结果表明,小波方法可以准确高效的应用于流动涡旋的演化预测,进而,讨论了小波方法在湍流数值模拟中的应用．     

Poche de tourbillon pour Euler 2D dans un ouvert à bord

We consider Euler equations for an homogeneous incompressible non viscous fluid inside a smooth bounded domain of the plane. For an initial data of smooth vortex patch type, we obtain existence and uniqueness of a solution of the same type, locally in time if the initial patch is tangent to the boundary of the domain, and globally in time if the patch is far away from the boundary. We use pseudo-differential calculus to take care of the boundary condition. For the tangent limit case, we show that the velocity gradient of a vortex patch is Hölder continuous up to the boundary of the patch, using singular integrals. Our method provide also a result for several mutually tangent vortex patches in the plane.     

Streamwise oscillations of a cylinder beneath a free surface: Free surface effects on vortex formation modes

A computational study of a viscous incompressible two-fluid model with an oscillating cylinder is investigated at a Reynolds number of 200 and at a dimensionless displacement amplitude of A=0.13 and for the dimensionless forcing cylinder oscillation frequency-to-natural vortex shedding frequency ratios, f/f₀=1.5,2.5,3.5. Specifically, two-dimensional flow past a circular cylinder subject to forced in-line oscillations beneath a free surface is considered. The method is based on a finite volume discretization of the two-dimensional continuity and unsteady Navier-Stokes equations (when a solid body is present) on a fixed Cartesian grid. Two-fluid model based on improved volume-of-fluid method is used to discretize the free surface interface. The study focuses on the laminar asymmetric flow structure in the near wake region and lock-on phenomena at a Froude number of 0.2 and for the dimensionless cylinder submergence depths, h=0.25, 0.5 and 0.75. The equivorticity patterns and pressure distribution contours are used for the numerical flow visualization. The code validations in special cases show good comparisons with previous numerical results.     

Homogenization of the Equations of Filtration of a Viscous Fluid in Two Porous Media

A homogenized model of filtration of a viscous fluid in two domains with common boundary is deduced on the basis of the method of two-scale convergence. The domains represent an elastic medium with perforated pores. The fluid, filling the pores, is the same in both domains, and the properties of the solid skeleton are distinct.     

An estimate of the radius of the viscous core of a vortex

A formula is derived which relates the radius of a vortex in a viscous fluid to its kinetic energy of turbulence.     

The Existence of Görtler Vortices in Separated Boundary Layers

The linear stability properties of Görtler vortices within a general separated boundary layer flow are addressed. There has been little previous theoretical work directed toward this problem and here we are able to characterize the important features of vortices over the complete wavenumber spectrum. This investigation complements earlier studies of vortices within an attached flow which demonstrated that there are three distinctive wavenumber régimes which together describe the most relevant possibilities for vortex behavior. In the first of these, at relatively small wavenumbers, the mode is inviscid in character; as the vortex wavenumber increases so the spatial amplification rate of the vortices increases until viscous effects become significant and the growth rate begins to diminish. As the wavenumber increases yet further so the vortex is completely stabilized. Here we discuss the corresponding structures which may exist within a separated flow and the most significant result we find is that the maximum growth rate of a mode in this type of flow is actually greater than that which occurs when the flow has not separated. In addition, the inviscid modes are shown to have a far more complicated configuration than within an attached boundary layer and, indeed, their structure can only be completely determined by implementation of numerical procedures.     

A new stabilizing technique for boundary integral methods for water waves

Boundary integral methods to compute interfacial flows are very sensitive to numerical instabilities. A previous stability analysis by Beale, Hou and Lowengrub reveals that a very delicate balance among terms with singular integrals and derivatives must be preserved at the discrete level in order to maintain numerical stability. Such balance can be preserved by applying suitable numerical filtering at certain places of the discretization. While this filtering technique is effective for two-dimensional (2-D) periodic fluid interfaces, it does not apply to nonperiodic fluid interfaces. Moreover, using the filtering technique alone does not seem to be sufficient to stabilize 3-D fluid interfaces.

Here we introduce a new stabilizing technique for boundary integral methods for water waves which applies to nonperiodic and 3-D interfaces. A stabilizing term is added to the boundary integral method which exactly cancels the destabilizing term produced by the point vortex method approximation to the leading order. This modified boundary integral method still has the same order of accuracy as the point vortex method. A detailed stability analysis is presented for the point vortex method for 2-D water waves. The effect of various stabilizing terms is illustrated through careful numerical experiments.