Features of the supercritical transition in the wake of a circular cylinder

The features of the wake behind a uniform circular cylinder at Re = 200, which is just beyond the critical Reynolds number of 3-D transition, are investigated in detail by direct numerical simulations by solving 3-D incompressible Navier-Stokes equations using mixed spectral-spectral-element method. The high-order splitting algorithm based on the mixed stiffly stable scheme is employed in the time discretization. Due to the nonlinear evolution of the secondary instability of the wake, the spanwise modes with different wavelengths emerge. The spanwise characteristic length determines the transition features and global properties of the wake. The existence of the spanwise phase difference of the primary vortices shedding is confirmed by Fourier analysis of the time series of the spanwise vorticity and attributed. to the dominant spanwise mode. The spatial energy distributions of various modes and the velocity profiles in the near wake are obtained. The numerical results indicate that the near wake is in 3-D quasi-periodic laminar state with transitional behaviors at this supercritical Reynolds number.

再入尾迹电磁特性的湍流效应研究

本文旨意在于通过探讨高超声速再入尾迹中的湍流等离子体与电磁波相互作用的机理，以及建立能反映此机理的应用性理论模型，从而提供一套可进行目标特性分析的方法，以便为工程部门的突防技术服务。本题目在再入气动物理现象研究中具有重要意义。综合分析指出，地面雷达观测到的非相干散射信号主要来源于再入尾迹的亚密湍流区产生的体积散射。因此，电磁散射特性分析主要针对尾迹亚密湍流等离子体。并且，这里所有的分析都是根据在工程应用中最成熟的一阶畸变波Born近似理论模型。再入尾迹电磁特性的湍流效应研究，着眼点就在于湍流等离子体场的研究。对湍流等离子体场理论模型，本文试图通过模式理论来表达，即求解平均化的全Navier-Stokes方程及其封闭方程k-ε-g模型，从而准确获得流动平均场和脉动场信息。这种表达方式较以前有了较大改进。注意到高超声速流动具有强烈可压缩性的特点，故使用的N-S平均方程由质量加权平均过程产生，湍流模型方程也经过可压缩性修正。方程的离散求解方法，都是运用带矢通量分裂的二阶TVD格式的有限体积法。再入尾迹湍流场的初始条件由近尾迹（底部）流动经N-S方程求解给定，初始值更加准确可靠。尾迹从层流到湍流的转捩过程采用相对成熟的半经验公式确定。飞行器的高超声速再入过程必然导致它周围的空气温度升高，使得流动表现出真实气体效应。对重点考察的湍流流动而言，真实气体效应主要表现为气体处于热化学平衡状态。就工程部门面临的实际问题，把一阶畸变波Born近似的解算方法做些改进，使其能够处理的范围从轴对称尾迹扩展到三维湍流等离子体场是必要的。这为深入的理论分析提供了有力的保障。在能够准确模拟湍流流动的刻划雷达散射截面的基础上，考察亚密湍流等离子体对电磁散射的影响。通过选择的几个有代表性的因素进行讨论，初步结果表明：湍流转捩方式、湍流尺度对尾迹雷达散射截面值计算影响不大，而电子组份脉动能初始值影响较明显，且在特定条件下湍流模型的影响亦不大。但由于湍流模型涉及脉动初始值，其影响需进一步确定。同时，一些今后开展继续此项研究工作的有益建议也提了出来。

Generation of long, laminar plasma jets at atmospheric pressure and effects of flow turbulence

Long, laminar plasma jets at atmospheric pressure of pure argon and a mixture of argon and nitrogen with jet length up to 45 fi,Hes its diameter could be generated with a DC are torch by! restricting the movement of arc root in the torch channel. Effects of torch structure, gas feeding, and characteristics of power supply on the length of plasma jets were experimentally examined. Plasma jets of considerable length and excellent stability could be obtained by regulating the generating parameters, including are channel geometry gas flow I ate, and feeding methods, etc. Influence of flow turbulence at the torch,nozzle exit on the temperature distribution of plasma jets was numerically simulated. The analysis indicated that laminar flow plasma with very low initial turbulent kinetic energy will produce a long jet, with low axial temperature gradient. This kind of long laminar plasma jet could greatly improve the controllability for materials processing, compared with a short turbulent are let.

Starting Nozzle Flow Simulation Using K-G Two-Equation Turbulence Model

