Large eddy simulation of particle deposition and resuspension in turbulent duct flows

Particle deposition and resuspension in a horizontal, fully developed turbulent square duct flow at four flow bulk Reynolds numbers (10,320, 83k, 215k and 250k) is simulated by applying large eddy simulation coupled with Lagrangian particle tracking technique. Forces acting on particles includes drag, lift, buoyancy and gravity. Four particle sizes are considered with the diameters of 5?μm, 50?μm, 100?μm and 500?μm. Results obtained for the fluid phase are in good agreement with the available experimental and numerical data. Predictions for particles show that particle size, flow Reynolds number and the duct (celling, floor and vertical) walls play important roles in near-wall particle deposition and resuspension. For the smallest particle (5?μm), the particle deposition rates in duct ceiling, floor and vertical walls are found to be similar with each other and all increase with the flow Reynolds number, while the particle resuspension tends to occur in the middle wall regions and corners of the duct with less influenced by the flow Reynolds number. The ceiling deposition rate gradually decreases with particle size while the floor and vertical wall deposition rates both increase with particle size. The ceiling particle deposition rate increases with Reynolds number while the floor deposition rate decreases with it. The vertical deposition rate for the small particles (5–50?μm) increases with the flow Reynolds number obviously, while for the large particles (100–500?μm) that becomes insensitive. In addition, the flow Reynolds number is found to have an obvious effect on particle resuspension while the effect of particle size on particle resuspension decreases with Reynolds number. Eventually, a dynamic analysis was conducted for particles deposition and resuspension in turbulent duct flows.     

Visualization and analysis of longitudinal vortices at curved walls of 2D laminar and turbulent channel flows

The three-dimensional structure of longitudinal vortices at the curved walls of both laminar and turbulent channel flows is visualized by the hydrogen bubble technique. Together with the conditional sampling of the turbulent characteristics at the wall processed by thevita method the dynamics of the near-wall structure is discussed.     

颗粒弥散强化复相陶瓷临界颗粒尺寸的优化设计

本文应用弹性力学和Ｇｒｉｆｆｉｔｈ微裂纹理论，对颗粒弥散强化陶瓷材料内的热应力进行了分析，导出了临界粒尺寸的计算公式，为颗粒弥散强化复相陶瓷材料的界面设计和颗粒尺寸选择提供了依据。     

高压喷嘴空化初生的大涡模拟

考虑瞬时压力脉动对空化初生的影响,采用动态亚格子应力模型对某高压柴油喷嘴内部瞬时流动进行大涡模拟。计算结果表明：喷孔入口分离形成的局部低压区内时均湍流脉动压力可达时均静压的2倍,两相流场时均空泡位置及形态与实验结果十分吻合,并与基于时均静压及湍流脉动压力预测的空化初生区域基本一致。分离剪切层失稳形成柱状涡并激发了转捩过程,在分离转捩区域柱状涡发生扭曲变形并产生强烈的瞬时压力脉动,从而导致了空化。喷孔入口倒角处理使得分离减弱,对分离转捩过程和空化初生具有抑制作用。     

Application of flamelet/progress-variable approach to the large eddy simulation of a turbulent jet flame of pulverized coals

In this study, the flamelet/progress-variable (FPV) approach was applied to a large eddy simulation of a pulverized coal jet flame. The FPV approach considers the characteristics of the pulverized coal flame, e.g., non-adiabatic system and several types of fuel streams, via additional representative variables. First, the applicability of the FPV approach to a turbulent flame with pulverized coals was confirmed through a comparison of the numerical solutions and experimental data. In this study, the pure pilot case was also investigated to clarify the effects of pulverized coals on the flame. The flame structure changes significantly upon the injection of pulverized coals, and the flame index suggests the coexistence of premixed and diffusion combustion modes even in the downstream region. In particular, the combustion mode fluctuates with time in the middle region of the flame. The fuel gas released from the pulverized coals should increase in this region; therefore, the release and combustion behavior of the volatile matter must be involved in the combustion mode variation. The evaluation of the combustion modes of fuel gas in the coal flame is useful for the design and optimization of pulverized coal combustors with next-generation technologies.     

Note on eddy formation for stokes flows in circular cavities

Abstract

