带有横向射流的三维超声速湍流流场分析

应用显式的五阶WENO格式,结合k-ε湍流模型,求解三维Favre平均N-S方程,计算了从方孔横向喷出的声速气流与马赫数为30的超声速气流的干扰流场。结果表明,在射流上游,射流的阻碍使超声速气流产生分离,形成两个主要的回流区域,主回流导致在方孔射流两侧形成马蹄涡区域,射流下游存在低压区域,形成较小的回流以及一对螺旋形旋涡。     

The critical cross-section of a vortex

Summary This paper discusses the problem of critical-flow cross-sections in vortex flows. It is shown that there are two different types of vortex flows, A-type and B-type vortices (say). An A-type vortex approaches its critical flow state as its cross-sectional area increases and departs from the critical state as the cross-sectional area is decreased. This property is associated with the particular dependence of total pressure and circulation on the stream function, and it holds for both subcritical and supercritical A-type vortices. On the other hand, both subcritical and supercritical B-type vortices approach their critical flow states as their cross-sectional areas are decreased and depart from their critical states for increasing cross-sectional area. As was shown by Benjamin, setting the first variation of the flow force with respect to the stream function equal to zero leads to Euler's equation of motion. The second variation also vanishes if the corresponding flow state is critical. In this case the sign of the third variation decides whether the flow is an A-type or a B-type vortex. Within the framework of inviscid-fluid flow theory the type of a vortex is preserved unless vortex breakdown occurs. Making use of the knowledge that vortex flows are controlled by two different types of critical-flow cross-sections a variety of vortex flow phenomena are investigated, including the two types of inlet vortices that are observed upstream of jet engines, the behavior of vortex valves, the flow characteristics of liquid-fuel atomizers and the bath tub vortex.     

A Class of Solutions for the Equations Governing the Steady Two-Dimensional Motion of a Viscous Incompressible Liquid

A general viscosity dependent solution for the stream function is found satisfying the simplest nondegenerate form of the steady flow Navier-Stokes equations for a viscous incompressible liquid. The solution is two-dimensional and is expressed in terms of arbitrary analytic functions in the fluid domain. This class of flows is generated by complex stream functions, and the region of definition is restricted by an inequality containing these analytic functions. A general potential flow, and degenerate Stokes or creeping flows are recovered as particular solutions in limiting cases.     

Exact solutions of problems of the three-dimensional flow of ideal viscous incompressible fluids near cylindrical surfaces

Three-dimensional flows of an incompressible fluid, the parameters of which depend on two coordinates and time, are considered. The stream surfaces of such flows are cylindrical. The equations of continuity and the Navier-Stokes equations can be transformed to relations, one of which is the equation for the stream function the other is the integral of the equations relating the pressure and the stream function, and the third is a linear equation for the projection of the velocity vector onto the axis parallel to the generatrix of the cylindrical surfaces. The problems of modelling the flows are considered on the basis of the exact solutions of the Navier-Stokes equations and Euler's equations using examples. Relations for the distribution of the flow parameters in the channel created by hyperbolical cylinders are derived for the case of unsteady inviscid flow. The streamlines of these flows are situated on the side surfaces of the hyperbolical cylinders and intercept the generatrices of the cylinders at certain indirect angles. The flow around a circular cylinder and the flow of fluid inside an elliptic cylinder are considered in the case of steady inviscid flow. The streamlines on the circular cylinder are arranged transverse to the cylinder (the projection of the velocity vector onto the coordinate axis, parallel to the generatrix of the cylinder, is equal to zero). Far from the cylinder the streamlines are also situated on a cylindrical surfaces, but not transverse to the cylinder, making certain indirect angles with the generatrix. Viscous three-dimensional flows, possessing a certain symmetry, are considered. In the case of radial symmetry the streamlines are helical lines. The non-planar Couette flow between parallel moving planes is characterized by the fact that the velocity vectors, being situated in the same plane, are collinear, while the velocity vectors in parallel planes are not collinear. Relations for viscous steady three-dimensional flows, using well-known relations, obtained for the stream function of two-dimensional flows, are given.     

Vortex breakdown in turbulent pipe flows

In many technical applications turbulent flows with embedded slender vortices exist. Depending on the boundary conditions vortex breakdown can occur. The purpose of this work is to develop and implement a solution scheme for large‐eddy simulations of vortex breakdown in turbulent pipe flows. One of the main problems in this simulation is the formulation of the inflow boundary condition for a fully developed turbulent flow with an embedded vortex. For that purpose a rescaling technique is developed in which a solution at a downstream location is inserted at the inflow boundary after an appropriate rescaling. To determine rescaling laws for pipe flows with an embedded vortex, analytical velocity profiles of swirling flows are first prescribed in a laminar flow. From the spatial development of the vortex a scaling law is deduced. In a next step this procedure is to be transferred to turbulent flows.     

The optimal convergence rates of non-isentropic subsonic Euler flows through the infinitely long three-dimensional axisymmetric nozzles