The starting process of two-dimensional nozzle flows has been simulated with Euler, laminar and k - g two-equation turbulence Navier-Stokes equations. The flow solver is based on a combination of LUSGS subiteration implicit method and five spatial discretized schemes, which are Roe, HLLE, MHLLE upwind schemes and AUSM+, AUSMPW schemes. In the paper, special attention is for the flow differences of the nozzle starting process obtained from different governing equations and different schemes. Two nozzle flows, previously investigated experimentally and numerically by other researchers, are chosen as our examples. The calculated results indicate the carbuncle phenomenon and unphysical oscillations appear more or less near a wall or behind strong shock wave except using HLLE scheme, and these unphysical phenomena become more seriously with the increase of Mach number. Comparing the turbulence calculation, inviscid solution cannot simulate the wall flow separation and the laminar solution shows some different flow characteristics in the regions of flow separation and near wall.

Direct numerical simulation of transition and turbulence in compressible mixing layer

The three-dimensional compressible Navier-Stokes equations are approximated by a fifth order upwind compact and a sixth order symmetrical compact difference relations combined with three-stage Ronge-Kutta method. The computed results are presented for convective Mach number Mc = 0.8 and Re = 200 with initial data which have equal and opposite oblique waves. From the computed results we can see the variation of coherent structures with time integration and full process of instability, formation of Lambda-vortices, double horseshoe vortices and mushroom structures. The large structures break into small and smaller vortex structures. Finally, the movement of small structure becomes dominant, and flow field turns into turbulence. It is noted that production of small vortex structures is combined with turning of symmetrical structures to unsymmetrical ones. It is shown in the present computation that the flow field turns into turbulence directly from initial instability and there is not vortex pairing in process of transition. It means that for large convective Mach number the transition mechanism for compressible mixing layer differs from that in incompressible mixing layer.

Turbulence in laminar premixed V-flames

Strong velocity fluctuations had been found in the laminar premixed V-flames. These velocity fluctuations are closely related to the chemical reaction. But the effects of the upstream combustible mixture velocity on the velocity fluctuations inside the flame are quite weak. The probability distribution function (PDF) of the velocity in the centre region of the flame appears "flat top" shaped. By analyzing the experiment results the flame-flow interactions are found to affect the flame not only at large scale in the flow field but also at small scale inside the flame. These effects will give rise to flame generated small scale turbulences.

Optimized Group Velocity Control Scheme and DNS of Decaying Compressible Turbulence of Relative High Turbulent Mach Number

To overcome the difficulty in the DNS of compressible turbulence at high turbulent Mach number, a new difference scheme called GVC8 is developed. We have succeeded in the direct numerical simulation of decaying compressible turbulence up to turbulent Mach number 0.95. The statistical quantities thus obtained at lower turbulent Mach number agree well with those from previous authors with the same initial conditions, but they are limited to simulate at lower turbulent Mach numbers due to the so-called start-up problem. The energy spectrum and coherent structure of compressible turbulent flow are analysed. The scaling law of compressible turbulence is studied. The computed results indicate that the extended self-similarity holds in decaying compressible turbulence despite the occurrence of shocklets, and compressibility has little effects on relative scaling exponents when turbulent Mach number is not very high.

Influence of Obstacle-Produced Turbulence on Development of Premixed Flames

An investigation into influence of obstructions on premixed flame propagation has been carried out in a semi-open tube. It is found that there exists flame acceleration and rising overpressure along the path of flame due to obstacles. According to the magnitude of flame speeds, the propagation of flame in the tube can be classified into three regimes: the quenching, the choking and the detonation regimes. In premixed flames near the flammability limits, the flame is observed first to accelerate and then to quench itself after propagating past a certain number of obstacles. In the choking regime, the maximum flame speeds are somewhat below the combustion product sound speeds, and insensitive to the blockage ratio. In the more sensitive mixtures, the transition to detonation (DDT) occurs when the equivalence ratio increases. The transition is not observed for the less sensitive mixtures. The dependence of overpressure on blockage ratio is not monotonous. Furthermore, a numerical study of flame acceleration and overpressure with the unsteady compressible flow model is performed, and the agreement between the simulation and measurements is good.

Simulation of Coal Combustion by AUSM Turbulence-Chemistry Char Combustion Model and a Full Two-Fluid Model

