Direct Numerical Simulation of Passive Scalar in Decaying Compressible Turbulence

The passive scalars in the decaying compressible turbulence with the initial Reynolds number (defined by Taylor scale and RMS velocity) Re=72, the initial turbulent Mach numbers (defined by RMS velocity and mean sound speed) Mt=0.2-0.9, and the Schmidt numbers of passive scalar Sc=2-10 are numerically simulated by using a 7th order upwind difference scheme and 8th order group velocity control scheme. The computed results are validated with different numerical methods and different mesh sizes. The Batchelor scaling with k(-1) range is found in scalar spectra. The passive scalar spectra decay faster with the increasing turbulent Mach number. The extended self-similarity (ESS) is found in the passive scalar of compressible turbulence.

不可压N-S方程高效算法及二维槽道湍流分析

构造了基于非等距网格的迎风紧致格式,并将其与三阶精度的Adams半隐方法相结合,构造了求解不可压N-S方程高效算法。在数值模拟的基础上对二维槽道流动进行了分析,得到了Reynolds数从6 000到15 000的二维流动饱和态解(所谓"二维槽道湍流");定性及定量结果均与他人的数值计算结果吻合十分理想。对流场进行了统计分析,指出了"二维湍流"与三维湍流统计特性的区别。

Unsteady Model of Drop Marangoni Migration in Microgravity

The experiments of drop Marangoni migration have been performed by the drop shift facility of short period of 4.5 s, and the drop accelerates gradually to an asymptotic velocity during the free fall. The unsteady and axisymmetric model is developed to study the drop migration for the case of moderate Reynolds number Re = O(1), and the results are compared with the experimental ones in the present paper. Both numerical and experimental results show that the migration velocity for moderate Reynolds number is several times smaller than that given by the linear YGB theory.

Experimental Study on the Thermal Argon Plasma Generation and Jet Length Change Characteristics at Atmospheric Pressure

The generation, jet length and flow-regime change characteristics of argon plasma issuing into ambient air have been experimentally examined. Different torch structures have been used in the tests. Laminar plasma jets can be generated within a rather wide range of working-gas flow rates, and an unsteady transitional flow state exists between the laminar and turbulent flow regimes. The high-temperature region length of the laminar plasma jet can be over an order longer than that of the turbulent plasma jet and increases with increasing argon flow rate or arc current, while the jet length of the turbulent plasma is less influenced by the generating parameters. The flow field of the plasma jet has very high radial gradients of plasma parameters, and a Reynolds number alone calculated in the ordinary manner may not adequately serve as a criterion for transition. The laminar plasma jet can have a higher velocity than that of an unsteady or turbulent jet. The long laminar plasma jet has good stiffness to withstand the impact of laterally injected cold gas and particulate matter. It could be used as a rather ideal object for fundamental studies and be applied to novel materials processing due to its attractive stable and adjustable properties.

Planar Thermocapillary Migration of Two Bubbles in Microgravity Environment

A theoretical investigation is performed on the thermocapillary motion of two bubbles in arbitrary configuration in microgravity environment under the assumption that the surface tension is high enough to keep the bubbles spherical. The two bubbles are dr

Extended self-similar scaling law of multi-scale eddy structure in wall turbulence

The longitudinal fluctuating velocity of a turbulent boundary layer was measured in a water channel at a moderate Reynolds number. The extended self-similar scaling law of structure function proposed by Benzi was verified. The longitudinal fluctuating velocity, in the turbulent boundary layer was decomposed into many multi-scale eddy structures by wavelet transform. The extended self-similar scaling law of structure function for each scale eddy velocity was investigated. The conclusions are I) The statistical properties of turbulence could be self-similar not only at high Reynolds number, but also at moderate and low Reynolds number, and they could be characterized by the same set of scaling exponents xi (1)(n) = n/3 and xi (2)(n) = n/3 of the fully developed regime. 2) The range of scales where the extended self-similarity valid is much larger than the inertial range and extends far deep into the dissipation range,vith the same set of scaling exponents. 3) The extended selfsimilarity is applicable not only for homogeneous turbulence, but also for shear turbulence such as turbulent boundary layers.

