Large Eddy Simulation of Complex Turbulent Flows: Physical Aspects and Research Trends

In the current paper, we have primarily addressed one powerful simulation tool developed during the last decades-Large Eddy Simulation (LES), which is most suitable for unsteady three-dimensional complex turbulent flows in industry and natural environment. The main point in LES is that the large-scale motion is resolved while the small-scale motion is modeled or, in geophysical terminology, parameterized. With a view to devising a subgrid-scale(SGS) model of high quality, we have highlighted analyzing physical aspects in scale interaction and-energy transfer such as dissipation, backscatter, local and non-local interaction, anisotropy and resolution requirement. They are the factors responsible for where the advantages and disadvantages in existing SGS models come from. A case study on LES of turbulence in vegetative canopy is presented to illustrate that LES model is more based on physical arguments. Then, varieties of challenging complex turbulent flows in both industry and geophysical fields in the near future-are presented. In conclusion; we may say with confidence that new century shall see the flourish in the research of turbulence with the aid of LES combined with other approaches.

Diffusive Wave Solutions for Open Channel Flows with Uniform and Concentrated Lateral Inflow

The diffusive wave equation with inhomogeneous terms representing hydraulics with uniform or concentrated lateral inflow intoa river is theoretically investigated in the current paper. All the solutions have been systematically expressed in a unified form interms of response function or so called K-function. The integration of K-function obtained by using Laplace transform becomesS-function, which is examined in detail to improve the understanding of flood routing characters. The backwater effects usuallyresulting in the discharge reductions and water surface elevations upstream due to both the downstream boundary and lateral infloware analyzed. With a pulse discharge in upstream boundary inflow, downstream boundary outflow and lateral inflow respectively,hydrographs of a channel are routed by using the S-functions. Moreover, the comparisons of hydrographs in infinite, semi-infiniteand finite channels are pursued to exhibit the different backwater effects due to a concentrated lateral inflow for various channeltypes.

Micro-Scale Gas Flows

了解微尺度气体流动特点是微机电系统设计和优化的基础.有关的研究可以上溯到20世纪初Knudsen的平面槽道流动质量流量的测量和Millikan的小球阻力系数的测量,实验结果揭示了稀薄气体效应即尺度效应对气体运动的重要影响.由于流动特征长度很小,微尺度气流经常处于滑流区甚至过渡领域,流动的相似参数为Knudsen数和Mach数.因此可以考虑利用相似准则,通过增大几何尺寸、减小压力的途径,解决微机电系统实验观测遇到的困难.为解决直接模拟MonteCarlo方法分析微机电系统中低速稀薄气流遇到的统计涨落困难,我们提出了信息保存法(IP),该方法能够有效克服统计散布,并已成功用于多种微尺度气流.

Fluid-Solid Interaction in Particle-Laden Flows

Liu Qingquan, Singh V.P

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.

Kinetic Assessment of Measured Mass Flow Rates and Streamwise Pressure Distributions in Microchannel Gas Flows

Measured mass flow rates and streamwise pressure distributions of gas flowing through microchannels were reported by many researchers. Assessment of these data is crucial before they are used in the examination of slip models and numerical schemes, and in the design of microchannel elements in various MEMS devices. On the basis of kinetic solutions of the mass flow rates and pressure distributions in microchannel gas flows, the measured data available are properly normalized and then are compared with each other. The 69 normalized data of measured pressure distributions are in excellent agreement, and 67 of them are within 1 +/- 0.05. The normalized data of mass flow-rates ranging between 0.95 and 1 agree well with each other as the inlet Knudsen number Kn (i) < 0.02, but they scatter between 0.85 and 1.15 as Kn (i) > 0.02 with, to some extent, a very interesting bifurcation trend.

Sand-Dust Storm and Two-Phase Flows

从稀相气固两相流理论出发,针对沙尘暴问题的特点,给出描述固相拟流体本构关系和气固相间相互作用的方法,探讨确定沙尘悬浮临界判据的途径,导出多场耦合下含尘大气运动的基本方程,从而可为定量预报沙尘暴系统结构特征和长距离输送传播提供理论基础。

Influence of Coriolis force on bulk flows during horizontal Bridgman growth on a centrifuge: a numerical study

Three-dimensional and time-dependent numerical simulations are performed For melt convection in horizontal Bridgman crystal growth tinder high gravity conditions by means of a centrifuge. The numerical results show that Coriolis Force can cause a stabilizing effect on the fluctuations of the melt flow under a specific relation direction and relation rates of the centrifuge as reported in previous experiments (Ma et al., Materials Processing in High Gravity, Plenum Press, New York, 1994, p. 61). The present simulation provides details of the now features associated with the effect of the Coriolis force. There are also some differences between the present three-dimensional and former two-dimensional numerical solutions particularly in the prediction of the critical conditions and flow patterns.