This paper is concerned about the optimal convergence rates of non-isentropic subsonic flows at far fields in three-dimensional infinitely long axisymmetric nozzles. By using the stream function formulation for the compressible Euler equations, the subsonic Euler flows are equivalent to a quasilinear elliptic equation of the stream function. The key points to prove the convergence rates of subsonic flows at far fields are the choice of compared functions and the maximum principles.     

Algorithm for solving the Navier–Stokes equations for the modeling of creeping flows

We study a coupled algorithm for solving the two-dimensional Navier–Stokes equations in the stream function–vorticity variables. The algorithm is based on a finite-difference scheme in which the inertial terms in the vortex transport equation are taken from the lower time layer and the dissipative terms, from the upper time layer. In the linear approximation, we study the stability of this algorithm and use test computations to show its advantages when modeling flows at moderate Reynolds numbers.     

Tollmien-Schlichting/vortex interactions in compressible boundary-layer flows

A weakly nonlinear interaction of oblique Tollmien-Schlichtingwaves and longitudinal vortices in compressible, high Reynoldsnumber, boundary-layer flow over a flat plate is consideredfor all ranges of the Mach number. The interaction equationsconsist of equations for the vortex which is indirectly forcedby the waves via a boundary condition, whereas a vortex termappears in the amplitude equation for the wave pressure. Thedownstream solution properties of interaction equations arefound to depend on the sign of an interaction coefficient. Thisparticular type of weakly nonlinear interaction was first proposedby Hall & Smith (1989), who considered incompressible flows;however, there are some errors in their formulation. Correctedresults for the incompressible regime are presented for comparisonwith those calculated for compressible flows. Compressibilityis found to have a significant effect on the interaction properties,principally through its impact on the waves and their governingmechanism, the ‘triple-deck’ structure. It is foundthat, in general, the flow quantities will grow slowly withincreasing downstream coordinate. However, for flows with Machnumber values below 2, there exists a small band of wave obliquenessangles for which the solutions terminate in abrupt, finite-distance‘break-ups’.     

A lattice Boltzmann model for the eddy–stream equations in two-dimensional incompressible flows

A lattice Boltzmann model for two-dimensional incompressible flows with eddy–stream equations is proposed. By using two kinds of distribution functions and employing several higher-order moments of equilibrium distribution functions, the eddy equation and stream function equation with the second-order truncation error are obtained. In the numerical examples, we compared the numerical results of this scheme with those obtained by other classical method. The numerical results agree well with the classical ones.     

Three-Dimensional Nonlinear Dispersive Waves on Shear Flows

The Green–Naghdi equations describing three-dimensional water waves are considered. Assuming that transverse variations of the flow occur at a much shorter lengthscale than variations along the wave propagation direction, we derive simplified asymptotic equations from the Green–Naghdi model. For steady flows, we show that the approximate model reduces to a one-dimensional Hamiltonian system along each stream line. Exact solutions describing a wide class of free-boundary flows depending on several arbitrary functions of one argument are found. The numerical results showing different patterns of steady three-dimensional waves are presented.     

Three-dimensional sand ripples as the product of vortex instability

Three-dimensional sand ripples can be observed under steady liquid flows in both nature and industry. Some examples are the ripples observed on the bed of rivers and in petroleum pipelines conveying sand. Although of importance, the formation of these patterns is not completely understood. There are theoretical and experimental evidence that aquatic ripples grow from two-dimensional bed instabilities, so that a straight vortex is formed just downstream of their crests. The proposition of Raudkivi (2006) [18], that three-dimensionality has its origin in a vortex instability, is employed here. This paper presents a linear stability analysis of the downstream vortex in order to obtain the transverse scales of three-dimensional ripples. The obtained wavelength is compared with experimentally observed ripples.     

Numerical analysis of the dynamics of distributed vortex configurations

A numerical algorithm is proposed for analyzing the dynamics of distributed plane vortex configurations in an inviscid incompressible fluid. At every time step, the algorithm involves the computation of unsteady vortex flows, an analysis of the configuration structure with the help of heuristic criteria, the visualization of the distribution of marked particles and vorticity, the construction of streamlines of fluid particles, and the computation of the field of local Lyapunov exponents. The inviscid incompressible fluid dynamic equations are solved by applying a meshless vortex method. The algorithm is used to investigate the interaction of two and three identical distributed vortices with various initial positions in the flow region with and without the Coriolis force.     

Point vortices and classical orthogonal polynomials

Stationary equilibria of point vortices in the plane and on the cylinder in the presence of a background flow are studied. Vortex systems with an arbitrary choice of circulations are considered. Differential equations satisfied by generating polynomials of vortex configurations are derived. It is shown that these equations can be reduced to a single one. It is found that polynomials that are Wronskians of classical orthogonal polynomials solve the latter equation. As a consequence vortex equilibria at a certain choice of background flows can be described with the help of Wronskians of classical orthogonal polynomials.     