In this note, a simple model of a circular cavity with two arcs rotated in opposite directions on the boundary is employed to demonstrate flow patterns of recirculation. A particular streamline “cutting” the domain into two separating eddies is formed as a circular arc. When the rotating angular speed of the arc is proportional to the length of the arc, the streamline degenerates to a straight line.     

Velocity Fluctuations in a Particle-Laden Turbulent Flow over a Backward-Facing Step

Dilute gas-particle turbulent flow over a backward-facing step is numerically simulated. Large Eddy Simulation (LES) is used for the continuous phase and a Lagrangian trajectory method is adopted for the particle phase. Four typical locations in the flow field are chosen to investigate the two-phase velocity fluctuations. Time-series velocities of the gas phase with particles of different sizes are obtained. Velocity of the small particles is found to be similar to that of the gas phase, while high frequency noise exists in the velocity of the large particles. While the mean and rms velocities of the gas phase and small particles are correlated, the rms velocities of large particles have no correlation with the gas phase. The frequency spectrum of the velocity of the gas phase and the small particle phase show the -5/3 decay for higher wave number, as expected in a turbulent flow. However, there is a "rising tail' in the high frequency end of the spectrum for larger particles. It is shown that large particles behave differently in the flow field, while small particles behave similarly and dominated by the local gas phase flow.     

Convection of heat and thermodynamic irreversibilities in two-phase,turbulent nanofluid flows in solar heaters by corrugated absorber plates

The effects of simultaneous implementation of corrugated walls and nanoparticles upon the performance of solar heaters are investigated. Triangular and sinusoidal wall profiles along with varying concentration of nanoparticles are analyzed. The multi-phase mixture and the SST κ-ω models are used to simulate turbulent nanofluid flows inside the corrugated channels. The staggered computational grid is employed for storing the velocity and pressure terms at cell faces and cell center, respectively. The governing equations are first discretized by employing a second-order upwind differencing technique and are then solved by means of pressure-based finite volume approach. The convergence criterion is also presented for the validation of obtained results. The effects of wall profiles and nanoparticle concentration on the pertinent parameters including Nusselt number, pressure drop, performance evaluation criterion (PEC), and thermal and frictional irreversibilities are studied. This reveals that, in general, the triangular duct features superior heat transfer and inferior hydraulic characteristics in comparison with the sinusoidal duct. It is demonstrated that as long as the base fluid (water) is used the highest value of PEC corresponds to the straight duct. Yet, by introducing nanofluids the PEC values of the corrugated ducts exceed those of the straight duct. The analysis further shows that on the basis of the performance evaluation criterion, the sinusoidal duct appears to be a better choice in comparison with the triangular duct. However, the situation is reversed when thermodynamic irreversibilities are considered. It is argued that vortex formation in the two investigated wavy walls and shear layer developed in the triangular case are the essential physical reasons for the observed thermal, hydraulic and entropic behaviors.     

Capturing non-Brownian particle transport in periodic flows

In this paper we computationally examine the motion of a dilute suspension of slightly non-neutrally buoyant solid spheres as they migrate across the curved fluid streamlines of a viscous cellular flow. This is done by incorporating particle-fluid interactions into a continuum-based Lagrangian advection model derived from the Basset–Boussinesq–Oseen (BBO) equation, where the flow field is mimicked by using a perturbed streamfunction. Although the purely regular cellular flow is able to capture maximum velocity and particle diameter effects that are observed experimentally, it has several shortcomings. Most significantly, it is unable to capture the secondary island structures that exist in many rotating flow systems, nor the impact that these structures are observed to have on particle migration. Our results in this work demonstrate significant interplay between the underlying fluid structure and the non-trivial equilibrium locations of the non-Brownian particles, in agreement with previous experimental work. We also evaluate the effect of the Saffman lift force on the lateral migration of the solid spheres.     

Large eddy simulation of pollutant gas dispersion with buoyancy ejected from building into an urban street canyon