Flow Characteristics of Liquids in Microtubes Driven by a High Pressure

The flow characteristics of liquids in microtubes driven by a high pressure ranging from 1 MPa to 30 MPa are studied in this paper. The diameter of the microtube is from 3 μm to 10 μm and liquids composed of simple small molecules are chosen as the working fluids. The Reynolds number ranges from 0. 1 to 24. The behavior of isopropanol and carbon tetrachloride under high pressure is found different from the prediction from conventional Hagen-Poiseuille (HP) equation. The normalized friction coefficient C* increases significantly with the pressure. From an analysis of the microtube deformation, liquid compressibility, viscous heating and wall slip, it may be seen that the viscosity at high pressure plays an important role here. An exponential function of viscosity vs pressure is introduced into the HP equation to counteract the difference between experimental and theoretical values. However, this difference is not so marked for di-water.

Calculations of longitudinal and transverse velocity structure functions using a vortex model of isotropic turbulence

The longitudinal structure function (LSF) and the transverse structure function (TSF) in isotropic turbulence are calculated using a vortex model. The vortex model is composed of the Rankine and Burgers vortices which have the exponential distributions in the vortex Reynolds number and vortex radii. This model exhibits a power law in the inertial range and satisfies the minimal condition of isotropy that the second-order exponent of the LSF in the inertial range is equal to that of the TSF. Also observed are differences between longitudinal and transverse structure functions caused by intermittency. These differences are related to their scaling differences which have been previously observed in experiments and numerical simulations.

高雷诺数流动的控制方程体系和扩散抛物化Navier-Stokes方程组的意义和用途

在计算机发达的时代,高雷诺(Re)数绕流计算中有无必要使用简化NS方程组,本文讨论这个问题.主要内容如下:(1)高Re数绕流包含3种基本流动:所有方向对流占优流动、所有方向对流扩散竞争流动和部分方向对流占优部分方向对流扩散竞争流动(简称干扰剪切流动),3个基本流动的特征彼此不同且在流场中所占领域大小彼此相差悬殊,NS方程区域很小,它们的最简单控制方程组Euler、Navier-Stokes(NS)和扩散抛物化(DP)NS方程组的数学性质彼此不同,因此利用Euler-DPNS-NS方程组体系分析计算高Re数绕流流动就是一个合乎逻辑的选择,该法与利用单一NS方程组的常用方法可以彼此检验和补充.(2)流体之间以及流体与外界的动量、能量和质量交换,流态从层流到湍流的演化主要发生在干扰剪切流动中,干扰剪切流及其最简单控制方程--DPNS方程组具有基础意义;DPNS方程组笔者在1967年已提出.(3)诸简化NS方程组:DPNS、抛物化(P)NS、薄层(TL)NS、黏性层(VL)NS方程组的发展、相互关系,它们的历史贡献和今后的用途;它们的数学性质均为扩散抛物型,但它们包含的黏性项彼此有所不同;从流体力学角度来看,它们中只有DPNS方程组能够准确描述干扰剪切流动.提出把诸简化NS方程组统一为DPNS方程组的建议.(4)干扰剪切流--DPNS方程组与无干扰剪切流--边界层方程组之间的关系以及进一步研究干扰剪切流的意义.

Effects of the Spanwise Characteristic Length on the Transition Feature in the Wake of a Cylinder

The transition 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 of 3-D incompressible Navier-Stokes equations. The spanwise characteris-tic length determines the transition features and global properties of the wake.

Direct Numerical Simulation of Supersonic Turbulent Boundary Layer Flow