在微重力条件下无限长方柱中的自然对流*

为了研究在宇宙空间微重力环境中.自然对流对流体运动的影响,将变量展成Grashof的摄动级数,使用摄动理论将Navier-Stokes方程组简化成:关于温度T的Poisson方程,关于流函数ψ的非齐次biharmonic方程.选取一无限长封闭方柱体,假定在柱体边界上预先给定一种线性温度分布,使用数值计算方法求解上述简化方程组,得到各阶流函数和各阶温度值,进而详细地研究了方柱中流体的运动状况,分析和讨论了某些参数,如Grashof数和Prandtl数对流体运动的影响,最后将计算结果与由未简化方程推算的结果进行比较,证实近似方法正确地简化了复杂的流体运动过程,并且可以推广、运用到三维问题上.     

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.     

Equilibrium statistical theory for nearly parallel vortex filaments

The first mathematically rigorous equilibrium statistical theory for three‐dimensional vortex filaments is developed here in the context of the simplified asymptotic equations for nearly parallel vortex filaments, which have been derived recently by Klein, Majda, and Damodaran. These simplified equations arise from a systematic asymptotic expansion of the Navier‐Stokes equation and involve the motion of families of curves, representing the vortex filaments, under linearized self‐induction and mutual potential vortex interaction. We consider here the equilibrium statistical mechanics of arbitrarily large numbers of nearly parallel filaments with equal circulations. First, the equilibrium Gibbs ensemble is written down exactly through function space integrals; then a suitably scaled mean field statistical theory is developed in the limit of infinitely many interacting filaments. The mean field equations involve a novel Hartree‐like problem with a two‐body logarithmic interaction potential and an inverse temperature given by the normalized length of the filaments. We analyze the mean field problem and show various equivalent variational formulations of it. The mean field statistical theory for nearly parallel vortex filaments is compared and contrasted with the well‐known mean field statistical theory for two‐dimensional point vortices. The main ideas are first introduced through heuristic reasoning and then are confirmed by a mathematically rigorous analysis. A potential application of this statistical theory to rapidly rotating convection in geophysical flows is also discussed briefly. © 2000 John Wiley & Sons, Inc.     

Transonic nozzle flows and free boundary problems for the full Euler equations

We establish the existence and uniqueness of transonic flows with a transonic shock through a two-dimensional nozzle of slowly varying cross-sections. The transonic flow is governed by the steady, full Euler equations. Given an incoming smooth flow that is close to a constant supersonic state (i.e., smooth Cauchy data) at the entrance and the subsonic condition with nearly horizontal velocity at the exit of the nozzle, we prove that there exists a transonic flow whose downstream smooth subsonic region is separated by a smooth transonic shock from the upstream supersonic flow. This problem is approached by a one-phase free boundary problem in which the transonic shock is formulated as a free boundary. The full Euler equations are decomposed into an elliptic equation and a system of transport equations for the free boundary problem. An iteration scheme is developed and its fixed point is shown to exist, which is a solution of the free boundary problem, by combining some delicate estimates for the elliptic equation and the system of transport equations with the Schauder fixed point argument. The uniqueness of transonic nozzle flows is also established by employing the coordinate transformation of Euler-Lagrange type and detailed estimates of the solutions.     

Explicit Solutions of the Steady Two-Dimensional Navier-Stokes Equations

A method is described to determine exact general solutions of the steady two-dimensional Navier-Stokes equations. The solutions can also be utilized to construct the stream function for a class of unsteady flows, which in turn contain two arbitrary analytic functions of a complex variable defined in the fluid region.     

Analyzing a Two-Stage Queueing System with Many Point Process Arrivals at Upstream Queue

We consider a two-stage queueing system where the first (upstream) queue serves many flows, of which a fixed set of flows arrive to the second (downstream) queue. We show that as the capacity and the number of flows aggregated at the upstream queue increases, the overflow probability at the downstream queue converges to that of a simplified single queue obtained by removing the upstream queue from the original two-stage queueing system. Earlier work shows such convergence for fluid traffic, by exploiting the large deviation result that the workload goes to zero almost surely, as the number of flows and capacity is scaled. However, the analysis is quite different and more difficult for the point process traffic considered in this paper. The reason is that for point process traffic the large deviation rate function need not be strictly positive (i.e., I(0)=0), hence the workload at the upstream queue may not go to zero even though the number of flows and capacity go to infinity. The results in this paper thus make it possible to decompose the original two-stage queueing system into a simple single-stage queueing system.     

Possibility of explaining the existence of vortexlike hydrodynamic structures based on the theory of stationary kinetic equations

The possibility of describing vortex structures in quasi-one-dimensional plane flows by applying kinetic equations and bifurcation theory is examined. The Lyapunov-Schmidt method is used to obtain a system of Riccati-type generalized bifurcation equations. An analysis of its properties leads to conditions for the existence of vortex structures.