非线性流动稳定性理论

<正> 流体的湍流运动普遍存在于大气、海洋、飞行器周围、推进装置和流体机械中,探讨湍流形成的条件和过程是流动稳定性理论的研究对象。 早在1843年,Stokes就预见到流体状态转捩的原因是失稳。整整100年前(1883年),Reynolds让液体流入不同口径的圆管,并在对称轴上注入一股纤细的颜色水,以便明显、敏感地观察到湍流的发生。他发现: (1)流速较低时,可以看到层次分明的层流,流速增加到一定程度后,就转变成高度无序的湍流状态;

一个推广Virial定理的数学方法和应用

本文在提出广义Gauss定理的思路下,提出了广义Stokes定理和一个推广Virial定理的新方法,后者应用于天体磁流体力学和引力平衡问题时得到的结果有:(1)气体具有运动时的平衡系统的判据,(2)磁场对气团形态的影响。应用于实验室等离子体平衡问题时,其结果有:(1)包围在气体中孤立磁场的特性,(2)发现内包无力场必须外包一有力场,(3)无力场的形态,应用于Tokmak等离子体环时的结果有:(1)环的胖瘦对环表面磁压的影响,(2)两个外加磁场分量分别和气压、环的胖瘦,截面形态、环电流分布和逆磁或顺磁的平衡关系,搞清楚了外加磁场约束等离子体总体平衡的物理机制。

探测气体平衡及非平衡态的相干反斯托克斯喇曼散射技术

<正> 一、引言 这篇综述的目的在于介绍最近发展起来的相干反斯托克斯喇曼散射(Coherent Anti-Stokes Raman Scattering(CARS))方法在诊断平衡及非平衡态气体中的应用,包括实验及理论两部分.文内还扼要介绍了我们正在筹备的测试系统的方案。 在气体中,特别在稀薄气体非平衡流动中,探测粒子数密度分布;混合含气中各组分,各能态的分布;平动、振动、转动温度;粒子的激发和弛豫过程,都是比较困难的,而这些

稀薄气流在麦克斯韦类型边界条件下对球的绕流

本文讨论低速稀薄气流(K_n《1)中圆球阻力问题,在Knudsen层中取分子在表面作Maxwell类型反射的壁面条件求解线性B-K-W方程,而以Stokes区的匹配渐近展开解为无穷远处的边界条件,在表面与均匀来流间温差为小量的假设下得到了小K_n数下圆球绕流的解。除稀薄气体效应外,还考虑了惯性效应和分子在表面的反射机制,所得的球阻表达为K_n、Re和热适应系数α的函数。

Effect of nonlinear wave-current interaction on flow fields and hydrodynamic forces

A fifth-order theory for solving the problem of interaction between Stokes waves and exponential profile currents is proposed. The calculated flow fields are compared with measurements. Then the errors caused by the linear superposition method and approximate theory are discussed. It is found that the total wave-current field consists of pure wave, pure current and interaction components. The shear current not only directly changes the flow field, but also indirectly does sx, by changing the wave parameters due to wave-current interaction. The present theory can predict the wave kinematics on shear currents satisfactorily. The linear superposition method may give rise to more than 40% loading error in extreme conditions. When the apparent wave period is used and the Wheeler stretching method is adopted to extrapolate the current, application of the approximate theory is the best.

Numerical simulation of coherent structure in two-dimensional compressible mixing layers

The coherent structure in two-dimensional mixing layers is simulated numerically with the compressible Navier-Stokes equations. The Navier-Stokes equations are discretized with high-order accurate upwind compact schemes. The process of development of flow structure is presented: loss of stability, development of Kelvin-Helmholtz instability, rolling up and pairing. The time and space development of the plane mixing layer and influence of the compressibility are investigated.

Numerical study of bifurcation solutions of spherical Taylor-Couette flow

The steady bifurcation flows in a spherical gap (gap ratio sigma=0.18) with rotating inner and stationary outer spheres are simulated numerically for Re(c1)less than or equal to Re less than or equal to 1 500 by solving steady axisymmetric incompressible Navier-Stokes equations using a finite difference method. The simulation shows that there exist two steady stable flows with 1 or 2 vortices per hemisphere for 775 less than or equal to Re less than or equal to 1 220 and three steady stable flows with 0, 1, or 2 vortices for 1 220

Numerical simulation of axisymmetric unsteady incompressible flow by a vorticity-velocity method

A new numerical method for solving the axisymmetric unsteady incompressible Navier-Stokes equations using vorticity-velocity variables and a staggered grid is presented. The solution is advanced in time with an explicit two-stage Runge-Kutta method. At each stage a vector Poisson equation for velocity is solved. Some important aspects of staggering of the variable location, divergence-free correction to the velocity held by means of a suitably chosen scalar potential and numerical treatment of the vorticity boundary condition are examined. The axisymmetric spherical Couette flow between two concentric differentially rotating spheres is computed as an initial value problem. Comparison of the computational results using a staggered grid with those using a non-staggered grid shows that the staggered grid is superior to the non-staggered grid. The computed scenario of the transition from zero-vortex to two-vortex flow at moderate Reynolds number agrees with that simulated using a pseudospectral method, thus validating the temporal accuracy of our method.

Triple-region structure for turbulent-flow in a square duct - a finite-element approach