The dispersion of buoyancy driven smoke soot and carbon monoxide (CO) gas, which was ejected out from side building into an urban street canyon with aspect ratio of 1 was investigated by large eddy simulation (LES) under a perpendicular wind flow. Strong buoyancy effect, which has not been revealed before, on such pollution dispersion in the street canyon was studied. The buoyancy release rate was 5 MW. The wind speed concerned ranged from 1 to 7.5 m/s. The characteristics of flow pattern, distribution of smoke soot and temperature, CO concentration were revealed by the LES simulation. Dimensionless Froude number (Fr) was firstly introduced here to characterize the pollutant dispersion with buoyancy effect counteracting the wind. It was found that the flow pattern can be well categorized into three regimes. A regular characteristic large vortex was shown for the CO concentration contour when the wind velocity was higher than the critical re-entrainment value. A new formula was theoretically developed to show quantitatively that the critical re-entrainment wind velocities, u_c, for buoyancy source at different floors, were proportional to −1/3 power of the characteristic height. LES simulation results agreed well with theoretical analysis. The critical Froude number was found to be constant of 0.7.     

On the longitudinal dispersion of heavy particles in a horizontal circular pipe

The problem of longitudinal dispersion of heavy particles of finite size in a horizontal circular pipe carrying a steady turbulent flow was incompletely dealt with in the classical work of Elder [J. Fluid Mech. 5 (1959) 544]. In this paper, the problem is re-examined and a formal expression in terms of a cell function for the dispersion coefficient is deduced using an asymptotic analysis. The dispersion coefficient can be calculated upon finding numerically the first few components of the Fourier expansion of the cell function. Prediction by the present theory is found to agree better than Elder's theory with some published experimental data.     

Behaviour of organised disturbances in fully developed turbulent channel flow

In our earlier work we have shown the relevance of stability theory in understanding the sustenance of turbulence in turbulent boundary layers. Here we adopt the same model to study the evolution of organised disturbances in turbulent channel flow. Since the dominant modes are wall modes we find that the stability characteristics in the two flows are nearly identical although the boundary conditions (at the edge of the boundary layer and at the centre of the channel) are different. Comparisons with the experiments of Hussain and Reynolds are also presented.     

Wavelet diagnostics for detection of coherent structures in instantaneous turbulent flow imagery: A review

A review of work over the last decade shows that 2D wavelet techniques applied on flow imagery can provide powerful insights into the nature and lifecycle of coherent structures (the latter through wavelet movies) in turbulent shear flows. The advantage of wavelet techniques in often being able to infer the nature of coherent motion from a single image is emphasized. The techniques are first calibrated by using them on well-known results in the turbulent mixing layer. They are then applied to jets and plumes, and it is shown how off-source heating in such flows can disrupt the coherent structures in the unheated flow. A suitably reduced version of the present method, using discrete wavelet transforms on signals from a finite array of sensors, could be a useful diagnostic tool in near-real-time detection of coherent structures or patterns for the purpose of selecting appropriate control signals to the actuators in a flow-control system.     

Investigation of turbulent heat transfer and nanofluid flow in a double pipe heat exchanger

In present study, heat transfer and turbulent flow of water/alumina nanofluid in a parallel as well as counter flow double pipe heat exchanger have been investigated. The governing equations have been solved using an in-house FORTRAN code, based on finite volume method. Single-phase and standard k-ε models have been used for nanofluid and turbulent modeling, respectively. The internal fluid has been considered as hot fluid (nanofluid) and the external fluid, cold fluid (base fluid). The effects of nanoparticles volume fraction, flow direction and Reynolds number on base fluid, nanofluid and wall temperatures, thermal efficiency, Nusselt number and convection heat transfer coefficient have been studied. The results indicated that increasing the nanoparticles volume fraction or Reynolds number causes enhancement of Nusselt number and convection heat transfer coefficient. Maximum rate of average Nusselt number and thermal efficiency enhancement are 32.7% and 30%, respectively. Also, by nanoparticles volume fraction increment, the outlet temperature of fluid and wall temperature increase. Study the minimum temperature in the solid wall of heat exchangers, it can be observed that the minimum temperature in counter flow has significantly reduced, compared to parallel flow. However, by increasing Reynolds number, the slope of thermal efficiency enhancement of heat exchanger gradually tends to a constant amount. This behavior is more obvious in parallel flow heat exchangers. Therefore, using of counter flow heat exchangers is recommended in higher Reynolds numbers.     

颗粒流动模型研究进展