An algebraic unified second-order moment (AUSM) turbulence-chemistry model of char combustion is introduced in this paper, to calculate the effect of particle temperature fluctuation on char combustion. The AUSM model is used to simulate gas-particle flows, in coal combustion in a pulverized coal combustor, together with a full two-fluid model for reacting gas-particle flows and coal combustion, including the sub-models as the k-epsilon-k(p) two-phase turbulence niodel, the EBU-Arrhenius volatile and CO combustion model, and the six-flux radiation model. A new method for calculating particle mass flow rate is also used in this model to correct particle outflow rate and mass flow rate for inside sections, which can obey the principle of mass conservation for the particle phase and can also speed up the iterating convergence of the computation procedure effectively. The simulation results indicate that, the AUSM char combustion model is more preferable to the old char combustion model, since the later totally eliminate the influence of particle temperature fluctuation on char combustion rate.

Numerical Study Of Heat Transfer Mechanism In Turbulent Supercritical Co2 Channel Flow

Direct numerical simulation (DNS) of supercritical CO2 turbulent channel flow has been performed to investigate the heat transfer mechanism of supercritical fluid. In the present DNS, full compressible Navier-Stokes equations and Peng-Robison state equation are solved. Due to effects of the mean density variation in the wall normal direction, mean velocity in the cooling region becomes high compared with that in the heating region. The mean width between high-and low-speed streaks near the wall decreases in the cooling region, which means that turbulence in the cooling region is enhanced and lots of fine scale eddies are created due to the local high Reynolds number effects. From the turbulent kinetic energy budget, it is found that compressibility effects related with pressure fluctuation and dilatation of velocity fluctuation can be ignored even for supercritical condition. However, the effect of density fluctuation on turbulent kinetic energy cannot be ignored. In the cooling region, low kinematic viscosity and high thermal conductivity in the low speed streaks modify fine scale structure and turbulent transport of temperature, which results in high Nusselt number in the cooling condition of the supercritical CO2.

Oil-Water Two-Phase Flow Inside T-Junction

The oil/water two-phase flow inside T-junctions was numerically simulated with a 3-D two-fluid model, and the turbulence was described using the mixture k - epsilon model. Some experiments of oil/water flow inside a single T-junction were conducted in the laboratory. The results show that the separating performance of T-junction largely depends oil the inlet volumetric fraction and flow patterns. A reasonable agreement is reached between the numerical simulation and the experiments for both the oil fraction distribution and the separation efficiency.

Sediment incipience in turbulence generated in a square tank by a vertically oscillating grid

The failure of hydraulic structures in many estuaries and coastal regions around the world has been attributed to sediment transport and local scour. The sediment incipience in homogenous turbulence generated by oscillating grid is studied in this paper. The turbulent flow is measured by particle tracer velocimetry (PTV) technique. The integral length scale and time scale of turbulence are obtained. The turbulent flow near the wall is measured by local optical magnification. The sediment incipience is described by static theory. The relationship of probability of sediment incipience and the turbulent kinetic energy were obtained experimentally and theoretically. The distribution of the turbulent kinetic energy near the wall is found to obey the power law and the turbulent energy is further identified as the dynamic mechanism of sediment incipience.

Attractor crisis and bursting in a fluid flow with two no-slip directions

Characteristic burtsing behavior is observed in a driven, two-dimensional viscous flow, confined to a square domain and subject to no-slip boundaries. Passing a critical parameter value, an existing chaotic attractor undergoes a crisis, after which the flow initially enters a transient bursting regime. Bursting is caused by ejections from and return to a limited subdomain of the phase space, whereas the precrisis chaotic set forms the asymptotic attractor of the flow. For increasing values of the control parameter the length of the bursting regime increases progressively. Passing another critical parameter value, a second crisis leads to the appearance of a secondary type of bursting, of very large dynamical range. Within the bursting regime the flow then switches in irregular intervals from the primary to the secondary type of bursting. Peak enstrophy levels for both types of bursting are associated to the collapse of a primary vortex into a quadrupolar state.

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.

Some features of the turbulence produced by rising air bubbles in water

Presented in this paper is an experimental study on the characteristics of the turbulence produced by rising air bubbles in water. The measurements of turbulent velocities were made by using visualization technique of particle streak and computer image processing of the flow field. The turbulence features have been examined, showing that the rising bubble-produced turbulence can be approximately modeled by homogeneous turbulence as in the case of grid turbulence in air.