First-passage failure and its feedback minimization of quasi-partially integrable Hamiltonian systems

An n degree-of-freedom Hamiltonian system with r (1¡r¡n) independent 0rst integrals which are in involution is calledpartially integrable Hamiltonian system. A partially integrable Hamiltonian system subject to light dampings andweak stochastic excitations is called quasi-partially integrable Hamiltonian system. In the present paper, the procedures for studying the 0rst-passage failure and its feedback minimization of quasi-partially integrable Hamiltonian systems are proposed. First, the stochastic averaging methodfor quasi-partially integrable Hamiltonian systems is brie4y reviewed. Then, basedon the averagedIt ˆo equations, a backwardKolmogorov equation governing the conditional reliability function, a set of generalized Pontryagin equations governing the conditional moments of 0rst-passage time and their boundary and initial conditions are established. After that, the dynamical programming equations and their associated boundary and 0nal time conditions for the control problems of maximization of reliability andof maximization of mean 0rst-passage time are formulated. The relationship between the backwardKolmogorov equation andthe dynamical programming equation for reliability maximization, andthat between the Pontryagin equation andthe dynamical programming equation for maximization of mean 0rst-passage time are discussed. Finally, an example is worked out to illustrate the proposed procedures and the e9ectiveness of feedback control in reducing 0rst-passage failure.

Numerical Simulation and Experimental Investigation about Internal and External Flows

In this paper, TASCflow3D is used to solve inner and outer 3D viscous incompressible turbulent flow (R-e = 5.6 X 10(6)) around axisymmetric body with duct. The governing equation is a RANS equation with standard k-epsilon turbulence model. The discrete method used is a finite volume method based on the finite element approach. In this method, the description of geometry is very flexible and at the same time important conservative properties are retained. The multi-block and algebraic multi-grid techniques are used for the convergence acceleration. Agreement between experimental results and calculation is good. It indicates that this novel approach can be used to simulate complex flow such as the interaction between rotor and stator or propulsion systems containing tip clearance and cavitation.

Shocked Flows Induced by Supersonic Projectiles Moving in Tubes

A numerical study on shocked flows induced by a supersonic projectile moving in tubes is described in this paper. The dispersion-controlled scheme was adopted to solve the Euler equations implemented with moving boundary conditions. Four test cases were carried out in the present study: the first two cases are for validation of numerical algorithms and verification of moving boundary conditions, and the last two cases are for investigation into wave dynamic processes induced by the projectile moving at Mach numbers of M-p = 2.0 and 2.4, respectively, in a short time duration after the projectile was released from a shock tube into a big chamber. It was found that complex shock phenomena exist in the shocked flow, resulting from shock-wave/projectile interaction, shock-wave focusing, shock-wave reflection and shock-wave/contact-surface interactions, from which turbulence and vortices may be generated. This is a fundamental study on complex shock phenomena, and is also a useful investigation for understanding on shocked flows in the ram accelerator that may provide a highly efficient facility for launching hypersonic projectiles.

Statistical Simulation of Rarefied Gas Flows in Micro-Channels

Rarefied gas flows through micro-channels are simulated using particle approaches, named as the information preservation (IP) method and the direct simulation Monte Carlo (DSMC) method. In simulating the low speed flows in long micro-channels the DSMC method encounters the problem of large sample size demand and the difficulty of regulating boundary conditions at the inlet and outlet. Some important computational issues in the calculation of long micro-channel flows by using the IP method, such as the use the conservative form of the mass conservation equation to guarantee the adjustment of the inlet and outlet boundary conditions and the super-relaxation scheme to accelerate the convergence process, are addressed. Stream-wise pressure distributions and mass fluxes through micro-channels given by the IP method agree well with experimental data measured in long micro-channels by Pong et al. (with a height to length ratio of 1.2:3000), Shih et al. (l.2:4800), Arkilic et al. and Arkilic (l.3:7500), respectively. The famous Knudsen minimum of normalized mass flux is observed in IP and DSMC calculations of a short micro-channel over the entire flow regime from continuum to free molecular, whereas the slip Navier-Stokes solution fails to predict it.

Probe for Measurements of Density/Conductivity in Flows of Conducting Fluids

A probe utilizing the bipolar pulse method to measure the density of a conducting fluid has been developed. The probe is specially designed such that the concentration of a stream tube can be sampled continuously. The density was determined indirectly from the measurement of solution conductivity. The probe was calibrated using standard NaCl solutions of varying molarity and was able to rapidly determine the density of a fluid with continuously varying conductance. Measurements of the conductivity profiles, corresponding density profiles, and their fluctuation levels are demonstrated in a channel flow with an electrolyte injected from a slot in one wall.