Direct numerical simulations of a spatially evolving supersonic flat-plate turbulent boundary layer flow with free Mach number M = 2.25 and Reynolds number Re = 365000/in are performed. The transition process from laminar to turbulent flow is obtained by solving the three-dimensional compressible Navier-Stokes, equations, using high-order accurate difference schemes. The obtained statistical results agree well with the experimental and theoretical data. From the numerical results it can be seen that the transition process under the considered conditions is the process which skips the Tolimien-Schlichting instability and the second instability through the instability of high gradient shear layer and becomes of laminar flow breakdown. This means that the transition process is a bypass-type transition process. The spanwise asymmetry of the disturbance locally upstream imposed is important to induce the bypass-type transition. Furthermore, with increasing the time disturbance frequency the transition will delay. When the time disturbance frequency is large enough, the transition will disappear.

Shock-Induced Near-Wall Two-Phase Flow Structure Over a Micron-Sized Particles Bed

The present paper studies numerical modelling of near-wall two-phase flows induced by a normal shock wave moving at a constant speed, over a micronsized particles bed. In this two-fluid model, the possibility of particle trajectory intersection is considered and a full Lagrangian formulation of the dispersed phase is introduced. The finiteness of the Reynolds and Mach numbers of the flow around a particle as well as the fineness of the particle sizes are taken into account in describing the interactions between the carrier- and dispersed- phases. For the small mass-loading ratio case, the numerical simulation of flow structure of the two phases is implemented and the profiles of the particle number density are obtained under the constant-flux condition on the wall. The effects of the shock Mach number and the particle size and material density on particle entrainment motion are discussed in detail.The obtained results indicate that interphase non-equilibrium in the velocity and temperature is a common feature for this type of flows and a local particle accumulation zone may form near the envelope of the particle trajectory family.

Squeeze-film Effects in MEMS Devices with Perforated Plates for Small Amplitude Vibration

Squeeze-film effects of perforated plates for small amplitude vibration are analyzed through modified Reynolds equation (MRE). The analytical analysis reckons in most important influential factors: compressibility of the air, border effects, and the resistance caused by vertical air flow passing through perforated holes. It is found that consideration of air compressibility is necessary for high operating frequency and small ratio of the plate width to the attenuation length. The analytical results presented in this paper agree with ANSYS simulation results better than that under the air incompressibility assumption. The analytical analysis can be used to estimate the squeeze-film effects causing damping and stiffness added to the system. Since the value of Reynolds number involved in this paper is low (< 1), inertial effects are neglected.

Flows with Inertia in a Three-Dimensional Random Fiber Network

A fiber web is modeled as a three-dimensional random cylindrical fiber network. Nonlinear behavior of fluid flowing through the fiber network is numerically simulated by using the lattice Boltzmann (LB) method. A nonlinear relationship between the friction factor and the modified Reynolds number is clearly observed and analyzed by using the Fochheimer equation, which includes the quadratic term of velocity. We obtain a transition from linear to nonlinear region when the Reynolds numbers are sufficiently high, reflecting the inertial effect of the flows. The simulated permeability of such fiber network has relatively good agreement with the experimental results and finite element simulations.

Heating characteristics of blunt swept fin-induced shock wave turbulent boundary layer interaction

An experimental study was conducted on shock wave turbulent boundary layer interactions caused by a blunt swept fin-plate configuration at Mach numbers of 5.0, 7.8, 9.9 for a Reynolds number range of (1.0.similar to 4.7) x 10(7)/m. Detailed heat transfer and pressure distributions were measured at fin deflection angles of up to 30 degrees for a sweepback angle of 67.6 degrees. Surface oil flow patterns and liquid crystal thermograms as well as schlieren pictures of fin shock shape were taken. The study shows that the flow was separated at deflection of 10 degrees and secondary separation were detected at deflection of theta greater than or equal to 20 degrees. The heat transfer and pressure distributions on flat plate showed an extensive plateau region followed by a distinct dip and local peak close to the fin foot. Measurements of the plateau pressure and heat transfer were in good agreement with existing prediction methods, but pressure and heating peak measurements at M greater than or equal to 6 were significantly lower than predicted by the simple prediction techniques at lower Mach numbers.