颗粒流动在自然界中和各种工业过程中广泛存在,但人们对于其机理认识的还不深入。描述颗粒运动的模型有很多,连续介质模型应用简单但准确性比较低,离散微粒学模型是近来人们研究的一个热点,以每个颗粒为考察对象,能够更准确地反应颗粒系统的性质。本文介绍了描述颗粒流动的模型,概述了各模型的理论和应用,通过对多种模型的比较可以看到,每个模型都有一定的使用范围,要更准确、更方便地描述颗粒系统的运动,还要进行深入地研究。     

Particle electrification and levitation in a continuous particle feed and dispersion system with vibration and external electric fields

Electrification and levitation of particles in a continuous particle feed and dispersion system have been studied both theoretically and experimentally. This system consisted of a vibrator and inclined parallel electrodes. A mesh and a vibrating plate were used for the upper and lower electrodes, respectively. A dc voltage was applied to one of the electrodes and the other electrode was grounded. Particles fed to the lower electrode were charged by induction and levitated upward by the Coulomb forces. When the applied voltage was high enough, the particles passed through the mesh electrode. The charge of the particles was measured with a Faraday cup, and the particle behavior was observed with a high-speed microscope camera. The particle charges were also analyzed from experimentally obtained particle trajectories and numerically calculated electric fields. Finally, the conditions for the effective levitation and dispersion of the charged particles and their mechanisms were studied and have been described in detail.     

Particle screening phenomena in an oblique multi-level tumbling reservoir: a numerical study using discrete element simulation

Due to their wide usage in industrial and technological processes, granular materials have captured great interest in recent research. The related studies are often based on numerical simulations and it is challenging to investigate computational phenomena of granular systems. Particle screening is an essential technology of particle separation in many industrial fields. This paper presents a numerical model for studying the particle screening process using the discrete element method that considers the motion of each particle individually. Dynamical quantities like particle positions, velocities and orientations are tracked at each time step of the simulation. The particular problem of interest is the separation of round shape particles of different sizes using a rotating tumbling vertical cylinder while the particulate material is continuously fed into its interior. This rotating cylinder can be designed as a uniform or stepped multi level obliqued vertical vessel and is considered as a big reservoir for the mixture of particulate material. The finer particles usually fall through the sieve openings while the oversized particles are rebounded and ejected through outlets located around the machine body. Particle–particle and particle–boundary collisions will appear under the tumbling motion of the rotating structure. A penalty method, which employs spring-damper models, will be applied to calculate the normal and frictional forces. As a result of collisions, the particles will dissipate kinetic energy due to the normal and frictional contact losses. The particle distribution, sifting rate of the separated particles and the efficiency of the segregation process have been studied. It is recognized that the screening phenomenon is very sensitive to the machines geometrical parameters, i.e. plate inclinations, shaft eccentricities and aperture sizes in the sieving plates at different levels of the structure. The rotational speed of the machine and the feeding rate of the particles flow have also a great influence on the transportation and segregation rates of the particles. In an attempt to better understand the mechanism of the particle transport between the different layers of the sifting system, different computational studies for achieving optimal operation have been performed.     

A priori evaluation of dynamic subgrid models of turbulence in square duct flow

A priori tests of two dynamic subgrid-scale (SGS) turbulence models have been performed using a highly resolved direct numerical simulation database for the case of turbulent incompressible flow in a straight duct of square cross-section. The model testing is applied only to the homogeneous flow direction where grid filtering can be applied without the introduction of commutation errors. The first SGS model is the dynamic (Smagorinsky/eddy viscosity) SGS model (DSM) developed by Germano et al. [1] while the second is the dynamic two parameter (mixed) model (DTM) developed by Salvetti and Banerjee [2]. For the Smagorinsky model we have used both the Fourier cut-off filter and a modified Gaussian filter which has the property that it removes aliasing errors in consistent a priori model-testing for spectral-based datasets. Results largely consistent with those found for plane channel flow are observed but with some slight differences in the corner regions. As found in prior studies of this sort, there is a very poor correlation of the modelled and exact subgrid-scale dissipation in the case of the DSM. The DSM over-predicts subgrid-scale dissipation on average. Instantaneously, the model provides an inaccurate representation of subgrid-scale dissipation, in general underestimating the magnitude by approximately one order of magnitude. On the other hand, the DTM shows excellent agreement with the exact SGS dissipation over most of the duct cross-section with a correlation coefficient of approximately 09.