The Reynolds-averaged Navier-Stokes equations for describing the turbulent flow in a straight square duct are formulated with two different turbulence models. The governing equations are then expanded as a multi-deck structure in a plane perpendicular to the streamwise direction, with each deck characterized by its dominant physical forces as commonly carried out in analytical work using triple-deck expansion. The resulting equations are numerically integrated using higher polynomial (H-P) finite element technique for each cross-sectional plane to be followed by finite difference representation in the streamwise direction until a fully developed state is reached. The computed results using the two different turbulence models show fair agreement with each other, and concur with the vast body of available experimental data. There is also general agreement between our results and the recent numerical works anisotropic k-epsilon turbulence model.

Computation of turbulent-flow in a square duct - aspects of the secondary flow and its origin

A numerical study of turbulent flow in a straight duct of square cross-section is made. An order-of-magnitude analysis of the 3-D, time-averaged Navier-Stokes equations resulted in a parabolic form of the Navier-Stokes equations. The governing equations, expressed in terms of a new vector-potential formulation, are expanded as a multi-deck structure with each deck characterized by its dominant physical forces. The resulting equations are solved using a finite-element approach with a bicubic element representation on each cross-sectional plane. The numerical integration along the streamwise direction is carried out with finite-difference approximations until a fully-developed state is reached. The computed results agree well with other numerical studies and compare very favorably with the available experimental data. One important outcome of the current investigation is the interpretation analytically that the driving force of the secondary flow in a square duct comes mainly from the second-order terms of the difference in the gradients of the normal and transverse Reynolds stresses in the axial vorticity equation.

Upwind compact method and direct numerical-simulation of driven flow in a square cavity

A high-order accurate finite-difference scheme, the upwind compact method, is proposed. The 2-D unsteady incompressible Navier-Stokes equations are solved in primitive variables. The nonlinear convection terms in the governing equations are approximated by using upwind biased compact difference, and other spatial derivative terms are discretized by using the fourth-order compact difference. The upwind compact method is used to solve the driven flow in a square cavity. Solutions are obtained for Reynolds numbers as high as 10000. When Re less than or equal to 5000, the results agree well with those in literature. When Re = 7500 and Re = 10000, there is no convergence to a steady laminar solution, and the flow becomes unsteady and periodic.

Approximate theory of natural-convection with small grashof number

To gain some insight into the behaviour of low-gravity flows in the material processing in space, an approximate theory has been developed for the convective motion of fluids with a small Grashof number Gr. The expansion of the variables into a series of Gr reduces the Boussinesq equation to a system of weakly coupled linearly inhomogeneous equations. Moreover, the analogy concept is proposed and utilized in the study of the plate bending problems in solid mechanics. Two examples are investigated in detail, i. e. the 2-dimensional steady flows in either circular or square infinite closed cylinder, which is horizontally imposed at a specified temperature of linear distribution on the boundaries. The results for stream function ψ, velocity u and temperature T are provided. The analysis of the influences of some parameters such as the Grashof number Gr and the Prandtl number Pr, on motions will lead to several interesting conclusions. The theory seems to be useful for seeking for an analytical solutions. At least, it will greatly simplify the complicated problems originally governed by the Navier-Stokes equation including buoyancy. It is our hope that the theory might be applicable to unsteady or 3-dimensional cases in future.

Compressible laminar boundary-layer flows of a dusty gas over a semi-infinite flat-plate

The compressible laminar boundary-layer flows of a dilute gas-particle mixture over a semi-infinite flat plate are investigated analytically. The governing equations are presented in a general form where more reasonable relations for the two-phase interaction and the gas viscosity are included. The detailed flow structures of the gas and particle phases are given in three distinct regions : the large-slip region near the leading edge, the moderate-slip region and the small-slip region far downstream. The asymptotic solutions for the two limiting regions are obtained by using a seriesexpansion method. The finite-difference solutions along the whole length of the plate are obtained by using implicit four-point and six-point schemes. The results from these two methods are compared and very good agreement is achieved. The characteristic quantities of the boundary layer are calculated and the effects on the flow produced by the particles are discussed. It is found that in the case of laminar boundary-layer flows, the skin friction and wall heat-transfer are higher and the displacement thickness is lower than in the pure-gas case alone. The results indicate that the Stokes-interaction relation is reasonable qualitatively but not correct quantitatively and a relevant non-Stokes relation of the interaction between the two phases should be specified when the particle Reynolds number is higher than unity.

Theory of local mode excitation in polyatomics by frequency-modulated lasers

that the Stokes-interaction relation is reasonable qualitatively but not correct

Characteristics of a CW HF chemical-laser calculated from a simplified two-dimensional model

A two-dimensional simplified model of an HF chemical laser is introduced. Using an implicit finite difference scheme, the solution of two adjacent parallel streams with diffusion mixing and chemical reaction is generated. A contour of mixing and reaction boundary is obtained without presupposition. The distribution of the HF(v) concentrations, gas temperature and the optical small signal gain (alpha sub V, J) on the flowing plane (X, Y) are presented. Compared with the solution solved directly from a set of Navier-Stokes equations, the results of these two methods agree with each other qualitatively. The influences of the different velocity, temperature (T sub 0) and composition of the two streams on the small signal gain after the nozzle exit are investigated. It is interesting that for larger J with a fixed v, the peaks of alpha sub v-T sub 0 profiles move towards higher T sub 0. The computing method is simple and only a short computing time is needed.