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.

Variational approach to the closure problem of turbulence theory

A new method is proposed to solve the closure problem of turbulence theory and to drive the Kolmogorov law in an Eulerian framework. Instead of using complex Fourier components of velocity field as modal parameters, a complete set of independent real parameters and dynamic equations are worked out to describe the dynamic states of a turbulence. Classical statistical mechanics is used to study the statistical behavior of the turbulence. An approximate stationary solution of the Liouville equation is obtained by a perturbation method based on a Langevin-Fokker-Planck (LFP) model. The dynamic damping coefficient eta of the LFP model is treated as an optimum control parameter to minimize the error of the perturbation solution; this leads to a convergent integral equation for eta to replace the divergent response equation of Kraichnan's direct-interaction (DI) approximation, thereby solving the closure problem without appealing to a Lagrangian formulation. The Kolmogorov constant Ko is evaluated numerically, obtaining Ko = 1.2, which is compatible with the experimental data given by Gibson and Schwartz, (1963).

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.

Subgrid scale fluid velocity timescales seen by inertial particles in large-eddy simulation of particle-laden turbulence

Large-eddy simulation (LES) has emerged as a promising tool for simulating turbulent flows in general and, in recent years,has also been applied to the particle-laden turbulence with some success (Kassinos et al., 2007). The motion of inertial particles is much more complicated than fluid elements, and therefore, LES of turbulent flow laden with inertial particles encounters new challenges. In the conventional LES, only large-scale eddies are explicitly resolved and the effects of unresolved, small or subgrid scale (SGS) eddies on the large-scale eddies are modeled. The SGS turbulent flow field is not available. The effects of SGS turbulent velocity field on particle motion have been studied by Wang and Squires (1996), Armenio et al. (1999), Yamamoto et al. (2001), Shotorban and Mashayek (2006a,b), Fede and Simonin (2006), Berrouk et al. (2007), Bini and Jones (2008), and Pozorski and Apte (2009), amongst others. One contemporary method to include the effects of SGS eddies on inertial particle motions is to introduce a stochastic differential equation (SDE), that is, a Langevin stochastic equation to model the SGS fluid velocity seen by inertial particles (Fede et al., 2006; Shotorban and Mashayek, 2006a; Shotorban and Mashayek, 2006b; Berrouk et al., 2007; Bini and Jones, 2008; Pozorski and Apte, 2009).However, the accuracy of such a Langevin equation model depends primarily on the prescription of the SGS fluid velocity autocorrelation time seen by an inertial particle or the inertial particle–SGS eddy interaction timescale (denoted by $\delt T_{Lp}$ and a second model constant in the diffusion term which controls the intensity of the random force received by an inertial particle (denoted by C_0, see Eq. (7)). From the theoretical point of view, dTLp differs significantly from the Lagrangian fluid velocity correlation time (Reeks, 1977; Wang and Stock, 1993), and this carries the essential nonlinearity in the statistical modeling of particle motion. dTLp and C0 may depend on the filter width and particle Stokes number even for a given turbulent flow. In previous studies, dTLp is modeled either by the fluid SGS Lagrangian timescale (Fede et al., 2006; Shotorban and Mashayek, 2006b; Pozorski and Apte, 2009; Bini and Jones, 2008) or by a simple extension of the timescale obtained from the full flow field (Berrouk et al., 2007). In this work, we shall study the subtle and on-monotonic dependence of $\delt T_{Lp}$ on the filter width and particle Stokes number using a flow field obtained from Direct Numerical Simulation (DNS). We then propose an empirical closure model for $\delta T_{Lp}$. Finally, the model is validated against LES of particle-laden turbulence in predicting single-particle statistics such as particle kinetic energy. As a first step, we consider the particle motion under the one-way coupling assumption in isotropic turbulent flow and neglect the gravitational settling effect. The one-way coupling assumption is only valid for low particle mass loading.

海底沉积层中水合物成藏过程模拟

水合物成藏是一种动态演变过程，从生成和分解速度相对大小角度提出了天然气水合物的动态成藏理论；分析指出天然气水合物储层状态主要存在三种类型：成长型（包括渗透型与扩散型）、成熟型和消退型，必须结合经济型指标对这三种类型水合物藏进行开采可行性研究。 通过在自制的透明鼓泡反应釜内模拟气体迁移与水合物形成实验，解释了气体迁移对体系内传热与传质、溶解性质、过压等的影响因，进一步指出需要研究新的水合物相图。 从传热角度建立了水合物生长模型，计算并比较了南海北部陆坡琼东南盆地里甲烷水合物在渗漏系统和扩散系统下的生长速度，指出前者比后者快约20-40倍。 建立了水-气（游离气）-水合物-盐反应的多相流模型，演绎了在四个不同时刻随着水合物的形成，溶解度、渗透率、毛细压力、饱和度及盐度等的联动演变关系，并分析和比较了在沉积层水合物稳定区域内不同游离气饱和度和甲烷渗漏通量情况下水合物的形成过程和三种成藏类别。 反演了神狐海域水合物的形成，指出该区很可能是属于发育阶段的成长型渗漏系统水合物储藏，按照甲烷渗漏通量为0.5kg/m2·a计算，该地区已演化了约4800a，还需要约2700ka孔隙中水合物可到最大体积饱和度，达约75%，并在此后游离气将溢出沉积层进入海水中；计算得最终的甲烷水合物储量约为1.47×109m3，即该区最终将储存着约2200亿立方米甲烷气，是现在预测的约13倍！

黄土坡面水力学特征参数与土壤侵蚀量间关系研究

通过室内人工降雨模拟试验,研究黄土坡面水流的水力学特征参数与土壤侵蚀量间的关系。结果表明:随着放水流量的增加,出口断面平均流速、雷诺数在逐渐增大,弗劳德数在逐渐减小,坡面流流态也从层流状态逐渐过渡到紊流状态,当放水流量在10.1~19.9 L/min范围内时,坡面水流的流态为层流与紊流之间的混和状态;出口断面的平均流速随着放水流量的增加而增大,二者之间呈良好的幂函数关系;黄土坡面的累计泥沙侵蚀总量随着雷诺数的增加而增大,二者呈幂函数关系;单宽土壤侵蚀产沙量与雷诺数之间拟合关系的好坏受雷诺数取值的影响,当雷诺数的值在临界雷诺数(取500)附近时,单宽土壤侵蚀产沙量与雷诺数之间的负相关性良好。研究结果不仅可以揭示土壤侵蚀水动力过程的内在机制,而且对黄土高原地区生态环境建设具有重要的现实意义。

Instabilities in a non-transferred direct current plasma torch operated at reduced pressure

Using an oscilloscope, a high-speed video camera and a double-electrostatic probe system, the periodicity and amplitude of the fluctuations in arc voltage, jet luminance and ion saturation current of a plasma jet were monitored to investigate various sources of instabilities and their effects in a non-transferred dc plasma torch operated at reduced pressure. The results show that besides a 300 Hz main fluctuation inherited from the power supply, arc voltage fluctuation of 3–4 kHz with an amplitude less than 5% of the mean voltage was mainly affected by the total gas flow rate. The arc voltage fluctuation can affect the energy distribution of the plasma jet which is detectable by electrostatic probes and a high-speed video camera. The steadiness of energy transfer is also affected by the laminar or turbulent flow state of the plasma.

CFD study of a sudden-expanding coal combustor using Euler-Euler and Euler-Lagrange models

In this study, the Euler-Euler (E-E) and Euler-Lagrange (E-L) models designed for the same chemical mechanism of heterogeneous reactions were used to predict the performance of a typical sudden-expanding coal combustor. The results showed that the current E-E model underestimated the coal burnout rate because the particle temperature fluctuation on char combustion is not adequately considered. A comparison of the E-E and E-L simulations showed the underestimation of heterogeneous chemical reaction rates by the E-E model. (C) 2010 Elsevier Ltd. All rights reserved.

Air flow over a mountain and the convective boundary-layer

In this paper, the analytical model coupling the convective boundary layer (CBL) with the free atmosphere developed by Qi and Fu (1992) is improved. And by this improved model, the interaction between airflow over a mountain and the CBL is further discussed. The conclusions demonstrate: (1) The perturbation potential temperatures in the free atmosphere can counteract the effect of orographic thermal forcing through entraining and mixing in the CBL. If u(M)BAR > u(F)BAR, the feedback of the perturbation potential temperatures in the free atmosphere is more important than orographic thermal forcing, which promotes the effect of interfacial waves. If u(M)BAR < u(F)BAR, orographic thermal forcing is more important, which makes the interfacial height and the topographic height identical in phase, and the horizontal speeds are a maximum at the top of the mountain. (2) The internal gravity waves propagating vertically in the free atmosphere cause a strong downslope wind to become established above the lee slope in the CBL and result in the hydraulic jump at the top of the CBL. (3) With the CBL deepening, the interfacial gravity waves induced by the potential temperature jump at the top of the CBL cause the airflow in the CBL to be subcritical.

An experimental study on turbulent coherent structures near a sheared air-water interface

The turbulence structures near a sheared air-water interface were experimentally investigated with the hydrogen bubble visualization technique. Surface shear was imposed by an airflow over the water flow which was kept free from surface waves. Results show that the wind shear has the main influence on coherent structures under air-water interfaces. Low- and high- speed streaks form in the region close to the interface as a result of the imposed shear stress. When a certain airflow velocity is reached, "turbulent spots" appear randomly at low-speed streaks with some characteristics of hairpin vortices. At even higher shear rates, the flow near the interface is dominated primarily by intermittent bursting events. The coherent structures observed neat sheared air-water interfaces show qualitative similarities with those occurring in near-wall turbulence. However, a few distinctive phenomena were also observed, including the fluctuating thickness of the instantaneous boundary layer and vertical vortices in bursting processes, which appear to be associated with the characteristics of air-water interfaces.

Simulation of Free-Surface Flow in a Tank Using the Navier-Stokes Model and Unstructured Finite Volume Method

Modelling free-surface flow has very important applications in many engineering areas such as oil transportation and offshore structures. Current research focuses on the modelling of free surface flow in a tank by solving the Navier-Stokes equation. An unstructured finite volume method is used to discretize the governing equations. The free surface is tracked by dynamically adapting the mesh and making it always surface conforming. A mesh-smoothing scheme based on the spring analogy is also implemented to ensure mesh quality throughout the computaiton. Studies are performed on the sloshing response of a liquid in an elastic container subjected to various excitation frequencies. Further investigations are also carried out on the critical frequency that leads to large deformation of the tank walls. Another numerical simulation involves the free-surface flow past as submerged obstacle placed in the tank to show the flow separation and vortices. All these cases demonstrate the capability of this numerical method in modelling complicated practical problems.

Flow properties of a dusty-gas point source in a supersonic free stream

By using Lagrangian method, the flow properties of a dusty-gas point source in a supersonic free stream were studied and the particle parameters in the near-symmetry-axis region were obtained. It is demonstrated that fairly inertial particles travel along oscillating and intersecting trajectories between the bow and termination shock waves. In this region,formation of "multi-layer structure" in particle distribution with alternating low- and highdensity layers is revealed. Moreover, sharp accumulation of particles occurs near the envelopes of particle trajectories.

Modeling Study on the Entrainment of Ambient Air into Subsonic Laminar and Turbulent Argon Plasma Jets

Modeling study is performed to reveal the special features of the entrainment of ambient air into subsonic laminar and turbulent argon plasma jets. Two different types of jet flows are considered, i.e., the argon plasma jet is impinging normally upon a flat substrate located in atmospheric air surroundings or is freely issuing into the ambient air. It is found that the existence of the substrate not only changes the plasma temperature, velocity and species concentration distributions in the near-substrate region, but also significantly enhances the mass flow rate of the ambient air entrained into the jet due to the additional contribution to the gas entrainment of the wall jet formed along the substrate surface. The fraction of the additional entrainment of the wall jet in the total entrained-air flow rate is especially high for the laminar impinging plasma jet and for the case with shorter substrate standoff distances. Similarly to the case of cold-gas free jets, the maximum mass flow-rate of ambient gas entrained into the turbulent impinging or free plasma jet is approximately directly proportional to the mass flow rate at the jet inlet. The maximum mass flow-rate of ambient gas entrained into the laminar impinging plasma jet slightly increases with increasing jet-inlet velocity but decreases with increasing jet-inlet temperature.

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.

Receptivity To Free-Stream Disturbance Waves For Hypersonic Flow Over A Blunt Cone

A high-order shock-fitting finite difference scheme is studied and used to do direction numerical simulation (DNS) of hypersonic unsteady flow over a blunt cone with fast acoustic waves in the free stream, and the receptivity problem in the blunt cone hypersonic boundary layers is studied. The results show that the acoustic waves are the strongest disturbance in the blunt cone hypersonic boundary layers. The wave modes of disturbance in the blunt cone boundary layers are first, second, and third modes which are generated and propagated downstream along the wall. The results also show that as the frequency decreases, the amplitudes of wave modes of disturbance increase, but there is a critical value. When frequency is over the critial value, the amplitudes decrease. Because of the discontinuity of curvature along the blunt cone body, the maximum amplitudes as a function of frequencies are not monotone.